Abstract:SUMMARY The receiver function (RF) method is one of the most widely used passive seismic approaches for imaging subsurface structures from shallow sedimentary basins to the deep Earth’s interior. RF is typically computed with a deconvolution operation that enables separating structural response from source–time function embedded in teleseismic wavefields. However, seismic recordings of temporary arrays are often contaminated with strong random noise and even amplitude bias due to poor installation conditions. Both factors can cause instability to deconvolution and severely degrade the accuracy, thereby the imaging quality, of the RF method. This underscores the importance of effective denoising algorithms in RF processing. In this study, we explore the application of high-resolution Radon transform (RT) to improve the conventional RF imaging workflow. Contrary to the commonly implemented post-processing (i.e. after deconvolution) noise suppression schemes, we introduce the RT to data pre-processing (i.e. before deconvolution). This method seeks a sparse representation of teleseismic wavefields in the τ–p domain by iteratively solving a least-squares minimization problem with the conjugate gradient algorithm. Synthetic test with a 2-D step-Moho model shows that non-linear phase arrivals including incoherent noise and diffraction energy are effectively removed in resulting RFs, with the signal-to-noise ratio increased by as much as ∼8 dB. Real data experiments using the Hi-CLIMB network in the Tibetan Plateau demonstrate the superior performance of the proposed workflow in regularizing the wavefield and improving the coherence of converted phases across the recording array. Consequently, imaging results of common conversion point stacking using a single and a group of teleseismic events both recover subtle converted phases from the Moho and potential lithospheric discontinuities that are otherwise obscured by noise arising from the conventional processing workflow. This study highlights the necessity of wavefield regularization in the RF method and calls for improved data processing techniques in array-based seismic imaging.

Upside down Rayleigh-Marchenko: a practical, yet exact redatuming scheme for seabed seismic acquisitions

Reflected wave least squares reverse time migration with angle illumination compensation

Marchenko Equation for Time Reversal Imaging in an Underwater Acoustic Waveguide

Elastic Reverse Time Migration Without Wavefield Decomposition

Sequential Filtering Applied to Multi-Beam Echo-Sounding

Downward continuation of marine seismic reflection data: An undervalued tool to improve velocity models

Reflection Waveform Inversion Based on Full-Band Seismic Data Reconstruction for Salt Structure Inversion

Waveform Imaging Based on Linear Forward Representations for Scalar Wave Seismic Data

Virtual plane-wave imaging via Marchenko redatuming

Target-oriented least-squares reverse-time migration with Marchenko redatuming and double-focusing: Field data application

Accurate Source-Receiver Positioning Method for a High-Resolution Deep-Towed Multichannel Seismic Exploration System

Improving Receiver Function Imaging with High-Resolution Radon Transform

Target-Oriented Prestack Correlative Least Squares Reverse Time Migration Using Marchenko Redatumed Data

First-order multiples imaging aided by water bottom

Virtual sources and receivers in the real Earth -- considerations for practical applications

Generating Azimuth–Reflection Angle Gathers from Reverse Time Migration Using the High-Dimensional Local Phase Space Approximation of Seismic Wavefields

Marchenko Imaging from Rugged Topography in Mountainous Areas

Converted Wave Reverse Time Migration and Demigration Based on Forward Expression of Elastic Reflected Waves

Mesh-free Radial-Basis-function-generated Finite Differences and Their Application in Reverse Time Migration

Full Reciprocity-Gap Waveform Inversion in the frequency domain, enabling sparse-source acquisition