Accelerating Seismic Imaging Using Novel Boundary Condition Handling

In this paper, we present a load balancing method that attempts to avoid typical load balancing problems by removing the need for absorbing boundary condition calculations. We accomplish this by introducing an invertible, nonlinear transformation that maps the entire infinite model space into a finite space. This transform creates a finite halo of points around the modeling area of interest that models the infinite region of space beyond the limits of the model. The nonlinearity results in a slowing of wave propagation in the halo region and thereby provides a buffer zone that shields the model from unphysical reflections. Since the calculation is the same at every point of the model, load balancing is guaranteed.