本文收到日期:2018-12-25
基金项目:2019年中国地震局震情跟踪定向工作任务(项目编号:2019010109)
Henan Earthquake Agency, Zhengzhou 450016, China
0 引言
地震活动在许多中强地震发生前存在明显平静异常,是国内外地震学界普遍关注的地球物理异常现象之一。在同一个研究区域内,与背景平静的地震活动速率相比,主震前活动速率在统计上显著减小,而在时间上该速率的减小或持续到主震发生。在地震学分析预测方法中,4级地震活动异常平静是比较基本的时间进程分析方法,反映了地震活动平静地区在一定时间间隔内地震应变的积累,具有明显的物理意义。
蒋海昆等(2000)、平建军等(2001)对华北地区中强地震前平静异常进行研究,表明ML≥4.0地震平静是华北地区强震前的一个重要异常特征。目前,对于4级地震平静异常与MS≥4.5地震的相关性研究不多,仅郑先进等(2010)对华东地区做过类似研究。中原地区(30°—38°N,109°—118°E)位于华北地区南部,地质构造复杂,太行山前断裂带、汤东断裂、聊城—兰考断裂带等数条活动构造沿NE向、NNE向排列,历史上发生过多次大震,如1830年河北磁县7½级地震、1937年山东菏泽7级地震、1966年河北邢台7.2级地震等。现代中小地震仍频繁发生,据统计,1970年1月至2018年10月共发生27次MS≥4.5地震,其中8次MS≥5.0地震。中原地区地处平原,人口稠密,交通发达,防震减灾形势严峻,1996年至今被列为全国地震重点监视防御区之一,2003年以来又被确定为全国可能发生5级地震值得关注地区。因此,有必要对该区ML≥4.0地震平静异常与MS≥4.5地震的相关性进行研究。
1 资料选取
收集整理中国地震局地震预测研究所编录的全国ML≥1.0地震目录,选取中原地区(30°—38°N,109°—118°E)1970年1月—2018年10月223次ML≥4.0地震数据,采用时间进程分析方法,计算相邻2次4.0级以上地震的平静间隔,求其平均值,并选取该区MS≥4.5地震,回溯检验4.0级地震平静异常与MS≥4.5地震的对应关系。通过研究中原地区地震平静异常现象,将ML≥4.0地震平静大于等于90天判定为异常,认为地震平静异常出现后,中原地区有61.8%的概率发生MS≥4.5地震。
2 地震平静异常判别
对于选取的223次中原地区1970年以来ML≥4.0地震,分析震级—时间关系,统计相邻2次地震平静间隔(图 1),结果如下:①154次地震平静间隔小于90天,8次地震平静间隔在90—100天,23次地震平静间隔在101—150天,12次地震平静间隔在151—200天,10次地震平静间隔在201—250天,9次地震平静间隔在251—300天,6次地震平静间隔在大于300天;②平静间隔200天以内的地震占88.7%,最大平静间隔为1 537天;③223次ML≥4.0地震平均平静间隔约为89.4天,平均复发周期为87.9天。
综合分析认为,中原地区ML≥4.0地震平静异常判别标准为大于等于90天。
3 4级地震平静异常与MS≥4.5地震对应分析
3.1 MS≥4.5地震统计
据统计,中原地区1970年1月至2018年10月发生MS≥4.5地震27次(同一地震序列只计一次最大地震),其中MS≥5.0地震8次。MS≥4.5地震分布见图 2,具体参数见表 1。
表 1(Table 1
表 1 中原地区1970年至今MS≥4.5地震目录Table 1 The basic parameters of earthquakes with MS≥4.5 since 1970 in Central China
序号 |
时间 |
φN/° |
λE/° |
MS |
深度/km |
地点 |
1 |
1970-08-10 |
35.70 |
116.88 |
4.6 |
34 |
山东曲阜 |
2 |
1971-06-05 |
37.43 |
113.42 |
4.8 |
- |
山西和顺 |
3 |
1972-10-12 |
36.97 |
114.32 |
4.8 |
16 |
河北沙河 |
4 |
1973-03-11 |
31.37 |
116.17 |
4.6 |
6 |
安徽霍山 |
5 |
1973-11-29 |
32.88 |
111.52 |
4.7 |
7 |
河南淅川 |
6 |
1974-06-06 |
37.60 |
115.13 |
4.5 |
24 |
河北宁晋 |
7 |
1978-06-05 |
35.35 |
113.92 |
4.5 |
10 |
河南新乡 |
8 |
1979-03-02 |
33.18 |
117.42 |
5.0 |
11 |
安徽固镇 |
9 |
1979-05-22 |
31.10 |
110.47 |
5.1 |
16 |
湖北秭归 |
10 |
1979-06-19 |
37.10 |
111.87 |
5.1 |
- |
山西介休 |
11 |
1980-03-09 |
37.28 |
112.07 |
4.6 |
22 |
山西平遥 |
12 |
1980-08-02 |
36.05 |
113.95 |
4.7 |
19 |
河南林州 |
13 |
1981-11-09 |
37.43 |
114.97 |
5.8 |
20 |
河北隆尧 |
14 |
1982-03-11 |
33.20 |
110.28 |
4.5 |
25 |
陕西山阳 |
15 |
1983-11-07 |
35.27 |
115.28 |
6.0 |
12 |
山东菏泽 |
16 |
1985-11-30 |
37.23 |
114.82 |
5.2 |
15 |
河北任县 |
17 |
1989-12-25 |
35.58 |
111.30 |
4.9 |
11 |
山西侯马 |
18 |
1998-07-11 |
34.95 |
110.63 |
4.6 |
22 |
山西永济 |
19 |
2002-04-22 |
37.27 |
114.42 |
5.0 |
10 |
河北内丘 |
20 |
2002-09-03 |
37.85 |
112.58 |
4.6 |
18 |
山西太原 |
21 |
2003-11-25 |
36.17 |
111.63 |
4.6 |
19 |
山西洪洞 |
22 |
2010-01-24 |
35.57 |
110.77 |
4.8 |
8 |
山西河津 |
23 |
2010-10-24 |
34.07 |
114.65 |
4.6 |
6 |
河南太康 |
24 |
2013-12-16 |
31.10 |
110.45 |
5.1 |
10 |
湖北巴东 |
25 |
2014-03-30 |
30.92 |
110.78 |
4.6 |
5 |
湖北秭归 |
26 |
2016-03-12 |
34.93 |
110.83 |
4.5 |
10 |
山西盐湖 |
27 |
2018-10-11 |
31.03 |
110.47 |
4.5 |
6 |
湖北秭归 |
|
表 1 中原地区1970年至今MS≥4.5地震目录
Table 1 The basic parameters of earthquakes with MS≥4.5 since 1970 in Central China |
3.2 MS≥4.5地震映震统计
中原地区1970年1月—2018年10月ML≥4.0地震平静异常与MS≥4.5地震对应统计结果见表 2。由表 2可知,4级地震平静大于等于90天平静间隔异常标准的共68次,其后1年半内对应MS≥4.5地震的有42次,占61.8%,其中有13次被MS≥4.5地震直接打破;虚报26次,占38.2%;而研究区内共发生27次MS≥4.5地震,震前有平静异常的共25次,占92.6%,漏报2次,占7.4%。2次漏报中的一次地震发震时间为1970年8月10日山东曲阜MS 4.6,漏报原因可能与无1970年前地震数据有关。
表 2(Table 2
表 2 中原地区1970年以来ML≥4.0地震平静异常与MS≥4.5地震对应统计Table 2 Relation statistics between the quiescence of earthquake with ML≥4.0 and earthquake with MS≥4.5 since 1970 in Central China
序号 |
ML≥4.0地震平静异常参数 |
| 地震“三要素” |
发震间隔 T3-T2/月 |
地震状态 |
开始时间T1 |
结束时间T2 |
持续天数 |
发震时间T3 |
参考地名 |
MS |
1 |
1972-04-03 |
1972-09-07 |
157.4 |
| 1972-10-12 |
河北沙河 |
4.8 |
1.1 |
地震活跃 |
2 |
1972-11-17 |
1973-03-11 |
114.3 |
1973-03-11 |
安徽霍山 |
4.5 |
0 |
3 |
1973-03-31 |
1973-09-22 |
193.8 |
1973-11-29 |
河南淅川 |
4.7 |
2.2 |
4 |
1973-12-04 |
1974-03-07 |
92.4 |
1974-06-06 |
河北宁晋 |
4.5 |
3.0 |
5 |
1974-03-07 |
1974-06-06 |
91.8 |
1974-06-06 |
河北宁晋 |
4.5 |
0 |
6 |
1974-07-18 |
1975-04-22 |
277.6 |
|
|
|
| 地震平静 |
7 |
1975-04-22 |
1976-02-14 |
298.2 |
|
|
|
|
8 |
1976-02-14 |
1976-08-30 |
197.8 |
|
|
|
|
9 |
1976-08-30 |
1977-01-01 |
124.6 |
1978-06-05 |
河南新乡 |
4.5 |
17.2 |
10 |
1977-03-10 |
1977-07-09 |
121.0 |
1978-06-05 |
河南新乡 |
4.5 |
10.9 |
11 |
1977-10-14 |
1978-05-24 |
221.6 |
1978-06-05 |
河南新乡 |
4.5 |
0.4 |
12 |
1978-06-05 |
1978-10-01 |
117.5 |
1979-03-02 |
安徽固镇 |
5.0 |
5.0 |
13 |
1978-10-01 |
1979-03-02 |
152.2 |
1979-03-02 |
安徽固镇 |
5.0 |
0 |
|
|
|
| 1979-05-22 |
湖北秭归 |
5.1 |
2.7 |
地震活跃 |
|
|
|
| 1979-06-19 |
山西介休 |
5.1 |
3.6 |
14 |
1979-06-19 |
1980-03-09 |
263.7 |
1980-03-09 |
山西平遥 |
4.6 |
0 |
15 |
1980-03-09 |
1980-08-02 |
145.9 |
1980-08-02 |
河南林州 |
4.7 |
0 |
16 |
1980-08-07 |
1980-12-15 |
130.2 |
1981-11-09 |
河北隆尧 |
5.8 |
10.8 |
17 |
1981-01-22 |
1981-07-05 |
164.9 |
1981-11-09 |
河北隆尧 |
5.8 |
4.2 |
18 |
1981-07-05 |
1981-11-09 |
126.1 |
1981-11-09 |
河北隆尧 |
5.8 |
0 |
|
|
|
| 1982-03-11 |
陕西山阳 |
4.5 |
4.0 |
19 |
1983-03-09 |
1983-09-15 |
189.7 |
1983-11-07 |
山东菏泽 |
6.0 |
1.7 |
20 |
1983-11-07 |
1984-09-25 |
323.0 |
1985-11-30 |
河北任县 |
5.2 |
14.1 |
21 |
1984-09-25 |
1985-01-26 |
123.7 |
1985-11-30 |
河北任县 |
5.2 |
10.1 |
22 |
1985-01-26 |
1985-05-02 |
95.6 |
1985-11-30 |
河北任县 |
5.2 |
6.9 |
23 |
1985-06-25 |
1985-11-30 |
158.1 |
1985-11-30 |
河北任县 |
5.2 |
0 |
24 |
1986-02-18 |
1986-11-17 |
272.4 |
|
|
|
|
25 |
1986-12-16 |
1988-02-12 |
423.3 |
|
|
|
|
26 |
1988-04-15 |
1988-12-16 |
244.9 |
1989-12-25 |
山西侯马 |
4.9 |
12.3 |
27 |
1989-02-03 |
1989-06-02 |
118.1 |
1989-12-25 |
山西侯马 |
4.9 |
6.7 |
28 |
1989-07-23 |
1989-12-14 |
144.4 |
1989-12-25 |
山西侯马 |
4.9 |
0.3 |
29 |
1989-12-25 |
1990-07-14 |
201.3 |
|
|
|
| 地震平静 |
30 |
1990-08-03 |
1990-12-24 |
142.8 |
|
|
|
|
31 |
1990-12-24 |
1991-04-17 |
114.1 |
|
|
|
|
32 |
1991-04-25 |
1991-08-21 |
118.0 |
|
|
|
|
33 |
1992-02-26 |
1992-09-02 |
188.9 |
|
|
|
|
34 |
1992-09-02 |
1993-01-01 |
121.2 |
|
|
|
|
35 |
1993-02-19 |
1993-06-26 |
127.1 |
|
|
|
|
36 |
1993-08-06 |
1994-04-26 |
262.8 |
|
|
|
|
37 |
1994-04-26 |
1998-07-11 |
1537.6 |
1998-07-11 |
山西永济 |
4.6 |
0 |
38 |
1998-07-11 |
1998-12-23 |
164.4 |
|
|
|
| 地震活跃 |
39 |
1998-12-26 |
1999-05-20 |
145.5 |
|
|
|
|
40 |
1999-05-20 |
1999-12-30 |
223.3 |
|
|
|
|
41 |
2000-01-28 |
2000-04-29 |
91.6 |
|
|
|
|
42 |
2000-04-29 |
2001-04-27 |
363.4 |
2002-04-22 |
河北内丘 |
5.0 |
11.8 |
43 |
2001-07-25 |
2001-12-13 |
141.2 |
2002-04-22 |
河北内丘 |
5.0 |
4.3 |
44 |
2001-12-13 |
2002-04-22 |
129.6 |
2002-04-22 |
河北内丘 |
5.0 |
0 |
|
|
|
| 2002-09-03 |
山西太原 |
4.6 |
4.3 |
45 |
2002-12-25 |
2003-04-27 |
122.8 |
2003-11-25 |
山西洪洞 |
4.6 |
6.9 |
46 |
2003-04-27 |
2003-11-25 |
212.0 |
2003-11-25 |
山西洪洞 |
4.6 |
0 |
47 |
2003-11-25 |
2004-05-27 |
184.0 |
|
|
|
| 地震平静 |
48 |
2004-05-27 |
2005-02-14 |
263.0 |
|
|
|
|
49 |
2005-07-21 |
2006-04-09 |
262.6 |
|
|
|
|
50 |
2006-04-09 |
2006-07-26 |
107.5 |
|
|
|
|
51 |
2006-07-26 |
2006-10-27 |
93.6 |
|
|
|
|
52 |
2006-10-28 |
2007-06-03 |
217.7 |
|
|
|
|
53 |
2007-06-03 |
2008-01-03 |
214.2 |
|
|
|
|
54 |
2008-03-24 |
2008-11-22 |
242.7 |
2010-01-24 |
山西河津 |
4.8 |
14.1 |
55 |
2008-11-22 |
2009-04-06 |
135.3 |
2010-01-24 |
山西河津 |
4.8 |
9.6 |
56 |
2009-04-06 |
2008-11-05 |
212.4 |
2010-01-24 |
山西河津 |
4.8 |
2.6 |
57 |
2010-04-07 |
2010-10-24 |
200.1 |
2010-10-24 |
河南太康 |
4.7 |
0 |
地震活跃 |
58 |
2011-03-08 |
2011-06-17 |
101.4 |
|
|
|
|
59 |
2011-10-13 |
2013-04-04 |
539.2 |
2013-12-16 |
湖北巴东 |
5.1 |
8.3 |
60 |
2013-04-04 |
2013-12-16 |
255.8 |
2013-12-16 |
湖北巴东 |
5.1 |
0 |
61 |
2013-12-16 |
2014-03-27 |
100.5 |
2014-03-30 |
湖北秭归 |
4.6 |
0.1 |
62 |
2014-05-27 |
2015-01-18 |
235.5 |
2016-03-12 |
山西盐湖 |
4.5 |
13.8 |
63 |
2016-03-23 |
2016-01-20 |
303.2 |
2016-03-12 |
山西盐湖 |
4.5 |
1.7 |
64 |
2016-03-27 |
2016-12-18 |
266.3 |
|
|
|
|
65 |
2017-02-23 |
2017-05-30 |
96.0 |
2018-10-11 |
湖北秭归 |
4.5 |
16.4 |
66 |
2017-09-10 |
2018-02-09 |
151.8 |
2018-10-11 |
湖北秭归 |
4.5 |
8.0 |
67 |
2018-04-06 |
2018-07-10 |
94.4 |
2018-10-11 |
湖北秭归 |
4.5 |
3.0 |
68 |
2018-07-10 |
2018-10-11 |
93.3 |
2018-10-11 |
湖北秭归 |
4.5 |
0 |
|
表 2 中原地区1970年以来ML≥4.0地震平静异常与MS≥4.5地震对应统计
Table 2 Relation statistics between the quiescence of earthquake with ML≥4.0 and earthquake with MS≥4.5 since 1970 in Central China |
根据公式R值=报对地震数/应报地震数-预报占用时间/研究总时间,计算其地震预报效能评分,得到R = 25/27 - 221.2/586 = 0.547,中原地区ML≥4.0地震平静异常大于90天后,统计有61.8%的平静异常与中原地区MS≥4.5地震存在一定对应关系。经查97.5%置信水平R0 = 0.171,可见R>R0,表明中原地区4级地震平静与其后对应MS ≥ 4.5地震至少有97.5%的置信度。上述研究表明,ML≥4.0地震平静异常90天指标具有一定预测价值。
由表 2可见,映震距4级地震平静异常结束时间:有13次4级地震平静被MS≥4.5地震直接打破,占28%;4级地震平静异常结束后6个月内发生MS ≥ 4.5地震30次,占65%;9个月内发震35次,占76%;12个月内发震40次,占87%;18个月内发震46次,占100%。因此出现ML ≥ 4.0地震平静异常后,即可预测未来1.5年内存在发生MS≥4.5地震的可能,但优势发震时间为1年内(87%)。
4 地震虚报分析
使用中原地区ML≥4.0地震平静异常大于等于90天,预测其后该区存在发生MS≥4.5地震的可能虽具有一定价值,但虚报率较高,达38.2%。如能减少虚报率,可使该方法具有更大的应用价值。
根据研究区域1970—2018年MS≥4.5地震M—t图,中原地区MS≥4.5地震活动存在较明显的地震活跃时段和地震平静时段(图 3,表 3)。由表 2和表 3可知,对中原地区26次虚报的地震中,有8次发生在该区域地震活跃时段,18次发生在该区域地震平静时段,统计结果见表 4。
表 3(Table 3
表 3 中原地区地震周期划分Table 3 The seismic period division in Central China
地震活跃时段 |
| 地震平静时段 |
序号 |
起止时间 |
持时/月 |
N4.5/N5 |
对应最大地震 |
序号 |
起止时间 |
持时/月 |
时间 |
地点 |
M |
1 |
1970-08-10—1974-06-06 |
45.9 |
6 / 0 |
1972-10-12 |
河北沙河 |
4.8 |
| 1 |
1974-06-07—978-06-04 |
47.9 |
2 |
1978-06-05—1989-12-25 |
138.7 |
11 / 6 |
1983-11-07 |
山东菏泽 |
6.0 |
2 |
1989-12-26—998-07-10 |
102.5 |
3 |
1998-07-11—2003-11-25 |
64.5 |
4 / 1 |
2002-04-22 |
河北内丘 |
5.0 |
3 |
2003-11-26—2010-01-23 |
73.9 |
4 |
2010-01-24—2018-10-31 |
105.3 |
6 / 1 |
2013-12-16 |
湖北巴东 |
5.1 |
|
|
|
注:N4.5指4.5级地震发生次数,N5指5级地震发生次数 |
|
表 3 中原地区地震周期划分
Table 3 The seismic period division in Central China |
表 4(Table 4
表 4 4级地震平静异常分布状况Table 4 The distribution of quiet anomalies for ML 4 earthquakes
地震活跃时段 |
| 地震平静时段 |
映震数 |
映震率 |
虚报数 |
虚报率 |
映震数 |
映震率 |
虚报数 |
虚报率 |
35 |
81.4% |
8 |
18.6% |
| 7 |
28.0% |
18 |
72.0% |
|
表 4 4级地震平静异常分布状况
Table 4 The distribution of quiet anomalies for ML 4 earthquakes |
由表 4明显可见,地震活跃时段映震率明显高于虚报率,而地震平静时段虚报率则明显高于映震率。若把地震平静时段中出现的18次虚报异常进行适当的技术处理,则可明显减少地震虚报率。在地震预测实践中,或可进行以下操作:在区域地震平静时段,如4级地震平静异常后1年内仍无MS≥4.5地震对应,则不再对该异常继续跟踪。但根据前文提及的1年内发震率占87%,可能还要承担13%的失败风险。由于震例太少,如何更有效地减小地震虚报带来的影响,有待积累资料进行深入探讨。
5 结论
综上所述,可以得到以下结论:①中原地区ML≥4.0地震平静异常与异常结束后MS≥4.5地震存在较好的相关性,映震率达61.8%;②中原地区出现ML≥4.0地震平静异常,可预测该区未来1.5年内MS≥4.5地震的发生,优势发震时间在1年内(占87%);③地震预测效能评分R值=0.547,大于97.5%置信度对应的最低值R0(0.171),说明R值至少有97.5%的置信度,表明该方法具有一定预测价值;④中原地区ML≥4.0地震平静异常大于90天后,统计有61.8%的平静异常与中原地区MS≥4.5地震相关,但地震虚报较多(占38.2%),特别是在中原地区地震平静时段内,地震虚报率高达72%,如结合研究区地震活动平静时段分析,可明显减少地震虚报率。