地下核试验通常在严格保密的条件下进行，核试验场通常为人烟稀少的偏远地区，获取地下核试验的近场观测数据十分困难，从某种客观程度上讲，利用区域少量或单台观测数据分析和识别地下核试验具有十分重要的现实意义[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]：

(1)

式中：θ代表观测台站到震源的方位角；G为各个滑断层产生的Green函数，下标ss代表纯走滑断层(倾角90°，滑动角0°)，ds代表纯倾滑断层(倾角90°，滑动角90°)，dd代表45°的斜滑断层(倾角45°，滑动角90°)[11]，ep是纯爆炸源的Green函数；z,r,t分别代表垂向、径向和切向。式(1)的优点在于对所反演矩张量不设任何约束条件，不仅能反演双力偶成分(DC)，还能准确地反演对角线成分(ISO)和线性偶极补偿源成分(CLVD) 。

分析下文数据，发现输血后不良反应组患者的各项指标均低于无不良反应组，差异有统计学意义（P＜0.05）；两组患者输血后的血常规指标比对，差异无统计学意义（P＞0.05）。

2 实际应用

2.1 地下核爆炸震源矩张量反演

由表1可以看出，7次地下核试验的震源矩张量反演结果中，均包含EXP源及矩张量迹为“零”的成分，地下核爆炸震源的球对称机制可以解释EXP源成分的存在。将迹为“零”的成分进一步分解为DC源和CLVD源，明显地看出， CLVD成分普遍较大，与Dreger等[18]给出的统计结果是一致的。爆炸导致的震源区应力破裂及层裂介质回落拍击的物理过程[19-22]，可以很好地解释CLVD成分。震源区介质的各向异性以及路径中的构造应力释放等， 均可以导致剪切位错成分的存在，从结果来看，DC源比重最小，而天然地震的震源机制主要为DC源，说明地下核试验的震源机制与天然地震是不同的。成分较大CLVD源，必然会体现在地下核爆炸激发的地震波波形数据中，其对角分量的偏差，可以直接激发S波，该波速度普遍低于上地幔S波速度，这样会在Moho面反射，形成能量很大的导波，继而导致地下核试验的波形数据中含有丰富的低频Lg波，实际地下核爆炸地震波数据中，可以普遍的观察到这一现象。

图1 88044事件波形处理 Fig.1 Processing procedure of 88044 event

想想看，未来你不需要去学校上学，只要坐在家里戴上VR眼镜，就可以畅游在知识的海洋，不会有老师枯燥的讲解让你昏昏欲睡，你可以潜入古墓拜访木乃伊，你可以到罗马斗兽场目睹角斗士的殊死拼搏，你可以在长平战场上体验一回被活埋的滋味……

(2)

项目工程造价主要是指整个工程建设需要投入的费用总额，也可将其进一步细化为2个部分：（1）工程项目投资。基于投资者视角进行分析，造价即一个工程建设应花费的总额。（2）基于承包者视角进行分析，造价是建设工程发包承包实施结果作为商品进入到流通领域，能够体现出的货币量，可简单理解为是工程承发包价格[1]。

88044事件波形拟合结果如图2所示。

图2 88044事件观测数据与理论数据比较 Fig.2 Comparison between data and synthesis of 88044 event

景德镇先辈们留下辉煌的陶瓷产业需要我们延续她的光芒，保留这一伟大的陶瓷产业。去年，中共中央办公厅、国务院办公厅联合印发了《关于实施中华优秀传统文化传承发展工程的意见》，意见指出：“对于传统中华文脉，全力提升人民群众文化素养，维护国家文化安全，增强国家文化软实力，推进国家治理体系和统治能力的现代化，具有重要的意义，”陶瓷文化也是中华传统文化的重要组成部分。中国的陶瓷史可以说就是中国的文化史，欧美各国一度就是通过了解中国的瓷器来了解中国的文化，所以把中国叫做“chian”，也就是景德镇昌南。前1000年景德镇陶瓷发展就是在柴窑窑火的烧制下走过，景德镇陶瓷伴随着柴窑窑火也必将走向下一个1000年。

M=MISO+MDC+MCLVD

(3)

矩张量M的特征值及其分解结果为[17]：

(4)

对88094事件进行矩张量反演，过程与处理88044事件过程相同，反演得到的震源矩张量为：

(5)

88094事件波形拟合结果如图3所示。矩张量M的特征值及其分解结果为：

(6)

图3 88094事件观测数据与理论数据比较 Fig.3 Comparison between data and synthesis of 88094 event

对88125事件反演得到的震源矩张量为:

(7)

88125事件波形拟合结果如图4所示。矩张量M的特征值及其分解结果为：

采用Bowers等[16]的方法对矩张量进行分解，该分解方法的优点在于可以直接计算出ISO、DC、CLVD源成分的比重：

(8)

图4 88125事件观测数据与理论数据比较 Fig.4 Comparison between data and synthesis of 88125 event

采用同样的方法，对IRIS中乌鲁木齐台站记录到质量较好的苏联地下核试验场附近发生的3次天然地震事件的震源矩张量进行了反演，3次事件的基本参数如表2所示，计算理论Green函数的震源深度采用Global CMT Project给出的结果。

转眼上元在即。今年不同于往常，只因灯节毗邻皇后生辰，圣上恩准举办灯会。灯会专辟禁区，供皇族和品级高的皇亲国戚、股肱之臣观赏。

表1 苏联7次地下核试验的震源矩张量解 Table 1 Moment tensor of seven underground nuclear tests

试验编号EXPDCCLVDVR8812524.89%10.22%64.89%87.77%8809422.71%13.56%63.73%89.06%8825819.95%13.83%66.22%83.72%8918928.16%13.74%58.10%73.97%8804422.97%12.98%64.05%85.26%8835235.09%14.40%50.51%81.36%8831738.67%7.98%53.35%78.91%

对Incorporated Research Institutions for Seismology (IRIS)乌鲁木齐站记录到质量较好的苏联的7次地下核试验进行分析，采用适合东哈萨克斯坦地区的Steven模型来计算理论Green函数[12]，震中距取970 km，震源时间函数采用光滑的阶跃函数[13]。图1是对代号为88044事件记录数据进行处理的结果，处理顺序依次为：去均值并采用反卷积(图1(a))、去倾校正(图1(b))，同时均采用反卷积仪器传递函数；将数据旋转到z、r、t分量(图1(c))；利用Butterworth滤波器进行滤波，滤波周期为50 s到100 s (图1(d))；图中：BHE、BHN和BHZ分别代表东、北和竖直方向。

2.2 天然地震震源矩张量反演

同理，对代号分别为88258、89189、88352、88317的事件数据进行了震源矩张量反演，对7次地下核试验的矩张量结果进行分析，分别得到EXP源(数值同ISO成分)、DC源和CLVD源在震源中所占比重，结果见表1。

对基本Green函数进行相同周期范围的Butterworth滤波，利用式(1)对88044事件进行矩张量反演，与实际观测数据进行方差缩减(VR)计算，取方差缩减最大结果为最佳解[14-15]。一般情况下，地下核试验的爆炸装置埋深在1.5 km以内，理论Green函数的震源深度为0.2～1.5 km，步长为0.1 km，地震矩统一取为M0≈1017 N·m。震源矩张量反演结果为：

表2 苏联地下核试验场中天然地震的相关信息 Table 2 Information about the three natural earthquakes

事件代号事件日期震中位置震源深度震级/Mb080292008-01-2949.74N,83.48E24.8km4.6081112008-04-2045.97N,82.73E21.9km4.9981931998-07-1247.79N,82.78E35.4km5.2

对天然地震进行相同步骤的反演，08029事件反演结果为：

(9)

08029事件波形拟合结果如图5所示。矩张量M的特征值及其分解结果为：

(10)

图5 08029事件观测数据与理论数据比较 Fig.5 Comparison between data and synthesis of 08029 event

图6为事件08111的波形拟合结果，震源矩张量反演结果为：

(11)

矩张量M的特征值及其分解结果为：

(12)

图6 08111事件观测数据与理论数据比较 Fig.6 Comparison between data and synthesis of 08111 event

事件98193的矩张量反演结果为：

(13)

98193事件波形拟合结果如图7所示，其矩张量M的特征值及其分解结果为：

针对百米林带植物群落中封闭草地型空间占比较大、乔木种植过密的问题。建议梳理群落空间层次,对一些植物长势过于杂乱、种植密度过高的群落进行适当的抽稀处理,去除长势较差的植物,使群落层次更加清晰,在提升景观视觉效果的同时也增加可游憩面积。

(14)

图7 98193事件观测数据与理论数据比较 Fig.7 Comparison between data and synthesis of 98193 event

对于这3次天然地震的矩张量结果进行分析，同样得到EXP源、DC源和CLVD源在震源中所占比重，结果见表3。

表3 苏联地下核试验场附近天然地震震源矩张量解 Table 3 Moment tensor of underground nuclear tests

试验编号EXPDCCLVDVR080291.8%75.72%22.48%86.69%081116.1%79.26%14.64%85.21%9819323.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台站距试验场有近1 000 km，地震波的传播路径和地壳速度模型对反演结果会有影响，同时台站的记录数据可能受到较大的干扰因素，数据质量会受到影响，因此由一维地壳速度模型反演得到准确的结果是比较困难的。

参考文献:

[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.

[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.

HE 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.

[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.

[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.

[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.

[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.

LIN 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.

[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.

[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.

[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.

HE Yongfeng, CHEN Xiaofei, ZHANG Haiming. The excitation of Lg wave by underground nuclear explosions[J]. Chinese Journal of Geophysics, 2005,48(2):367-372.

[22] 何永锋,赵克常,张献兵.地下核爆炸地震波二次源特征[J].地球物理学报,2012,55(5):1742-1748.

HE 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.

[23] 郑建常,陈运泰.基于Langston分解和Hilbert变换约束的区域偏量矩张量反演方法及应用[J].地震学报,2012,34(2):171-190.

ZHENG 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.

[24] 许力生,陈运泰.震源深度误差对矩张量反演的影响[J].地震学报,1997,19(5):462-470.