2020年5月18日云南巧家Ms5.0地震总结

admin管理员组

文章数量:1530842

2024年6月8日发(作者：)

第41卷 第5期

2020年 10月

地震地磁观测与研究

SEISMOLOGICAL AND GEOMAGNETIC

OBSERVATION AND RESEARCH

Vol.41 No.5

Oct. 2020

doi: 10. 3969/j. issn. 1003-3246. 2020. 05. 021

2020年5月18日云南巧家

M

S

5.0地震总结

卢 显 解孟雨 张小涛 王 月 姜祥华 杨 文

姚 琪 张凌空 解 滔 田 雷 薛 艳

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

摘要

系统总结2020年5月18日云南巧家

M

S

5.0地震发生前出现的地震活动和地球

物理观测异常，其中：

①

地震活动异常：震前存在震群活动、Benioff应变比和地震发

生率异常等中短期异常；

②

地球物理观测异常：震中200 km范围内分布重力、地磁、

水位、水温、断层氢等测项，其中异常测项15项，大部分为趋势异常，异常台项比较

低（5%），短临异常不突出，巧家电磁波为地震的短期异常；

③

综合方法：多参数概率

谱计算结果在震前出现异常。对巧家

M

S

5.0地震序列进行跟踪研究，发现：主震与最

，主震释放能量占序列总能量的99.64%，属主余型

大余震的震级差为1.6（

M

L

震级）

地震序列；主震震源机制显示为走滑型，根据余震分布，认为发震断层位于莲峰断裂

和昭通—鲁甸断裂之间的NS向次级断裂。综合分析认为，在巧家

M

S

5.0地震发生前，

地震活动具有背景异常逐渐向中短期转变的过程，地球物理观测异常以趋势异常为主，

短期异常不突出。

关键词

巧家

M

S

5.0地震；地震活动异常；地球物理观测异常；序列跟踪

0 引言

据中国地震台网测定，2020年5月18日21时47分在云南省昭通市巧家县（27.18°N，

103.16°E）发生

M

S

5.0地震，截至2020年5月19日6时，造成4人死亡（巧家县3人、

鲁甸县1人）、23人受伤（/sh/2020/05-19/）。

此次巧家5.0级地震震源深度8 km，震中位于安宁河—则木河断裂带，与震中最

近的断层为荞麦地—大箐断裂，距离约14 km，与2020年度“四川道孚至云南永善7.0

级左右年度地震重点危险区”相距约51 km，同时位于2020年5—7月全国地震概率预

测图所示西南地区相对高概率区附近。中国地震台网中心给出的震源机制解显示，该地

震具有走滑型破裂性质，其中NS向节面走向与大凉山断裂至包谷垴—小河断裂的连线

一致，显示此次地震可能是一次NS向次级断裂活动。

分析此次巧家5.0级地震发生前出现的各项异常，发现地震活动具有一定预测意义

的中期和短期震前异常；震中附近的地球物理观测以趋势异常为主，短期异常不突出；

在综合计算方法中，多参数概率谱显示出短期异常。文中就2020年巧家

M

S

5.0地震发

生前的震群、Benioff应变比和地震发生率变化及形变、电磁、流体等地球物理观测异常

作者简介：卢显（1985—），女，博士，高级工程师，主要从事数字地震学和地震遥感应用的研究。

E

-

mail：luxian@

基金项目：国家自然科学基金青年基金（项目编号：41704062）

本文投稿日期：2020

-

08

-

30

154地震地磁观测与研究41卷

与多参数概率谱方法在震前体现的异常特征进行了分析，以期通过对此次地震的系统回

顾和总结，为滇东北地区5级地震活动积累资料，为震前预报分析提供基础数据和震例

资料，进一步为防灾减灾救灾服务。

1 构造背景和历史地震

2020年巧家

M

S

5.0地震发生在巴颜喀拉块体、川滇块体与华南块体三者之间，以挤

压和左旋走滑为主要活动特征的大凉山次级块体内部。该次级块体吸收了来自川滇块体

和巴颜喀拉块体的挤压作用，以各边界断裂带的挤压和内部大凉山断裂带、峨边断裂带

等NNW向左旋走滑次级断裂为主要活动特征（程佳等，2014）。

此次巧家5.0级地震震中位于小江断裂以东、莲峰断裂与西鱼河—昭通断裂之间，

是发生在大凉山次级块体南缘的一次中等强度的破坏性地震。该地区除发育安宁河断裂、

则木河断裂、小江断裂、大凉山断裂、马边—盐津断裂等主干活动断层外，还存在NE

向昭通—莲峰断裂带和昭通—鲁甸断裂带（闻学泽等，2013）。该区域近年来中强地震

多发，先后发生了2005年8月5日会泽

M

S

5.3地震、2012年9月7日彝良

M

S

5.7、5.6

地震、2014年8月3日鲁甸

M

S

6.5地震（李文涛等，2020）。1900年以来该区历史地震

分布见图1，图中沙滩球为巧家

M

S

5.0地震震源机制解（据中国地震台网中心）。

102°103°104°E

甘

洛

—

竹

核

断

裂

莲

峰

断

裂

28°N

断

裂

小

江

2020-05-18

M

S

5.0

通

断

裂

鱼

河

—

昭

27°

2014-08-03

M

S

6.5

西

1

0

0

k

m

巧家地震

震源机制解

M

S

5.0—5.9

M

S

6.0—6.9

M

S

7.0—7.9

26°

图1 2020年巧家

M

S

5.0地震区域历史地震分布

Fig.1 Historical earthquakes in the study region of 2020 Qiaojia

M

S

5.0 earthquake

第5期

卢 显等：2020年5月18日云南巧家

M

S

5.0地震总结155

1900年以来，巧家地震震中100 km

103.1°103.2°103.3°103.4°103.5°E

27.4°N

范围内发生

M

S

≥

5.0地震27次，其

中5.0—5.9级地震25次，6.0—6.9

27.3°

级地震2次，距今时间最短、震级最

大、空间距离最近的地震是2014年8

27.2°

月3日云南省昭通市鲁甸县

M

S

6.5地

27.1°

震，与此次巧家5.0级地震相距约

20 km（图1）。2014年鲁甸

M

S

6.5地震

27.0°

即发生在大凉山块体内部，莲峰断裂与西

鱼河—昭通断裂之间（图1）。房立华等

断裂

26.9°

（2014）利用双差定位方法，根据震中距

巧家地震

鲁甸地震序列

200 km范围内29个地震台站记录的地震

26.8°

波形资料，对鲁甸地震序列进行精定位，

结果显示，余震分布呈SE、SW两个优势

图2 鲁甸

M

S

6.5地震序列与巧家

M

S

5.0地震分布

Fig.2 Earthquake sequence of Ludian

M

S

6.5

方向，具有不对称共轭特征。2个余震条

earthquake and the Qiaojia

M

S

5.0 earthquake

带夹角约100°，展布长度均约16 km。鲁

甸地震序列早期余震分布主要沿NW—SE向垂直于昭通—鲁甸断裂，而此次巧家5.0

级地震距其发震断层及余震分布集中区较远，分析认为，此次巧家地震不是2014年鲁甸

（图2）。

M

S

6.5地震的余震

蓬

峰

断

裂

2 巧家

M

S

5.0地震前地震活动异常

2020年巧家

M

S

5.0地震发生前，震中附近先后出现巧家震群、Benioff应变比和地

震发生率异常，这些异常在震前9个月至震前3个月依次出现，异常空间位置集中在此

次地震震中周边，预测的震级水平多为5级。

2.1 巧家震群

2020年巧家5.0级地震前的突出现象是巧家震群的发生，该震群自2020年1月17

日开始活动，至5月7日共发生

M

L

≥

1.0地震82次，其中最大地震为3月18日

M

L

3.9

地震，次大地震为2月7日

M

L

3.8地震，最大日频次12次（图3），符合震群定义（国家

地震局预测预防司，1997）。进一步计算震群参数

U

、

F

、

K

、

ρ

、

b

、

h

，其中

U

= 0.601 7，

F

= 1.592 6，

K

= 0.684 2，

ρ

= 0.435 8，

b

= 0.570 1，

h

= 0.500 0。按照前兆震群的判

定标准（国家地震局预测预防司，1997），以上6个参数中有4个显示为前兆震群，因

此判定该震群为前兆震群（表1），对后

表1 巧家

M

S

5.0震群前兆震群性质判定指标

续震群周边中强地震有预测意义（黄耘

Table 1 Determination index of precursory swarm

properties of Qiaojia

M

S

5.0 earthquake swarm

等，2011；陈玉华等，2013；曹凤娟等，

参数名称参数值 判定标准判定结果

2015）。此次巧家5.0级地震震中距震群

0.601 7 ＞0.50前兆震群

U

值

约12 km，且发生在震群最大日频次后

1.592 6＞0.70前兆震群

F

值

0.435 8 ＜0.55前兆震群

ρ

值

101天，结合前兆震群计算参数，分析

0.684 2

≤

0.70非前兆震群

k

值

认为，此次巧家震群为2020年5月18

0.500 0 ＜1.00前兆震群

h

值

0.570 1

≤

0.65非前兆震群

b

值

日巧家地震前的短期异常。

西

鱼

河

—

昭

通

断

裂

156

4.0

3.0

震

级

M

L

2.0

1.0

01-10

15

(b)

10

频

度

01-30

地震地磁观测与研究41卷

(a)

02-1903-1003-3004-1905-09

5

0

01-1001-30

02-1903-10

t/月-日

03-3004-1905-09

图3 巧家震群

M—t

及频次图

Fig.3

M

-

t

and frequency figures of Qiaojia earthquake swarm

2.2 Benioff应变比

Benioff应变比是区域内某一特定日期后3个月内地震能量的开方及与该日期前相同

时间窗长地震能量的开方和之比（杨文等，2014）。应变比的高低反映了区域内滑动速率

的快慢。对我国西部Benioff应变比进行空间扫描，结果显示，2020年1月，在此次巧

家地震震中附近存在显著高值异常[图4(a)]。据以往震例分析该高值区域附近1年内5

级地震发生的可能，发现2020年5月18日巧家5.0级地震即发生在异常区内部，而在

此次地震发生后，震中附近的Benioff应变比高值异常消失[图4(b)]。因此，认为此次

Benioff应变比异常应为巧家5.0级地震映震异常。

92°96°100°104°108°E

92°96°100°104°

108°E

R

0.8

0.6

0.4

28°28°

0.2

0.0

24°24°

-0.2

-0.4

20°

2019-08-01—2020-01-31

(a)

2019-12-01—2020-05-31

(b)

20°

-0.6

-0.8

32°N32°N

(a)2019年8月1日至2020年1月31日阶段应变比；(b)2019年12月1月至2020年5月31日阶段应变比

图4 巧家

M

S

5.0地震前后Benioff应变比异常

Fig.4 Benioff strain ratio anomalies before and after Qiaojia

M

S

5.0 earthquake

第5期

卢 显等：2020年5月18日云南巧家

M

S

5.0地震总结157

2.3 地震发生率

地震发生率异常是基于泊松分布得到的实际地震发生率相对于背景地震发生率的异

常偏离程度，参数值接近1反映地震活动显著增强，接近0反映地震活动显著减弱。利

用G

-

C法（陈凌等，1998），将云南巧家

M

S

5.0地震删除余震后

M

L

3.0以上地震目录作

为计算资料，以5年为背景窗长、3个月为计算窗长，进行地震发生率异常时空扫描，

结果见图5。由图5可见，2020年5月18日云南巧家

M

S

5.0地震前3个月，在震中附

近出现地震发生率高值异常，震前2个月异常面积进一步扩大，震前1个月异常明显减

弱随后消失，之后在原异常区域发生巧家

M

S

5.0地震。

98°100°102°104°106°108°E98°100°102°104°106°108°E98°100°102°104°106°108°E

32°N32°N32°N

(a)(b)(c)

30°

2020-05-18

5.0

30°

2020-05-18

5.0

30°

2020-05-18

5.0

28°

26°

24°

22°

28°

26°

24°

22°

28°

26°

24°

22°

1

0.9

0.8

0.7

0.6

0.5

0.4

0.3

0.2

0.1

图5 巧家5.0级地震前地震发生率异常演化

Fig.5 Anomaly evolution of earthquake rate before the Qiaojia

M

S

5.0 earthquake

(a) 震前3个月异常分布；(b)震前2个月异常分布；(c)震前1个月异常分布

3 地球物理观测异常

30°N

巧家地震震中200 km范围内地球物理

形变

电磁

观测台站（国家测震台网）主要布设钻孔应

流体

蓝色：趋势

变、垂直摆、水平摆、重力、地磁、地电

绿色：短期

灰色：背景

阻率、水位、水温、断层氢等测项，共279

28°

个观测台项（图6），除个别台站外，大多数

台站的观测资料连续、可靠。

200 km

100 km

此次地震发生前，震中200 km范围内

26°

地球物理异常测项15项（表2），其中：震

中100 km范围内（测项119项）异常测项7

项，异常测项比（异常测项数/观测测项数）

24°

较低（5%）；震中101—200 km范围内（测

图6 云南巧家

M

S

5.0地震震中300 km范围内

地球物理观测台及异常测项分布

项161项）异常测项8项，异常测项比为

Fig.6 Fixed geophysical observation stations and

8/161（5%）。在15项异常中有14项为趋势

abnormal distribution of Qiaojia

M

S

5.0 earthquake

异常，均非此次巧家5.0级地震异常，而云

南省昭通市地震局地球物理观测网的巧家电磁波NS、EW分量作为短期异常，应为该地

震的映震异常。中国地震台网中心形变研究室研究认为，位于鲜水河断裂带的侏倭跨断

层测距短期异常应为此次地震映震异常。

100°102°104°E106°

158地震地磁观测与研究41卷

表2 云南巧家

M

S

5.0地震震中200 km范围内地球物理测项异常

Table 2 Fixed geophysical anomalies within 200 km of epicenter of Qiaojia

M

S

5.0 earthquake

台站

汤家坪

汤家坪

宁南

尔乌

金河

金河

金河

金河

测项

(A←B)水准

(A←C)水准

(5←4)水准

(3←A)水准

YRY

-

4钻孔线应变NE

YRY

-

4钻孔线应变EW

YRY

-

4钻孔线应变NW

YRY

-

4钻孔线应变NS

震中距/km异常性质

43

43

46

111

137

137

137

137

趋势

趋势

趋势

趋势

趋势

趋势

趋势

趋势

台站

西昌

南山

平地

巧家

昭通茨院

大关

会理31井

（原06井）

测项

地磁GM4—

Z

地磁GM4—

Z

地磁GM4—

Z

电磁波

水位

水位

静水位

震中距/km

99

160

170

待填

46

95

147

异常性质

趋势

趋势

趋势

短期

趋势

趋势

趋势

3.1 电磁学科异常分析

震中200 km范围内共有3个电磁异常测项，均为趋势异常，分别为西昌、南山和

平地的地磁异常（表2）。此次巧家地震震级较低，未达到3个地磁异常测项预测的目标

震级，非映震异常。云南省昭通市地震局地球物理观测网的巧家电磁波NS和EW分量

（图7）于2020年5月6日同步出现扰动幅度增大现象，5月18日发生巧家

M

S

5.0地震。

该台电磁波观测资料在2019年下半年出现较为突出的异常信息，表现在2019年6月17

日长宁

M

S

6.0地震序列（震中距约250 km）及之后几次5级以上地震发生前，电磁曲线均

出现明显的高值异常现象，故此次5月6日出现的电磁扰动幅度增大异常应为巧家

M

S

5.0

地震映震异常。

磁

场

N

S

分

量

/

n

T

6185

兴文6.0

珙县5.4

珙县5.6

内江5.4

资中5.2

巧家5.0

6183

(a)

6181

2019-06

5420

兴文6.0

珙县5.4

珙县5.6

070809

内江5.4

101112205

巧家5.0

磁

场

N

S

分

量

/

n

T

资中5.2

5418

(b)

5416

2019-06

t/年-月

205

图7 巧家电磁波日均值曲线

Fig.7 The daily average curve of Qiaojia electromagnetic wave

3.2 形变学科异常分析

（1）8项形变异常。震中200 km范围内共分布8个形变异常测项，均为趋势异常，

且均非此次巧家5.0级地震映震异常。此次地震发生前还存在以下趋势异常：汤家坪水

准测线A—B/A—C、宁南水准测线5—4和尔乌水准测线3—A记录数据异常，且震后

异常仍持续，尚不能判定是否为此次

M

S

5.0地震异常，需继续跟踪。定点形变趋势异常

测项为金河钻孔应变NE、EW、NW、NS分量，异常特征表现为：2019年6月15日起，

第5期

卢 显等：2020年5月18日云南巧家

M

S

5.0地震总结159

4分量呈快速上升的张性变化，且幅度较大，12月1日呈转折下降的压性变化，2020

年6月1日后应变曲线再次转折上升，而此次巧家5.0级震中周边300 km范围内曾发

生13次

M

S

5.0以上地震，震中距最近的一次地震为2018年西昌

M

S

5.1地震，震中距为

19 km，以上震前均无此类应变异常出现，故认为该异常信度不高，不作为本次巧家5.0

级地震的震前异常。

（2）侏倭跨断层异常。四川省地震局测绘工程院于2018年初在鲜水河断裂北段侏倭

新建跨断层测量场地，该测点距此次巧家地震震中548 km，虽然间距较远，但通过对以

往震例的研究，认为地震往往发生在异常区外围，而非异常区内。中国地震台网中心形

变预报研究室在此次巧家5.0级地震发生前，根据侏倭跨断层测距拉张异常，提出川滇

交界东部有发生5.5级左右地震的危险性，并于2020年4月3日提出短期预测意见，且

在之后的周月地震会商报告中坚持该预测意见，直至巧家地震发生。

侏倭场地每月测量1期，测线1—3和2—3的测距在2018年全年变化不大，呈

平稳趋势，2019年1月，2条测线的测距曲线突然分别快速上升（拉张）2.80 mm和

2.73 mm，2月11日达峰值（图8），后转折下降（收缩）。其间，于6月17日发生长宁

M

S

6.0

地震，并于珙县等地发生多次5级以上地震。2020年2月，测线1—3和2—3的测距

曲线再次同步上升，4月10日达峰值，后转折下降，于4月1日发生石渠

M

S

5.6地震，

5月18日发生巧家

M

S

5.0地震。2组震例的共同特征是曲线上升（拉张）—转折下降（收

缩）—发震，等待时间约2—6个月。因此，侏倭跨断层测距的转折下降变化应为巧家5.0

级地震的震前异常。

250902.0

250901.5

测

距

/

m

m

250901.0

250900.5

250900.0

250899.5

250899.0

250898.5

2018-03-052019-03-05

t/年-月-日

2020-03-052018-03-05

2019-03-05

t/年-月-日

2020-03-05

珙县5.6

兴文5.7

(a)

威远5.4

长宁6.0

石渠5.6

巧家5.0

(b)

162082

测

距

/

m

m

兴文5.7

长宁6.0

威远5.4

石渠5.6

巧家5.0

162081

162080

162079

珙县5.6

图8 侏倭跨断层测距曲线与周边600 km范围内5级以上地震对应

Fig.8 The corresponding situation of Zhuwo cross fault curve and

M

S

≥

5.0 earthquakes within 600 km of Zhuwo

(a)测线1—3；(b)测线2—3

3.3 流体学科异常分析

震中200 km范围内，流体学科测项在震前存在3项异常，均为趋势背景异常，分

别为昭通茨院水位、大关水位和会理水位测项，其中茨院水位2020年1月出现破年变高

值异常，大关水位2019年9月出现破年变低值异常，会理水位2019年3月出现快速上

升后年变幅增大异常。此次巧家5.0级地震发生前，未发现3项异常出现明显短临变化。

昭通茨院水位、大关水位主要对应川滇交界东部地震，据以往震例，目前2项水位

变化异常非此次巧家5.0级地震异常；会理水位异常此前无震例对应，且目前该异常并

160地震地磁观测与研究41卷

未恢复，仍需跟踪分析，不能判定是否为此次5.0级地震异常。

4 综合计算方法

根据王海涛等（2008）提出的基于多种地震学参数的地震对应相关谱识别地震前兆异

常的方法，截取2020年5月9日前巧家地震序列记录，采用多参数概率谱方法进行计

算，结果显示，巧家地震发生在预测的危险区域和时间范围内。多参数概率谱方法主要

选择以下7个参数进行综合计算：

①

描述地震活动震级和频度关系的

b

值、缺震、

η

值；

②

描述地震危险性的

D

值；

③

描述地震活动自身强度变化的

M

f

值、频度

N

值；

④

环境

因子调制参数

R

m

值。

2019年6月17日四川长宁

M

S

6.0地震（104.90°E，28.34°N）发生后，为了进行跟

踪分析，以震中为中心，提取半径200 km范围内2006年1月至2018年12月发生的

M

L

5.5

地震，利用震前小震活动的多种地震学参数综合异常特征，外推预测2019年1月以来的

震情形势。时间跟踪分析表明：目标地震在震前均出现多参数综合高值异常；研究区域

在2016—2017年未发生目标地震，多参数概率谱曲线出现高值虚报，见图9虚框内所

示地震。研究发现，2016—2017年，该区域多次发生

M

L

5.0以上地震，因此容易出现

不太突出的多参数概率谱高值现象。

6.0

5.9

5.6

5.5

6.8

5.3

5.3

5.3

5.15.2

5.65.9

5.3

5.8

6.0

5.5

6.3

巧家地震

5.7

5.6

0.4

综

合

概

率

指

数

0.3

0.2

0.1

2017-12

2010-122013-12

t/年-月

2016-122019-12

图9 多参数概率谱变化曲线

Fig.9 The change of earthquake corresponding probability spectrum of multi-parameter

2019年10月多参数数值结果开始转折上升，云南巧家地震发生时间位于综合指数

上升阶段，震中位于研究区域边缘，综合分析认为，2019年10月至2020年4月，在长

宁地震震中200 km范围内出现的多参数概率谱异常，应为2020年巧家5.0级地震的震

前异常。

5 巧家

M

S

5.0地震震源参数与序列跟踪

5.1 震源物理参数

2020年巧家

M

S

5.0地震发生后，不同研究机构利用不同方法给出该地震的震源机

制解，结果见表3，可见各机构所得震源机制解较为一致，滑动角反演结果均显示此次

地震为走滑型地震。地震序列余震空间分布显示，余震呈近NS向展布，认为其发震断

第5期

卢 显等：2020年5月18日云南巧家

M

S

5.0地震总结161

层应为位于莲峰断裂与昭通—鲁甸断裂之间的NS向断裂。结合余震空间展布，选取近

NS向节面作为破裂面，即选取各机构反演所得节面

Ⅰ

为地震破裂面。各机构反演所得

断层面均较陡，倾角在60°—70°，且滑动角在

-

6°—34°，表明地震破裂以左旋走滑为

主兼正断分量。

表3 不同机构给出的云南巧家

M

S

5.0地震震源机制解

Table 3 Different focal mechanism solutions by different institutions

研究机构

中国地震台网中心

中国地震局地球物理研究所

中国地震局地震预测研究所

走向/(°)

170

180

174

节面

Ⅰ

倾角/(°)

70

60

72

滑动角/(°)

-

10

-

24

-

6

走向/(°)

263

283

266

节面

Ⅱ

倾角/(°)

81

69

84

滑动角/(°)

-

160

-

148

-

162

5.2 余震序列跟踪

表4 巧家

M

S

5.0地震序列目录（

M

L

≥

3.0）

截至2020年6月30日，巧家

M

S

5.0地

Table 4 The catalog of Qiaojia

M

S

5.0

震余震区共记录

M

L

≥

0.0地震780次，其

earthquake sequence (

M

L

≥

3.0)

中

M

L

0.0—0.9地震485次，

M

L

1.0—1.9

编

震中位置/(°)

发震日期震级深度/km

φ

N

λ

E

地震238次，

M

L

2.0—2.9地震53次，

M

L

3.0—

号

8

12020

-

05

-

1827.20 103.17

M

S

5.0

3.9地震4次，最大余震为5月19日

M

L

3.8

722020

-

05

-

1827.18 103.13

M

L

3.6

1032020

-

05

-

1927.18 103.15

M

S

3.5

（

M

S

3.5）地震（表4，图10）。

542020

-

05

-

2027.18 103.15

M

L

3.2

由

M

—

t

及日频次图（图10）可见，余震

3.0952020

-

06

-

0127.18 103.13

M

L

集中发生在震后5天内，共发生

M

L

≥

0.0

地震339次，其中5月19日和5月20日发生余震相对较多，共197个

M

L

≥

0.0地震。

5月22日以后，余震日频次持续衰减，但衰减速度较慢，共发生

M

L

≥

0.0地震441次，

占现阶段记录余震总数的56.5%。

6

5

震

级

4

3

2

1

0

100

80

地

震

数

60

40

20

0

05-1805-2205-2605-3006-0306-0706-11

t/年-月-日

06-1506-1906-2306-2707-01

(b)

(a)

图10 巧家

M

S

5.0地震序列

M

—

t

及频次图

Fig.10

M

-

t

diagram and the daily cumulative rate of Qiaojia

M

S

5.0 earthquake sequence

162地震地磁观测与研究41卷

巧家地震序列主震与最大余震的震级差为1.6（

M

L

震级），主震释放能量占序列总能

量的99.64%，属主余型地震序列。从空间分布看，余震活动集中在主震附近13 km范围内，

近似NS向分布，长轴约23 km，短轴约13 km。最大余震位于主震西南方，与主震较近，

间距约3 km（图11）。

基于巧家地震发生前后余震区的地震目录，利用ZMAP软件中最小完整性震级

M

c

的计算方法（Wyss et al，1999；Wiemer et al，2000；Mignan et al，2012），即结合拟合优

度测试（the Goodness-of-Fit Test，简写为GFT）和修正最大曲率法（简称为修正MAXC），

计算得到巧家

M

S

5.0地震前后余震区的

M

c

值分别为

M

L

1.7 ± 0.2和

M

L

0.5 ± 0.3（图

12），震前

M

c

值大于震后

M

c

值。此结果与一般中强地震后，由于小震级余震缺失造成

的

M

c

值变大现象不相符。可能是因为，巧家

M

S

5.0地震发生后，云南省地震局于5月

19日在震中附近架设3个流动台站，增强了余震区地震监测能力，可完整记录

M

L

≥

1.0

地震（党文杰等，2020）。同时，由于GFT和修正MAXC均存在对于实际最小完整性震

级的低估现象（Mignan et al，2012），故选定

M

L

1.0为余震序列的最小完整性震级。

102.8°103.0°103.2°103.4°103.6°E

7

(a)

6

最

小

完

整

性

震

级

M

L

5

4

3

2

1

M

c

=1.7±0.1

2012-10-162016-08-01

最小完整性震级M

L

最小完整性震级不确定M

L

实测地震

7

6

5

4

3

2

1

M

L

M

L

金

阳

断

裂

27.6°N

27.4°

2008-12-31

27.2°

最

小

完

整

性

震

级

M

L

断

2020-05-18

7

6

5

4

3

2

1

7

(b)

6

5

4

3

2

1

M

c

= 0.6 ± 0.3

2012-10-162016-08-01

t/年-月-日

基于上述最小完整性震级，利用最大似然法（Marzocchi et al，2003），对巧家

M

S

5.0

地震序列震级—频度曲线进行G—R关系拟合，所得序列

a

值和

b

值分别为3.15 ± 0.08

和0.68 ± 0.03[图13(a)]，进而推测序列最大余震震级为

M

L

4.6，大于实际发生的最大

余震震级

M

L

3.8。同时，由图13(a)可见，在震级

M

L

2.5—4.6范围内，利用G—R关系

拟合曲线给出的地震累积数明显较多，表明地震存在缺失现象，未来可能发生该震级范

围的余震。此外，为刻画地震序列的衰减特征，计算巧家序列的

h

值。据刘正荣等（刘

正荣等，1979，1986；刘正荣，1984）、蒋海昆等（2007）的研究，若

h

＞ 1，地震序列一般为

则

木

河

裂

昭

通

甸

鲁

—

断

裂

27.0°

图11 巧家

M

S

5.0地震序列震中分布

Fig.11 The epicenter distribution of Qiaojia

M

S

5.0 earthquake sequence

北

段

断

裂

小

江

26.8°

2008-12-312020-05-18

图12 巧家

M

S

5.0地震震前(a)和震后

(b)余震区最小完整性震级

Fig.12 The complete magnitude of the aftershock area

before (a) and after (b) Qiaojia

M

S

5.0 earthquake

第5期

卢 显等：2020年5月18日云南巧家

M

S

5.0地震总结163

主余型地震序列，若

h

≤

1，地震序列一般为前震序列，后续会发生更大地震。截至6月30日，

巧家

M

S

5.0地震序列

h

值为1.497，大于1，表明该序列属主余型地震序列[图13(b)]。

10

3

(a)

18

实际累积频度—震级分布

拟合G—R关系

最小完整性震级：M

1.0

估计最大余震震级：M

4.6

a = 3.15 ± 0.08; b = 0 .68 ± 0.03

M

= M

1.0

dt = 24 h

h = 1.497

h 值参考线

拟合曲线

实测数据

(b)

0.30.40.50.6

16

14

12

N

R

(

t

)

10

8

0.7

10

2

地

震

数

2015-05-18—07-01

0.8

0.9

1.0

1.1

1.2

1.3

1.4

1.6

2.0

3.0

4.0

10

1

6

4

2

10

0

00.511.522.5

M

L

33.544.55

00

t/d

图13 巧家

M

S

5.0地震序列

b

值(a)、

h

值(b)计算结果

Fig.13 The

b

-value (a) and

h

-value (b) of Qiaojia

M

S

5.0 earthquake sequence

6 讨论

6.1 地震活动和地球物理观测异常

（1）地震活动。2020年巧家

M

S

5.0地震的前突出地震活动异常为巧家震群，该震群

自2020年1月17日开始活动，分析判定为巧家地震序列的前兆震群。2019年8月至

2020年1月的应变比显示川滇交界地区存在高值异常，2020年5月18日云南巧家

M

S

5.0

地震即发生在异常区域，地震后异常逐渐消失。此外，震前3个月在震中地区出现地震

发生率高值异常，震前2个月异常面积进一步扩大，震前1个月异常明显减弱随后消失，

之后在原异常区域发生巧家

M

S

5.0地震。根据地震序列动态跟踪过程，巧家

M

S

5.0地震

序列属主余型地震序列，主震震源机制为走滑型。

（2）地球物理观测。震中200 km范围内共布设279个地球物理观测台项，其中异常

测项15项，异常测项比（异常测项数/观测测项数）为5%。15项异常多为趋势异常，短

期异常不突出，而巧家电磁波震前出现扰动幅度增大的异常，属短期异常，应为本次地

震异常；形变预报研究室提出，侏倭跨断层在震前出现测距曲线上升后转折下降的现象，

并于下降过程中对应巧家地震。利用综合方法中多参数概率谱进行计算，在长宁余震区

的跟踪研究中，出现曲线上升异常现象，本次巧家地震发震时间位于综合指数上升阶段，

震中位于跟踪区域边缘，应为震前短期异常。

6.2 震中附近应力变化

强震活动通常发生在活动断裂带上具有高应力积累的凹凸体（Aki，1984）或者闭锁

段（Wyss et al，2000）。

b

值与应力水平成反比，

b

值越低，应力水平越高（Urbancic et al，

1992）。巧家以北至东川以南段位于小江断裂带最北段，1900—1966年5级地震频发，

1966年2月5日东川

M

S

6.5地震发生后处于5级以上地震平静状态，直至2005年8月5

164地震地磁观测与研究41卷

日于该段南部发生会泽

M

S

5.3地震。根据M7专项工作组（2012）的研究结果，自1733年

东川7级地震发生后，巧家以北至东川以南段处于7级强震空区，且该段大地震破裂的原

地平均复发间隔一般为100—200年。毛燕等（2016）计算发现，小江断裂带最北段

b

值在

中偏低的水平，表明此段当前应力水平较高，中小地震频发，未来发生强震的可能性较大。

地震发生时，断层破裂引起的同震位错会导致破裂面周边应力场变化，基于弹性位错理

论，计算所得同震应力场为二阶张量，投影到目标断层面上得到其正应力和剪应力变化，最

终得到地震所导致的库仑破裂应力变化（Barka，1999；Shan et al，2009；万永革等，2010）。

付芮等（2015）利用弹性位错理论和分层岩石圈模型，计算2014年鲁甸

M

S

6.5地震导致的同

震库仑应力变化场，讨论主震对余震分布及周边主要活动断层的影响（图14，图中红色五角

星为此次巧家地震），认为鲁甸地震加大了则木河断裂和小江断裂交界处的应力积累。

103°E103.5°E104°E

28°N

1.00

0.10

0.05

昭

通

断

裂

0.003

0.001

-0.001

-0.003

-0.01

27°

-0.03

-0.05

昭通

鲁甸

巧家

断

小

江

裂

-0.10

-1.00

图14 鲁甸

M

S

6.5地震对周边主要活动断层应力积累的影响（据付芮等，2015）

Fig.14 The influence of stress accumulation of Ludian

M

S

6.5 earthquake on main

active faults nearby (according to Fu Rui et al, 2015)

闻学泽等（2008）通过对古地震及鲜水河—小江断裂带应力转移计算结果的研究指

出，在则木河断裂和小江断裂处的巧家附近，存在一个直径20 km以上、应力积累较

强的地震空区，具有发生中强以上地震的风险，而2014年鲁甸

M

S

6.5地震的发生，进

一步增强了该区域今后中强以上地震的危险性。另外，鲁甸地震提升了莲峰断裂NE段

0.003 MPa以下应力的积累（Shan et al，2013），虽然小于应力触发阈值（0.01 MPa）（Harris

et al，1998），应力增量变化不明显，但作为现今的地震空区之一，当时认为需关注该段

未来发生地震的危险性。2020年巧家

M

S

5.0地震恰发生在则木河断裂与小江断裂处应力

积累较强的地震空区，反映了该区域应力水平的持续累积状态。

6.3 三维地壳速度结构

区域地壳速度结构成像结果显示，川滇地块的下地壳自北向南存在2条低速带，一

应

力

变

化

/

M

P

a

则

木

河

断

裂

莲

峰

断

裂

27.5°

0.03

0.01

第5期

卢 显等：2020年5月18日云南巧家

M

S

5.0地震总结165

条向西南延伸至丽江—小金河断裂和红河断裂，另一条延伸至东南部，位于小江断裂带

下方。这2条地壳通道流和峨眉山大火成岩省刚性块体，通过走滑断层的挤压作用，在

青藏高原东南缘晚中新世至今的区域动力学演变中起着重要作用（Bao et al，2015；Yang

et al，2020）。2020年巧家

M

S

5.0地震发生在高速异常的峨眉山大火成岩省中带，震中

10 km深度为高速体和低速体交界区域，地震活动频繁（徐涛等，2015）。

李大虎等（2019）应用区域地震和远震联合反演的方法，得到川滇地区三维地壳速度

结构图像，并重点剖析了昭通、莲峰断裂带P波速度结构（图15，图中黑线为区内主要

断裂，白色五角星为2014年鲁甸

M

S

6.5地震，红色五角星为2020年巧家

M

S

5.0地震），

研究结果表明，在10 km深度层上，昭通、莲峰断裂带及周边的上地壳速度结构呈现明

显的横向不均匀分布，昭通、莲峰断裂带处于高波速异常区，其间存在局部相对低速异

常体。2020年巧家

M

S

5.0地震和2014年鲁甸

M

S

6.5地震均发生在该高低速异常体分界

线附近，略偏向高速体一侧，随着深度增加，巧家和鲁甸地震附近中下地壳速度异常分

布特征出现一定趋势性变化，即存在大范围低速异常分布，这种中下地壳物质相对较为

软弱的结构特征有利于应力在脆性的上地壳内积累和集中，加之横向挤压的构造应力场

作用，易发生破裂继而引发地震，从而对巧家

M

S

5.0地震的孕育和发生起到重要作用。

102°

MG

JY

ZT

LD

MY

HD

QJ

HD

深度=1 km

MY

HD

QJ

HD

深度=10 km

104°E102°

MG

104°E102°

MG

104°E

断裂

YJ

YL

ZT

LD

YJ

YL

28°N

YJ

YL

ZT

LD

JY

28°N

JY

28°N

巧家地震

鲁甸地震

MY

HD

QJ

HD

深度=15 km

6

26°

4

2

0

28°N

-2

-4

-6

深度=50 km

102°

MG

JY

104°E

YJ

YL

ZT

LD

26°

102°

MG

104°E

YJ

YL

ZT

LD

26°

102°

MG

104°E

YJ

YL

ZT

LD

28°N

JY

28°N

JY

MY

HD

QJ

HD

深度=30 km

MY

HD

QJ

HD

深度=40 km

MY

HD

QJ

HD

26°26°26°

图15 三维P波速度异常分布（据李大虎等， 2019）

Fig.15 The images of 3D P wave velocity perturbation (according to Li Dahu et al, 2019)

7 结论

全面梳理、系统总结每一次显著地震的地震地质环境和震源参数及地震序列跟踪过

程，尤其是震前出现并发现的地震活动及地球物理观测异常变化，对地震预测的科学研究

具有重要意义。2020年巧家

M

S

5.0地震震中与2020年度“四川道孚至云南永善7.0级左右

年度地震重点危险区”相距约51 km，同时位于2020年5—7月全国地震概率预测图所示

西南地区的相对高概率区附近。通过对巧家

M

S

5.0地震震例的系统总结，得出以下结论。

d

v

/

v

(

%

)

166地震地磁观测与研究41卷

（1）近30年来，巧家

M

S

5.0地震震中附近发生多次5级以上地震，其中2014年8

月3日鲁甸

M

S

6.5地震距离最近，约20 km，但此次地震不是其余震。对巧家地震余震

序列跟踪分析发现，主震与最大余震的震级差较大，且主震释放能量较高，判定该序列

属主余型地震序列。

（2）巧家

M

S

5.0地震发生在莲峰断裂和昭通—鲁甸断裂之间，余震序列呈NS向展

布。多家研究机构利用不同方法给出此次5.0级主震的震源机制解，结果较为一致，且

滑动角反演结果均显示此为走滑型地震。

（3）通过对此次巧家

M

S

5.0地震前异常的梳理，发现映震异常有巧家震群、Benioff

应变比、地震发生率、巧家电磁波、侏倭跨断层及多参数概率谱，且以中期异常为主，

而巧家电磁波和侏倭跨断层可作为此次地震的短期异常。

本文撰写得到王海涛研究员的指导和鼓励，蒋海昆研究员、晏锐研究员、孟令媛研

究员和闫伟高级工程师亦给予帮助，中国地震台网中心国家地震科学数据中心（

）提供数据支撑，在此对他们及中国地震台网中心预报部同事的辛苦工

作，一并表示衷心感谢。

参考文献

曹凤娟，翟丽娜，殷轶娜，等. 辽宁地区震群序列统计特征及预测意义研究[J]. 中国地震，2015，31（1）：45

-

55.

陈玉华，姚家骏，李戈云，等. 1980年以来青海地区震群活动与前震序列特征分析[J]. 中国地震，2013，29（4）：489

-

500.

程佳，刘杰，徐锡伟，甘卫军. 大凉山次级块体内强震发生的构造特征与2014年鲁甸6.5级地震对周边断层的影响[J]. 地震地质，

2014，36（4）：1 228

-

1 243.

（3）：

党文杰，李文涛，曹舸斌，等. 应急流动监测台网在云南巧家

M

S

5.0地震中的应用[J]. 地震地磁观测与研究，2020，41

193

-

199.

（4）：1 173

-

1 185.

房立华，吴建平，王未来，等. 云南鲁甸

M

S

6.5地震余震重定位及其发震构造[J]. 地震地质，2014，36

付芮，单斌，熊熊，等. 2014年云南鲁甸地震同震库仑应力对余震分布及周边断层的影响[J]. 地震地质，2015，37（4）：

1 084

-

1 095.

国家地震局预测预防司. 测震学分析预报方法[M]. 北京：地震出版社，1997：105

-

174.

黄耘，孙业君，杨云，叶碧文. 江苏及邻区震群活动特征分析[J]. 中国地震，2011，27（1）：72

-

82.

蒋海昆，傅征祥，刘杰，等. 中国大陆地震序列研究[M]. 北京：地震出版社，2007：38

-

44.

李大虎，丁志峰，吴萍萍，等. 川滇交界东段昭通、莲峰断裂带的深部结构特征与2014年鲁甸

M

S

6.5地震[J]. 地球物理学报，

2019，62（12）：4 571

-

4 587.

2020，41（3）：

李文涛，党文杰，金明培，等. 2020年云南巧家

M

S

5.0地震震源机制解与发震断层归属[J]. 地震地磁观测与研究，

184

-

192.

刘正荣，钱兆霞，王维清. 前震的一个标志—地震频度的衰减[J]. 地震研究，1979，（4）：1

-

9.

刘正荣. 根据地震频度衰减预报地震的工作细则[J]. 地震，1984，（1）：35

-

37.

刘正荣，孔绍麟. 地震频度衰减与地震预报[J]. 地震研究，1986，9（1）：1

-

10.

毛燕，刘自凤，叶建庆，等. 小江断裂带强震危险性分析[J]. 地震研究，2016，39（2）：213

-

217.

M7专项工作组. 中国大陆大地震中

-

长期危险性研究[M]. 北京：地震出版社，2012：35

-

36.

陈凌，刘杰，陈颙，等. 地震活动性分析中余震的删除[J]. 地球物理学报，1998，41（Z1）：244

-

252.

万永革，沈正康，盛书中，等. 2008年新疆于田7.3级地震对周围断层的影响及其正断层机制的区域构造解释[J]. 地球物理学报，

2010，53（2）：280

-

289.

王海涛，王琼. 基于地震对应概率谱分析的前兆异常识别研究[J]. 地震研究，2008，31（4）：330

-

334.

闻学泽，范军，易桂喜，等. 川西安宁河断裂上的地震空区[J]. 中国科学：地球科学，2008，38（7）：797

-

807.

闻学泽，杜方，易桂喜，等. 川滇交界东段昭通、莲峰断裂带的地震危险背景[J]. 地球物理学报，2013，56（10）：3 361

-

3 372.

徐涛，张忠杰，刘宝峰，等. 峨眉山大火成岩省地壳速度结构与古地幔柱活动遗迹：来自丽江

-

清镇宽角地震资料的约束[J].

中国科学：地球科学，2015，45（5）：561

-

576.

杨文，周龙泉，刘杰，等. 新疆于田两次7.3级地震之后南北地震带Benioff应变比分布特征研究[J]. 中国地震，2014，30（2）：

219

-

225.

Aki K. Asperities, barriers, characteristic earthquakes and strong motion prediction[J]. Journal of Geophysical Research: Solid

Earth, 1984, 89(B7): 5 867

-

5 872.

第5期

卢 显等：2020年5月18日云南巧家

M

S

5.0地震总结167

Bao X W, Sun X X, Xu M J, et al. Two crustal low

-

velocity channels beneath SE Tibet revealed by joint inversion of Rayleigh wave

dispersion and receiver functions[J]. Earth and Planetary Science Letters, 2015, 415: 16

-

24.

Barka A. The 17 August 1999 Izmit earthquake[J]. Science, 1999, 285(5 435): 1 858

-

1 859.

Harris R A, Simpson R W. Suppression of large earthquakes by stress shadows: a comparison of coulomb and rate

-

and

-

state

failure[J]. Journal of Geophysical Research: Solid Earth, 1998, 103(B10): 24 439

-

24 451.

Marzocchi W, Sandri L. A review and new insights on the estimation of the

b

-

value and its uncertainty[J]. Annals of Geophysics,

2003, 46(6): 1 271

-

1 282.

Mignan A, Woessner J. Estimating the magnitude of completeness for earthquake catalogs[M/OL]. Community Online Resource for

Statistical Seismicity Analysis, 2012 [2020

-

08

-

21]. /export/sites/corssa/.galleries/articles

-

pdf/Mignan

-

Woessner

-

2012

-

CORSSA

-

Magnitude

-

of

-

.

Shan B, Xiong X, Zheng Y, et al. Stress changes on major faults caused by

M

W

7.9 Wenchuan earthquake, May 12, 2008[J]. Science

in China Series D：Earth Sciences, 2009, 52(5): 593

-

601.

Shan B, Xiong X, Zheng Y, et al. Stress changes on major faults caused by 2013 Lushan earthquake, and its relationship with 2008

Wenchuan earthquake[J]. Science China Earth Sciences, 2013, 56(7): 1 169

-

1 176.

Urbancic T I, Trifu C I, Long J M, et al. Space

-

time correlations of

b

values with stress release[J]. Pure and Applied Geophysics,

1992, 139(3): 449

-

462.

Wiemer S, Wyss M. Minimum magnitude of completeness in earthquake catalogs: examples from Alaska, the Western United States,

and Japan[J]. Bulletin of the Seismological Society of America, 2000, 90(4): 859

-

869.

Wyss M, Hasegawa A, Wiemer S, et al. Quantitative mapping of precursory seismic quiescence before the 1989,

M

7.1 off

-

Sanriku

earthquake, Japan[J]. Annali di Geofisica, 1999, 42(5): 851

-

869.

Wyss M, Schorlemmer D, Wiemer S. Mapping asperities by minima of local recurrence time: San Jacinto

-

Elasinore fault zone[J].

Journal of Geophysical Research, 2000, 105(B4): 7 829

-

7 844.

Yang Y, Yao H J, Wu H X, et al. A new crustal shear

-

velocity model in Southwest China from joint seismological inversion and its

implications for regional crustal dynamics[J]. Geophysical Journal International, 2020, 220(2): 1 379

-

1 393.

About Qiaojia

M

S

5.0 earthquake in Yunnan on May 18, 2020

LU Xian，XIE Mengyu，ZHANG Xiaotao，WANG Yue，JIANG Xianghua，YANG Wen，

YAO Qi，ZHANG Lingkong，XIE Tao，TIAN Lei and XUE Yan

(China Earthquake Networks Center, Beijing 100045, China)

Abstract

The observed seismicity and geophysical anomalies before the Qiaojia M

S

5.0 earthquake in Yunnan

on May 18, 2020, were systematically summarized. In terms of seismicity, some medium-term and

short-term anomalies such as earthquake swarm, Benioff strain ratio, and earthquake occurrence rate

emerged before the Qiaojia M

S

5.0 earthquake. In terms of geophysical observations, measurements

such as gravity, geomagnetism, water level, water temperature, fault hydrogen, and so on were available

within 200 km of the epicenter of the Qiaojia earthquake, and most of the 15 measurements which

show abnormal before the mainshock was long-term anomalies, the anomaly ratio was relatively low

(5%), and short-term information is not obvious before the Qiaojia earthquake. Qiaojia electromagnetic

wave as a short-term anomaly could be explained to this earthquake. The anomaly was identified by the

earthquake corresponding probability spectrum of multi-parameter of the comprehensive method before

the Qiaojia earthquake. The dynamic tracking process of the Qiaojia M

S

5.0 earthquake sequence was

analyzed and summarized. The magnitude difference between the mainshock and the largest aftershock

was 1.6 (M

L

), and the energy released by the mainshock accounts for 99.64% in the whole sequence,

therefore, the earthquake sequence of the Qiaojia M

S

5.0 earthquake was a mainshock-aftershock type.

The focal mechanism for the mainshock was the strike-slip type. The seismogenic fault was located on

the NS trending secondary fault between the Lianfeng fault and Zhaotong-Ludian fault. In summary,

the process of seismic anomalies was changing gradually from long-term to short-term. Geophysical

anomalies are mainly trending type and short term anomalies did not appear or were not founded before

the Qiaojia M

S

5.0 earthquake.

Keywords: Qiaojia M

S

5.0 earthquake

，

seismicity anomaly

，

observed geophysical anomaly

，

earthquake sequence tracking

本文标签： 地震巧家发生序列活动