降雨對大柏舍地電臺地電阻率觀測的影響

吳鶴帥 羅 娜 趙志遠 張 波 尹小兵

（中國河北 054000 紅山基準地震臺）

0 引言

1966 年邢臺地震發生后，我國將地球物理勘探電阻率方法用于地震監測預報（杜學彬，2010；解滔等，2016；張國苓等，2019）。在1976年唐山MS7.8、1988年瀾滄—耿馬MS7.6、2008 年汶川MS8.0 等中強地震前多次記錄到地電阻率異常，且異常通常表現為趨勢性下降、破年變異常變化（錢復業等，1980，1991；錢家棟，1993；汪志亮等，1995；王志賢等，1996；張學民等，2009；錢家棟等，2013；杜學彬等，2015）。地電阻率觀測已成為研究地震預報的有效手段之一。

地電阻率觀測易受到觀測系統、電流、局部異常體、水位、降雨等因素的影響（李菊珍等，2004；陳遠東等，2007；張國苓等，2013；羅娜等，2016）。降雨對地電阻率的干擾比較復雜，既存在即時效應，又存在滯后效應，與地電阻率季節性變化相關較好，需綜合考慮滲透過程、蒸發作用、電性結構、降雨方式、地形等因素的影響（金安忠，1981；李飛等，2004；張偉峰等，2015）。如：昌黎地電臺地電阻率在降雨當天出現大幅下降變化；嘉峪關地電臺地電阻率在降雨時數據反而上升（趙和云等，1986；喬子云等，2019）；代縣臺地電阻率年變動態主要由降雨造成，且表層介質影響系數為負，具有“反向年變”特征；周至臺地電阻率年變主要是因為地表薄層受降雨、溫度等季節因素影響（石富強等，2014；郭文峰等，2019）。可以建立滲透模型，分析降雨及滲透過程對地電阻率的影響機制，利用褶積濾波法去除降雨對地電阻率的影響（張學民等，1996；宋曉磊等，2006）。

河北省大柏舍地電臺（下文簡稱大柏舍臺）始建于1966 年邢臺地震之后，當時由蘭州地震大隊建立和管理，屬于國家基本臺，是我國建設較早的地電臺站。文中以大柏舍地電臺地電阻率觀測數據為基礎，結合觀測日志、降雨量數據及電性結構，分析不同降雨方式對地電阻率觀測的影響，為地電阻率異常識別提供依據。

1 大柏舍臺地電阻率觀測概況

大柏舍臺臺址位于邢臺市隆堯縣柏舍村，地處隆堯斷裂和新河斷裂交叉部位（圖1），所在區域屬邢臺老震區，周圍斷裂較為發育，為地震監測敏感點，觀測數據較為可靠，映震效能較好（王志賢，1987；王志賢等，1997）。多年來監測預報評估結果為A，震例R值平均評分約為0.24。

圖1 大柏舍臺地質構造Fig.1 Geological structure graph at Dabaishe station

該臺地電阻率測區地勢平坦，采用對稱四極法進行觀測。“九五”時期，使用ZD8B 儀進行觀測，整點值采樣，布設NS向測道和EW向測道，兩測道相互垂直，裝置系數均為3.142，供電、測量電極布設在農田中，采用60 cm × 50 cm × 2 cm 方形鉛板電極，均埋深2 m，供電極距A1B1=A2B2=1.5 km，測量極距M1N1=M2N2=0.5 km（圖2），外線路采用橫擔絕緣子架空方式。2016 年ZD8M 儀器入庫觀測，外線路采用地埋方式。2017 年11 月，大柏舍臺利用DUK-2A 電法測量系統重新進行電測深勘探，采用水平層狀介質模型，得到電測深電性結構與影響系數，見表1。

表1 大柏舍電測深曲線反演的電性結構Table 1 The electrical structure obtained by inverting electric sounding data at Dabaishe station

圖2 大柏舍臺地電阻率測區環境Fig.2 The geoelectric resistivity observation environment around Dabaishe station

2 影響系數理論

將視電阻率測區劃分為任意N塊區域，每塊區域視為介質均勻，設電阻率為ρi，i=1，2，…，N，在測區電性結構確定、觀測裝置和極距以及布極位置確定時，視電阻率ρa為各分區介質電阻率的函數（錢家棟等，1985；Park et al，1991），公式如下

一般，各區域介質電阻率在一定時間內相對變化量較小，可對式（1）進行Taylor 級數展開，其中二階及高階項遠小于一階項，可略去不計，則視電阻率相對變化可表示為各區域介質電阻率相對變化加權和，公式如下

同時，影響系數Bi滿足如下關系

3 降雨對地電阻率觀測的影響

選取大柏舍臺2011—2019 年地電阻率和降雨量資料，分析單日降雨、連續降雨、雷雨、強降雨對地電阻率觀測的影響。

3.1 單日降雨

2014 年4 月25 日，大柏舍地電測區出現降雨天氣，臺站記錄降雨量為27.3 mm。觀測數據顯示，該臺地電阻率無明顯變化，自然電位差受到影響（圖3），其中：NS 測道自然電位差出現上升—下降變化，上升幅度11.7 mV；EW 測道自然電位差出現下降—上升變化，下降幅度17.2 mV。降雨停止后，隨著天氣放晴地面慢慢變干，兩測道自然電位差在幾天內逐步恢復。

七是打好精準扶貧攻堅戰，走中國特色減貧之路。黨的十八大以來，黨中央以前所未有的政策力度向貧困宣戰，脫貧攻堅取得了舉世矚目的偉大成就。解決好脫貧問題的“最后一公里”，需要充分發揮政治和制度優勢，精準施策，把提高脫貧質量放在首位，將發展農業產業、易地扶貧搬遷、生態補償、發展教育、社會保障兜底的舉措做細做實，為建成全面小康社會打贏這場具有決定意義的攻堅戰。

圖3 單日降雨對地電阻率和自然電位差影響Fig.3 The effect of single-day rainfall on earth resistivity and spontaneous potential difference observations

3.2 連續降雨

2015 年5 月1 日至7 日，大柏舍臺當地出現連續降雨，降雨總量達45 mm，其中5 月1 日降雨量25.1 mm。此次降雨對地電觀測的影響見圖4。由圖4 可見，大柏舍臺地電阻率未受降雨影響，NS 測道自然電位差5 月1 日起出現快速上升—下降變化，上升幅度9.5 mV，EW 測道自然電位差出現快速下降—上升變化，下降幅度4.3 mV；在自然電位差數據恢復過程中，5 日再次出現降雨，降雨量12.9 mm，NS、EW 測道數據再次出現快速上升、下降變化，降雨結束后，觀測數據緩慢恢復正常。

圖4 連續降雨對地電阻率和自然電位差影響Fig.4 The effect of continuous rainfall on earth resistivity and spontaneous potential difference observations

3.3 雷雨

2017 年10 月8 日至10 日，大柏舍地電測區內出現降雨天氣，其中8 日為雷雨天氣，降雨總量為53.3 mm。此次降雨對該臺地電觀測的影響見圖5。由圖5 可見，NS 測道自然電位差8 日出現快速大幅上升—小幅下降變化，18 時至19 時上升幅度41.8 mV，19 時至20 時突降16.4 mV，21 時后開始緩慢下降；EW 測道自然電位差18 時至19 時快速小幅上升6.5 mV，19 時至20 時大幅下降15.3 mV，隨后緩慢上升；9 日、10 日再次降雨，兩測道自然電位差雖有小幅波動，但整體呈恢復狀態。

圖5 雷雨對地電阻率和自然電位差影響Fig.5 The effect of thunderstorms on earth resistivity and spontaneous potential difference observations

3.4 強降雨

2016 年7 月18 日—21 日，河北省在低渦及副高外圍暖濕氣流共同影響下，出現當年最強一次強暴雨過程，全省平均降水量達154.6 mm，其中邢臺市平均降水204 mm，柏鄉、寧晉、隆堯等4 縣相繼滯洪進水0.2—1.5 m。大柏舍臺19 日、20 日單日降雨量為64.5 mm、50.5 mm，當月累計降雨總量高達173.4 mm，較2015 年同期偏多5 倍（邢臺市水務局，2016；吳麗英等，2018）。

此次降雨過程突發性強、范圍廣、強度大，大柏舍臺不僅自然電位差受到暴雨干擾，地電阻率觀測也受到影響，見圖6。由圖6 可見，7 月19 日至20 日，大柏舍臺自然電位差出現數次大幅階變，EW 測道最大變化幅度為43 mV，NS 測道最大變化幅度為20 mV，21 日后逐步恢復正常變化；7 月19 日以來，地電阻率NS 測道僅7 月20 日出現階變，最大變化幅度為0.37 Ω·m，而EW 測道出現大幅臺階變化，且持續至8 月14 日，最大變化幅度為1.27 Ω·m，8 月15 日起逐漸恢復正常形態。

圖6 強降雨對地電阻率和自然電位差影響Fig.6 The effect of heavy rainfall on earth resistivity and spontaneous potential difference observations

4 討論

降雨當天或者第二天，地表形成薄飽水層，表層介質含水率升高，表層真電阻率下降。若臺站表層介質影響系數為正，地電阻率觀測值上升；若臺站表層介質影響系數為負，地電阻率觀測值將下降。隨著雨水下滲、蒸發等，雨水影響深度逐漸增加，水飽和程度逐步減小，至恢復原有程度（趙和云等，1986）。

根據大柏舍電測深曲線反演的電性結構（表1），計算AB/2=750 m 時各層影響系數，結果見表2。由表2 可見，B3是大柏舍臺地下第三層的影響系數，所占比重較大，表明該臺地電阻率觀測主要反映第三層介質變化，而B1系數相對較小，小幅降雨對大柏舍地電阻率短期影響不明顯，僅自然電位差受到干擾。

雷電是一種強大的電磁干擾，近距離雷電電壓較高，在放電瞬間引起觀測數據大幅度改變，嚴重時對地充電，可導致測區自然電位急劇增長，且數日不穩定（錢家棟等，1985；王燚坤等，2011）。降雨伴隨雷電現象時，大柏舍臺自然電位差變化幅度較大。

表層影響系數B1較小，降雨量較少時，可以看出地電阻率觀測不受影響，但2016 年7 月19 日至20 日降雨量高達115 mm，大柏舍臺表層介質電阻率發生改變，從而引起地電阻率變化。同時，ZD8M 儀器EW 向地電阻率出現階升變化，幅度約0.1 Ω·m，證明此次降雨確實對地電阻率觀測造成了影響（圖7）。

圖7 2016 年7 月5 日—8 月23 日大柏舍地電阻率測值（ZD8M）Fig.7 The earth resistivity observations by ZD8M at Dabaishe station from July 5 to August 23,2016

持續幾天的強降雨造成空氣濕度增大，同時，因ZD8B 地電阻率觀測系統外線路使用年限較長，出現老化現象，線路絕緣性能明顯降低。7 月21 日臨時檢查結果顯示，EW 方向M 極測量線路絕緣電阻為80 MΩ，N 極測量線路絕緣電阻為200 MΩ（表3），線路出現漏電現象，線路漏電導致ZD8B 儀地電阻率觀測數據出現大幅階降變化。而用于背景場觀測的ZD8M 儀采用地埋線路，因表層介質影響系數B1為負，強降雨造成地電阻率觀測值上升，階變幅度比ZD8B 儀數據記錄幅度小。

表3 EW 向測量線路絕緣電阻Table 3 EW measuring component insulation resistance

5 結論

文中總結了不同降雨方式及降雨量對大柏舍地電阻率觀測的影響，并通過影響系數理論對結果進行了探討，得到以下結論。

（1）降雨對地電觀測的影響。在降雨較小或者雷雨天氣條件下，自然電位差變化較大。大柏舍臺表層介質系數較小，單日降雨且降雨量較小時，地電阻率觀測不受影響，自然電位差受到影響，并于降雨結束后緩慢恢復；出現連續降雨時，在自然電位差恢復過程中再次出現上升、下降變化；受雷雨影響時，自然電位差NS 測道出現快速大幅上升—小幅下降變化，EW 測道出現小幅下降—大幅上升變化，變化幅度比單純降雨時偏大。大范圍強暴雨時，大柏舍臺自然電位差和地電阻率觀測同時受到影響。因測區表層介質對地電阻率相對變化影響系數為負值，降雨后地表形成薄飽水層，表層真電阻率降低，地電阻率觀測值升高。

（2）觀測裝置老化對地電阻率觀測的影響。由于ZD8B 儀外線路采用架空方式，線路老化現象嚴重，降雨時線路漏電影響程度大于降雨對地電阻率觀測的影響，導致ZD8B 東西向地電阻率呈現階降變化；而背景場觀測的ZD8M 儀采用地埋線路，階變幅度比ZD8B儀小。

（3）地電臺站可采取一些措施降低降雨影響。一般，地電阻率觀測受降雨影響較明顯，自然電位差則無變化，例如：昌黎臺受降雨影響，地電阻率測值出現大幅度下降；嘉峪關臺大雨當日，地電阻率測值出現大幅度突升。而大柏舍臺在小幅降雨時，僅自然電位差受到影響，大范圍強降雨才會對地電阻率產生影響。為此，在日常觀測、臺站建設中，可增加電極距或采用深井地電阻率觀測，抑制降雨對地電觀測產生的干擾。

降雨對地電阻率的影響機制較復雜，文中僅定性分析降雨對大柏舍臺地電阻率觀測的影響，至于引起該臺地電阻率發生變化的降雨量大小及相應的變化幅度，尚需大量數值模擬工作來進行量化分析。