地震接收函数.docVIP

地震接收函数 篇一：接收函数法基本原理 深部速度结构反演的接收函数法 3.1远震P波波形接收函数的求取方法 接收函数法是利用远震P波波形的单台记录来反演台站下方一维S波速度结构的波形反演方法。远震P波波形含有关于震源时间函数、源区介质结构、上地幔传播路径以及接收区介质结构的丰富信息。远震P波波形与这些影响机制的关系可表示成： D(t)?S(t)*MS(t)*MRay(t)*MR(t)*I(t) （6） 其中：D(t)为所记录的远震P波波形数据； S(t)为震源时间函数； MS(t)为近源介质结 I(t)构响应；MR(t)为台站下方接收介质的响应；MRay(t)为P波在地幔中传播的透射响应； 为仪器响应。 在以上因素中，除了仪器响应外，其它因素都是难以一一加以确定的。而只有台站下方介质的响应才是我们所感兴趣的、可用来反演台站下方地壳、上地幔速度结构的波形信息。因此要有一种方法将接收介质的响应从整个P波波形中分离出来，而接收函数法就是这样一种行之有效的方法。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)分别代表介质结构脉冲响应的垂直分量、径向分量和切向分量。 对于许许多多波形简单的远震事件的观测表明，深源远震地表位移的垂直分量表现为尖脉冲的时间函数与仪器响应的褶积，紧随其后的续至震相非常小（Burdick and Helmberger,1974）。理论计算也表明，即使地壳内存在角度适中的强速度界面，陡角度入射P波所产生的转换波及地壳内部的鸣震震相的垂直分量也是非常小的（Burdick and Langston,1977）。因此，可以认为介质结构响应的垂直分量近似为Dirac函数，即： EV(t)??(t) （8） 显然在（8）的假设条件下，地表位移的垂直分量可作为与接收介质响应无关的远震P波波形的影响因素， 也就是说，地表位移的垂直分量可以近似为仪器响应和有效震源时间函数 的褶积： I(t)*S(t)?DV(t)（9） 如果三分量地震仪的脉冲响应都一致的话，那么用DV(t)对DR(t)、DT(t)分别作反褶积处理就可以得到ER(t)、ET(t)。反褶积在频率域可表示成： ?ER(?)?????ET(?)???DR(?)D(?)?RI(?)S(?)DV(?)（10） DT(?)DT(?)?I(?)S(?)DV(?) 将ER(?)、ET(?)分别反变换回时间域，就可得到介质结构响应的径向分量ER(t)和切向分量ET(t)，也就是所谓的径向接收函数和切向接收函数。由于实际地震资料是有限带宽的， 且包含随机噪声，直接用（10）在频率域作除法运算往往是不稳定的。为了确保频率域反褶积的稳定性，常采用Helmnerger and Wiggins(1971)提出的频率域反褶积稳定算法。 为了克服频率域反褶积的固有缺点，需要发展数据长度和滤波因子长度灵活可变的时间域反褶积方法，以获得稳定而精确的接收函数。Wiener 滤波反褶积以远震P波波形的垂直分量作为输入，以接收函数作为滤波因子，以远震P波波形的水平分量（径向和切向）作为期望输出， 通过远震P波波形垂直分量与接收函数的褶积得到Wiener 滤波器的实际输出，以期望输出与实际输出的均方误差取极小，作为求取接收函数的准则。Wiener 滤波反褶积可归结为Toeplitz方程的求解； Levinson递推算法大大提高了Wiener 滤波的计算效率； Toeplitz方程的非奇异性保证了Wiener 滤波反褶积的稳定性。此外还可引入镶窗技术来压制噪声，提高接收函数的精度。另外还可以用最大熵谱反褶积法（Tselentis,1990）来求取接收函数。 3.2模型接收函数的计算方法 为了反演速度结构，需要将从实际资料得到的接收函数与理论模型的接收函数相拟合，从而确定速度参数随深度的变化。而要得到理论接收函数，就需要计算某一模型下的理论地震图。而为了在反演中求解参数变化对理论结果的影响，就需要求解微分地震图。这是一个十分复杂的工作。对此Kennett（1979,1980）已有过详细的描述和推导，这里就不再贅述。需要说明的是，由于理论地震图是在频率域合成的，用地表位移的径向分量直接除以垂直分量，最后反变换回时间域，就可以得到接收函数理论地震图。这与实际