Seismic scattering originated from structural heterogeneity covers a wide range of scales within the Earth's interior. Conventionally, stochastic approaches are employed to study high-frequency scattering process from random heterogeneity (Sato et al., 2012). One typical example is the characterization of P and S coda waves from local and regional earthquakes (Aki, 1969). On the other hand, deterministic imaging of subsurface structures has long been undertaken using either backward or forward scattering. For example, in earthquake seismology, receiver functions rely on the forward P-to-S scattering to image seismic discontinuities in the crust and mantle (Langston, 1979). Seismic reflection surveys utilize backward reflected waves to characterize petroleum reservoirs in exploration seismology and the Earth's crust in controlled-source crustal seismology (Prodehl & Mooney, 2012; Sheriff & Geldart, 1995). Unlike subhorizontal structures, strong lateral heterogeneity associated with near-vertical structure poses a significant challenge in deterministic scattered-wave imaging. In exploration seismology, complex structures can be imaged using densely distributed sources and receivers, in combination with sophisticated imaging or inversion techniques, such as reverse time migration or full-waveform inversion (Sheriff & Geldart, 1995; Virieux & Operto, 2009). However, such frameworks do not generally apply to earthquake seismology. The sparse distribution of earthquake sources and seismic stations greatly limits the detection ability of subsurface lateral heterogeneities, such as basin edges, fault zones, and Moho offsets, despite them having strong effects on seismic waveforms. Furthermore, many seismic imaging methods, such as common conversion point stacking of conventional receiver functions, assume subhorizontal structures, which cannot readily be applied to image strong lateral heterogeneities. Body-to-surface wave conversion is a special case of seismic scattering originated from strong lateral heterogeneity. Compared to body-to-body wave scattering, scattered surface waves propagate horizontally along