The Spatial and Temporal Scanning of the b-value in the Western and Mid-western Section of the Zhangjiakou-Bohai Fault Belt
-
摘要: 张家口-渤海断裂带作为华北平原地区重要的活动断裂带,地震活动频繁,是我国地震监测预测重点区域之一。本文选取该断裂带西段及中西段1970-2016年的地震目录,采用最大似然法进行时间扫描,分析显示研究区b值为0.28-1.52,其随时间变化的特点是在大地震发生前降至最低,震后逐渐恢复;研究区空间扫描结果显示,该区b值的平均值为0.93,其中怀安-万全盆地北缘断裂和蓟运河断裂平均b值较低,反映该区域应力水平较高。综合以上结果,本次研究揭示出研究区地震危险性的时间和空间差异,为对研究区地震危险性评价提供基础数据。
-
关键词:
- b值 /
- 张渤断裂带西段及中西段 /
- 地震危险性
Abstract: The Zhangjiakou-Bohai fault belt, as one of the important active faults in North China, as well as earthquake-prone area, becomes one of the crucial region to monitor the earthquake activities in China. In this paper we select the western and mid-western section of the Zhangjiakou-Bohai fault belt as the study region, and apply the maximum likelihood as calculating method to conduct analysis of earthquake catalogue during 1970-2016. The temporal analyses indicate that the seismic b-value ranges from 0.28 to 1.52 and drops to the lowest before the great earthquake and gradually recovers after the earthquake. The spatial analyses illustrate that the b-value shows 0.93 in average, and that in the lower average b-value region, such as north rim fault of the Huai'an-Wanquan basin and Jiyunhe fault shows stronger regional stress level. The spatiotemporal characteristic analyses of the b-value, indicates that the degree of seismic risk varies in time and space in the study area. Our results provide basic data for seismic hazard assessment for the Zhangjiakou-Bohai fault belt. -
引言
对地震动特征规律的分析是地震工程学科研究的主要问题 (刘启方等,2006)。近几十年来国内外发生了多次破坏性地震,对所获得地震动记录进行分析,不难看出位于发震断层附近区域内的地震动具有许多独特性质,例如竖向效应。竖向效应是指近断层区域内竖向地震作用远远超过规范所规定的值的现象,国内外众多地震的强震记录显示竖向地震峰值是水平向的1/2—2/3(Ambraseys等,2003),我国《建筑抗震设计规范》(GB 50011—2010) 规定竖向地震影响系数为水平向的65%(中华人民共和国国家标准,2010),但在许多地震灾害中的近断层区域内出现了竖向地震动高于水平向地震动的现象 (冉志杰等,2012;赵国辉等,2008;周锡元等,2006;Niazi等,1991)。
中国大陆强震台网于2008年3月投入运行,在2008年5月12日汶川特大地震中获得了丰富的近断层强震记录,为相关问题的研究提供了重要的基础数据。通过对强震记录的分析发现近断层地震动竖向与水平分量之比最高达到了1.4(谢俊举等,2010;于海英等,2008)。但是已有的对于竖向效应问题的研究主要集中于对地震动动力特性的分析,关于地震动对工程结构影响的分析较少,尤其缺少对影响效应的定量分析。
本文基于此研究背景,选取汶川地震近断层强震记录为基础数据,以单自由度体系P-Δ效应为研究目标,对近断层竖向地震动所产生的P-Δ效应进行分析研究,对于工程结构的抗震问题具有一定的意义。
1. P-Δ效应的计算方法
已有研究表明,P-Δ效应产生的放大作用是竖向地震导致结构破坏的主要原因之一 (贺秋梅等,2014;刘启方等,2006)。根据结构抗震思想,大多数工程结构在进行抗震设计时都需要转化为等效的单自由度体系的叠加进行分析 (胡聿贤,2006;李宏男,2013;梁炯丰等,2013),所以本文对单自由度体系竖向地震作用下的P-Δ效应进行分析。
1.1 P-Δ效应概念描述
P-Δ效应是指体系在动力荷载作用下,由于其竖向作用使体系结构产生动力附加弯矩的过程,相当于在体系上附加了一个水平地震作用,其原理如图 1所示 (胡聿贤,2006;袁一凡等,2012)。图中m表示质量,P(t) 表示体系在某一时刻承受的水平动力荷载,F(t) 表示体系在某一时刻承受的竖向荷载,包括重力以及竖向动力荷载等。
根据其基本原理,对于地震作用下的单自由度体系,在考虑了竖向地震作用后就变成了具有2个方向自由度的双自由度体系。由于地震所产生的运动以水平运动为主,故本文分析时不考虑竖向地震作用造成的竖向相对运动,因此该体系在某一时刻承受的竖向荷载作用为重力与竖向地震作用引起的惯性力。
强震仪所获得的地震加速度以av(t) 表示,同时竖直方向以加速度向上为正、向下为负,根据其基本原理P-Δ效应等效出的水平地震作用如下式 (1) 所示:
$$ {P_0}(t) = m\left[ {g + {a_v}(t)} \right] \cdot \frac{{u(t)}}{H} $$ (1) 式中,u(t) 为体系在某时刻的位移大小,H为体系高度,P0(t) 为与竖向地震作用等效的水平地震作用。
1.2 P-Δ效应的计算方法
将P-Δ效应所产生的附加水平动力作用带入动平衡方程,则得到考虑P-Δ效应的动力平衡方程,表达如下:
$$ m\ddot u{\rm{(}}t{\rm{)}} + c\dot u{\rm{(}}t{\rm{)}} + ku{\rm{(}}t{\rm{)}} = m\; \cdot \;\left[ {g + {a_v}{\rm{(}}t{\rm{)}}} \right]\; \cdot \;\frac{{u{\rm{(}}t{\rm{)}}}}{H} + m{a_h}{\rm{(}}t{\rm{)}} $$ (2) 式中,c是单自由度体系阻尼系数,k是回复力系数,ah(t) 是水平地震加速度。
式 (2) 可以根据杜哈密积分进行简化求解,进而求得考虑P-Δ效应的单自由度体系的地震反应,从而进一步求得P-Δ效应的放大作用,具体方法为:
(1) 按照地震动力时间步利用杜哈密积分求得每一时间步时刻的水平地震作用动力反应,即得到u(ti) 与ü(ti);
(2) 根据每一时间步的水平位移按照公式 (3) 计算竖向荷载所产生的等效水平地震作用,利用杜哈密积分求解等效水平地震作用所产生的附加地震反应:ü0(ti)
$$ {\ddot u_0}({t_i}) = \left[ {g + {a_v}(t)} \right]\; \cdot \;\frac{{u(t)}}{H} $$ (3) (3) 据所计算出的水平地震反应ü(t) 与附加地震反应ü0(t),按照公式 (4) 计算P-Δ效应所产生的放大作用,式中β为效应放大系数,越大表征所产生放大作用越强。
$$ \beta = \frac{{\ddot u(t) + {{\ddot u}_0}(t)}}{{\ddot u(t)}} $$ (4) 2. 竖向地震动P-Δ效应放大作用的分析
为分析近断层竖向地震动P-Δ效应放大作用的特点,本文分别选取汶川地震中近断层与中远场强震记录,按上述方法分析计算实际地震动所产生P-Δ效应的放大系数,并与中远程强震记录对比分析。
2.1 强震数据的选取
按近断层定义,本文先选取汶川地震中发震断层附近20个强震台站的强震记录,台站的断层距小于60km,具体信息见表 1。再选取位于中远场的20个强震台站的强震记录,台站的断层距大于100km,具体信息见表 2。强震记录来源于中国地震局工程力学研究所下属的中国强震台网中心数据库。
表 1 所选取强震记录的近场台站信息Table 1. The information of near-site stations selected in the study编号 51MZQ 51JYH 51PXZ 51AXT 51JYD 51SFB 51MXN 51WCW 51MXT 51JYC 断层距/km 7.17 13.6 21.05 25.89 26.43 26.73 27.09 27.72 27.77 30.54 编号 51DXY 51DYB 51LXT 51QLY 51LXM 51PWM 51LXS 51GYZ 51BXZ 51XJD 断层距/km 31.19 34.33 46.4 49.17 49.18 51.28 51.58 55.15 57.75 59.09 表 2 所选取强震记录的中远场台站信息Table 2. The information of remote stations selected in the study编号 51HSD 51CXQ 51YAS 51JZB 62SHW 51JZZ 51TQL 51HYQ 51HYJ 51LDD 断层距/km 104.5 104.7 111.1 114.2 121.1 125.4 137.2 155.4 163 177.1 编号 51SMW 51LDL 62TSH 51SMM 51MBD 51YXX 51YXZ 51SMC 51LDS 51MNL 断层距/km 188.7 190.9 191.4 211.5 222.2 249.2 262.2 272.9 278.5 300.4 为充分表达研究目的,给出了汶川地震发震断层与所选取近断层强震台站的空间分布,如图 2所示。
2.2 P-Δ效应放大系数计算结果
对于任意台站,所获得强震记录分为水平向 (EW、NS) 与竖直向 (UD),所以任意台站可以获得2个P-Δ效应放大系数。分别计算所选取近断层与中远场强震台站记录的放大系数,并将计算结果绘制成频率直方图,如图 3所示。计算中,结合大多数工程结构的动力特点,其单自由度体系阻尼比选择为0.05,自振周期分别取1s、2s与3s。
对比近断层与中远场强震台站记录所计算出的放大系数,从整体上可以看出:对于近断层所获得的强震记录其P-Δ效应放大系数值较大,而中远场的P-Δ效应放大系数值较小;对于中远场地震动,其地震作用的水平分量本身就较弱,故其竖向地震动所产生的P-Δ效应可以忽略。因此竖向地震动所产生P-Δ效应的放大作用是近断层地震动所产生动力作用的主要特点,对于近断层地震动应分析其所产生的P-Δ效应。
2.3 放大系数的统计分析
地震动具有很强的不确定性 (Niazi等,1991),为对建筑抗震设计提供具体参考,在完成对所选取强震记录初步分析后,应对多条强震记录的初步分析结果进行统计规律分析。
基于所得到的20个近断层台站获得的强震记录,计算得到不同自振周期条件下P-Δ效应的放大系数。利用概率图工具分析放大系数服从的概率分布,如图 4所示 (自振周期为1s)。根据分析,放大系数服从正态分布。进一步计算不同自振周期单自由度体系放大系数的统计参数,计算结果如表 3所示。
表 3 放大系数的统计参数Table 3. Statistical parameter of amplification factor自振周期/s 1 2 3 均值 1.012 1.229 1.701 方差 0.002 0.014 0.063 从表 3可以看出,对于同一自振周期的单自由度体系,其放大系数的离散性较小,放大系数主要分布于均值附近,且主要受单自由度体系自振周期的影响,自振周期越大,其P-Δ效应放大系数越高。为此,可以认为,在近断层区域内,竖向地震动所产生P-Δ效应的放大作用具有普遍性,其放大作用主要受自振周期大小的影响。
3. P-Δ效应放大系数谱的建立
3.1 放大系数谱的建立思路
在近断层范围内竖向地震引起的P-Δ效应会产生普遍放大作用,其放大系数主要受自振周期的影响,因此可以参考地震动反应谱的基本思想建立P-Δ效应放大系数谱,以获得竖向地震动P-Δ效应的放大作用与体系自振周期的关系。以上分析可知,当体系自振周期一定时P-Δ效应放大系数的离散性较小,所以可以用自振周期为一定时不同台站强震记录放大系数的均值作为P-Δ效应放大系数的代表值,从而建立放大系数谱,建立步骤如下:
(1) 设定不同的自振周期,分别为Tk=0.02×k(k=1,2,…,200)。
(2) 对于所设定的不同的自振周期,分别计算每个台站地震记录的P-Δ效应放大系数,即得到βi(Tk),其表示第i个台站的强震记录在自振周期为Tk时的P-Δ效应放大系数。
(3) 当体系自振周期为Tk时,求得不同台站强震记录放大系数的均值,即u[β(Tk)],进而求得不同自振周期条件下,其放大系数的均值u[β(Tk)]。
(4) 绘制放大系数曲线,横坐标为Tk,纵坐标为u[β(Tk)],基于所得曲线利用最小二乘法拟合,得到规准化后的放大系数谱。
3.2 放大系数谱的建立
建立近断层竖向地震P-Δ效应的放大系数谱,需先求得Tk-u[β(Tk)]曲线,然后进行最小二乘拟合。为了方便工程应用,拟合函数选用线性函数进行。从Tk-u[β(Tk)]曲线中可以看出,在自振周期Tk < 2s时放大系数几乎都小于1.1,曲线趋近于一条斜率为零的直线,在自振周期Tk > 2s时放大系数β > 1.1,P-Δ效应具有明显的放大作用,β随Tk的变化趋近于单调上升的线性函数曲线,所以对放大系数规准谱曲线以Tk=2s为分界点进行分段拟合。利用最小二乘法进行分段拟合并经过简化处理得到规准化后的放大系数谱,如图 5所示。
进一步拟合放大系数谱的数学关系,其关系表达式为式 (5):
$$ \beta = \left\{ \begin{array}{l} {\rm{1}}.{\rm{05}},\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; {\rm{0s}}< {T_k} \le {\rm{1}}.7{\rm{s}}\\ {\rm{0}}.{\rm{6}} \times {T_k} + {\rm{0}}.{\rm{03}},\;\;\;\;{\rm{1}}.{\rm{7s}}\; < {T_k} \le {\rm{4s}} \end{array} \right. $$ (5) 此公式可为近断层区域内抗震设计问题中考虑竖向地震作用P-Δ效应的参考公式,在结构抗震设计中可以与反应谱结合使用。
4. 结语
本文以汶川大地震近断层的强震数据为基础,对近断层竖向地震动所产生的P-Δ效应进行分析,结果如下:
(1) 汶川地震中相比于中远场区域,近断层区域内竖向地震作用所产生的P-Δ效应具有明显的放大作用,并且具有普遍性。P-Δ效应放大系数主要受体系自振周期的影响,当自振周期一定时,不同地震动的放大系数服从正态分布并且离散性较小,主要分布于均值附近。
(2) 建立了放大系数随体系自振周期变化的放大系数谱,从而为结构抗震设计提供了参考依据;认为在可能发生强烈地震的活断层附近区域内,结构抗震设计过程中应该考虑竖向地震动所产生P-Δ效应引起的放大作用,放大系数按照规准后的放大系数谱曲线确定,在结构抗震设计中,放大系数谱可以与反应谱结合使用。
(3) 由于中国大陆地区强震记录数量的限制,本文仅仅选择了汶川地震的近断层强震记录;随着近断层强震数据的增多,应根据发震断层与地质构造特点,建立适用于不同区域、不同地质场地条件的放大系数谱,从而使其在工程应用中日益完善。
-
表 1 研究区地震数据统计结果(据国家地震前兆台网中心测定)
Table 1. The statistical seismic data of the study area
震级/ML 1.0-1.9 2.0-2.9 3.0-3.9 4.0-4.9 5.0-5.9 6.0-6.9 7.0-7.9 8.0-8.9 频数/个 4315 2688 378 188 23 3 1 0 震源深度/km 0-0.9 1-5 6-10 11-15 16-20 21-25 26-30 31-61 频数/个 2200 593 2631 1132 603 289 130 18 表 2 研究区ML≥5.0的地震序列(据国家地震前兆台网中心测定)
Table 2. The earthquake sequence data (ML≥5.0)of the study area
日期(年-月-日) 时间 北纬/° 东经/° 震源深度/km ML 参考地点 2006-07-04 11:56 38.89 116.28 20 5.5 河北文安县 1999-03-11 21:18 41.16 114.39 5.5 河北张北县 1998-01-10 11:50 41.10 114.31 10 6.2 河北尚义县 1990-07-21 08:41 40.58 115.83 8 5.0 北京延庆县 1988-07-23 13:51 40.08 114.22 5.0 河北阳原县 1983-04-03 10:16 40.75 114.78 20 5.1 河北万全县 1980-02-07 13:31 39.52 117.90 5.2 天津宁河县 表 3 研究区主要活动断裂和平均b值
Table 3. The main active faults and the average b-value of the study area
编号 断裂名称 平均b值 长度
/km产状 断裂性质 平均滑动速率/mm·a-1 最新活动时代 走向 倾向 倾角 F1 张家口断裂 0.80 70 NNW SE 60° 正断层 >0.1 Q3 F2 新保安-沙城断裂 0.76 32 NWW SW 50°—70° 正断层 0.2 Q4 F3 南口-孙河断裂 0.78 58 NWW SW,NE 70° 正断层 0.3 Q4 F4 夏垫断裂 0.78 23 NE SE 50°—70° 正走滑 Q4 F5 永定河断裂 0.78 26 NW SW,NE 70°—75° 正走滑 Q2 F6 廊坊-武清断裂 1.02 50 NW SW 正断层 Q3 F7 蓟运河断裂 0.68 50 NW SW 70° 正断层 Q4 -
陈凌, 刘杰, 陈颙等, 1998.地震活动性分析中余震的删除.地球物理学报, 41(S1):244-252. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx1998s1025&dbname=CJFD&dbcode=CJFQ 陈培善, 白彤霞, 李保昆, 2003.b值和地震复发周期.地球物理学报, 46(4):510-519. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx200304012&dbname=CJFD&dbcode=CJFQ 方颖, 江在森, 张晶等, 2008.张家口-渤海断裂带现今运动状态分析.大地测量与地球动力学, 28(1):11-15. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dkxb200801005&dbname=CJFD&dbcode=CJFQ 高战武, 徐杰, 宋长青等, 2001.张家口-蓬莱断裂带的分段特征.华北地震科学, 19(1):35-42, 54. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=hdkd200101005&dbname=CJFD&dbcode=CJFQ 韩孔艳, 2009.张家口-渤海构造带的分段性与地震活动特征研究.北京:中国地震局地质研究所. 何付兵, 白凌燕, 王继明等, 2013.夏垫断裂带深部构造特征与第四纪活动性讨论.地震地质, 35(3):490-505. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzdz201303004&dbname=CJFD&dbcode=CJFQ 侯贵廷, 叶良新, 杜庆娥, 1999.渤张断裂带的构造机制及其地质意义.地质科学, 34(3):375-380. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzkx199903012&dbname=CJFD&dbcode=CJFQ 赖晓玲, 李松林, 孙译, 2007.渤海及邻区3次7级以上地震的深部构造背景.大地测量与地球动力学, 27(1):31-33, 54. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dkxb200701006&dbname=CJFD&dbcode=CJFQ 李纪汉, 1987.b值影响因素的初步研究.地震学刊, 7(2):49-53. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzxk198702005&dbname=CJFD&dbcode=CJFQ 李全林, 陈锦标, 于渌等, 1978.b值时空扫描——监视破坏性地震孕育过程的一种手段.地球物理学报, 21(2):101-125. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx197802001&dbname=CJFD&dbcode=CJFQ 李涛, 吴小平, 王莹, 2010.基于分区空间扫描研究云南地区b值分布特征.云南大学学报(自然科学版), 32(6):670-675. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=yndz201006013&dbname=CJFD&dbcode=CJFQ 刘静伟, 吕悦军, 谢富仁, 2016.利用b值空间分布解析2014年11月22日康定MS 6.3地震的发生.地球物理学进展, 31(2):553-558. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwj201602007&dbname=CJFD&dbcode=CJFQ 刘艳辉, 赵根模, 吴中海等, 2015.青藏高原东南缘及邻区近年来地震b值特征.地质通报, 34(1):58-70. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=zqyd201501005&dbname=CJFD&dbcode=CJFQ 刘志宏, 梁一鸿, 徐仲元等, 2011.构造地质学.2版.北京:地质出版社. 任雪梅, 高孟潭, 冯静, 2011.地震目录的完整性对b值计算的影响.震灾防御技术, 6(3):257-268. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20110305&journal_id=zzfyjs 沈繁銮, 1999.地震活动性参数b值和年发生率的一种取值方法.华南地震, 19(4):41-45. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=hndi199904006&dbname=CJFD&dbcode=CJFQ 沈小七, 姚大全, 陈安国等, 2010.江淮地区地震精定位及b值随深度的变化研究.地震地磁观测与研究, 31(6):1-7. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzgj201006002&dbname=CJFD&dbcode=CJFQ 盛艳蕊, 张子广, 周晓成等, 2015.新保安-沙城断裂带土壤气地球化学特征分析.地震, 35(4):90-98. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dizn201504010&dbname=CJFD&dbcode=CJFQ 孙文福, 顾浩鼎, 1992.怎样正确计算b值.东北地震研究, 8(4):13-27. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ddyj199204002&dbname=CJFD&dbcode=CJFQ 索艳慧, 李三忠, 刘鑫等, 2013.中国东部NWW向活动断裂带构造特征:以张家口-蓬莱断裂带为例.岩石学报, 29(3):953-966. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=ysxb201303018&dbname=CJFD&dbcode=CJFQ 王辉, 曹建玲, 荆燕等, 2012.川滇地区强震活动前b值的时空分布特征.地震地质, 34(3):531-543. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzdz201203016&dbname=CJFD&dbcode=CJFQ 王帅军, 张先康, 刘宝峰等, 2005.张渤地震构造带中西段及邻区深部构造探测.大地测量与地球动力学, 25(3):110-115. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dkxb200503021&dbname=CJFD&dbcode=CJFQ 汪素云, 俞言祥, 2009.震级转换关系及其对地震活动性参数的影响研究.震灾防御技术, 4(2):141-149. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20090202&journal_id=zzfyjs 王熠熙, 张辉, 刘双庆等, 2015.河北平原地震带b值时空变化特征.地震工程学报, 37(1):188-195. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=zbdz201501031&dbname=CJFD&dbcode=CJFQ 武敏捷, 朱红彬, 岳晓媛等, 2013.基于多地震活动参数的首都圈地区地震危险性分析.华北地震科学, 31(1):25-30. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=hdkd201301007&dbname=CJFD&dbcode=CJFQ 谢卓娟, 吕悦军, 彭艳菊等, 2012.东北地震区小震资料完整性分析及其对地震活动性参数的影响研究.中国地震, 28(3):256-265. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=zgzd201203004&dbname=CJFD&dbcode=CJFQ 谢卓娟, 李山有, 吕悦军, 2015.滇西南地区主要活动断裂的b值空间分布特征.地球科学——中国地质大学学报, 40(10):1755-1766. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqkx201510016&dbname=CJFD&dbcode=CJFQ 张建中, 宋良玉, 1981.地震b值的估计方法及其标准误差——应用蒙特卡罗方法估计b值精度.地震学报, 3(3):292-301. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzxb198103007&dbname=CJFD&dbcode=CJFQ 张磊, 白凌燕, 蔡向民等, 2014.北京平原南口-孙河断裂带北西段活动性分析.中国地质, 41(3):902-911. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dizi201403017&dbname=CJFD&dbcode=CJFQ 张琳琳, 唐兰兰, 2015.新疆天山地区b值时空变化特征.地震, 35(3):93-101. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dizn201503010&dbname=CJFD&dbcode=CJFQ 周月玲, 尤惠川, 2010.张家口断裂第四纪构造变形与活动性研究.震灾防御技术, 5(2):157-166. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20100202&journal_id=zzfyjs El-Isa Z. H., Eaton D. W., 2014. Spatiotemporal variations in the b-value of earthquake magnitude-frequency distributions:classification and causes. Tectonophysics, 615-616:1-11. doi: 10.1016/j.tecto.2013.12.001 Gutenberg B., Richter C. F., 1954. Seismicity of the earth and associated phenomenon. 2nd ed. Princeton:Princeton University Press. Ishimoto M., Iida K., 1939. Observations of earthquakes registered with the microseismograph constructed recently. Bulletin of the Earthquake Research Institute, 17:443-478. Keilis-Borok V. I., Knopoff L., Rotvain I. M., 1980. Bursts of aftershocks, long-term precursors of strong earthquakes. Nature, 283(5744):259-263. doi: 10.1038/283259a0 Khan P. K., Chakraborty P. P., 2007. The seismic b-value and its correlation with Bouguer gravity anomaly over the Shillong Plateau area:tectonic implications. Journal of Asian Earth Sciences, 29(1):136-147. doi: 10.1016/j.jseaes.2006.02.007 Khan P. K., Ghosh M., Chakraborty P. P., et al., 2011. Seismic b-value and the assessment of ambient stress in Northeast India. Pure and Applied Geophysics, 168(10):1693-1706. doi: 10.1007/s00024-010-0194-x Maden N., Öztürk S., 2015. Seismic b-values, Bouguer Gravity and heat flow data beneath eastern Anatolia, Turkey:tectonic implications. Surveys in Geophysics, 36(4):549-570. doi: 10.1007/s10712-015-9327-1 Shi Y. L., Bolt B. A., 1982. The standard error of the Magnitude-frequency b value. Bulletin of the Seismological Society of America, 72(5):1677-1687. https://pubs.geoscienceworld.org/ssa/bssa/article-abstract/72/5/1677/102113/the-standard-error-of-the-magnitude-frequency-b -