Interest in synthetic modeling of seismic wavefields scattered in the subsurface is growing today due to its applicability in various forward and inverse problems of geophysics. It has been used extensively for general evaluation of the subsurface structure, in survey design and illumination studies, and also as the basis for imaging and inversion algorithms. Models with complex geological structures that contain strong-contrast or irregular reflectors and shadow zones and where conventional algorithms fail to simulate realistic wavefields present a particular challenge.
