E-mail Alert Rss
 

物探与化探, 2020, 44(3): 480-488 doi: 10.11720/wtyht.2020.1430

地质调查·资源勘查

基于重力值标准差研究长宁6.0级地震活动性

陈挺1, 严迪1, 杨剑2, 廖国忠2, 武斌1, 冉中禹1

1. 四川省地质矿产勘查开发局物探队,四川 成都 610072

2. 中国地质调查局 成都地质调查中心,四川 成都 610082

A study of seismic activity of Changning Ms 6.0 earthquake based on standard deviation of gravity

CHEN Ting1, YAN Di1, YANG Jian2, LIAO Guo-Zhong2, WU Bin1, RAN Zhong-Yu1

1. Geophysical Prospecting Party of SBGMR,Chengdu 610072, China

2. Chengdu Geological Survey Center, China Geological Survey,Chengdu 610082, China

责任编辑: 王萌

收稿日期: 2019-09-3   修回日期: 2020-01-9   网络出版日期: 2020-06-20

基金资助: 中国地质调查局成都地质调查中心项目“四川会理—盐源及邻区1:25万区域重力调查”.  DD20190033

Received: 2019-09-3   Revised: 2020-01-9   Online: 2020-06-20

作者简介 About authors

陈挺(1984-),男,四川绵阳人,高级工程师,研究方向为重磁正反演。Email:48009735@qq.com 。

摘要

通过对长宁6.0级地震前后经固体潮和零点位移改正后的重力静态观测数据、重力仪实时观测标准差(S.D.)和舍弃数量(REJ)数据的分析,认为发震时检测到的重力值变化往往是重力仪震动滤波功能失效,导致重力仪读数错误。长宁6.0级地震在地震波到达之前未能观测到有效的重力异常变化信号,这是因为当重力仪震动滤波功能开启时,大于4~6倍标准差的高频噪声将被自动拒绝,重力场的瞬时变化在重力静态观测数据中反映不出来,但是这种变化可以通过S.D.指标来反映。由于S.D.值包括了重力仪自身噪声和环境震动等重力高频干扰信号,它可以用来衡量信号噪声的大小。通过对地震发生时重力仪在不同时刻记录S.D.值与地震震级相关性分析认为,地震震级与重力仪记录S.D.值呈明显的正相关关系。在记录到的21次地震中,在发震前有2次重力仪记录S.D.值有明显突变,有1次在发震前8 min观测到REJ值的突变,因此可以根据S.D.值的突变并结合重力值和REJ值的变化来预报地震和判定单次地震活动持续时间。

关键词： 长宁地震 ; 标准差 ; 重力异常 ; 地震预报

Abstract

Through an analysis of the static observation data of gravity, the standard deviation (S.D.) and the rejection (REJ) of real-time observation data of gravimeter before and after the Ms 6.0 earthquake in Changning, the authors hold that the change of gravity value detected during vibration is usually the mistaken reading of gravimeter caused by the failure of vibration filtering function of gravimeter. In this earthquake of Ms 6.0, no effective signal of gravity anomaly was observed before the arrival of seismic wave, and this is because high-frequency noise greater than 4 to 6 standard deviations is automatically rejected when the gravimeter vibration filtering function is turned on, transient changes in the gravitational field are not reflected in static gravity observations, but can be reflected by the S.D. indicator. Since the value of S.D. includes high frequency interference signals such as gravimeter's own noise and environment vibration, it can be used to measure the size of signal noise. By analyzing the correlation between the S.D. value recorded by gravimeter and the earthquake magnitude in different time spans when the earthquake occurs, it is concluded that the earthquake magnitude and the S.D. value recorded by gravimeter have a significant positive correlation. Among the 21 earthquakes recorded, 2 had significant changes in S.D values recorded by gravimeters prior to the earthquake and a sudden change of REJ value was observed 8 minutes before the earthquake. It can therefore be held that the earthquake can be predicted in advance and the duration of single earthquake activity can be judged according to the sudden change of S.D. value, the change of gravity value and the REJ value.

Keywords： Changning earthquake ; standard deviation ; gravity anomaly ; earthquake prediction

PDF (2209KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

陈挺, 严迪, 杨剑, 廖国忠, 武斌, 冉中禹. 基于重力值标准差研究长宁6.0级地震活动性. 物探与化探[J], 2020, 44(3): 480-488 doi:10.11720/wtyht.2020.1430

CHEN Ting, YAN Di, YANG Jian, LIAO Guo-Zhong, WU Bin, RAN Zhong-Yu. A study of seismic activity of Changning Ms 6.0 earthquake based on standard deviation of gravity. Geophysical and Geochemical Exploration[J], 2020, 44(3): 480-488 doi:10.11720/wtyht.2020.1430

0 引言

CG-5相对重力仪是加拿大Scintrex公司制造生产的的新型数字重力仪[1],分辨率为0.001×10-5 m/s2,重复率优于0.005×10-5 m/s2。广泛应用于重力研究、地震预报等领域[2]。当地震活动突然爆发时,地下岩石的空间位置以及密度分布形态发生改变,从而引起观测重力值的变化。吴云成等使用拉科斯特重力仪在辽宁西部地区进行重力复测,分析了朝阳4.9级地震前后重力异常变化情况[3]。王武星分析了昆仑山口西8.1级地震和印尼8.9级地震前后部分重力观测的变化,并对这些采样间隔在1 min以下的重力数据进行了高通滤波和频谱分析[4]。车兆宏分析了张北6.2级地震前后华北地区水准、GPS、重力、定点形变、流动形变、地应变观测资料及前兆的演化过程[5]。祝意青等通过流动重力测量,观测到较可靠的重力随时间变化数据[6,7,8,9,10]。梁伟锋等基于河西地区流动重力观测资料获得的重力场时空变化结果,研究发现重力场在活动断裂附近变化较大。这些观测记录表明重力变化能较好地反映地下深部物质运移与地壳密度变化等构造活动信息[11]。

长宁6.0级地震发生前,笔者正在云南省盐津县普洱镇进行重力静态观测,成功记录到长宁6.0级地震发生前后的重力静态观测数据、重力仪实时观测标准差(S.D.)以及舍弃数量(REJ)数据。截止2019年7月23日,在盐津县普洱镇和珙县上罗镇重力仪成功记录到地震活动21次,说明长宁地震具有主震突出、余震丰富的特点。笔者在对长宁 6.0级地震前后重力仪观测标准差数据和舍弃数量值分析的基础上,认为标准差值包括了重力仪自身噪声和环境震动等重力高频干扰信号,可以用来衡量信号噪声的大小,同时舍弃数量(REJ)的突变也指示了高频强震动干扰的存在。通过对比地震发生时重力仪在不同时刻记录的标准差值,认为标准差结合舍弃数量值可用于地震活动不确定性测量,对长宁地区发震识别与记录等具有科学与实际意义。

1 长宁震区地质背景

根据1:20万区域地质资料,长宁地区主要出露中、古生界地层(见图1),将区内地层系统列于表1。岩浆岩主要出露晚二叠系玄武岩。玄武岩主要由致密状、斑状、杏仁状玄武岩、凝灰角砾岩、凝灰岩组成,矿物成分主要为拉长石、单斜辉石,主要分布于测区西部巡场镇—坝底村一带。长宁6.0级地震震中主要出露寒武系地层向奥陶系地层过渡的白云岩和灰岩,珙县5.1级地震震中主要出露志留系罗惹坪组灰岩。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   长宁地区区域地质图(地质底图来源于91卫图助手)

Fig.1   Regional geological map of Changning area(the geological base map is derived from the assistant of 91 satellite map)


表1   长宁地区地层系统

Table 1  Stratigraphic system table of Changning area

界系统地层名称符号厚度/m主要岩性
中生界白垩系下统嘉定群K1j200~587.5长石石英砂岩
上统蓬莱镇组J3p0~452.7砂岩、钙质泥岩
侏罗系中统遂宁组J2sn192~1050钙质泥岩夹长石石英砂岩
下沙溪庙组J2s1380~915砂岩、暗紫色泥岩
三叠系上统香溪群T3~J1x330.8~680粉砂岩、石英砂岩
中统雷口坡组T2l0~336.5泥质白云岩夹粉砂岩
下统嘉陵江组T1j108.6~375.4白云岩、灰岩
飞仙关、铜街子组T1f+t16~713凝灰质页岩、含铜砂岩
古生界二叠系上统乐平组P2l0~338砂岩、泥质页岩
峨眉山玄武岩P2β0~1374玄武岩、凝灰角砾岩、凝灰岩
下统未分组P10~129粉砂岩、灰岩
志留系中统罗惹坪组S2l379.6~847.3灰岩
下统龙马溪组S1l20~388砂质灰岩、页岩
奥陶系中、上统未分组O2+30~439砂质页岩、灰岩
下统未分组O121~386页岩、灰岩
寒武系中统娄山关群∈2-3l316.5~364白云质灰岩、白云岩
高台组∈2g>104钙质粉砂岩

新窗口打开| 下载CSV


长宁地区主要地质构造为长宁背斜,该背斜东起叙永县,西至珙县天星桥一带,呈不对称的复式大背斜,背斜轴宽阔平缓,地表上在双河乡一带。长宁背斜是一个多断裂构造的背斜,断裂多发育在由古生界地层组成的背斜核部或两翼,以小断裂为主,呈高角度压性逆冲断层,震中附近断层倾角最高达85°(图1)。断裂与褶皱紧密相伴,主要断裂往往伴随褶皱而发育,构造的延伸方向与褶皱线形迹一致[12,13]。

2 观测S.D.值数据可靠性分析

标准差在重力仪读数记录里用“S.D.”表示,重力仪进行静态观测时,S.D.值每秒钟更新一次。S.D.值的计算是根据已采集到的观测重力值求取标准差,假设重力仪单点观测持续时间设置为60 s,则重力仪记录的S.D.值为60 s内所有瞬时重力值读数记录的标准差。

为了保证重力仪读数的准确性和可重复性,当打开“Auto Rejection”选项时,大于4~6倍标准差的高频噪声将被自动拒绝。因此当激活震动滤波功能时,得到的重力观测值是滤波处理后的结果。在外业观测时,S.D.指标主要用来衡量单点观测质量,当S.D.较小时,表示环境震动噪声处于较低水平,设置观测时间内的数值较接近平均值,观测质量合格;当S.D.较大时,表示环境震动噪声处于较高水平,观测时间内出现很多高频重力信号,数据质量不合格,需要重新观测。正是由于S.D.值包括了重力仪自身噪声和环境震动等重力高频干扰信号,它可以用来衡量信号噪声的大小。通过对比地震发生时重力仪在不同时刻记录的S.D.值,认为S.D.值可用于地震活动不确定性测量。

为了说明观测S.D.值数据可靠性,在珙县使用两台CG-5重力仪在同一地点进行两次观测,重力仪编号分别是40242#、41036#。两台CG-5重力仪在使用前均进行了格值标定及各项性能试验,重力仪的一致性试验精度为0.022×10-5 m/s2,试验精度结果满足规范要求。

2019年7月22日16时26分在四川宜宾市长宁县发生4.1级地震,震源深度9 km,震中位于北纬28.35°,东经104.91°。两台CG-5重力仪布置在距离震中25.5 km的珙县进行静态数据采集,表2是发震时两台重力仪记录的不同时刻静态观测重力值数据和实时观测重力值数据的标准差值(S.D.),以及高频噪声舍弃数量(REJ)统计。

表2   长宁4.1级地震重力仪观测数据统计

Table 2  Statistical table of observation data of the gravimeter of Ms 4.1 earthquake in Changning

序号观测日期观测开始
时间		观测持续时间
DUR值
/s		经固体潮和零漂改正
后的重力值
/(10-5 m·s-2)		仪器记录S.D.
/(10-5 m·s-2)		舍弃数量
REJ值/%		仪器号
12019/7/2216:24:14604845.0020.02941036#
22019/7/2216:25:34604845.0010.01919
32019/7/2216:26:54604844.9941.4044
42019/7/2216:28:14604844.9990.1991
52019/7/2216:29:34604845.0000.0581
62019/7/2216:30:54604845.0020.0359
72019/7/2216:32:14604845.0050.0230
82019/7/2216:33:34604845.0030.0164
92019/7/2216:34:54604845.0090.0111
102019/7/2216:36:14604845.0070.013
112019/7/2216:24:13604815.7360.012040242#
122019/7/2216:25:33604815.7370.0181
132019/7/2216:26:53604815.7281.3431
142019/7/2216:28:13604815.7360.1960
152019/7/2216:29:33604815.7360.0521
162019/7/2216:30:53604815.7360.0230
172019/7/2216:32:13604815.7380.0190
182019/7/2216:33:33604815.7380.0110
192019/7/2216:34:53604815.7390.010
202019/7/2216:36:13604815.7400.0080

新窗口打开| 下载CSV


重力仪观测持续时间DUR值均为60 s,为了对观测S.D.值数据可靠性进行分析,根据Excel软件中的CORREL函数计算两台重力仪观测S.D.值数据的相关系数为0.999 949,说明两台重力仪的观测S.D.值数据具有较大的相关性。另外两台重力仪观测开始时间仅相差1 s,从观测S.D.值的数值大小变化来分析,40242#重力仪相对41036#重力仪在16:26:53、16:28:13、16:29:33三个不同时刻的相对变化率分别为4.34%、1.51%、10.34%。说明利用S.D.值来分析地震活动性是可靠的。2019年7月23日22时25分在四川宜宾市长宁县发生3.0级地震,震源深度6 km,震中位于北纬28.38°,东经104.94°。两台CG-5重力仪布置在距离震中30.2公里的珙县进行静态数据采集,表3是发震时两台重力仪记录的不同时刻静态观测重力值数据和实时观测重力值数据的标准差值以及高频噪声舍弃数量值统计表。重力仪观测持续时间DUR值均为60 s,对两台重力仪的观测S.D.值数据进行了相关分析,相关系数为0.997 503,说明两台重力仪的观测S.D.值数据具有较大的相关性。表3观测S.D.值数据的相关系数小于表2,这主要是由于表3两台重力仪的观测开始时间相差20 s,即在各自的60 s观测时段内仅同步观测40 s所致。

表3   长宁3.0级地震重力仪观测数据统计

Table 3  Statistical table of observation data of the gravimeter of Ms 3.0 earthquake in Changning

序号观测日期观测开始
时间		观测持续时间
DUR值
/s		经固体潮和零漂改正
后的重力值
/(10-5 m·s-2)		仪器记录S.D.
/(10-5 m·s-2)		舍弃数量
REJ值/%		仪器号
12019/7/2322:19:42604844.9250.022641036#
22019/7/2322:21:02604844.9310.0217
32019/7/2322:22:22604844.930.01414
42019/7/2322:23:42604844.9320.0239
52019/7/2322:25:02604844.9520.38429
62019/7/2322:26:22604844.9340.0660
72019/7/2322:27:42604844.9330.0170
82019/7/2322:29:02604844.9330.0090
92019/7/2322:30:22604844.9320.0148
102019/7/2322:31:42604844.9350.00814
112019/7/2322:20:02604815.7370.013040242#
122019/7/2322:21:22604815.7370.012
132019/7/2322:22:42604815.7360.0070
142019/7/2322:24:02604815.7390.012
152019/7/2322:25:22604815.7620.49926
162019/7/2322:26:42604815.740.0480
172019/7/2322:28:02604815.7380.0110
182019/7/2322:29:22604815.7370.0080
192019/7/2322:30:42604815.7380.011
202019/7/2322:32:02604815.7370.0114

新窗口打开| 下载CSV


表2的S.D.值变化幅值较表3大,最大达 1.404×10-5 m/s2,两台重力仪发震时刻经固体潮和零漂改正后的相对重力值较发震前的变化量分别为-0.007×10-5 m/s2、-0.009×10-5 m/s2,说明观测时误差较大的高频重力数据已经被重力仪自带的滤波器过滤掉。表3的S.D.值变化幅值较小,发震时刻经固体潮和零漂改正后的相对重力值较发震前的变化量分别为0.020×10-5 m/s2、0.023×10-5 m/s2,说明部分由小震动引起的小于4倍标准差的重力观测数据未被过滤掉。据表2和表3的舍弃数量(REJ)值统计结果知,发震时刻大致与REJ极大值相对应,说明发震时的高频震动干扰最强。41036#仪器在同一时刻的REJ值均大于或等于40242#仪器,且表2和表3中41036#仪器REJ值在地震发生前有明显变化,而40242#仪器REJ值在地震发生前无明显变化,说明REJ值与重力仪的弹簧弹性系数相关,从REJ值的统计结果知,41036#仪器对高频干扰噪声更敏感,弹簧灵敏度明显高于40242#仪器,因此用于地震监测的重力仪应选择弹簧灵敏度较高的仪器。由于REJ值与重力仪自身的弹簧弹性系数相关,如需使用REJ值来预测地震,需针对单台重力仪的大量试验数据来选择合适的REJ值判断指标。

3 观测S.D.值与地震震级相关性研究

截止2019年7月23日,41036#重力仪共记录到21次地震,见图2和表4。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   长宁地震震中分布位置与地形叠合

Fig.2   The overlap map of topography and epicenter distribution of Changning earthquake


表4   发震时刻重力仪记录观测数据统计

Table 4  Statistical table of observation data recorded by gravimeter at seismic moment

序号震中经纬度地震级数震源深度
/km		仪器布设
地点		震中距
/km		重力仪记录时间观测持续时间
DUR值
/s		仪器记录
S.D.
/(10-5 m·s-2)		舍弃数量
REJ值/%		震级与S.D.的
相关系数
1北纬28.34°东经104.90°616盐津县72.906月17日22时55分6023.066100.90
2北纬28.43°东经104.77°5.116盐津县62.706月17日23时36分6011.70819
3北纬28.41°东经104.80°2.99盐津县65.106月18日01时16分600.0341
4北纬28.36°东经104.86°4.412盐津县69.506月18日05时06分600.1190
5北纬28.40°东经104.82°310盐津县66.706月19日03时00分600.0573
6北纬28.38°东经104.89°38盐津县73.206月19日04时40分600.0310
7北纬28.41°东经104.77°3.48盐津县62.606月20日01时48分600.0560
8北纬28.39°东经104.84°3.17盐津县68.306月20日10时58分600.0410
9北纬28.41°东经104.76°3.25盐津县61.306月20日13时18分600.0322
10北纬28.41°东经104.76°3.114盐津县61.306月20日14时28分600.20824
11北纬28.40°东经104.75°3.415盐津县60.106月20日21时00分600.0391
12北纬28.41°东经104.75°3.79盐津县60.506月20日22时00分600.0933
13北纬28.44°东经104.81°3.210盐津县66.906月24日04时47分600.18817
14北纬28.43°东经104.76°3.48盐津县61.906月29日22时08分600.16116
15北纬28.44°东经104.81°3.18盐津县66.906月30日21时33分600.0750
16北纬28.41°东经104.86°3.24珙县30.807月08日00时38分600.14150.94
17北纬28.38°东经104.84°38珙县27.107月22日16时18分600.29825
18北纬28.35°东经104.91°4.19珙县25.507月22日16时26分601.4044
19北纬28.38°东经104.94°36珙县30.207月23日22时25分600.38429
20北纬28.20°东经104.77°2.98珙县7.807月15日16时48分600.935471
21北纬28.14°东经104.74°3.211珙县6.207月19日01时21分6011.00622

新窗口打开| 下载CSV


表4是发震时重力仪记录静态观测数据统计表。在盐津县观测的15次记录中,当时震中距介于60.1~73.2 km时,为了研究地震震级与重力仪记录S.D值的关系,根据Excel软件中的CORREL函数计算了地震震级与对应S.D值的相关系数,计算得出相关系数为0.90,由于相关系数是研究变量之间线性相关程度的量,相关系数值越大,说明地震震级与S.D值之间线性相关程度越高。在珙县观测的6次记录中,当震中距介于25.5~30.8 km,相关系数为0.94;当震中距介于6.2~7.8 km,相关系数为1。需要指出的是,相关系数有一个明显的缺点,即当震中距介于6.2~7.8 km,仅有7月15日和7月19日两组统计样本数据,相关系数的绝对值为1,因此在观测记录较少时,仅凭相关系数较大就判定震级与重力仪记录S.D值之间有密切的线性关系是不妥当的。通过以上相关分析认为,当震中距介于60.1~73.2 km、25.5~30.8 km时,地震震级与重力仪记录S.D值呈明显的正相关关系,即在震中距相近的情况下,S.D值越大,地震震级越高。当震中距介于6.2~7.8 km时,由于统计样本数据不足,未能准确判定地震震级与重力仪记录S.D值的相关关系。表4的舍弃数量(REJ)值与震级和重力仪记录S.D值无明显的相关关系,这可能是由于当观测持续时间选择60 s时,REJ值为60 s内的统计值,无法对瞬时高频干扰进行更精细的反映。REJ值的突变显示了高频干扰的存在,在今后工作中,可以将观测持续时间选择为30 s,来进一步统计分析REJ值与震级和重力仪记录S.D值的相关关系。

4 利用S.D.值进行地震预测可行性分析

2019年6月17日22时55分在四川宜宾市长宁县发生6.0级地震, 06月17日23时36分在四川宜宾市珙县发生5.1级地震。从表5知,22时55分重力仪记录S.D.值高达23.066×10-5 m/s2,REJ值为10%。需注意的是,发震前十分钟即22时45分,S.D.值仅为0.021×10-5 m/s2,但REJ值却高达19%,REJ值的突变能否作为地震发生前兆的判断依据还需更多的试验记录来支撑。23时36分珙县5.1级地震继长宁6.0级地震后爆发,S.D.值突变为11.708×10-5 m/s2,直到06月18日00时06分,S.D.值降为0.013×10-5 m/s2,表明该次地震活动已经平息。另外从观测重力值来分析,长宁县6.0级地震和珙县5.1级地震突然爆发引起的强烈震动,导致重力仪自带的地震滤波功能失效,观测重力值与地震爆发前观测值的差值分别为4.466×10-5 m/s2、-0.055×10-5 m/s2。23时46分S.D.值为0.104×10-5 m/s2,近似为平静场中观测值的5倍,但观测重力值与地震爆发前的观测值相差仅为-0.001×10-5 m/s2,可见重力仪自带的地震滤波功能已经滤掉了重力高频干扰信号。剔除发震时刻的畸变重力数据,在云南省盐津县观测到的长宁县6.0级地震、珙县5.1级地震并未引起观测重力观测值的明显变化。直至06月18日00时06分S.D.值变为0.013×10-5 m/s2,说明两次地震震动叠加在一起的持续活动时间约为71分钟。

表5   地震前后重力仪记录静态观测数据统计

Table 5  Statistical table of static observation data recorded by gravimeter before and after earthquake

序
号		发震时间震中经纬度地震
级数		震源
深度
/km		仪器布设
地点		震中距
/km		重力仪记录时间观测持续时间
DUR值
/s		经固体潮和零漂
改正后的重力值
/(10-5 m·s-2)		仪器记录S.D.
/(10-5 m·s-2)		舍弃数量
REJ值
/%
106月17日
22时55分		北纬28.34°
东经104.90°		616盐津县普
洱镇		72.906月17日22时35分604756.9490.0180
06月17日22时45分604756.9530.02119
06月17日22时55分604761.41923.06610
06月17日23时05分604756.9510.7660
06月17日23时15分604756.9480.1550
06月17日23时25分604756.9520.0382
206月17日
23时36分		北纬28.43°
东经104.77°		5.116盐津县普
洱镇		62.706月17日23时36分604756.89711.70819
06月17日23时46分604756.9510.1040
06月17日23时56分604756.9530.02426
06月18日00时06分604756.9510.0130
306月20日
01时44分		北纬28.41°
东经104.77°		3.48盐津县普
洱镇		62.606月20日01时28分604757.3960.0190
06月20日01时38分604757.3940.0430
06月20日01时48分604757.3910.0560
06月20日01时58分604757.390.0260
06月20日02时08分604757.3910.0280
06月20日02时18分604757.3890.0150
407月15日
16时48分		北纬28.20°
东经104.77°		2.98珙县上
罗镇		7.807月15日16时29分604830.0310.0060
07月15日16时30分604830.0310.0480
07月15日16时32分604830.0430.040
07月15日16时33分604830.030.0230
07月15日16时34分604830.030.0170
07月15日16时36分604830.0310.0150
07月15日16时37分604830.030.0131
07月15日16时38分604830.0290.021
07月15日16时40分604829.9960.02242
07月15日16时41分604830.0290.020
07月15日16时42分604830.0310.0160
07月15日16时44分604830.030.0120
07月15日16时45分604830.0310.0161
07月15日16时46分604830.0330.0130
07月15日16时48分604830.0340.93547
07月15日16时49分604830.0320.0581
07月15日16时50分604830.0280.0562
07月15日16时52分604830.030.0250
07月15日16时53分604830.0290.0360
07月15日16时54分604830.030.0140

新窗口打开| 下载CSV


在记录到的21次地震中,有2次在发震前重力仪记录S.D.值有明显突变。第1次是2019年6月20日01时44分在四川宜宾市珙县发生的3.4级地震(见图3),在地震发生前的01时38分,S.D.值变为0.043×10-5 m/s2,该时刻经固体潮和零点位移改正后的静态观测重力值为4757.394×10-5 m/s2。从表5知,S.D.值与上一记录时刻相差0.024×10-5 m/s2,重力值与上一记录时刻相差-0.002×10-5 m/s2,重力值无明显变化,REJ值也无变化。由于静态观测时设置的读数时间间隔是10 min,未能记录到01时27分到01时37分之间的S.D.值,因此S.D.值可能在发震6 min之前就已发生突变。直至02时18分S.D.值变为0.015×10-5 m/s2,说明此次地震活动从发震到停止用时约34 min。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   珙县3.4级地震标准差随时间变化

Fig.3   The change of S.D. value over time in Gongxian Ms 3.4 earthquake


第2次是2019年7月15日16时48分在四川宜宾市珙县发生的2.9级地震。为了更详细地研究重力场和S.D.值随时间的变化规律,重力仪的读数时间间隔在盐津县观测时设置为10 min,在珙县观测时间设置为1 min。从图4知,在发震时刻S.D.值变为0.935×10-5 m/s2,由于发震时刻S.D.值幅值变化太大,在图4中很难观测到发震前后S.D.值的微小变化。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   珙县2.9级地震标准差随时间变化

Fig.4   The change of S.D. value over time in Gongxian Ms 2.9 earthquake


图5中剔除了16时48分的S.D.值后,可以明显观测到,在地震发生前的16时30分和16时32分,S.D.值发生了明显改变,其值分别为0.048×10-5 m/s2、0.040×10-5 m/s2,经固体潮和零点位移

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   珙县2.9级地震标准差随时间变化(剔除发震时刻S.D.突变点)

Fig.5   The change of S.D. value over time in Gongxian Ms 2.9 earthquake(the mutation point of S.D. at the time of shock was removed)


改正后的静态观测重力值分别为4 830.031×10-5 m/s2、4 830.043×10-5 m/s2,重力值与上一记录时刻相差0.012×10-5 m/s2,说明在此时刻由震动产生的部分高频重力异常未被重力仪自带的滤波器滤除掉。在16时30分,S.D.值近似为平静场中观测值的2.5倍,即在发震前18 min重力仪已记录到S.D.值的突变。值得注意的是16时40分,即在发震前8 min,S.D.值仅为0.022×10-5 m/s2,重力值与上一记录时刻相差-0.033×10-5 m/s2,REJ值高达42%,根据重力值和REJ值的变化可对地震进行提前预报。在直至16时54分S.D.值变为0.014×10-5 m/s2,说明此次地震活动从发震到停止用时约6 min。

通过以上实际记录S.D.值数据分析认为,在某些特定条件下,根据S.D.值的突变并结合重力值和REJ值的变化来提前预报地震是可行的。

5 结论与建议

通过对地震发生时重力仪在不同时刻记录S.D.值和REJ值的分析,认为S.D.值结合REJ值可用于地震活动不确定性测量。首先,对两台重力仪两次实测S.D.值数据进行了相关分析,相关系数为0.99,说明两台重力仪的观测S.D.值数据具有较大的相关性,利用S.D.值来分析地震活动性是可靠的。然后,通过对地震震级与发震时刻重力仪记录S.D值的分析,认为当震中距介于60.1~73.2 km、25.5~30.8 km时,地震震级与重力仪记录S.D值呈正相关关系。REJ值与重力仪自身的弹簧弹性系数相关,如需使用REJ值来预测地震,需针对单台重力仪的大量试验数据来选择合适的REJ值判断指标。最后,在记录到的21次地震中,有2次在发震前重力仪记录S.D值有明显突变,S.D.值近似为平静场中观测值的2倍;有1次在发震前8 min,观测到REJ值的突变。因此在某些特定条件下,根据S.D.值的突变并结合重力值和REJ值的变化来预报地震和研究单次地震从发震到结束的持续时间是可行的。

本次长宁6.0级地震,在地震波到达之前未能观测到有效的重力异常变化信号,地震前后未观测到明显的重力值变化。这可能是由于本次地震未引起地下地质体空间分布形态的较大改变且测站距离震中距离较远所致。在盐津县观测时,重力仪观测时间间隔设置为10 min,数据采集密度不够,未能在地震发生前记录到S.D.值的突变,但在发震前10 min记录到REJ值的突变。在以后的工作中,重力仪观测时间间隔设置为30~60 s,可以观测到有效的重力值和S.D.值,从而进一步研究地震活动不确定性规律。另外在震区进行小震中距重力静态观测时,地震震级与重力仪记录S.D.值是否也呈正相关关系,在进一步分析中是否要加入震源深度作为约束?这需要通过更多的采样数据来证实。

参考文献

肖凡, 何志堂, 张宏伟, 等.

CG-5型相对重力仪测量精度分析

[J]. 测绘技术装备, 2011,13(2):6-8.

[本文引用: 1]

Xiao F, He Z T, Zhang H W, et al.

Accuracy analysis of cg-5 relative gravimeter

[J] . Surveying and Mapping Equipment, 2011,13(2):6-8.

[本文引用: 1]

曹金国, 李永明, 王存军, 等.

CG-5重力仪传感器信号应用解析

[J]. 北京测绘, 2018,32(1):62-65.

[本文引用: 1]

Cao J G, Li Y M, Wang C J, et al.

The application analysis of CG-5 gravimeter sensor signal

[J]. Beijing Surveying and Mapping, 2018,32(1):62-65.

[本文引用: 1]

吴云成, 李子涛, 石作亭, 等.

朝阳4.9级地震前后重力变化

[J]. 东北地震研究, 1995,12(4):48-53.

[本文引用: 1]

Wu Y C, Li Z T, Shi Z T, et al.

Gravity variation bdfore and after Chaoyang Ms 4.9 earthquake

[J]. Seismological Research of Northeast China, 1995,12(4):48-53.

[本文引用: 1]

王武星, 马丽, 黄建平.

强地震前后重力观测中异常变化现象的研究

[J]. 地震, 2007,27(4):53-63.

[本文引用: 1]

Wang W X, Ma L, Huang J P.

Analysis of anomaly in gravity observation before and after strong earthquakes

[J]. Earthquake, 2007,27(4):53-63.

[本文引用: 1]

车兆宏, 谢觉民, 张晶, 等.

张北6.2级地震地形变前兆演化特征

[J]. 地震, 1999,19(4):315-322.

[本文引用: 1]

Che Z H, Xie J M, Zhang J, et al.

Precursory evolution characteristics of ground deformation before Zhangbei Ms 6.2 earthquake

[J]. Earthquake, 1999,19(4):315-322.

[本文引用: 1]

祝意青, 梁伟锋, 刘练.

北祁连河西地区流动重力监测与地震预测

[J]. 地震地磁观测与研究, 2002,23(1):39-44.

[本文引用: 1]

Zhu Y Q, Liang W F, Liu L.

The mobile gravity surveilance and earthquake prediction in the northern Qilian-Hexi area

[J]. Seismological and Geomagnetic Observation and Research, 2002,23(1):39-44.

[本文引用: 1]

祝意青, 胡斌, 李辉, 等.

新疆地区重力变化与伽师6.8级地震

[J]. 大地测量与地球动力学, 2003,23(3):66-69.

[本文引用: 1]

Zhu Y Q, Hu B, Li H, et al.

Variation of gravity field of Xinjiang area and Jiashi Ms 6.8 earthquake

[J]. Journal of Geodesy and Geodynamics, 2003,23(3):66-69.

[本文引用: 1]

祝意青, 胡斌, 张永志.

永登5.8级地震前后的重力场动态图像特征研究

[J]. 地壳形变与地震, 1999,19(1):71-77.

[本文引用: 1]

Zhu Y Q, Hu B, Zhang Y Z.

Study of the dynamic pattern characteristics of gravity field around Yongdeng earthquake Ms 5.8

[J]. Crustal Deformation and Earthquake, 1999,19(1):71-77.

[本文引用: 1]

祝意青, 陈兵, 张希, 等.

景泰5.9级地震前后的重力变化研究

[J]. 中国地震, 2001,17(4):356-363.

[本文引用: 1]

Zhu Y Q, Chen B, Zhang X, et al.

Gravity variation research before and after Jingtai Ms 5.9 earthquake

[J]. Earthquake Research in China, 2001,17(4):356-363.

[本文引用: 1]

祝意青, 梁伟锋, 徐云马, 等.

汶川Ms 8.0地震前后的重力场动态变化

[J]. 地震学报, 2010,32(6):633-640.

[本文引用: 1]

Zhu Y Q, Liang W F, Xu Y M, et al.

Dynamic variation of gravity field before and after Wenchuan Ms 8.0 earthquake

[J]. Acta Seismologica Sinica, 2010,32(6):633-640.

[本文引用: 1]

梁伟锋, 赵云峰, 徐云马, 等.

河西地区重力变化与2013年门源Ms 5.1地震研究

[J]. 地震工程学报, 2014,36(4):1008-1013.

[本文引用: 1]

Liang W F, Zhao Y F, Xu Y M, et al.

Study on the gravity variation in Hexi area and the Menyuan Ms 5.1 earthquake

[J]. China Earthquake Engineering Journal, 2014,36(4):1008-1013.

[本文引用: 1]

蔡一川, 王宇航.

四川长宁地区宽频带数字近震记录特征

[J]. 四川地震, 2015,154(1):1-5.

[本文引用: 1]

Cai Y C, Wang Y H.

The broadband digitalized recorded characteristic of near earthquakes in Changning area

[J]. Earthquake Reserch in Sichuan, 2015,154(1):1-5.

[本文引用: 1]

阮祥, 程万正, 张永久, 等.

四川长宁盐矿井注水诱发地震研究

[J]. 中国地震, 2008,24(3):226-234.

[本文引用: 1]

Ruan X, Cheng W Z, Zhang Y J, et al.

Research of earthquakes induced by water injections in salt mines in Changning,Sichuan

[J]. Earthquake Research in China, 2008,24(3):226-234.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址：北京市学院路29号航遥中心 邮编：100083
电话：010-62060192；62060193 E-mail：whtbjb@sina.com