接收函数波动方程叠后偏移方法.doc

接收函数波动方程叠后偏移方法介绍 陈凌，张耀阳 中国科学院地质与地球物理研究所 接收函数波动方程叠后偏移方法是我们借鉴勘探地震学中发展成熟的反射波偏移成像技术，2005年以来发展的一种新的天然地震接收函数偏移成像方法（Chen et al., 2005a; 2005b）。该方法包括两个基本步骤：时间域CCP 叠加和深度域波场反向延拓。与反射地震学中处理反射波资料的CMP 叠加类似，CCP 叠加用以提高资料的信噪比。波场反向延拓则是一个将产生Ps 转换波（图1d）的转换波源偏移至其真实位置的过程。在CCP 叠加中，用一个1D 参考模型将来自单个地震事件的接收函数按垂直入射情况做时间校正后，根据其转换点位置进行组合叠加。所得到的CCP 叠加道集可以近似看作是零偏移距资料，即转换波从速度间断面的转换点处垂直向上传播至地表而被记录到的波场。对于复杂结构，只采用CCP 叠加将会给成像结果带来假像，而偏移对于转换界面的正确成像是必不可少的（Ryberg and Weber, 2000; Poppeliersand Pavlis, 2003a）。我们采用的波场偏移原理与反射地震学中普遍采用的爆炸反射面原理（Claerbout, 1985; Sheriff and Geldart, 1995）类似。在爆炸反射面模型中：每一个反射界面都被当作是布满了密度正比于界面反射强度的爆炸源。将这些爆炸源同时在时间t = 0 激发，地震波以实际速度的二分之一向上传播以满足双向传播走时（图1a, b）。这样在地表上接收的波场可以模拟由CMP 叠加得到的零偏移距（零炮点-接收点距离）道集。倒转上述过程，即将地表观测的零偏移距波场反向传播至t = 0 时刻，就可以得到一个使所有反射体正确归位的深度偏移成像结果。在接收函数偏移中，通过CCP 叠加产生的零偏移距接收函数道集可以用来模拟地表记录到的将所有地下转换波源同时在t = 0 时刻激发得到的波场，而源的强度正比于Ps 的转换系数（图1c）。这种模型称为爆炸转换波源模型，是由Pavlis（2003）首先提出并应用于接收函数成像的。将接收函数的CCP 叠加道集反向延拓至整个空间并在t = 0 时刻成像，就能使所有转换波偏移至产生转换波的相应地下速度间断面的真实位置处。Ryberg和Weber（2000）提出了同样的接收函数偏移方案。我们所提出的???移方法与Ryberg和Weber 方法的本质区别在于采用的波场延拓技术不同。在他们的研究中采用基于射线的Kirchhoff 高频近似方法，而我们则应用波动方程波场传播技术。 1．CCP 叠加 传统的接收函数CCP 叠加过程包括接收函数的走时校正和组合叠加（Dueker andSheehan, 1997, Zhu, 2000; Ai et al., 2003）。这种数据处理过程与反射地震学中的CMP 叠加类似，但两者之间仍存在着显著差别。首先，在接收函数处理过程中，时间校正是为了补偿由不同震中距入射P 波造成的P 波转换成S 波的走时差。因此，通常是以一个1D 参考模型计算得到的参考震中距的走时作为标准，对所有接收函数进行时间校正。参考震中距一般选取为所有资料的平均震中距，以使时间校正量达到最小。由于接收函数的典型震中距是30°到90°，参考震中距也应落在该范围内，因此必然对应于P 波倾斜入射的情况（入射平面波的水平慢度p ≠ 0）。另一方面，在反射资料分析中，反射波的走时均按垂直入射（p = 0）进行校正。如此得到的CMP 叠加道集才能近似作为叠后偏移所需要的零偏移距资料。为了应用叠后偏移对接收函数进行处理，定义相对于假设零偏移距情况（震中距为180°，p = 0）的Ps 转换波走时校正量为 图1反射地震学中获得零偏移距资料的探测系统(a)、叠后偏移采用的爆炸反射面模型(b) ；接收函数叠后深度偏移采用的爆炸转换波源模型(c)，以及Ps 转换波与地表多次反射波示意图(d)。 (1) 其中接收函数中Ps 转换波走时TPS 是基于展平地球模型中平面波入射假定得到的直达P 波与转换S 波的走时差（Gurrola et al., 1994; Dueker and Sheehan, 1997） (2) 上式中D, VP, VS 分别为Ps 转换深度，随深度变化的P 波速度和S 波速度。在走时校正之后，根据各个深度处的转换点位置将所有接收函数重新组合到各自的共转换点单元中，以达到最佳聚焦的目的。对共转换点单元中的接收函数分别进行叠加而构成的CCP 道集近似相当于P 波垂直入射的观测资料，因而将作为下一步叠后深度偏移的初始波场。 2．基于波动方程的波场反向延拓 为了由接收函数的 CCP 叠加道集对产生Ps 转换波的地下间断面真实成像，需要采用一个适当的波场