接收函数方法软件

1 接收函数研究概况：

转换波的地壳测深方法自70年代被介绍到我国，并曾经成为除人工地展测深以外研究地壳和上地幔结构的重要方法(邵学钟和张家茹，1978；刘启元和邵学钟，1985；张家茹和邵学钟，1994)。它利用远震p波入射到台站下方时在介质间断面上产生的ps转换震相与透射p 波的相对到时差研究地下介质间断面的深度分布。

转换波测深的一些主要思想在进一步的接收函数研究中得到了极大发展。

langston (1979)利用远震p波波形的这个特点提出了等效震源假定，并提出了从长周期远震体波波形数据中分离接收台站下地球介质对入射p波的脉冲响应(即接收函数)的方法。 owens et al. (1984) 将接收函数的方法进一步扩展到宽频带记录的情况，并发展了相应的远震体波接收函数的线性波形反演方法。利用远震接收函数反演方法，人们可以根据宽频带远震p波的波形数据获得台站下方岩石圈的s波速度结构。

其理论和方法也获得了不断的改进和发展.其中，randall(1989)提出了计算微分地震图的高效率方法，ammon et al. (1990) 针对接收函数反演的非唯一性提出了保留接收函数径向分量绝对振幅的接收函数分离方法。刘启元等(1996)提出了从宽频带地震台阵资料获取三分量接收函数的方法并实现了基于tarantola矢量反演理论的接收函数非线性反演方法，接收函数的反演方法在国内外己获得了日益广泛的实际应用。

在研究基于一维介质假设的接收函数及其反演方法的同时，针对接收函数切向分量上地震波能量的研究也在同时进行。主要是研究介质的非均匀性，各向异性。

zandt & ammon (1995)以及zhu & kanamori(2000) 利用接收函数ps转换震相和多次震相研究了地壳厚度和地壳平均poisson比结构。由于地壳poisson比结构包含着比地壳p波和s波速度结构更多的地壳介质成分和动力学演化信息，这个方法得到了越来越多的应用

近年来，接收函数方法的一个重要发展方向就是将地震勘探中广为应用的地震偏移技术移植到天然地震台阵观测数据的解释，用以研究地壳和上地慢速度间断面的横向变化。

2 接收函数的提取方法：

所谓远震体波接收函数即地震台站下方(接收区)介质对以近垂直角度入射到接收区的地震体波的脉冲响应。

传统的接收函数，也可以称为ｐ波接收函数，主要利用坐标旋转和反卷积的方法从远震ｐ波中分离ｐｓ震相，并且目前已经成为一种常规的数据分析工具来研究区域性的壳幔结构。

首先，可以利用后方位角把原始的ｚｎｅ三分量地震记录旋转到ｚｒｔ（垂向，径向和切向）坐标系，后方位角根据震中和台站的位置由理论计算得到。

langston(1979)提出用震源等效化方法来消除有效震源时间函数对远震p波波形的影响，得到了所谓的接收函数。他认为从一系列水平分层或倾斜分层介质底部入射的平面p波产生的地表位移响应在时间域可表示为：

?dv(t)?i(t)*s(t)*ev(t)??dr(t)?i(t)*s(t)*er(t) （7）

?d(t)?i(t)*s(t)*e(t)t?t

其中，s(t)代表入射平面波的有效震源时间函数，i(t)代表仪器的脉冲响应，ev(t)、er(t)、et(t)分别代表接收区介质结构脉冲响应的垂直分量、径向分量和切向分量。

频率域表示为：

?dv(?)?i(?)s(?)ev(?)??dr(?)?i(?)s(?)er(?)

?d(?)?i(?)s(?)e(?)t?t

理论计算与实际观测表明，近垂直入射的远震p波波形的垂直分量主要由近似脉冲的直达波构成，尾随波列能量较弱，可忽略不计，于是可做如下近似：

ev(t)??(t),ev(?)?1.

于是，地震p波记录的垂直分量可以近似为：

dv(t)?i(t)*s(t)

将频率域稍作变换，径向和切向的介质结构相应可表示为：

er(?)?

et(?)?dr(?)d(?)?ri(?)s(?)dv(?)dt(?)d(?)?t................(1)i(?)s(?)dv(?) 再将er(?)和et(?)分别反变换到时间域，即可得到介质结构的径向接收函数er(t)和切向接收函数et(t)。

由此可看到，接收函数的提取实际是一个反褶积运算的问题，按照计算域不同可大致分为两种：频率域反褶积和时间域反褶积。各种反褶积方法各有优点，没有一种方法具有绝对的优势；在数据质量较好时，它们之间的差异很小。

2.1 频率域接收函数提取

目前最常用的是频率域的水准量反褶积方法。在处理实际资料时，（1）式中垂直分量的频谱有时会有零点或者很小的值出现，从而导致运算的不稳定。通常引入“水准量”来压制靠近零点值的频谱（clayton & wiggins, 1976）:

*dr(?)dr(?)?dv(?)er(?)??g(?)dv(?)?(?)

et(?)?dt(?)dt(?)?d(?)?g(?)dv(?)?(?)*v

**其中，c为水准量常数，0<??<1，?(?)?max?dv(?)dv(?),c?max?dv(?)dv(?)}}，由经验选择。g(?)为频率高斯滤波函数，表达式为g(?)?a0exp(??2/4a2)，a0为归一化常数，a为高斯滤波器的宽度因子，a越大滤波后波形中含有的高频成分越多，实际中a的选取要根据垂向波形的脉冲宽度。

频率域内水准量反褶积方法提取接收函数简单便捷，但固有缺点就在于引入水准量的同时也降低了精度。

2.2 时间域接收函数提取

时间域内提取接收函数的方法多种多样，比如wiener滤波反褶积（taylor et al., 1982）、最大熵谱反褶积（tselentis, 1990）、最大或然性反褶积（shumway et al., 1982; 刘启元等，1996）、同步时间域反褶积（gurrola et al., 1995），以及迭代反褶积方法（kikuchi and kanamori, 1982; ligorria and ammon, 1999）等。

wiener 滤波反褶积以远震p波波形的垂直分量作为输入，以接收函数作为滤波因子，以远震p波波形的水平分量（径向和切向）作为期望输出，通过远震p波波形垂直分量与接收函数的褶积得到wiener 滤波器的实际输出，以期望输出与实际输出的均方误差取极小，作为求取接收函数的准则。wiener 滤波反褶积可归结为toeplitz方程的求解； levinson递推算法大大提高了wiener 滤波的计算效率； toeplitz方程的非奇异性保证了wiener 滤波反褶积的稳定性。

也可以利用入射角再把ｚｒｔ坐标系选转到ｌｑｔ射线坐标系坐标系：

以脉冲反褶积为例，在时间域中利用最小化最小二乘意义上的观测的l分量和期望的具有归一化幅度的δ脉冲函数之间的差异的方法获得一个逆滤波器，然后把这个逆滤波器分别与ｌｑｔ三分量进行卷积，分别得到期望的零相位脉冲（ｌ分量）

、ｑ分量的接收函数和ｔ

分量接收函数．最后，反卷积后得到的ｌｑｔ三分量根据l分量包含的脉冲函数的最大幅度进行归一化处理。q分量的接收函数包含远震ｐ波穿透台站下方的地震间断面产生的ｐｓ转换波以及源自间断面和地表之间多次反射的多次波如ｐｐｐｓ和ｐｐｓｓ，通常所说的ｐ波接收函数指的是q分量的接收函数（图１）。t分量的接收函数则显示了台站下方倾斜和各向异性结构．

使用q分量代替r分量的优点主要体现在两个方面：一个方面是q分量上直达ｐ波能量