徐纪人李海兵曾祥芝赵志新
1） 中国江苏东海 222300 江苏东海大陆深孔地壳活动国家野外科学观测研究站
2） 中国北京 100037 自然资源部深部地球动力学重点实验室
3） 中国北京 100037 中国地质科学院地质研究所
微小地震活动观测和研究对于了解地壳的原位实时地应力状态、研究地球内部结构和地壳及断裂活动、理解地球动力学过程和地震预测等均具有极为重要的意义（Liuet al，1992；
Ito，Zoback，2000）.利用宽频地震仪进行地震观测已取得了显著的成效，远距离传输地震观测也可以在非地震噪声干扰小的场所进行.然而，远距离传输地震观测难以观测到远离地震仪设置场所的研究区域的微小地震，且利用远距离传播的地震学方法解释地球构造特征的能力会随距离的增加而减弱.远距离传播的面波也会导致高频成分的丢失，而高频成分正是推测地壳和上地幔结构不可或缺的信息.
在研究区域安装地震仪进行微小地震活动监测，可在一定程度上弥补上述远距离地震观测的不足，然而许多研究区域却是观测地震的背景噪声高值区域.近年来的高频地震学研究也需要在人类活动频繁的都市背景噪声高值地区监测地震活动，因此地震学者一直在艰难地探索着如何在高噪声背景的地震活动区域观测地震才能获取高信噪比的地震波形.深井下安装地震仪进行地震观测可以减少地表的非地震噪声干扰影响，实现背景噪声高值地区微小地震活动的近距离观测，从而有效地提高观测波形信噪比（徐纪人，赵志新，2009；
Zeng，Yang，2021），因而深井地震学与地球物理学相结合的观测研究逐渐发展成为多学科的前沿性科学研究（涂毅敏，陈运泰，2002；
Baischet al，2002；
Ishiiet al，2003；
Okuboet al，2004）.深部地震观测方兴未艾（Tembeet al，2009；
龚正等，2020；
滕吉文，2021），其成果已被广泛地应用于地球科学诸多领域，极大地推进了现代地震学和地球物理学以及深井钻孔板块构造学等多学科的研究（李海兵等，2018）.近年来，断层井下地震观测研究在断层应力状态、震源过程及地震预报研究领域取得了前沿性成果（Chavarriaet al，2003）.钻孔地震观测有利于建立三维地表地震波速度模型（Layet al，2020），从而推进速度结构与地震关系等的研究（Jeppson，Tobin，2015）.
深井地震仪记录到大量M−0.5——0.5的微小地震事件，其中包括与中强地震相关的前震和余震，这些微小地震的记录促进了地震序列活动及其震源位置的研究（Lellouchet al，2019），其精确定位的地震活动对于构造研究具有重大意义 （Maet al，2012；
Shashidharet al，2020；
Hunget al，2022）.区域性深井长期地震观测广泛地应用于地震预测研究（Asaiet al，2005；
Fujimotoet al，2007），且在地震断层活动机制分析中发挥了重要作用（Ikariet al，2015）.此外，利用井下地震观测资料进行震源机制反演求得微小地震群震源的断层平面解，有助于建立地壳应力不均匀性模型（Langetet al，2020；
Luttrell，Hardebeck，2021）.
利用井下地震波测量资料结合岩石物理模型可为地震断裂特性的探索提供理论基础（Xuet al，2020）.井孔内微震波研究显示，应力场的方向与岩石脆性破裂特征有关（Pierdominiciet al，2020）.使用钻井地震资料可进行裂谷构造的地球动力学研究（Leiet al，2020）.深井地震观测数据结合岩芯信息的地震波数值模拟研究，可以提高复杂地下结构非均质性解释的准确度（Zhaoet al，2012；
Jeonget al，2020）.深井地震观测记录因受场地和衰减效应的影响比较小，能够为震源谱分析提供更直接的约束（Shearer，Abercrombie，2021）.根据井中入射波与反射波之间的传递函数与频率相关的衰减进行场地近地表剪切波衰减特性研究，有利于预测地震动特征（Rigaet al，2019）.井下地震应变波的研究也将促进高频地震学的发展（Zhanget al，2019）.
综上，高精度微小地震观测是未来地球科学发展的必然趋势.本文将主要探索高地震背景噪声的江苏东海深井观测站处于不同深度的地震仪记录微小地震的观测能力，以及近距离监测区域微小地震活动的有效途径，藉此分析地下不同深度地震仪观测到的地震波形信噪比的差异，研究不同深度地震波形的保真度，以期通过深井地震观测资料研究结果来促进地球科学的发展.
东海深井地壳活动观测站邻近郯庐断裂带（图1）.在东海中国大陆科学钻探的5 158 m深井中的不同深度安装地震仪进行观测，以便于中国大陆构造运动的研究（徐纪人等，2016）.造山运动时期形成的左旋走滑的郯庐断裂带不仅古地震活跃，当代地震活动性也很强（徐纪人等，2008）.1668年山东郯城M8.5地震就发生在郯庐断裂带上.现代震源机制解结果显示郯庐断裂带的当代构造活动呈右旋走滑运动（周翠英等，2013；
Xuet al，2021）.在复杂构造活动带进行地壳活动的监测研究，对于提高郯庐断裂带附近地震活动及我国东部地壳构造断裂带活动的监视能力具有重要意义.
图1 江苏东海地壳活动国家野外科学观测研究站周缘地质构造图（引自徐纪人等，2016）Fig.1 Geological structure around the Jiangsu Donghai Crustal Activity in Continental Deep Hole National Observation and Research Station in China （after Xu et al，2016）
江苏东海地壳活动观测研究站周围多有水晶采石活动及水晶岩石加工厂等，位于观测深井东南约3 km处的采石场时有爆破工程，附近的公路交通等也很繁忙.这些均会对地震观测造成背景噪声高值干扰，属于常规地面地震观测的不利因素.在该地区开展地面和井下地震活动观测，探索在背景噪声高值干扰地区地震观测的有效途径，有助于研究地壳活动和郯庐断裂带构造运动的特征.
江苏东海地壳活动观测研究站位于5 000 m深井中，井下最高温度达142.4 ℃ （何丽娟等，2006），井孔中的压力高达 50 MPa.第一期观测分别在井下 544.5 m，1 559.5 m，2 545.5 m三处不同深度安装了美国Sondi公司制造的主频为 4.5 Hz的高精度数字三分量地震仪，分别记为L1，L2和L3.该仪器无供电要求，可以在深井高温高压环境下连续工作.同时，在地表浅处设置了一台英国Guralp公司制造的3ESP型宽频三分量地震仪L0，该仪器的频宽为50 Hz——60 s，连同深处的三台地震仪组成了深井地震观测网（Xuet al，2016）.基于此观测网进行地面与井下不同深度地震仪记录到的地震波形特征及其各自的地震波形信噪比特征的研究.表1列出了各深度处地震仪参数及其围岩岩性等物理参数，其中地震波速为实验室岩心测量结果.
表1 深井中的地震仪参数及工作环境岩性和地球物理参数Table 1 Seismograph parameters and lithologic and geophysical parameters around the seismographs in the borehole
本文分析研究地震波形信噪比，采用了如下公式：
式中，r为地震波形信噪比，S为信号能量，N为噪声能量.噪声能量计算是截取初至波到时之前100 s的数据计算平均值，信号能量是截取横波到时之后记录中的最大振幅能量，噪声N或信号S各自的能量计算公式为
式中：E表示式（1）中计算的噪声N或信号S的能量；
AEW，ANS和AU分别为噪声N或者信号S的东西、南北和垂直分量的振幅.
3.1 地面高噪声背景下的微小地震观测
图2为东海观测站的地面和深井地震仪记录到的2013年2月3日江苏灌云ML0.8地震的波形图.该震震中位于（34.4°N，119.7°E），震中距为 97.2 km，震源深度为 5 km.由图2a的原始地震图记录可见东海观测站周围地面干扰有时相当严重，以至于从地面宽频地震仪L0记录中难以识别地震波形.干扰造成了微小地震的漏记现象，显示了地面设置的地震仪在高噪声背景条件下记录微小地震的局限性.为了与井下地震仪记录波形进行比较，图2b给出了图2a中原始地震图记录经过中通滤波的波形，可见滤波后的波形中滤掉了低频干扰，在时间轴大约07 ∶ 27处隐约可辨认出地震事件，但是P波和S波震相皆不清晰，整个波形序列中非地震信号干扰仍十分严重.而图2c中深井下2 554 m处的L3地震仪在同样噪声背景下能够记录到低噪声地震波形图，P波初动和S波波形清晰，显示出井下地震仪在高噪声地区近距离记录微小地震的优越性.
图2 东海观测站地面与井下地震仪记录的2013年2月3日江苏灌云ML0.8地震（震中距 97.2 km，方位角 80°）的三分量波形图（a） 地面宽频地震仪 L0记录波形；
（b） 图（a）波形经中通滤波后的波形；
（c） 井下地震仪 L3记录波形Fig.2 Three-component waveforms of the ML0.8 earthquake （epicentral distance is 97 km，azimuth is 80°） in Guanyun，Jiangsu on February 3，2013 recorded by the seismographs on the ground and in the borehole of the Donghai station（a） Waveforms recorded by the broadband seismograph L0 on the ground；
（b） Waveforms after medium-pass filtering to the waveforms in Fig.（a）；
（c） Waveforms recorded by the seismograph L3 in the borehole underground
3.2 地面与井下地震仪记录波形保真度分析
图3 显示了2013年5月24日安徽利辛附近ML2.7地震地面和井下不同深处地震仪记录波形.该地震位于（33.0°N，116.3°E），震中距为 271.5 km.图3 中井下地震仪L1和L3记录中基本无噪声干扰，地震波形三分量的P波初动及P波和S波波形震相都清晰可辨，各分量P波波形和S波波形分别相似.图3a中地面宽频地震仪L0的记录在震前和震后都有明显的长周期噪声干扰，这反映了东海观测站不时处于地面高噪声状态.本文观测资料显示，东海观测站非地震干扰强弱程度是无规则的.地面地震仪记录的垂直分量和近南北分量的P波初动及S波虽可辨，但波形与井下两地震仪的相应波形差别较大，尤其是其近南北分量振幅很小，近于畸变.该地震震中位于东海观测站西南方向，方位角为56°.畸变原因除与观测方位有关外，可能是受到了场地效应或者路径衰减的影响，也可能是受到地面噪声干扰的影响.这就意味着井下地震观测记录受场地或者浅层噪声的影响较地表地震仪小.为了进一步与井下主频为4.5 Hz的地震仪记录波形比较，图3b给出了图3a地面宽频地震仪波形经中通滤波后的波形，可见图3a中的低频干扰虽被滤掉，但与井下地震仪记录波形（图3c和3d）相较，其近南北分量的振幅仍然很小，近于畸变.结合上述井下不同深度地震仪波形的相似性分析可知，井下地震仪观测的微小地震波形具有良好的保真度，其观测质量明显优于地面台站，这可能给地震观测学带来启示.综合分析井下和地面地震观测也有益于推进地壳浅层结构特征和地面震动动力学方面的研究.
图3 东海观测站地面与井下不同深度地震仪记录的2013年5月24日安徽利辛ML2.7地震（震中距 271.5 km，方位角 56°）的三分量波形图（a） 地面 L0 地震仪记录波形；
（b） 图（a）中波形经中通滤波后所得波形；
（c） 井下 544.5 m 深处地震仪 L1 记录波形；
（d） 井下 2 545.5 m 深处地震仪 L3 记录波形Fig.3 Three-component seismic waveforms of the ML2.7 event in Lixin，Anhui on May 24，2013 recorded by seismometers on the ground and at different depths in the borehole underground of the Donghai station（a） Waveforms recorded by seismograph L0 on the ground；
（b） Waveforms after medium-pass filtering to the waveforms in Fig.（a）；
（c） Waveforms recorded by seismograph L1 at 544.5 m depth in the borehole；
（d） Waveforms recorded by seismograph L3 at 2 545.5 m depth in the borehole
3.3 井下地震仪微小地震观测能力探索
微小地震是深井地震学的重要观测研究对象.图4给出了东海站深井地震仪L3观测到的三个微小地震事件的波形，其中：图4a 为 2013年2月4日江苏徐圩台附近（34.2°N，119.3°E）ML−1.3地震（震中距为 64 km）的波形，其信噪比为 55.08 dB；
图4b为 2013年 1月21 日江苏灌云（34.2°N，119.4°E）ML−0.5 地震（震中距为 71 km）的波形，其信噪比为 52.84 dB；
图4c为2013年6月12日江苏涟水ML0.3地震（震中距为112 km）的波形，其信噪比为67.13 dB.从图中均可见清晰的P波和S波波形，无明显的噪声扰动；
这两个负震级地震波形图中的波形虽然存在不同程度的噪声，但都观测到了可供分析的清楚的P波初动和S波波形；
在背景噪声水平比较大的东海深井地震观测站能够高质量地观测到周围0震级或者负震级的微小地震.因此，东海深井地震观测站延展了高噪声地区观测微小地震的能力，可促进深井地震学的发展.目前进行中的深井长期观测可以进一步探索东海观测站观测微小地震下限的能力.
图4 东海观测站深井地震仪L3观测的微小地震三分量波形Fig.4 Three-component waveforms of the earthquakes recorded by seismograph L3 in the borehole underground of the Donghai station
统计分析地面和井下不同深度地震仪记录的地震波形信噪比，是探索避免地面噪声干扰、获取零噪声波形的深井地震学观测研究的新途径.本文分析了2012年上半年期间L0，L1，L2和L3四套地震仪记录到的116个地震事件的地震波形信噪比，如图5所示.可见，深井下不同深度地震仪观测到的地震波形信噪比大多在40 dB以上变化，深井地震仪L1，L2和L3的地震波形信噪比的平均值分别为65.33，69.20和74.15 dB，相较地面地震仪L0，其观测的近震地震波形的平均信噪比有明显的提高.图5中不同深度的地震仪观测到的波形信噪比变化范围虽然不同，但是井下三处不同深度地震仪的波形信噪比曲线的变化趋势基本一致，仅变化幅度大小有别.深井中三组地震仪信噪比平均值以及观测到优良地震波形的数量也均随仪器安装深度的增加而增加，这意味着地震噪声源主要来自地面，表明井下地震观测具有较好的抗地面噪声干扰的能力，是在地面高噪声环境中实现高精度观测微小地震活动的可行途径.因此深井中地震观测对避免地面噪声具有显著的积极效应.如图5中的观测结果所示，在东海观测站安装的深度大于1 500 m的地震仪L2和L3观测到的地震波形平均信噪比都接近或大于70 dB，观测结果比较理想.
图5 东海站地面L0与深井下L1，L2和L3地震仪观测地震信噪比计算结果Fig.5 Signal-to-noise ratio curves of seismic waveforms observed by seismometers L0 on the surface and L1，L2 and L3 in the borehole underground in Donghai station
本文分析的地震均为中小地震事件，震中距均在1 000 km以内，其优良地震波形几乎全部是由500 m以下的深井地震仪所记录，其中深度为2 545.5 m的地震仪L3记录到的优良地震波形最多.图5中地面安装的地震仪L0所记录到的1 000 km之内的中小事件的地震波形信噪比大多不足40 dB，平均值也较低，对于2012年3月16日江苏新沂ML1.1地震（震中距为34 km），地震仪L0观测到的地震波形的信噪比仅为11.27 dB.此外，图5中地面地震仪L0记录的地震波形信噪比变化与深井地震仪L1，L2和L3三者的信噪比曲线变化形态差别也很大.虽然地面宽频地震仪可以记录到远处大震的波形，但是在避免地面噪声干扰能力方面远不如深井地震仪.
本文分析了东海地壳活动观测研究站在地面和井下不同深度安装的四套地震仪所记录波形的特征及其信噪比随深度的变化，以探索深井地震仪观测微小地震的能力.由于东海深井观测站周围常出现强烈的非地震干扰，地面地震仪记录中无法识别微小地震的波形，而井下地震仪却可以记录到良好的地震波形，深井地震观测是创新深井地震学发展的必然趋势（Oyeet al，2004）.东海观测站井下地震仪可以记录到零震级或负震级微小地震的波形，记录中也可检测到明显的P波和S波震相.相对于地表观测，井下地震观测延拓了高噪声背景区域观测微小地震的下限，不仅提高了高噪声背景区域地震活动性的研究水平，而且为地壳和上地幔结构研究提供了重要资料（沈伟森等，2010；
王芳等，2017）.深井地震台网观测可以有效地避开地面噪声，接收更多高精度的地球内部信息，因此可以促进郯庐断裂及其周围地区的地壳活动监测研究.
东海深井地震观测获得了高保真波形.观测期间不同深度井下地震仪记录到的三分量P波和S波波形都清晰，背景噪声小，信噪比高，而且相对应的三分量波形彼此相似，这意味着深井地震观测受地面噪声干扰以及场地效应的影响小，波形保真度高.地面地震仪记录到的地震波形经常出现噪声扰动，甚至有的分量波形畸变，振幅异常，且与深井地震仪记录的相应三分量波形相似性差，这可能是受到地面噪声干扰或场地衰减的影响所致（罗诚等，2018）.因此，地面及深井下记录波形可用于场地效应等研究，与地面地震观测记录相比，基于井下地震仪记录的高保真波形所获得的介质运动学与动力学参数能更好地反映基岩的物理特征（王俊国等，1988）.深井观测地震波形为震源研究和震源谱分析提供了更直接的约束，这也是深井地震学的优势之一.
本文地震波形信噪比的统计结果表明，深井地震仪观测波形的信噪比明显高于地面地震仪的信噪比，且信噪比平均值随地震仪安装深度的增大而增加.观测结果进一步表明，深井地震仪观测可有效地获取高信噪比、高保真度的观测波形.深井中1559.5 m深处地震仪观测波形的平均信噪比为69.2dB，接近70 dB，2545.5 m深处观测波形的信噪比达到74.15 dB，二者都基本达到了良好地震波形的保真度和信噪比值.这表明，1500 m深处进行的深井地震观测基本可以避开地面噪声的干扰.深井地震学观测研究将有力地促进地球科学的基础理论研究和实际应用.
此外，观测资料显示，对于少数事件而言，井下较浅地震仪记录的地震波形信噪比高于深处的地震波形信噪比，这可能暗示了东海地区及周边地区岩石的非线性效应.地面宽频地震仪L0有时能够记录到可识别的大的远震波形，显示了在环境高噪声背景地区观测远处大地震的优势；
井下4.5 Hz主频地震仪记录到的远震的P波或S波初动清晰.因此可以考虑进行深井与地面地震仪联合观测研究，这将会取得更好的观测结果.
十分感谢三位审稿人对本文提出的宝贵意见和建议.
猜你喜欢深井信噪比波形两种64排GE CT冠脉成像信噪比与剂量对比分析研究现代仪器与医疗(2022年1期)2022-04-19基于Halbach阵列磁钢的PMSM气隙磁密波形优化防爆电机(2020年4期)2020-12-14绳短不能汲深井——保持“本领恐慌”当代陕西(2020年24期)2020-02-01赵楼煤矿：千米深井开采进入“智时代”山东煤炭科技(2019年12期)2019-12-27基于深度学习的无人机数据链信噪比估计算法北京航空航天大学学报(2019年9期)2019-10-26用于SAR与通信一体化系统的滤波器组多载波波形雷达学报(2018年5期)2018-12-05低信噪比下基于Hough变换的前视阵列SAR稀疏三维成像雷达学报(2017年3期)2018-01-19集水明排与深井井点联合降水工程应用中国港湾建设(2017年11期)2017-12-19基于ARM的任意波形电源设计通信电源技术(2016年3期)2016-03-26深井高应力软岩巷道支护参数优化铜业工程(2015年4期)2015-12-29