降雨對(duì)大柏舍地電臺(tái)地電阻率觀測(cè)的影響

吳鶴帥 羅 娜 趙志遠(yuǎn) 張 波 尹小兵

（中國(guó)河北 054000 紅山基準(zhǔn)地震臺(tái)）

0 引言

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

地電阻率觀測(cè)易受到觀測(cè)系統(tǒng)、電流、局部異常體、水位、降雨等因素的影響（李菊珍等，2004；陳遠(yuǎn)東等，2007；張國(guó)苓等，2013；羅娜等，2016）。降雨對(duì)地電阻率的干擾比較復(fù)雜，既存在即時(shí)效應(yīng)，又存在滯后效應(yīng)，與地電阻率季節(jié)性變化相關(guān)較好，需綜合考慮滲透過程、蒸發(fā)作用、電性結(jié)構(gòu)、降雨方式、地形等因素的影響（金安忠，1981；李飛等，2004；張偉峰等，2015）。如：昌黎地電臺(tái)地電阻率在降雨當(dāng)天出現(xiàn)大幅下降變化；嘉峪關(guān)地電臺(tái)地電阻率在降雨時(shí)數(shù)據(jù)反而上升（趙和云等，1986；喬子云等，2019）；代縣臺(tái)地電阻率年變動(dòng)態(tài)主要由降雨造成，且表層介質(zhì)影響系數(shù)為負(fù)，具有“反向年變”特征；周至臺(tái)地電阻率年變主要是因?yàn)榈乇肀邮芙涤?、溫度等季?jié)因素影響（石富強(qiáng)等，2014；郭文峰等，2019）?？梢越B透模型，分析降雨及滲透過程對(duì)地電阻率的影響機(jī)制，利用褶積濾波法去除降雨對(duì)地電阻率的影響（張學(xué)民等，1996；宋曉磊等，2006）。

河北省大柏舍地電臺(tái)（下文簡(jiǎn)稱大柏舍臺(tái)）始建于1966 年邢臺(tái)地震之后，當(dāng)時(shí)由蘭州地震大隊(duì)建立和管理，屬于國(guó)家基本臺(tái)，是我國(guó)建設(shè)較早的地電臺(tái)站。文中以大柏舍地電臺(tái)地電阻率觀測(cè)數(shù)據(jù)為基礎(chǔ)，結(jié)合觀測(cè)日志、降雨量數(shù)據(jù)及電性結(jié)構(gòu)，分析不同降雨方式對(duì)地電阻率觀測(cè)的影響，為地電阻率異常識(shí)別提供依據(jù)。

1 大柏舍臺(tái)地電阻率觀測(cè)概況

大柏舍臺(tái)臺(tái)址位于邢臺(tái)市隆堯縣柏舍村，地處隆堯斷裂和新河斷裂交叉部位（圖1），所在區(qū)域?qū)傩吓_(tái)老震區(qū)，周圍斷裂較為發(fā)育，為地震監(jiān)測(cè)敏感點(diǎn)，觀測(cè)數(shù)據(jù)較為可靠，映震效能較好（王志賢，1987；王志賢等，1997）。多年來監(jiān)測(cè)預(yù)報(bào)評(píng)估結(jié)果為A，震例R值平均評(píng)分約為0.24。

圖1 大柏舍臺(tái)地質(zhì)構(gòu)造Fig.1 Geological structure graph at Dabaishe station

該臺(tái)地電阻率測(cè)區(qū)地勢(shì)平坦，采用對(duì)稱四極法進(jìn)行觀測(cè)?！熬盼濉睍r(shí)期，使用ZD8B 儀進(jìn)行觀測(cè)，整點(diǎn)值采樣，布設(shè)NS向測(cè)道和EW向測(cè)道，兩測(cè)道相互垂直，裝置系數(shù)均為3.142，供電、測(cè)量電極布設(shè)在農(nóng)田中，采用60 cm × 50 cm × 2 cm 方形鉛板電極，均埋深2 m，供電極距A1B1=A2B2=1.5 km，測(cè)量極距M1N1=M2N2=0.5 km（圖2），外線路采用橫擔(dān)絕緣子架空方式。2016 年ZD8M 儀器入庫觀測(cè)，外線路采用地埋方式。2017 年11 月，大柏舍臺(tái)利用DUK-2A 電法測(cè)量系統(tǒng)重新進(jìn)行電測(cè)深勘探，采用水平層狀介質(zhì)模型，得到電測(cè)深電性結(jié)構(gòu)與影響系數(shù)，見表1。

表1 大柏舍電測(cè)深曲線反演的電性結(jié)構(gòu)Table 1 The electrical structure obtained by inverting electric sounding data at Dabaishe station

圖2 大柏舍臺(tái)地電阻率測(cè)區(qū)環(huán)境Fig.2 The geoelectric resistivity observation environment around Dabaishe station

2 影響系數(shù)理論

將視電阻率測(cè)區(qū)劃分為任意N塊區(qū)域，每塊區(qū)域視為介質(zhì)均勻，設(shè)電阻率為ρi，i=1，2，…，N，在測(cè)區(qū)電性結(jié)構(gòu)確定、觀測(cè)裝置和極距以及布極位置確定時(shí)，視電阻率ρa(bǔ)為各分區(qū)介質(zhì)電阻率的函數(shù)（錢家棟等，1985；Park et al，1991），公式如下

一般，各區(qū)域介質(zhì)電阻率在一定時(shí)間內(nèi)相對(duì)變化量較小，可對(duì)式（1）進(jìn)行Taylor 級(jí)數(shù)展開，其中二階及高階項(xiàng)遠(yuǎn)小于一階項(xiàng)，可略去不計(jì)，則視電阻率相對(duì)變化可表示為各區(qū)域介質(zhì)電阻率相對(duì)變化加權(quán)和，公式如下

同時(shí)，影響系數(shù)Bi滿足如下關(guān)系

3 降雨對(duì)地電阻率觀測(cè)的影響

選取大柏舍臺(tái)2011—2019 年地電阻率和降雨量資料，分析單日降雨、連續(xù)降雨、雷雨、強(qiáng)降雨對(duì)地電阻率觀測(cè)的影響。

3.1 單日降雨

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

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

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

3.2 連續(xù)降雨

2015 年5 月1 日至7 日，大柏舍臺(tái)當(dāng)?shù)爻霈F(xiàn)連續(xù)降雨，降雨總量達(dá)45 mm，其中5 月1 日降雨量25.1 mm。此次降雨對(duì)地電觀測(cè)的影響見圖4。由圖4 可見，大柏舍臺(tái)地電阻率未受降雨影響，NS 測(cè)道自然電位差5 月1 日起出現(xiàn)快速上升—下降變化，上升幅度9.5 mV，EW 測(cè)道自然電位差出現(xiàn)快速下降—上升變化，下降幅度4.3 mV；在自然電位差數(shù)據(jù)恢復(fù)過程中，5 日再次出現(xiàn)降雨，降雨量12.9 mm，NS、EW 測(cè)道數(shù)據(jù)再次出現(xiàn)快速上升、下降變化，降雨結(jié)束后，觀測(cè)數(shù)據(jù)緩慢恢復(fù)正常。

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

3.3 雷雨

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

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

3.4 強(qiáng)降雨

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

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

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

4 討論

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

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

雷電是一種強(qiáng)大的電磁干擾，近距離雷電電壓較高，在放電瞬間引起觀測(cè)數(shù)據(jù)大幅度改變，嚴(yán)重時(shí)對(duì)地充電，可導(dǎo)致測(cè)區(qū)自然電位急劇增長(zhǎng)，且數(shù)日不穩(wěn)定（錢家棟等，1985；王燚坤等，2011）。降雨伴隨雷電現(xiàn)象時(shí)，大柏舍臺(tái)自然電位差變化幅度較大。

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

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

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

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

5 結(jié)論

文中總結(jié)了不同降雨方式及降雨量對(duì)大柏舍地電阻率觀測(cè)的影響，并通過影響系數(shù)理論對(duì)結(jié)果進(jìn)行了探討，得到以下結(jié)論。

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

（2）觀測(cè)裝置老化對(duì)地電阻率觀測(cè)的影響。由于ZD8B 儀外線路采用架空方式，線路老化現(xiàn)象嚴(yán)重，降雨時(shí)線路漏電影響程度大于降雨對(duì)地電阻率觀測(cè)的影響，導(dǎo)致ZD8B 東西向地電阻率呈現(xiàn)階降變化；而背景場(chǎng)觀測(cè)的ZD8M 儀采用地埋線路，階變幅度比ZD8B儀小。

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

降雨對(duì)地電阻率的影響機(jī)制較復(fù)雜，文中僅定性分析降雨對(duì)大柏舍臺(tái)地電阻率觀測(cè)的影響，至于引起該臺(tái)地電阻率發(fā)生變化的降雨量大小及相應(yīng)的變化幅度，尚需大量數(shù)值模擬工作來進(jìn)行量化分析。