共查询到20条相似文献，搜索用时 171 毫秒
1.
用Wiener滤波方法提取台站接收函数 总被引：10，自引：1，他引：10
本文提出了一种在时间域用Wiener滤波方法提取台站接收函数的方法，用远震P波波形的垂直分量为输入，接收函数作为滤波因子，远震P波波形的径向和切向分量作为期望输出，通过期望输出与实际输出的均方误差达极小，来提取接收函数。接收函数的计算可归结为Toeplitz方程的求解，可以采用Levinson递推算法。Toeplitz方程的非奇异性保证了Wiener滤波反褶积方法的稳定性。合成地震图与观测地震图的检验表明，用Wiener滤波方法测定台站接收函数是一种有效的时间域反褶积方法。 相似文献
2.
Receiver Functions from Autoregressive Deconvolution 总被引：4，自引：0，他引：4
Qingju Wu Yonghua Li Ruiqing Zhang Rongsheng Zeng 《Pure and Applied Geophysics》2007,164(11):2175-2192
Summary Receiver functions can be estimated by minimizing the square errors of Wiener filter in time-domain or spectrum division in
frequency domain. To avoid the direct calculation of auto-correlation and cross-correlation coefficients in Toeplitz equation
or of auto-spectrum and cross-spectrum in spectrum division equation as well as empirically choosing a damping parameter,
autoregressive deconvolution is presented to isolate receiver function from three-component teleseismic P waveforms. The vertical
component of teleseismic P waveform is modeled by an autoregressive model, which can be forward and backward, predicted respectively.
The optimum length of the autoregressive model is determined by the Akaike criterion. By minimizing the square errors of forward
and backward predicting filters, autoregressive filter coefficients can be recursively solved, and receiver function is also
estimated in the similar procedure. Both synthetic and real data tests show that autoregressive deconvolution is an effective
method to isolate receiver function from teleseismic P waveforms in time-domain. 相似文献
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.
提出一种在时间域采用最大熵谱反褶积提取台站接收函数的方法，以最大熵作为自相关函数和互相关函数的递推准则，利用Toeplitz方程及Levinson递推算法，得到预测误差滤波系数的递推公式，从而计算台站接收函数．外推运算过程中，反射系数总是小于1，保证了最大熵谱反褶积的稳定性．时窗外数据熵极大提高了接收函数的分辨率．合成地震图与实测地震图的检验表明, 最大熵谱反褶积是一种在时间域测定台站接收函数的有效方法． 相似文献
6.
用接收函数方法研究上海地震台阵下地壳结构 总被引：1，自引：1，他引：0
利用上海地震台阵16个台站记录的远震资料,采用接收函数线性反演方法,对台阵下的地壳速度结构进行研究,获得了研究区域内地壳厚度和地壳速度的分布特征。研究结果表明,研究区域Moho面深度约为33±2 km,Moho面深度基本不变,地幔顶部S波速度约4.4 km/s,地壳内没有发现明显的低速层。 相似文献
7.
8.
Mapping crustal S—wave velocity structure with SV—component receiver function method 总被引：1，自引：0，他引：1
In this article,we analyze the characters of SV-component receiver function of teleseismic body waves and its advantages in mapping the S-wave velocity structure of crust in detail.Similar to radial receiver function,SV-component receiver function can be obtained by directly deconvolving the P-component from the SV-component of teleseismic recordings.Our analyses indicate that the change of amplitude of SV-component receiver function against the change of epicentral distance is less than that of radial receiver function.Moreover,the waveform of SV-component receiver function is simpler than the radial receiver function and gives prominence to the PS converted phases that are the most sensitive to the shear wave velocity structure in the inversion.The synthetic tests show that the convergence of SV-component receiver function inversion is faster than tnat of the radial receiver function inversion.As an example,we investigate the S-wave velocity structure beneath HIA sta-tion by using the SV-component receiver function inversion method. 相似文献
9.
详细讨论了远震体波SV分量接收函数的特点及其在反演地壳S波速度结构中的优势．与径向接收函数类似，SV分量接收函数可通过对远震体波的SV分量直接反褶积P分量获得．研究分析表明:与径向接收函数相比，SV分量接收函数的振幅随震中距的变化更加稳定，波形简单且突出了对结构最敏感的PS转换波信息．理论数值实验显示:在反演地壳S波速度结构时，SV分量接收函数比径向接收函数具有更好的收敛性．作为实例，利用SV分量接收函数反演方法反演了海拉尔台下的S波速度结构． 相似文献
10.
11.
利用新疆和田地震台阵3 km孔径范围内架设的9个子台(包括1套宽频带和9套短周期地震仪)记录的3年远震波形数据,对比研究和分析了短周期地震仪接收函数的稳定性和可靠性.通过比较和分析短周期和长周期地震仪获取的接收函数波形,结果发现:(1)短周期地震仪记录与宽频带地震仪记录得到的接收函数有很好的一致性,且具有较高的稳定性,无论α取值为1.5还是2.5,短周期地震仪接收函数与宽频带地震仪接收函数都具有较强的线性相关性(相关系数0.9),但Ps震相均存在小幅的振幅差(约20%);(2)采用接收函数震相到时信息的方法(如H-κ叠加搜索),短周期地震仪可以代替宽频带地震仪;(3)由于短周期地震仪缺乏0.155 Hz以下的低频信号和在1 Hz以下频段非线性的振幅响应,仅仅采用短周期地震仪接收函数波形反演台站下方S波速度结构,获得下地壳到上地幔顶部的速度偏差较大(约0.3 km/s),可能会造成错误解释(如下地壳低速层),因此需要和其它对波速值敏感的数据(如面波频散)进行联合分析. 相似文献
12.
One of the problems in signal processing is estimating the impulse response function of an unknown system. The well-known Wiener filter theory has been a powerful method in attacking this problem. In comparison, the use of stochastic approximation method as an adaptive signal processor is relatively new. This adaptive scheme can often be described by a recursive equation in which the estimated impulse response parameters are adjusted according to the gradient of a predetermined error function. This paper illustrates by means of simple examples the application of stochastic approximation method as a single-channel adaptive processor. Under some conditions the expected value of its weight sequence converges to the corresponding Wiener optimum filter when the least-mean-square error criterion is used. 相似文献
13.
鄂尔多斯地块东南缘地带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值分布具有明显的分块特征,块体内部结构比较稳定,东缘东段地壳结构相对一致,东缘东侧与西侧地壳深度具有明显的差异性,从山西断陷以东向西地壳厚度逐渐增厚,很好地对应了其地质构造特点. 相似文献
14.
利用天然地震记录直达P波后续20s内的波形信息,研究了地壳、土地慢速度结构反演的方法、唯一性、精度及应注意的问题,并将求解有条件极值的惩罚数法和求解无条件极值的单纯形法引入到波形反演.数值计算证明,波形反演方法对记录的误差有压制作用,在深源、远震、各向同性水平层状近似合理的前提下,可用一个地震事件在射线平面内的二分量记录反演台站区地壳、上地幔P、S波速度结构.在新疆天山地区转换波流动台站观测中,选择出较好的记录进行波形反演,给出了该转换被测深剖面上6个台站下方的地壳、上地慢P、S波速度结构.通过塔里木盆地内一测点上波形反演和地震勘探结果的对比证明,波形反演方法具有满意的精度. 相似文献
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.
分析了相邻井下、地表地震计记录的背景噪声及近震、远震波形特征,以及相关的振幅谱及近震、远震接收函数波形特征。结果表明,井下地震计记录的背景噪声比地表地震计小1个数量级,其近震、远震波形记录较清晰、简洁;井下、地表地震计记录的远震振幅谱和径向接收函数具有较好的一致性,但背景噪声、近震振幅谱、近震接收函数、远震切向接收函数波形存在一定差异;由地表地震计得到的近震震级比井下地震计的大0.12,地表地震计到井下地震计之间的P波速度约为2 km/s。 相似文献
17.
冰盖厚度是研究南极冰盖质量、建立冰盖动力学模型的基本参数,对于冰川均衡调整、冰盖物质平衡及全球气候变化研究具有重要意义.基于地震学的远震接收函数和H-Kappa格网搜索方法可以用于地震台站下方冰盖厚度的可靠探测,不仅能与冰雷达获得的冰盖厚度进行独立对比,还可以与冰雷达方法相互补充,进一步填补南极大陆冰盖厚度探测空白区.本文利用布设于南极大陆冰盖上方的流动地震台阵记录到的远震波形数据,基于接收函数方法对台阵下方的冰盖厚度进行了研究.结果显示:基于远震接收函数方法的冰盖厚度与Bedmap2冰厚格网模型相比,二者差别大多在200 m以内;少数台站差值达到600 m左右,这一差别可能与Bedmap2测线分布空区、冰雷达测深不确定性以及冰盖内部复杂波速结构等因素有关.本文研究结果表明:利用南极大陆冰盖上方的流动地震台阵,基于远震接收函数方法可以获得比较可靠的南极冰盖厚度,为独立验证冰雷达的探测结果并弥补冰雷达探测空白区提供了有效方法.同时,部分台站接收函数波形的复杂性可能暗示了南极大陆数千米厚的冰盖内部结构不是均一的,仍然存在比较复杂的内部结构变化.因此,有必要进一步利用包括接收函数波形拟合、地震面波反演等方法对南极大陆冰盖厚度及其内部精细结构进行更为深入的研究. 相似文献
18.
Twenty broadband seismographs were deployed along Hongyuan, Sichuan to Wuwei, Gansu. 81 teleseismic events were recorded in one year. We computed receiver functions from teleseismic waveform data and obtained S wave velocity structure beneath each station along the profile by using receiver function inversion method. The results revealed that the crustal structure is very complex and crustal average S wave velocity is to be on the low side. Low velocity structure generally exists in the depth range of 10～40 km in the crust between Aba arc fault and northern edge fault of Qinling earth＇s axis and it is a tectonic feature of complex geological process such as ancient A＇nyemaqen Tethys ocean from closing and side colliding to subducted plate exhumed or thrust rock slice lifted. The Moho is about 50 km depth along the profile and is slightly deeper in the south than in the north. 相似文献
19.
2004年东乌珠穆沁旗地震震源参数研究 总被引：5，自引：0，他引：5
利用波形反演方法, 研究并确定了2004年3月24日发生在内蒙古东乌珠穆沁旗与西乌珠穆沁旗交界处地震的震源机制及深度。 近震波形对震源机制有很好的约束, 而远震体波对于震源深度分辨率很高, 因此综合利用了近震和远震波形数据。 对于近震宽频带波形记录, 用CAP方法反演震源机制解。 利用远震直达P波、 pP波和sP波的到时能较好的确定震源深度。 综合利用近震宽频带波形和远震体波波形记录, 反演得出2004年03月24日地震为逆冲型, 震级为MW5.3。 其两个节面分别为: 节面I的走向、 倾角、 滑动角分别为147°、 22°、 87°, 节面Ⅱ的走向、 倾角、 滑动角分别为 330°、 68°、 91°; 震源深度为(12±2) km, 是一个典型的浅源上地壳地震。 与前震序列对比, 节面Ⅰ是可能的发震构造。 相似文献
20.
To determine the crustal structure in central Tibet, we used teleseismic waveform data recorded by 18 stations in the INDEPTH-Ⅲ seismic array across the central Tibet from the central Lhasa terrane to the central Qiangtang terrane. The S-wave velocity structures beneath stations are determined by inverting the stacked radial receiver function using the GA method. The first order features in the receiver function are modeled. Our results show that the Moho in Qiangtang is about 8 km shallower than that in Lhasa terrane along the INDEPTH-Ⅲ profile. It maybe suggests the northward subduction of the Lhasa mantle lid beneath the Qiangtang terrane is affected by the India-Asia collision. We conclude that there exist low velocity zone in the middle crust across the northern Lhasa and Qiangtang terrane, which can be related to the high temperature upper mantle beneath that. 相似文献