共查询到20条相似文献，搜索用时 312 毫秒
1.
用Wiener滤波方法提取台站接收函数 总被引：10，自引：1，他引：10
本文提出了一种在时间域用Wiener滤波方法提取台站接收函数的方法，用远震P波波形的垂直分量为输入，接收函数作为滤波因子，远震P波波形的径向和切向分量作为期望输出，通过期望输出与实际输出的均方误差达极小，来提取接收函数。接收函数的计算可归结为Toeplitz方程的求解，可以采用Levinson递推算法。Toeplitz方程的非奇异性保证了Wiener滤波反褶积方法的稳定性。合成地震图与观测地震图的检验表明，用Wiener滤波方法测定台站接收函数是一种有效的时间域反褶积方法。 相似文献
2.
Wiener filtering is used to estimate receiver function in a time-domain. With the vertical component of 3-component teleseismic P waveform as the input of a Wiener filter, receiver function as the filter response, and radial and tangential components as the expected output, receiver function is estimated by minimizing the error between expected and actual outputs. Receiver function can be obtained by solving the Toeplitz equation using the Leviuson algorithm. The non-singularity of the Toeplitz equation ensures the stability of Wiener Deconvolution. Both synthetic and observational seismogram checks show that Wiener Deconvolution is an effective time-domain method to estimate receiver function from teleseismic P waveform. 相似文献
3.
4.
Receiver function estimated by maximum entropy deconvolution 总被引：2，自引：0，他引：2
Maximum entropy deconvolution is presented to estimate receiver function, with the maximum entropy as the rule to determine auto-correlation and cross-correlation functions. The Toeplitz equation and Levinson algorithm are used to calculate the iterative formula of error-predicting filter, and receiver function is then estimated. During extrapolation, reflective coefficient is always less than I, which keeps maximum entropy deconvolution stable. The maximum entropy of the data outside window increases the resolution of receiver function. Both synthetic and real seismograms show that maximum entropy deconvolution is an effective method to measure receiver function in time-domain. 相似文献
5.
分析了相邻井下、地表地震计记录的背景噪声及近震、远震波形特征,以及相关的振幅谱及近震、远震接收函数波形特征。结果表明,井下地震计记录的背景噪声比地表地震计小1个数量级,其近震、远震波形记录较清晰、简洁;井下、地表地震计记录的远震振幅谱和径向接收函数具有较好的一致性,但背景噪声、近震振幅谱、近震接收函数、远震切向接收函数波形存在一定差异;由地表地震计得到的近震震级比井下地震计的大0.12,地表地震计到井下地震计之间的P波速度约为2 km/s。 相似文献
6.
利用Ｓｈａｍｗａｙ和Ｄｅｒ的多道最大或然性反褶积原理，讨论了从单台三分量远震Ｐ波波形中分离接收函数径向与垂向分量复谱比的方法．根据Ｔａｒａｎｔｏｌａ的波形反演理论发展了接收函数复谱比的非线性反演方法．本文的反演方法除了拟合接收函数的复谱比之外，还需拟合时间域中接收函数垂向与径向分量的初至振幅比．合成波形数据反演的结果表明本文方法的反演结果不依赖初始模型．利用该方法和ＣＤＳＮ台网兰州台记录的宽频带远震Ｐ波波形数据，研究了该台站的接收函数及其随方位角的变化．接收函数复谱比非线性反演给出了该台站下方１４０ｋｍ深度岩石层的Ｓ波速度结构，得到兰州台下方地壳厚度为５４ｋｓ，岩石层厚度为９４ｋｍ． 相似文献
7.
提出一种在时间域采用最大熵谱反褶积提取台站接收函数的方法，以最大熵作为自相关函数和互相关函数的递推准则，利用Toeplitz方程及Levinson递推算法，得到预测误差滤波系数的递推公式，从而计算台站接收函数．外推运算过程中，反射系数总是小于1，保证了最大熵谱反褶积的稳定性．时窗外数据熵极大提高了接收函数的分辨率．合成地震图与实测地震图的检验表明, 最大熵谱反褶积是一种在时间域测定台站接收函数的有效方法． 相似文献
8.
井下地震计波形记录的P波垂向分量存在频谱极小（spectrum null）现象,导致接收函数的结果不稳定.本文以首都圈地区为例,基于平面波入射的传播矩阵理论,发展了用于计算井下地震计的接收函数正演方法.在此基础上,分析了井下地震计波形垂向分量频谱极小现象,研究其对接收函数稳定性的影响.结果表明,井下地震计波形记录垂向分量的频谱极小开始出现的频率和地震计的埋深相关.该现象可造成反卷积提取的接收函数不稳定,且不稳定情况出现在频谱极小附近的频段,可通过选择合适的高斯因子压制其对接收函数的影响. 相似文献
9.
为研究青藏高原的岩石层构造及其动力学过程,根据记录到的来自台站东北方向的大量宽频带远震林波波形资料,应用时间域的最大熵谱反褶积算法,得到了11个（全部）PASSCAL（地壳与岩石层的地震台阵研究计划）台站的接收函数．利用时间域广义线性反演的Jumping（跳动）算法,引入模型光滑度约束,并将合成地震图的Kennett完全算法及微分地震图的Randall快速算法用于接收函数的正演计算,由台站接收函数获得了各台站下方的一维S波速度分布．反演结果表明,青藏高原的Moho界面在班公-怒江缝合带附近存在明显的深度错断;在日喀则、拉萨、桑雄和安多等地的地壳内部,可以连续观测到三个显著的速度界面h1,h2和h3,其中h1和h2可以连续追踪到温泉、二道沟和不冻泉等地,而h3在班公-怒江缝合带以北消失;在日喀则、拉萨、桑雄、安多、二道沟和不冻泉等地有壳内低速层．关于班公-怒江缝合带附近Moho界面的错断现象,一个可能的解释是班公-怒江缝合带是印度地壳向欧亚下地壳挤入的前沿． 相似文献
10.
11.
12.
使用相匹配滤波技术提取基阶面波振型的相位,然后再利用时变滤波进行 振幅谱校正,以此分离出基阶振型面波信号。对合成面波理论地震图的处理结果表 明,该方法可精确地从含高阶（最高阶数为５）的综合地震图中分离出基阶面波,与利 用同一模型合成的基阶面波相比,分离出的基阶面波的振幅谱畸变很小,根据中国数 字化地震台网（ＣＤＳＮ）记录的长周期面波资料,对与震源位于同一大圆弧上的北京。 海拉尔台观测的瑞雷波进行处理,分离出基阶振型信号,然后利用两台的基阶振型信 号进行反褶积得到台间格林函数,由格林函数计算得到台间周期为１３～１４２ ｓ的面 波衰减系数和群速度,所得结果也较为平稳。 相似文献
13.
14.
Ricker‐compliant deconvolution spikes at the center lobe of the Ricker wavelet. It enables deconvolution to preserve and enhance seismogram polarities. Expressing the phase spectrum as a function of lag, it works by suppressing the phase at small lags. A by‐product of this decon is a pseudo‐unitary (very clean) debubble filter where bubbles are lifted off the data while onset waveforms (usually Ricker) are untouched. 相似文献
15.
Bohuslav R??ek Jaroslava Plomerová Vladislav Babu?ka 《Studia Geophysica et Geodaetica》2012,56(1):107-140
Joint inversion of teleseismic P-waveforms and local group velocities of surface waves retrieved from ambient seismic noise
has been performed to model velocity structure of the crust and uppermost mantle of the Bohemian Massif. We analysed P-waveforms
of 381 teleseismic earthquakes recorded at 54 broadband seismic stations located on the territory of the Czech Republic and
in its close surroundings. Group velocities of Rayleigh and Love surface waves were obtained by cross-correlating long-term
recordings of seismic noise. The basis for waveform inversion is the well-known methodology of P-to-S receiver functions constructed
from converted phases. Due to instabilities in direct inversion of receiver functions caused by the necessity of applying
deconvolution, we propose an alternative formulation to fit observed and calculated radial components of P waveforms. The
joint inversion is transformed into a search for the minimum of the cost function defined as a weighted sum of waveform and
group velocity misfits. With the use of the robust stochastic optimizer (Differential Evolution Algorithm), neither derivatives
nor a starting model are needed. The task was solved for 1D layered isotropic models of the crust and the uppermost mantle.
We have performed a sequence of inversions with models containing one, two, three and four layers above a half-space. By using
statistical criteria (F-test) we were able to select the simplest velocity models satisfying data and representing local geological
structures. Complex crustal models are typical for stations located close to boundaries of major tectonic units. The relatively
low average P to S wave-velocity ratio is in agreement with the generally accepted view that the BM crust is predominantly
felsic. 相似文献
16.
《Geofísica Internacional》2014,53(2):211-220
We apply a single-step, finite-fault analysis procedure to derive a coseismic slip model for the large MW 7.4 Ometepec-Pinotepa Nacional, Mexico earthquake of 20 March 2012, using teleseismic P waveforms recorded by the Global Seismographic Network. The inversion is conducted in near-realtime using source parameters available from the USGS/NEIC and the Global Centroid Moment Tensor (gCMT) project. The fault orientation and slip angle are obtained from the gCMT mechanism assuming that the fault coincides with the shallow-dipping nodal plane. The fault dimensions and maximum rise time are based on the magnitude reported for the event. Teleseismic data from the USGS/NEIC Continuous Waveform Buffer database are used in the inversion with record start times set to the P-wave arrivals used to compute the earthquake hypocenter. The inversion is stabilized by requiring a smooth transition of slip across the fault while minimizing the seismic moment. These constraints are applied using a smoothing weight that is estimated from the inverse problem, allowing the recovery of the least-complicated rupture history in a single step. Inversion of the deconvolved, ground-displacement waveforms reveals a simple, circular rupture similar in extent to the source identified by the USGS/NEIC using body-and surface-wave data, indicating that the teleseismic P waves can provide a first-order source model for the event in near-realtime. Additional inversions conducted using velocity records identify a more-detailed rupture model characterized by an elliptical 2500 km2 source region extending updip and downdip from the hypocenter. This elliptical source preserves the orientation and overall dimensions of a dual-source slip model obtained recently by other investigators using local strong motions and global seismic waveforms. The results indicate that velocity waveforms could provide additional details of the earthquake rupture in near-realtime, finite-fault inversions using teleseismic P waves. 相似文献
17.
鄂尔多斯地块东南缘地带Moho深度变化特征研究 总被引：7，自引：2，他引：5
鄂尔多斯地块东南缘是主要的历史强震活跃区,曾经多次发生6级或以上的强烈地震,其边缘边界具有较强的地震活动性.本文利用该区域内分布的固定台站数据记录的大量远震体波波形资料,应用频率域反褶积方法提取远震P波接收函数,由H-κ方法测定了各台站下方的Moho深度和Vp/Vs值.研究结果表明:鄂尔多斯地块东南缘的Vp/Vs值介于1.6~1.9之间.东缘的Moho深度介于33.4~45 km之间,太原断陷盆地附近的Moho深度较浅,最浅处为33.4 km;东部北段的延怀盆地、蔚县盆地、阳原盆地和南段的临汾盆地附近Moho深度变化不大,平均深度为40 km.而在东缘东侧,因存在着山西断陷带,导致块体边缘的Moho深度要小于块体内部的Moho深度.块体南缘的Moho深度介于31.0~53.1 km之间,自东段向西段Moho深度逐渐变大,从渭河盆地附近的31.0 km增厚至秦岭造山带地段的53.1 km.总之,鄂尔多斯地块东南缘地带的Moho深度和Vp/Vs值分布具有明显的分块特征,块体内部结构比较稳定,东缘东段地壳结构相对一致,东缘东侧与西侧地壳深度具有明显的差异性,从山西断陷以东向西地壳厚度逐渐增厚,很好地对应了其地质构造特点. 相似文献
18.
A method is proposed to obviate the shortcomings of conventional deconvolution approaches applied to vibroseis data. The vibroseis wavelet reduces the time domain resolution of the earth's impulse response by restricting its passband. The spectrum of the wavelet is assumed to be a “low quefrency”phenomenon, and hence it can be estimated by low cut cepstral filtering. The wavelet's amplitude spectrum can then be removed by spectral division. By using an approach which is consistent with the principle of maximum entropy, the undetermined portions of the seismogram's Fourier transform can be filled in by autoregressive prediction. The process of initially deconvolving in a restricted passband reduces the enhancement of noise contaminated parts of the spectrum, and the spectral extension scheme increases the time domain resolution of the process. 相似文献
19.
考虑到面波频散对介质S波速度、接收函数对界面深度的各自敏感性优势,综合利用面波和接收函数资料实现联合反演,求取滇西地区壳幔速度结构. 本文利用适配滤波频时分析技术处理覆盖滇西地区的长周期面波资料,获得105～1050s周期范围内的面波群速度频散,进而利用分格反演方法提取研究区内1°×1°网格纯路径频散;基于滇西地区宽频带三分量远震记录,经反褶积后得到台站下方的远震P波接收函数. 联立面波纯路径频散信息和接收函数资料建立系统方程,利用阻尼最小二乘法实现联合反演,从而获得滇西地区壳幔S波速度结构. 结果表明,滇西地区以红河断裂为界,东西两侧壳幔结构存在明显差异,断裂西侧约20km深度处存在一厚度为10km左右的低速层,而东侧并不明显;滇缅泰块体上的畹町、沧源一带属于上地幔低速区,而另一个地幔低速区则位于滇中块体上的康滇古隆起上,两处地幔低速区与大地高热流分布、强震活动具有较好的对应关系. 相似文献
20.
S. M. DEREGOWSKI 《Geophysical Prospecting》1978,26(1):252-290
In a previous paper the author showed how, by computing an inverse filter in the frequency domain, an automatic compromise could be made between the conflicting requirements to spike a wavelet and to keep the attendant noise amplification within bounds. This paper extends the technique to take account of errors in the estimated shape of the wavelet defined to the deconvolution process. The drastic effects which such errors can have if they are ignored are demonstrated. A novel form of filter–called the “self-matching filter”–is defined which allows the user to limit not only the noise amplification but also the sensitivity of the filter to random uncertainties in the estimated wavelet. This is achieved by whitening the spectrum only within automatically selected pass bands whilst suppressing other noise-dominated or uncertainly defined frequency components. Conventional Wiener filtering is shown to be a special case of this more general filter, namely one in which the wavelet uncertainty is completely ignored. The type of phase spectrum which the output pulse should be designed to possess (e.g. zero phase or minimum phase) is briefly discussed. 相似文献