2018-05-28松原MS5.7地震地電場(chǎng)變化特征研究

李夢(mèng)瑩 張志宏 焦明若 方禹心

1 遼寧省地震局，沈陽(yáng)市黃河北大街44號(hào)， 110034 2 沈陽(yáng)地震監(jiān)測(cè)中心站，沈陽(yáng)市東陵路208巷32號(hào)，110161

根據(jù)中國(guó)地震臺(tái)網(wǎng)(http:∥www.ceic.ac.cn)正式測(cè)定，2018-05-28 01:50吉林松原發(fā)生MS5.7地震。本次地震為2013-11-23吉林前郭MS5.8地震后的首次較大地震，震中區(qū)也是中國(guó)東部1975-02-04遼寧海城MS7.3地震和1976-07-28河北唐山MS7.8地震后地震活動(dòng)性最強(qiáng)的地區(qū)，對(duì)研究地震前地球物理場(chǎng)中短期異常有著重要意義。

我國(guó)地電場(chǎng)觀測(cè)與地震預(yù)測(cè)研究始于1966年邢臺(tái)MS7.2地震后，曾記錄到不少震前異常變化[1]，Varotsos等[2-3]用多極距觀測(cè)法解決了地電場(chǎng)觀測(cè)系統(tǒng)噪聲問(wèn)題，自此國(guó)內(nèi)地電場(chǎng)地震預(yù)測(cè)預(yù)報(bào)方法得到快速發(fā)展[4-8]。本文使用2018-05-28吉林松原MS5.7地震震前1 a以上及震后6個(gè)月震中500 km范圍內(nèi)13個(gè)地電場(chǎng)臺(tái)站分鐘值數(shù)據(jù)，分析地震前后自然電場(chǎng)變化，得到震前自然電場(chǎng)ESP異常特征，并結(jié)合巖體裂隙水電荷滲流模型計(jì)算異常臺(tái)站地電場(chǎng)優(yōu)勢(shì)方位角a，討論地電場(chǎng)優(yōu)勢(shì)方位角在地震預(yù)測(cè)方面的優(yōu)勢(shì)。

1 原理及方法介紹

1.1 自然電場(chǎng)和大地電場(chǎng)

地電場(chǎng)是指固體地球內(nèi)部和外部的各種非人工電流系與地球介質(zhì)相互作用，在地球表面產(chǎn)生的電場(chǎng)，包括大地電場(chǎng)ET和自然電場(chǎng)ESP。ET與磁層和電離層中電流系的運(yùn)動(dòng)有關(guān)，ESP則是地殼內(nèi)部各類物理化學(xué)作用引起的正負(fù)電荷分離產(chǎn)生的地電場(chǎng)[1]。ET和ESP通過(guò)測(cè)量地表每個(gè)測(cè)向兩點(diǎn)之間的電位差，得到其隨時(shí)間的變化，如圖1所示，通過(guò)埋設(shè)在地表下的測(cè)量電極A和B，以一定時(shí)間間隔測(cè)量2個(gè)電極的自然電位差序列(VAB)i(i=N-1，N-2，…，2，1，0)，其中i=0表示當(dāng)前時(shí)間，i=1表示當(dāng)前時(shí)間之前第1個(gè)時(shí)間點(diǎn)，i=N-1表示當(dāng)前時(shí)間之前第N-1個(gè)時(shí)間點(diǎn)。按式(1)計(jì)算得出的自然電位差序列(VAB)i算術(shù)平均值，即為當(dāng)前時(shí)間的自然電場(chǎng)ESP；按式(2)計(jì)算此方向上的大地電場(chǎng)ET分量值，其中負(fù)號(hào)表示大地電場(chǎng)的方向由高電位指向低電位。

(1)

式中，(VAB)i(i=N-1，N-2，…，2，1，0)為電極A和B之間的自然電位差值序列，單位mV；N為參與計(jì)算的自然電位差值序列數(shù)據(jù)總數(shù)，不少于24 h之內(nèi)的所有測(cè)量數(shù)據(jù)個(gè)數(shù)。

(2)

圖1 地電場(chǎng)觀測(cè)基本原理示意圖Fig.1 Schematic diagram of basic principle of geoelectric field observation

大地電場(chǎng)矢量由2個(gè)正交方向的大地電場(chǎng)分量值合成，2個(gè)正交方向分別為x軸(東向分量)和y軸(北向分量)，由2個(gè)大地電場(chǎng)分量值Ex和Ey按式(3)和式(4)計(jì)算大地電場(chǎng)的幅值E和方位角a：

(3)

(4)

式中，Ex為大地電場(chǎng)東向分量，Ey為大地電場(chǎng)北向分量，單位均為mV/km。

地電場(chǎng)的觀測(cè)數(shù)據(jù)中，ESP為自然電場(chǎng)，ET為大地電場(chǎng)，Er為干擾成分，則地電場(chǎng)E可表達(dá)為：

E=ESP+ET+Er

(5)

設(shè)地電場(chǎng)觀測(cè)分鐘值數(shù)據(jù)為Ei，日均值計(jì)算自然電場(chǎng)ESP的表達(dá)式可簡(jiǎn)化為：

(6)

應(yīng)用式(6)能夠基本消除大地電場(chǎng)ET和干擾成分Er。

1.2 地電場(chǎng)優(yōu)勢(shì)方位角a

中國(guó)學(xué)者基于大地電場(chǎng)ET日變波形的時(shí)域和頻域特征，提出了大地電場(chǎng)ET日變波源于電離層Sq電流和潮汐力學(xué)說(shuō)，即大地電場(chǎng)的潮汐機(jī)理[5-6]。譚大誠(chéng)等[8]在潮汐機(jī)理的基礎(chǔ)上建立了大地電場(chǎng)ET的巖體裂隙水(電荷)滲流(移動(dòng))模型，解釋了大地電場(chǎng)ET強(qiáng)度和方向的變化是地下介質(zhì)應(yīng)力的積累改變巖體裂隙結(jié)構(gòu)的過(guò)程。圖2(b)為裂隙水電荷日周期移動(dòng)的理想模型，巖體結(jié)構(gòu)發(fā)育程度決定了實(shí)際的電荷移動(dòng)方向。用a表示場(chǎng)地巖體裂隙水電荷移動(dòng)的主滲流方位，也可稱之為地電場(chǎng)優(yōu)勢(shì)方位角a。

圖2 ET巖體裂隙水(電荷)滲流(移動(dòng))模型Fig.2 ET model of fissure water in rock with seepage(movement)or charges

當(dāng)臺(tái)站地電場(chǎng)NS和NW測(cè)向之間的相關(guān)性高時(shí)，地電場(chǎng)ET的優(yōu)勢(shì)方位角a(北偏東)計(jì)算公式如下：

(7)

式中，ANW(i)、ANS(i)分別為NW、NS測(cè)向第i階潮汐諧波振幅，Ai的計(jì)算公式如下：

(8)

對(duì)于一個(gè)數(shù)據(jù)序列yt(時(shí)間序列總數(shù)為n)，數(shù)學(xué)上可表示為：

(9)

2 臺(tái)站地質(zhì)構(gòu)造概況

本文選取2018年松原MS5.7地震震中500 km范圍內(nèi)采樣時(shí)間間隔為1 min、測(cè)量頻段為0～0.005 Hz、“雙L”6測(cè)道長(zhǎng)短極距布設(shè)、不極化電極的13套ZD9A-Ⅱ型儀器和2套GEF-2型儀器地電場(chǎng)臺(tái)站[9](表1)，圖3為臺(tái)站所處現(xiàn)今地質(zhì)構(gòu)造區(qū)域。鐘以璋等[10]通過(guò)研究衛(wèi)星影像數(shù)據(jù)和地球物理場(chǎng)觀測(cè)數(shù)據(jù)發(fā)現(xiàn)，研究區(qū)地殼構(gòu)造由斷塊組成，地震震中沿?cái)嗔褞Х植迹艺共挤较蚺c區(qū)域內(nèi)主要斷裂帶的走向大體一致。研究區(qū)地屬松遼盆地，該盆地為中國(guó)東北部的一個(gè)大型中-新生代NE向沉積盆地，地跨中國(guó)東北三省和內(nèi)蒙古自治區(qū)[11]。21世紀(jì)以來(lái)，松遼盆地及周邊地區(qū)地震活動(dòng)明顯增強(qiáng)，2002年至今已連續(xù)發(fā)生5級(jí)以上天然地震16次，高立新等[12-13]認(rèn)為2011年日本MW9.0地震對(duì)中國(guó)東北地區(qū)地球動(dòng)力學(xué)環(huán)境影響顯著，地殼運(yùn)動(dòng)狀態(tài)由擠壓轉(zhuǎn)為拉張。2011年以來(lái)中國(guó)東北地區(qū)相繼發(fā)生遼寧燈塔MS5.1地震、內(nèi)蒙古科爾沁M(jìn)S5.3地震、吉林前郭MS5.8震群等7次5.0級(jí)以上淺源地震。

圖3 2018年松原MS5.7地震震中及500 km內(nèi)地電場(chǎng)臺(tái)站分布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

表1 研究區(qū)地電臺(tái)站概況

3 分析研究

3.1 自然電場(chǎng)變化

根據(jù)場(chǎng)源的不同，可將地電場(chǎng)E分為大地電場(chǎng)ET和自然電場(chǎng)ESP，大地電場(chǎng)ET的起源是空間Sq電流系和固體潮，自然電場(chǎng)ESP則源于地下介質(zhì)的物理、化學(xué)作用。自然電場(chǎng)ESP包括地下金屬礦體與周圍溶液經(jīng)過(guò)氧化還原反應(yīng)的氧化還原場(chǎng)、地下電介質(zhì)在巖體裂隙中流動(dòng)的過(guò)濾電場(chǎng)及由2種電解質(zhì)濃度差異引起的接觸-擴(kuò)散電場(chǎng)。膨脹和壓力實(shí)驗(yàn)表明，巖體在應(yīng)力變化作用下，能夠激發(fā)自然電場(chǎng)的變化[1]。因此，監(jiān)測(cè)孕震過(guò)程的應(yīng)力變化可能觀測(cè)到自然電場(chǎng)短期異常[14]。

圖4為2017～2018年本文所選臺(tái)站中自然電場(chǎng)ESP變化較大，且觀測(cè)系統(tǒng)和場(chǎng)地較好的4個(gè)場(chǎng)地?cái)?shù)據(jù)。圖4(a)表明，2017-03林甸臺(tái)ESP的EW、NW測(cè)向開(kāi)始出現(xiàn)小幅度異常躍變，持續(xù)至 7 月初躍變幅度峰值約為400 mV·km-1，2018-03開(kāi)始第2次異常躍變，且異常測(cè)向變幅方向相反、幅度差較??；5月初ESP的EW、NW測(cè)向恢復(fù)背景值，松原MS5.7地震后ESP沒(méi)有觀測(cè)到同震響應(yīng)。圖4(b)表明，2018-02綏化臺(tái)ESP三測(cè)向出現(xiàn)同步異常，躍變幅度峰值較小，小于 50 mV·km-1。圖4(c)表明，2018-03初望奎臺(tái)ESP的NS、N45°W測(cè)向出現(xiàn)躍變異常，躍變峰值約為500 mV·km-1，且兩異常測(cè)向的變幅同步性較高，地震發(fā)生在ESP異常幅度恢復(fù)期，震后 6 個(gè)月ESP恢復(fù)正常。圖4(d)表明，2018-03初肇東臺(tái)ESP的EW、NS測(cè)向觀測(cè)值突增，異常起始幅度已達(dá)峰值500 mV·km-1。

由表2可知，綏化臺(tái)ESP的異常幅度相對(duì)于背景值變化并不顯著，異常信度最低；肇東臺(tái)、望奎臺(tái)和林甸臺(tái)異常變化顯著，異常信度較高，其中林甸臺(tái)的異常幅度相對(duì)較小，異常起始時(shí)間較早。從圖4也可看出，肇東臺(tái)、望奎臺(tái)和林甸臺(tái)ESP異常測(cè)向表現(xiàn)為兩道，其中望奎臺(tái)和林甸臺(tái)同屬富裕-明水?dāng)嗔?圖3)，而肇東臺(tái)ESP異常變化速率顯著，短時(shí)期內(nèi)達(dá)到了峰值，發(fā)震構(gòu)造研究結(jié)果[15-16]也支持肇東臺(tái)的異常特征。

表2 松原MS5.7 地震地電場(chǎng)異常臺(tái)站數(shù)據(jù)統(tǒng)計(jì)

3.2 大地電場(chǎng)巖體裂隙優(yōu)勢(shì)方位

圖5為基于大地電場(chǎng)巖體裂隙水(電荷)滲流(移動(dòng))模型計(jì)算的4個(gè)異常場(chǎng)地巖體裂隙水主滲流方位。圖5(a)為林甸臺(tái)大地電場(chǎng)優(yōu)勢(shì)方位角α的變化散點(diǎn)分布，可以看出，2017-03林甸臺(tái)方位角α的突變范圍出現(xiàn)了大幅度的發(fā)散，Δα約為90°，對(duì)比圖4(a)中林甸臺(tái)自然電場(chǎng)ESP的變化開(kāi)始時(shí)間發(fā)現(xiàn)，兩者的變化初始時(shí)間明顯對(duì)應(yīng)。由此可見(jiàn)，2017-03林甸臺(tái)自然電場(chǎng)ESP大幅變化的物理機(jī)理可能是場(chǎng)地巖體裂隙結(jié)構(gòu)在應(yīng)力加卸載作用下發(fā)生了較大改變，表現(xiàn)為間斷的巖體裂隙偏轉(zhuǎn)或剪裂現(xiàn)象。圖5(c)和5(d)為望奎臺(tái)和肇東臺(tái)大地電場(chǎng)優(yōu)勢(shì)方位角α的變化散點(diǎn)分布，可以看出，兩者方位角α的突變范圍呈大幅度的擠壓狀態(tài)，Δα在45°左右。對(duì)比圖4(c)和4(d)可知，望奎臺(tái)自然電場(chǎng)ESP同步發(fā)生異常變化，松原MS5.7地震后望奎臺(tái)大地電場(chǎng)優(yōu)勢(shì)方位角α的擠壓狀態(tài)持續(xù)至9月，后恢復(fù)發(fā)散狀態(tài)，而肇東臺(tái)震前由擠壓恢復(fù)至擴(kuò)張。

圖4 異常臺(tái)站在地震前后地電場(chǎng)變化(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)

圖5 異常臺(tái)站大地電場(chǎng)優(yōu)勢(shì)方位角α的變化散點(diǎn)分布Fig.5 Variation scatter diagram of the dominant azimuth angle α of the telluric field at anomalous stations

4 結(jié) 語(yǔ)

李永生等[15]的全波形矩張量反演結(jié)果表明，扶余-肇東斷裂東段為2018年松原MS5.7地震的發(fā)震構(gòu)造；李君等[16]利用雙差定位法對(duì)地震序列進(jìn)行重新定位，并利用MSATSI軟件反演得到松原地區(qū)應(yīng)力場(chǎng)，結(jié)果表明，松原地區(qū)應(yīng)力場(chǎng)既受太平洋板塊對(duì)北東亞板塊的西向俯沖，又受區(qū)域構(gòu)造運(yùn)動(dòng)的控制，在中國(guó)東北地區(qū)產(chǎn)生了近 EW向的主壓應(yīng)力。因此，自然電場(chǎng)變化量ΔESP最大的臺(tái)站應(yīng)位于扶余-肇東斷裂，且EW測(cè)向異常顯著，這就解釋了肇東臺(tái)ESP的EW測(cè)向異常顯著和巖體裂隙優(yōu)勢(shì)方位角出現(xiàn)高達(dá)45°擠壓的異常特征。震中周邊場(chǎng)地的孕震過(guò)程是應(yīng)力不斷加卸載的過(guò)程，巖體裂隙結(jié)構(gòu)會(huì)因應(yīng)力加卸載的變化而變化。但實(shí)際結(jié)果是，場(chǎng)地巖體結(jié)構(gòu)的差異導(dǎo)致不同場(chǎng)地觀測(cè)到的地電場(chǎng)ESP變化差異較大，這可能是異常臺(tái)站只有兩道顯著異常的原因。

本文通過(guò)研究近年來(lái)我國(guó)東北地區(qū)地震活動(dòng)最為活躍的吉林松原地區(qū)地電場(chǎng)觀測(cè)資料，結(jié)合2018-05-28松原MS5.7地震，分別對(duì)地震前后自然電場(chǎng)ESP和大地電場(chǎng)ET的背景、臨震、震后響應(yīng)時(shí)序變化進(jìn)行對(duì)比分析，并應(yīng)用巖體裂隙水電荷移動(dòng)滲流模型計(jì)算了研究區(qū)所有地電場(chǎng)臺(tái)站的電荷移動(dòng)優(yōu)勢(shì)方位，得到如下結(jié)論：

1)多極距、多方向的地電場(chǎng)布極方式可有效識(shí)別干擾，但不易提取前兆異常信息；

2)2018-05-28松原 MS5.7 地震前后，綏化臺(tái)、林甸臺(tái)、肇東臺(tái)、望奎臺(tái)自然電場(chǎng)ESP的異常起始時(shí)間存在準(zhǔn)同步特征；

3)基于巖體裂隙水電荷滲流模型及大地電場(chǎng)優(yōu)勢(shì)方位的結(jié)果，可得出巖體結(jié)構(gòu)發(fā)育狀況；

4)松原MS5.7 地震前后異常臺(tái)站自然電場(chǎng)ESP的大幅度變異現(xiàn)象，可能是松原地區(qū)的區(qū)域應(yīng)力場(chǎng)受到構(gòu)造板塊應(yīng)力的影響,望奎臺(tái)及肇東臺(tái)的大地電場(chǎng)優(yōu)勢(shì)方位角α呈擠壓異常。

利用巖體裂隙水電荷滲流模型計(jì)算得到的電荷移動(dòng)方向可有效捕捉場(chǎng)地孕震巖體結(jié)構(gòu)的微變化，由于孕震過(guò)程復(fù)雜，對(duì)該研究方法的適用可能存在區(qū)域性，后續(xù)需要更多震例支持。

致謝:感謝甘肅省地震局譚大誠(chéng)研究員提供地電分析程序，黑龍江省地震局和吉林省地震局提供地電場(chǎng)觀測(cè)數(shù)據(jù)。