禁核试核查国际监测系统中的USRK及KSRS台阵对朝鲜地下核试验场的地震信号检测能力评估

薛方正, 徐雄, 靳平, 卢娜, 王红春

薛方正, 徐雄, 靳平, 卢娜, 王红春. 2018: 禁核试核查国际监测系统中的USRK及KSRS台阵对朝鲜地下核试验场的地震信号检测能力评估. 地震学报, 40(4): 448-460. DOI: 10.11939/jass.20170223
引用本文: 薛方正, 徐雄, 靳平, 卢娜, 王红春. 2018: 禁核试核查国际监测系统中的USRK及KSRS台阵对朝鲜地下核试验场的地震信号检测能力评估. 地震学报, 40(4): 448-460. DOI: 10.11939/jass.20170223
Xue Fangzheng, Xu Xiong, Jin Ping, Lu Na, Wang Hongchun. 2018: Evaluation of capability for seismic signal detection of the USRK and KSRS arrays of IMS to the North Korean underground nuclear test site. Acta Seismologica Sinica, 40(4): 448-460. DOI: 10.11939/jass.20170223
Citation: Xue Fangzheng, Xu Xiong, Jin Ping, Lu Na, Wang Hongchun. 2018: Evaluation of capability for seismic signal detection of the USRK and KSRS arrays of IMS to the North Korean underground nuclear test site. Acta Seismologica Sinica, 40(4): 448-460. DOI: 10.11939/jass.20170223

禁核试核查国际监测系统中的USRK及KSRS台阵对朝鲜地下核试验场的地震信号检测能力评估

基金项目: 西北核技术研究所长线课题(14101503)资助
详细信息
    通讯作者:

    王红春: e-mail: wanghongchun@nint.ac.cn

  • 中图分类号: P315.63

Evaluation of capability for seismic signal detection of the USRK and KSRS arrays of IMS to the North Korean underground nuclear test site

  • 摘要: 根据禁核试核查国际监测系统中的USRK及KSRS台阵的实际记录数据,统计建立了这两个台阵的本底噪声水平分布概率模型。基于地下爆炸的源频谱模型及比例关系,对2009年5月25日朝鲜核试验在这两个地震台阵上的实际记录信号幅值进行了比例折算,得到了更小震级的“朝鲜地下核试验”在这两个地震台阵上的多频带信号幅值。在3倍信噪比检测条件下,利用概率模型法分别得到了这两个地震台阵对朝鲜地下核试验场小震级地下核试验的多频带信号检测能力。结果表明,震中距在400 km左右的两台阵对朝鲜地下核试验地震信号的检测阈值均随着频率的升高而逐渐降低,即高频段的检测能力强于低频段。这两个地震台阵的最高信号检测能力均出现在6.0—9.0 Hz频带,90%检测概率条件下,USRK和KSRS在该频带的检测阈值分别为mb1.5和mb2.2。
    Abstract: The noise level distribution probability models of the USRK and KSRS arrays are derived from the observed data of the two arrays which belongs to the International Monitoring System (IMS) Arrays. Based on the scaling model of source spectrum of underground explosion, amplitudes of seismic wave signal excited by the second North Korean nuclear test on May 25, 2009 and recorded by these two arrays are scaled to amplitudes of smaller underground explosions with lower magnitude. Thus, the amplitudes of the North Korean nuclear tests with lower magnitude are obtained. By the probabilistic model method, the detection capabilities of these two arrays for seismic signal from Korean nuclear tests with different magnitudes are given in the multi-band on the condition that signal-to-noise ratio is greater than three. The result indicates that the detection threshold is reduced with the frequency increasing at these two arrays on the condition epicentral distance within 400−500 km. Capability in the high frequency is much stronger than that in the low frequency. In the band of 6.0−9.0 Hz, the detection thresholds of the USRK and KSRS arrays aremb1.5 and mb2.2 respectively in the 90% detection probability.
  • 井下地震观测可以减少地表的噪声干扰影响。井底背景噪声的频谱幅度约为地表的数百分之一,甚至千分之一(Takahashi,Hamada,1975),可有效地提高观测波形的信噪比(徐纪人,赵志新,2009)。井下地震仪可以在高噪声背景区域观测到更多的微小近震,近距离观测的地震波更多地保留了其中的高频成分,而高频成分正是解析地球内部物理、地壳和上地幔结构不可缺少的信息,对研究断裂活动、理解地质过程等都起着极为重要的作用(Aster,Shearer,1991Yamauchi et al,2005)。井下观测地球物理状态的变化,有利于地震活动性的分析及震源机制和地震预报研究(Oye et al,2004)。目前井下地球物理观测研究已逐渐发展成为多学科的前沿性科学研究课题(Baisch et al,2002),其成果推进了现代地震学、地球物理学以及钻孔板块构造学等学科的发展(Zoback et al,2011)。由于强地面运动受到厚沉积物的显著影响,井下观测也是研究沉积平原地区强地面运动的有效途径(Zhao et al,2004)。我国井下地震观测研究始于20世纪70年代,至上世纪80年代中期已有百余个井下地震观测台,遍布全国(修济刚,1988)。

    近年来井下地震观测台站在数量及深度上都得到迅速发展,观测多使用国产地震仪,且观测质量优良。在我国1 000多个数字化测震台站中,井下观测台约200个(张明等,2019)。2001年建成的首都圈测震台网的107个台站中,58个观测台站配备有井下地震仪(张尉等,2009),已发展成为先进的井下地震观测台网,促进了我国地震事业的发展。图1给出了我国大陆地区表1中部分井下观测台站的分布。其分布与地震构造带关系密切,主要分布在我国东部沉积层覆盖比较厚的地区,目的在于监控地震活动及开展地震研究工作。增设井下观测台站,形成科学的全方位观测台网,可以提高地震观测研究能力。为了监视大断裂的构造活动和地震活动,我国在西部一些地区也增设了井下地震观测研究台站。

    图  1  中国大陆地区部分井下地震台站分布简图
    Figure  1.  A distribution diagram of some underground seismic stations in Chinese mainland
    表  1  我国部分井下观测台站基本情况
    Table  1.  Basic information of some underground observation stations in China
    台站名 井深/m 地震仪型 参考文献 台站名 井深/m 地震仪型 参考文献
    首都圈台网 150—480 宽、甚频 短周期 刘渊源等(2 011 上海东滩 421 FSS-3DBH 裴晓等(2 012
    北京市台网 宽、甚频 短周期 兰从欣等(2 005 上海金泽 305 FSS-3DBH
    首都圈白家疃 257 韦士忠和李玉萍(1 990 上海上戏 370 FSS-3DBH
    首都圈文安 266 上海南汇 280 JDF-2
    首都圈东三旗 250 上海大新中学 375 FSS-3DBH
    首都圈大兴 110 上海竹园 317 JDF-1
    首都圈雄县 358 上海虹桥 651 768
    首都圈龙门庄 447 上海八角厅 780 JDF-1
    吉林松原 384 FSS-3DBH 韦庆海等(2 015 上海张江 350 KS-2000M 裴晓等(2 013
    吉林松原 243 FSS-3DBH 陈闯等(2 022 上海崇明 463 李伟等(2 013
    大庆新台 708 JD-2 韦庆海等(2 015 江苏宝应 460 CMG-3TB 仇中阳等(2 014
    河北赵县 260 BBVS-60 郑德高等(2 018 江苏高邮 440 CMG-3TB
    河北唐海 480 短周期 郑德高等(2 018 江苏淮安 315 JDF-2
    河北涿县 320 768 李彦林和郑淑兰(1 989 江苏涟水 400 CMG-3TB
    河北邯郸 400 JD-2 张新东(2 002 江苏射阳 380 CMG-3TB
    河北肥乡 400 JD-2 江苏盐城 445 CMG-3TB
    河北临漳 440 JD-2 江苏金湖 447 GL-S60B 宫杰等(2 019
    天津静海 371 768 赵惠君等(1 991 江苏滨海 470 GL-S60B
    天津芦台 276 768 江苏丹阳 175 GL-S60B
    天津武清 450 768 江苏响水 410 GL-S60B
    内蒙赤峰 90 GL-S120B 郭延杰等(2 020 江苏高邮 458 GL-S60B
    甘肃天水 337 短周期 蔡耐芳(1 990 江苏建湖 443 GL-S60B
    新疆喀什 283 GL-S60B 赵瑞胜等(2 021 江苏启东 410 GL-S60B
    山西太原 500 JD-2 张少泉等(1 988 江苏盐城 436 GL-S60B
    陕西定边 300 BBSV-60BH 李少睿等(2 016 江苏泰兴 425 GL-S60B
    宁夏灵武 248 JDF-2 江苏东台 458 GL-S60B
    宁夏陶乐 245 JDF-2 江苏溧阳 83 CMG-DM24 mk3 胡米东(2 014
    河南安阳 393 FSS-3DBH 江苏盐城 445 CMG-DM24 mk3
    河南清丰 308 FSS-3DBH 江苏南通 105 CMG-DM24 mk3
    四川泸州 95 CMG-3TB 江苏大丰 366 短周期 徐元耀(1 994
    山东荷泽 370 768 周焕鹏(1 986 江苏淮阴 325 井下摆
    安徽六安 126 GL-S60B 石英杰等(2 021 江苏海安 420 JD-2
    安徽霍邱 150 GL-S60B 石英杰等(2 021 云南昆明 2 02 GL-S60B 李雷等(2 018
    浙江景宁 68 FSS-3DBH 张明等(2 019 云南昆明 452 JD-2 修济刚(1 988
    浙江北仑 86 FSS-3DBH 云南大寨 375 井下仪 王芳等(2 017
    浙江南麂岛 110 GL-60DBH 广东汕头 2 00 TBG-60B 郭德顺等(2 014
    上海普陀 564 768 叶世元和柳国华(1 987 江苏东海 4050 宽频 短周期 Xu等(2 016
    上海海运 600 768 叶世元和柳国华(1 987
    下载: 导出CSV 
    | 显示表格

    我国井下台网的观测技术已日趋完善成熟。在井下密封、井下定向、方位补偿、电缆应力解除等观测技术领域也取得了令人满意的进步。而适合不同井孔孔径的锁固装置的成功研制及利用反馈控制技术调节井下长周期地震仪的零位漂移、振幅效应等技术,可保持井下地震仪长期稳定工作(高龙生,1986)。通过编码遥控装置调整井下的地震仪摆体的夹紧、放松状态以及标定等多种功能,可改进井下地震仪的工作状态(陈宏水,付迺珍,1987)。规范井下地震仪安装过程以保证其良好的观测能力(胡米东,2014);规范安装观测房和井口之间的防雷设施及线路布设,使井下地震台站稳定可靠地运行(陈吉锋等,2014)。根据条件选择落底安装井下地震仪,可以提高地震仪的观测能力(宫杰等,2019)。利用波形相关方法,分析地表与井下地震仪方位角检测以及一致性对比测试,使井下地震仪方位角检测精度优于4° (李少睿等,2016)。在缆线应力释放完成后固定井下地震仪,同地表地震仪进行数据相关性分析,以确定方位角的可靠性(张明等,2019)。

    井下观测多使用我国研制的先进的短周期、宽频和甚宽频带地震观测仪器(刘瑞丰等,2008)。井下观测资料的分析方法也在不断精益求精。对比分析各种地震观测参数,研究人员提出了高效、可靠的识别与提取井下微地震信号的方法(李稳等,2016)。利用我国井下台网已经可以进行全频域、高精度、长期稳定的地震监测研究,井下台网也成为获取可靠资料的保障,为地震科学的发展提供了更坚实的基础。表1列出了本文涉及的井下观测台站的深度、观测地震仪的频段特性以及观测成果来源。

    上述先进的观测仪器和井下观测技术提高了地震观测能力和科学研究水平。地面背景噪声常导致地震波形很难被地面地震台识别,采用井下摆观测方式就可以有效抑制各种人为活动噪声的干扰。宽频地震仪安装在地下10 m的中微风化层基岩上,可有效减小噪声(郭德顺等,2014)。研究平原地面、井下与山区地震的纵、横波的运动学和动力学特征得知,井下记录比较接近设置于基岩上的地震台的记录。井下记录的地震波形的噪声功率均值较小(张阿瑶等,2018),可以更精确地记录到地脉动(裴晓等,2013)。上海张江井下台的背景噪声的均方根比地表台站要小1个数量级,而且近震、远震波形记录都很清晰。井下摆的放大倍数亦可达地面摆的十倍多(李彦林,郑淑兰,1989)。井下观测的数据精度也能提高1—2个数量级,并且井下观测系统的有效动态范围超出地面系统约10%—30%,能记录到的地震动信号的范围更宽(叶世元等,1988)。井下90 m深地震仪也能观测到比山洞中深处花岗岩台址上地震仪更远的同等大小的地震(郭延杰等,2020)。

    东海台阵深井地震仪可以检测几十至几百千米外的微小地震的清晰波形,而地表记录却识别不出这些地震的迹象(表2)。图2示出了东海台阵的日常监测屏上不同深度地震仪显示的2019年9月22日山东烟台ML3.1地震纪录图。图中地表记录识别不出地震迹象,井下400 m深处记录可看出明显地震波形。而且,地下3 km深处记录也可看出清晰的地震波形,表明监控记录可靠。井下观测填补了在高噪声区域获得高精度地震波形的空白,在提高地震监测分辨能力及观测地震波精度等方面具有明显优势(徐元耀,1994徐纪人,赵志新,2006)。

    表  2  江苏东海垂直地震台阵地面与井下观测能力比较
    Table  2.  Comparison of observation ability between surface and underground of vertical seismic array in Jiangsu Donghai
    ML 震中距
    Δ/km
    地面观测宽频
    带地震仪
    井下观测
    井深/m 仪器 观测效果
    −1.3 64 无法识别 2545 短周期 图像清晰
    −0.5 71 无法识别 2545 短周期 图像清晰
    0.8 97 无法识别 3500 短周期 图像清晰
    3.1 492 无法识别 400 宽频 图像清晰
    下载: 导出CSV 
    | 显示表格
    图  2  东海垂直地震台阵的地面及井下仪器记录的2019年9月22日烟台ML3.1地震
    G0:地表宽频地震仪;G1:400 m深宽频地震仪;G3:3 km深短周期地震仪
    Figure  2.  The September 22,2019,Yantai ML3.1 earthquake recorded by surface and underground instruments of Donghai
    vertical seismic array G0:Surface broadband seismograph;G1:400 m deep broadband seismograph; G3:3 km deep short period seismometer

    井下地震观测台网记录到的地震波形,噪声小且精度高,这些可靠的高质量资料较精确地反映了地震波路径及震源信息。在远离地面干扰环境中,长期观测地球内部的信息,所得结果客观真实。用其高分辨率记录测定地震参数,研究震源过程、地震成像和地球内部结构,可以改善研究成果的精度,有利于创新地震波理论研究,且目前已取得了可喜的新成果。

    井下观测研究避免了地表噪声和场地效应的影响,可全方位布局地震观测台,这样增用井下观测资料得到的地震定位结果更加准确(仇中阳等,2014),进而可推进对地震序列活动及震源位置的研究,以及断层构造运动等的研究。分析井下摆与地面摆测定的地震参数差异成因,并给出两者之差的适当校正值(朱雅山,李永勤,1992),可提高地震观测成果的客观性与可靠性。

    井下观测地震波形噪声小,有利于高精度测定地震P波和S波到时、振幅、周期等参数(赵惠君等,1991),推进了地震学的研究。上海南汇台用统计方法分析了井下记录地震尾波衰减与震级之间的关系,求得了垂直向地震波持续时间与对数震级Mv在近震范围内(Δ<5°)呈线性关系(李慧民等,1992)。用井下地震仪记录的最大振幅测定近震震级,误差分析表明,结果合理可靠(章纯等,1992)。吉林地区台站由井下地震仪得到的近震震级比地表地震仪的小0.12 (陈闯等,2022)。河北涿县井深350 m的短周期摆与地面同型摆记录的地震波形对比研究表明,地面S波两水平向最大地动位移平均为井下摆的4.7倍。地面摆垂直向最大地动位移平均为井下摆的2.7倍。因此,井下地震摆观测所得震级比地面摆平均低0.7。如采用波形振动持续时间测定震级,井下记录可得到与地面记录较为一致的结果。涿县地震台P波首波到时,地面地震仪摆要比井下摆滞后0.3 s (张寿康等,1986)。大庆井下地震台网不同深度各子台观测结果表明,近震震级偏差和校正与井下摆深度和介质性质等有关(袁卫红等,2012)。

    新疆喀什麦盖提台站井下宽频地震仪记录的近、远震波形比地面的更为清晰。井下地震仪得到的近震和远震震级分别比地表地震仪得到的震级小0.518和0.025,且功率谱也相对较低(赵瑞胜等,2021)。昆明基准地震台山洞内地面上的宽频带地震仪观测结果ML比井下结果平均大0.18 (李雷等,2018)。宽频带地震仪井下记录的地震波的运动学参数、地震波走时比地面的小,统计得出地面地震台的动力学参数地动振幅大约是井下振幅的10倍(王俊国等,1988),与上述短周期地震仪研究结果类似。地震矩和矩震级观测也表明井下震级较地面要小(仇中阳等,2014)。这可能是由于地面记录存在较大的地震波振幅效应,某些情形下地面振幅相对于井下振幅的放大因子可能大于1。井下观测结果是相对客观而科学的。图3显示出部分井下与地面观测的同一地震事件的震级比较。图中红线为横轴与纵轴等值参考线,横轴和纵轴数值分别表示地面地震仪与井下地震仪观测到同一事件的震级。可看出井下震级明显小于地面震级,特别是对于M4.0以下的小地震(毛华锋等,2014仇中阳等,2014赵瑞胜等,2021)。目前井下观测资料在分析地震参数、进行地震速报及监测地震活动等方面日益发挥出重要作用(修济刚,1988),推进了深井地震学研究的发展。

    图  3  同一地震井下与地面观测的震级比较
    Figure  3.  Comparison of seismic magnitude between underground and surface observation for the same earthquake

    低噪音的井下地震观测资料已用于地球内部结构研究,有益于建立高精度三维地震波速度模型(Lay et al,2020)。地表和井下台站波形资料的联合利用提高了监控区域地震波速研究精度,揭示了京津唐地区地震波速度结构明显的横向不均匀性。并且这些资料也被用于探索地震成因及地震预报,发现许多地震均发生在P波低速异常体与高速异常体的交界地带(于湘伟等,2003)。唐山震区的三次强震均发生在低速体与高速体之间(齐诚等,2006)。地震波速时间的变化显示2006年文安MS5.1大地震前一年,文安和唐山地区波速比明显出现了可信的前兆性降低过程(王林瑛等,2008)。

    井下地震观测记录到了体波的地表反射震相。这一发现促进了地震学观测的发展。井下地震仪记录的地震波震相较之地面记录更为丰富(乌统昱,赵惠君,1992)。大庆台井下708 m地震仪记录到2005年7月25日林甸M5.1地震和2003年6月26日松原M2.3地震的地表反射震相。用井下接收的地面反射震相PP,SP等可以确定井下观测站与地面之间的速度值和界面埋深(林云松,1989)。静海台井下地震仪记录到周围微震波形的4个主要震相,即P,S,地面反射波PP和SS,其中SS震相的水平分量通常较强。松原台384 m深处井下摆记录到许多地方震的地表反射波(韦庆海等,2015)。用地表反射波资料求得松原台地表到井下地震仪之间的P波速度约为2 km/s (陈闯等,2022)。新疆喀什台站用井下资料求得麦盖提下方浅层的平均S波速度为0.59 km/s (赵瑞胜等,2021)。由此结果构成的地表沉积层速度模型可有效改善地震定位精度。利用反卷积方法识别近震入射波与地表反射波,可建立高精度浅层地震波速度模型(王芳等,2017Riga et al,2019)。

    井下摆反射波形分析是一种测量深达数百米S波速度的有效方法。分析首都圈台网中58个井下摆记录的100多次M≥3.0地震的直达波及对应地表反射波的波形可知,这些台站上方厚约300 m浅地表土层P波和S波平均速度分别约为1.6—2.0 km/s和0.34—0.48 km/s,波速比约为4.0—5.3。该区域P波速度和波速比均无显著横向变化。清晰的震相图也表明,土层没有明显的速度断层(刘渊源等,2011)。首都圈地表S波平均速度梯度为0.8 m/(s·m),浅部100 m的平均S波速小于300 m/s,深度500 m处S波速度为800 m/s (沈伟森等,2010)。井下摆近震记录研究表明,华北盆地近地表P波平均速度约为1.98 km/s,S波平均速度约为0.46 km/s (朱音杰等,2020)。上述地震波速度成果对计算地球内部精细结构、地震波传播在地球内部的衰减特征以及资源评估均具有重要意义,同时也可为研究场地效应提供直接证据(罗诚等,2018)。

    研究地震记录的背景噪声功率谱可知,井下地震台在全频范围内记录数据质量明显优于地表地震仪(石英杰等,2021)。井下台网高质量的资料为研究震源和地球内部物理学奠定了坚实的基础,取得了大量可信赖的高精度成果。

    利用井下和地面观测资料求得了京津唐地区1 132次ML1.7—6.2地震的精确震源参数(于湘伟等,2003),促进了震源过程研究。用最小二乘法拟合北京地区震源参数与震级的定量表达式,证实拟合所使用资料的台网规模对震源参数结果影响程度的差异。场地响应对震源参数的计算结果影响也比较大(兰从欣等,2005)。江苏部分地震台井下观测波形拐角频率比地面观测的小0.397 Hz,井下观测的地震矩和矩震级分别比地面的小8.642×1014 N·m和0.32,断层错距和应力降分别比地面的小2.268 cm和2.033 MPa (毛华锋等,2014)。仇中阳等(2014)的结果显示,井下记录地震波拐角频率较地面约小0.28 Hz;地震矩和矩震级数值较地面分别小3.104×1014 N·m和0.32。井下测定的地震波谱高频较弱,拐角频率较低,地震矩、应力降和平均位错都比地面基岩测定的结果偏低,而震源等效半径则偏高,井下记录高频成份没有地面的丰富(韦士忠,李玉萍,1990)。接收函数正演计算和井下理论波形研究发现,井下地震仪记录P波波形垂向分量的频谱极小值随深度的增加而减小(郑德高等,2018)。首都圈西北部地区快剪切波平均偏振方向为NE方向,它暗示了该区域的主压应力水平方向;慢剪切波时间延迟急剧的梯度变化可能与地壳深部的温度变化有关(吴晶等,2007)。上述结果对于探索地震发生机制有参考意义。加之钻孔井下记录受场地效应的影响比地表地震仪小,可为研究震源提供更直接的约束(Shearer,Abercrombie,2021),其成果丰富了现代地震学研究宝库。

    强地面运动的振幅和持续时间受到厚沉积物速度结构的显著影响。研究地震波传播在沉积物中的非线性效应,以及近地表地震波衰减和场地效应对于精确预测强地面运动至关重要(Stephenson et al,2005)。

    将短周期井下摆与地面仪记录对比可知,地表土层对地震波有放大作用。北京西部地区地表基岩台的场地响应总体上比东部平原区井下台的响应大,其低频端场地响应变化平稳。东部平原区井下台的场地响应也与频率有关(兰从欣等,2005)。品质因子Q值与频率的定量关系显示,台站场地响应在低频段幅度变化亦保持稳定(蔡耐芳,1990卞真付等,2005)。将地面和井下观测场地效应资料用于波场理论计算,可减少波场解的不确定性,为工程地震提供更可靠的依据(张少泉等,1992)。

    井下强震动观测台阵中加速度仪记录是研究强地面运动的重要方法(Hollender et al,2023)。地表与井下实测数据可有效地用于研究各种土层的地震响应(丁海平等,2014)。进一步反演结果指出,土层的S波和P波速度及厚度是影响场地响应的水平与竖向谱比的重要因素(李小军等,2020)。由井下基岩记录的地震波作为输入研究土层地震动响应可知,土层对地震波的放大效应与波形主频率有关,放大倍数最高可达约20倍。该结论对软土地区的烈度区划、建筑的抗震设计有现实意义(徐永林等,2002)。

    深井地震垂直台阵观测研究有益于地壳活动性的监测及地震成因的研究,开创了高精度地震学研究的新途径(Ma,2021)。台湾车笼埔断层带上的井下垂直地震台阵中,多台波形水平分量之间的尾波互相关函数的研究结果显示,该区地震各向异性呈现出季节性变化(Hung et al,2022)。

    江苏东海大陆深孔地壳活动国家野外科学观测研究站在国际大陆科学钻探项目(ICDP)之一的中国大陆科学钻探的5000 m深井中,建设了我国大陆第一个深井垂直地震台阵,即东海大陆深井垂直地震台阵(简称东海垂直地震台阵)。除设立地面观测地震仪之外,还在井下不同深度设置了4层地震观测仪。在相距深井约500 m的另一个450 m井中进行了综合地球物理观测。2008年汶川MS8.0地震后,在汶川地震区跨断层建立了6口200—3000 m不同深度的井下观测台,深入研究大地震发生机制,监测断层活动。

    东海大陆深井垂直地震台阵及其周围处于高噪声背景区域,致使地表地震仪记录中难于识别ML0.8以下地方震波形。深井地震台阵能够观测到比地面地震仪更多的零级或负震级微小地震的清晰波形,观测资料清晰可靠。井下不同深度地震台观测的地震序列的波形信噪比随时间变化的趋势基本一致。其地震波形基本达到了70 dB以上的高信噪比。在高噪声背景区域获得了高保真度的微小地震波形(徐纪人等,2022),提高了监测地壳活动和地震活动的能力。不同深度处的记录波形的信噪比可为井下观测台站建设提供科学参考依据。低噪声高质量的地震波形研究成果可以加深对地震运动学和动力学的理解;建立精细波速模型,有助于全面清晰地揭示地球内部结构,创新地震波传播理论研究(Ma,2021)。分析低噪声的微小地震至大地震的波相特性,可为震源和震源谱的研究提供直接约束,更深刻地解释震源破裂过程和机制。利用深井地震波形研究场地效应可以提高强地面运动预测水平。未来深井垂直地震台阵观测研究或可在高精度地震学研究中发挥巨大作用。

    本文研究表明,井下观测可以在高噪声背景区域近距离观测地震活动,应科学地布局观测台站,以取得高质量地震资料,进而推进地震科学的发展。用井下资料研究波速空间分布的非均匀性,可以缩小因地面台站台基的不同而导致的各台站间的P,S波的平均速度差异,提高所建立的区域波速模型的可信度和科学性。井下观测与地面观测的震级、矩震级比地面基岩测定的结果偏低(图3),其差异可能与井下仪器上层介质的波形非线性增幅效应与波的频率有关(Zhao et al,2004)。井下观测拐角频率低于地表可能与井下仪器上层介质对不同频率波分量吸收与放大有关(周焕鹏,1986)。场地响应对震源谱参数影响也较大,在低频域,绝大多数井下台的场地响应都大于1,高频端,绝大多数井下台的场地响应都小于1(兰从欣等,2005),这些因素都可能导致井下台站观测的震级、拐角频率都低于地面观测结果。上述差异的成因也具有多重复杂性(Abercrombie,1997),仍是今后科学探索的重要课题。井下地震观测开创了我国在地下深处观测研究地球内部的新途径,提高了地震观测能力,为地球科学研究作出了巨大贡献;对资源评估、地震预报、防震减灾也具有重大的现实意义(徐纪人等,2022)。本文未能引用全国的数百个井下观测台的全部累累成果,这些论文和成果可以在相关的文献中见到。

    深井地震观测研究是现今国内外的前沿课题,提高其观测研究水平可实现高精度地震学研究,而高精度地球物理学正是创新未来的必然发展轨迹(Zeng,Yang,2021)。我国还将要沿安宁河断裂带建设12个井深300 m和两个井深1000 m的综合地震观测研究井(王晓蕾等,2021);目前已经建成8个井深300 m和两个井深1000 m的综合地震观测台阵,其中湖南2000 m深井宽频地震观测站已取得重大成果。东海大陆深孔地壳活动国家野外科学观测研究站将继续增强汶川及东海地区井下综合地球物理长期观测研究工作。随着井下台网观测研究工作的发展,我国的高精度科学研究成果将更上一层楼。

  • 图  1   一年内(a)和一天内(b)的噪声样本分布

    Figure  1.   The distribution of noise samples in a year (a) and a day (b)

    图  2   USRK台阵(a)和KSRS台阵(b)在各频带的累积概率分布曲线

    Figure  2.   Cumulative probability distribution curves of the KSRS array (a) and the KSRS array (b) in the multi-band

    图  3   Pn波位移谱残差随源频谱拐角频率fc的变化曲线

    Figure  3.   Residual variation of displacement spectral of Pn wave with corner frequency fc of source frequency spectrum

    图  4   理论Pn波频谱与实际观测频谱的比较

    Figure  4.   The comparison of spectrum between theoretical Pn waveform and observed waveform

    图  5   部分小当量爆炸的源频谱(a)以及相对于参考事件的谱比值(b)

    Figure  5.   Source spectrum of explosion events (a) and spectral ratio to the reference event (b) with small yield

    图  6   USRK (左)和KSRS (右)台阵的多频带波形数据

    Figure  6.   Waveforms in the multi-band for the USRK (left) and KSRS (right) arrays

    图  7   USRK台阵(a)和KSRK台阵(b)各频带检测概率

    Figure  7.   Detection probability of the USRK array (a) and the KSRS array (b) in the multi-band

    图  8   不同震级事件USRK(a)和KSRS(b)台阵的检测阈值

    Figure  8.   Detection threshold with different magnitudes for USRK array (a) and KSRS array (b)

    图  9   2010年5月12日朝鲜核试验场小震事件的波形

    图中时间轴零点处为Pn波信号到时位置

    Figure  9.   Waveforms of the small earthquake on May 12,2010 from the North Korea nuclear test site

    The zero position of the time axis is the arrival time of Pn wave

    表  1   不同震级朝鲜地下核试验在USRK和KSRS台阵上不同频带Pn波信号的实测和预测幅值

    Table  1   Observed and predicted amplitudes of Pn waves with different frequencies at USRK and KSRS arrays generated by North Korea nuclear test with different magnitudes

    mb 台阵 信号幅值/(μm·s–1
    0.75—1.5 Hz 1.0—2.0 Hz 2.0—4.0 Hz 3.0—6.0 Hz 4.0—8.0 Hz 6.0—9.0 Hz
    4.5 USRK 1.37 2.68 6.66 4.41 2.59 2.25
    KSRS 1.05 1.86 5.16 4.24 8.42 27.64
    4.0 USRK 0.43 0.84 2.79 2.57 1.72 1.56
    KSRS 0.34 0.56 2.13 2.45 5.33 19.16
    3.0 USRK 0.04 0.08 0.28 0.35 0.35 0.49
    KSRS 0.03 0.05 0.22 0.33 1.09 5.99
    2.0 USRK 0.004 0.008 0.03 0.03 0.04 0.05
    KSRS 0.003 0.005 0.02 0.03 0.11 0.65
    下载: 导出CSV
  • 冯建刚, 蒋长胜, 韩立波, 陈继峰. 2012. 甘肃测震台网监测能力及地震目录完整性分析[J]. 地震学报, 34(5): 646-658.

    Feng J G, Jiang C S, Han L B, Chen J F. 2012. Analysis on the monitoring capability of seismic networks and completeness of earthquake catalogues in Gansu region[J]. Acta Seismologica Sinica, 34(5): 646-658 (in Chinese).

    何少林. 2003. 甘肃省地震台网监测能力分析[J]. 地震地磁观测与研究, 24(6): 103-108.

    He S L. 2003. The analysis of monitoring ability of seismic network of Gansu province[J]. Seismological and Geomagnetic Observation and Research, 24(6): 103-108 (in Chinese).

    李志海, 蒋长胜, 黄瑜, 王宝柱. 2011. 新疆地区地震目录最小完整性震级和台网科学布局研究[J]. 地震学报, 33(6): 763-775.

    Li Z H, Jiang C S, Huang Y, Wang B Z. 2011. Minimum magnitude of completeness for earthquake catalogue and scientific layout of seismic stations in Xinjiang region[J]. Acta Seismologica Sinica, 33(6): 763-775 (in Chinese).

    王燕, 刘俊民, 王海军, 唐恒专, 李靓, 唐伟, 刘哲函. 2011. 阈值监测技术在禁核试地震台网监测能力评估中的应用[J]. 核电子学与探测技术, 31(2): 163-166.

    Wang Y, Liu J M, Wang H J, Tang H Z, Li J, Tang W, Liu Z H. 2011. The application of threshold monitoring which is used for evaluating seismic network[J]. Nuclear Electronics & Detection Technology, 31(2): 163-166 (in Chinese).

    Bahavar M, North R. 2002. Estimation of background noise for international monitoring system seismic stations[J]. Pure Appl Geophys, 159(5): 911-944.

    CTBTO. 2017. Station profiles[EB/OL]. [2017–07–18]. http://www.ctbto.org/verification-regime/station-profiles/.

    Denny M D, Johnson L R. 1991. The explosion seismic source function: Models and scaling laws reviewed[G]//Explosion Source Phenomenology. Washington DC: American Geophysical Union: 1–24.

    Gerlach G S, Goforth T T, Seymour F F, Swanson J G. 1966. Estimates of the Detection Capability of Four VELA-UNIFORM Seismological Observatories[R]. Garland, Texas: Technical Report: 1–65.

    Koper K D, Herrmann R B, Benz H M. 2008. Overview of open seismic data from the north Korean event of 9 October 2006[J]. Seismol Res Lett, 79(2): 178-185.

    Kværna T, Ringdal F, Baadshaug U. 2007. North Korea’s nuclear test: The capability for seismic monitoring of the North Korean test site[J]. Seismol Res Lett, 78(5): 487-497.

    Rougier E P, Howard J K, Earl E B, Christopher R. 2011. Constraints on burial depth and yield of the 25 May 2009 North Korean test from hydrodynamic simulations in a granite medium[J]. Geophys Res Lett, 38(16): L16316.

    Sereno Thomas J Jr, Bratt S R, Bache T C. 1988. Simultaneous inversion of regional wave spectra for attenuation and seismic moment in Scandinavia[J]. J Geophys Res: Solid Earth, 93(B3): 2019-2035.

    Yang X, Lay T, Xie X B, Thorne M S. 2007. Geometric Spreading of Pn and Sn in a Spherical Earth Model[J]. Bull Seismol Soc Am, 97(6): 2053-2065.

    Zhang M, Wen L X. 2015. Seismological evidence for a low-yield nuclear test on 12 May 2010 in North Korea[J]. Seismol Res Lett, 86(1): 138-145.

    Zhao L F, Xie X B, Wang W M, Yao Z X. 2012. Yield Estimation of the 25 May 2009 North Korean Nuclear Explosion[J]. Bull Seismol Soc Am, 102(2): 467-478.

    Zhao L F, Xie X B, Tian B F, Chen Q F, Hao T Y, Yao Z X. 2015. Pn wave geometrical spreading and attenuation in Northeast China and the Korean Peninsula constrained by observations from North Korean nuclear explosions[J]. J Geophys Res: Solid Earth, 120(11): 7558-7571.

图(9)  /  表(1)
计量
  • 文章访问数:  1514
  • HTML全文浏览量:  689
  • PDF下载量:  40
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-12-20
  • 修回日期:  2018-03-26
  • 网络出版日期:  2018-06-21
  • 发布日期:  2018-06-30

目录

/

返回文章
返回