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.

Estimation and Removal of the Residual Normal Moveout Based on the S- Transform

Elastic Reverse Time Migration Without Wavefield Decomposition

Seismic Data Two-Step Recovery Approach Combining Sparsity-Promoting and Hyperbolic Radon Transform Methods

Isotropic Elastic Reverse Time Migration

An Effective Scheme for Residual Moveout Picking Using a Markov Decision Process

Application of Seismic Gather Flattening Technique Based on Time-Frequency Analysis

Multiple-Eliminated Technique Based on Curvelet Transform

Seismic Linear Noise Attenuation Based on the Rotate-Time-Shift FK Transform

Accurate Elevation and Normal Moveout Corrections of Seismic Reflection Data on Rugged Topography

Adaptive Stabilization Forq-Compensated Reverse Time Migration

Substituting Smoothing with Low-Rank Decomposition - Applications to Least-Squares Reverse Time Migration of Simultaneous Source and Incomplete Seismic Data

Seismic Data Enhancement Based on Common-Reflection-Surface-Based Local Slope and Trimmed Mean Filter

A stable Q-compensated reverse-time migration method using a modified fractional viscoacoustic wave equation

Improving Receiver Function Imaging with High-Resolution Radon Transform

Structure-oriented prestack waveform inversion

Warped-mapping Based Multi-Gather Joint Prestack Q Estimation

Random Noise Attenuation Using a New Spectral Decomposition Method

Seismic trace reconstruction approach based on noise immunity weighted parabolic radon transform

A Multitask Learning-Based Dynamic Wavelet Amplitude Spectra Extraction Method and Its Application in Q Estimation

Application of spectral decomposition using regularized non-stationary autoregression to random noise attenuation

An Efficient and Effective Common Reflection Surface Stacking Approach Using Local Similarity and Plane-Wave Flattening