引 言

地震波(P，S，pP，sP等)在传播过程中，遇到间断面会发生反射或转换，所产生的新震相(P₄₁₀P，P₆₆₀P，S₄₁₀S，S₆₆₀S，S₄₁₀P，S₆₆₀P，s₄₁₀P等)为间断面(“410”，“660”等)的深入研究提供了地震学基础(Flanagan，Shearer，1998a； Deuss，2009； Schmerr，Thomas，2011)． 这些前驱震相在单条记录中的能量(幅度)往往很弱，常淹没于噪声之中； 但是通过叠加大量的观测波形数据，则可以压制噪声，更加有效地提取与间断面相关的有用信号(Rost，Thomas，2002； 臧绍先，周元泽，2002)．

很多研究者利用下行的转换震相在间断面的转换点深度确定该间断面的深度． 例如： Richards和Wicks(1990)使用S_dP转换波，利用非线性叠加方法对汤加地区下方的“670”间断面的深度和性质进行了深入研究； 谢彩霞等(2012)利用4次根倾斜叠加方法有效提取了离源下行的S_dP次生转换震相，进一步确认了汤加—斐济下方300 km附近速度间断面的存在； Li等(2008)通过叠加S--P转换波，计算了中国东北地区“660”间断面的深度，并探讨了俯冲带对“660”间断面的影响．

发生在俯冲带中的深震会在近源区速度间断面产生反射震相(p₄₁₀P，s₄₁₀P等)和转换震相(S₄₁₀P，S₆₆₀P等)(Flanagan，Shearer，1998b)． 长周期震相(PP，SS)前驱波方法对狭窄俯冲带间断面的横向分辨具有很大局限性(Flanagan，Shearer，1998a； Schmerr，Thomas，2011)． 接收函数方法是研究地球内部间断面的一种有效方法，可以通过P波与其上行转换波P_ds的到时差来确定间断面的深度(Langston，1979)，但其主要适用于三分量地震台站下方壳、 幔间断面的研究．

在深源地震(震源深度h＞300 km)分布密集的俯冲带(如汤加—斐济俯冲带)，使用深震震相sP在近源区速度间断面底界面反射的震相s₄₁₀P确定狭窄俯冲带间断面的横向变化具有很大的优势． 本文将通过倾斜叠加大量的观测波形数据，有效提取sP在近源区“410”速度间断面底界面反射的弱前驱波s₄₁₀P，利用sP震相与其前驱波s₄₁₀P震相的到时差计算汤加—斐济俯冲带“410”速度间断面的深度．

1.   台网资料

甘肃地震台网由甘肃数字测震台网和甘东南野外观测流动台阵构成，二者均处于青藏高原东北缘． 甘肃“十五”数字测震台网自2008年6月正式运行，由兰州、 高台、 安西、 嘉峪关和天水等5个有人值守的国家数字地震台和39个区域遥测数字地震台组成(冯建刚等，2012)，目前拥有包括邻省(宁夏、 陕西、 四川、 青海、 内蒙)在内的76个数字测震台站． 甘东南野外观测流动台阵是中国地震局地质研究所和中国地震局兰州地震研究所共同合作在甘东南地区勘址架设，共设有7条测线，150个流动台站，平均台间距约为10 km(图 1)． 所有台站均使用REF TEK-130型数据采集器和Guralp公司生产的CMG-3ESPC宽频带地震计． 甘东南野外观测流动台阵运行期间(2009年11月—2011年12月)积累了连续、 可靠、 高质量的地震观测波形数据，为深入研究地球内部物理结构提供了宝贵资料(秦满忠等，2015)．

图  1  甘肃数字测震台网和甘东南野外观测流动台阵分布

Figure  1.  Distribution of the Gansu Digital Seismic Network (blue triangles) and southeastern Gansu temporary observation array (black triangles)

下载: 全尺寸图片 幻灯片

本文选取了甘肃数字测震台网和甘东南野外观测流动台阵共同记录的2011年9月15日发生在汤加—斐济俯冲带(21.61°S，179.53°W，h=644 km，m_B=7.3)的深远地震数字观测波形数据，定位结果参考了USGS(2011)地震目录．

2.   数据处理

深震震相sP与其前驱波s₄₁₀P在地幔中的传播路径基本相同． sP是上行的S波在震中附近地表反射转换的P波； s₄₁₀P是上行的S波在震中附近“410”间断面底界面反射转换的P波(图 2)，该震相在初至震相P与sP之间到达(图 3)．

图  2  射线(sP， s₄₁₀P和P)传播示意图

Figure  2.  Ray paths of sP， s₄₁₀P and P

下载: 全尺寸图片 幻灯片

图  3  甘肃地震台网记录到的2011年9月 15日发生在汤加—斐济俯冲带的深远 地震的数字观测波形

s_MohoP， s_LABP和s₄₁₀P震相分别表示sP在 近源区莫霍面， 岩石圈和“410”的前驱波 The seismic phases s_MohoP， s_LABP and s₄₁₀P are the precursors from near-source underside reflection off Moho， lithosphere-asthenosphere boundary， and “410”, respectively

Figure  3.  Waveforms of m_B=7.3 event occurred in Tonga-Fiji on September 15， 2011， recorded by the Gansu Seismic Network

下载: 全尺寸图片 幻灯片

首先对所选取的汤加—斐济深远地震观测波形数据作重采样、 去倾斜、 去均值处理; 然后作带通滤波处理(0.1—0.01 Hz)，同时去掉记录畸形、 信噪比差的观测波形； 最终选取了甘肃数字测震台网和甘东南野外观测流动台阵共同记录到的170条垂直分量观测波形数据． 如图 3所示，汤加—斐济深远震在台网记录中的震中距为88°—102°．

通过倾斜叠加大量的观测波形数据，可以压制噪声，有效地提取与间断面相关的弱信号(Ritsema et al，1995； 臧绍先，周元泽，2002)． 在进行倾斜叠加的过程中，选取震中距为94°的台站(红崖山台)作为参考台，sP为参考震相，并以sP震相到时为零时刻(图 3)． 信号i(s₄₁₀P)在任意台s的到时为Δt_ir+ΔpΔD_rs，以红崖山台为参考台的倾斜叠加结果为

式中，Δt_ir为参考震中距记录的信号i(s₄₁₀P)与参考震相sP之间的走时差，Δp为信号i在任意台s的水平慢度与信号i(震中距为参考台震中距)水平慢度之差，ΔD_rs为任意台s与参考台的震中距之差， m为使用地震记录的台站数．

以sP为参考震相，使用倾斜叠加方法获得了在震中附近“410”间断面底界面反射的前驱波信号s₄₁₀P，其反射点位置为(21.25°S，179.99°W)． 该信号在倾斜叠加波形(图 3中红色线条)中成像清晰，与参考震相sP的观测走时差为137.43 s，使用TauP软件结合理论IASP91全球速度模型得到的s₄₁₀P与震相sP理论走时差为137.04 s； 通过使用基于CRUST1.0模型(Laske et al，2013)校正后的汤加—斐济地区速度模型IASP91_Tonga(图 4)，反算出“410”间断面底部反射点的深度约为398.5 km，“410”间断面抬升约11.5 km，这与俯冲带与近源区“410”间断面的相互作用有关．

图  4  全球速度模型IASP91以及修改后的汤加—斐济地区速度模型IASP91_Tonga

Figure  4.  IASP91 velocity model and modified velocity model (IASP91_Tonga) of the Tonga-Fiji zone in this paper

下载: 全尺寸图片 幻灯片

3.   讨论与结论

甘肃数字测震台网和甘东南野外观测流动台阵记录的汤加—斐济深震震相pP和sP均具有尖锐、 清晰和较大振幅的记录特征，sP记录振幅甚至超过了直达P波，这为我们使用前驱波(s_MohoP，s_LABP和s₄₁₀P等)对近源区间断面(莫霍面，岩石圈，“410”等)的深入研究提供了基础资料．

鉴于前驱波震相在单条记录中的能量(振幅)往往很弱，本文通过倾斜叠加大量的观测波形数据，获得了在近源区“410”间断面底界面反射的前驱波震相s₄₁₀P，并使用CRUST1.0模型校正后的汤加—斐济地区速度模型IASP91_Tonga，反算出“410”间断面底部反射点的深度约为398.5 km．

图 3中震相s_LABP和s_MohoP相对于参考震相sP的走时差分别为23.16 s和10.62 s. 我们同样对震源区的地壳及岩石圈结构作了修改(图 4)，使用TauP软件计算得到近源区岩石圈的深度约为46 km，莫霍面深度约为10 km．

俯冲带对“410”和“660”的影响是目前研究的一个热点(Li et al，2008)． 该研究涉及到间断面的性质，又能反映地幔对流的性质和形式． 俯冲带是冷的下沉物质，其通过间断面时会引起间断面的形态变化． 若“410”为相变界面，俯冲带会使它抬升(Vidale，Benz，1992； Collier，Helffrich，1997); 若“660”为相变界面，俯冲板块中的冷物质将使它下沉． 因此，汤加—斐济俯冲带中冷的物质与“410”作用时，俯冲带及相邻地幔中的橄榄石-尖晶石的相变界面会上升(蒋志勇等，2003)． 本文中得到的汤加—斐济俯冲带“410”间断面抬升约11.5 km与该结论一致．

本文使用倾斜叠加技术成功提取了sP前驱波(s₄₁₀P，s_LABP和s_MohoP)，获得了俯冲带间断面的深度，给出了汤加—斐济俯冲带新的地震波速度模型(IASP91_Tonga)． 本文结果将为深入研究该区域的地球结构提供基础资料，对正确认识俯冲带的复杂结构以及深入理解地球深部的动力学过程具有重要的参考价值．