Study on Characteristics of Seismic Sequence Types and Early Aftershock Prediction Method in Hebei and Its Surrounding Regions
-
摘要: 本文梳理了河北及邻区（京津冀地区）1970年至2022年MS≥3.5地震目录，对MS≥3.5主震进行地震类型判定。研究发现，在该地震目录中，主余型地震66次，多震型地震24次，孤立型地震58次。本文对研究区内地震类型及其空间分布特征进行分析，并利用震级均值特性计算地震序列的最大余震震级，结果表明：（1）研究区域内的地震类型以主余型为主（45%），孤立型次之（39%），多震型最少（16%）。9次MS≥5.0主震中，主余型3次，孤立型6次。随着震级增大，主余型地震所占比例也随之增加。（2）主余型地震序列和多震型地震序列中，最大余震在震后1 d内发生的占比最大。（3）地震类型的空间分布显示，3种类型的地震，沿着张家口—渤海地震带和河北平原地震带分布。孤立型地震在这2个条带上分布比较均匀；主余型地震主要分布在张家口、唐山、邢台一带；多震型地震与主余型地震分布情况基本相同。（4）利用震级均值特性，计算该区MS≥5.0主震（主余型）的最大余震震级，发现3个研究震例计算得到的最大余震震级与实际记录比较相近，可以用此方法预估河北地区未来中强震的最大余震，为震后强余震的判断提供参考方法。Abstract: This study reorganizes the earthquake catalog of MS≥3.5 events from 1970 to 2022 in Hebei and its surrounding regions (Beijing-Tianjin-Hebei) and classifies the earthquake types of MS≥3.5 main earthquakes. The analysis identifies 66 mainshock-aftershock earthquakes, 24 multiple-mainshock earthquakes, and 58 isolated earthquakes. The study examines the types and spatial distribution of earthquakes and estimates the maximum aftershock magnitude of earthquake sequences using the average magnitudes of aftershock sequences. The key findings are as follows: (1) Mainshock-aftershock earthquakes are the most prevalent (45%), followed by isolated earthquakes (39%) and multiple-mainshock earthquakes (16%). Among the MS≥5.0 earthquakes (nine in total), three are mainshock-aftershock events, while six are isolated earthquakes. Additionally, the proportion of mainshock-aftershock earthquakes increases with earthquake magnitude. (2) In both mainshock-aftershock and multiple-mainshock sequences, the largest aftershocks predominantly occur within one day after the main event. (3) The spatial distribution of earthquake sequences reveals that all three earthquake types are concentrated along the Zhangjiakou-Bohai and Hebei Plain seismic belts. Isolated earthquakes are evenly distributed across these belts, while mainshock-aftershock and multiple-mainshock earthquakes are mainly concentrated in Zhangjiakou, Tangshan, and Xingtai. (4) The maximum aftershock magnitude for MS≥5.0 mainshock-aftershock earthquakes is estimated using the average magnitudes of aftershock sequences. The calculated values closely match observed data, indicating that this method can serve as a predictive tool for assessing the maximum aftershock magnitude of moderate to strong earthquakes in Hebei. This approach provides a valuable reference for evaluating the likelihood of strong aftershocks following seismic events.
-
引言
微震检测在地震学诸多研究中均十分重要，如大地震前震活动性（Kato等，2012；Brodsky等，2014）、地震触发（Hill等，1993；Kilb等，2000；Gomberg等，2001；Peng等，2010；Meng等，2013, 2014；Wang等，2015）、大地震成核（Kato等，2012；Brodsky等，2014；Ellsworth等，2018）、断层形态刻画（Mori等，1990；Dreger等，2000；Kuge，2003；Kao等，2007；Yang等，2009）、水力压裂诱发地震检测（House，1987；Rutledge等，2003；Eisner等，2008；Anikiev等，2014）、小当量核试验检测（Waldhauser等，2004；Wen等，2010；Zhang等，2015a）等，同时，微震检测研究具有一定挑战性。国内外学者提出了多种地震事件识别方法，如应用较多的长短时窗比检测法（Stevenson，1976；Earle等，1994；Kao等，2004, 2007；Liao等，2012；Gharti等，2010；Grigoli等，2013a, 2013b；Drew等，2013；Langet等，2014）、利用互相关技术的波形模板匹配检测方法、震源扫描叠加法和波形自相关检测技术。由于微震记录信噪比较低，长短时窗比检测法虽较简单、计算速度快，但无法有效应用于微震检测中。震源扫描叠加法检测能力较低，限制了对微震震级的检测，对震源位置进行搜索时需利用绝对走时，该过程十分依赖速度模型。因此，本文对应用较广的波形模板匹配检测方法、匹配定位技术、基于图形处理器（Graphics Processing Unit，简称GPU）加速的匹配定位技术（Liu等，2020）与波形自相关检测技术进行研究，应用基于GPU加速的匹配定位技术与双差地震定位法（Waldhauser等，2000）对北京地区19个台站记录的2015年连续地震资料进行分析，获得完备的地震目录，并进行精定位，从而为北京地区地震活动性与发震断层深部构造研究提供基础资料。
1. 波形模板匹配检测方法
存在噪声的情况下，将模板事件波形与输入的待识别连续数据进行互相关是检测已知信号的主要手段，可应用于多个台站的情况（Gibbons等，2006），以此为基础的波形模板匹配检测方法基本检测步骤为：（1）对模板事件波形与连续波形进行相同参数滤波；（2）选定滑动时间窗口与移动步长，进行模板事件波形与连续波形之间的互相关，搜索与模板事件相似的事件，并得到相应的互相关系数（Shelly等，2007；Peng等，2009；Meng等，2013）；接着对多个台站分量记录到的模板事件波形与记录的连续波形进行互相关叠加，并计算平均互相关系数，以此衡量检测到的事件与模板的相似度；（3）将检测到的高于一定阈值的事件作为结果，如果检测结果中出现重复事件，则选取窗口长度内互相关系数最高的事件。判断检测到的事件是否作为地震事件，取决于检测阈值，因此阈值的设定尤为重要。通常对每个互相关系数进行中位数绝对偏差计算（Shelly等，2007）。检测结果中的震中位置代替了相应模板检测事件的震中位置，同时，根据模板事件和检测事件之间的震级比与振幅对数比线性关系计算检测事件的震级（Peng等，2009；Wang等，2015）。在数据处理过程中，Meng等（2013）设置了数据自检环节，为使所有模板事件均实现完美自检（得到的互相关系数＞0.98），在叠加互相关系数时允许1个数据点存在偏移，即在每个时间点，其自身及2个相邻点间的最高互相关系数均用于叠加。进行偏移处理后，可叠加所有分量最佳相关窗口。
通过波形模板匹配检测方法可得到更完整的地震目录，降低地震目录的完备震级。Shelly等（2007）通过波形模板匹配检测方法，将先前记录和定位的低频地震作为模板事件，检测出了日本四国岛震颤活动，同时认为此次震颤活动为低频地震。该方法可在时间和空间上以高分辨率寻找震颤来源，表明震颤和缓慢异化是单个过程的不同表现。Peng等（2009）应用该方法，以发生在帕克菲尔德断层附近的3647个重定位地震事件作为模板，检测2004年9月28日帕克菲尔德Mw6.0地震发生后3d内漏掉的余震事件，捕捉到了比北加利福尼亚地震台网标准目录给出的余震事件多10倍的地震事件，并在此基础上研究了该次地震发生后早期余震的迁移特性。Meng等（2013）应用波形模板匹配检测方法检测了2003年圣西蒙Mw6.5地震发生前46d到发生后31d圣安德烈斯断层帕克菲尔德段的地震，得到的地震事件数量为北加利福尼亚地震台网地震目录数量的8倍，且发现沿圣安德烈斯断裂蠕动部分的地震事件数量在圣西蒙主震后呈显著下降趋势，这与主震产生的负静态应力变化（即应力阴影）密切相关。Kato等（2014）应用波形模板匹配检测方法对2014年智利伊基克Mw8.1地震震源附近记录到的连续波形数据进行了检测，得到的地震事件数量为美国地质勘探局地震台网地震目录数量的10倍，并确定了多个地震迁移序列，还识别到了重复地震，推测为前震序列中沿板块边界断层的地震滑动情况。Walter等（2015）应用波形模板匹配检测方法分析了哥斯达黎加尼科亚半岛2012年9月5日7.6级地震发生前的地震活动，大大提高了地震目录的完备性，确定了主震前约35min、距震源15km以内的小地震（震级 < 2.2），并分析了地震事件触发的可能性。Huang等（2017）通过波形模板匹配检测方法恢复2015年Mw8.0戈尔卡地震发生前80d到发生后30d的地震，检测到的余震数量是国际地震中心地震目录所列余震数量的2倍，发现2015年5月12日7.3级最大余震发生前，发生了1次沿主震破裂断层走向的余震扩展。
波形模板匹配检测方法可有效检测信噪比较低的地震信号，但进行模板事件波形与连续波形互相关时耗时较长。目前，GPU计算已从固定功能的图形设备演变为高度可编程的并行处理器。波形模板匹配检测方法在多个级别上表现出并行性：（1）存在多个模板，每个模板均需与连续波形进行匹配；（2）连续波形本身由一连串的数据点窗口组成；（3）模板和连续波形均包含多个台站记录的多个分量数据；（4）在并行结构的更底层，每个分量均包含多个数据点，这些数据点为互相关分析输入文件，可实现并行。
为降低时间成本，Meng等（2012）提出了基于GPU显著提速的检测算法。为探索这种层次结构的并行性，研究人员将计算分解为多个任务，利用分解后的任务计算模板和滑动窗口内连续波形的相关性，其中每个任务均将1个模板通过滑动窗口与相同长度的连续波形进行匹配。该检测方法的设计集中在以下算法方面：（1）将总任务分组并分配给GPU进行计算；（2）将每个任务的GPU代码优化。与串行的CPU代码相比，Meng等（2012）通过GPU处理将1张英伟达（Nvidia）显卡的速度提高了约40倍。
对于任务分配，需注意以下方面：（1）每组任务不应超过GPU的内存容量，使GPU有足够的内存容纳数据（包括输入、输出和临时处理数据）；（2）分配时应涵盖所有任务，并在各小组间保持均衡，以确保将所有任务分配给GPU，各任务GPU计算时间相对一致；（3）分配任务时应将传送和冗余的计算成本最小化。
Meng等（2012）通过使用GPU提速代码，仅用约48h完成了以下计算：使用2000个模板波形、80d的连续波形，共检测到约24000个地震事件，约为南加州地震台网地震目录数量的70倍，并根据新检测到的事件和南加州地震台网地震目录分析计算了地震活动率。Meng等（2014）利用该方法，在2010年Mw7.2El Mayor-Cucapah主震附近的索尔顿海地热田和圣哈辛托断裂带得到了更完整的地震目录，新检测到的目录揭示了2个研究区域在主震后的地震活动率均显著提高。Liu等（2019）利用波形模板匹配检测方法和双差地震定位法研究门源主震后余震的时空演变，将中国地震台网中心提供的地震目录和重定位后的地震目录中的事件作为模板，扫描了2016年连续波形数据，得到的事件较中国地震台网中心提供的地震目录和重定位后的地震目录数量分别多3倍和11倍。利用更完整的地震目录，对研究区域余震时空演变进行分析。
波形模板匹配检测方法使用的互相关技术对于微弱信号极其敏感，检测到的地震震级远小于其他方法，但仅能检测模板事件附近的地震，并假设待检测地震与模板事件发生在同一位置，不能提供检测事件的位置信息（张淼，2015）。
2. 匹配定位技术
针对波形模板匹配检测方法的局限性，Zhang等（2015a）研发了微震检测新方法，即匹配定位技术，该技术与波形模板匹配检测方法最大不同之处在于，在相关波形叠加前考虑了微震与模板之间的位置差，从而在检测的同时对检测事件进行了定位。
匹配定位技术首先根据已有台网地震目录确定扫描所需的模板事件，并确定模板事件的位置信息；然后在每个台站的所有分量上选取模板事件振幅最大的震相作为参考震相，并将其与连续波形进行滑动互相关；接着以模板事件所在位置为搜索中心，在经度、纬度和深度方向上对可能存在微震的位置进行网格化搜索，搜索过程中首先计算微震每个可能存在的格点位置与参考模板位置到同一台站的到时差，通过到时差对每个台站各分量滑动互相关值进行位置校正；最后计算各台站所有分量叠加后的滑动互相关系数平均值与信噪比，当上述2个参数大于设定的阈值时，即认为捕捉到了1个地震事件，并将该网格点的空间位置确定为相应的震源位置。捕捉到的地震事件震级依靠于所有台站分量中参考震相的振幅比中位数(张淼，2015；Peng等，2009；Meng等，2013)。
匹配定位技术将滑动互相关系数平均值与信噪比作为阈值控制输出结果，通过降低检测标准增加检测数量，降低误检概率。另外，为避免重复检测，在特定时间窗口中选取相关系数最大的地震作为检测到的地震事件。
Zhang等（2015a）应用匹配定位技术检测并定位了朝鲜于2010年5月12日进行的小当量核试验，将朝鲜2009年和2013年核试验产生的地震作为模板事件，对2010年4月1日至2010年5月31日距朝鲜核试验基地200km内的7个地震台站记录到的连续波形进行扫描，得到2010年核试验位置，进一步确认该事件的频谱比Pg/Lg，最终认定其为爆炸类型，为朝鲜2010年5月12日小当量核试验的存在提供了地震学证据，展示了利用地震数据对小当量核试验进行科学检测的能力。Zhang等（2015b）利用匹配定位技术对2007年和2014年日本Ontake火山2次活动过程中产生的连续波形进行扫描，检测到2次活动中的地震数量分别为日本气象厅提供的地震目录数量的37倍和30倍。根据新的完备地震目录，发现2014年日本Ontake火山活动中存在地震数量增多和减少交替出现的规律，该规律出现在大规模地震活动初期（张淼，2015）。
应用匹配定位技术时，首先进行相对走时校正，然后对互相关系数进行最大化叠加，随后在模板周围三维范围内搜索新识别事件的最佳位置，最后提供检测事件的位置信息。该方法虽已使用共享内存并行编程（OpenMP）进行加速计算，但由于计算机CPU内核有限，随着地震数据量的增加和搜索网格的减小，处理数据产生的时间成本仍不可忽略。计算时间主要由匹配定位技术中的波形互相关及校正时间偏移后的互相关过程决定，将任务进行合理分配是亟待解决的问题。Liu等（2020）对匹配定位技术进行改进，得到基于GPU加速的匹配定位技术，该技术为模板事件波形添加了权重因子，并进一步完善了GPU计算问题，首先将模板事件波形和连续波形加载到CPU随机存取存储器；然后计算所有模板事件位置与潜在事件位置的偏移量；接着将模板事件波形和连续波形从CPU随机存取存储器传输到GPU随机存取存储器，再传输到线程块共享内存，每个模板有2个GPU内核循环，以计算滑动相关系数，并根据偏移量将滑动相关系数叠加，将叠加后的滑动相关系数传回CPU随机存取存储器；最后合格的检测结果形成了新目录。基于GPU加速的匹配定位技术增强了检测能力，扩大了检测范围，显著提高了计算效率，且降低了可检测地震的震级，可对检测到的地震事件进行定位。Liu等（2020）利用基于GPU加速的匹配定位技术对2015年阿拉善左旗地震余震进行检测，以中国地震台网中心提供的38次余震作为模板事件波形，扫描主震后7天的连续波形，检测到的余震数量为中国地震台网地震目录数量的20倍左右，并依据余震空间分布与震源机制解，判断此次5.8级地震发生在东西向隐伏走滑断层上。
3. 波形自相关检测技术
Brown等（2008）在解决地颤动中低频地震搜索问题时，提出了波形自相关检测技术，该技术无需先验的地震信号信息，可认为是彻底的盲搜索技术。首先对连续波形进行适当的预处理，如滤波等；然后将各台站各分量连续波形数据按一定时窗长度间隔分成若干部分，由于无先验的模板信息，相对于Shelly等（2007）使用的时间窗口，波形自相关检测技术时窗长度更长；接着将所有可能成对的时窗分别进行互相关并叠加，选取超过一定阈值的事件作为候选事件；最后将所有候选事件对应在同一台站的波形进行额外互相关运算，并按互相关系数高低进行分组并叠加形成高信噪比模板，运用波形模板匹配检测方法进行补充搜索（刘翰林等，2017）。
Brown等（2008）运用波形自相关检测技术对6个台站数据进行处理，选取6s的时间窗进行计算，将大于5倍中位数绝对偏差的事件确定为候选事件，并利用这些事件进行模板匹配滤波扫描，得到287个低频地震事件，通过对低频地震事件进行分析，认为日本西南部的深部震颤是重复的低频地震群，发生在30—35km、地震发生带板块界面下倾面处。
波形自相关检测技术可在模板准备阶段提供帮助，得到更多高信噪比模板事件，但该技术依赖于大量相关计算，对于数据量庞大的研究而言，仍存在计算量大和时间成本过高等问题。
4. 基于GPU加速的匹配定位技术在2015年北京地区微震检测中的应用
4.1 台站分布与研究数据
本文收集了中国地震台网中心与北京市地震局在北京地区布设的19个固定地震台站记录到的2015年连续波形，研究区域位置和台站分布如图 1所示。在台网中心提供的地震目录中筛选245个地震事件作为模板事件，利用基于GPU加速的匹配定位技术对连续波形进行扫描检测。
图 1 研究区域台站分布图Figure 1. Distribution of stations in the study area注：红色虚线表示研究区域断裂带，蓝色三角表示台站。4.2 检测过程
（1）模板事件和连续波形均进行相同参数的减采样率处理，将采样率由100Hz降至20Hz，以提升计算速度。
（2）对模板事件和连续波形进行2—8Hz带通滤波，根据区域一维速度模型，使用TauP（Crotwell等，1999)计算模板事件P波和S波理论到时，并对到时进行人工核验与校正。
（3）模板扫描窗口为S波到时前1s至后5s，对模板事件每个分量波形数据质量赋予权重因子。考虑研究区域范围、模板分布与时间成本，以模板事件位置为中心，设置搜索区域为8km×8km，搜索步长为2km（图 2）。
图 2 研究区域模板事件分布与扫描范围图Figure 2. The distribution of templates and scanning range for each event注：蓝色三角表示台站，黑色虚线表示断层，红点表示2015年台网地震目录中筛选出的模板事件，灰色实线框表示对应模板事件的扫描范围。（4）利用双差地震定位法对地震事件进行精定位（杨智娴等，2003；朱艾斓等，2005；张广伟等，2011；胡幸平等，2013），共获得922个精确的地震位置。
4.3 结果与分析
本文利用245个模板事件，经扫描检测得到高达1229个地震事件，约为地震台网地震目录数量的5倍。模板事件为2015年1月15日2.0级地震事件，检测到的事件为2015年1月16日1.01级地震事件（图 3）。精定位后的地震分布如图 4所示，由图 4可知，地震主要集中在黄庄-高丽营断裂（图 4中白线）附近。地震数量-震级关系如图 5所示，由图 5可知，检测到的地震事件以1级以下微震为主。
图 4 模板事件（黄色）与检测到的地震事件（红色）精定位后震中分布图Figure 4. The distribution of the catalog events (yellow) and detected events (red) after relocation检测到地震在潭柘寺地区密集分布，震源深度主要为5—10km，分布趋势为北西向，如图 6所示，这补充了台网中心地震目录在该范围的空缺。由于该区域附近的模板事件数量较少，导致检测到的结果数量少，可见南东侧地震分布较深，北西侧地震分布较浅，提示断裂倾向为南东。在黄庄-高丽营断裂带北端，震源深度分布较广泛，以10—20km为主，震源深度较潭柘寺地区深，总体分布情况与台网地震目录相似。在垂直于断裂带走向的剖面（图 6（e）），地震集中分布情况显示出断裂倾向为南东，且断层倾角较大。检测到的地震目录完善了地震台网原始地震目录，更好地揭示了该地区断层形态。
图 6 扫描重定位后的台网地震目录事件（黄色）与检测到的地震事件（红色）震源深度沿不同剖面方向的垂直分布图Figure 6. The distributions of catalog events (yellow) and detected events (red) focal depth along the AB, AC, IK, and IJ sections after relocation5. 总结与讨论
完备地震目录在地震学研究中十分关键且重要，近年来，微震检测与识别是地震学研究的热门课题。本文主要介绍国内外有效且应用广泛的微震识别方法，包括波形模板匹配检测方法、匹配定位技术和波形自相关检测技术，其中波形模板匹配检测方法可在信噪比较低的连续波形数据中识别出与模板相似度较高的地震事件，但不能同时进行地震定位；匹配定位技术弥补了波形模板匹配检测方法的不足，可在模板周围进行全方位检测，将识别到与模板事件相似度高的事件作为结果，且可对检测到的地震进行定位，但时间成本随着搜索网格的缩小而增加；基于GPU加速的匹配定位技术克服了匹配定位技术的缺点，通过添加权重因子，可更有效地进行地震检测和定位，更好地丰富地震目录；波形自相关检测技术虽可在模板准备阶段提供一定帮助，但由于计算量过大，时间成本过高，不便应用于数据量较大的研究。
通过对微震识别方法进行总结与对比，可知匹配定位技术适用于各类微弱信号的检测，在地震迁移过程中可提供精准的位置信息，也可较好地刻画断层走向，对于地震活动性的研究具有十分重要的意义。在今后的研究工作中，可利用该方法获取更完备的地震目录，得到更多的地下断层信息，为组织开展地震区划与各类基础工程建设提供充实的地震资料，方便进行相应的风险评估工作，提高地震防御能力。引入GPU并行计算可提升计算速度，有效提高微震识别效率与精度，不同学科的交叉与结合将成为地震学研究趋势。
利用基于GPU加速的匹配定位技术对北京地区2015年地震资料进行遗漏地震检测和识别，共识别1229个地震事件，约为地震台网地震目录数量的5倍，获得了较精确的地震位置、震级等信息。对检测到的地震事件进行精定位，可知地震在潭柘寺地区密集分布，震源深度主要为5—10km，南东侧地震分布较深，北西侧地震分布较浅，提示黄庄-高丽营断裂倾向为南东；黄庄-高丽营断裂北端震源深度以10—20km为主，揭示出断层倾向为南东，且断层倾角较大。本文研究结果为研究北京地区地震活动性与黄庄-高丽营断裂周边构造提供了重要的地震事件基础数据支撑。基于GPU加速的匹配定位技术能检测出大量台网地震目录中未记录的地震事件，遗漏事件以小于1级的微震为主，从而补充了台网地震目录。
致谢: 感谢北京市地震局为本研究提供的连续波形数据。本研究得到北京市自然科学基金（8212041）与上海佘山地球物理国家野外科学观测研究站开发基金（2020K02）的资助。 -
图 2 唐山MS7.8地震后MS≥5.0地震的M-T和N-T图（1976-07-28—2023-12-31）
Figure 2. The M-T and N-T plots of earthquakes with MS≥5.0 after the Tangshan MS7.8 earthquake
图 3 最大余震与主震时间间隔-频次图
Figure 3. Time frequency plot of the interval time between the maximum aftershock and the main shock
图 4 河北及邻区地震类型分区及震中分布图
Figure 4. Division of sub-regions and epicenters of earthquake types in Hebei and its surrounding regions
表 1 河北及邻区MS≥3.5主余型地震（1970—2022）
Table 1. Mainshock-aftershock earthquakes of MS≥3.5 in Hebei and its surrounding regions （1970—2022）
序号 发震时间
（年-月-日）纬度/（°N） 经度/（°E） 震级MS 最大余震
震级MS最大余震与主震
时间间隔/d主震与最大
余震震级差主震与最大
余震的震中距/km震中位置 1 1970-01-14 37.30 114.88 3.6 3.2 54.0 0.4 3.3 河北隆尧 2 1970-06-05 37.55 114.57 3.6 2.1 0.1 1.5 0.9 河北临城 3 1971-02-15 37.63 115.07 3.6 3.1 29.0 0.5 24.7 河北宁晋 4 1971-02-22 36.50 113.75 3.6 2.3 1.0 1.3 25.4 河北涉县 5 1971-12-27 38.42 114.42 4.6 2.7 0.0 1.9 0.0 河北行唐 6 1972-09-07 37.30 114.82 4.1 3.2 91.0 0.9 24.4 河北隆尧 7 1974-06-06 37.60 115.13 4.4 3.8 42.0 0.6 19.0 河北宁晋 8 1974-07-18 37.43 115.15 3.8 2.0 6.0 1.8 26.6 河北新河 9 1974-10-03 39.68 118.77 3.8 2.2 0.0 1.6 0.0 河北滦县 10 1976-07-28 39.63 118.18 7.8 7.1 0.6 0.7 46.9 河北唐山 11 1977-10-14 36.82 114.32 3.7 2.2 11.0 1.5 9.7 河北武安 12 1978-05-19 40.50 116.57 4.0 2.6 6.0 1.4 4.8 北京怀柔 13 1980-03-28 40.10 119.40 4.0 3.8 7.0 0.2 14.6 河北青龙满族自治县 14 1981-02-08 39.70 118.78 4.2 3.9 35.0 0.3 51.5 河北昌黎 15 1981-11-09 37.43 114.97 5.8 4.3 78.0 1.5 9.5 河北隆尧 16 1981-12-06 36.87 114.07 4.1 3.0 27.0 1.1 21.6 河北武安 17 1982-01-26 37.40 114.87 4.3 2.7 52.0 1.6 14.9 河北隆尧 18 1982-03-08 39.87 118.67 4.4 3.7 22.0 0.7 16.4 河北迁安 19 1982-05-29 36.58 114.57 4.3 2.7 3.0 1.6 22.3 河北邯郸 20 1982-10-19 39.88 118.98 4.9 2.6 0.0 2.3 1.1 河北卢龙 21 1982-12-10 40.47 116.55 4.4 3.1 0.0 1.3 2.1 北京怀柔 22 1983-10-05 39.85 118.90 4.2 2.7 10.0 1.5 18.9 河北卢龙 23 1984-01-07 39.72 118.75 4.8 3.4 64.0 1.4 34.1 河北滦县 24 1985-01-26 37.37 114.77 3.7 2.9 89.0 0.8 11.1 河北隆尧 25 1985-05-22 39.83 118.53 4.2 2.6 0.2 1.6 6.7 河北滦县 26 1985-10-05 39.78 118.45 4.6 2.7 8.0 1.9 55.4 河北唐山 27 1985-11-21 40.08 115.83 4.2 2.9 0.0 1.3 2.8 北京门头沟 28 1986-01-11 37.43 114.90 3.6 2.4 29.0 1.2 35.9 河北隆尧 29 1986-02-15 37.75 115.23 4.3 4.0 3.0 0.3 7.1 河北辛集 30 1986-12-16 37.27 114.75 3.7 2.6 7.0 1.1 3.7 河北隆尧 31 1987-06-07 39.53 118.08 4.0 3.2 33.0 0.6 48.2 河北丰南 32 1987-07-16 39.78 118.63 4.0 3.1 13.0 0.9 26.8 河北滦县 33 1988-10-26 39.87 118.48 3.8 3.0 78.0 0.8 25.7 河北迁安 34 1989-06-18 39.67 118.32 3.8 2.3 1.0 1.5 0.8 河北唐山 35 1989-07-23 37.37 115.03 3.8 2.9 80.0 0.9 5.3 河北巨鹿 36 1989-10-05 39.62 118.25 3.7 2.0 1.0 1.7 2.1 河北唐山 37 1989-12-14 37.60 115.28 4.3 2.7 68.0 1.6 32.1 河北冀县 38 1990-08-03 37.88 115.03 3.7 2.3 58.0 1.4 41.3 河北晋县 39 1991-05-07 39.80 118.70 3.8 2.9 0.0 0.9 5.9 河北滦县 40 1991-07-11 39.68 118.38 3.8 2.6 2.0 1.2 3.3 河北唐山 41 1991-07-27 39.85 118.80 3.9 2.3 10.0 1.6 14.9 河北滦县 42 1991-10-07 37.93 115.07 3.6 2.4 1.0 1.2 3.7 河北晋县 43 1991-10-17 39.73 118.42 3.8 3.3 79.0 0.5 5.6 河北唐山 44 1994-12-23 40.47 115.55 3.8 1.9 0.0 1.9 0.0 河北怀来 45 1995-07-20 40.27 115.43 3.6 2.9 1.0 0.7 3.5 河北涿鹿 46 1998-01-10 41.10 114.30 6.2 4.1 0.1 2.1 17.1 河北张北 47 1998-04-14 39.68 118.47 4.6 3.9 0.0 0.7 0.0 河北滦县 48 1999-05-12 40.18 118.78 3.7 3.1 74.0 0.6 50.5 河北迁安 49 1999-11-29 41.17 114.57 3.6 2.7 7.0 0.9 0.8 河北张北 50 2000-04-10 41.12 114.52 3.8 2.7 68.0 1.1 7.4 河北张北 51 2000-06-25 39.82 118.33 3.9 2.9 11.0 1.0 14.2 河北滦县 52 2001-12-28 39.75 118.68 3.7 2.0 17.0 1.7 16.9 河北滦县 53 2002-12-25 37.60 115.10 3.7 3.3 25.0 0.4 3.5 河北宁晋 54 2003-11-15 39.70 118.80 3.6 2.2 35.0 1.4 8.2 河北昌黎 55 2004-01-20 39.73 118.78 4.6 3.3 53.0 1.3 38.6 河北昌黎 56 2008-03-11 39.98 118.90 3.9 2.0 0.0 1.9 1.6 河北卢龙 57 2010-03-06 39.70 118.50 4.3 3.0 0.3 1.3 0.0 河北滦县 58 2010-04-09 39.47 118.07 4.1 2.4 0.1 1.7 1.1 河北丰南 59 2012-05-28 39.72 118.47 4.7 3.1 0.3 1.6 6.8 河北唐山 60 2012-06-18 39.62 117.57 4.0 3.4 69.0 1.6 11.9 天津宝坻 61 2015-09-14 39.73 118.82 4.2 3.3 0.0 0.9 1.8 河北昌黎 62 2016-09-10 39.68 118.33 4.1 2.8 25.0 1.3 2.8 河北开平 63 2017-09-04 37.50 114.35 3.8 2.0 0.0 1.8 6.6 河北临城 64 2021-04-16 39.75 118.72 4.3 2.9 69.0 1.4 9.2 河北滦州 65 2021-06-22 41.17 114.40 3.9 2.4 0.3 1.5 0.0 河北张北 66 2022-12-04 39.90 118.88 3.7 2.2 0.4 1.5 2.2 河北卢龙 
下载: 导出CSV
表 2 河北及邻区MS≥3.5多震型地震（1970—2022）
Table 2. Multiple mainshock earthquakes of MS≥3.5 in Hebei and its surrounding regions （1970—2022）
序号 发震时间
（年-月-日）纬度/（°N） 经度/（°E） 震级MS 最大余震
震级MS最大余震与
主震时间间隔/d主震与最大
余震震级差主震与最大余震的
震中距/km震中位置 1 1970-06-25 37.17 114.82 3.6 3.6 0.04 0 0 河北任县 2 1972-08-24 40.47 115.42 3.6 3.1 39.0 0.5 15.1 河北怀来 3 1980-08-07 37.52 114.42 3.8 3.6 0.04 0.2 7.9 河北赞皇 4 1981-09-06 40.32 115.85 3.9 3.3 1.0 0.6 3.4 河北怀来 5 1981-09-19 39.80 118.73 3.6 3.3 15.0 0.3 45.4 河北滦县 6 1982-01-27 39.75 118.77 3.9 4.4 40.0 0.5 15.6 河北昌黎 7 1983-01-31 37.43 115.02 3.6 4.0 6.0 0.4 11.5 河北宁晋 8 1983-08-21 39.72 118.32 3.6 3.0 0.42 0.6 4.2 河北唐山 9 1983-09-15 37.52 114.27 3.6 3.1 1.0 0.5 25.5 河北赞皇 10 1984-02-27 39.43 118.02 3.8 3.4 13.0 0.4 37.3 河北丰南 11 1984-12-13 39.65 118.22 3.6 3.3 31.0 0.3 21.0 河北唐山 12 1985-04-22 39.72 118.77 4.6 4.2 0.3 0.4 3.4 河北昌黎 13 1987-08-08 39.33 117.88 4.2 3.6 0.4 0.6 1.7 天津宁河 14 1988-07-25 39.52 118.10 3.9 4.4 1.0 0.5 5.5 河北丰南 15 1988-07-30 39.77 118.72 3.8 4.1 4.0 0.3 19.8 河北滦县 16 1990-03-16 39.77 118.45 3.7 3.2 35.0 0.5 7.8 河北唐山 17 1991-05-29 39.72 118.30 4.8 5.2 1.0 0.4 5.1 河北唐山 18 1991-08-21 37.33 114.70 3.8 3.6 1.0 0.2 15.8 河北隆尧 19 1996-12-16 40.17 116.50 4.0 3.4 0.2 0.6 10.3 北京顺义 20 2002-08-03 40.48 115.53 3.9 3.3 0.3 0.6 1.4 河北沙城 21 2003-04-23 39.33 117.63 3.6 3.8 0.2 0.2 2.9 天津宁河 22 2003-05-22 39.73 118.43 3.6 3.2 30.0 0.4 8.6 河北唐山 23 2006-05-03 39.77 118.45 3.6 3.8 0.01 0.2 0 河北唐山 24 2019-12-03 40.43 114.55 3.6 3.1 0.01 0.5 0 河北怀安
下载: 导出CSV
表 3 河北及邻区MS≥3.5孤立型地震（1970—2022）
Table 3. Isolated earthquakes of MS≥3.5 in Hebei and its surrounding regions （1970—2022）
序号 发震时间（年-月-日） 纬度/（°N） 经度/（°E） 震级MS 震中位置 1 1970-05-25 39.55 118.12 4.2 河北丰南 2 1971-08-05 38.37 114.52 4.3 河北行唐 3 1972-03-25 40.43 116.60 4.5 北京怀柔 4 1972-10-12 36.97 114.32 4.8 河北沙河 5 1973-05-13 39.63 117.37 3.8 天津宝坻 6 1973-09-21 39.05 116.33 4.5 河北霸州 7 1973-12-31 38.47 116.55 5.3 河北河间 8 1975-12-01 39.25 118.00 3.7 天津滨海新区 9 1976-04-22 38.68 116.67 4.3 河北大城 10 1977-02-04 39.97 117.13 4.0 河北三河 11 1977-03-10 36.30 114.30 4.0 河北磁县 12 1978-10-03 40.17 116.50 4.0 北京顺义 13 1979-05-31 40.20 119.13 3.7 河北青龙满族自治县 14 1980-11-05 39.32 116.78 3.6 河北廊坊 15 1980-11-10 40.65 115.63 3.7 河北赤城 16 1981-07-24 39.72 118.35 4.0 河北唐山 17 1981-08-03 39.50 117.90 3.8 天津宁河 18 1981-08-25 40.55 115.00 4.2 河北宣化 19 1981-10-16 39.30 117.85 3.6 天津滨海新区 20 1983-04-03 40.75 114.78 4.7 河北万全 21 1983-07-21 40.82 114.70 3.9 河北万全 22 1983-08-08 40.67 115.38 3.6 河北张家口 23 1983-08-26 39.80 117.22 3.8 天津宝坻 24 1985-10-17 39.28 114.65 3.6 河北涞源 25 1985-11-30 37.23 114.82 5.2 河北任县 26 1986-11-10 40.05 116.72 4.2 北京顺义 27 1987-03-21 38.28 114.25 4.0 河北平山 28 1987-11-11 40.28 114.80 4.2 河北宣化 29 1988-05-09 39.30 118.10 3.7 河北丰南 30 1988-07-23 40.08 114.22 4.6 河北阳原 31 1989-12-25 40.33 118.95 4.2 河北青龙 32 1990-05-23 40.22 116.47 3.8 北京昌平 33 1990-07-21 40.58 115.83 4.6 北京延庆 34 1990-07-23 39.75 118.48 4.5 河北唐山 35 1990-09-22 40.08 116.37 4.0 北京昌平 36 1990-12-24 37.90 115.02 3.6 河北晋县 37 1991-08-22 37.47 114.73 3.6 河北柏乡 38 1992-02-14 39.77 118.43 3.9 河北唐山 39 1992-07-22 39.28 117.93 4.5 天津滨海新区 40 1993-11-18 39.60 117.45 3.9 天津宝坻 41 1995-10-06 39.67 118.33 5.0 河北唐山 42 1997-05-25 40.70 114.87 4.2 河北宣化 43 2001-03-02 40.37 114.92 4.2 河北宣化 44 2001-07-16 38.32 116.47 3.8 河北沧县 45 2002-04-22 37.27 114.42 5.0 河北隆尧 46 2002-05-19 39.37 117.98 4.2 河北丰南 47 2006-07-04 38.90 116.28 5.1 河北文安 48 2014-09-06 40.27 115.42 4.3 河北涿鹿 49 2016-06-23 40.97 114.18 4.0 河北尚义 50 2017-01-02 40.77 114.15 3.6 河北怀安 51 2018-02-12 39.37 116.67 4.2 河北永清 52 2019-04-21 36.33 114.43 3.7 河北临漳 53 2019-12-5 39.31 118.04 4.5 河北丰南 54 2020-07-12 39.78 118.43 5.1 河北古冶 55 2021-07-15 37.57 115.08 3.6 河北宁晋 56 2022-02-06 36.97 115.43 3.6 河北威县 57 2022-05-21 39.97 119.00 3.6 河北卢龙 58 2022-10-03 38.37 113.73 4.3 河北平山
下载: 导出CSV
表 4 河北及邻区不同类型地震统计结果（1970—2022）
Table 4. Statistical results of different sequence types of earthquakes in Hebei and its surrounding regions （1970—2022）
序号 震级分档MS 主余型地震次数（占比） 多震型地震次数（占比） 孤立型地震次数（占比） 总次数 1 3.5～3.9 33（44%） 20（26%） 23（30%） 76 2 4.0～4.9 30（48%） 4（6%） 29（46%） 63 3 5.0～5.9 1（14%） 0（0） 6（86%） 7 4 6.0～6.9 1（100%） 0（0） 0（0） 1 5 7.0～7.9 1（100%） 0（0） 0（0） 1
下载: 导出CSV
表 5 河北及邻区各分区地震类型统计结果
Table 5. Statistical results of different sequence types of earthquakes in each partition of Hebei and its surrounding regions
区域 震级分档MS 主余型次数（占比） 多震型次数（占比） 孤立型次数（占比） 总次数 张家口区（Ⅰ） 3.5～4.9 5（42%） 3（25%） 4（33%） 12 4.0～4.9 0（0） 0（0） 9（100%） 9 5.0～5.9 0（0） 0（0） 0（0） 0（0） 6.0～6.9 1（100%） 0（0） 0（0） 1 7.0～7.9 0（0） 0（0） 0（0） 0（0） 北京区（Ⅱ） 3.5～4.9 0（0） 1（17%） 5（83%） 6 4.0～4.9 4（36%） 1（9%） 6（55%） 11 5.0～5.9 0（0） 0（0） 0（0） 0（0） 6.0～6.9 0（0） 0（0） 0（0） 0（0） 7.0～7.9 0（0） 0（0） 0（0） 0（0） 津唐区（Ⅲ） 3.5～4.9 14（44%） 11（34%） 7（22%） 32 4.0～4.9 18（64%） 3（11%） 7（25%） 28 5.0～5.9 0（0） 0（0） 2（100%） 2 6.0～6.9 0（0） 0（0） 0（0） 0（0） 7.0～7.9 1（100%） 0（0） 0（0） 1 冀中区（Ⅳ） 3.5～4.9 0（0） 0（0） 2（100%） 2 4.0～4.9 0（0） 0（0） 2（100%） 2 5.0～5.9 0（0） 0（0） 2（100%） 2 6.0～6.9 0（0） 0（0） 0（0） 0（0） 7.0～7.9 0（0） 0（0） 0（0） 0（0） 冀西区（Ⅴ） 3.5～4.9 0（0） 0（0） 0（0） 0（0） 4.0～4.9 1（25%） 0（0） 3（75%） 4 5.0～5.9 0（0） 0（0） 0（0） 0（0） 6.0～6.9 0（0） 0（0） 0（0） 0（0） 7.0～7.9 0（0） 0（0） 0（0） 0（0） 冀南区（Ⅵ） 3.5～4.9 14（58%） 5（21%） 5（21%） 24 4.0～4.9 7（78%） 0（0） 2（22%） 9 5.0～5.9 1（33%） 0（0） 2（67%） 3 6.0～6.9 0（0） 0（0） 0（0） 0（0） 7.0～7.9 0（0） 0（0） 0（0） 0（0）
下载: 导出CSV
表 6 河北及邻区MS≥5.0主余型地震（1970—2022）
Table 6. Mainshock-aftershock earthquake types of MS≥5.0 in Hebei and its surrounding regions （1970—2022）
序号 发震时间（年-月-日） 纬度/（°N） 经度/（°E） 震级MS 发震位置 最大余震震级MS 1 1976-07-28 39.63 118.18 7.8 河北唐山 7.1 2 1981-11-09 37.43 114.97 5.8 河北隆尧 4.3 3 1998-01-10 41.10 114.30 6.2 河北张北 4.1
下载: 导出CSV
表 7 研究震例的最大余震震级计算参数及结果
Table 7. Aftershock magnitude calculation parameters and results of studied earthquakes
序号 发震时间（年-月-日） 震级MS 发震位置 c q k MS（63%） MS（10%） MS（2%） 实际最大余震MS 1 1976-07-28 7.8 河北唐山 0.5833 0.0694 7.6556 5.6 6.3 6.8 7.1 2 1981-11-09 5.8 河北隆尧 1.1550 0.1750 8.46 3.9 4.4 4.8 4.3 3 1998-01-10 6.2 河北张北 2.6972 0.4026 9.8089 3.6 4.1 4.4 4.1
下载: 导出CSV
-
毕金孟,蒋长胜,曹付阳,2023. 震后早期阶段余震预测研究进展. 地球与行星物理论评(中英文),54(5):498−511.Bi J. M., Jiang C. S., Cao F. Y., 2023. Research progress of aftershock forecasting in the early stage after the mainshock. Reviews of Geophysics and Planetary Physics, 54(5): 498−511. (in Chinese) 崔效锋,谢富仁,李瑞莎等,2010. 华北地区构造应力场非均匀特征与煤田深部应力状态. 岩石力学与工程学报,29(S1):2755−2761.Cui X. F., Xie F. R., Li R. S., et al., 2010. Heterogeneous features of state of tectonic stress filed in North China and deep stress in coal mine. Chinese Journal of Rock Mechanics and Engineering, 29(S1): 2755−2761. (in Chinese) 国家地震局预测预防司,1997. 测震学分析预报方法. 北京:地震出版社,152−153. 郭蕾,孙丽娜,周依等,2021. 唐山地区近十年地震的重定位研究. 地震学报,43(6):706−715.Guo L., Sun L. N., Zhou Y., et al., 2021. Relocation of earthquakes in Tangshan area in the recent decade. Acta Seismologica Sinica, 43(6): 706−715. (in Chinese) 蒋海昆,李永莉,曲延军等,2006. 中国大陆中强地震序列类型的空间分布特征. 地震学报,28(4):389−398. doi: 10.3321/j.issn:0253-3782.2006.04.006Jiang H. K., Li Y. L., Qu Y. J., et al., 2006. Spatial distribution features of sequence types of moderate and strong earthquakes in Chinese mainland. Acta Seismologica Sinica, 28(4): 389−398. (in Chinese) doi: 10.3321/j.issn:0253-3782.2006.04.006 蒋海昆,王锦红,2023. 适用于机器学习的地震序列类型判定特征重要性讨论. 地震研究,46(2):155−172.Jiang H. K., Wang J. H., 2023. Discussion on the importance of the features for the judgement of earthquake sequence types applicable to machine learning. Journal of Seismological Research, 46(2): 155−172. (in Chinese) 宋美琴,李斌,吕芳等,2011. 山西及周边地区地震类型分区特征. 地震学报,33(5):691−696. doi: 10.3969/j.issn.0253-3782.2011.05.013Song M. Q., Li B., Lü F., et al., 2011. A study on sub-region division of earthquake types in Shanxi and its surrounding regions. Acta Seismologica Sinica, 33(5): 691−696. (in Chinese) doi: 10.3969/j.issn.0253-3782.2011.05.013 王华林,周翠英,耿杰,1997. 中国大陆及邻区地震序列类型的分区特征和震源环境讨论. 地震,17(1):34−42.Wang H. L., Zhou C. Y., Geng J., 1997. Discussion on zoning features of earthquake sequence type and source environment in Chinese mainland and its adjacent area. Earthquake, 17(1): 34−42. (in Chinese) 吴开统,1971. 地震序列的基本类型及其在地震预报中的应用. 地震站线,7(11):45−51. 许晓庆,盛菊琴,马禾青,2016. 中国西北地区东部地震序列参数与预测余震发生的关系分析. 地震工程学报,38(4):598−605. doi: 10.3969/j.issn.1000-0844.2016.04.0598Xu X. Q., Sheng J. Q., Ma H. Q., 2016. Relation between earthquake sequence parameters and the prediction of aftershocks in the eastern part of northwest China. China Earthquake Engineering Journal, 38(4): 598−605. (in Chinese) doi: 10.3969/j.issn.1000-0844.2016.04.0598 叶友清,苏金蓉,易桂喜等,2019. 利用震级均值特性建立有上限震级概率模型对2017年九寨沟7.0级地震的余震概率预测. 地球物理学进展,34(2):489−499. doi: 10.6038/pg2019CC0089Ye Y. Q., Su J. R., Yi G. X., et al., 2019. Probabilistic prediction for the earthquake sequence of 2017 Jiuzhaigou, Sichuan, M S7.0 by using limited-magnitude probability distribution based on the average magnitudes of aftershock sequence. Progress in Geophysics, 34(2): 489−499. (in Chinese) doi: 10.6038/pg2019CC0089 周翠英,侯海峰,华爱军等,2001. 山东及其近海区地震序列类型地理分布特征. 华北地震科学,19(3):15−27. doi: 10.3969/j.issn.1003-1375.2001.03.002Zhou C. Y., Hou H. F., Hua A. J., et al., 2001. The geographic distribution feature of earthquake sequence type in Shandong and its offshore area. North China Earthquake Sciences, 19(3): 15−27. (in Chinese) doi: 10.3969/j.issn.1003-1375.2001.03.002 Ogata Y., 1989. Statistical model for standard seismicity and detection of anomalies by residual analysis. Tectonophysics, 169(1-3): 159−174. doi: 10.1016/0040-1951(89)90191-1 Page M. T., van der Elst N., Hardebeck J., et al., 2016. Three ingredients for improved global aftershock forecasts: tectonic region, time-dependent catalog incompleteness, and intersequence variability. Bulletin of the Seismological Society of America, 106(5): 2290−2301. doi: 10.1785/0120160073 Reasenberg P. A., Jones L. M., 1989. Earthquake hazard after a mainshock in California. Science, 243(4895): 1173−1176. doi: 10.1126/science.243.4895.1173 期刊类型引用(0)
其他类型引用(4)
-
下载:
点击查看大图
计量
- 文章访问数: 86
- HTML全文浏览量: 29
- PDF下载量: 15
- 被引次数: 4
