荐言献策

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

2016年8月20日唐山震群小震震源机制

王宁 蔡璐 王莉婵 王亚玲 郭垚嘉 刘丽

王宁,蔡璐,王莉婵, 等. 2016年8月20日唐山震群小震震源机制[J]. 华北地震科学,2020, 38(3):49-55. doi:10.3969/j.issn.1003−1375.2020.03.009
引用本文: 王宁,蔡璐,王莉婵, 等. 2016年8月20日唐山震群小震震源机制[J]. 华北地震科学,2020, 38(3):49-55. doi:10.3969/ j.issn.1003−1375.2020.03.009
WANG Ning,CAI Lu,WANG Lichan, et al. Focal Mechanism of Small Earthquakes in Tangshan Earthquake Swarm on August 20, 2016[J]. North China Earthqauke Sciences,2020, 38(3):49-55. doi:10.3969/j.issn.1003−1375.2020.03.009
Citation: WANG Ning,CAI Lu,WANG Lichan, et al. Focal Mechanism of Small Earthquakes in Tangshan Earthquake Swarm on August 20, 2016[J]. North China Earthqauke Sciences,2020, 38(3):49-55. doi:10.3969/j.issn.1003−1375.2020.03.009

2016年8月20日唐山震群小震震源机制

doi: 10.3969/j.issn.1003-1375.2020.03.009
基金项目: 河北省地震局三结合课题(DZ20190422039);星火计划面上项目(DZ20190422038);中国地震局震情跟踪课题(2019010105);河北省地震局星火项目(DZ20180319008)
详细信息
    作者简介:

    王宁(1988—),女,河北邯郸人,工程师,主要从事地震监测方面研究. E-mail:819087955@qq.com

  • 中图分类号: P315.6

Focal Mechanism of Small Earthquakes in Tangshan Earthquake Swarm on August 20, 2016

  • 摘要: 运用新近发展的适用于求解小震震源机制的广义极性振幅技术(GPAT),求解2016年8月20日唐山震群ML2.0以上地震事件的震源机制,结果显示:①震群的震源机制以正断型为主,占总数的69.5%;②对于震群中小于ML3.5的地震,GPAT得到的矩震级MW大于台网观测报告常用的ML震级,事件震级越小,差距越明显,最大震级偏差为1.0;③GPAT计算深度能够反映震源的真实深度。总之,GPAT方法在唐山震群震源机制反演中效果良好。
  • 图  1  震群主震震中与台站分布图

    注:▲为台站,显示台站代码标注的为参与计算的台站;★为主震震中位置

    图  2  震群震中分布图

    图  3  2016-08-21 17:46:51.7GPAT震源机制反演结果

    注:值为正表示截取的震相周期内绝对值最大点的符号为正,反之为负

    图  4  震源机制解分布图

    注:红色为正断型和正断兼走滑分量类型;黑色为走滑型;绿色为逆断型和逆断兼走滑分量类型;蓝色为无法确定型。

    图  5  矩震级MW与地方震级ML对比图

    图  6  深度对比图

    表  1  华北速度模型

    VP/(km/s)VS/(km/s)层厚/km
    16.013.51515
    26.513.80718
    38.014.68
    下载: 导出CSV

    表  2  震群震源机制解及断层类型

    编号发震时刻经度/(ºE)纬度/(ºN)深度/kmML节面Ⅰ/(°)节面Ⅱ/(°)类型
    走向倾向滑动角走向倾向滑动角
    12016-08-21 08:55:16.3118.35239.70262.130182−11320323−9U
    22016-08-21 17:15:35.1118.34439.69973.626638−11111355−74NF
    32016-08-21 17:46:51.7118.34739.70072.325551−11711546−60NF
    42016-08-21 18:09:59.3118.34639.69573.324357−12411346−50NF
    52016-08-21 19:45:02.4118.35139.69562.923180−10310416−38NF
    62016-08-21 19:46:04.2118.35439.69962.62829−12713982−87U
    72016-08-21 19:53:26.7118.35939.70362.525125−15814181−67U
    82016-08-21 20:46:12.6118.36339.68962.526251−11211544−65NF
    92016-08-21 20:47:59.0118.34039.67662.127433−16017082−60U
    102016-08-22 04:07:30.1118.35539.69973.522254−16012074−38SS
    112016-08-22 05:50:32.3118.34739.70062.927840−10712052−76NF
    122016-08-22 15:15:29.1118.35139.69453.326359−12313444−48NF
    132016-08-22 15:17:53.3118.34239.67762.121664−11210823−19NS
    142016-08-22 22:22:00.5118.34739.69363.025043−13812763−55NF
    152016-08-23 02:24:54.2118.34539.68062.226051−13413857−49NF
    162016-08-23 12:51:46.1118.35139.70262.222576−13412146−20NS
    172016-08-23 16:21:39.4118.34939.69252.129051−13116454−51NF
    182016-08-23 21:10:49.1118.33539.687172.623878−16214472−13SS
    192016-08-23 22:22:54.7118.34239.69672.223531−1138162−77NF
    202016-08-24 10:08:25.6118.36539.70153.124832−1018159−83NF
    212016-08-25 11:59:23.7118.36339.69672.926355−10210337−74NF
    222016-08-27 15:52:35.2118.33839.70772.127226−747465−98NF
    232016-08-29 13:02:54.0118.34539.68562.326149−13313657−52NF
    242016-08-29 13:12:55.3118.34739.67652.325938−907952−90NF
    252016-08-29 17:41:38.7118.35239.68972.626346−10810847−72NF
    262016-08-30 20:35:04.1118.35539.68562.429536−13817071−66NF
    272016-08-31 18:15:27.7118.34739.68972.626446−11511849−66NF
    282016-09-01 01:01:15.9118.34339.69172.23024−558087−93U
    292016-09-01 08:51:20.7118.34939.69672.326427−847763−93NF
    302016-09-05 03:39:19.9118.35739.69152.624263−1049130−64NF
    312016-09-05 04:33:22.3118.36239.69762.526443−11011350−72NF
    322016-09-07 19:58:31.6118.35739.69752.128129−14515974−66NF
    332016-09-10 18:19:15.4118.35739.69663.025986−11316123−9U
    342016-09-10 23:35:59.2118.32939.67962.135255901724690TS
    352016-09-10 18:09:37.3118.35139.70464.326733−13013265−67NF
    362016-09-10 18:13:00.8118.35639.69352.628540−9611350−85NF
    372016-09-11 06:57:02.9118.35739.68572.0338703423558156SS
    382016-09-11 07:20:44.1118.32439.68372.031376−204871−165SS
    392016-09-11 07:41:31.7118.33139.68672.12656615656826SS
    402016-09-11 14:41:39.6118.33139.68672.023568−13614384−6SS
    412016-09-11 16:44:13.2118.34739.70072.026025−25167−153U
    422016-09-11 22:05:05.2118.34639.69472.11347714489167SS
    432016-09-13 09:59:10.3118.34839.68262.924358−1069235−66NF
    442016-09-13 10:50:09.0118.34739.69152.426320−494075−103NF
    452016-09-14 20:00:17.4118.34539.68062.524063−957127−80NF
    462016-09-16 03:34:42.5118.35639.69482.726264−12413942−41NF
    472016-09-17 12:06:37.3118.35139.68072.126628−17317288−62U
    482016-09-17 16:18:33.2118.34939.69672.829277−845424−146NF
    492016-09-17 16:20:04.1118.34639.67972.123549−1118645−67NF
    502016-09-19 07:18:06.4118.34539.70172.923266−1007626−68NF
    512016-09-19 07:31:38.4118.35139.68352.723547−1077946−73NF
    522016-09-20 11:28:38.4118.34439.68362.823869−1039025−60NF
    532016-09-21 21:57:55.0118.34839.70063.226442−11912154−66NF
    542016-09-24 09:09:35.6118.35339.67982.224664−17620486−26SS
    552016-09-29 14:22:36.6118.35339.68653.025954−13714057−45NF
    562016-10-01 12:39:37.2118.35039.68582.129260−14917187−85NF
    572016-10-05 23:27:40.2118.35239.70573.526632−12612765−70NF
    582016-10-07 16:32:04.7118.34839.68852.127643−10211248−79NF
    592016-10-19 04:49:01.1118.34939.67372.22957−14617186−84U
    下载: 导出CSV
  • [1] Helmberger D V, Engen G R. Modeling the long-period body waves from shallow earthquakes at regional ranges[J]. Bulletin of the Seismological Society of America, 1980, 70(5): 1699-1714.
    [2] Zhao L S, Helmberger D V. Source estimation from broadband regional seismograms[J]. Bulletin of the Seismological Society of America, 1994, 84(1): 91-104.
    [3] Kisslinger C. Evaluation of S to P amplitude rations for determining focal mechanisms from regional network observations[J]. Bulletin of the Seismological Society of America, 1980, 70(4): 999-1014.
    [4] 梁尚鸿, 李幼铭, 束沛镒, 等. 利用区域地震台网、振幅比资料测定小震震源参数[J]. 地球物理学报, 1984, 27(3): 249-257. doi:  10.3321/j.issn:0001-5733.1984.03.005
    [5] Hardebeck J L, Shearer P M. Using S/P amplitude ratios to constrain the focal mechanisms of small earthquakes[J]. Bulletin of the Seismological Society of America, 2003, 93(6): 2434-2444. doi:  10.1785/0120020236
    [6] 严川, 许力生. 一种地方与区域地震震源机制反演技术: 广义极性振幅技术(一)——原理与数值实验[J]. 地球物理学报, 2014, 57(8): 2555-2572. doi:  10.6038/cjg20140816
    [7] 严川, 许力生, 张旭, 等. 一种地方与区域地震震源机制反演技术: 广义极性振幅技术(二)——对实际震例的应用[J]. 地球物理学报, 2015, 58(10): 3601-3614. doi:  10.6038/cjg20151014
    [8] 惠少兴, 严川, 王平, 等. 陕西地区小微震震源机制研究[J]. 地震工程学报, 2018, 40(2): 288-293. doi:  10.3969/j.issn.1000-0844.2018.02.288
    [9] Wang R J. A simple orthonormalization method for stable and efficient computation of Green's functions[J]. Bulletin of the Seismological Society of America, 1999, 89(3): 733-741.
    [10] Zoback M L. First-and second-order patterns of stress in the lithosphere: The World Stress Map Project[J]. Journal of Geophysical Research, 1992, 97(B8): 11703-11728. doi:  10.1029/92JB00132
    [11] Kanamori H. The energy release in great earthquakes[J]. Journal of Geophysical Research, 1977, 82(20): 2981-2987. doi:  10.1029/JB082i020p02981
    [12] Hanks T C, Kanamori H. A moment magnitude scale[J]. Journal of Geophysical Research, 1979, 84(B5): 2348-2350. doi:  10.1029/JB084iB05p02348
    [13] 陈运泰, 刘瑞丰. 地震的震级[J]. 地震地磁观测与研究, 2004, 25(6): 1-12. doi:  10.3969/j.issn.1003-3246.2004.06.001
    [14] 王宁, 蔡璐, 李雪英, 等. 2016年8月唐山震群视应力变化特征[J]. 地震地磁观测与研究, 2018, 39(3): 52-58. doi:  10.3969/j.issn.1003-3246.2018.03.008
  • [1] 朱音杰, 刘檀, 赵英萍, 刘新, 李冬圣, 丁成.  CAP方法反演2016年唐山ML4.3地震震源机制解 . 华北地震科学, 2017, 35(1): 50-55. doi: 10.3969/j.issn.1003-1375.2017.01.008
    [2] 李铂, 崔鑫, 苗庆杰, 穆娟.  CAP方法反演乳山震群3.0级以上地震震源机制解 . 华北地震科学, 2016, 34(3): 14-19. doi: 10.3969/j.issn.1003-1375.2016.03.003
    [3] 孙晴, 李守勇.  利用小震震源机制解反演濮阳地区应力场 . 华北地震科学, 2015, 33(4): 14-19. doi: 10.3969/j.issn.1003-1375.2015.04.003
    [4] 夏乱保, 孙佩卿.  红山地震台近震震级ML的台基校正 . 华北地震科学, 1989, 7(2): 71-78.
    [5] 韦士忠, 石汝斌, 李玉萍, 赵英芳.  用选频谱研究唐山余震序列较大地震前小震波谱和震源参数的某些特征 . 华北地震科学, 1989, 7(1): 43-51.
    [6] 平建军.  单台波速比对邢台5.8级中强余震异常显示的分析 . 华北地震科学, 1987, 5(4): 72-78.
    [7] 靳雅敏, 于新昌, 扈志刚.  由小震应力场的动态变化估计河北及邻近地区的地震危险性 . 华北地震科学, 1987, 5(1): 83-89.
    [8] 秦保燕, 徐纪人, 翟文杰.  南北地震带北段现今地壳应力场讨论 . 华北地震科学, 1986, 4(1): 1-15.
    [9] 吕培苓, 刘连柱.  1981年邢台5.8级强余震前后小地震的应力降估计 . 华北地震科学, 1986, 4(2): 26-32.
    [10] 漆家福.  由微(小)震综合断层面解推测林县地区的现代应力场 . 华北地震科学, 1985, 3(2): 41-45.
  • 加载中
图(6) / 表 (2)
计量
  • 文章访问数:  1775
  • HTML全文浏览量:  929
  • PDF下载量:  23
出版历程
  • 收稿日期:  2019-08-07
  • 网络出版日期:  2020-08-05
  • 刊出日期:  2020-07-01

2016年8月20日唐山震群小震震源机制

doi: 10.3969/j.issn.1003-1375.2020.03.009
    基金项目:  河北省地震局三结合课题(DZ20190422039);星火计划面上项目(DZ20190422038);中国地震局震情跟踪课题(2019010105);河北省地震局星火项目(DZ20180319008)
    作者简介:

    王宁(1988—),女,河北邯郸人,工程师,主要从事地震监测方面研究. E-mail:819087955@qq.com

  • 中图分类号: P315.6

摘要: 运用新近发展的适用于求解小震震源机制的广义极性振幅技术(GPAT),求解2016年8月20日唐山震群ML2.0以上地震事件的震源机制,结果显示:①震群的震源机制以正断型为主,占总数的69.5%;②对于震群中小于ML3.5的地震,GPAT得到的矩震级MW大于台网观测报告常用的ML震级,事件震级越小,差距越明显,最大震级偏差为1.0;③GPAT计算深度能够反映震源的真实深度。总之,GPAT方法在唐山震群震源机制反演中效果良好。

English Abstract

王宁,蔡璐,王莉婵, 等. 2016年8月20日唐山震群小震震源机制[J]. 华北地震科学,2020, 38(3):49-55. doi:10.3969/j.issn.1003−1375.2020.03.009
引用本文: 王宁,蔡璐,王莉婵, 等. 2016年8月20日唐山震群小震震源机制[J]. 华北地震科学,2020, 38(3):49-55. doi:10.3969/ j.issn.1003−1375.2020.03.009
WANG Ning,CAI Lu,WANG Lichan, et al. Focal Mechanism of Small Earthquakes in Tangshan Earthquake Swarm on August 20, 2016[J]. North China Earthqauke Sciences,2020, 38(3):49-55. doi:10.3969/j.issn.1003−1375.2020.03.009
Citation: WANG Ning,CAI Lu,WANG Lichan, et al. Focal Mechanism of Small Earthquakes in Tangshan Earthquake Swarm on August 20, 2016[J]. North China Earthqauke Sciences,2020, 38(3):49-55. doi:10.3969/j.issn.1003−1375.2020.03.009
    • 震源机制可以用来反映震源区域在地震发生前后的动力学运动过程,准确求解的震源机制可以用来判断震源类型,揭示地震区域应力场运动特征,并判定地震断层破裂面性质及其孕震构造。因此,震源机制的求解极为重要。震源机制反演方法很多,从最早1980年Helmberger等[1]发现Pnl波可以用于确定震源机制开始,国内外学者开展了很多研究工作,并取得了丰富的成果。关于中强震的震源机制的确定,现阶段运用较多并被接受的是改进的CAP方法[2]。而小震震源机制的确定,一直是学术界的难点。小震发生能量小,高频信号衰减快速,波形触发台站较少,其本身的局限性决定了其震源机制求解困难。即便如此,在小震求解方面的研究一直没有停止过。如利用P波初动和振幅比方法[3-5]确定小震震源机制等,虽然取得了一些成果,但现阶段并没有被广泛的认可。针对小震研究现状,本文采用了近几年被成功运用[6-8]的广义极性振幅技术(GPAT)对2016年唐山震群小震震源机制反演进行尝试,以期为唐山地区小震震源机制的研究提供一条途径。

      唐山地区,地震构造复杂,活动断裂发育,是华北地震高发区,尤其是1976年唐山大地震之后,它一直是河北省重点地震监测区域。近几年来唐山地区地震发生尤其活跃,并以2016年8月唐山震群最为典型。本文以此震群ML2.0以上地震事件波形资料为基础,利用广义极性振幅技术(GPAT)反演这些事件的震源机制解,通过分析其相关特征,对判定震群的孕震机理、区域应力场相关特征等具有重要意义。

    • 广义极性振幅法(GPAT)是严川提出来的一种求解小震震源机制的非线性反演方法,是现有初动-振幅类方法的一般化[6]。此方法以震相初动极性和振幅为研究对象,并运用网格搜索技术来进行震源机制反演。

      一般地,地震产生的地表位移可以表示为:

      $$ {\rm{u}}\left( {x,\;t} \right) = {u_j}{e_j}, $$ (1)
      $$ {u_j} = u_j^{\rm{P}} + u_j^{\rm{S}} + u_j^{\rm{F}} + u_j^{\rm{O}},\;\;{{j}} = 1,2,3 $$ (2)

      式中:j=1,2,3,分别代表EW、NS和UD三分向;P、S、F、O分别表示P波、S波、面波和其他震相。假设2个行矢量:

      $$ {v_1} = \left[ {{a_{ijk}}} \right],\;{v_2} = \left[ {{b_{ijk}}} \right];\;i = 1,2, \cdots N,\;j = 1,2,3,\;k = 1,2,3 $$ (3)

      式中:i表示不同的台站,j为不同的分量,k分别代表直达P波、S波和面波;行矢量v1的元素a是带有极性的振幅绝对值最大值,行矢量v2的元素b代表P波初动极性。

      按照式(3),对于给定的震源产生的波场可以写出类似的矢量:

      $$ {v_1}^{\prime} = \left[ {a_{ijk}^{'}} \right],\;{v_2}^{\prime} = \left[ {b_{ijk}^{'}} \right] $$ (4)
      $$ {\text{令}},\;\eta = \left[ {\begin{array}{*{20}{l}} {{v_1}}&{\omega {v_2}} \end{array}} \right] = {\left[ {\begin{array}{*{20}{l}} {{a_{ijk}}}&{\omega {b_{ijk}}} \end{array}} \right]^{\rm{T}}} = \left[ {\begin{array}{*{20}{l}} {{a_m}} \end{array}} \right] $$ (5)
      $$ {\eta ^\prime } = \left[ {\begin{array}{*{20}{l}} {v_1^{'}}&{\omega v_2^{'}} \end{array}} \right] = {\left[ {\begin{array}{*{20}{l}} {a_{ijk}^{'}}&{\omega b_{ijk}^{'}} \end{array}} \right]^{\rm{T}}} = \left[ {{\beta _m}} \right] $$ (6)

      式中:ɑm为观测矢量元素;βm为合成矢量元素;T为行矢量向列矢量的转置;w为振幅信息与P波初动信息的相对权重值。

      假定震源的震源机制标量、地震矩、震源位置与实际震源的完全相同,则有$\eta $=$\eta '$;而如果假定只有标量地震矩不同但其他参数与实际震源的相同,那么$\eta $$\eta '$线性相关。设定ρ为理论震源与实际震源参数的相关系数,则

      $$ \rho = \frac{{\displaystyle\sum\limits_m^M {\left( {{\alpha _m} - \bar \alpha } \right)} \left( {{\beta _m} - \bar \beta } \right)}}{{\sqrt {\displaystyle\sum\limits_m^M {{{\left( {{\alpha _m} - \bar \alpha } \right)}^2}} } \sqrt {\displaystyle\sum\limits_m^M {{{\left( {{\beta _m} - \bar \beta } \right)}^2}} } }} $$ (7)

      式中:$\overline a $$\overline b $分别表示ɑmβm的平均值。

      另外,关于w相对权重的取值问题定义2个列向量XY

      $$ X = {\left[ {\begin{array}{*{20}{l}} {{a_{ijk}}}&{\omega {b_{ij}}} \end{array}} \right]^{\rm{T}}} = \left[ {\begin{array}{*{20}{l}} {{x_m}} \end{array}} \right] $$ (8)
      $$ Y = {\left[ {\begin{array}{*{20}{l}} {{a_{ijk}}}&{ - \omega {b_{ij}}} \end{array}} \right]^{\rm{T}}} = \left[ {\begin{array}{*{20}{l}} {{y_m}} \end{array}} \right] $$ (9)

      两者的相关系数为:

      $$ \kappa = \frac{{\displaystyle\sum\limits_m^M {\left( {{x_m} - \bar x} \right)} \left( {{y_m} - \bar y} \right)}}{{\sqrt {\displaystyle\sum\limits_m^M {{{\left( {{x_m} - \bar x} \right)}^2}} } \sqrt {\displaystyle\sum\limits_m^M {{{\left( {{y_m} - \bar y} \right)}^2}} } }} $$ (10)

      如果w=0,则k=1,表明极性不参与计算,只有振幅参与计算;如果w→+∞,则k→-1,表明振幅对k值的影响很小。由此,根据w的取值不同,振幅和初动极性发挥不同的作用。

      根据以上公式可知,βm的确定受参与计算的事件震源机制和震源位置的影响。如果震中位置确定,则受震源机制与震源深度影响,此时ρ是关于震源机制、震源深度的函数:

      $$\rho = \rho (\psi ,\theta ,\lambda ,d)$$ (11)

      式中:ψθλd表示断层走向、倾角、滑动角、震源深度。此时震源机制反演变成非线性求解ρ$ \to $1的问题。

    • 通过首都圈和“十五”测震台网的建设和改造,现河北省测震台网由1个台网中心、71个区域数字地震台站组成,台网中心汇集了省内所有台站数据,以及国家测震台网中心接收周边98个台站(北京29个、天津31个、山东14个、山西9个、内蒙古5个、河南5个、辽宁5个)的数据,共有168个台站的实时波形数据。基于台站布设分布情况,全省范围内地震监测能力达到ML2.0,有些地区已高达ML1.0。

      本文使用了河北台网记录的2016年8月唐山震群资料(图1),并对每一次参与计算的地震事件进行了仔细的筛选。反演震源机制时按照震中距<100 km、波形清晰、较高信噪比、P波震相记录应多于6个台站、排除重叠地震的条件进行筛选,选出59个(图2)符合条件的天然地震事件,其中ML2.0~2.9地震49个,ML3.0以上地震10个。

      图  1  震群主震震中与台站分布图

      图  2  震群震中分布图

    • 唐山震群波形资料震中距范围选择在0~100 km。此次唐山震群震级较小,波形记录的频率较高,故将震群的滤波范围设置在1~4 Hz之间,采样频率确定为滤波频率上限的5倍。格林函数库计算的方法为反射-折射率法[9],深度范围为0~35 km,深度间隔为1 km。计算格林函数时选用的速度模型为河北台网现阶段在用的地壳速度模型—华北速度模型(表1)。

      表 1  华北速度模型

      VP/(km/s)VS/(km/s)层厚/km
      16.013.51515
      26.513.80718
      38.014.68
    • 根据华北速度模型和近震台站记录,利用GPAT方法反演得到震群的最佳震源机制解(表2)。图3为地震2016—08—21 17:46:51.7实际观测数据与理论合成数据的对比结果,反演结果中观测数据初动符号与合成数据初动符号极性完全一致,观测数据振幅与合成数据振幅相关性很好,可见震源机制解的可信赖度很高。

      表 2  震群震源机制解及断层类型

      编号发震时刻经度/(ºE)纬度/(ºN)深度/kmML节面Ⅰ/(°)节面Ⅱ/(°)类型
      走向倾向滑动角走向倾向滑动角
      12016-08-21 08:55:16.3118.35239.70262.130182−11320323−9U
      22016-08-21 17:15:35.1118.34439.69973.626638−11111355−74NF
      32016-08-21 17:46:51.7118.34739.70072.325551−11711546−60NF
      42016-08-21 18:09:59.3118.34639.69573.324357−12411346−50NF
      52016-08-21 19:45:02.4118.35139.69562.923180−10310416−38NF
      62016-08-21 19:46:04.2118.35439.69962.62829−12713982−87U
      72016-08-21 19:53:26.7118.35939.70362.525125−15814181−67U
      82016-08-21 20:46:12.6118.36339.68962.526251−11211544−65NF
      92016-08-21 20:47:59.0118.34039.67662.127433−16017082−60U
      102016-08-22 04:07:30.1118.35539.69973.522254−16012074−38SS
      112016-08-22 05:50:32.3118.34739.70062.927840−10712052−76NF
      122016-08-22 15:15:29.1118.35139.69453.326359−12313444−48NF
      132016-08-22 15:17:53.3118.34239.67762.121664−11210823−19NS
      142016-08-22 22:22:00.5118.34739.69363.025043−13812763−55NF
      152016-08-23 02:24:54.2118.34539.68062.226051−13413857−49NF
      162016-08-23 12:51:46.1118.35139.70262.222576−13412146−20NS
      172016-08-23 16:21:39.4118.34939.69252.129051−13116454−51NF
      182016-08-23 21:10:49.1118.33539.687172.623878−16214472−13SS
      192016-08-23 22:22:54.7118.34239.69672.223531−1138162−77NF
      202016-08-24 10:08:25.6118.36539.70153.124832−1018159−83NF
      212016-08-25 11:59:23.7118.36339.69672.926355−10210337−74NF
      222016-08-27 15:52:35.2118.33839.70772.127226−747465−98NF
      232016-08-29 13:02:54.0118.34539.68562.326149−13313657−52NF
      242016-08-29 13:12:55.3118.34739.67652.325938−907952−90NF
      252016-08-29 17:41:38.7118.35239.68972.626346−10810847−72NF
      262016-08-30 20:35:04.1118.35539.68562.429536−13817071−66NF
      272016-08-31 18:15:27.7118.34739.68972.626446−11511849−66NF
      282016-09-01 01:01:15.9118.34339.69172.23024−558087−93U
      292016-09-01 08:51:20.7118.34939.69672.326427−847763−93NF
      302016-09-05 03:39:19.9118.35739.69152.624263−1049130−64NF
      312016-09-05 04:33:22.3118.36239.69762.526443−11011350−72NF
      322016-09-07 19:58:31.6118.35739.69752.128129−14515974−66NF
      332016-09-10 18:19:15.4118.35739.69663.025986−11316123−9U
      342016-09-10 23:35:59.2118.32939.67962.135255901724690TS
      352016-09-10 18:09:37.3118.35139.70464.326733−13013265−67NF
      362016-09-10 18:13:00.8118.35639.69352.628540−9611350−85NF
      372016-09-11 06:57:02.9118.35739.68572.0338703423558156SS
      382016-09-11 07:20:44.1118.32439.68372.031376−204871−165SS
      392016-09-11 07:41:31.7118.33139.68672.12656615656826SS
      402016-09-11 14:41:39.6118.33139.68672.023568−13614384−6SS
      412016-09-11 16:44:13.2118.34739.70072.026025−25167−153U
      422016-09-11 22:05:05.2118.34639.69472.11347714489167SS
      432016-09-13 09:59:10.3118.34839.68262.924358−1069235−66NF
      442016-09-13 10:50:09.0118.34739.69152.426320−494075−103NF
      452016-09-14 20:00:17.4118.34539.68062.524063−957127−80NF
      462016-09-16 03:34:42.5118.35639.69482.726264−12413942−41NF
      472016-09-17 12:06:37.3118.35139.68072.126628−17317288−62U
      482016-09-17 16:18:33.2118.34939.69672.829277−845424−146NF
      492016-09-17 16:20:04.1118.34639.67972.123549−1118645−67NF
      502016-09-19 07:18:06.4118.34539.70172.923266−1007626−68NF
      512016-09-19 07:31:38.4118.35139.68352.723547−1077946−73NF
      522016-09-20 11:28:38.4118.34439.68362.823869−1039025−60NF
      532016-09-21 21:57:55.0118.34839.70063.226442−11912154−66NF
      542016-09-24 09:09:35.6118.35339.67982.224664−17620486−26SS
      552016-09-29 14:22:36.6118.35339.68653.025954−13714057−45NF
      562016-10-01 12:39:37.2118.35039.68582.129260−14917187−85NF
      572016-10-05 23:27:40.2118.35239.70573.526632−12612765−70NF
      582016-10-07 16:32:04.7118.34839.68852.127643−10211248−79NF
      592016-10-19 04:49:01.1118.34939.67372.22957−14617186−84U

      图  3  2016-08-21 17:46:51.7GPAT震源机制反演结果

      据震源机制解的3个应力轴倾角大小分类,Zoback将震源机制解分为正断型(NF)、正断为主兼走滑型(NS)、走滑型(SS)、逆断型(TF)、逆断为主兼走滑型(TS)、无法确定型(U)共6种[10]。据此标准,对本文所获取的59个震源机制解(表2)进行分类统计分析,不同类型的震源机制解用不同的颜色表示(图4),正断型和正断兼走滑分量类型为41次,占69.5%;走滑型地震为8次,占总数的13.6%;逆断型和逆断兼走滑分量类型为1次,占1.7%;无法确定型为9次,占15.3%。由此可见,本次震群以正断型为主,应力状态主要为剪切拉张。

      图  4  震源机制解分布图

      GPAT震级计算采用的是矩震级MW[11-12],台网观测报告采用近震震级ML[13],二者之间定义不同,必然存在一些差异。图5表明:震群矩震级MW和近震震级ML有从小到大一致性越来越好的趋势,对于ML4.5以上地震是否一致性继续变好,受震群震级范围限制,没有具体数据支持,有待后续研究。另外对于ML3.6以下的地震事件,矩震级一般大于近震震级,且事件震级越小,二者差距越明显,最大差距可达到1级,这和先前研究结果基本一致[7]

      图  5  矩震级MW与地方震级ML对比图

      为了验证震源深度的正确性,把GPAT计算产出的所有事件的G深度与观测报告产出的震源深度(L深度)进行对比(图6a),并把相同事件的2种震源深度的平均值(M深度)假设为震源的正确深度,可以看出G深度和L深度有很好的相关性。为了做出更加精确的定量评价,计算了G深度与M深度的深度偏差(图6b),其中最小偏差为0 km,最大偏差为2.5 km,平均偏差为0.19 km,其中偏差小于2.0 km的占91.5%,整体偏差较小,而0.19 km的偏差可能是由于起始破裂点与“矩心”之间的差别产生的[6],也可能是人工拾取震相误差或者选择台站不同导致的。总而言之,对于小震,GPAT深度与观测报告产出的震源深度具有相当好的一致性,它能够反映震源的真实深度。

      图  6  深度对比图

    • 结合地震发震背景分析,震群发生于唐山断裂带中NE向的唐山-古冶断裂附近,主震为9月10日ML4.3地震。GPAT确定的震源机制节面Ⅰ走向为132°,倾向65°,滑动角−67°;节面Ⅱ走向为267°,倾向33°,滑动角−130°;T轴方位角205°,仰角17°;B轴方位角302°,仰角21°;P轴方位角79°,仰角63°。节面Ⅰ与唐山-古冶断裂带东北段走向基本一致,判断主破裂面为节面Ⅰ,类型为正断型,兼右旋走滑活动,这与CAP和振幅比方法计算的结果总体一致[14]。另外,震群内部事件震源机制结果相似性很高,也为震源机制解的正确性提供了可靠依据。

      以2016年8月河北台网产出的事件波形资料为研究对象,利用GPAT方法计算唐山震群震源机制,得出以下结论:

      1)震群整体震源机制解结果相似性好,结果可信度很高;震源机制解结果以正断型为主,应力状态以剪切拉张为主。

      2)矩震级一般大于近震震级,且事件震级越小,二者差距越明显,最大差距可达到1级;受震群震级范围限制,对于较大地震震级相关性暂时无法判断,有待后续研究。

      3)GPAT深度与台网日常观测报告产出的震源深度具有相当好的一致性,GPAT深度能够反映震源的真实深度。

      总之,GPAT方法可以作为目前台网反演小震震源机制、矩震级和震源深度的一种有效的、实用的手段。

参考文献 (14)

目录

    /

    返回文章
    返回