文章快速检索     高级检索
  大地测量与地球动力学  2022, Vol. 42 Issue (12): 1276-1280, 1299  DOI: 10.14075/j.jgg.2022.12.012

引用本文  

李梦莹, 张志宏, 焦明若, 等. 2018-05-28松原MS5.7地震地电场变化特征研究[J]. 大地测量与地球动力学, 2022, 42(12): 1276-1280, 1299.
LI Mengying, ZHANG Zhihong, JIAO Mingruo, et al. Study on Variation Characteristics of Natural Electric Field in Songyuan MS5.7 Earthquake on May 28, 2018[J]. Journal of Geodesy and Geodynamics, 2022, 42(12): 1276-1280, 1299.

项目来源

中国地震局地球物理研究所基本科研业务费专项(DQJB21B34);中国地震局震情跟踪定向工作任务(2022010412,2022010301);辽宁省地震局科技项目(2022025, 2022035)。

Foundation support

Special Fund of the Institute of Geophysics, CEA, No. DQJB21B34;The Earthquake Tracking Task of CEA, No. 2022010412, 2022010301; Science and Technology Project of Liaoning Earthquake Agency, No.2022025, 2022035.

通讯作者

张志宏,高级工程师,主要从事地震电磁研究,E-mail: 568950612@qq.com

Corresponding author

ZHANG Zhihong, senior engineer, majors in earthquake electromagnetic, E-mail: 568950612@qq.com.

第一作者简介

李梦莹,工程师,主要从事地震电磁研究,E-mail: 1593919580@qq.com

About the first author

LI Mengying, engineer, majors in earthquake electromagnetic, E-mail: 1593919580@qq.com.

文章历史

收稿日期:2022-02-21
2018-05-28松原MS5.7地震地电场变化特征研究
李梦莹1     张志宏1     焦明若1     方禹心2     
1. 辽宁省地震局,沈阳市黄河北大街44号,110034;
2. 沈阳地震监测中心站,沈阳市东陵路208巷32号,110161
摘要:通过研究东北地区地电场观测资料,结合2018-05-28吉林松原MS5.7地震,分别对震中周边自然电场ESP和大地电场ET地震前后背景、临震、震后响应时序变化进行对比分析,并应用岩体裂隙水电荷移动渗流模型计算研究区电荷移动的优势方位。结果表明,震前绥化、林甸、肇东、望奎自然电场ESP出现明显的小幅度突跳或大幅度跃变,电荷移动的优势方位震前3个月内出现45°挤压或拉张异常特征。
关键词松原MS5.7地震自然电场大地电场裂隙结构异常机理

根据中国地震台网(http://www.ceic.ac.cn)正式测定,2018-05-28 01:50吉林松原发生MS5.7地震。本次地震为2013-11-23吉林前郭MS5.8地震后的首次较大地震,震中区也是中国东部1975-02-04辽宁海城MS7.3地震和1976-07-28河北唐山MS7.8地震后地震活动性最强的地区,对研究地震前地球物理场中短期异常有着重要意义。

我国地电场观测与地震预测研究始于1966年邢台MS7.2地震后,曾记录到不少震前异常变化[1],Varotsos等[2-3]用多极距观测法解决了地电场观测系统噪声问题,自此国内地电场地震预测预报方法得到快速发展[4-8]。本文使用2018-05-28吉林松原MS5.7地震震前1 a以上及震后6个月震中500 km范围内13个地电场台站分钟值数据,分析地震前后自然电场变化,得到震前自然电场ESP异常特征,并结合岩体裂隙水电荷渗流模型计算异常台站地电场优势方位角a,讨论地电场优势方位角在地震预测方面的优势。

1 原理及方法介绍 1.1 自然电场和大地电场

地电场是指固体地球内部和外部的各种非人工电流系与地球介质相互作用,在地球表面产生的电场,包括大地电场ET和自然电场ESPET与磁层和电离层中电流系的运动有关,ESP则是地壳内部各类物理化学作用引起的正负电荷分离产生的地电场[1]ETESP通过测量地表每个测向两点之间的电位差,得到其随时间的变化,如图 1所示,通过埋设在地表下的测量电极AB,以一定时间间隔测量2个电极的自然电位差序列(VAB)i(i=N-1,N-2,⋯,2,1,0),其中i=0表示当前时间,i=1表示当前时间之前第1个时间点,i=N-1表示当前时间之前第N-1个时间点。按式(1)计算得出的自然电位差序列(VAB)i算术平均值,即为当前时间的自然电场ESP;按式(2)计算此方向上的大地电场ET分量值,其中负号表示大地电场的方向由高电位指向低电位。

$ V_{\mathrm{SP}}=\frac{1}{N} \sum\limits_{i=0}^{N-1}\left(V_{A B}\right)_i $ (1)
图 1 地电场观测基本原理示意图 Fig. 1 Schematic diagram of basic principle of geoelectric field observation

式中,(VAB)i(i=N-1,N-2,⋯,2,1,0)为电极AB之间的自然电位差值序列,单位mV;N为参与计算的自然电位差值序列数据总数,不少于24 h之内的所有测量数据个数。

$ E_{\mathrm{T}}=-\frac{\left(V_{A B}\right)_0-V_{\mathrm{SP}}}{\overline{A B}} $ (2)

式中,(VAB)0为当前时间点的电位差值,单位mV;VSP为式(1)计算得到的自然电场值,单位mV;AB为2个测量电极之间的距离,单位km。

大地电场矢量由2个正交方向的大地电场分量值合成,2个正交方向分别为x轴(东向分量)和y轴(北向分量),由2个大地电场分量值ExEy按式(3)和式(4)计算大地电场的幅值E和方位角a

$ E=\sqrt{E_x^2+E_y^2} $ (3)
$ a=\tan ^{-1} \frac{E_x}{E_y} $ (4)

式中,Ex为大地电场东向分量,Ey为大地电场北向分量,单位均为mV/km。

地电场的观测数据中,ESP为自然电场,ET为大地电场,Er为干扰成分,则地电场E可表达为:

$ E=E_{\mathrm{SP}}+E_{\mathrm{T}}+E_r $ (5)

设地电场观测分钟值数据为Ei,日均值计算自然电场ESP的表达式可简化为:

$ E_{\mathrm{SP}} \approx \frac{1}{1\;440} \sum\limits_{i=0}^{1\;439} E_i $ (6)

应用式(6)能够基本消除大地电场ET和干扰成分Er

1.2 地电场优势方位角a

中国学者基于大地电场ET日变波形的时域和频域特征,提出了大地电场ET日变波源于电离层Sq电流和潮汐力学说,即大地电场的潮汐机理[5-6]。谭大诚等[8]在潮汐机理的基础上建立了大地电场ET的岩体裂隙水(电荷)渗流(移动)模型,解释了大地电场ET强度和方向的变化是地下介质应力的积累改变岩体裂隙结构的过程。图 2(b)为裂隙水电荷日周期移动的理想模型,岩体结构发育程度决定了实际的电荷移动方向。用a表示场地岩体裂隙水电荷移动的主渗流方位,也可称之为地电场优势方位角a

图 2 ET岩体裂隙水(电荷)渗流(移动)模型 Fig. 2 ET model of fissure water in rock with seepage(movement)or charges

当台站地电场NS和NW测向之间的相关性高时,地电场ET的优势方位角a(北偏东)计算公式如下:

$ a \approx 180^{\circ}-\left(\frac{180}{\pi}\right) \tan ^{-1}\left(\sqrt{2} \frac{\sum\limits_{i=1}^{10} A_{\mathrm{NW}(i)}}{\sum\limits_{i=1}^{10} A_{\mathrm{NS}(i)}}-1\right) $ (7)

式中,ANW(i)ANS(i)分别为NW、NS测向第i阶潮汐谐波振幅,Ai的计算公式如下:

$ A_i=\sqrt{a_i^2+b_i^2} $ (8)

其中,ai=$ \frac{2}{n} \sum\limits_{t=1}^n y_t \cos \left(\frac{2 \pi i t}{n}\right)$bi=$\frac{2}{n} \sum\limits_{i=1}^n \sin \left(\frac{2 \pi i t}{n}\right) $

对于一个数据序列yt(时间序列总数为n),数学上可表示为:

$ y_t=\bar{y}+\sum\limits_{i=1}^{n / 2}\left[a_i \cos \left(\frac{2 \pi i t}{n}\right)+b_i \sin \left(\frac{2 \pi i t}{n}\right)\right] $ (9)
2 台站地质构造概况

本文选取2018年松原MS5.7地震震中500 km范围内采样时间间隔为1 min、测量频段为0~0.005 Hz、“双L”6测道长短极距布设、不极化电极的13套ZD9A-Ⅱ型仪器和2套GEF-2型仪器地电场台站[9](表 1),图 3为台站所处现今地质构造区域。钟以璋等[10]通过研究卫星影像数据和地球物理场观测数据发现,研究区地壳构造由断块组成,地震震中沿断裂带分布,且展布方向与区域内主要断裂带的走向大体一致。研究区地属松辽盆地,该盆地为中国东北部的一个大型中-新生代NE向沉积盆地,地跨中国东北三省和内蒙古自治区[11]。21世纪以来,松辽盆地及周边地区地震活动明显增强,2002年至今已连续发生5级以上天然地震16次,高立新等[12-13]认为2011年日本MW9.0地震对中国东北地区地球动力学环境影响显著,地壳运动状态由挤压转为拉张。2011年以来中国东北地区相继发生辽宁灯塔MS5.1地震、内蒙古科尔沁MS5.3地震、吉林前郭MS5.8震群等7次5.0级以上浅源地震。

表 1 研究区地电台站概况 Tab. 1 Overview of local stations in the study area

图 3 2018年松原MS5.7地震震中及500 km内地电场台站分布 Fig. 3 Schematic diagram of the distribution of the 2018 Songyuan MS5.7 earthquake epicenter and the distribution of electric field stations within 500 km
3 分析研究 3.1 自然电场变化

根据场源的不同,可将地电场E分为大地电场ET和自然电场ESP,大地电场ET的起源是空间Sq电流系和固体潮,自然电场ESP则源于地下介质的物理、化学作用。自然电场ESP包括地下金属矿体与周围溶液经过氧化还原反应的氧化还原场、地下电介质在岩体裂隙中流动的过滤电场及由2种电解质浓度差异引起的接触-扩散电场。膨胀和压力实验表明,岩体在应力变化作用下,能够激发自然电场的变化[1]。因此,监测孕震过程的应力变化可能观测到自然电场短期异常[14]

图 4为2017~2018年本文所选台站中自然电场ESP变化较大,且观测系统和场地较好的4个场地数据。图 4(a)表明,2017-03林甸台ESP的EW、NW测向开始出现小幅度异常跃变,持续至7月初跃变幅度峰值约为400 mV ·km-1,2018-03开始第2次异常跃变,且异常测向变幅方向相反、幅度差较小;5月初ESP的EW、NW测向恢复背景值,松原MS5.7地震后ESP没有观测到同震响应。图 4(b)表明,2018-02绥化台ESP三测向出现同步异常,跃变幅度峰值较小,小于50 mV ·km-1图 4(c)表明,2018-03初望奎台ESP的NS、N45°W测向出现跃变异常,跃变峰值约为500 mV ·km-1,且两异常测向的变幅同步性较高,地震发生在ESP异常幅度恢复期,震后6个月ESP恢复正常。图 4(d)表明,2018-03初肇东台ESP的EW、NS测向观测值突增,异常起始幅度已达峰值500 mV ·km-1

图 4 异常台站在地震前后地电场变化(2017-1-30~2019-01-30) Fig. 4 Variation of electric field of abnormal stations before and after earthquake (January 30, 2017 to January 30, 2019)

表 2可知,绥化台ESP的异常幅度相对于背景值变化并不显著,异常信度最低;肇东台、望奎台和林甸台异常变化显著,异常信度较高,其中林甸台的异常幅度相对较小,异常起始时间较早。从图 4也可看出,肇东台、望奎台和林甸台ESP异常测向表现为两道,其中望奎台和林甸台同属富裕-明水断裂(图 3),而肇东台ESP异常变化速率显著,短时期内达到了峰值,发震构造研究结果[15-16]也支持肇东台的异常特征。

表 2 松原MS5.7地震地电场异常台站数据统计 Tab. 2 Statistics of electric field in anomalous stations of the Songyuan MS5.7 earthquake
3.2 大地电场岩体裂隙优势方位

图 5为基于大地电场岩体裂隙水(电荷)渗流(移动)模型计算的4个异常场地岩体裂隙水主渗流方位。图 5(a)为林甸台大地电场优势方位角α的变化散点分布,可以看出,2017-03林甸台方位角α的突变范围出现了大幅度的发散,Δα约为90°,对比图 4(a)中林甸台自然电场ESP的变化开始时间发现,两者的变化初始时间明显对应。由此可见,2017-03林甸台自然电场ESP大幅变化的物理机理可能是场地岩体裂隙结构在应力加卸载作用下发生了较大改变,表现为间断的岩体裂隙偏转或剪裂现象。图 5(c)5(d)为望奎台和肇东台大地电场优势方位角α的变化散点分布,可以看出,两者方位角α的突变范围呈大幅度的挤压状态,Δα在45°左右。对比图 4(c)4(d)可知,望奎台自然电场ESP同步发生异常变化,松原MS5.7地震后望奎台大地电场优势方位角α的挤压状态持续至9月,后恢复发散状态,而肇东台震前由挤压恢复至扩张。

图 5 异常台站大地电场优势方位角α的变化散点分布 Fig. 5 Variation scatter diagram of the dominant azimuth angle α of the telluric field at anomalous stations
4 结语

李永生等[15]的全波形矩张量反演结果表明,扶余-肇东断裂东段为2018年松原MS5.7地震的发震构造;李君等[16]利用双差定位法对地震序列进行重新定位,并利用MSATSI软件反演得到松原地区应力场,结果表明,松原地区应力场既受太平洋板块对北东亚板块的西向俯冲,又受区域构造运动的控制,在中国东北地区产生了近EW向的主压应力。因此,自然电场变化量ΔESP最大的台站应位于扶余-肇东断裂,且EW测向异常显著,这就解释了肇东台ESP的EW测向异常显著和岩体裂隙优势方位角出现高达45°挤压的异常特征。震中周边场地的孕震过程是应力不断加卸载的过程,岩体裂隙结构会因应力加卸载的变化而变化。但实际结果是,场地岩体结构的差异导致不同场地观测到的地电场ESP变化差异较大,这可能是异常台站只有两道显著异常的原因。

本文通过研究近年来我国东北地区地震活动最为活跃的吉林松原地区地电场观测资料,结合2018-05-28松原MS5.7地震,分别对地震前后自然电场ESP和大地电场ET的背景、临震、震后响应时序变化进行对比分析,并应用岩体裂隙水电荷移动渗流模型计算了研究区所有地电场台站的电荷移动优势方位,得到如下结论:

1) 多极距、多方向的地电场布极方式可有效识别干扰,但不易提取前兆异常信息;

2) 2018-05-28松原MS5.7地震前后,绥化台、林甸台、肇东台、望奎台自然电场ESP的异常起始时间存在准同步特征;

3) 基于岩体裂隙水电荷渗流模型及大地电场优势方位的结果,可得出岩体结构发育状况;

4) 松原MS5.7地震前后异常台站自然电场ESP的大幅度变异现象,可能是松原地区的区域应力场受到构造板块应力的影响, 望奎台及肇东台的大地电场优势方位角α呈挤压异常。

利用岩体裂隙水电荷渗流模型计算得到的电荷移动方向可有效捕捉场地孕震岩体结构的微变化,由于孕震过程复杂,对该研究方法的适用可能存在区域性,后续需要更多震例支持。

致谢: 感谢甘肃省地震局谭大诚研究员提供地电分析程序,黑龙江省地震局和吉林省地震局提供地电场观测数据。

参考文献
[1]
钱复业, 赵玉林. 地电场短临预报方法研究[J]. 地震, 2005, 25(2): 33-40 (Qian Fuye, Zhao Yulin. Study on Geoelectric Field Method for Short-Term and Impending Earthquake Prediction[J]. Earthquake, 2005, 25(2): 33-40) (0)
[2]
Varotsos P, Alexopoulos K. Physical Properties of the Variations of the Electric Field of the Earth Preceding Earthquakes[J]. Tectonophysics, 1984, 110(1-2): 73-98 DOI:10.1016/0040-1951(84)90059-3 (0)
[3]
Varotsos P, Alexopoulos K. Physical Properties of the Variations of the Electric Field of the Earth Preceding Earthquakes, Ⅱ, Determination of Epicenter and Magnitude[J]. Tectonophysics, 1984, 110(1-2): 99-125 DOI:10.1016/0040-1951(84)90060-X (0)
[4]
席继楼, 陈敏, 董蕾, 等. 长宁MS6.0地震前后地电场秒数据变化特征研究[J]. 地震, 2020, 40(4): 129-143 (Xi Jilou, Chen Min, Dong Lei, et al. Study on the Second Data Variation Characteristics of Geoelectric Field before and after Changning MS6.0 Earthquake[J]. Earthquake, 2020, 40(4): 129-143 DOI:10.3969/j.issn.1003-3246.2020.04.018) (0)
[5]
黄清华, 刘涛. 新岛台地电场的潮汐响应与地震[J]. 地球物理学报, 2006, 49(6): 1 745-1 754 (Huang Qinghua, Liu Tao. Earthquakes and Tide Response of Geoelectric Potential Field at the Niijima Station[J]. Chinese Journal of Geophysics, 2006, 49(6): 1 745-1 754) (0)
[6]
谭大诚, 赵家骝, 席继楼, 等. 潮汐地电场特征及机理研究[J]. 地球物理学报, 2010, 53(3): 544-555 (Tan Dacheng, Zhao Jialiu, Xi Jilou, et al. A Study on Feature and Mechanism of the Tidal Geoelectrical Field[J]. Chinese Journal of Geophysics, 2010, 53(3): 544-555) (0)
[7]
谭大诚, 王兰炜, 赵家骝, 等. 潮汐地电场谐波和各向波形的影响要素[J]. 地球物理学报, 2011, 54(7): 1 842-1 853 (Tan Dacheng, Wang Lanwei, Zhao Jialiu, et al. Influence Factors of Harmonic Waves and Directional Waveforms of Tidal Geoelectrical Field[J]. Chinese Journal of Geophysics, 2011, 54(7): 1 842-1 853) (0)
[8]
谭大诚, 辛建村, 王建军, 等. 大地电场岩体裂隙模型的应用基础与震例解析[J]. 地球物理学报, 2019, 62(2): 558-571 (Tan Dacheng, Xin Jiancun, Wang Jianjun, et al. Application Foundation and Earthquake Case Analysis of the Telluric Field Rock Crack Model[J]. Chinese Journal of Geophysics, 2019, 62(2): 558-571) (0)
[9]
刘长生, 张思萌, 杨维辉, 等. 黑龙江地电场方位角异常与中强地震的关系探讨[J]. 防灾减灾学报, 2020, 36(2): 33-39 (Liu Changsheng, Zhang Simeng, Yang Weihui, et al. The Relationship between Azimuth Anomaly of Geoelectric Field and Moderate Strong Earthquakes in Heilongjiang Province[J]. Journal of Disaster Prevention and Reduction, 2020, 36(2): 33-39) (0)
[10]
钟以璋, 卢造勋, 刘天革. 东北地区地壳结构特征与地震分布[J]. 地震地质, 1979, 1(2): 13-21 (Zhong Yizhang, Lu Zaoxun, Liu Tiange. Characteristics of Crustal Structures in Northeast China and Their Relations to the Distribution of Strong Earthquakes[J]. Seismology and Geology, 1979, 1(2): 13-21) (0)
[11]
周建波, 张兴洲, 马志红, 等. 中国东北地区的构造格局与盆地演化[J]. 石油与天然气地质, 2009, 30(5): 530-538 (Zhou Jianbo, Zhang Xingzhou, Ma Zhihong, et al. Tectonic Framework and Basin Evolution in Northeast China[J]. Oil and Gas Geology, 2009, 30(5): 530-538) (0)
[12]
高立新. 中国东北地区地震活动的动力背景及其时空特征分析[J]. 地震, 2011, 31(1): 41-51 (Gao Lixin. Dynamic Background and Space-Time Characteristics Analysis of Earthquake Activities in Northeast China[J]. Earthquake, 2011, 31(1): 41-51) (0)
[13]
高立新, 戴勇. 现今中国东北地区地球动力学环境[J]. 大地测量与地球动力学, 2020, 40(11): 1 101-1 107 (Gao Lixin, Dai Yong. The Present Geodynamic Environment of Northeast China[J]. Journal of Geodesty and Geodynamics, 2020, 40(11): 1 101-1 107) (0)
[14]
谭大诚, 赵家骝, 席继楼, 等. 青藏高原中强地震前的地电场变异及构成解析[J]. 地球物理学报, 2012, 55(3): 875-885 (Tan Dacheng, Zhao Jialiu, Xi Jilou, et al. The Variation of Waveform and Analysis of Composition for the Geoelectrical Field before Moderate or Strong Earthquake in Qinghai-Tibetan Plateau Regions[J]. Chinese Journal of Geophysics, 2012, 55(3): 875-885) (0)
[15]
李永生, 赵谊, 李继业, 等. 2018年5月28日吉林松原MS5.7地震发震构造分析[J]. 地震学报, 2018, 42(1): 12-23 (Li Yongsheng, Zhao Yi, Li Jiye, et al. Seismogenic Structure Analysis of the MS5.7 Songyuan, Jilin, Earthquake on May 28, 2018[J]. Acta Seismologica Sinica, 2018, 42(1): 12-23) (0)
[16]
李君, 王勤彩, 郑国栋, 等. 2018年5月松原MS5.7地震序列发震断层及应力场特征[J]. 地震学报, 2019, 41(2): 207-218 (Li Jun, Wang Qincai, Zheng Guodong, et al. Characteristics of Seismogenic Faults and Stress Fields of the Songyuan MS5.7 Earthquake Sequence in May 2018[J]. Acta Seismologica Sinica, 2019, 41(2): 207-218) (0)
Study on Variation Characteristics of Natural Electric Field in Songyuan MS5.7 Earthquake on May 28, 2018
LI Mengying1     ZHANG Zhihong1     JIAO Mingruo1     FANG Yuxin2     
1. Liaoning Earthquake Agency, 44 North-Huanghe Street, Shenyang, 110034, China;
2. Shenyang Earthquake Monitoring Center Station, 32 Lane of 208 Dongling Road, Shenyang, 110161, China
Abstract: We study the observational data of the geoelectric field in northeast China and combined with the Songyuan MS5.7 earthquake on May 28, 2018, in Jilin. At the same time, we comparatively analyze the background, imminent and post-earthquake response time series changes of the natural electric field ESP and the telluric field ET around the epicenter. In the study area, we calculate the dominant direction of charge movement by applying the charge movement seepage model of rock fissure water. The results show that the natural electric field ESP in Suihua, Lindian, Zhaodong, and Wangkui show obvious small-amplitude sudden jumps or large-scale jumps before the earthquake, and 45° extrusion or tension anomalies feature occurred in the dominant direction of charge movement within 3 months before the earthquake.
Key words: Songyuan earthquake; natural electric field; telluric field; fracture structure; mechanism