“…The search for more reliable seismic images and additional subsurface information, such as fracture distribution, drives advances in seismic acquisition, such as larger offsets, wider azimuths, and multicomponent recording. All of these advances facilitate incorporating anisotropy and elasticity into wavefield extrapolation and reverse time migration (RTM), which is the state-of-art wavefield imaging algorithm suitable for complex geologic structures (Baysal et al, 1983;Lailly, 1983;McMechan, 1983;Levin, 1984;Chang and McMechan, 1987;Hokstad et al, 1998;Farmer et al, 2009;Zhang and Sun, 2009). Although seismic acquisition improves with such advances, it always involves practical limitations, such as finite and irregular data sampling, which negatively impact wavefield im-aging methods.…”