从2006年10月9日朝鲜第一次地下核试验以来，分别于2009年5月25日，2013年2月13日和2016年1月6日相继进行了三次规模较大的核爆试验(以下简称NKT1，2，3和4).四次朝鲜核爆的试验地点非常接近，在中国东北的地震记录中能够清楚地识别出区域震相中的Pn，Pg，Lg和Rayleigh波(例如，Zhao et al.，2008，2012，2014).由于NKT2和NKT3的震级较大，地震记录的信噪比相对较高，甚至在跨过日本海的日本F-NET台网的台站也记录到清楚的P波震相.朝鲜核试验所产生的地震数据，因为震源简单，记录质量高，能够为估算爆炸当量、分析朝鲜半岛核爆事件与地震事件之间的震源差异提供重要的信息，还可以用于研究核爆记录覆盖区域的地壳和上地幔的地震速度和衰减结构.因此，在全世界范围内激发了许多研究者的极大兴趣(Kim and Richards，2007；Kværna et al.，2007；Salzberg and Marshall，2007；Bonner et al.，2008；Koper et al.，2008；Hong et al.，2008；Patton and Taylor，2008；Salzberg，2008；Tibuleac et al.，2008；Zhao et al.，2008；Hong and Rhie，2009；Kim et al.，2009；Chi et al.，2010；Murphy et al.，2010，2011；Ni et al.，2010；Schlittenhardt et al.，2010；Shin et al.，2010；Rougier等，2011；Zhang and Wen，2013).

朝鲜核爆发生以后，很多研究者给出了地震学当量的估计.但是，因为使用资料和方法不同，当量估计结果差异较大.例如，对于2006年10月9日朝鲜地下核试验，Salzberg和Marshall(2007)利用同一台站不同事件地震记录的能量包络比值法，参照中国东北地区用于地震测深的已知当量的化爆事件，获得NKT1的当量估计为0.45 kt.基于最小埋藏深度假定，Zhao等(2008)利用中国东北及邻近地区的地震Lg波资料得到的当量估计是0.48 kt.但是，Chun等(2011)使用中朝边界附近的流动台网记录的区域Lg波资料得到的结果是1.22 kt.Murphy等(2013)利用远震P波振幅谱得到的核爆当量与埋深分别为0.9 kt和200 m.

早期核试验的核爆当量较高，通常是计算远震P波震级，利用震级与当量之间的经验关系获得爆炸当量的估计(例如，Ringdal et al.，1992;Douglas and Marshall，1996;Murphy and Barker，2001).也可以通过远震P波波形反演求解核爆等效震源参数，利用相对静力学强度给出核试验的当量估计(赵连锋等，2005).随着核武器小型化，地下核试验的当量逐渐减小，区域地震记录成为研究核爆当量的首选资料.例如，Vergino和Mensing(1990)发展了使用区域震相Pn波计算体波震级的方法；Nuttli(1973，1986a，1986b)以Lg波列中第三峰值振幅为基础发展了一种区域地震级系统，并建立了Lg波震级与爆炸当量之间的关系.Patton(1988)将Nuttli的方法推广应用到Nevada测试点(NTS)的更大的数据体.Hansen等(1990)和Ringdal等(1992)发现基于Lg波均方根(rms)振幅的震级估计更加稳定，即使台站很少，也能给出令人满意的当量估计.

最近，一些研究者尝试使用区域地震波形的高频波列，建立了计算震级和估计爆炸当量的尾波法(Mayeda and Walter，1996；Mayeda et al.，2003；Mayeda，1993；Murphy et al.，2009).尾波是地方震记录图中体波或面波之后的某些震荡，通常被认为是由于分布于地壳和上地幔的横向不均匀性引起的散射波，尾波的激发依赖于台站的地质情况(Aki，1969；Aki and Chouet，1975).在地方性地震范围内，体波震相Pn，Pg，Sn和短周期面波Lg较容易识别，但波形较为复杂且容易受到源和路径不均匀性的影响，但其尾波对辐射花样、传播路径等不敏感，波形特征稳定.因此，地方性地震范围内单个地震台站体波震相的尾波有可能用于测量地下核试验的当量和埋深.Pasyanos等(2012)将这种使用单一震相尾波方法扩展到多个震相，以等效震源函数、传播效应、台基响应和传递函数的物理模型为基础，建立了一种通过单一地震台站拟合区域波形全波列的包络，同时获得地下核试验当量和埋藏深度的方法.

我们收集中国东北地区距离朝鲜核试验场较近的11个地震台站的地方及区域地震波形资料，参考Pasyanos等(2012)的全波列尾波方法，调查四次朝鲜地下核试验的爆炸当量和埋藏深度.

选取中国东北地区距离朝鲜核试验场(表 1)震中距范围在200 km到500 km之间的11个地震台站以及2002年到2013年间该区域发生的27个天然地震记录(M_L≥2.5)来进行数据分析(表 2).

表 1(Table 1)

表 1 朝鲜地下核试验事件目录 Table 1 Event parameters of the North Korean underground nuclear tests 朝鲜 核试验 时间 (yyyy-mm-dd) 发震时刻 (hh:mm:ss) UTC 纬度(°N) 经度(°E) NKT1 2006-10-09 01:35:28.00001 41.28742 129.10832 NKT2 2009-05-25 00:54:43.1142 41.2940 129.0775 NKT3 2013-02-12 02:57:51.2741 41.2918 129.0733 NKT4 2016-01-06 01:30:00.9706 41.3003 129.0678 注：数据源1.USGS(Wen and Long,2010;Zhao et al.,2014).2.卫星图像(Wen and Long,2010;Zhao et al.,2014).

表 1 朝鲜地下核试验事件目录 Table 1 Event parameters of the North Korean underground nuclear tests

表 2(Table 2)

表 2 区域地震目录 Table 2 Earthquake parameters used in this study 时间 (yyyy-mm-dd) 发震时刻 (hh:mm:ss)UTC 纬度(°N) 经度(°E) 2002-04-16 22:52:38 40.70 128.67 2002-05-05 10:02:24 40.06 127.24 2002-11-25 11:19:38 41.80 128.17 2002-12-24 15:57:56 41.95 128.15 2003-03-03 08:40:00 41.49 126.91 2003-03-19 10:15:12 42.02 127.43 2003-07-24 09:54:22 41.87 128.19 2003-10-25 11:03:08 41.92 128.10 2004-09-08 14:24:21 42.10 128.02 2004-12-16 18:59:12 41.80 127.98 2005-04-15 06:34:15 42.15 127.71 2005-05-30 03:40:31 41.32 128.25 2006-02-27 15:32:47 40.09 127.32 2006-11-11 14:11:08 41.74 128.03 2007-12-28 07:14:40 41.59 126.83 2008-05-30 00:17:11 41.34 123.36 2008-06-12 06:26:00 42.93 124.74 2008-09-03 13:43:44 42.75 124.16 2009-02-24 18:05:42 41.84 123.74 2009-08-05 12:08:12 42.35 127.22 2009-12-27 17:36:25 41.86 123.56 2010-10-09 05:45:14 42.37 128.42 2010-10-09 06:07:09 42.35 128.39 2011-02-18 15:25:58 41.73 126.89 2011-06-08 22:03:30 44.82 129.88 2013-08-11 07:46:38 42.90 124.40 2013-08-17 06:56:36 41.90 123.90

表 2 区域地震目录 Table 2 Earthquake parameters used in this study

基于尾波单一散射模型(Aki，1969；Aki and Chouet，1975)，Mayeda等(2003)给出在中心频率f和震中距r的单个区域震相尾波包络(图 2)计算公式：

图 1

Fig. 1

图 1 研究区域背景图.朝鲜地下核试验点(五角星)，天然地震记录(黑色十字形)，区域台站(黑色三角形) Fig. 1 Map showing the background of study area and North Korean underground nuclear tests(star),earthquakes(black cross),and stations(black triangles)

图 2

Fig. 2

图 2 单个震相理论包络形状 Fig. 2 The shape of single phase theory envelope

(1)

其中S(f)为震源谱，G(r)为几何扩展项，B(f，r)为衰减项，P(f)为台基响应，T(f)为传递函数，E(t，f，r)为描述尾波形状的函数，

(2)

H为单位阶跃函数；t_a为到时；γ(f，r)和b(f，r)为尾波形状参数.

地方性或区域地震范围地震全记录的包络(图 3a)为该区域噪声水平A_noise与所有区域震相包络之和(Pasyanos et al.，2012)：

图 3

Fig. 3

图 3 区域事件理想包络形状(a)和实际波形包络形状(b).(b)为MDJ台站记录的2002年4月16日天然地震的不同频带的包络，横轴为时间(取发震时刻后250 s的记录)，纵轴为速度振幅的对数. Fig. 3 (a)Expected envelope shapes for regional events;(b)Shapes of real waveforms envelopes in different frequency bands of the earthquake event on 16 April 2002 recorded at MDJ station.X-axes are times(250 s after the origin time).Y-axes are the log amplitude of velocity.

(3)

应用于实际观测资料，首先将两个水平分量记录经过反褶积去掉仪器响应，获得地面速度记录，然后，利用八阶零相(四个极点，双通)巴特沃兹滤波器得到窄频带波形，应用窄带包络公式(f为中心频率，v(t)为一个改正后的带通滤波地震图地面速度，h(t)为v(t)对应的希尔伯特变换)：

(4)

求出两水平分量的窄带包络，再求取两个水平包络的平均值，经过平滑处理，最后获得该频带的尾波包络(图 3b).这个过程也能用单分量(如垂直分量)来完成，但比起单分量，两个水平分量能提供一个更平稳光滑的包络.对于天然地震事件和地下核试验事件，其尾波包络(1)式中的各参数项只有震源项不同，其他参数项都可以由天然地震记录资料获得.应用于地下核试验事件中，通过天然地震记录构建地下核爆的理论包络，与实际结果比较，可以给出地下核试验的爆炸当量和埋藏深度(Pasyanos et al.，2012).

我们用事件前噪声水平来经验地估计A_noise，取Pn到达之前的噪声窗，进行统计处理，得到均值.对于壳幔模型，地壳用单层均匀介质层表示，壳幔边界水平，忽略坡度变化，地幔水平方向分块，每个方块内速度均匀，射线在地幔顶部走直线，忽略速度横向变化引起的实际射线弯曲，忽略地球曲率和地幔顶部垂向速度梯度的影响(裴顺平，2002).由于所用数据中天然地震和四次核试验事件的震源非常浅，震中距离较近，我们将Pg和Lg波的传播路径近似为震源到台站的直线(图 4).Pn和Sn波的走时方程可描述为：

图 4

Fig. 4

图 4 区域震相射线路径图及壳幔模型 Fig. 4 Ray paths of regional phases and the crust-mantle model

(5)

其中t为走时，T₁为到时，T₀为发震时刻，d_s、d_m和d_r分别为震源路径、地幔路径和接收路径的长度，s_c为地壳慢度，s为地幔慢度，dl为对路径的积分.假设单层水平地壳，其内及地幔顶部的速度不变，根据Snell定律可将(5)式转化为：

(6)

其中H_s为震源处地壳厚度，h_s为震源深度，H_r为台站处地壳厚度，h_r为台站高程，r为震中距.

中国东北地区Pn波平均波速度约为8.0~8.05 km·s^-1，单层水平地壳P波平均速度约为6.3 km·s^-1(贾军等，2011；裴顺平等，2004a；张国民等，2012)，Sn波平均速度约为4.55 km·s^-1，地壳S波速度约为3.6 km·s^-1(裴顺平等，2004b)，Pg波平均速度约为6.0 km·s^-1，Lg波群速度为3.5 km·s^-1(赵连锋，2008).

天然地震震源采用点源模型，用Brune(1970)ω²震源模型表示，其中S₀(f)为零频震源谱，f_c为源拐角频率，M₀为源的标量地震矩，ρ_s为源介质密度，ρ_r为接收处介质密度，c_s为源介质速度，c_r为接收处介质速度，R_θφ为平均辐射花样，在研究区域内相关参数值R_θφS=0.63(Sn，Lg)，R_θφP=0.43(Pn，Pg).根据Walter和Taylor(2001)的MDAC2模型得到拐角频率、地震矩及视应力间的关系：

(7)

(8)

(9)

其中σ_a为视应力，M′₀为参考地震矩，σ′_a为在M′₀处的应力，ψ为比例系数(对于恒定的视应力ψ=0)，ζ为P，S波拐角频率比例因子，P波和S波的拐角频率相近时ζ趋于1.α_s和β_s分别为源P，S波速度.

采用Mueller-Murphy模型(Mueller and Murphy，1971)来描述爆炸震源的震源谱，爆炸源P波部分震源谱的表达式为：

(10)

其中f_c为拐角频率；ρ_s为源介质密度；，v_s(P)和v_s(S)分别为近源P，S波速；R_e为弹性半径，R_e=kY^1/3(ρgh)^-0.417，k为比例系数，Y为当量，ρ为介质密度，g为重力常数，h为埋藏深度；f₀=，f₁与R_e成反比；P_p=1.5ρgh为峰值压力；在不可压缩下，为稳态压力；R_c=cY^0.29h^-0.11为在花岗岩中的空腔半径，c为与介质有关的参数.部分参数值ρ_s=2700 kg·m^-3，v_s(P)=6.0 km·s^-1，v_s(S)=3.5 km·s^-1，R_e=257.7 m，c=24.6619，P₀=2.76×10⁶ N·m²，f₁=4.35.爆炸源S波震源谱，S波拐角频率与P，S波速比有关，用相同的P，S波速比来设定低频S波水平，并用于与P波相同的函数形式模拟(Fisk，2006).

尾波形状参数b(f，r)和γ(f，r)代表尾波振幅的衰减，分别依赖于介质内在吸收引起的尾波几何扩展和衰减速率，b(f，r)控制尾波包络后半部分的斜率，γ(f，r)控制着紧跟直达波之后的尾波初期的衰减，γ(f，r)值越大表明尾波初期更陡峭或是有更大的直达峰值.根据Mayeda等(2003)给出的在地方性及区域地震范围内b(f，r)和γ(f，r)的函数关系：

(11)

(12)

其中γ₀，γ₁，γ₂和b₀，b₁，b₂为常量，可用网格搜索法确定.我们利用区域内27个天然地震的地震记录资料，在选定的频带内为每个地震事件-台站对用网格搜索找出γ的最佳参数值，γ的变化范围从0到1.6每次增加0.1.当γ固定之后用最小二乘法来反演求解b(Yoo et al.，2011).通过匹配E(t，f，r)与实际波形的衰减速率，由S波尾波和Lg波尾波校准得到我们所用区域台站的尾波衰减参数b(f，r)和γ(f，r)的数值(图 5).

图 5

Fig. 5

图 5 MDJ台站尾波衰减参数γ(f,r)和b(f,r)在频带2~4 Hz内所得结果 Fig. 5 Coda-decay parameters γ(f,r) and b(f,r) in the 2~4 Hz frequency band at MDJ station

不同震相的几何扩展项G(r)的函数形式为(Zhu et al.，1991；Walter and Taylor，2001)：对于Pn、Pg和Sn震相，G(r)=(1/r)^1.1；对于Lg震相，.衰减项B(f，r)=，品质因子Q=Q₀f^η，Q₀为f=1时Q值，η为频率关系指数.Xie等(2006)调查了欧亚大陆的Lg波Q值并给出横向的Q₀值图像，区域内Q₀均值约为450；赵连锋(2008)给出中国东北及邻近地区地震Lg波Q值层析成像的结果，平均Q₀值为429；Hong和Rhie(2009)通过区域地震台站2002年4月16日天然地震记录，给出Pn、Pg、Sn和Lg波衰减参数Q₀和η值.根据上述资料，MDJ台站所用到的衰减参数如表 3所示.

表 3(Table 3)

表 3 MDJ台站各震相的Q0值和η值 Table 3 Values of Q0 and η at MDJ station Q0 η Pn 230 0.48 Pg 500 0.18 Sn 180 0.49 Lg 429 0.15

表 3 MDJ台站各震相的Q0值和η值 Table 3 Values of Q0 and η at MDJ station

尾波包络(1)式中几何扩展、衰减和尾波形状函数项确定后，将S(f)归一化，两边取对数log₁₀，应用L1范数改正并通过网格搜索比较MDJ台站记录的2002年4月16日和2004年12月16日两次天然地震事件来得到MDJ台站台基响应和传递函数，即对理论包络做出改正，再应用于地下核试验事件中(图 6).

图 6

Fig. 6

图 6 MDJ台站理论包络改正 2002年(a)和2004年(b)地震事件MDJ台站的拟合结果(2~4 Hz频带内)，实际波形包络(黑色波形)，未改正的地震理论包络(蓝色虚线)，改正后的地震理论包络(绿色实线).横轴为时间(取发震时刻后250s的记录)，纵轴为速度振幅的对数. Fig. 6 Corrected theoretical envelopes at MDJ station (a)Seismic event of 2002(2~4 Hz);(b)Seismic event of 2004(2~4 Hz)with envelopes of real data(black waveforms),uncorrected theoretical seismic envelopes(blue dash lines),corrected theoretical envelopes(green lines).X-axes are times(250 s after the origin time).Y-axes are log amplitudes of velocity.

MDJ台站记录到朝鲜四次地下核试验的垂直分量速度记录如图 7所示.对于朝鲜的四次地下核试验事件，根据中国东北地区地震台网的区域地震台站全波列尾波的理论包络，采用爆炸源模型，介质类型为花岗岩，通过网格搜索比较不同当量、不同埋深事件的理论包络与实际结果，确定四次地下核试验事件的最佳拟合理论包络(图 8，9).从图 8中可以看出，尾波包络的振幅对于当量和埋深的变化是比较敏感的.当量1 kt、埋深100 m的理论包络(蓝色实线)比2006年核试验事件实际包络的当量和埋深要大，但比2009年核试验事件实际包络的当量和埋深要小得多.对于2013年和2016年两次核试验事件，当量15 kt、埋深500 m的理论包络(蓝色实线)大于实际包络的当量和埋深，而当量5 kt、埋深500 m的理论包络(黑色虚线)小于实际包络的当量和埋深，对比两次核试验事件实际波形包络可以看出两次核试验事件的实际当量和埋深差别较小.我们将当量和深度分别从0.01 kt到100 kt、10 m到1000 m的范围内进行搜索，找出最佳拟合范围(图 9)获得2006年核试验事件的当量为0.6 kt(0.4~0.8 kt)，埋藏深度为150 m(50~250 m)，2009年核试验事件的当量为3 kt(1.5~4.5 kt)，埋藏深度为350 m(250~450 m)，2013年核试验事件的当量为10 kt(8~12 kt)，埋藏深度为500 m(300~700 m)，而2016年核试验事件的当量为8 kt(6~10 kt)，埋藏深度为500 m(300~700 m).

图 7

Fig. 7

图 7 MDJ台站四次朝鲜地下核试验事件垂直分量速度记录 波形上标出了视群速度值. Fig. 7 Vertical velocities of four North Korean underground nuclear tests recorded at MDJ station Apparent group velocity values are marked on waveforms.

图 8

Fig. 8

图 8 MDJ台站四次朝鲜地下核试验事件拟合结果(2~4 Hz频带内) 实际波形包络(黑色波形)，最佳拟合理论包络(红色实线).2006年(a)和2009年(b)图中蓝色实线为当量1 kt、埋深100 m的理论包络，2013年(c)和2016年(d)图中蓝色实线和黑色虚线分别为当量15 kt、埋深500 m及当量5 kt、埋深500 m的理论包络.横轴为时间(取发震时刻后250 s的记录)，纵轴为速度振幅的对数. Fig. 8 Results of four North Korean underground nuclear tests in 2006(a),2009(b),2013(c)and 2016(d)at MDJ station(2~4 Hz) The envelopes of real data(black waveforms),best fitting theoretical envelopes(red lines)are plotted.Blue lines in 2006(a)and 2009(b)represent the theoretical envelopes of 1 kt with 100 m depth,while for 2013(c)and 2016(d)the blue lines and black dash lines represent the theory envelopes of 15 kt with 500 m and 5 kt with 500 m depths,respectively.X-axes show times(250 s after the origin time).Y-axes show the log amplitudes of velocity.

图 9

Fig. 9

图 9 2006年(a)、2009年(b)、2013年(c)和2016(d)年四次朝鲜地下核试验事件的MDJ台站网格搜索结果 黑色五角星为最佳拟合，色标范围表示不拟合度由最佳拟合处向四周延伸越来越大. Fig. 9 Results of grid search for 2006(a),2009(b),2013(c)and 2016(d)North Korean underground nuclear tests at MDJ station The best fit is indicated by the star.The color scale from the best fit extending all around means misfit getting greater.

我们选用除MDJ台站外其他10个区域地震台站采用爆炸源模型来研究2013年地下核试验事件，结果如图 10所示.从图 10中可以看出，当量10 kt、埋藏深度500 m时，10个区域台站尾波的理论包络与实际资料都拟合较好，表明整个区域地震台网内的台站结果具有高度的相似性，说明地方性及区域地震范围内单个地震台站体波震相的尾波有可能用于测量地下核试验的当量和埋深.

图 10

Fig. 10

图 10 2013年朝鲜地下核试验事件10个区域台站拟合结果(2~4 Hz频带内) 实际波形包络(黑色波形)，最佳拟合理论包络(红色实线).横轴为时间(取发震时刻后250 s的记录)，纵轴为速度振幅的对数. Fig. 10 Fitting results of North Korean underground nuclear test in 2013 at 10 regional stations(2~4 Hz) Envelopes of real data(black waveforms),best fitting theoretical envelopes(red lines)are displayed.X-axes show the time(250 s after the origin time).Y-axes show log amplitudes of velocity.

Chi等(2010)，Murphy等(2010)，Salzberg(2008)，Zhang和Wen(2013)，Pasyanos等(2012)及Zhao等(2008，2012，2014，2015)都曾给出朝鲜地下核试验核爆当量和深度的研究结果.表 4给出了本文与前人调查结果的比较，可以看出本文结果与表中几位国内外专家的研究一致性较高，说明地方性及区域地震范围内单个地震台站体波震相的尾波是可能同时约束地下核试验的爆炸当量和埋藏深度的.

表 4(Table 4)

表 4 四次朝鲜地下核试验事件的研究结果 Table 4 Research results of four North Korean underground nuclear tests 2006-05-25 2009-10-09 2013-02-12 2016-01-06 当量(kt) 深度(m) 当量(kt) 深度(m) 当量(kt) 深度(m) 当量(kt) 深度(m) Chi 0.4(0.2~0.8) - 3.2(1.6~6.4) - - - - - Murphy 0.9 200 4.6 550 - - - - Salzberg 0.85±0.1 - - - - - - - Zhang - - 7.0±1.9 610 12.2±3.8 430 11.13±4.0* - Pasyanos 0.5±0.3 100(20~300) 2(1~5) 200(70~600) - - - - Zhao 0.48 - 2.35 - 7.47 - 3.58** - 本文 0.6±0.2 150±100 3.0±1.5 350±100 10.0±2.0 500±200 8.0±2.0 500±200 注：*中国科学院大学精确测定朝鲜核爆( http://news.cri.cn/201617/14d96c7a-fe28-4307-fcee-d4b4e2ed8a5a.html).**研究快报：2016年1月6日朝鲜核试验定位、识别和当量估计的初步调查结果(http://www.igg.cas.cn/xwzx/kyjz/201601/t20160107_4515052.html).

表 4 四次朝鲜地下核试验事件的研究结果 Table 4 Research results of four North Korean underground nuclear tests

尾波对辐射花样、传播路径等不敏感，且波形特征稳定，本文通过构建一个理论的区域震相尾波振幅包络，来拟合地方性及区域地震范围内单一地震台站实际波形记录，调查朝鲜四次地下核试验的爆炸当量和埋藏深度.所有区域震相的时间域振幅包络由其震源谱函数、传播特征、台基响应、传递函数及尾波形状函数组成，我们通过研究区域内天然地震记录资料对各参数项进行实际分析，得到计算尾波理论包络所需各参数值，例如27个天然地震记录资料来获得MDJ台站所选用频带内的尾波形状参数b(f，r)和γ(f，r)值.我们对当量和埋藏深度进行网格搜索来拟合台站水平分量记录的尾波包络，得到2006年核试验事件的当量为0.6 kt±0.2 kt，埋藏深度为150 m±100 m，2009年核试验事件的当量为3 kt±1.5 kt，埋藏深度为350 m±100 m，2013年核试验事件的当量为10 kt±2 kt，埋藏深度范围为500 m±200 m，而2016年核试验事件的当量为8 kt±2 kt，埋藏深度范围为500 m±200 m.

由于核爆当量与埋藏深度依赖于爆炸源模型和传递函数，而传递函数如前所述是由该区域内天然地震记录资料中经验地推导出的，需要在未来能更加透彻地理解爆炸源模型和传递函数，并希望能给出更好、更合理的解释.另外，区域范围内四个区域震相Pn，Pg，Sn，Lg在实际波形中能较清楚地观测到，而在地方震距离内(Δ<100 km)的地震台站上只有Pg和Lg两个震相，所以我们下一步的工作是研究从近场到区域地震范围内，该方法在更多更密集的台站上的具体表现.

感谢美国地震学研究机构联合会(IRIS)提供的四次朝鲜地下核试验的数据资料，感谢赵连锋副研究员、郝金来副研究员和赵旭副研究员提供的宝贵建议和帮助.