平武地電觀測實驗

劉學華

(四川省平武縣防震減災局，四川平武 622550)

平武地處四川北部，縣境西北有1976年曾經發生兩次松平7.2級地震的近南北向長約60 km的虎牙斷裂，北東向龍門山大斷裂斜貫平武縣境東南。“5·12”汶川特大地震時，平武地區因此斷裂的劇烈活動遭受了巨大破壞，造成數千人傷亡，地震烈度全境均為Ⅷ度，斷裂經過處局部烈度達到Ⅺ度。平武自有史記載以來曾經多次遭受強烈地震破壞，是一個存在較大地震危險性的潛在震源區。在這樣特殊的地震地質背景條件下，進行相關的地震前兆觀測和地震預測預報的探索研究是一件非常有意義的工作。觀測地電變化是一種很成熟且較早運用于地震預測預報的地震前兆方法之一，在國內和國際上都有相應的理論和實驗研究成果。地電觀測方法原理簡單，設備制作和安裝也較為簡便［1］。自20世紀70年代以來，有關單位在平武境內曾經開展了多種形式的地震前兆觀測工作，取得了大量的觀測數據。在場地地質條件良好和觀測方法運用得當的情況下，地電觀測方法可獲得有成效的觀測結果，在1976年松平大地震的成功預報中發揮過重要作用。但是近年來，地電觀測手段在一些地方逐步被放棄被淘汰，到2008年汶川地震前，平武地震臺只有測震和地下水氡觀測兩種監測手段。因儀器設備老化和“5·12”地震的破壞，水氡觀測不得不停止。僅剩下數字地震觀測，而無其他任何前兆觀測裝備，這種現狀對有著近40年地震監測歷史的臺站來說無疑是相當遺憾的。2009年7月上旬筆者接受甘肅隴南地震局趙衛國的建議和指導，在平武縣城西制作安裝了一套地下水電流觀測設備，利用這里有利的自然條件開始地電流觀測實驗工作。截至2014年底已經連續實驗觀測了5年多，其間于2011年5月對其進行了半自動化觀測改造，共取得實驗觀測數據30 000多個。通過分析發現，這些觀測數據具有良好的地震前兆對應，尤其是對較小范圍的有感地震對應較為明顯，是一種值得進一步繼續加強觀測和推廣運用的非常有利的前兆補充手段。

1 實驗觀測簡介

1.1 水點選取

在平武縣城涪江西岸馬園的一處被稱作龍洞的山石裂隙中，有一處常年徑流不息的泉水，其水文地質環境穩定，環境干擾極小，加之筆者正居家在該泉眼近處，所以觀測和管理維護也極為方便，因此該泉點是進行地電流觀測實驗的理想場所。該泉點在地質構造上位于龍門山后山斷裂帶北側與虎牙斷裂東側的交匯處。距離最近的龍門山后山斷裂之一的青溪斷裂約12 km，距南壩斷裂約40 km，距虎牙斷裂約25 km。這幾條主要斷裂帶上都曾經有過強烈地震活動，尤其是1976年虎牙斷裂帶松平7.2級地震群，2008年龍門山斷裂帶汶川8級地震，均在此區域產生了巨大的破壞，這里是一個研究地震與地電變化對應關系的好地方。區域地震地質及震中分布圖參見圖1。泉點附近出露地層的地質年代為志留系，泉點出水處巖性為志留系茂縣群灰色、灰綠色絹云英千枚巖夾變質石英砂巖和透鏡狀砂質結晶灰巖。泉水類型為基巖裂隙水，出水方式為裂隙自流，泉水清澈透明、白色、無臭、無味，流量大于1 000 ml/s，水溫變化在13.5～17.5℃之間。水質類型:HCO3-，Mg++，Ca++。水點電場性質應屬于上升泉電場，泉水涌出處具有電位極大值［2］。

圖1 區域地震地質及主要地震震中分布

1.2 地電觀測

平武地震臺已多年不進行地電等前兆觀測，所以沒有現成的儀器設備可供使用，選用部分部件只能自制加工。利用鉛塊和銅板制作兩個電極，微安表上加裝部分電阻元件、開關和電池構成測量電路。初期實驗觀測采用無源加有源測量相結合電路，實驗裝置測量電路如圖2所示。測量電路由室外和室內兩部分組成。室內測量部分由微安表和限流電阻及調整電位器構成無源電路，可以測量自然電場兩組數據;外加一組電阻和直流電池構成有源電路，測量自然和人工相結合電場兩項目。2011年5月上旬為了實現半自動化觀測，減輕人工負擔，將電路進行了簡化，去除有源部件，安裝了由視頻攝像頭和一臺老舊筆記本電腦構成的測量記錄系統，同時加大觀測密度，固定測量每小時實時記錄。安裝時，鉛電極(Pn)埋設在洞穴泉水出口深處，用以探測來自巖層深部由地下水攜帶的微電場信息，銅電極(Cu)安裝在距離洞口約10 m遠處一棵生長良好的大樹的南面，距離地面約4 m高處，電極方位為N30°E。兩個電極通過電纜連接至室內與測量電路相連接構成回路。使用電信電話鋼絲自承重銅質雙芯線電纜，線路電阻經測試小于10 Ω，信號傳輸距離約200 m。現場安裝時直接測量兩電極間電流，一般約為120 μA，由電纜引入室內后測量電流最大值為100 μA。

圖2 實驗原理圖

實驗觀測自2009年7月11日開始，采用每天人工多次讀取數據的方法進行。觀測初期數據變化平緩，原因是限流電阻取值不合理，經過對電路進行改進，7月19日開始記錄數據。讀數時，轉換測量檔位開關，分別讀取1、2、3、4檔位數據，填寫原始數據。原始觀測數據及記錄格式如表1所示。1、2檔數據為不加電直接測量地下水中的電流值，加入測量電路后1檔測量值在80μA左右。3、4檔數據為加上人工電場的電流值，外加人工電場后在4擋測量值約為120μA。設置第3、4檔的目的是檢測測量電路工作狀態是否正常。每次測量完畢恢復檔位到第1檔，使測量電路保持為自然電場待測狀態。2011年5月12日對觀測系統進行半自動化改造后，每小時自動通過攝像頭拍攝微安表頭盤面，然后及時由人工讀出每日24個整點數據。

表1 原始觀測數據及記錄格式

2 觀測數據及處理

根據實驗觀測方法，每天測定1至3個時點數據，截至2011年5月初，觀測到實驗電流值4組共計600多個數據。2011年5月12日后改為整點測量第1檔，數據量迅速增加，至2014年底，共取得觀測實驗數據32 000多個。觀測數據使用Excel電子表格錄入計算機進行存儲和處理，也可方便地轉存儲為其它格式。在Excel中設置為動態曲線，以便隨時檢測和糾正數據錄入過程中的人為錯誤，并方便觀察原始數據的動態變化。觀測實驗數據的原始曲線見圖3所示。由圖3可看出，觀測數據曲線可分為兩個階段。第一個階段自開始觀測到2010年底可看作是觀測系統與地質體結合逐步趨于穩定的時間段。第二階段是2011年初以后的時段是比較穩定和正常的，數據曲線可以看出明顯的年變形態(2010年已能看出年變)，盡管觀測過程中發現數據變化與大氣溫度無直接的相關性，但總體趨勢仍然與年溫差正相關，只是數據變化與溫度上升并不同步。觀測數據反映出和降雨降水關系密切，非常容易識別。降雨時電流值增加迅速，降雨結束后數據迅速回落到正常水平，與前兆異常持續時間較長有很大區別。使用Matlab等相關數據處理程序，也很容易濾除降雨干擾［3］。

圖3 原始觀測曲線

為了分析該點觀測數據與地震的對應關系，對觀測數據做了較為簡單的差分處理，如圖4所示。同時將觀測時段內影響較大的四川及周邊地區尤其是汶川地震余震區內ML≥4.5級地震與觀測數據曲線的變化情況進行對比。對于ML4級以上的地震，在觀測的原始數據和經過處理的數據曲線上大多能在地震前較早時間找到前兆異常的正突變反應，充分證明了該泉點的地電流觀測具有較好的映震效應，且震前異常提前量從1天到幾十天，異常反映良好。2011年5月12日改進為整點觀測后，觀測數據的異常變化對應地震的效果更明顯，圖5所示為部分時間段數據變化情況。根據已經取得的觀測數據，可以直觀地看出無論是無源的自然電場數據，還是有源觀測數據，全部4組數據完全存在線性相關的變化形態，與所使用的由線性電阻元件組成的線性電子測量電路系統的特性相吻合。

圖4 觀測數據差分曲線

3 觀測數據與地震對應的討論

圖5 改進觀測方法后的地震對應效果

通過5年多時間的觀測和與地震的對比，可以看到該水點的電流測值有明顯的地震對應效果，尤其是對ML≥4.5級以上的地震具有非常明顯震前對應，震級越大越明顯。但對于ML＜4.5級的地震的對應情況較差。異常形態均為正突跳。從2010年1月的最低值開始起伏上升，3月22日上升到最大值，至3月23日在1天之內迅速下降到最低，相對變化幅度達到40%，在起伏變化23天后青海玉樹發生7.1級地震［4］。從開始觀測半年后的2010年1月起，將幾次地震前異常非常明顯的數據進行了統計，結果見表2。由表2可見，地電流異常相對變化幅度越大，持續時間越長，未來發生的地震震級越大，震中距離也越遠，但單從此并不能得出距離的實際判斷結論。當繼續觀測積累的震例樣本增加到一定數量后，應可以進行異常時間、幅度、震中距和震級等之間的量化關系擬合計算。

表2 地震前異常數據統計

該水點經過近幾年的觀測證明，地電的數據變化非常穩定，無工業或其它人為等不良環境干擾，依據取得的觀測數據可以較為直觀地對地震異常與否做出初步判定［5］。在觀測過程中，當出現較為明顯的測值突跳，而且其變化幅度與背景觀測值相比至少達到40%的變幅時，可初步判定為地震前兆異常，但需要排除雷雨雷電降雨天氣或其他干擾因素(除降雨以外無其他干擾出現的記錄)，若需要可進一步加密觀測，以確定異常變化趨勢。再觀察數據曲線形態形成突跳狀尖峰并回落到背景值附近，即可完全確定本次地震前兆異常確實存在，可預示將有地震在附近發生。如2010年5月24日1檔測值變化幅度達到54%，25日測值基本回落到背景值附近，即為一突跳，后于25日都江堰發生ML5.2級、29日玉樹發生Ms6.0級地震。在這種測值回落非常明顯的形態下，即使未來發生的地震震級沒有超過5級，但在周邊地區尤其是龍門山斷裂帶及其附近大都會有4級左右地震發生。若觀測數據在高值附近起伏變化，逐漸形成緩慢變化的凸峰并且持續時間明顯較長時，則為地電場大背景調整，異常為假，可撤銷預報警戒。如2010年4月27日數據變化幅度達到40%，測值一直在高值點附近起伏，到5月5日緩慢恢復到正常值附近，形成一個凸起的峰值，此即為一個正常調整狀態，而非地震前兆異常。

由于觀測采用人工讀數的方式，因此觀測的精度不高。同時由于制作的電極不夠標準和規范，例如電極尺寸過小，可能造成不能完全拾取地電流信號，還有靈敏度太低等缺陷。考慮到業余人工觀測作業工作量負擔，無法保證不間斷觀測。而非連續的觀測，數據離散性可能過大，每天觀測一個(或幾個)數據，進行濾波(采樣率為1～1/24)以后會遺漏地震前兆信息。

已經取得的觀測數據和震例對應結果，還不能有效進行數據與地震信息的相關性分析，有待繼續觀測并進行優化總結，以利應用于實際的地震監測預報工作。使用基于GIS的地震分析預報系統等專業軟件進行研究和處理，可能會得到更準確、更科學的結論。采用更先進的電極和使用全自動采樣的方式提高數據觀測精度，例如采樣頻率提高，減少人為誤差以利于充分拾取地震前兆信息，觀察其演變過程。該泉點水量充沛，冬暖夏涼，水質良好，還可輔助水流量、水溫、水離子組分等多種前兆觀測，以利于多方法、多途徑地充分獲取地震前兆信息。理想的途徑是在一定的項目支撐下采用國家或行業標準設備進行自動化數據采集［6］。

該試驗是筆者工作之余進行的。該實驗能夠得以開展和堅持，特別感謝甘肅隴南地震局趙衛國副局長的鼓勵和大力支持，2011年5月1日在隴南交流學習時獲贈一臺SE-11型地電觀測儀，也安裝在本實驗點近處同步進行觀測(有關觀測情況擬另文討論)。在本文寫作成稿過程中得到中國地震局地殼研究所楊選輝研究員的傾力指導，修改時還得到了《地震地磁觀測與研究》編輯部和《四川地震》編輯部的老師和專家的指導，在此一并致以誠摯的謝意。

［1］何世根，沈啟興，等.地電觀測技術［M］.北京:地震出版社，2000.50-67.

［2］國家地震局《地震觀測技術》編寫組.地震觀測技術［M］.北京:地震出版社，1994:217-222.

［3］馬君釗，等.大地電場觀測各類干擾源的調研與分析［J］.地震地磁觀測與研究，2010(5):65-72.

［4］張小濤，等.地電場方法在川滇地區地震觀測中的應用［J］.地震地磁觀測與研究，2010(6):93-98.

［5］萬迪堃，汪成民，等.地下水異常與地震短臨預報［M］.北京:地震出版社，1993:147-172.

［6］孔令昌，王志敏，王桂清，等.大地微電流檢測儀的研制和觀測［J］.地震地磁觀測與研究，2002(5):69-71.