Abstract:Receiver function techniques are widely used in imaging crustal and mantle structure beneath a seismic station. The weak P-to-S conversions at deep seismic structures are usually masked by strong shallow reverberations when unconsolidated sediments are present below the station, making it nearly impossible to utilize receiver function techniques. We develop a method to estimate sediment and crustal structures beneath a seismic station based on wavefield downward continuation and decomposition method. The method parametrizes velocity structure beneath the station with a stack of constant velocity layers overlying a homogeneous half-space, and approximates the teleseismic P wave and its coda by the structural response to an incoming plane P wave. Our method is based on the principle that the upgoing S wavefield is absent in the half-space, and searches for the optimum velocity and thickness of the layers that give the minimum S-wave energy flux from the half-space to the layers. An iterative grid-search algorithm from the top to the bottom layers is employed to implement the search. In this study, we only use models comprising either only one crust layer or two layers (sediment + crust) with a half-space mantle, although models with more layers are also implementable. The method is especially useful in resolving seismic structure beneath a station sitting on unconsolidated sediments. It not only can be used to determine the sediment thickness and velocity structure, but also provides an effective way to generate subsurface receiver functions, which are formed by deconvolving the upgoing P wavefield from the upgoing S waves at the top of hardrock crust, and thus are free from shallow reverberations. The technique is applied to various synthetic data generated with different types of velocity model and noise levels, and appears to have good capability in recovering the input models. We further applied this method to teleseismic data recorded at a station inside the Songliao Basin in northeast China. The estimated sediment thickness and velocity agrees well with the results of previous active-source studies. The subsurface receiver functions also show a superior power in exposing the Moho Ps conversions, resulting in a well-defined peak in the H-kappa domain, which are absent in the regular receiver function data.

Simultaneous Structure Enhancing Filtering and Fault Detection Using Plane Wave Prediction

Omnidirectional Plane-Wave Destruction

Enhancing seismic reflections using empirical mode decomposition in the flattened domain

Nonlocal Structure Constrained Neighborhood Preserving Embedding Model and Its Application for Fault Detection

Empirical Mode Decomposition Based Smoothing in the Flattened Domain

Large Dip Calculation via Robust Nonstationary Plane-Wave Destruction

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

3-D Structural Complexity-Guided Predictive Filtering: A Comparison Between Different Non-Stationary Strategies

Accelerated Plane-Wave Destruction

Frequency Adaptive Fault Detection By Feature Pyramid Network With Wavelet Transform

An Optimal Spatio-Temporal Hybrid Model Based on Wavelet Transform for Early Fault Detection

Robust Nonstationary Local Slope Estimation

Physics-informed deep filtering of ultrasonic guided waves for incipient defect inspection of large-scale square tube structures

Warped P-SV wavelet distortion correction using a time-frequency adaptive shaping filter

Nonstationary local slope estimation via forward-backward space derivative calculation

Deblending of Simultaneous Source Data Using a Structure-Oriented Space-Varying Median Filter

3-D Time-Space Joint Deconvolution for Enhancing Recognition Accuracy of Seismic Microstructure

Geometry-informed phase space warping for reliable fatigue damage monitoring

Estimating Sedimentary and Crustal Structure Using Wavefield Continuation: Theory, Techniques and Applications

Structure-Oriented Singular Value Decomposition for Random Noise Attenuation of Seismic Data

An Adaptive Deconvolution Method with Improve Enhanced Envelope Spectrum and Its Application for Bearing Fault Feature Extraction