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. -
引言
全面禁止核试验条约组织（Comprehensive Nuclear-Test-Ban Treaty Organization，简写为CTBTO）为了核查世界各签约国遵守条约的情况，建立了以地震、次声、水声和放射性核素为主要监测手段的国际监测系统（International Monitoring System，简写为IMS），并在奥地利维也纳建立了接收、存储和处理监测系统实时监测数据的国际数据中心（International Data Center，简写为IDC）。其中地震监测作为地下核试验的一种有效远区监测技术手段，一直倍受CTBTO重视。截至2017年7月，已经完成了由50个主要台站和120个辅助地震台站所组成的庞大地震监测系统89%的核准与95%的建设（CTBTO，2017）。IMS的地震监测系统具有对全球1 kt以上当量的地下核试验的监测能力。由于国际监测系统的地震台站分布不均匀，对部分地区的实际监测阈值可能已经低于1 kt。地震台站的信号检测能力是地震台网事件监测能力的一个重要基础指标。地震台网的事件监测能力事实上主要是由组成地震台网的单个地震台站的信号检测能力所决定的，因此准确评估地震台站对重点地区地震事件，特别是地下爆炸地震信号的检测能力，对于进一步优化台网布局、提高对敏感地区可能地下核试验的监测能力以及快速响应能力具有重要的指导意义。
朝鲜是自IMS投入运行后，多次进行地下核试验的国家，因此IMS对该地区的信号监测能力也备受国内外核监测地震学家的关注。朝鲜第一次地下核试验后不久，Kværna等（2007）运用阈值监测法得到IMS台网在相应时段内对朝鲜地下核试验场地震信号的实时动态检测阈值及事件监测阈值。他们在计算检测阈值时以朝鲜第一次地下核试验作为校准事件，以4倍信号信噪比作为信号检测条件，以3个以上IMS地震台站上的信号信噪比同时达到4倍信号信噪比以上作为形成事件监测能力的条件。评估结果表明，在2006年朝鲜进行第一次地下核试验期间，IMS对朝鲜地下核试验场的实时信号检测能力大约为mb3.0，而相应的地震事件监测能力大约为mb3.8。不过Kværna在开展IMS监测系统对朝鲜地下核试验场的地震监测能力评估工作时，属于IMS的位于韩国原州的地震台阵KSRS并未向IDC实时传送数据，同时，位于俄罗斯乌苏里斯克的地震台阵USRK也并未建成。Kværna等利用长度约为3 min的KSRS台阵的监测数据进行分析后，认为增加KSRS台阵后，IMS系统对朝鲜核试验场的信号检测能力还可以下降0.5级，达到2.5级，而相应的事件监测能力基本上没有变化。Kværna等的研究结果还表明，如果将中国牡丹江的地震台站MDJ的记录数据与IMS的监测数据一起处理，则IMS对朝鲜核试验场地下爆炸的地震信号的检测能力平均可以达到2.0级，而事件的检测能力则可以达到3.5级（Kværna et al，2007 ）。因此可以判断，KSRS向IDC实时传输数据，特别是USRK台阵建成后，IMS地震台网对朝鲜地下核试验场的地震信号检测能力将会有很大的提高。Kværna等在进行阈值监测能力评估时，只在低于4 Hz的较低频带进行计算，而决定地震台站信号检测能力的台站本底噪声在高频段的强度要远低于低频段，同时地下爆炸相对于同震级的天然地震，震源源频谱的拐角频率较高，导致高频地震波激发更强。因此在区域震震中距范围内，地震台站在中高频段的信号检测能力可能要强于低频段（Bahavar，North，2002）。阈值监测方法作为一种地震台网实时监测能力的评估方法，被广泛地应用于地震台站（网）的监测能力实时计算中（王燕等，2011）。然而，利用阈值监测方法虽然可以得到监测系统在运行过程中的瞬时监测能力，但其计算结果受当前地震台站的非平稳噪声和非目标区域地震信号的干扰较大。因此，阈值监测方法只能确定监测目标区域是否发生震级大于当前阈值的地震事件，其较短时间内的评估结果并不能反映地震台站的真实信号检测能力以及地震台网的事件监测能力。此外，基于地震观测报告和区域地震目录，应用“震级-序号”，“最大曲率”及“完整性震级范围”等方法分析区域台网的目录最小完整性也是评估地震台网监测能力的常用方法（李志海等，2011；冯建刚等，2012）。相比较而言，概率模型法是在对每个地震台站的本底噪声进行统计、建立噪声水平分布概率模型的基础上，对地震台站的信号检测能力进行评估，因此其评估结果更能反映单个地震台站的信号检测能力。
为研究当前全面禁止核试验条约组织对朝鲜核试验场的监测能力，本文将以地下爆炸的源强度比例修正模型为基础，应用概率模型法评估属于国际监测系统的USRK和KSRS台阵对朝鲜地下核试验场地下爆炸地震信号的检测能力。
1. 方法原理
地震台站的信号检测能力主要取决于台站处的背景噪声水平、传播路径的衰减及地震设备的性能指标（何少林，2003）。
根据Gerlach等（1966）的研究，台站噪声幅值的对数值符合正态分布，即
$\lg {A_{\rm{n}}} {\text{～}} N(\mu{\text{，}}{\gamma ^2}){\text{，}}$
(1) 式中，An为噪声的速度幅值，μ为噪声幅值的期望，γ为该分布的标准差。对于幅值为As的信号，其在地震台站上的信噪比大于K的概率Pd可以表示为
${P_{\mathop{\rm d}\nolimits} } \text{＝} P(K{A_{\rm n}} \text{＜} {A_{\rm s}}) \text{＝} P\left({\rm lg}\frac{{K{A_{\rm n}}}}{{{A_{\rm s}}}} \text{＜} 0\right){\text{，}}$
(2) 式中P表示概率运算。因此，在给定信号信噪比的条件下，地震台站能检测到信号的概率就是随机变量X=lg（KAn/As）＜0的概率。如果将信号幅值As看作一个确定值，由式（1）可知，X的分布形式为
$X {\text{～}} N({\rm lg}K \text{＋} \mu \text{－} {\rm lg}{A_{\rm s}}{\text{，}}{\gamma ^2}){\text{，}}$
(3) 求X＜0的概率并将其转换为标准正态分布形式，即可得到对不同信号幅值的检测概率
$P(X \text{＜} 0) \text{＝} \varPhi \left(\frac{{\lg {A_{\rm s}} \text{－} \mu \text{－} \lg K}}{\gamma }{\rm{ }}\right){\text{，}}$
(4) 式中Φ为标准正态分布概率。通过查阅标准正态分布表，可以得到不同信号幅值所对应的检测概率，也可反推一定检测概率对应的阈值。对于特定台站对特定场地的监测，可以利用台站记录到的历史地震事件的信号幅值得到信号的幅值-震级关系，从而利用信号幅值反算震级。朝鲜历次核试验的震级集中在mb4.0—5.0，根据地下核试验的震源频谱模型，在拐角频率fc以下，地震台站记录到的不同核试验的信号幅值对数差与试验之间的震级差相同，因此对于拐角频率以下的信号可以利用已知震级核试验的实测信号幅值，根据比例关系折算得到不同震级地下核试验的信号幅值；但当信号的频率大于拐角频率时，震级差与信号的幅值对数差之间不再满足比例关系，则不能直接利用已知震级核试验的高频实测信号幅值折算更小震级核试验的高频信号幅值。为此，本文运用地下爆炸的源频谱模型，利用已知较大震级核试验实测的高频信号得到较小震级地下核试验的高频信号幅值，从而实现地震台站高频信号检测能力的评估。
对于球对称地下爆炸源，假设作用在弹性边界上的应力时间函数为阶跃函数，则距离弹性边界r处的远场径向位移的频谱可以近似表示为
$u_r {\text{≈}} \frac{{{{\varPhi }_\infty }}}{{r\left[ \text{－} {{\left(\displaystyle\frac{\omega}{\omega _{\rm{e}}}\right)}^2} \text{＋} {{2\eta }{\displaystyle\frac{\omega} {\omega _{\rm{e}}}{\rm i}} \text{＋} 1} \right]}}{\text{，}}$
(5) 式中：ωe为源本征频率，ωe＝2β/re＝2π fc；re为弹性半径；η＝β/α为源区附近介质的阻尼系数；β和α分别为介质的S波和P波速度；Φ∞为爆炸源的折合位移势的稳态值。根据Denny和Johnson （1991）的结果，在固定埋深的情况下，稳态折合位移势正比于爆炸当量W，即
${{\varPhi }_\infty } \propto W{\text{，}}$
(6) 同时式（5）中的源本征频率与爆炸当量之间存在比例关系
${\omega _{\rm{e}}} \propto {W^{ - {{1/ 3}}}}{\text{.}}$
(7) 当距离r大于爆炸的弹性边界且满足远场条件时，可以将式（5）所表示的频谱称为爆炸的源时间函数S（ω）。假如从r处到台站的地震波传播可以用传播函数T（ω）表示，则台站记录到的位移频谱可以表示为
$A\left( \omega \right) \text{＝} S\left( \omega \right)T\left( \omega \right){\text{.}}$
(8) 当两次爆炸的爆心位置的差异非常小时，可以认为从r处到台站的地震波传播函数T（ω）相同，并且假设两次爆炸源区附近的介质完全相同，这时两次爆炸在同一地震台站上记录到的地震信号之间的谱比值为
$\frac{{{A_1}\left( \omega \right)}}{{{A_2}\left( \omega \right)}} \text{＝} \frac{{{S\!_1}\left( \omega \right)}}{{{S\!_2}\left( \omega \right)}} \text{＝} \frac{{{W\!_1}}}{{{W\!_2}}} \frac{{{{\left({{{{\displaystyle\frac{\omega }{{{\omega _{{\rm e}2}}}}}}}}\right)}^2} \text{－} 2\eta {{{{\displaystyle\frac{\omega }{{{\omega _{{\rm e}2}}}}}}}}{\mathop{\rm i}\nolimits} \text{－} 1}}{{{{\left( {{{{{\displaystyle\frac{\omega }{{{\omega _{{\rm e}1}}}}}}}}} \right)}^2} \text{－} 2\eta {{{{\displaystyle\frac{\omega }{{{\omega _{{\rm e}1}}}}}}}}{\mathop{\rm i}\nolimits} \text{－} 1}}{\text{，}}$
(9) 式中，W1和W2为两次爆炸的当量，ωe1和ωe2分别为两次爆炸的源本征频率。因此，当已知某次当量为W1的地下爆炸的源本征频率为ωe1，其相应的mb震级为mb 1 ，其在某台站上的地震信号频谱为A1（ω）。根据式（9）可以得到当量为W2的地下爆炸在该台站上的地震信号频谱A2（ω）。同时在固定埋深的情况下，当量为W2时地下爆炸的震级mb 2可由
${{m}}{_{\rm{b}}^2} \text{＝} {{m}}{_{\rm{b}}^1} \text{－} \lg \frac{{{W_1}}}{{{W_2}}}.$
(10) 计算得到。通过以上比例折算，我们可以根据当量为W1，震级为mb 1的爆炸在频率ω1处的信号幅值，得到震级为mb 2的爆炸在相应频率处的信号幅值As。
2. 台站检测能力评估
2.1 台阵基本情况
本文所研究的USRK和KSRS台阵分别位于俄罗斯的乌苏里斯克和韩国原州，是目前IMS系统中距离朝鲜地下核试验场最近的两个台阵。KSRS台阵距朝鲜丰溪里核试验场约430 km，由28个子台构成，由于数据连续性问题，本文实际上仅使用中心区域的16个台站参与聚束。USRK台阵相距朝鲜地下核试验场约400 km，由10个子台构成孔径约为3.6 km的3层同心圆小孔径台阵。
2.2 台站本底噪声水平建模
概率模型法是通过统计建立噪声幅值样本的概率模型来确定对不同震级事件的检测概率。噪声统计过程中要避免地震事件信号及信号尾波的干扰。本文选取IDC公报中2014年至2016年全球范围内mb≥4.5地震共计约3 500次，截取每个事件中P波震相到时前15 s至前5 s的波形作为噪声波形。事件选取同时考虑不同季节、一天内不同时段可能存在的噪声强弱差异。为验证所选噪声信号样本在时间上是均匀分布的，分别统计一年和24小时的噪声样本数，其分布如图1所示。从统计结果可以看出，噪声样本在一天内不同时刻基本平均分布，按月份分布未表现出明显的季节差异性，该噪声样本可以反映地震台站的平均本底噪声水平，可用于后续计算。
![]()
图 1 一年内（a）和一天内（b）的噪声样本分布Figure 1. The distribution of noise samples in a year （a） and a day （b）对噪声样本数据按以下流程进行处理：首先，对台阵各子台记录到的波形数据进行质量控制，去除毛刺、记录为零或靠摆的数据段，再对质量合格的波形进行聚束和去仪器响应得到地面的振动速度波形（单位：μm/s）。由于存在连续波形缺失、数据质量不合格、数据间断等问题，最终两个台阵分别形成了大约3 200个有效噪声样本数据；其次，测量不同频带内的噪声幅值。为避免滤波产生的边界效应，以事件信号到时为中心前后各取100 s波形数据进行滤波。滤波频带分别为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，共计6个频带。为避免走时计算不准确导致事件信号本底进入噪声幅值测量时间窗，截取滤波后事件信号理论到时前15 s至前5 s共计10 s的波形作为噪声波形，测量每段噪声波形的最大峰-峰值作为其噪声幅值；最后，统计噪声峰-峰幅值对数的均值和方差。
得到噪声峰-峰幅值对数的均值和方差后，即可利用正态分布模型对不同频带内的噪声水平概率分布进行拟合。图2给出了不同频带内USRK台阵和KSRS台阵的本底噪声幅值对数的累积概率分布。可以看出，USRK台阵各频带的噪声幅值对数的均值随着频率的升高逐渐降低，而KSRS台阵噪声幅值对数的均值从0.75—1.5 Hz至2.0—4.0 Hz逐渐降低，但随后又升高。这可能是由于KSRS台阵的仪器响应校正曲线在高频段不准而引起。两个台阵噪声幅值对数的标准差在各频带无明显变化。
![]()
图 2 USRK台阵（a）和KSRS台阵（b）在各频带的累积概率分布曲线Figure 2. Cumulative probability distribution curves of the KSRS array （a） and the KSRS array （b） in the multi-band2.3 多频带信号幅值折算
根据式（9）和（10），只要利用一次已知当量地下核试验在地震台站上的实际观测信号幅值，即可预测出该爆心附近其它当量地下核试验在该台站上的信号幅值。朝鲜核试验之后，国内外有大量研究人员均尝试对其试验当量进行估算（Rougieret al，2011 ；Zhao et al，2012 ）。其中Rougier等（2011）以地表未出现明显破坏作用为边界条件，采用流体动力学数值模拟的方法对核试验的比例埋深进行约束，给出了2009年5月25日朝鲜第二次地下核试验的当量为5.7 kt，埋深为375 m。根据Denny和Johnson （1991），当量为Y，埋深为h的地下核试验的拐角频率可以表示为
${f_{\rm c}} \text{＝} \frac{{0.64{{\left( {\rho {\rm g}h} \right)}^{0.5522}}{{10}^{0.0025P_{\rm G}}}}}{{\pi \; {\rho ^{0.7245}}{\beta ^{0.0642}}{Y^{0.3333}}}}{\text{，}}$
(11) 式中，ρ为源区介质密度，h为试验装置的埋深，PG为介质的干燥孔隙率。根据Rougier等（2011），取密度为2 680 kg/m3，P波速度为5 495 m/s，S波速度为3 269 m/s；干燥花岗岩介质的孔隙率大约为0.5% （Koper et al，2008 ）。根据2009年5月25日朝鲜地下核试验的当量和埋深，可以估算其拐角频率为2.1 Hz。为了提高拐角频率估算的准确性，本文在1.1—3.1 Hz范围内对最优拐角频率进行搜索，即通过
$\min \left( {\sum\limits_{k = 1}^K {d_k^2} } \right)$
(12) 来搜索一个最优拐角频率fcopt。式中dk为地震台站上记录到的Pn波的位移谱与理论位移谱之间的残差，表示为
${d_k} \text{＝} A_{{\rm{o}}}^{{\rm{Pn}}}\left( {{f_k}} \right) \text{－} A_{{\rm{t}}}^{{\rm{Pn}}}\left( {{f_k}，{f_{\rm c}}} \right){\text{，}}$
(13) 式中Ao Pn（fk）和At Pn（fk）分别为Pn波在第k个频点fk上的观测位移谱和理论位移谱。根据Sereno等（1988），台站上记录到的Pn波位移理论振幅谱可以表示为
$A_{{\rm{t}}}^{{\rm{Pn}}}\left( {f{\text{，}}\!\!\!{f_{\rm c}}} \right) \text{＝} S\!\left({f{\text{，}}{f_{\rm c}}} \right)G\left( {\varDelta{\text{，}}f} \right)\exp \left( { \text{－} \frac{{{\rm{\pi }}f\tau }}{{Q\left( f \right)}}} \right){\text{，}}$
(14) 式中，S（ f，fc）为源频谱，G（Δ，f ）为Pn波从爆心到地震台站的几何扩散因子，Δ为爆心到台站的距离，Q（ f ）为传播路径的平均非弹性衰减，τ为Pn波的传播走时。根据Yang等（2007）对Pn波几何扩散因子的研究结果，G（Δ，f ）可以表示为
$G\left( {{\rm{\varDelta }}{\text{，}}f} \right) \text{＝} \frac{{{{10}^{{n_3}\left( f \right)}}}}{{{{\rm{\varDelta }}_0}}}{\left( {\frac{{{{\rm{\varDelta }}_0}}}{{\rm{\varDelta }}}} \right)^{{n_1}\left( f \right)\lg (\frac{{{{\rm{\varDelta }}_0}}}{{\rm{\varDelta }}}) \text{＋} {n_2}\left( f \right)}}{\text{，}}$
(15) ${n_i}\left( f \right) \text{＝} {n_{i1}}{\left[ {\lg^2 \left( {\frac{f}{{{f_0}}}} \right)} \right]} \text{＋} {n_{i2}}\lg \left( {\frac{f}{{{f_0}}}} \right) \text{＋} {n_{i3}}{\text{，}}\quad i \text{＝} 1{\text{，}}2{\text{，}}3{\text{，}}$
(16) 式中Δ0和f0分别为参考距离和参考频率。利用朝鲜3次地下核试验激发的Pn波在297个宽频带地震台站上的记录信号，Zhao等（2015）对适用于朝鲜半岛及我国东北地区的Pn波几何扩散因子系数进行了搜索反演，同时给出了朝鲜半岛和我国东北地区的Pn波视衰减Q值为237，视衰减Q值的频率指数为0.361。根据Denny和Johnson （1991），当量为Y，埋深为h的地下核爆炸的稳态折合位移势可以表示为
${\varPhi _\infty } \text{＝} \frac{{3.4 \text{×} {{10}^9}Y}}{{{\beta ^{1.1544}}{{\left( {\rho {\rm g}h} \right)}^{0.4385}}{{10}^{0.0344P_{\rm G}}}}}.$
(17) 利用式（17）和式（5），可以计算出1 km距离处的位移谱，将其作为式（14）中的S（f，fc），再利用几何扩散因子和视非弹性衰减系数即可对最优源拐角频率进行搜索。由于本文研究的两个地震台站只有USRK台站装有宽频带地震计，因此我们只使用USRK台站宽频带仪器记录到的2009年5月25日朝鲜地下核试验的Pn波位移谱对最优拐角频率进行搜索。图3给出了最优频率的搜索结果，可以看出，当源频谱拐角频率为2.7 Hz时，位移谱拟合残差取得最小值。图4给出了取最优拐角频率时理论Pn波位移谱对观测Pn波位移谱的拟合结果，可以看出，理论Pn波位移谱与实际观测谱基本吻合，表明2009年5月25日朝鲜地下核试验的源频谱拐角频率应该在2.7 Hz附近。
![]()
图 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本文以2009年5月25日的核试验作为参考爆炸事件，根据其当量、震级以及通过对实际观测数据拟合得到的源频谱拐角频率，利用比例关系得到了其它更小当量试验的震级及其与参考事件之间的源频谱比值。图5给出了部分小当量爆炸对参考事件的稳态强度进行归一化后的源频谱及相应的源频谱比值。
![]()
图 5 部分小当量爆炸的源频谱（a）以及相对于参考事件的谱比值（b）Figure 5. Source spectrum of explosion events （a） and spectral ratio to the reference event （b） with small yield在时间域内折算小震级核试验的信号幅值时，首选利用相应的Pn波慢度对USRK和KSRS台阵各台站记录到的2009年5月25日朝鲜核试验的信号进行聚束并去除仪器响应，然后再利用通带分别为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的3阶巴特沃斯带通滤波器进行滤波。图6给出了经过上述处理后的观测波形。
![]()
图 6 USRK （左）和KSRS （右）台阵的多频带波形数据Figure 6. Waveforms in the multi-band for the USRK (left) and KSRS (right) arrays从Pn波到时前0.5 s起，在宽为5.5 s的时间窗内测量信号的峰−峰幅值。表1给出了2009年5月25日朝鲜核试验在上述两个台阵上的实测Pn波幅值。不同震级地下核试验在上述两个地震台站上的Pn波预测幅值A2可以由
${A^2} \text{＝} {A^1} {\mathop{\rm mean}\nolimits} \left[ {\frac{{{S_2}\left( {{f_1}:{f_2}} \right)}}{{{S_1}\left( {{f_1}:{f_2}} \right)}}} \right] \text{＝} {A^1} {R_{21}}\left( {{f_1}:{f_2}} \right)$
(18) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
表 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 magnitudesmb 台阵 信号幅值/（μ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 2.4 台站多频带信号检测概率计算
利用USRK和KSRS台阵的本底噪声幅值对数均值和方差，将不同震级对应的各频带幅值带入式（4）并取信噪比K为3.0，可以计算得到这两个台阵对朝鲜核试验场不同震级地下核试验的3倍信噪比的信号检测概率，结果如图7所示。
![]()
图 7 USRK台阵（a）和KSRK台阵（b）各频带检测概率Figure 7. Detection probability of the USRK array （a） and the KSRS array （b） in the multi-band从图7中分别读取0.9，0.5，0.3的检测概率所对应的震级，所得到的值即为对应的检测阈值，结果如图8所示。可以看出，UKSK和KSRS台阵的检测阈值随频率的升高逐渐下降，0.9的检测概率条件下，在0.75—1.5 Hz频带内，两台阵的检测阈值分别约为mb3.5和mb3.4；在6.0—9.0 Hz频带内，两台阵的检测阈值分别为mb1.5和mb2.2。
![]()
图 8 不同震级事件USRK（a）和KSRS（b）台阵的检测阈值Figure 8. Detection threshold with different magnitudes for USRK array (a) and KSRS array (b)3. 讨论与结论
北京时间2010年5月12日8时8分45.067秒在朝鲜核试验场附近发生了一次震级约为mb（1.44±0.13）的事件（Zhang，Wen，2015）。根据本文的评估结果，USRK台阵在6.0—9.0 Hz频带可能记录到了该次小事件的信号。截取该时刻的波形，进行聚束处理，并按上文频带进行滤波，结果如图9所示。可以看出，在4.0—8.0 Hz和6.0—9.0 Hz频段内均可以清晰地看到该次事件的Pn震相，其信噪比分别约为3.4和2.5。f-k分析的结果表明，该信号的方位角和水平慢度与来自朝鲜核试验场附近的地震信号的方位角和水平慢度相一致，表明图9中4.0—8.0 Hz和6.0—9.0 Hz的微弱地震信号来自于2010年5月12日朝鲜核试验场附近的小震级地震。由此可以验证USRK台阵对朝鲜核试验场附近的地震事件的信号检测能力可以达到mb（Lg）1.44。考虑到2010年5月12日朝鲜核试验场地震事件的震级测量不确定范围及mb（Lg）震级与mb震级之间的等效性，这一观测结果与本文得到的USRK台阵对朝鲜核试验场具有mb1.5的信号检测能力的结果相吻合。同时，KSRS台阵在上述频带并未检测到该次小地震的信号，说明KSRS台阵的检测阈值要大于mb1.44。
![]()
图 9 2010年5月12日朝鲜核试验场小震事件的波形图中时间轴零点处为Pn波信号到时位置Figure 9. Waveforms of the small earthquake on May 12，2010 from the North Korea nuclear test siteThe zero position of the time axis is the arrival time of Pn wave此外，从评估结果可以看出，在0.75—1.5 Hz和2.0—4.0 Hz频带内，KSRS的信号检测阈值比USRK分别弱0.1和0.2个震级单位，这与KSRS台阵距离朝鲜核试验场更远有关。但在6.0—9.0 Hz频带内，KSRS的检测能力要比USRK弱0.7个震级单位。从表1和图2b可以看出，在该频带内，KSRS台阵的信号幅值要比USRK强13倍，同时本底噪声水平也高50倍。通过分析这两个台阵的仪器响应特性，我们认为可能是由于IDC提供的KSRS台阵仪器响应曲线在6.0—9.0 Hz频带内存在过校正使得两个台阵之间的信号幅值和噪声水平有较大差异。不过由于过校正同时作用于信号和噪声，因此并不会影响信号检测能力的评估结果。KSRS台阵的卫星影像图表明，该台阵距离居民区和河流的距离较近，这可能是导致其高频本底噪声水平偏高、高频信号检测能力下降的主要原因。
本文通过台站本底噪声建模及多频带信号幅值折算，计算得到USRK台阵和KSRS台阵对朝鲜核试验场不同震级地下核试验的Pn波信号检测概率，进而得到不同信号检测概率对应的震级阈值。在0.9的检测概率条件下，USRK台阵和KSRS台阵在0.75—1.5 Hz，2.0—4.0 Hz，6.0—9.0 Hz频带对朝鲜核试验场地下核试验的检测阈值分别为mb3.5，mb2.3，mb1.5和mb3.4，mb2.5，mb2.2。结果表明：震中距为400 km左右的两台阵均随着频率的升高，检测阈值逐渐降低，即信号检测能力增强；在高频段，USRK的监测能力强于KSRS。
-
![]()
图 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.
计量
- 文章访问数: 1514
- HTML全文浏览量: 689
- PDF下载量: 40
