黑龍江地電場方位角異常與中強地震的關(guān)系探討

劉長生，張思萌，楊維輝，康 健，高雙玲

（1.黑龍江省地震局，黑龍江 哈爾濱 150090；2.德都地震臺，黑龍江 德都 164199；3.牡丹江地震臺，黑龍江 牡丹江 157009）

0 引言

地震孕育和發(fā)生必然會導致地下介質(zhì)電性結(jié)構(gòu)發(fā)生改變，從而導致地表大地電場的變化。地電場是地表大地電磁場的一部分，由較穩(wěn)定的自然電場和變化的大地電場構(gòu)成[1]。1984 年Varotsos，Alexopoulos 和Nomicos 提出利用地電場觀測資料提取地震電信號預報地震的方法，即VAN 法。近些年，國內(nèi)地震學者利用多種方法對地震前地電場中可能包含的異常信息進行了分析和提取，毛桐恩[2]提取了河北張北M6.2地震前大地電場矢量方位的顯著變化；馬欽忠[3]基于地電場多極距原理論證了可能與汶川地震有關(guān)的地電場異常前兆，發(fā)現(xiàn)震前脈沖信號并非源于同一地點，不同的應力積累區(qū)可能是異常信號的源區(qū)。譚大城[4]應用大地電場巖體裂隙模型對國內(nèi)典型震例進行分析，發(fā)現(xiàn)震前巖體裂隙優(yōu)勢方位出現(xiàn)偏轉(zhuǎn)和大幅度突跳等結(jié)構(gòu)變異。艾薩·伊斯馬伊力等[5]分析新疆北天山地區(qū)兩次6 級地震前巖體裂隙方位角也發(fā)現(xiàn)類似變化。上述研究表明，地電場可以很好的反應地下介質(zhì)的電性變化，開展孕震過程中大地電場的變化特征研究對地震預測預報具有非常重要的意義。

中國“十五”數(shù)字化臺網(wǎng)建設完成至今，黑龍江省及鄰區(qū)發(fā)生多次M>4.0 中強地震，大多地震發(fā)生在松遼盆地內(nèi)部，省內(nèi)有六個的電場臺布設在其中，本文使用大地電場方位角法分析和討論裂隙方位角異常變化特征與中強地震之間的關(guān)系。

1 地質(zhì)構(gòu)造背景

黑龍江省位于中國東北部，北向和東向皆與俄羅斯接壤，西向與內(nèi)蒙古自治區(qū)、南向與吉林省相連[6]本地區(qū)構(gòu)造運動較為活躍，地震構(gòu)造背景復雜，省內(nèi)及鄰區(qū)多次發(fā)生M4.5～6.0 的中強地震。黑龍江地區(qū)分布有北東東向的嫩江斷裂帶、依蘭—伊通斷裂帶、敦化—密山斷裂帶共三條一級大型斷裂帶；北西向的濱州斷裂和近北南的肇東—扶余斷裂帶等多條次一級斷裂。它們組成了黑龍江地區(qū)活動斷裂構(gòu)造的基本框架（圖1）。

圖1 黑龍江省地電場臺與斷裂分布圖Fig.1 Distribution of electric field stations and faults in the Heilongjiang region

表1 黑龍江省地電場觀測臺站基本參數(shù)

2 觀測資料

黑龍江省地震局在中國地震局進行“十五”數(shù)字觀測網(wǎng)絡改造過程中新增了大地電場觀測項目，該項目在黑龍江省共建立了7 個大地電場觀測臺站，它們分別為綏化地震臺、德都地震臺、通河地震臺、林甸地震臺、肇東地震臺、望奎地震臺、密山地震臺，具體位置如圖1 所示。觀測外線路均采用雙“L”形布極方式，長極距均為300m，短極距為150m，布設采用地埋的方式，電極的埋藏深度為3.0±0.5m，中心點在同一坑內(nèi)埋設兩個電極，從表1 中可以看到各個地震臺站觀測系統(tǒng)的配置參數(shù)。地電臺網(wǎng)于2007 年6 月開始產(chǎn)出數(shù)據(jù)。

3 方法原理和計算過程

3.1 方法原理

在實際場地中，利用巖石裂隙水主體滲流模型原理計算方位角α 時，之前需先進行相關(guān)分析，通常會得出三條ET結(jié)果曲線，如果EW、NE 之間的ET相關(guān)性高，就表明裂隙優(yōu)勢走向介于EW 和NE 之間，即選取該相關(guān)性高的α 曲線作為巖石裂隙狀態(tài)變化的判定標準曲線。但也可能出現(xiàn)NS、NE 或NS、EW 相關(guān)性高，則相應進行選?。槐疚慕Y(jié)合各地電臺站電極布設方式，給出大地電場優(yōu)勢方位角α（北偏東）計算公式如下：

公式（1）為綏化臺滲流方位角計算公式。式中，α 為巖石裂隙優(yōu)勢方位走向；AEW和ANE表示地電場EW 和NE 向觀測數(shù)據(jù)的潮汐諧波振幅值；n 為階數(shù)，表示各臺優(yōu)勢周期個數(shù)，可通過傅里葉變化得到，一般n 計算到n 階即可[4]。

上面公式僅是應用EW、NE 測向數(shù)據(jù)的計算公式。當應用其它測向計算α 時，其表達式需相應調(diào)整。

3.2 計算過程

具體計算過程如下：首先，利用中國地震前兆預處理數(shù)據(jù)庫中每日大地電場三個方向長極距分鐘值采樣數(shù)據(jù)，不同方位之間兩兩計算，共計算出三個相關(guān)系數(shù)，以其中兩個方位之間諧相關(guān)系數(shù)最大的結(jié)果計算方位角；其次，計算出電場中潮汐成分所占比例較高的前10 階振幅之和作為其諧波振幅；最后，基于諧波振幅計算最終方位角。

4 松遼盆地內(nèi)六個地電場波形特征、巖石裂隙優(yōu)勢走向分析

4.1 波形特征分析

譚大城[7]將全國地電場觀測波形分為兩類：一是TGF-A 型，其與固體潮密切相關(guān)的正弦型波形，二是TGF-B 型，其與氣體潮汐有關(guān)的空間Sq電流產(chǎn)生的近梯形波形。認為在沒有或離地表水域稍遠處，地下巖石含水度高、透水性好的區(qū)域Sq電流系在地表的感應引起裂隙水周期性滲流產(chǎn)生大地電場可能是主要原因[8]。黑龍江省七個電場臺站中除了密山臺波形為無序型外，松遼盆地內(nèi)六個臺均表現(xiàn)出較好的TGF-B型，如圖2，表2 所示。

圖2 黑龍江省大地電場臺日變化波形對比（20180501）Fig.2 Comparison of waveform characteristics of platform diurnal variation in Heilongjiang electric field

表2 黑龍江省地電場臺觀測波形分類

基于巖石物理學[9]，巖體內(nèi)部的裂隙發(fā)展過程分為生成、發(fā)育、長大、破碎，而后[4]利用中國大陸地電場臺網(wǎng)中的40 個觀測環(huán)境較好及裝置相對穩(wěn)定的臺站，計算了2017 年上半年的場地巖體裂隙方位角正常變化背景，結(jié)果顯示，80%的場地巖體裂隙處于發(fā)育階段，15%的場地處于長大階段，5%的場地處于破碎階段。本文共計算了黑龍江省六個臺站方位角結(jié)果，其中德都臺、望奎臺、肇東臺、綏化臺和林甸臺是發(fā)育階段、通河臺是長大階段，處于發(fā)育階段的臺站占多數(shù)。

4.2 松遼盆地內(nèi)六個電場臺站巖石裂隙優(yōu)勢走向分析

依據(jù)珧立珣[10]、張鳳鳴[11]等的構(gòu)造應力場研究資料表明，黑龍江地區(qū)的構(gòu)造應力場主壓應力方向σ1約為54°，接近整個東北地區(qū)地殼應力場N58°E，在研究區(qū)內(nèi)不同的臺站場地內(nèi)主壓應力σ1變化明顯，本地區(qū)主要以水平應力為主（圖3）。巖石力學理論認為：兩組共軛剪裂面與主壓應力方向的理論夾角為45°，圖4a 左圖是由σ1方位繪出的兩組共軛剪裂面L1 和L2，走向分別為N9°E、N99°E，這是最大剪應力作用的兩個理論走向，α 是巖石裂隙的大致走向。

圖3 黑龍江及鄰區(qū)斷裂帶的斷層錯動及大地電場臺分布（據(jù)姚立珣，1992）Fig.3 Fault movement and distribution of geoelectric station in Heilongjiang and its adjacent area

圖4 黑龍江省六個地電場臺裂隙優(yōu)勢走向分析Fig.4 Analysis of the preferential crack orientation of the geoelectrical field at six earthquake stations in Heilongjiang province

巖石裂隙實際走向會小于理論計算值[4]，從表3 中可知，望奎、肇東、綏化臺和通河臺的值與巖石學理論得出的可能走向更加接近，其差的形成可能是因為主壓應力與剪裂面的夾角稍小于45°所引起；德都和林甸臺方位角度偏小，與理論走向相差較大，可能與各臺站局部場地、裂隙度、含水度和滲透率等因素差異明顯有關(guān)。

表3 2018 年5 月1 日松遼盆地內(nèi)六臺站巖石裂隙水主體滲流方向

5 典型震例分析

2007 年國家“十五”數(shù)字地震臺網(wǎng)建設以來，黑龍江省及鄰區(qū)發(fā)生多次M≥4.5 地震，選取其中較為典型，同時震中周圍地電臺站數(shù)量較多的地震進行分析，其中震前原始分鐘值數(shù)據(jù)曲線出現(xiàn)異常的唯一典型震例為2009 年5 月10 日安達M4.5 地震，而分析方位角結(jié)果與地震關(guān)系則出現(xiàn)三次典型震例，分別為：2009 年5 月10 日安達M4.5、2013 年6 月20 日內(nèi)蒙古莫力達瓦M5.0，2017 年7 月至2018 年3 月吉林松原M4.0 級叢集地震等共3 組（6 個）地震，如表4 所示。

表4 黑龍江及鄰區(qū)M＞4.0 地震目錄

5.1 安達M4.5 地震綏化臺大地電場分鐘值曲線異常

2009 年5 月10 日安達發(fā)生M4.5 地震，發(fā)震斷層為濱州斷裂。綏化地震臺距離震中128km，從地震發(fā)生前的4 月20 日左右開始，綏化臺大地電場六條曲線均出現(xiàn)明顯的短期異常變化，如圖5 所示。異常結(jié)束的時間基本與發(fā)生地震時間一致。經(jīng)過核實，異常存續(xù)期間臺站周邊并未發(fā)現(xiàn)干擾源。

圖5 2009 年5 月10 日安達M4.5 地震前綏化臺地電場異常變化（2009.03.01—2009.05.31）Fig.5 Abnormal change of electric field at Suihua Seismic Station before Anda M4.5 earthquake on May 10th，2009

5.2 吉林寧江M4.0 叢集地震前望奎方位角異常

2017 年7 月至2018 年3 月間，吉林松原寧江出現(xiàn)M4.0 地震叢集現(xiàn)象，震級分別為M4.9、M4.5、M4.3、M4.2，距 離 該 震 中 約230km 的望奎臺大地電場出現(xiàn)方位角α 突跳的異常現(xiàn)象。由圖6 可以看到2017 年望奎方位角出現(xiàn)明顯的由“長大”階段轉(zhuǎn)變?yōu)椤捌扑椤彪A段，由正常的30°以內(nèi)變化變幅轉(zhuǎn)變?yōu)榇笥?5°變化，變化顯著。

圖6 吉林松原M4.0 叢集地震發(fā)生前望奎大地電場方位角變化特征Fig.6 Azimuth variation characteristics of Heilongjiang Wangkui geoelectric field before Jilin Songyuan M4.0 cluster earthquake

5.3 2009 年5 月10 日安達4.5 級地震

2008 年11 月至2009 年1 月綏化臺電場方位角數(shù)據(jù)出現(xiàn)大幅度改變，如圖7 所示，方位角從之前的穩(wěn)定的30°～50°范圍波動，表現(xiàn)為“發(fā)育”階段，變幅為20°左右，逐漸演變成為幅度為45°上下波動，即從“長大”階段轉(zhuǎn)變成“破碎”階段，在此期間發(fā)生了震中距120km 的安達4.5 地震，而后方位角恢復至正常狀態(tài)。

圖7 2009 年黑龍江安達M4.5 地震前綏化大地電場方位角變化特征Fig.7 Azimuth variation characteristics of Suihua geoelectric field before Heilongjiang Anda M4.5 earthquake in 2009

6 結(jié)論

通過對黑龍江省七個大地電場臺分鐘值數(shù)據(jù)進行梳理和方位角計算結(jié)果與中強地震的對應關(guān)系進行分析，可以得出以下結(jié)論：

（1）黑龍江地區(qū)地震臺站地電場和地磁場日變化具有同步性和相關(guān)性。地電場觀測的日變化波形形態(tài)屬于“峰-谷-峰”型，其中綏化、德都、通河屬于峰-谷-峰型；林甸屬于峰-谷型；望奎北南向和北西屬于雙峰-雙谷型，東西向是雙峰-單谷型；所有臺站觀測波形都屬于TGF-B 型，波形形態(tài)受空間Sq電流系產(chǎn)生的感生電流控制。

（2）與前人對本地區(qū)構(gòu)造應力場方位走向進行對比，得到望奎臺、肇東臺、綏化臺和通河臺等臺的值更符合巖石學理論的可能走向，其差可能是因為主壓應力與剪裂面的夾角稍小于45°所引起；德都和林甸臺方位角度偏小，與理論走向相差較大，可能與各臺站局部場地、裂隙度、含水度和滲透率等因素差異明顯有關(guān)。

（3）省內(nèi)七個大地電場臺原始數(shù)據(jù)受干擾（如電極引線破損、外界環(huán)境、降雨等）變化幅度超100mv/km 以上時，方位角數(shù)值會受到顯著影響，幅度在30mv/km 以下時，方位角基本不受干擾。

（4）通過對大地電場滲流方位角結(jié)果與歷史中強地震之間分析，發(fā)現(xiàn)黑龍江及鄰區(qū)顯著地震前會出現(xiàn)巖體裂隙剪裂變現(xiàn)象，兩者之間存在較好的對應關(guān)系。

總體上，對黑龍江及鄰區(qū)2 組（5 次）典型震例進行分析發(fā)現(xiàn)，在強震震中200km 范圍內(nèi)，巖體裂隙處于“長大”或“發(fā)育”好的場地較集中，并且這些場地在中強震前更容易發(fā)生巖體剪裂和大幅度突跳或集中等結(jié)構(gòu)變異的變化特征。

致謝：甘肅省地震局譚大城研究員提供方位角計算軟件和給予的指導和幫助，黑龍江省地震局監(jiān)測中心提供的觀測數(shù)據(jù)。