張國苓 喬子云 賈立峰 梁紅杰 韓和平 王艷榮

?

河北陽原臺地電阻率變化分析1

張國苓1）喬子云1）賈立峰1）梁紅杰1）韓和平2）王艷榮3）

1）河北省地震局，石家莊050021?2）陽原地電臺，張家口075800?3）三河地震臺，三河065200

河北陽原臺地電阻率NE測向2014年以來出現(xiàn)年變幅變小的年變畸變現(xiàn)象，經(jīng)核實(shí)，這時的觀測系統(tǒng)工作正常，而降雨量與地電阻率的異常變化具有一定的相關(guān)性，但通過褶積濾波法去除降雨量的影響后，這種異常變化依舊存在。同時，陽原臺NE測向的地電阻率與自然電位在2014年的變化具有同步性，也進(jìn)一步驗(yàn)證了異常變化的可靠性。陽原臺地電阻率觀測位于晉冀蒙交界地區(qū)，對地殼的應(yīng)力-應(yīng)變響應(yīng)比較敏感，其地電阻率異常變化與1989年大同地震和1998年張北地震有較好的對應(yīng)關(guān)系，并顯示出持續(xù)時間長、幅度大和恢復(fù)慢的特點(diǎn)。因此，陽原臺地電阻率異常變化的可信度較高，可為監(jiān)測晉冀蒙交界地區(qū)的地震活動提供參考。

地電阻率 陽原 異常分析

引言

在孕震過程中震中區(qū)及附近應(yīng)力場的改變，將會引起巖石的裂隙度及裂隙的連通性發(fā)生變化，此時地電阻率也會隨之發(fā)生改變。野外原地地電阻率的實(shí)驗(yàn)結(jié)果顯示，在低應(yīng)力下地電阻率會隨著巖石介質(zhì)體積被壓縮導(dǎo)致地電阻率的下降，相反，巖石介質(zhì)體積膨脹時地電阻率會上升（趙玉林等，1983）。近40年來地電阻率觀測在地震監(jiān)測預(yù)報(bào)工作中已經(jīng)積累了大量的實(shí)際案例，例如：錢復(fù)業(yè)等（1982；1987）和錢家棟（1993）發(fā)現(xiàn)，在大震發(fā)生前，震中區(qū)及其附近地區(qū)的地電阻率一般會表現(xiàn)出2—3年、幅度為百分之幾到百分之十幾的趨勢性變化，其中多以下降型負(fù)異常為主；同時趙玉林等（1978）也發(fā)現(xiàn)，在1976年唐山地震發(fā)生前，昌黎臺、寶坻臺、唐山臺等的地電阻率也出現(xiàn)了下降變化，另外，在1976年松潘地震發(fā)生前，武都臺的地電阻率也出現(xiàn)了下降變化；還有張學(xué)民等（2009）也發(fā)現(xiàn)，在2008年汶川地震發(fā)生前，郫縣臺的地電阻率在單測向下降的幅度達(dá)到7.2%，下降持續(xù)時間在2年左右，并在震前5個月出現(xiàn)異常轉(zhuǎn)折，而地震就發(fā)生在地電阻率恢復(fù)上升的階段。

本文基于陽原臺地電阻率觀測呈現(xiàn)出的有規(guī)則的冬高、夏低年周期變化特征，以及2014年2月出現(xiàn)的NE測向年變幅減?。曜兓儸F(xiàn)象）和NW測向的年變形態(tài)及幅度與往年相當(dāng)?shù)馁Y料（圖1），并結(jié)合觀測系統(tǒng)檢測和背景場環(huán)境以及各年度日均值數(shù)據(jù)的相關(guān)性分析結(jié)果，重點(diǎn)分析了降雨量對陽原臺地電阻率NE測向的趨勢性影響。同時還進(jìn)一步驗(yàn)證了陽原臺地電阻率NE測向出現(xiàn)的異常變化（年變畸變現(xiàn)象）的可靠性。

1 陽原臺地電阻率觀測概況

陽原臺位于陽原縣西城鎮(zhèn)北關(guān)村，處于晉冀蒙三省交界處，屬陰山—燕山和太行山兩大構(gòu)造交匯處。本地區(qū)主要活動斷裂為北東東向延展的南山山前斷裂（又稱六棱山山前斷裂）和北山山前斷裂（又稱八棱山山前斷裂），兩條斷裂均為傾向盆地中心的正斷層。由于陽原臺的地電阻率測點(diǎn)位于斷裂帶和沖積扇上，因此有利于監(jiān)測來自地殼深部的信息。

陽原臺地電阻率觀測有N22°E和N68°W兩個測向，使用的儀器為ZD8B型地電儀，供電極距和測量極距分別為1.2km和0.2km。測區(qū)的電測深曲線為QQ型，地下電性結(jié)構(gòu)可分為三層：第一層為第四紀(jì)黃土層，夾有礫石，電阻率為125.19—130.61Ω·m，層厚為18.0—19.8m；第二層為泥層，電阻率為11.29—15.27Ω·m，層厚為51.2—101.9m；第三層為粘土層，電阻率為4.12—4.21Ω·m，層厚為113.6—145.3m。第三層以下為砂質(zhì)灰?guī)r層，電阻率為2.75—3.24Ω·m。圖2為陽原臺巖性柱狀圖、電測深曲線和區(qū)域地質(zhì)構(gòu)造圖。由于測區(qū)內(nèi)有大秦鐵路穿過，同時還有小型液化氣站位于測區(qū)中央，所以對數(shù)據(jù)觀測有一定的影響。

2 陽原臺地電阻率變化分析

2.1 觀測系統(tǒng)檢測

根據(jù)《地震及前兆數(shù)字觀測技術(shù)規(guī)范（電場觀測）》（中國地震局科技監(jiān)測司，2001）的要求，陽原臺地電阻率觀測儀器和系統(tǒng)每個季度定期檢測一次，其結(jié)果顯示測量儀器校檢、穩(wěn)流電源性能、供電線、測量線絕緣等均符合規(guī)范的要求，同時各電極接地電阻穩(wěn)定。2014年6月16日，由河北省地震局預(yù)測研究中心電磁組和中國地震臺網(wǎng)中心組成的工作組，到陽原臺進(jìn)行了現(xiàn)場考察，給出了電極連接良好、外線路無破損、線桿及瓷瓶無損毀的測評意見。同時，工作組還用ZC-8型接地電阻儀測量了各電極的接地電阻，并利用ZC-7型兆歐表測量了各測量線的絕緣電阻，其結(jié)果與2013年、2014年的各季度校驗(yàn)基本一致，具體如表1所示。從表中可以看出，供電極接地電阻小于或等于30Ω，測量線路對地絕緣電阻大于或等于5MΩ，均符合規(guī)范要求。

表1 部分電極的接地電阻結(jié)果

2.2 背景場環(huán)境

工作組對陽原臺周邊一定范圍內(nèi)新建（或改建）工程以及其它環(huán)境變化情況進(jìn)行了專門的巡查和走訪。在測區(qū)內(nèi)，有一個1993年9月修建的小型液化氣站；1987年修建的大秦鐵路從測區(qū)通過，鐵路距NW道的東供電極（B3）只有20m；同時NW道的西測量電極（M3）附近有一個小型養(yǎng)羊場，其建于2014年4月位于NW測線的北邊，距M3 電極約100m；另外NW道的東測量電極（N3）附近有兩個蘑菇大棚，其建于2010年位于NW測線的南邊，距N3電極約200m。由于養(yǎng)羊場和蘑菇大棚采用鋼絲網(wǎng)接地，所以形成了一個低阻體，在雨季時其對NW道地電阻率有一定影響，但對NE道干擾不大，如圖3所示。根據(jù)筆者核實(shí)，周邊建筑物造成的環(huán)境影響在時間上與2014年陽原臺地電阻率NE測向出現(xiàn)的年變畸變并不同步。

2.3 地電阻率與降雨的關(guān)系

陽原臺的地電阻率觀測資料顯示，每逢雨季地電阻率就會出現(xiàn)突變，且N68°W測道比N22°E測道受降雨影響更明顯。筆者統(tǒng)計(jì)了2003年以來陽原地區(qū)的降雨量后發(fā)現(xiàn)，陽原的日地電阻率受降雨的影響比較明顯，降雨后地電阻率測值開始下降，而且表現(xiàn)為降雨越大，測值下降的幅度越大；降雨停止后，隨著地面慢慢變干，測值又逐漸恢復(fù)到變化前的狀態(tài)。陽原臺地電阻率在夏季呈現(xiàn)“山羊胡子”狀變化，其主要原因是受降雨影響（見圖1）。從月降雨量和地電阻率的關(guān)系（圖4）可以看出，每年的6—9月為降雨量較大的季節(jié)，地電阻率較低，存在冬高、夏低的年變規(guī)律。筆者通過統(tǒng)計(jì)地電阻率月均值與月降雨總量的相關(guān)關(guān)系后發(fā)現(xiàn)，其相關(guān)系數(shù)達(dá)0.7左右，說明地電阻率與降雨量有較好的相關(guān)性。另外從機(jī)理上看，這可能也與大秦鐵路穿過測區(qū)有關(guān)（圖5），降雨時地下電極與鐵軌形成通路，供電后就會有部分電流流向鐵軌，從而造成測值降低。

地電阻率觀測的是一定體積內(nèi)巖層電性結(jié)構(gòu)隨時間的變化，而區(qū)域內(nèi)表層電性可能受測區(qū)內(nèi)環(huán)境、氣象和水文條件變化的影響（劉允秀等，1999）。降雨對地電阻率的影響比較復(fù)雜，它不僅有一定的即時效應(yīng)，還存在一定的時間滯后效應(yīng)，利用褶積濾波法可以較好地去除降雨對地電阻率的影響（張學(xué)民等，1996）。

根據(jù)褶積濾波法原理，即把土層看成一個系統(tǒng)，降雨量C()作為輸入信號，當(dāng)降雨量對s的影響函數(shù)為()、地電阻率的趨勢變化為()=0+1時，則地電阻率，如果將降雨的即時及滯后效應(yīng)考慮進(jìn)去的話，則有：

例如選取滯后影響時間為5個月時，則上式可表示為：

張學(xué)民等（1996）曾利用上述方法對昌黎臺的地電阻率進(jìn)行了降雨量校正，計(jì)算結(jié)果表明其效果要比使用一般的回歸方法好。筆者利用Matlab軟件對2003—2014年陽原臺降雨量對NE向地電阻率的影響值進(jìn)行了多元回歸計(jì)算，其結(jié)果如圖6所示。從圖中可以看出，陽原臺NE向地電阻率年變主要是由降雨量引起的，當(dāng)去掉降雨量的影響后，地電阻率相對穩(wěn)定一些，年變不再清晰，但是還是存在趨勢性變化，也就是說這種趨勢性的地電阻率年變畸變下降變化，不是僅由降雨量造成的。

2.4 地電阻率與自然電位的對比

圖7是陽原臺NE向自然電位變化曲線，從圖中可以看出，NE測向的自然電位在2014年的幅值要比往年的高，出現(xiàn)了年變畸變，而NW測向的自然電位變化不明顯，這表明在時間和測向上自然電位與地電阻率有很好的對應(yīng)關(guān)系。尹相好等（2010）在總結(jié)了前人的研究后認(rèn)為，當(dāng)自然電位和地電阻率同時發(fā)生異常變化時，可提高異常的可信度。同時異常出現(xiàn)的時間與地震的發(fā)生也有比較好的對應(yīng)關(guān)系，例如：1998年12月9日山東莒縣陵陽L4.1級地震和2001年4月14日山東五蓮于里L(fēng)3.8級地震發(fā)生前，在山東莒縣深井的地電阻率和自然電位同時出現(xiàn)了異常變化。從物理成因上來說，在地震發(fā)生前由于應(yīng)力的積累，地下介質(zhì)會發(fā)生變化，裂隙的膨脹和發(fā)展以及地下流體的共同作用，可能會造成自然電位和地電阻率同時發(fā)生變化。

2.5 震例總結(jié)

張國民等（1995）的研究表明，在中國大陸發(fā)生的地震往往具有成組活動的特征，即地震的發(fā)生在時間和空間上具有相對集中的特點(diǎn)。同時，大量的地震案例也證實(shí)，地電阻率的異常變化與成組地震的發(fā)生有一定的對應(yīng)性，這種現(xiàn)象在地電阻率臺站密集的華北地區(qū)表現(xiàn)的尤為顯著。在1989年大同地震和1998年張北成組地震發(fā)生時，就有多個臺站記錄到了長時間的地電阻率異常，并且異常的時空變化對成組強(qiáng)震活動的遷移有一定的指示意義。這種異常多呈現(xiàn)出比單個地震更為復(fù)雜的形態(tài)特征，異常往往持續(xù)到成組地震中最后一個強(qiáng)震發(fā)生后才結(jié)束（王志賢等，1994；1996a；1996b；汪志亮等，2000；2002；張學(xué)民等，2003）?？傮w表現(xiàn)為趨勢性變化背景下的反復(fù)升降，并呈現(xiàn)出異常持續(xù)時間長、幅度大和恢復(fù)慢的特點(diǎn)（嚴(yán)玲琴等，2013）。

自1977年至今，陽原臺地電阻率觀測已經(jīng)連續(xù)進(jìn)行了37年，在此期間，距離臺站300km范圍內(nèi)發(fā)生的S5.0級以上地震共有8次（去除余震）。在張學(xué)民等（2003）總結(jié)震例的基礎(chǔ)上，筆者對陽原臺地電阻率觀測數(shù)據(jù)進(jìn)行了重新分析后發(fā)現(xiàn)，在幾次地震發(fā)生前陽原臺的地電阻率都曾出現(xiàn)過顯著的年變畸變異常變化，這種地電阻率的異常變化對距離陽原臺300km以內(nèi)的地震具有很好的映震能力，尤其是對晉冀蒙交界地區(qū)發(fā)生的中強(qiáng)地震，具體映震能力如表2所示。

表2 陽原臺地電阻率異常變化對地震的映震能力

3 結(jié)論與討論

通過對觀測系統(tǒng)的檢測，確認(rèn)陽原臺地電觀測儀器、外線路、電極工作狀態(tài)正常。測區(qū)背景場環(huán)境調(diào)查和輔助資料分析表明，降雨量和地下水位引起地電阻率數(shù)據(jù)異常的可能性較小，陽原臺的地電阻率變化反映了真實(shí)的應(yīng)力變化，可視為地震前兆異常。同時，通過地電阻率與自然電位的對比分析后發(fā)現(xiàn)，陽原臺地電阻率的NE測向與自然電位在2014年的變化具有同步性，驗(yàn)證了地電阻率異常的可靠性；陽原臺地電阻率的異常變化多以年變畸變?yōu)橹鳎?014年以來的異常形態(tài)與1998年大同地震和2006年文安地震發(fā)生前的異常形態(tài)較為類似。綜合分析認(rèn)為，陽原臺地電阻率異常變化的可信度較高，可視為地震前兆異常，類型為背景異常。

劉允秀，陳華靜，程瑞年等，1999．地電阻率與地下水位、大氣降水關(guān)系研究．中國地震，15（2）：184—189．

錢復(fù)業(yè)，趙玉林，于謀明等，1982．地震前地電阻率的異常變化．中國科學(xué)（B輯），（9）：831—839．

錢復(fù)業(yè)，趙玉林，許同春，1987．地電阻率季節(jié)干擾變化分析．地震學(xué)報(bào)，9（3）：289—302．

錢家棟，1993．與大震孕育過程有關(guān)的地電阻率變化研究．中國地震，9（4）：341—350．

汪志亮，王志賢，余素榮等，2000．地電前兆場特征及有關(guān)地震預(yù)報(bào)問題探討．地震，20（1）：147—153．

汪志亮，鄭大林，余素榮，2002．地震地電阻率前兆異?，F(xiàn)象．北京：地震出版社，3—47．

王志賢，彭遠(yuǎn)黔，王麗，1994．地電阻率短臨前兆異常的判定方法研究．地震，（1）：88—96．

王志賢，彭遠(yuǎn)黔，1996a. 地震地電阻率法的統(tǒng)計(jì)預(yù)報(bào)研究．山西地震，（2）：34—39．

王志賢，彭遠(yuǎn)黔，張學(xué)民，1996b．地電阻率異常變化與地震成組活動關(guān)系的研究．華北地震科學(xué)，14（1）：19—26．

嚴(yán)玲琴，鄭衛(wèi)平，張輝等，2013．臨夏臺地電阻率變化與震兆現(xiàn)象分析．中國地震，29（1）：168—176．

尹相好，田鳳東，趙慶紅等，2010．山東莒縣深井電阻率資料分析．防災(zāi)減災(zāi)學(xué)報(bào)，26（2）：70—74．

張學(xué)民，王志賢，臧明珍等，1996．降雨對地電阻率干擾的分析．華北地震科學(xué)，14（4）：71—75．

張學(xué)民，劉素英，2003．華北地區(qū)成組地震前后的地電阻率異常特征．華北地震科學(xué)，21（4）：10—18．

張學(xué)民，李美，關(guān)華平，2009．汶川8.0級地震前的地電阻率異常分析．地震，29（1）：108—115．

趙玉林，錢復(fù)業(yè)，楊體成，1983．原地電阻率變化的實(shí)驗(yàn)．地震學(xué)報(bào)，5（2）：217—225.

趙玉林，錢復(fù)業(yè)，1978．唐山7.8級強(qiáng)震前震中周圍形變電阻率的下降異常．地球物理學(xué)報(bào)，（3）：181—190.

張國民，耿魯民，張永仙等，1995．構(gòu)造塊體的成組孕震模型和前兆場某些特征的分析．地震學(xué)報(bào)，17（1）：1—10．

中國地震局科技監(jiān)測司，2001．地震及前兆數(shù)字觀測技術(shù)規(guī)范（電場觀測）．北京：地震出版社，21，26，34—37.

Analysis on the Anomaly of Resistivity at Yangyuan Station, Hebei Province

Zhang Guoling1), Qiao Ziyun1), Jia Lifeng1), Liang Hongjie1), Han Heping2)and Wang Yanrong3)

1）Earthquake Administration of Hebei Province, Shijiazhuang 050021, China?2）Yangyuan Seismic Station, Zhangjiakou 075800, China?3）Sanhe Seismic Station, Sanhe 066600, China

The Geoelectric resistivity monitoring from Yangyuan station detected that the annual amplitude of resistivity in NE direction became smaller since 2014. After verification, we found that the observation system is working properly and there was correlation between rainfall and ground resistivity. The abnormal is still exist after removing the effects of rain by convolution filtering method. The NE resistivity and natural potential changes in 2014 is synchronous, so the reliability of earth resistivity anomaly is improving. The ground resistivity variation has a good relationship with the history of earthquake group at Shanxi, Hebei and Inner Mongolia border area, such as 1989 Datong earthquake, 1998 Zhangbei earthquake, et al.. The anomaly of resistivity at Yangyuan station is sensitive to stress change of the earth crust. Therefore, the anomaly of the earth resistivity has a higher reliability to earthquake prediction, and we should pay close attention to the seismic activity at Shanxi, Hebei and Inner Mongolia border area.

Resistivity; Yangyuan Station; Anomaly analysis

國家自然科學(xué)基金（41274079）和震情跟蹤定向工作任務(wù)（2015010403）共同資助

2015-02-01

張國苓，女，生于1986年。河北省地震局工程師。主要研究方向：地震電磁學(xué)和地震預(yù)測研究。 E-mail：zhangguoling002@163.com．