The seismometer-detached ocean bottom seismograph and its experiments in the Gakkel Ridge，Arctic Ocean
-
摘要:
中国第十二次北极科学考察在密集浮冰覆盖的北冰洋加克洋中脊完成了国际首次大规模主动源海底地震仪探测，成功回收了43台海底地震仪中的42台，其中5台分体式海底地震仪均成功回收。本文介绍了针对该航次中面临的海底地震仪回收定位难题研发的一款新型分体式海底地震仪，其主要特点为：① 集成超短基线信标和声学应答系统，实现双重定位保障；② 使用分体式结构增强地震计与海底的耦合，有效提高信噪比；③ 选用国产地震计和浮力材料，实现核心部件国产化。分体式海底地震仪在北冰洋加克洋中脊的试验结果显示：地震背景噪声能量低，且记录的远震波形、两个近源微地震波形和三分量主动源地震记录均具有清晰可见的震相信息，表明该分体式海底地震仪能够满足冰下海底勘探需求。
Abstract:The ocean bottom seismograph （OBS） is an indispensable instrument for investigating the deep crust-mantle structure. China’s 12th Arctic scientific expedition marked a significant milestone in marine seismic exploration by achieving the first large-scale active exploration along the challenging Gakkel Ridge in the Arctic Ocean, where ice coverage poses formidable obstacles. Impressively, the expedition successfully recovered 42 out of 43 deployed OBSs, with all five seismometer-detached OBSs retrieved intact. This paper introduces a seismometer-detached OBS specifically developed to overcome the unique challenges encountered during the expedition, particularly in terms of OBS recovery and positioning. Its main features include:
1） The newly developed seismometer-detached OBS integrates cutting-edge technology to enhance its functionality and performance under harsh Arctic conditions. One of its key innovations is the combination of ultra-short baseline beacons and advanced acoustic response short baseline array positioning systems, which improves positioning accuracy through dual positioning technology. The positioning accuracy error of the short baseline array is within ±100 m, and the mature Sonardyne ultra-short baseline array loaded on the Xuelong-2 polar scientific research icebreaker can achieve a positioning accuracy up to one thousandth of the water depth. The use of such a dual positioning system can ensure accurate positioning of the instrument even in areas with dense ice layers. This ground-breaking design significantly improves the accuracy and efficiency of OBS recovery operations, enabling researchers to locate and recover instruments with unprecedented reliability and precision.
2） Referencing the current status of seismometer-detached OBSs both domestically and internationally, a titanium alloy seismic instrument chamber with a buoyancy of 17 kg in water has been designed. This chamber is incorporated into the seismic instrument separation structure, independently installed inside the OBS hull via flexible cables. Upon reaching the seafloor, it is timed to release, optimizing the coupling between the seismic instrument and the seabed. This structure also serves to reduce the impact of bottom currents on the seismometer. By enhancing this coupling, the OBS effectively improves the signal-to-noise ratio for recording seismic data, thus minimizing ambient noise and interference. Consequently, researchers can obtain clearer and more accurate seismic records, thereby facilitating a deeper understanding of the Earth’s geological processes and tectonic activities.
3） The OBS utilizes advanced domestic seismic instruments and buoyancy materials. On one hand, this provides new options for the materials used in OBS development, to some extent alleviating the problem of insufficient component production, which is advantageous for the large-scale production and industrialization. On the other hand, it enhances the flexibility of the OBS’s design, enabling the possibility of loading multi-functional modules onto the OBS. Additionally, using domestically produced seismometers makes it easier to optimize and develop hardware and software according to scientific research requirements. This signifies a significant step towards achieving self-sufficiency in marine instrument technology. This domestic production capacity not only enhances China’s scientific research capabilities, but also promotes innovation and technological advancement in the field of marine instruments, establishing independent intellectual property rights for key technologies of marine seismographs.
China’s 12th Arctic scientific expedition has yielded promising results, with the seismometer-detached OBS demonstrating exceptional performance in recording seismic signals. It was capable of collecting seismic data in the low-level horizontal seismic ambient noise environment along the Gakkel Ridge of the Arctic Ocean. Notably, the OBSs have exhibited low horizontal seismic ambient noise levels, underscoring their suitability for seismic exploration in ice-covered marine environments. Additionally, these OBSs successfully collected the waveform records of a small teleseismic, two micro-earthquakes and three-component active seismic exploration. This validates to some extent the effectiveness of the seismometer-detached structure in improving the signal-to-noise ratio, indicating that this type of OBS can meet the requirements for submarine exploration under the ice.
The successful development and application of the seismometer-detached OBS provide valuable experience for future submarine seismic exploration in extreme environments. Looking ahead, there are opportunities for further enhancement in various aspects such as real-time data transmission, longer battery life, and integration of more modularized sensors. The use of flexible buoyancy materials also liberates the submarine seismograph from the limitations imposed by traditional glass chambers. This enables the configuration of different types of modularized sensors, thus paving the way for the development of a versatile submarine observation platform with powerful functionality.
-
引言
地震预警系统具有减轻地震灾害的能力，在过去的数十年里，地震预警已经逐渐发展成为减轻地震灾害的一个重要手段（Nakamura，1988；Allen，Kanamori，2003；Zollo et al，2010 ；金星等，2012）。世界上许多地区建设了由实时监测地震的台站和中心处理系统构成的地震预警系统，包括日本（Odaka，2003；Horiuchi et al，2005 ；Kodera et al，2016 ）、美国加州（Allen，Kanamori，2003；Allen et al，2009 ；Böse et al，2009a ）、吉尔吉斯斯坦（Pittore et al，2014 ）、墨西哥（Espinosa-Aranda et al，1995 ，2009）、意大利（Zollo et al，2006 ）、罗马尼亚（Böse et al，2009b ）和伊比利亚半岛西南部（Pazos et al，2015 ）。我国北京、福建和台湾地区也陆续建成了地震预警系统，开展了一系列测试工作（Wu et al，1999 ，2013；马强，2008；Peng et al，2011 ；Chen，2015）。目前我国的国家地震烈度速报与预警工程正处于实施阶段，预计2020年投入测试工作，有望在全国范围内提供地震预警与烈度速报服务。
大地震往往沿发震断层主破裂方向产生较长距离的破裂，如MW7.0约60 km，MW7.5约130 km，MW8.0约300 km （Wells，Coppersmith，1994）。地震预警系统在对地震的震中、震级等参数进行估测时，通常将地震震源视为一个点源，这样的点源假设在实际应用中可以较好地表示较近距离范围内中小震的震源，但对于罕遇的大震（M＞7.5），则会出现明显的偏差。2011年日本东北（Tohoku） MW9.0地震中，在仅使用P波信息的条件下，地震预警系统给出的估测震级仅为M5.7，东京地区的烈度估测结果为日本气象厅地震烈度（Japan Meteorological Agency，简写为JMA）4度，由于预警系统未考虑断层破裂的长度，包括东京在内的日本关东部分地区均出现了烈度被低估的现象，东京地区的实际烈度达到了JMA5度（Hoshiba et al，2011 ）。地震预警系统在面对这样的大震时，如果能够迅速地估测出地震的破裂长度，则可通过估测的破裂长度修正震级和地震动参数，从而降低低估的可能性。
大地震破裂长度的快速估测近年来一直是地震预警系统研究的重要方向，而且相关研究也取得了一些突破。土耳其伊斯坦布尔的预警系统建设中采用PreSEIS方法（Böse，2006），该方法基于神经网络方法，将点源方法扩展到随机有限断层模型中，通过对历史有限断层模型数据集的训练，学习得到断层位置、破裂长度、地震震级等信息。这种方法能够在地震发生时给出地震破裂长度的估测值，且估测速度与现有的原地预警方法接近。Böse等（2012）开发了有限断层破裂探测器“FinDer”算法，该方法基于密集的地震台网和一组预先计算好的断层模板，在地震发生时实时测定地面运动的振幅，采用图像识别技术自动实时计算断层破裂（假设为线源）在地表投影的中心位置、长度和走向；Zhang等（2014）研究了利用近场强震动数据进行地震破裂过程实时成像的方法。近年来，随着近场台阵反演震源破裂过程的研究取得较好成果（Fletcher et al，2006 ；Vallée et al，2008 ；Meng et al，2012 ），研究人员开始尝试利用近场台阵来为地震预警系统提供信息。Meng等（2014）提出了利用近断层密集台阵波束形成技术反演震源方位角的方法来实时估测震源的破裂长度，用美国的UPSAR强震动台阵记录对多次地震进行破裂长度的估测，取得了较好的结果。该方法需要小孔径地震台阵提供精确的走时，UPSAR台阵能够较好地满足这一要求，但我国现有的一些强震动观测台阵（例如，在汶川地震中的四川自贡强震台阵）中各观测子台的观测仪器均为独立授时，各子台记录可能存在时钟同步误差问题，直接应用该算法还需要寻找适合每个台站的延时确定方法。因此，本文拟基于波束形成技术开发一套通过反演破裂主能量源方位角的震源破裂长度快速估算算法，同时提出调整各台站相位的方法，即在每个台站的延时无法确定的情况下，能够近似地得到每个台站之间的相位差；并在此基础上对2008年汶川地震四川自贡强震台阵记录到的强地震动记录进行分析，以探讨该方法对汶川地震的断层破裂长度进行实时估测的可行性。
1. 方法
利用近断层密集台阵记录反演破裂主能量源方位角的方法来实时估测震源的破裂长度，涉及波束形成技术，而其中的密集台阵中台站之间时间延迟（也就是相位差）的确定及其精度则是关键所在。
1.1 波束形成技术
波束形成技术是阵列信号处理的重要部分，该方法通过将各阵元输出进行加权求和，以达到增强期望信号、抑制干扰的目的。利用阵列信号处理技术分析相干信号的相位，可以确定波在阵列上的传播方向（张小飞等，2015），因此，可以通过密集地震台阵信号确定某一时刻地震波的传播方向，进而确定产生这部分地震波的破裂源方向。破裂源方向射线与预定义的线性断层的交点则可视为破裂的主能量源位置，也可将其视为地震的瞬时震中。基于这些瞬时震中即可估测断层破裂的尺度，即将每一时刻计算出的瞬时震中之间的最大距离视为断层破裂长度。
在波束形成过程中采用每两组波形之间的互相关系数进行叠加来代替原始波形的叠加，能够减少地震波在地壳浅部环境传播过程中的散射和多路径所造成的不利影响（Borcea et al，2002 ；Fletcher et al，2006 ）。因此，对于一个选定的时窗，将每两组地震动记录之间的互相关系数cij（θ）定义为
${c_{ij}}( \theta ) = \frac{{\sum\limits_t {{x_i}(t){x_j}(t - {\tau _{ij}}( \theta ))} }}{{\sqrt {(\sum\limits_t {x_i^2\sum\limits_t {x_j^2} } )} }},$
(1) 式中，xi为第i个台站记录到的地震信号，t为时间，τij（θ）为台站i与台站j之间的时间延迟，
${\tau _{ij}}( \theta ) = { \theta} \bullet ({r_i} - {r_j}) + \delta {t_i} + \delta {t_j},$
(2) 式中：θ为慢速矢量，对于单频波，其慢速矢量定义为大小为波速的倒数、方向与单频波传播方向相同的矢量；δti和δtj分别为台站i和j本身的记录延迟。将求得的τij（θ）舍入到最近的采样点上用于计算cij（θ）。
这样，按照式（1），在每一个时间窗下，求出每两组地震动记录之间的互相关系数cij（θ）并求和得到互相关系数的和cij s （θ）。由于慢速矢量包含着地震波传播视速度和反向方位角的信息，cij s （θ）是慢速矢量θ的函数，因此对于每一个时窗而言，当cij s （θ）的极值大于某个规定的阈值时，就认为此时θ所对应的方位角是地震主要辐射体所在的方向。Fletcher等（2006），Vallée等（2008）和Meng等（2012）均利用这种方法来推断破裂长度、断层分支及破裂速度，但其中涉及较多的人为推断，在实时应用中难以实现。
1.2 改进的相位调整方法
波束形成方法需要台阵中每个台站准确的延时数据δti。然而，现有的强震动观测台阵的强震动记录可能存在丢头或时钟同步误差较大等问题，因此在对各地震动时程进行分析之前，需要确定台阵中每个台站相对于标准台站的相位差。基于这个需求，Fletcher等（2006）曾提出利用大量余震、小震数据统计分析计算每个台站延时的方法，并取得了良好的效果；但对于本文研究所用的已获取汶川地震主震记录的自贡强震动观测台阵而言，能够获得的余震、小震记录少，无法通过该方法来计算台阵的延时。另外，Meng等（2012）曾利用10 s时窗下的相关性分析方法进行相位调整，但该方法旨在分析断层的破裂过程，10 s时间窗的设置对于分秒必争的地震预警系统显得尤为过长。
为解决上述台阵各台站延时确定的问题，本文提出相位调整的方法，具体步骤如下：
1） 利用Allen （1978，1982）的P波捡拾方法，对地震动时程的垂向分量进行首波初步捡拾；
2） 采用Akaike （1973）提出的赤池信息准则（Akaike information criterion，简写为AIC）进行首波的精确捡拾；
3） 给每一条地震动时程添加一个3 s的时窗，时窗起始点为首波的捡拾点；
4） 进行0.2—0.7 Hz带通滤波；
5） 利用相关性分析方法剔除丢头的数据；
6） 调整每个台站间水平向、垂向空间位置所造成的相位差，其中：水平向空间位置造成的相位差近似为台站空间位置矢量在点源震中到台阵的方位角上的投影长度除以莫霍面上的纵波波速；垂向空间位置造成的相位差近似为台站的高程差除以上地壳的纵波波速。
经过上述步骤，能够剔除丢头的地震动时程，获取满足波束形成技术的地震动时程。
图1给出了本文方法的技术路线，可以看出，该算法主要分为两个部分：第一部分为数据筛选与相位调整，该部分依照相位调整方法估算各台站相对于参考台站的相位差，并根据相位差调整每组数据的相位；第二部分为波束形成方法，按照1.1节所介绍的方法对断层破裂长度进行估算。
2. 基于汶川地震强震动台阵记录的应用分析
自贡地形影响强震动观测台阵是“十五”强震动台网系统建设项目的组成部分，于2007年建成并实施观测运行。该台阵位于我国四川省自贡市西山公园内的山体上，由8个台站组成，每个台站均配置了ETNA型记录器和ES-T型加速度计的数字强震动仪（杨宇等，2011；唐晖等，2012）。山顶台站（6号）与山底台站（0号或1号）的最大高程差达到72 m，4号台站与5号台站之间的高程差最小，为2 m。0号台站位于山底土层之上，1号台站位于山底基岩上，其它6个台站位于山脊不同高度的侏罗纪基岩上，各台站之间的空间相对距离处于52—400 m之间（阿布都瓦里斯，2013）。
![]()
图 2 自贡强震台阵记录到的汶川地震三分量加速度时程Figure 2. The three-component acceleration time histories of Wenchuan earthquake recorded by Zigong strong motion array自贡台阵位于汶川地震震中东南约200 km，在汶川地震中，该台阵获得了良好的主震加速度时程，各台站的加速度时程如图2所示，可见记录到的地震动时程的近震特征明显。
利用地震台阵估测地震主能量源反向方位角的方法是基于地震台阵信号的相位和相干性来追踪产生最强辐射源的区域。本文将2号台站作为参考台站来讨论台站之间的相位差。从图2中可以看出，1号台站的地震动时程相比其它台站数据存在明显的相位前移。为了筛选可用的地震动记录并对齐首波，需要估算地震波的到时。采用吴建平等（2009）研究中的汶川地区地下速度结构（图3b）以及刘瑞丰等（2014）的震相标识方法，粗略地估计出地震动记录各标志震相的到时；然后根据到时估算出P波与S波的到时差约为27 s。图3a为1号台站地震动时程的S波标志震相标记，可见Sn，Sb和Sg震相大概出现在15—21 s，考虑到估算的P波与S波的到时差，推断1号台站的地震动时程出现了丢头的现象，再加上波束形成方法需要对齐初至的Pn波相位，所以1号台站的地震动时程不能满足本研究需求。
台阵中的每个台站均存在各自的延时，其延时δt可以由大量的余震记录通过统计方法获得，由此精确地校正地震动时程的初始相位。由于可获得的自贡台阵数据量较少，我们无法通过余震记录的统计方法获得每个台站的δt，所以首先对齐首波，再根据各台站的空间相对位置确定每个台站与参考台站之间的相位差。本文采用P波捡拾方法（Allen，1978，1982），对地震动时程的垂向分量进行首波的初步捡拾，然后采用AIC方法（Akaike，1973）进行首波的精确捡拾，以此来对齐时程。
从各个台站之间Pn波到达后的初始10 s时间窗下波形的相关性（图3c）来看：由于1号台站记录到的地震动时程存在丢头的现象，其与其它台站之间的波形相关性非常小；而土层上的0号台站由于受到放大效应的影响，其地震动时程与除1号外其它6组地震动时程的相关性较差。在台阵中台站数量足够大的情况下，阵列数据处理中一般选择相关系数较高的地震动时程（Meng et al，2012 ）。在实时应用中，本文通过设置相关性系数的阈值，实现了自动择取相关性较好的、适合估测断层长度的数据的算法。考虑到本文可用的台站数据数量较少，计算中将阈值设置为0.3，剔除1号台站数据，选用了其余7个台站的数据。
![]()
图 3 强震动记录的筛选过程示意图（a） 1号台站前60 s三分量强震动记录及其S波标志震相；（b） 本文采用的汶川地区地下速度结构；（c） 自贡台阵记录到的汶川地震各台站加速度时程之间的相关性Figure 3. Diagram of the screening process for strong ground motion records（a） The pre 60 s three-component strong motion records of the station No.1 and the seismic phase of the S wave；（b） Underground velocity structure of seismic wave in Wenchuan area；（c） Correlation between acceleration time histories of each station recorded by Zigong array during Wenchuan earthquake由于台站的空间位置不同，对齐首波的时程并不能反映各台站之间的真实相位差，因此本文依据震中到台阵的方位角和各台站的空间位置计算各台站与参考台站之间的相位差，并根据相位差对各台站的波形相位进行调整。水平向的相位差由各台站的经纬度坐标、点源震中到台阵的方位角以及莫霍面的波速决定；垂向的相位差理论上可以通过多次余震的数据统计得到（Fletcher et al，2006 ），但由于自贡台阵缺少相关数据，故本文选用台站的高程差除以上地壳波速近似得到垂向的相位差。这样最终确定了各台站相对于参考台站（2号）的相位差，以此对各个台站记录的相位予以调整。
由于大地震中近场地震动峰值高且中低频成分丰富，同时含有较多的低频成分（周期1 s以上）（李小军，2001），再考虑到对于大地震的近场成像研究，采用低频成分往往能够较准确地估测震源尺度（Meng et al，2012 ），因此本文对地震动时程进行0.2—1 Hz的带通滤波，最终得到的加速度时程如图4所示。在3 s时间窗下，以1 s为滑移间隔，通过寻找每一秒的相关系数和的峰值位置来确定破裂主能量源方位角。预先定义龙门山断裂带为过震中、走向为225°的线型断层，则自贡台阵朝向该方位角的射线与这个预定义断层的交点即为破裂主能量源的位置。
![]()
图 4 经过震相对齐、滤波和相位调整后的地震动加速度时程Figure 4. The time histories of ground motion after seismic phase alignment，filtering and station delay processing图5给出了离线条件下汶川地震的破裂长度估测结果随时间的变化情况。图5a中，颜色代表断层各位置的相关系数和，圆点的大小代表每一时刻叠加的相关系数峰值，圆点越大说明能量越集中。地震波刚刚到达时，能量的来源集中，相关系数峰值较大，而加速度时程在100 s后很少出现大的相关系数峰值，这是由于这个阶段各种传播路径的地震波的能量叠加、成分和能量来源非常复杂所造成。
图5b为参考台站（2号）的垂向加速度时程，可以看到，地震动加速度时程显示在约18 s时Pn波传播至该台阵，约45 s时Sn波到达，约78 s时地震动时程达到了其加速度峰值（其它台站记录的地震动加速度峰值均在该时刻附近），此时图5a中估测出的破裂断层从震中向东北和西南方向各传播60 km、总长约为120 km的双向破裂，得到估测结果的时长为65 s。
![]()
图 5 离线条件下汶川地震瞬时主能量源实时估测结果（a）及同时刻2号参考台站的垂向加速度时程（b）图（a）中圆点的位置代表这一时刻的主能量源位置，圆点的大小代表这一时刻相关系数和的极值，虚线表示两阶段估测的破裂长度范围Figure 5. The outline result of the real-time estimation of the main energy source of Wenchuan earthquake （a） and the vertical acceleration record of the reference station No.2 at the same time （b）In Fig.（a），the dot position represents the location of main energy source at the moment，the dot size represents the maximum size of correlation coefficient at the moment，and the dotted lines represent range of the rupture从图5还可以看到，65 s时观测点的地震动远未达到该地的地震动最大值（位于70—80 s），说明在地震动加速度峰值到来之前便可获得较合适的破裂长度估测结果，这一先于地震动最大值出现之前的破裂长度估测对于地震预警系统非常重要。进一步可看到，120 km的断层破裂长度是由45—65 s时段的主能量估测结果所决定，反映了汶川大地震最主要的一次能量释放的破裂长度，该结果与张勇等（2008）的结论相似，表明本文方法能有效地给出最主要破裂过程的估测结果。根据张勇等（2008）的研究结果，汶川地震最主要的能量释放事件的发生时间介于震后14—34 s之间，该过程所释放的地震矩约占整个破裂过程地震矩的60%，相应滑移量超过3 m的断层长度约为120 km。考虑到地震波传播的时间因素，我们认为图5中45—65 s的估测结果较好地描述了发震后14—34 s这一阶段的断层破裂长度。此外，观察到该方法容易低估远离台阵方向上的破裂长度。120 s时利用本文方法得到的断层破裂长度为230 km，其中向南西延伸约170 km，向北东延伸约60 km，这个结果与其它研究结果相比有一定的差距，除了较张勇等（2008）的结果呈现较为明显的低估之外，还缺失一些震中东北方向上的重要破裂过程。从图5的结论来看，破裂过程到达42 s时，震中西南方向约150 km处出现了新的破裂，几乎在同一时间，虽然发生在北东方向上的新破裂的滑动量更大，但是这一过程未被探测出来。本文认为是震中北东方向上的破裂距离台阵较远，能量衰减较多，且该阶段中传播到达台阵的地震波成分复杂所致，所以北东方向上的破裂没有显现出来，而西南方向上的破裂距离台站较近，破裂过程被记录了下来。Meng等（2014）关于2012年El Mayor-Cucapah地震的研究中也出现了类似的现象。因此，这一方法对于汶川大地震主要破裂过程的破裂长度能够得到较好估测结果的同时，会在后期的估测过程中出现离台站较远的北东方向上破裂长度呈现明显的低估。这一现象与台阵相对于发震断层的距离较远有关。若在震中的北东、西南方向上有相应的台阵记录，将能够解决破裂长度低估的问题。这就是基于地震动观测台阵记录开展大地震断层破裂长度实时估测所面临的的限制条件。
本文所研究的断层破裂长度估测方法在地震预警原地、异地信息检测模式上有一定的局限性。原地地震预警系统是根据P波波速比更具破坏性的S波、面波波速快的物理原理，利用P波来估测未来S波、面波的影响程度。从图5中可以看到，本文的这种方法在S波到来后才能得到较为科学的估测结果，所以难以在原地地震预警系统中发挥作用。而异地地震预警系统是根据离断层较近的估测结果来为破坏性地震波还未到达的地区提供预警信息，本文方法能够基于近断层的地震动观测台阵记录计算得到断层破裂长度信息，以降低预警系统低估地震动强度和强烈地震动分布范围的风险。
如果要在S波到来前对震源长度进行有效的估测，可以从早期P波中寻找信息。如图5所示，在断层破裂长度的估测值达到120 km时，虽然地震动峰值加速度还未出现，但S波首波在约20 s前已经到达，这对于大地震来说是相当危险的，尤其是在近场，地震动峰值大，S波到达后很快会造成严重的破坏。根据早期的P波信息所得到的断层破裂长度为从震中向南西方向传播约50 km，由于这种方法对初始震中位置的估测出现了偏差，所以在破裂方向的估测上同样出现了偏差；但破裂过程得到了较为完整的记录，即从初始震中向北东方向传播，再从北东方向向南西方向传播。如果利用P波信息对汶川地震的主要破裂过程进行推断，结合初始震中的位置，可以得到汶川地震的初期破裂长度为从震中向北东方向延伸50 km，这虽然与利用S波进行估测的结果相比较小，但得到估测结果在时间上要比利用S波的早约40 s，而且这一结果也能够反映出地震破裂传播的主要方向。现阶段采用近场P波信息对震源破裂长度的估测研究还较少，这种估测方法仍需更多的大地震台阵记录进行准确性验证。
震源破裂长度的估算若要满足预警系统的需求，除了在危险来临前获取计算所需的数据，对算法本身的计算速度也有较高的要求。本文的时间窗漂移量设为1 s，那么，对于每一个时间窗的计算时间应在1 s以下，甚至要远小于1 s，以便为数据读取、通讯预留时间。
对算法本身的计算速度进行测试分析，具体的计算环境为：Intel i7-2600 k中央处理器，4核8线程，最大主频3.40 GHz；可用内存3.47 G；32位win7操作系统；Matlab2011a软件；并行算法的计算节点为8个。
算法速度的测试结果显示：数据准备过程用时4.149 7 s，利用断层实时估测的串行算法计算每1 s的断层破裂位置用时不超过2.4 s。对计算速度有显著影响的因素为：① 台阵的台站数量ns，用于计算每两个台站间的相关系数，时间复杂度为O（ns 2）；② 断层的分割数量nf，时间复杂度为O（nf）；③ 慢速矢量的间隔数量ls，时间复杂度为O（ls）；④ 时间窗长度nt。总的计算量为O（ns 2·nf·ls·nt），通过对整个断层分割进行8线程的并行计算，能够有效地提高计算速度。在并行计算初期，CPU需要经过一个加速的过程，之后该算法得到每一秒估测结果的时间均能够稳定在0.8 s附近，因此该算法理论上能够满足实时断层长度估测的应用。
3. 讨论与结论
本文采用近场地震动波束形成技术，利用自贡地形地震动观测台阵在汶川地震中得到的地震动时程，通过反演破裂主能量源方位角的方法开展了汶川地震发震断层破裂长度的实时估测，主要结果如下：
1） 完善了利用近场小口径地震动观测台阵实时估测断层破裂长度的技术路线，开发了用于估测断层破裂长度的算法，提出了调整各台站相位的方法，即在每个台站的延时无法确定的情况下，仍能够近似地得到每个台站之间的相位差。
2） 在地震动峰值加速度到来前约10 s，该方法得到的断层破裂长度为120 km，从震中向北东、南西各传播60 km，这一结论较好地反映了汶川大地震最主要的能量释放过程的断层破裂长度。
3） 受能量衰减因素的影响，该方法容易低估远离台阵方向上的破裂长度。最终的破裂长度估测结果约为230 km，主要从震中向南西方向延伸，这一结果低估了北东方向上的破裂，本文认为可通过在震中的北东、南西方向上各放置一个台阵的办法来降低破裂长度低估的可能性。
4） 讨论了使用近场P波估测破裂过程的可行性，通过该方法得到汶川地震的初期破裂为从震中向北东方向延伸50 km，且得到估测结果的时间比利用S波的早约40 s。基于P波信息的方法能够较准确地描述断层的初期破裂传播过程，但是对震中的估测存在一定的偏差。
5） 为提高实时估测震源破裂过程方法的实用性，分析了算法的计算速度。采用8线程并行计算算法，该算法下计算每秒的破裂位置用时约为0.8 s，能够为数据传输预留时间并满足实时破裂长度估测的应用。
本文的研究进一步表明，破裂长度估测方法仍然有其局限性，因此存在改进之处：第一，考虑到本文所利用的地震动观测台阵获得的强地震动记录数据较少，本文在处理台站之间高程带来的相位差时采取了近似的方法而不是利用大量数据进行拟合分析，这会给加速度时程的相位确定带来偏差；第二，相比美国的UPSAR台阵，自贡台阵的观测点数量少，台阵观测点口径也较小，且自贡台阵与龙门山断裂带相距较远，因此在计算主能量源位置时，空间不确定性较大，容易引起对破裂位置估测的偏移；第三，该方法虽然能够将主要的破裂过程估测出来，但是当一些局部破裂发生在远离台阵的断层走向上时，不容易被探测到，这会导致断层长度的低估。尽管如此，该方法仍然能够在大地震发生后第一时间内为其破裂长度的估测提供重要信息。
-
![]()
图 6 加克洋中脊分体式OBS15台站位置（a）及其记录的2021年海地MW5.8地震波形（b）
波形图使用0.04—0.2 Hz带通滤波进行了处理。下行小图分别为可识别的初至P波和PcP，PKiKP，初至S波和SKS和ScS的放大图。使用IASP91模型对上述地震波到时进行了预测，见图中红色和蓝色实线
Figure 6. The location of the seismometer-detached OBS15 station in the Gakkel Ridge （a） and the waveforms （b） of the Haiti MW5.8 earthquake in 2021
The waveforms were processed with bandpass filtering （0.04−0.2 Hz）. The figures in the second row are enlarged views of identifiable P-wave and PcP，PKiKP，S-wave and SKS，ScS waveforms，respectively. The arrival times of these seismic waves were predicted using the IASP91 model，indicated by the red and blue solid lines in the figures
![]()
图 1 新型宽频带分体式海底地震仪的结构模型
（a） 外视图；（b） 内部透视图
Figure 1. The structural models of the broadband seismometer-detached ocean bottom seismograph
（a） The external view；（b） The internal perspective view
![]()
图 2 短基线阵及其定位结果
（a） 短基线阵定位系统在雪龙2号月池车间的安装情况；（b） 短基线阵试验定位三维结果
Figure 2. The short baseline array and its positioning results
（a） The installation of the short baseline array positioning system in the moon pool workshop of the Xuelong-2；（b） The three-dimensional results of the short baseline array experimental positioning
![]()
图 3 超短基线在北冰洋加克洋中脊分体式海底地震仪回收过程中的三维实时定位结果
图（a）和（b）分别为分体式海底地震仪OBS9和OBS15的上浮轨迹
Figure 3. The three-dimensional real-time positioning results during the recovery process of the seismometer-detached OBS with ultra-short baseline in the Gakkel Ridge of Arctic Ocean
Figs. （a） and （b） show the ascending tracks of the seismometer-detached OBS9 and OBS15
![]()
图 4 分体式海底地震仪及其在加克洋中脊的试验
（a） 2021年中国第十二次北极科学考察布设和回收的海底地震仪位置和炸测位置示意图；（b） 分体式海底地震仪的投放；（c） 分体式海底地震仪的回收
Figure 4. The seismometer-detached OBS and its experiments in the Gakkel Ridge
（a） The locations of OBS deployment and recovery and air-gun shooting lines during the Chinese 12th Arctic scientific expedition in 2021；（b） Deployment of the seismometer-detached OBS；（c） Recovery of the seismometer-detached OBS
![]()
图 5 （a） 分体式OBS3垂直分量记录于2021年8月10—24日在北极加克洋中脊的功率谱密度分布；（b，c） 分体式OBS3两个水平分量的功率谱密度分布
图中彩色细实线代表每小时记录对应的功率谱密度，不同颜色代表分布在不同能量区间的概率，黑色粗实线为平均值，灰色粗实线为Peterson （1993）提出的全球背景噪音模型参考线
Figure 5. （a） The vertical component records of the OBS3 in the Gakkel Ridge of the Arctic collected from August 10 to 24，2021 presented in terms of power spectral density probability density functions；（b，c） The power spectral density distribution of the two horizontal components of the OBS3
In the figures color thin solid lines represent the power spectral density corresponding to each hourly record，with different colors indicating probability distribution in various energy ranges. The thick solid black line represents the average，and the thick solid gray lines serve as the reference lines for the global ambient noise models proposed by Peterson （1993）
![]()
图 7 加克洋中脊分体式OBS3 （a）和OBS17 （b）台站记录的微震信息
图中波形已经带通滤波（3—12 Hz）处理，红色和蓝色实线标注了使用IASP91模型对上述地震波到时进行的预测
Figure 7. Microseismic information recorded by the seismometer-detached OBS3 （a） and OBS17 （b） along the Gakkel Ridge
The waveforms have been bandpass filtered （3−12 Hz）. The arrival times of the seismic waves were predicted using the IASP91 model，indicated by the red and blue solid lines
![]()
图 8 布放在加克洋中脊的分体式OBS3台站主动源三分量地震记录的折合时间剖面（折合速度为6.0 km/s）
（a） x分量；（b） y分量；（c） z分量
Figure 8. Active source three-component seismic record sections with reduced time of the OBS3deployed along the Gakkel Ridge of Arctic Ocean （The reduced velocity is 6.0 km/s）
（a） x component；（b） y component；（c） z component
![]()
图 9 布放在加克洋中脊的OBS15台站主动源三分量地震记录折合时间剖面
（a） x分量；（b） y分量；（c） z分量
Figure 9. Active source three-component seismic record sections with reduced time of the OBS15 deployed along the Gakkel Ridge of Arctic Ocean
（a） x component；（b） y component；（c） z component
表 1 宽频带分体式海底地震仪组件及其技术参数
Table 1 Components and technical parameters of the broadband seismometer-detached OBS
地震计频带范围 60 s—50 Hz 标定信号DAC 16位 地震计满量程 10 mm/s （单峰值） 标定信号类型 脉冲，正弦波可选，参数可设置 地震计灵敏度 2000 V/（m/s）（双端） 标定信号启动方式 定时、指令 ADC 24位 数据记录 32 GB×3 数据采样率 50 sps，100 sps，200 sps，500 sps，
每个采集通道可单独设定采样率电子罗盘动态精度 ±1o 数字滤波器 线性相位FIR，最小相位FIR 电子罗盘调平范围 45° 数据采集动态范围 ＞135 dB 数据通信接口 LAN以太网 授时 GPS，北斗 数据通信协议 TCP，IP 守时 芯片级原子钟 最大工作水深 6 km 水听器频带范围 10—2000 Hz 工作时间 6个月—1年 方位角精度 俯仰范围±30°±0.1°，±（30°—45°）±0.2° 磁场测量精度 10 nT 
下载: 导出CSV
表 2 短基线阵和超短基线信标的技术参数
Table 2 The technical parameters of short baseline arrays and ultra-short baseline beacons
深度级 频段 收发机波束角 测距精度 短基线阵 4000 m12 kHz 半指向性 100 m 超短基线信标 7000 m19—34 kHz 半指向性 ＜15 mm
下载: 导出CSV
-
郝天珧,游庆瑜. 2011. 国产海底地震仪研制现状及其在海底结构探测中的应用[J]. 地球物理学报,54(12):3352–3361. Hao T Y,You Q Y. 2011. Progress of homemade OBS and its application on ocean bottom structure survey[J]. Chinese Journal of Geophysics,54(12):3352–3361 (in Chinese).
郝天珧,游庆瑜,王元,郭永刚,丘学林,黄松,徐亚,赵春蕾,张妍,徐锡强. 2022. 国产海底地震探测装备技术研发与应用[J]. 科学技术与工程,22(34):15020–15027. Hao T Y,You Q Y,Wang Y,Guo Y G,Qiu X L,Huang S,Xu Y,Zhao C L,Zhang Y,Xu X Q. 2022. Development technology and its application of Chinese ocean bottom seismometer[J]. Science Technology and Engineering,22(34):15020–15027 (in Chinese).
李江,庄灿涛,薛兵,朱小毅,陈阳,朱杰,彭朝勇,叶鹏,梁鸿森,刘明辉,杨桂存,周银兴,林湛,李建飞. 2010. 宽频带海底地震仪的研制[J]. 地震学报,32(5):610–618. Li J,Zhuang C T,Xue B,Zhu X Y,Chen Y,Zhu J,Peng C Y,Ye P,Liang H S,Liu M H,Yang G C,Zhou Y X,Lin Z,Li J F. 2010. Development of broadband ocean bottom seismograph (OBS)[J]. Acta Seismologica Sinica,32(5):610–618 (in Chinese).
李守军,包更生,吴水根. 2005. 水声定位技术的发展现状与展望[J]. 海洋技术,24(1):130–135. Li S J,Bao G S,Wu S G. 2005. A practical overview and prospect of acoustic positioning technology[J]. Journal of Ocean Technology,24(1):130–135 (in Chinese).
刘丹,杨挺,黎伯孟,吴越楚,王宜志,黄信锋,杜浩然,王建,陈永顺. 2022. 分体式宽频带海底地震仪的研制、测试和数据质量分析[J]. 地球物理学报,65(7):2560–2572. Liu D,Yang T,Li B M,Wu Y C,Wang Y Z,Huang X F,Du H R,Wang J,Chen Y S. 2022. Seismometer-detached broadband ocean bottom seismograph (OBS):Development,test,and data quality analysis[J]. Chinese Journal of Geophysics,65(7):2560–2572 (in Chinese).
牛雄伟,王力,丁巍伟,黄松华,孔凡圣,冯海泓,谭平川,史淼奇,傅建华,陈尚国. 2022. 一种适用于冰区的短基线承载结构:中国,202111670786.4[P]. [2022−06−21]. https://pss-system.cponline.cnipa.gov.cn/documents/detail?prevPageTit=changgui. Niu X W,Wang L,Ding W W,Huang S H,Kong F S,Feng H H,Tan P C,Shi M Q,Fu J H,Chen S G. 2022. Short base line bearing structure suitable for ice region:CN,202111670786.4[P]. [2022−06−21] https://pss-system.cponline.cnipa.gov.cn/documents/detail?prevPageTit=changgui (in Chinese).
阮爱国,李家彪,冯占英,吴振利. 2004. 海底地震仪及其国内外发展现状[J]. 东海海洋,22(2):19–27. Ruan A G,Li J B,Feng Z Y,Wu Z L. 2004. Ocean bottom seismometer and its development in the world[J]. Donghai Marine Science,22(2):19–27 (in Chinese).
阮爱国,牛雄伟,吴振利,吴招才,薛彬. 2009. 潮汕坳陷中生代沉积的折射波2D速度结构和密度[J]. 高校地质学报,15(4):522–528. Ruan A G,Niu X W,Wu Z L,Wu Z C,Xue B. 2009. The 2D velocity and density structure of the Mesozoic sediments in the Chaoshan depression[J]. Geological Journal of China Universities,15(4):522–528 (in Chinese).
阮爱国,李家彪,陈永顺,丘学林,吴振利,赵明辉,牛雄伟,王春龙,王显光. 2010. 国产I-4C型OBS在西南印度洋中脊的试验[J]. 地球物理学报,53(4):1015–1018. Ruan A G,Li J B,Chen Y S,Qiu X L,Wu Z L,Zhao M H,Niu X W,Wang C L,Wang X G. 2010. The experiment of broad band I-4C type OBS in the Southwest India ridge[J]. Chinese Journal of Geophysics,53(4):1015–1018 (in Chinese).
游庆瑜,刘福田,冉崇荣,王广福. 2003. 高频微功耗海底地震仪研制[J]. 地球物理学进展,18(1):173–176. You Q Y,Liu F T,Ran C R,Wang G F. 2003. High frequency micro-power ocean bottom seismograph[J]. Progress in Geophysics,18(1):173–176 (in Chinese).
张同伟,王向鑫,唐嘉陵,秦升杰. 2018. 深海超短基线定位系统现状及展望[J]. 舰船电子工程,38(10):1–6. doi: 10.3969/j.issn.1672-9730.2018.10.001 Zhang T W,Wang X X,Tang J L,Qin S J. 2018. Technical status and development trend of long range USBL[J]. Ship Electronic Engineering,38(10):1–6 (in Chinese).
Beyreuther M,Barsch R,Krischer L,Megies T,Behr Y,Wassermann J. 2010. ObsPy:A Python toolbox for seismology[J]. Seismol Res Lett,81(3):530–533. doi: 10.1785/gssrl.81.3.530
Christeson G L,Goff J A,Reece R S. 2019. Synthesis of oceanic crustal structure from two-dimensional seismic profiles[J]. Rev Geophys,57(2):504–529. doi: 10.1029/2019RG000641
Crawford W C,Webb S C. 2000. Identifying and removing tilt noise from low-frequency (＜0.1 Hz) seafloor vertical seismic data[J]. Bull Seismol Soc Am,90(4):952–963. doi: 10.1785/0119990121
Ding W W,Niu X W,Zhang T,Chen S G,Liu S H,Tan P C,Kong F S,Jin Z D,Huang S H,Wei C H,Fang Y X,Sun Q C,Li J B. 2022. Submarine wide-angle seismic experiments in the high Arctic:The JASMInE expedition in the slowest spreading Gakkel Ridge[J]. Geosyst Geoenviron,1(3):100076. doi: 10.1016/j.geogeo.2022.100076
Jokat W,Schmidt-Aursch M C. 2007. Geophysical characteristics of the ultraslow spreading Gakkel Ridge,Arctic Ocean[J]. Geophys J Int,168(3):983–998. doi: 10.1111/j.1365-246X.2006.03278.x
Jokat W,Ritzmann O,Schmidt-Aursch M C,Drachev S,Gauger S,Snow J. 2003. Geophysical evidence for reduced melt production on the Arctic ultraslow Gakkel mid-ocean ridge[J]. Nature,423(6943):962–965. doi: 10.1038/nature01706
Kennett B L N,Engdahl E R. 1991. Traveltimes for global earthquake location and phase identification[J]. Geophys J Int,105(2):429–465. doi: 10.1111/j.1365-246X.1991.tb06724.x
Korger E I M,Schlindwein V. 2014. Seismicity and structure of the 85°E volcanic complex at the ultraslow spreading Gakkel Ridge from local earthquake tomography[J]. Geophys J Int,196(1):539–551. doi: 10.1093/gji/ggt390
Koulakov I,Schlindwein V,Liu M Q,Gerya T,Jakovlev A,Ivanov A. 2022. Low-degree mantle melting controls the deep seismicity and explosive volcanism of the Gakkel Ridge[J]. Nat Commun,13(1):3122. doi: 10.1038/s41467-022-30797-4
Longuet-Higgins M S. 1950. A theory of the origin of microseisms[J]. Philos Trans Roy Soc A:Math Phys Eng Sci,243(857):1–35.
McGuire J J,Collins J A,Gouédard P,Roland E,Lizarralde D,Boettcher M S,Behn M D,van der Hilst R D. 2012. Variations in earthquake rupture properties along the Gofar transform fault,East Pacific Rise[J]. Nat Geosci,5(5):336–341. doi: 10.1038/ngeo1454
McNamara D E,Buland R P. 2004. Ambient noise levels in the continental United States[J]. Bull Seismol Soc Am,94(4):1517–1527. doi: 10.1785/012003001
Minshull T A,Muller M R,White R S. 2006. Crustal structure of the Southwest Indian Ridge at 66°E:Seismic constraints[J]. Geophys J Int,166(1):135–147. doi: 10.1111/j.1365-246X.2006.03001.x
Mosher D C,Chapman C B,Shimeld J,Jackson H R,Chian D,Verhoef J,Hutchinson D,Lebedeva-Ivanova N,Pederson R. 2013. High Arctic marine geophysical data acquisition[J]. Leading Edge,32(5):524–536. doi: 10.1190/tle32050524.1
Niu X W,Minshull T A,Li J B,Ruan A G,Wu Z L,Wei X D,Wang W,Li Y,Bayrakci G,Dong C Z,Ding W W,Fang Y X,Zhang J. 2023. Shear wave velocity structure and crustal lithology beneath the ultraslow spreading Southwest Indian Ridge at 50°E[J]. Geophys J Int,233(2):1416–1428. doi: 10.1093/gji/ggac516
Peterson J R. 1993. Observations and Modeling of Seismic Background Noise[R]. Albuquerque:U.S. Geol. Surv.:93−322.
Rekant P V,Gusev E A. 2016. Sediments in the Gakkel Ridge rift zone (Arctic Ocean):Structure and history[J]. Russ Geol Geophys,57(9):1283–1287. doi: 10.1016/j.rgg.2016.08.013
Schlindwein V,Müller C,Jokat W. 2007. Microseismicity of the ultraslow-spreading Gakkel Ridge,Arctic Ocean:A pilot study[J]. Geophys J Int,169(1):100–112. doi: 10.1111/j.1365-246X.2006.03308.x
Seher T,Crawford W C,Singh S C,Cannat M,Combier V,Dusunur D. 2010. Crustal velocity structure of the Lucky Strike segment of the Mid-Atlantic Ridge at 37°N from seismic refraction measurements[J]. J Geophys Res:Solid Earth,115(B3):B03103.
Stähler S C,Sigloch K,Hosseini K,Crawford W C,Barruol G,Schmidt-Aursch M C,Tsekhmistrenko M,Scholz J R,Mazzullo A,Deen M. 2016. Performance report of the RHUM-RUM ocean bottom seismometer network around La Réunion,western Indian Ocean[J]. Adv Geosci,41:43–63. doi: 10.5194/adgeo-41-43-2016
Sutton G H,Duennebier F K. 1987. Optimum design of ocean bottom seismometers[J]. Mar Geophys Res,9(1):47–65. doi: 10.1007/BF00338250
USGS. 2021. Earthquake lists,maps,and statistics[EB/OL]. [2023−08−18]. https://www.usgs.gov/natural-hazards/earthquake-hazards/lists-maps-and-statistics.
Webb S C. 1988. Long-period acoustic and seismic measurements and ocean floor currents[J]. IEEE J Oceanic Eng,13(4):263–270. doi: 10.1109/48.9239
计量
- 文章访问数: 165
- HTML全文浏览量: 17
- PDF下载量: 50
