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

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

?

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

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

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

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

地電阻率 陽原 異常分析

引言

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

本文基于陽原臺地電阻率觀測呈現出的有規則的冬高、夏低年周期變化特征，以及2014年2月出現的NE測向年變幅減小（年變畸變現象）和NW測向的年變形態及幅度與往年相當的資料（圖1），并結合觀測系統檢測和背景場環境以及各年度日均值數據的相關性分析結果，重點分析了降雨量對陽原臺地電阻率NE測向的趨勢性影響。同時還進一步驗證了陽原臺地電阻率NE測向出現的異常變化（年變畸變現象）的可靠性。

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

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

陽原臺地電阻率觀測有N22°E和N68°W兩個測向，使用的儀器為ZD8B型地電儀，供電極距和測量極距分別為1.2km和0.2km。測區的電測深曲線為QQ型，地下電性結構可分為三層：第一層為第四紀黃土層，夾有礫石，電阻率為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。第三層以下為砂質灰巖層，電阻率為2.75—3.24Ω·m。圖2為陽原臺巖性柱狀圖、電測深曲線和區域地質構造圖。由于測區內有大秦鐵路穿過，同時還有小型液化氣站位于測區中央，所以對數據觀測有一定的影響。

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

2.1 觀測系統檢測

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

表1 部分電極的接地電阻結果

2.2 背景場環境

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

2.3 地電阻率與降雨的關系

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

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

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

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

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

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

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

2.5 震例總結

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

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

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

3 結論與討論

通過對觀測系統的檢測，確認陽原臺地電觀測儀器、外線路、電極工作狀態正常。測區背景場環境調查和輔助資料分析表明，降雨量和地下水位引起地電阻率數據異常的可能性較小，陽原臺的地電阻率變化反映了真實的應力變化，可視為地震前兆異常。同時，通過地電阻率與自然電位的對比分析后發現，陽原臺地電阻率的NE測向與自然電位在2014年的變化具有同步性，驗證了地電阻率異常的可靠性；陽原臺地電阻率的異常變化多以年變畸變為主，2014年以來的異常形態與1998年大同地震和2006年文安地震發生前的異常形態較為類似。綜合分析認為，陽原臺地電阻率異常變化的可信度較高，可視為地震前兆異常，類型為背景異常。

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

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

錢復業，趙玉林，許同春，1987．地電阻率季節干擾變化分析．地震學報，9（3）：289—302．

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

汪志亮，王志賢，余素榮等，2000．地電前兆場特征及有關地震預報問題探討．地震，20（1）：147—153．

汪志亮，鄭大林，余素榮，2002．地震地電阻率前兆異常現象．北京：地震出版社，3—47．

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

王志賢，彭遠黔，1996a. 地震地電阻率法的統計預報研究．山西地震，（2）：34—39．

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

嚴玲琴，鄭衛平，張輝等，2013．臨夏臺地電阻率變化與震兆現象分析．中國地震，29（1）：168—176．

尹相好，田鳳東，趙慶紅等，2010．山東莒縣深井電阻率資料分析．防災減災學報，26（2）：70—74．

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

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

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

趙玉林，錢復業，楊體成，1983．原地電阻率變化的實驗．地震學報，5（2）：217—225.

趙玉林，錢復業，1978．唐山7.8級強震前震中周圍形變電阻率的下降異常．地球物理學報，（3）：181—190.

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

中國地震局科技監測司，2001．地震及前兆數字觀測技術規范（電場觀測）．北京：地震出版社，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

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

2015-02-01

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