Distinguishing underground nuclear test by full moment tensor inversion
-
摘要: 将震源矩张量中的各分量作为权重因子,利用基本Green函数的线性组合可以对地震波场进行描述,本文利用水平分层模型下的广义反射-透射系数方法的地震波场正演公式,反演了苏联东哈萨克斯坦地下核试验场的7次地下核爆炸、3次震中位置相近的天然地震的震源矩张量,反演结果表明地下核爆炸震源明显含有爆炸源成分,同时包含双力偶成分(DC)和线性偶极补偿源成分(CLVD),且CLVD成分比重大,CLVD的存在可用层裂机制来解释;天然地震震源矩张量的反演结果表明,DC源为主要成分,符合剪切位错震源模式。Abstract: The full moment tensor inversion is an important tool to monitor underground nuclear tests, thus has been employed in this paper to analyze the recorded data of seven underground nuclear tests and three natural earthquakes in the region of Eastern Kazakhstan. The results show that, except the obvious explosion source in underground nuclear tests, the components of double couple (DC) and compensate linear vector dipole (CLVD) also exist, in which the later plays a more important role and could be explained by the explosion-induced spall. In the case of natural earthquakes, the DC component is the main contribution, in consistent with the model of shear dislocation seismic source.
-
Key words:
- underground nuclear test /
- surface wave /
- Green function /
- full moment tensor inversion /
- CLVD source
-
地下核试验通常在严格保密的条件下进行,核试验场通常为人烟稀少的偏远地区,获取地下核试验的近场观测数据十分困难,从某种客观程度上讲,利用区域少量或单台观测数据分析和识别地下核试验具有十分重要的现实意义[1-2]。
通常情况下,可以采用忽略源区中的非线性效应等效力理论模型来描述震源[3-4],可以将震源定量描述为震源的断层面解(包括走向、滑动角、倾角)[5]。利用矩张量的形式分析震源机制,可以将观测资料、震源和地震波传播路径的相互关系归结为一种线性关系,随着理论地震图计算及有效分析能力的不断提高,利用复杂分层介质理论模型的震源机制反演方法取得许多重要进展。天然地震的反演研究以Jost等[6]的研究成果为代表。Wang等[1]通过进一步分析研究,得到了任意爆炸源和一个位错点源产生的地震波场表达式,以及其中包含的10个基本Green函数,其中包括爆炸源的基本Green函数。该地震波场表达式从理论上可以通过反演方式得到含有爆炸源成分的全矩张量解(full moment tensor)[5]。Minson等[7]和Dreger等[2]对Jost等[6]的包含爆炸Green函数的地震波场公式进一步修正,得到了能够准确反演爆炸源成分的公式,该公式采用Saikia等[8]定义的离散波数积分方法来研究理论Green函数。
本文中,采用基于Chen[9]和Yao等[10]的基于广义反、透射系数水平层状地球模中理论地震图的计算方法,对苏联的7次核爆炸及震中位置相近的天然地震进行反演,为探索核爆炸震源机制提供参考。
1. 方法原理
基于面波反演震源矩张量的理论地震波场可由矩张量M各分量作为权重的基本Green函数的线性组合来表示,早期工作中用于计算理论Green函数的介质模型比较简单,主要考虑地球表面的边界条件。本文中,采用文献[9-10]中理论地震图的计算方法进行分析,在圆柱坐标系中,理论位移公式为[5]:
{uz=Mxx[Gz,ss2cos(2θ)−Gz,dd3+Gz,ep3]+Myy[−Gz,ss2cos(2θ)−Gz,dd3+Gz,ep3]+Mzz[23Gz,dd+Gz,eq3]+Mxy[Gz,sssin(2θ)]+Mxz[−Gz,dscosθ]+Myz[−Gz,dssinθ]ur=Mxx[Gr,ss2cos(2θ)−Gr,dd3+Gr,ep3]+Myy[−Gr,ss2cos(2θ)−Gr,dd3+Gr,ep3]+Mzz[23Gr,dd+Gr,eq3]+Mxy[Gr,sssin(2θ)]+Mxz[−Gr,dscosθ]+Myz[−Gr,dssinθut=Mxx[Gt,ss2sin(2θ)]+Myy[−Gt,ss2sin(2θ)]+Mxy[−Gt,sscos(2θ)]+Mxz[Gt,dssinθ]+Myz[−Gt,dscosθ] (1) 式中:θ代表观测台站到震源的方位角;G为各个滑断层产生的Green函数,下标ss代表纯走滑断层(倾角90°,滑动角0°),ds代表纯倾滑断层(倾角90°,滑动角90°),dd代表 45°的斜滑断层(倾角45°,滑动角90°)[11],ep是纯爆炸源的Green函数;z, r, t分别代表垂向、径向和切向。式(1)的优点在于对所反演矩张量不设任何约束条件,不仅能反演双力偶成分(DC),还能准确地反演对角线成分(ISO)和线性偶极补偿源成分(CLVD)。
2. 实际应用
2.1 地下核爆炸震源矩张量反演
对Incorporated Research Institutions for Seismology (IRIS)乌鲁木齐站记录到质量较好的苏联的7次地下核试验进行分析,采用适合东哈萨克斯坦地区的Steven模型来计算理论Green函数[12],震中距取970km,震源时间函数采用光滑的阶跃函数[13]。图 1是对代号为88044事件记录数据进行处理的结果,处理顺序依次为:去均值并采用反卷积(图 1(a))、去倾校正(图 1(b)),同时均采用反卷积仪器传递函数;将数据旋转到z、r、t分量(图 1(c));利用Butterworth滤波器进行滤波,滤波周期为50s到100s (图 1(d));图中:BHE、BHN和BHZ分别代表东、北和竖直方向。
对基本Green函数进行相同周期范围的Butterworth滤波,利用式(1)对88044事件进行矩张量反演,与实际观测数据进行方差缩减(VR)计算,取方差缩减最大结果为最佳解[14-15]。一般情况下,地下核试验的爆炸装置埋深在1.5km以内,理论Green函数的震源深度为0.2~1.5km,步长为0.1km,地震矩统一取为M0≈1017N·m。震源矩张量反演结果为:
M=M0[0.7406−0.9489−0.1330−0.94890.84030.1197−0.13300.1197−0.3710] (2) 88044事件波形拟合结果如图 2所示。
采用Bowers等[16]的方法对矩张量进行分解,该分解方法的优点在于可以直接计算出ISO、DC、CLVD源成分的比重:
M=MISO+MDC+MCLVD (3) 矩张量M的特征值及其分解结果为[17]:
D=[−0.3868−0.15891.7556]=0.4033[111]+0.2279[−101]+1.1244[−0.5−0.51] (4) 对88094事件进行矩张量反演,过程与处理88044事件过程相同,反演得到的震源矩张量为:
M=M0[1.23461.6310−0.09231.63101.40600.0831−0.09230.0831−0.6308] (5) 88094事件波形拟合结果如图 3所示。矩张量M的特征值及其分解结果为:
D=[−0.67−0.272.95]=0.67[111]+0.40[−101]+1.88[−0.5−0.51] (6) 对88125事件反演得到的震源矩张量为:
M=M0[0.96411.2087−0.03781.20871.09110.0341−0.03780.0341−0.3848] (7) 88125事件波形拟合结果如图 4所示。矩张量M的特征值及其分解结果为:
D=[−0.40−0.172.25]=0.56[111]+0.23[−101]+1.46[−0.5−0.51] (8) 同理,对代号分别为88258、89189、88352、88317的事件数据进行了震源矩张量反演,对7次地下核试验的矩张量结果进行分析,分别得到EXP源(数值同ISO成分)、DC源和CLVD源在震源中所占比重,结果见表 1。
表 1 苏联7次地下核试验的震源矩张量解Table 1. Moment tensor of seven underground nuclear tests试验编号 EXP DC CLVD VR 88125 24.89% 10.22% 64.89% 87.77% 88094 22.71% 13.56% 63.73% 89.06% 88258 19.95% 13.83% 66.22% 83.72% 89189 28.16% 13.74% 58.10% 73.97% 88044 22.97% 12.98% 64.05% 85.26% 88352 35.09% 14.40% 50.51% 81.36% 88317 38.67% 7.98% 53.35% 78.91% 由表 1可以看出,7次地下核试验的震源矩张量反演结果中,均包含EXP源及矩张量迹为“零”的成分,地下核爆炸震源的球对称机制可以解释EXP源成分的存在。将迹为“零”的成分进一步分解为DC源和CLVD源,明显地看出,CLVD成分普遍较大,与Dreger等[18]给出的统计结果是一致的。爆炸导致的震源区应力破裂及层裂介质回落拍击的物理过程[19-22],可以很好地解释CLVD成分。震源区介质的各向异性以及路径中的构造应力释放等,均可以导致剪切位错成分的存在,从结果来看,DC源比重最小,而天然地震的震源机制主要为DC源,说明地下核试验的震源机制与天然地震是不同的。成分较大CLVD源,必然会体现在地下核爆炸激发的地震波波形数据中,其对角分量的偏差,可以直接激发S波,该波速度普遍低于上地幔S波速度,这样会在Moho面反射,形成能量很大的导波,继而导致地下核试验的波形数据中含有丰富的低频Lg波,实际地下核爆炸地震波数据中,可以普遍的观察到这一现象。
2.2 天然地震震源矩张量反演
采用同样的方法,对IRIS中乌鲁木齐台站记录到质量较好的苏联地下核试验场附近发生的3次天然地震事件的震源矩张量进行了反演,3次事件的基本参数如表 2所示,计算理论Green函数的震源深度采用Global CMT Project给出的结果。
表 2 苏联地下核试验场中天然地震的相关信息Table 2. Information about the three natural earthquakes事件代号 事件日期 震中位置 震源深度 震级/Mb 08029 2008-01-29 49.74N,83.48E 24.8km 4.6 08111 2008-04-20 45.97N,82.73E 21.9km 4.9 98193 1998-07-12 47.79N,82.78E 35.4km 5.2 对天然地震进行相同步骤的反演,08029事件反演结果为:
M=M0[−5.0908−6.67561.3793−6.67566.3099−1.26451.3793−1.2645−0.6965] (9) 08029事件波形拟合结果如图 5所示。矩张量M的特征值及其分解结果为:
D=[−8.2396−0.91359.6757]=0.1742[111]+7.3261[−101]+2.1754[−0.5−0.51] (10) 图 6为事件08111的波形拟合结果,震源矩张量反演结果为:
M=M0[0.0367−0.08031.5271−0.08030.0636−1.68231.5271−1.68230.3605] (11) 矩张量M的特征值及其分解结果为:
D=[−2.029−0.0312.5209]=0.1536[111]+1.998[−101]+0.3692[−0.5−0.51] (12) 事件98193的矩张量反演结果为:
M=M0[0.4216−0.00161.0429−0.00160.3435−1.07571.0429−1.07570.6834] (13) 98193事件波形拟合结果如图 7所示,其矩张量M的特征值及其分解结果为:
D=[−0.97320.38232.0393]=0.4828[111]+1.3555[−101]+0.201[−0.5−0.51] (14) 对于这3次天然地震的矩张量结果进行分析,同样得到EXP源、DC源和CLVD源在震源中所占比重,结果见表 3。
表 3 苏联地下核试验场附近天然地震震源矩张量解Table 3. Moment tensor of underground nuclear tests试验编号 EXP DC CLVD VR 08029 1.8% 75.72% 22.48% 86.69% 08111 6.1% 79.26% 14.64% 85.21% 98193 23.66% 66.47% 9.87% 83.33% 由上表的反演结果看出,相对地下核试验的震源矩张量反演结果来说,三次事件的震源中的DC成分是主要成分,天然地震的震源机制主要为剪切位错模式,这可以作为这3次事件中DC成分普遍大的合理解释。EXP、CLVD成分存在的原因是比较复杂的,对火山地震及深源地震震源的研究结果发现普遍存在CLVD成分,本文采用的一维模型,会忽略横向不均匀的效应。理论上来讲,算法也许会将EXP、DC以外的效应折合到CLVD成分中。
综上所述:在含有噪声、震相到时误差、地壳模型误差及震源深度误差等情况下,基于本文的公式的反演方法都能得到比较准确、稳定的结果[5, 23-24]。
3. 结论
对苏联实际地下核爆炸试验的反演结果表明,地下核爆炸的震源含有明显的爆炸源成分,非爆炸源成分主要是缘于构造应力释放和介质的不均匀性,CLVD源可看作是地下核爆炸中层裂震源机制的良好模型,其物理机制主要表现为层裂破裂介质的回落,此过程时间短,整体效果表现为无体积变化,CLVD源成分较大,稳定地分布在55%~65%之间,这种稳定的比重与地下核爆炸震源物理机制相关。3次发生在核试验场附近的天然地震震源机制反演结果表明,DC源占主要成分,符合天然地震震源为剪切位错模式的普遍结论。由于WMQ台站距试验场有近1000km,地震波的传播路径和地壳速度模型对反演结果会有影响,同时台站的记录数据可能受到较大的干扰因素,数据质量会受到影响,因此由一维地壳速度模型反演得到准确的结果是比较困难的。
-
表 1 苏联7次地下核试验的震源矩张量解
Table 1. Moment tensor of seven underground nuclear tests
试验编号 EXP DC CLVD VR 88125 24.89% 10.22% 64.89% 87.77% 88094 22.71% 13.56% 63.73% 89.06% 88258 19.95% 13.83% 66.22% 83.72% 89189 28.16% 13.74% 58.10% 73.97% 88044 22.97% 12.98% 64.05% 85.26% 88352 35.09% 14.40% 50.51% 81.36% 88317 38.67% 7.98% 53.35% 78.91% 表 2 苏联地下核试验场中天然地震的相关信息
Table 2. Information about the three natural earthquakes
事件代号 事件日期 震中位置 震源深度 震级/Mb 08029 2008-01-29 49.74N,83.48E 24.8km 4.6 08111 2008-04-20 45.97N,82.73E 21.9km 4.9 98193 1998-07-12 47.79N,82.78E 35.4km 5.2 表 3 苏联地下核试验场附近天然地震震源矩张量解
Table 3. Moment tensor of underground nuclear tests
试验编号 EXP DC CLVD VR 08029 1.8% 75.72% 22.48% 86.69% 08111 6.1% 79.26% 14.64% 85.21% 98193 23.66% 66.47% 9.87% 83.33% -
[1] WANG C Y, HERRMANN R B. A numerical study of P, SV, Sh-wave generation in a plane layered medium[J]. Bulletin of the Seismological Society of America, 1980, 70(4):1015-1036. [2] MINSON S E, DREGER D S. Stable inversions for complete moment tensors[J]. Geophysical Journal Internation, 2008, 174(2):585-592. doi: 10.1111/gji.2008.174.issue-2 [3] STUMP B W, JOHNSON L R. The determination of source properties by the linear inversion of seismograms[J]. Bulletin of the Seismological Society of America, 1977, 67:1489-1502. [4] AKI K, RICHARDS P G. Quantitative seismology, theory and methods[M]. San Francisco:W H Freeman, 1980. [5] 何永锋, 李锴, 刘炳灿, 等.基于面波数据的地下核爆炸的全元素矩张量反演方法[J].爆炸与冲击.2017, 37(5):945-950. http://www.bzycj.cn/CN/abstract/abstract9804.shtmlHE Yongfeng, LI Kai, LIU Bingcan, et al. Full moment tensor inversion method of underground nuclear explosions based on surface waves data[J]. Explosion and Shock Waves, 2017, 37(5):945-950. http://www.bzycj.cn/CN/abstract/abstract9804.shtml [6] JOST M L, HERRMANN R B. A student's guide to and review of moment tensor[J]. Seismological Research Letters, 1989, 60(2):37-57. [7] DREGER D, HELMBERGER D. Determination of source parameters at regional distances with three-components sparse network data[J]. Journal of Geophysical Research-Atmospheres, 1993, 98:8107-8125. doi: 10.1029/93JB00023 [8] SAIKIA C K. Modified frequency-wavenumber algorithm for regional seismograms using Filon's quadrature:modeling of Lg waves in eastern North America[J]. Geophysical Journal Internation, 1994, 118(1):142-158. doi: 10.1111/gji.1994.118.issue-1 [9] CHEN X F. A systematic and efficient method of computing normal modes for multilayered half-space[J]. Geophysical Journal Internation, 1993, 115(2):391-409. doi: 10.1111/gji.1993.115.issue-2 [10] YAO Z X, HARKRIDER D G. A generalized reflection-transmission coefficient matrix and discrete wavenumber method for synthetic seismograms[J]. Bulletin of the Seismological Society of America, 1983, 73(6):1685-1699. [11] LANGSTON C A. Source inversion of seismic waveforms:the Koyna, India, earthquakes of 13 September 1967[J]. Bulletin of the Seismological Society of America, 1981, 71(1):1-24. [12] MCLAUGHLIN K L, BARKER T G, DAY S M, et al. Effects of depth of burial on explosion and earthquake regional seismograms: regional discrimination and yield estimation: SSS_R_88_9844[R]. La Jolla, California, 1988. [13] DAY S M. Surface waves from underground explosions with spall:analysis of elastic and nolinear source models[J]. Bulletin of the Seismological Society of America, 1983, 73(1):247-264. [14] 林向东, 葛洪魁, 徐平, 等.近场全波形反演:芦山7.0级地震及余震矩张量解[J].地球物理学报, 2013, 56(12):4037-4047. doi: 10.6038/cjg20131209LIN Xiangdong, GE Hongkui, XU Ping, et al. Near field full waveform inversion:Lushan magnitude 7.0 earthquake and its aftershock moment tensor[J]. Chinese Journal of Geophysics, 2013, 56(12):4037-4047. doi: 10.6038/cjg20131209 [15] TEMPLETON D C, DREGER D S. Non-double couple earthquakes in the long valley volcanic region[J]. Bulletin of the Seismological Society of America, 2006, 96(1):69-79. doi: 10.1785/0120040206 [16] BOWERS D, HUDSON J A. Defining the scalar moment of a seismic source with a generak moment tensor[J]. Bulletin of the Seismological Society of America, 1999, 89(5):1390-1394. [17] FORD S R, DREGER D S, WALTER W R. Identifying isotropic events using an improved regional moment tensor inversion technique[J]. Journal of Geophysical Research, 2008, 114(B1). DOI: 10.1029/2008jb005743. [18] DREGER D, WOODS B. Regional distance seismic moment tensors of nuclear explosions[J]. Tectonophysics, 2002, 365(1):139-156. http://www.sciencedirect.com/science/article/pii/S0040195102003815 [19] MASSE R P. Review of seismic source models for underground nuclear explosions[J]. Bulletin of the Seismological Society of America, 1981, 71(4):1249-1268. [20] PATTON H J. Source models of the HARZER explosion from regional observations of fundamental-mode and higher mode surface waves[J]. Bulletin of the Seismological Society of America, 1988, 78(3), 1133-1157. [21] 何永锋, 陈晓非, 张海明.地下核爆炸Lg波的激发机制[J].地球物理学报, 2005, 48(2):367-372. doi: 10.3321/j.issn:0001-5733.2005.02.019HE Yongfeng, CHEN Xiaofei, ZHANG Haiming. The excitation of Lg wave by underground nuclear explosions[J]. Chinese Journal of Geophysics, 2005, 48(2):367-372. doi: 10.3321/j.issn:0001-5733.2005.02.019 [22] 何永锋, 赵克常, 张献兵.地下核爆炸地震波二次源特征[J].地球物理学报, 2012, 55(5):1742-1748. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201205031HE Yongfeng, ZHAO Kechang, ZHANG Xianbing, et al. The characteristic of the waveform from the second source induced by underground explosion[J]. Chinese Journal of Geophysics, 2012, 55(5):1742-1748. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201205031 [23] 郑建常, 陈运泰.基于Langston分解和Hilbert变换约束的区域偏量矩张量反演方法及应用[J].地震学报, 2012, 34(2):171-190. doi: 10.3969/j.issn.0253-3782.2012.02.005ZHENG Jianchang, CHEN Yuntai. Regional deviatoric moment tensor inversion based on Langston's decomposition and Hilbert transform constraints and its application[J]. Acta Seismologica Sinica, 2012, 34(2):171-190. doi: 10.3969/j.issn.0253-3782.2012.02.005 [24] 许力生, 陈运泰.震源深度误差对矩张量反演的影响[J].地震学报, 1997, 19(5):462-470. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199700067589 -