2020年3月23日新疆拜城5.0级和 7月13日霍城5.0级地震总结

2021-07-03苑争一宋治平姜祥华马亚伟解孟雨韩颜颜

地震地磁观测与研究 2021年2期

苑争一宋治平姜祥华马亚伟姚琪解孟雨于晨韩颜颜

（中国北京100045 中国地震台网中心）

0 引言

据中国地震台网测定，天山中部地区在2020 年3 月和7 月各发生1 次5.0 级地震，分别为3 月23 日拜城5.0 级（41.75°N，81.11°E）和7 月13 日霍城5.0 级地震（44.42°N，80.82°E）。其中，拜城5.0 级地震震中位于南天山地震带中段的拜城凹陷东北端，霍城5.0 级地震震中位于北天山地震带西段的伊犁盆地北缘断裂（喀什河断裂）附近。余震序列显示，3 月23 日拜城5.0 级地震余震较少，截至5 月23 日，共记录ML≥1.0 余震2 次，序列主震与最大余震的震级差为2.9，主震释放能量占序列的比例为99.99%，是一次孤立型事件；7 月13 日霍城5.0 级地震余震不多，截至9 月13 日，共记录ML≥1.0 余震9 次，序列主震与最大余震的震级差为2.5，主震释放能量占序列的比例为99.99%，也是一次孤立型事件。

文中介绍了拜城5.0 级和霍城5.0 级地震的区域构造背景和余震序列，详细梳理了2次地震发生前的活动异常，追踪了震前地震学参数演化特征，结合形变和电磁2 类地球物理观测手段震前异常现象，分析了各观测手段在时间、空间和强度方面的预测意义，以期为新疆地区5 级以上地震的判定积累震例资料。

1 历史地震和构造背景

1.1 历史地震

在2020 年3 月23 日拜城5.0 级地震震中100 km 范围内，1900 年以来共发生10 次5级以上地震，其中5.0—5.9 级地震7 次，6.0—6.9 级地震2 次，7.0—7.9 级地震1 次，1906 年3 月2 日新疆阿克苏地区温宿县7.0 级地震震级最大，距拜城5.0 级地震约90 km。在2020 年7 月13 日霍城5.0 级地震震中100 km 范围内，1900 年以来共发生10 次5 级以上地震（含余震），其中5.0—5.9 级地震8 次，6.0—6.9 级地震2 次，1958 年12 月21 日新疆温泉6.5 级地震震级最大，距霍城5.0 级地震约15 km（图1）。

图1 2020 年3 月23 日拜城5.0 级和7 月13 日霍城5.0 级地震震中附近历史地震分布Fig.1 The historical earthquake distribution near the epicenter of Baicheng MS5.0 earthquake on March 23 and Huocheng MS5.0 earthquake on July 13,2020,respectively

1.2 构造背景

天山地震带位于南部的塔里木盆地和帕米尔高原与北部的准噶尔盆地和哈萨克地台之间，是亚洲内部的特大造山带之一，也是世界上年轻、活跃的陆内造山带之一（蔡妍等，2019；龚正等，2020）。天山中部地区位于天山造山带中北部，展布于刚性的塔里木盆地和准噶尔盆地之间，受印度大陆和欧亚大陆碰撞远程效应的影响，区内分布数条近EW 向逆冲—褶皱活动构造带，是大陆内部变形和强震活动频繁的地区之一（Tapponnier et al，1979；沈军等，2003；Campbell et al，2013；张志斌等，2020）。

（1）2020 年拜城5.0 级地震。此次地震发生在南天山地震带中段拜城凹陷东北端。拜城凹陷是南天山中段与塔里木盆地北缘之间库车褶皱—冲断带构造单元的重要组成部分。库车褶皱—冲断带发育多种盐构造样式，盐构造空间上成排成带分布，自北向南主要形成克拉苏构造带、拜城凹陷和秋里塔格构造带（汪新等，2009）。克拉苏构造带为变形后缘，拜城凹陷为变形过渡带，而秋里塔克构造带为变形前锋带，盐构造主要发育于克拉苏和秋里塔格构造带，拜城凹陷受到微弱的构造挤压变形（唐鹏程等，2012）。库木格热木断裂与此次地震震中距离最近，相距约7 km，是克拉苏构造带的一部分，以逆冲为主（图2）。

图2 拜城5.0 级和霍城5.0 级地震和历史地震震源机制解（据GCMT，https://www.globalcmt.org/CMTsearch.html）Fig.2 The focal mechanisms for Baicheng MS5.0 earthquake,Huocheng MS5.0 earthquake,and historical earthquakes (according to GCMT)

（2）2020 年霍城5.0 级地震。此次地震发生在北天山地震带西段伊犁地块北缘。伊犁地块隶属于天山造山带伊犁—中天山微板块，主要由伊犁盆地构成（朱志新等，2013）。伊犁盆地是在塔里木板块和哈萨克斯坦板块南北对冲挤压应力作用下形成的大型内陆山间拗陷盆地，中生代受燕山运动和喜马拉雅运动改造作用明显，盆山接触地段发生显著逆冲推覆，新生代以来，印度板块向欧亚板块的持续俯冲、碰撞，引起青藏高原广泛强烈抬升和差异升降，对先期构造产生叠加改造。伊犁盆地北缘断裂（喀什河断裂带）作为博罗科努山与伊犁盆地分界断裂，控制了盆地北部边界，断裂地貌标志明显，北部为高山，南部为倾斜平原，相对高差达1 000 m 以上。野外地震地质勘查结果显示，晚第四纪以来，该断裂活动速率相对较低，垂直运动速率为0.9—1.2 mm/a（尹光华等，2003）。距霍城5.0 级地震最近的断层即为喀什河断裂，以右旋逆冲为主，距离约4 km（图2）。

震源机制解显示，拜城5.0 级地震为一次走滑型破裂地震，震中附近历史地震震源机制以近EW 向逆冲型破裂为主，优势主压应力轴方位近NS 向；霍城5.0 级地震为一次走滑型破裂地震，震中附近历史地震震源机制以NE 向具有一定逆冲分量的左旋走滑破裂为主，优势主压应力轴方位近NS 向。2020 年拜城5.0 级和霍城5.0 级地震所在的天山地震带中段，历史地震以近东向的逆冲型破裂为主，优势主压应力轴方位近NS 向（图2，表1）。

表1 拜城5.0 级和霍城5.0 级地震震源机制Table 1 The focal mechanism solutions for the Baicheng MS5.0 earthquake and Huocheng MS5.0 earthquake

2 地震活动特征

2020 年上半年，新疆地区5 级以上地震活跃，在1 月16 日库车5.6 级地震至7 月13日霍城5.0 级地震发生期间，天山地震带共发生9 次5 级以上地震（图1），地震活动较往年显著增强。本研究系统总结3 月23 日拜城5.0 级和7 月13 日霍城5.0 级地震发生前的活动特征，结果显示，拜城地震前具有明确预测意义的异常较少，仅有天山中段5 级以上地震成组1 项异常；霍城地震前具有3 级以上地震带状分布、地震发生率指数、地震危险度D值、震前小震调制比和b值等明确预测意义的异常。

2.1 地震活动异常

2.1.1 拜城地震。拜城地震前天山中部地区存在5 级以上地震成组异常。以3 个月时间内发生2 次5 级以上地震为成组标准，统计天山中部地区1970 年以来发生的5 级以上地震，结果表明，共发生50 次5 级以上地震，其中参与成组的地震有17 次，占比34%（17/50），共有8 组（表2）。2020 年1 月16 日新疆库车发生5.6 级地震，根据统计的成组规律，天山中部地区3 个月内存在发生5 级以上地震的可能，综合分析认为，2020 年拜城5.0 级地震为对应地震。

表2 天山中部5.0 级以上地震成组统计Table 2 Group list of MS≥5.0 earthquakes in the middle of Tianshan

2.1.2 霍城地震。（1）震前3 级以上地震带状分布。地震活动条带图像是地震预测的常用方法之一，表征了地震前区域地震活动由凌乱分散分布转为集中成带的现象（刘蒲雄，1992）。2020 年5 月10 日至7 月12 日，新疆伽师至克拉玛依出现3 级以上地震带状分布图像，根据条带起始震级与后续主震的一般关系（李莹甄等，2011），分析认为，该条带及周边后续存在发生5 级以上地震的可能，2020 年7 月13 日新疆霍城5.0 级地震即发生于此（图3）。

图3 霍城5.0 级地震前3 级以上地震带状分布Fig.3 Seismic belt with a magnitude of MS≥3.0 before Huocheng MS5.0 earthquake

（2）震前地震发生率指数异常演化。研究表明，强震发生前普遍存在地震活动频度增强与平静现象（陆远忠等，1997；梅世蓉等，1997），基于统计学模型构建的地震发生率指数，对地震活动显著增强与显著平静2 种典型异常具有定量识别的能力（姜祥华，2019）。采用5 年背景窗长，3个月计算窗长，进行地震发生率时空扫描。基于泊松分布，将地震发生率指数的显著增强或显著减弱变化统一映射为0—1 之间的概率值，并以色标显示，红色表示地震发生率指数显著增强，蓝色表示地震发生率指数显著减弱。2020 年7 月13 日霍城5.0 级地震发生前35—15 天，震源区附近地震发生率指数持续高值异常，并于震前10 天消失，而此次地震即发生在该显著高值异常区附近（图4）。

（3）震前地震危险度D值异常演化。本研究为考察地震危险度D值（王炜，1987）的时空演化特征，对D值计算进行一定改进，主要流程如下：①采用5 年背景窗长计算D值的均值及标准差；②采用1 年窗口计算当前时段的D值；③将当前D值减去背景均值后除以标准差，得到D值的相对变化（ΔD/σD）。空间扫描的网格大小取为1°×1°，计算此次5.0 级地震前的D值变化，结果见图5，可见震前5 个月，D值在震源区及附近出现高值异常，持续至7 月13 日霍城5.0 级地震发生（图5）。

（4）震前调制比异常。地震调制比异常可用于间接探测地壳介质的强弱分布，寻找高应力集中区，是地震预测研究的常用方法之一（韩颜颜等，2017；孟令媛等，2020）。霍城5.0 级地震发生前约10 个月，震中西部附近地区出现调制比高值异常，2019 年12 月底，该异常区向NE 发生调整变化，且异常值逐步升高，至2020 年2 月底达最高，后逐步减小，且该异常区向外扩展，至2020 年6 月底，在震中附近出现调制比高值异常，并于2020 年7 月13 日发生此次5.0 级地震。因此，该地震发生前，调制比异常持续存在，且呈现持续衰减过程，而在空间上，异常区表现出逐步向震中附近扩展的变化（图6）。

图6 霍城5.0 级地震前固体潮调制比时空演化过程(a) 2019-05-01—2019-10-31；(b) 2019-06-01—2019-11-30；(c) 2019-07-01—2019-12-31；(d) 2019-08-01—2020-01-31；(e) 2019-09-01—2020-02-29；(f) 2019-10-01—2020-03-31；(g) 2019-11-01—2020-04-30；(h) 2019-12-01—2020-05-31；(i) 2020-01-01—2020-06-30Fig.6 The spatio-temporal evolution of earth tidal modulation ratio before Huocheng MS5.0 earthquake

（5）震前b值异常。岩石力学实验显示，b值大小与应力水平呈反比关系（Scholz，1968；Schorlemmer et al，2005a），实际地震数据计算结果显示，在凹凸体区域的b值一般较低（Schorlemmer et al，2005b；易桂喜等，2008）。因此，计算区域b值，并分析相应的时空演化特征，是地震预测研究的常用方法之一（闻学泽等，2013；史海霞等，2018）。霍城5.0级地震发生前1 年，震中附近出现低b值异常，具体表现在：2019 年11 月至2020 年1 月，震中附近b值逐步升高，后转为降低，出现低值异常现象，同时异常区逐步扩大，2020 年4 月后，b值持续降低且异常区范围逐月扩大，至震前达最低，7 月13 日发生霍城5.0 级地震。可见，该地震发生前，b值表现为起伏变化，呈减小—增大—减小的变化特征，而相应的低b值异常则表现为增强—减弱—增强的循环变化特征，同时在空间上表现为异常区在震中附近的持续出现（图7）。

图7 霍城5.0 级地震前b值时空演化过程(a)2019-02-12—2019-08-12；(b)2019-03-12—2019-09-12；(c)2019-04-12—2019-10-12；(d) 2019-05-12—2019-11-12；(e)2019-06-12—2019-12-12；(f)2019-07-12—2020-01-12；(g)2019-08-12—2020-02-12；(h)2019-09-12—2020-03-12；(i)2019-10-12—2020-04-12；(j)2019-11-12—2020-05-12；(k)2019-12-12—2020-06-12；(l)2020-01-12—2020-07-12Fig.7 The spatio-temporal evolution of b-value before Huocheng MS5.0 earthquake

2.2 地球物理观测异常

统计结果显示，在2020 年3 月23 日拜城5.0 级地震和7 月13 日霍城5.0 级地震震中300 km 范围内，分布的地球物理观测台站分别为23 个和19 个，布设地倾斜、地应变、地磁、地电阻率、水位、水温、水化学等观测项。2 次震前共有10 个测项出现地球物理异常（图8，表3），含1 项重叠测项，其中：①拜城地震震中300 km 范围内共有4 项异常，分别为乌什伸缩仪EW 分量、尼勒克钻孔应变NE 和EW 分量、柯坪地电阻率NS 分量；②霍城地震震中300 km 范围内共有7 项异常，分别为尼勒克钻孔应变NE 分量、果子沟钻孔倾斜NS、EW 分量、温泉石英摆倾斜NS 分量、巩留钻孔倾斜NS、EW 分量、精河石英水平摆EW 分量。

表3 拜城5.0 级和霍城5.0 地震震中附近地球物理观测异常信息Table 3 Observed geophysical anomaly information of Baicheng MS5.0 earthquake and Huocheng MS5.0 earthquake,respectively

（1）乌什伸缩仪EW 分量。乌什伸缩EW 分量测项位于拜城5.0 级地震震中西南方向，相距约169 km。该测项自2012 年8 月21 日起由趋势压缩转为拉张，截至2020 年3 月23日拜城5.0 级地震发生前，异常持续2 700 多天，为长趋势异常（图9）。以往地震前该分量多次出现短期快速、大幅拉张变化，拜城5.0 级地震前无类似突出变化，且震后并无转向恢复迹象，故暂不作为此次地震异常，但仍具有预测意义。

图9 乌什伸缩仪EW 分量异常变化形态Fig.9 The EW component anomaly characteristic of cave strain at Wushi station

（2）尼勒克钻孔应变NE 和EW 分量。尼勒克钻孔应变测项位于拜城5.0 级地震震中NE 方向，相距约253 km；位于霍城5.0 级地震震中SE 方向，相距约162 km。NE 分量自2019 年10 月1 日出现反向拉张变化，至2020 年7 月底仍在持续，暂不作为拜城5.0 级和霍城5.0 级地震异常，需继续关注其后续变化。EW 分量在2019 年10 月23 日至2020 年2月15 日出现小幅度反向拉张的异常变化，后恢复压缩变化，异常结束（图10），分析认为，可以作为拜城5.0 级地震的映震异常。

图10 尼勒克钻孔应变异常变化形态(a) NE 分量；(b) EW 分量Fig.10 The anomaly characteristic of borehole strain at Nilka station

（3）果子沟钻孔倾斜NS 和EW 分量。果子沟钻孔倾斜测项位于霍城5.0 级地震震中SE方向，相距约14 km。震例回溯表明，2017 年8 月9 日精河6.6 级地震（震中距154 km）发生前，该测项NS 和EW 分量出现加速变化，异常持续时间分别为109 天和41 天。2020 年4 月中旬开始，NS 和EW 分量同步加速变化，并分别持续约15 天和5 天（图11）。中国地震台网中心预报部形变研究室周、月会商跟踪分析意见指出，异常结束后2—5 个月内，台站周边300 km 范围内有发生5.0 级以上地震的可能，此次霍城5.0 级地震即发生在异常结束后约2.5 个月。因此，NS 和EW 分量的同步加速变化应为此次霍城5.0 级地震异常。

图11 果子沟钻孔倾斜异常变化形态(a)NS 分量；(b)EW 分量Fig.11 The anomaly characteristic of borehole tilt at Guozigou station

（4）温泉石英摆倾斜NS 分量。温泉石英摆倾斜测项位于霍城5.0 级地震震中北侧，相距约55 km。2016 年开始，该测项由趋势南倾转折为趋势北倾，且年变形态发生较大变化，截至霍城地震发生前，异常已持续4 年多，震后无转折回返迹象（图12）。分析认为，此长趋势转折异常与霍城5.0 级地震无显著关联性，不作为该地震异常。

图12 温泉石英摆倾斜NS 分量异常变化形态Fig.12 The NS component anomaly characteristic of quartz pendulum tilt at Wenquan station

（5）巩留钻孔倾斜NS 和EW 分量。巩留钻孔倾斜测项位于霍城5.0 级地震震中东南侧，相距约170 km。2020 年6 月15 日起，NS、EW 分量同步快速变化，其中EW 分量变化幅度达5.5″，异常变化显著（图13）。经新疆维吾尔自治区地震局预报中心和伊犁州地震监测中心现场核实，排除周边环境、观测条件、气象因素等的影响，综合判定该变化真实可靠，并作为短期异常在周月地震会商中进行跟踪，发现异常结束25 天后发生霍城5.0 级地震。分析认为，巩留钻孔倾斜快速变化为此次5.0 级地震异常。

图13 巩留钻孔倾斜异常变化形态(a) NS 分量；(b) EW 分量Fig.13 The anomaly characteristic of the borehole tilt at Gongliu station

（6）精河石英摆倾斜EW 分量。精河石英摆倾斜测项位于霍城5.0 级地震震中东南侧，相距约170 km。2016 年9 月23 日起，该测项EW 分量由趋势西倾转为趋势东倾（图14），目前已持续4 年多，趋势东倾时段发生多次中强地震，包括2016 年12 月8 日呼图壁6.2 级（距离280 km）、2017 年8 月9 日精河6.6 级（距离39 km）、2018 年10 月16 日精河5.4 级（距离59 km）等地震。此次霍城5.0 级地震前后，该异常持续东倾，速率及趋势均无短期变化，分析认为非此次地震异常。

图14 精河石英摆倾斜EW 分量异常变化形态Fig.14 The EW component anomaly characteristic of the quartz pendulum tilt at Jinghe station

（7）柯坪地电阻率NS 分量。柯坪地电阻率NS 和EW 测道自2018 年出现同步下降变化，幅度分别为-1.4%和-1.7%，且NS 测道存在年变幅减小现象。根据以上地电阻率异常特征，2020 年度电磁学科会商会提出“南天山中段存在发生6 级左右地震可能”的预测意见，并将该区域划为年度危险区。2019 年下半年，柯坪地电阻率NS 测道年变幅度减小变化显著，2020 年初达最小，随后转折回升，于2020 年3 月23 日发生拜城5.0 级地震（震中距220 km），后观测数据开始正常的年变转向，近期无加速变化（图15），认为非此次地震异常。

图15 柯坪地电阻率NS 异常变化形态Fig.15 The NS component anomaly characteristic of apparent resistivity at Kalpin station

3 地震序列特征

（1）拜城5.0 级地震序列。2020 年3 月23 日至5 月23 日，此次地震震区共记录ML≥1.0 地震3 次，其中ML2.0—2.9 地震2 次，ML5.0—5.9 地震1 次，最大余震为4 月16 日ML2.5 地震（表4）。拜城地震后2 个月内发生ML≥2.0 余震2 次，与主震最小震级差为2.9，主震释放能量占全部序列的99.99%，为一次孤立型事件。从空间分布看，2 次余震位于主震西北侧，最大余震距主震约3.4 km（图16）。对于该地震序列，由于仅有3 次地震，无法计算相应的b值和h值参数。

表4 拜城5.0 级地震序列目录（ML≥2.0）Table 4 Catalogue of the Baicheng MS5.0 earthquake sequence (ML≥2.0)

图16 拜城5.0 级地震序列震中分布(a)、地震序列M—t图(b)及地震序列N—t图(c)Fig.16 Epicentral distribution (a),M-t(b) and N-tdiagram (c) of the Baicheng MS5.0 earthquake sequence

（2）霍城5.0 级地震序列。2020 年7 月13 日至9 月13 日，此次地震震区共记录ML≥1.0 地震10 次，其中ML1.0—1.9 地震4 次，ML2.0—2.9 地震5 次，ML5.0—5.9 地震1 次，最大余震为7 月13 日ML2.9 地震（表5）。霍城5.0 级地震序列主震与最大余震的震级差为2.5，主震释放能量占全部序列的99.99%，属于孤立型地震事件。从空间分布看，余震活动集中在主震附近20 km 范围内，呈NS 向展布，其长轴约30 km，短轴约13 km，最大余震距主震约3 km（图17）。对于该地震序列，因仅有10 次地震，无法计算相应的b值和h值参数。

表5 霍城5.0 级地震序列目录（ML≥2.0）Table 5 Catalogue of the Huocheng MS5.0 earthquake sequence (ML≥2.0)

4 讨论与结论

2020 年1 月1 日至7 月13 日，新疆天山地震带共发生5.0 级以上地震9 次，包括1 次 6.0级以上地震，5.0 级以上地震活动水平显著高于往年同期水平（1950 年以来，每年上半年5.0 级以上地震平均活动水平为3 次）。本文重点梳理了拜城5.0 级和霍城5.0 级地震前的中短期异常，在地球物理异常中，存在一项异常对应多个地震的情况，例如尼勒克钻孔应变EW 分量，该测项在2020 年1 月16 日库车5.6 级地震前出现“高频异常”和“反向变化异常”，地震发生后“高频异常”消失，“反向变化”持续，直至2 月底该异常变化基本恢复，异常结束后约1 个月于3 月23 日发生拜城5.0 级地震，并于7 月13 日发生天山地震带2020 年5.0级地震活跃过程中最后一个地震——霍城5.0 级地震，目前仍有若干异常尚未恢复或对应地震。因此，认为天山中部地区仍存在发生中强地震的可能。

地震序列演化结果显示，拜城5.0 级和霍城5.0 级地震均为孤立型地震事件，且2 次地震余震记录较少，均无法计算b值和h值等序列参数，综合分析得出以下结论。

（1）2 次5.0 级地震均发生在天山中部地区，前者位于拜城凹陷东北端，后者位于伊犁盆地北缘，其中：拜城5.0 级地震位于以逆冲为主的库木格热木断裂，震源机制解显示，该地震为走滑型破裂，与发震断裂性质存在一定差异；霍城5.0 级地震位于以右旋逆冲运动为主的喀什河断裂，震源机制解显示，该地震为走滑型破裂，与断裂性质存在一定差异。

（2）2 次5.0 级地震后余震记录均较少，其中：拜城5.0 级地震后2 个月仅发生ML≥2.0 地震2 次，最大余震为ML2.5 地震；霍城5.0 级地震后2 个月共记到ML≥1.0 地震10 次，最大余震为ML2.9 地震。2 次地震主震释放能量占比均超过99%，均属于孤立型地震事件。

（3）2 次5.0 级地震前的地震活动和地球物理观测异常差异较大，其中：拜城5.0 级地震前只存在1 项地震活动异常和4 项地球物理异常，分别为天山中部5 级以上地震成组、乌什伸缩仪EW 分量、尼勒克钻孔应变NE 和EW 分量、柯坪地电阻率NS 分量；霍城5.0级地震前具有中短期预测意义的地震活动和地球物理异常较为丰富，地震活动方面存在5项异常，包含3 级以上地震带状分布、地震发生率指数、D值、b值、调制比异常，而地球物理观测异常共7 项，有尼勒克钻孔应变NE 分量、果子沟钻孔倾斜NS 和EW 分量、温泉石英摆倾斜NS 分量、巩留钻孔倾斜NS 和EW 分量、精河石英水平摆EW 分量，趋势和短期变化显著，短期异常以形变为主，其中果子沟和巩留钻孔倾斜NS 和EW 分量在震前3 个月内均出现突出的短期变化，为霍城地震的发生提供了判定依据。

本文撰写得到王海涛研究员和刘杰研究员的指导和鼓励，蒋海昆研究员、晏锐研究员、孟令媛研究员和闫伟高级工程师亦给予帮助。中国地震台网中心国家地震科学数据中心（https://data.earthquake.cn）提供数据支撑，在此对他们及中国地震台网中心预报部同事的辛苦工作，一并表示衷心感谢。