大柏舍深井地電阻率觀測布極方式分析

寧亞靈，許家姝 ，解 滔 ，張國苓，盧 軍

1.山西省地震局,太原 030021 2.太原大陸裂谷動力學國家野外科學觀測研究站，太原 030021 3.吉林大學地球探測科學與技術學院，長春 130026 4.中國地震臺網中心,北京 100045 5.河北省地震局,石家莊 050021

0 引言

我國自1966年邢臺M7.2(M為震級)地震后開始了連續的定點地電阻率觀測，臺站通常布設2～3個測道，布極方式普遍采用對稱四極裝置，地表觀測電極埋設深度通常為1.5～2.0 m，觀測極距AB為500～2 400 m，地下探測深度為數百米或更深不等[1-2]。在近50年的連續監測中，在多次大地震前記錄到了顯著的中短期地電阻率異常[3-9]，且對發生在臺網內一些地震的三要素實施了1年時間尺度的中短期預測[10-11]，說明地電阻率是地震監測預報有效的方法之一。

隨著地區經濟的發展，蔬菜大棚、水管、鋼纜等金屬管線以及測區挖土蓄水、倉庫和道路等基礎設施建設改變了測區地表的電性結構，使地電阻率觀測數據質量嚴重下降，給地震監測預報數據分析帶來巨大困難[12-14]。為了有效地抑制來自測區地表淺層的干擾，我國地電工作者提出了深井觀測的方式，自20世紀80年代開始相繼進行了一些深井實驗觀測，并逐步開展了深井地電阻率觀測的相關理論研究，目前已有多個臺站在進行連續觀測[15-23]。已有的研究結果顯示，深井觀測能有效抑制地電阻率觀測中因地下水位變動、地表溫度變化和季節性降雨引起的年變化現象和地表局部電性結構變化引起的干擾。深井觀測在抑制地表淺層干擾的同時，更需要捕捉由孕震過程應力作用引起的測區深部介質變化信息。在電極埋深一定時，地電阻率觀測的探測范圍很大程度上取決于觀測極距。觀測極距越小，觀測值受地表的干擾越小，但所反映的深部信息也越少；由于測區地下介質的非均勻性，在區域應力作用下會呈現非均勻性變化，極距較小時容易遺漏相關信息。因此，深井觀測仍然需要采用較大的觀測極距，以反映更大三維空間范圍的介質變化信息。但隨著觀測極距的加大，觀測值所反映的空間范圍擴大，地表干擾對觀測的影響也相應地增加，深井觀測的作用也逐漸減弱。地電阻率影響系數理論可以定量地給出測區各區域介質電阻率變化對地電阻率觀測值變化的影響程度[24]，因此，通過計算在不同觀測極距和不同電極埋深時各區域介質對地電阻率的影響系數，可以評估深井觀測對地表干擾的抑制能力和對深部介質電阻率變化的響應能力[23]。

大柏舍深井地電阻率水平測道電極埋深100 m，觀測極距僅為60 m，能夠很好地抑制地表灌溉和季節性降雨等因素產生的干擾[25]，但能否較好地捕捉深部介質信息還需要進一步的分析。本文依據地電阻率影響系數理論，對大柏舍深井地電阻率觀測資料抑制地表干擾的原因和捕捉深部介質變化信息的能力進行分析，為今后深井地電阻率建設中電極埋深和觀測極距的選擇提供一定的參考。

1 臺站簡介

大柏舍地震臺位于河北省隆堯縣以南約5 km處，處于寧晉凸起和巨鹿凹陷交界的隆堯斷裂和新河斷裂的交叉部位(圖1)，臺站周圍斷裂較為發育，且多數近期仍處于活動階段。臺站所在區域為地震多發的邢臺老震區，是一個地震前兆監測的敏感點。大柏舍臺地電阻率觀測系統于1967 年建成并投入運行，是我國最早開展地震分析預報工作的地電阻率臺站，觀測數據具有較好的映震效能[26]。地電阻率測區主要為農業用地，地勢平坦，地形無明顯高差，地下水埋深約50 m。地表地電阻率觀測儀器目前為ZD8B型數字地電儀，采用對稱四極測量裝置(圖2a)，布設EW、NS向兩條測道，供電極距均為AB=1 500 m，測量極距均為MN=500 m，電極埋深為2.5～3.0 m。地表觀測數據質量較高，夏季降雨增加、溫度升高，表層介質電阻率降低時觀測值處于年變峰值，冬季則處于年變低值，呈現“夏高冬低”型年變形態。大柏舍臺于2010年1月1日增加深井地電阻率觀測，觀測儀器最初為深孔電極地電儀ATS-SR，2013年6月更換為ZD8B地電儀。深井觀測布設一條EW向水平測道和一條垂直向測道，各測道電極沿直線分布(圖2b)，供電極距為A1B1=A2B2=60 m，測量極距為M1N1=M2N2=20 m。水平測道所有電極埋深均為100 m，垂直測道4個電極埋深分別為40、60、80和100 m，垂直測道供電電極B2和水平測道供電電極A1共用一個電極。

圖1 大柏舍地區地質構造圖Fig.1 Geological structure map of the area in Dabaishe

a. 地表觀測布極圖；b. 深井觀測布極圖。圖2 大柏舍臺地表與深井地電阻率觀測電極布設示意圖Fig.2 Sketch map of the lay about the earth’s surface and deep-well apparent resistivity in Dabasishe seismic station

2 觀測資料分析

大柏舍深井地電阻率水平測道觀測值年變化形態不明顯(圖3a)，垂直測道則具有夏低冬高型年變形態，年變幅度為0.6 Ω·m(圖3b)。大柏舍臺位于河北平原地區，該區域農業灌溉用水量較大，地下水開采比較嚴重，地下水位埋深較大(圖3c)。對隆堯地區地下水位的調查顯示，地下水位受農業開采及降雨的影響有較為明顯的冬高夏低的年變化特征[27]。每年6月大量開采地下水用于冬小麥生長階段的灌溉，地下水得不到及時的補給而出現水位下降，期間垂直測道地電阻率觀測值呈現下降變化；7月雨季開始，地下水開采量也相應減少，同時受太行山山前側向徑流的補給作用，地下水位開始回升，垂直測道地電阻率則呈現上升變化；次年1—2月，由于停止農業灌溉地下水開采，地下水于2月底達到年內最高水位，期間垂直測道地電阻率出現上升變化。垂直測道地電阻率的年變化與水位變動有很好的對應關系：地下水位埋深較小時，地電阻率觀測值較高；地下水位埋深較大時，地電阻率觀測值則較低；在年變形態清晰的年份中，水位出現快速大幅度變動時，垂直測道地電阻率也出現相應的快速變化。大柏舍附近的地下水位資料表明，臺站所在區域地下水位屬于埋深較深的類型，受季節性降雨和抽水影響，近幾年地下水位埋深在50 m左右，并且逐年下降。垂直測道地電阻率A2、M2電極分別位于地下40和60 m，地下水位埋深在這兩個電極之間波動，且距離電極較近，水位的波動對垂直測道地電阻率觀測造成較大影響。

ρa.地電阻率觀測值；Hw.地下水位深度。圖3 深井地電阻率觀測資料和地下水位資料Fig.3 Data of the deep-well apparent resistivity observation and groundwater table

3 探測靈敏度理論

如果將測區地下三維空間按任意大小劃分為N塊三維子區域，每一塊子區域內視為均勻介質，電阻率為ρi，i=1，2，…，N，在觀測裝置、觀測極距和布極位置確定時，地電阻率觀測值ρa是各子區域介質電阻率的函數[28-29]：

(1)

通常情況下，各子區域介質電阻率變化緩慢，在一定時間段內電阻率相對變化量非常小，即Δρi/ρi<<1。因此，將式(1)作Taylor級數展開，略去二階項和高階項，地電阻率觀測值的相對變化可以表示為各子區域介質電阻率相對變化的加權和形式：

(2)

式中，Bi被稱之為地電阻率影響系數：

(3)

同時，所有子區域影響系數Bi滿足如下關系[30]：

(4)

從影響系數的定義式(2)中可知，在影響系數絕對值較大的區域，介質電阻率發生變化時將對地電阻率觀測值產生較大的影響，反之則影響很小。

4 布極方式分析

大柏舍地震臺地電阻率測區覆蓋層較厚，基巖埋深約700 m，地下介質導電性良好，電測深曲線顯示地層可視為QH型電性結構(圖4a)。依據電測深曲線采用水平層狀介質模型反演的電性結構示于圖4a：第一層介質為黃土層，厚度約5 m，電阻率約75 Ω·m；第二層介質為礫石和細砂互層，厚度約100 m，電阻率為25 Ω·m；第三層介質厚度約580 m，為亞砂土和亞黏土互層，電阻率較低，為12 Ω·m；其下為電阻率較高的基巖。

依據圖4a所示的水平層狀電性結構，計算水平測道在供電極距AB=60 m、測量極距MN=20 m時各層介質影響系數隨深度的變化，示于圖4b。用ρaw表示深井地電阻率觀測值，在電極埋深100 m時ρaw相對變化值則可以表示為如下形式：

(5)

井下觀測的目的在于抑制地表干擾和突出由孕震引起的深部巖層介質電阻率變化。觀測裝置埋深固定時，小極距對電極所在層位電阻率變化反映較大，隨著極距增加，對深部信息的反映能力增加。觀測極距足夠大后，各層影響系數與地表觀測時趨于相近，井下觀測的作用越來越小，失去意義。由于大柏舍深井地電阻率觀測電極埋深較深，觀測極距較小，基本能避免表層干擾并突出目標層介質電阻率變化。深井地電阻率水平測道主要體現了第二層和第三層介質的電阻率信息；第一層的影響系數為0.000 1，地表淺層介質對地電阻率觀測的影響可以忽略，因此水平測道地電阻率沒有觀測到明顯的年變化；第四層的影響系數也非常小，因此水平測道地電阻率對深部基巖電阻率變化的響應能力也是微乎其微的。地表觀測影響系數隨觀測極距變化示于圖4c，在極距AB=1 500 m時表層介質影響系數為負(由于是對數坐標，負數不能取對數，所以用絕對值，負值部分用虛線表示)，表層介質電阻率的升高和降低將引起地電阻率觀測值的下降和上升變化。春季降雨量開始增加、地表介質含水率增加、電阻率降低后地表視電阻率觀測值上升；進入秋季后降雨量減小，地表介質電阻率上升，地表視電阻率觀測值下降，因此大柏舍臺地表觀測出現“夏高冬低”型反常年變形態。電極埋深固定為H=100 m時，影響系數隨觀測極距的變化示于圖4d，可以看出：觀測極距較小時，視電阻率主要反映第二層介質電阻率變化；隨著觀測極距的增加，第二層介質影響系數減小、第三層介質影響系數增加，第一、四層基本為負。各層影響系數與地表觀測時逐漸趨于一致，說明對于固定的電極埋深，觀測極距超過一定范圍之后將失去深井觀測的作用。

圖4 大柏舍臺電性結構與影響系數Fig.4 Electrical structure and influence coefficient in Dabasishe seismic station

依據圖4a所示電性結構，計算水平測道觀測時各層介質影響系數隨觀測極距和裝置埋深的分布如圖5所示。在極距AB/2=150 m時，埋深100 m以上第一層介質影響系數變化梯度較大，100 m以下影響系數很小但變化較為平緩(圖5a)。第二層介質在小極距淺層觀測時影響系數較大，別的情況都很小(圖5b)。電極埋深在100 m以上，第三層介質影響系數占主導地位(圖5c)。極距AB/2在300 m以上、電極埋深在150 m以上才能探測到第四層介質的變化(圖5d)。分析圖5可知，小極距觀測主要反映裝置所在區域介質信息，在裝置埋深較淺時，主要反映淺層介質信息。大柏舍臺第三層介質厚度較大且電阻率很低，地電阻率觀測值主要反映該層介質信息，只有在觀測極距較大且裝置埋深也較大時，深井觀測才能對底層介質信息有所反映。淺層介質電阻率易受地下水位變動、含水飽和度季節性增減和溫度變化的影響，深井觀測應盡可能少地反映這部分干擾信息。通常認為，深層介質較之淺層松散沉積層更能有效地傳遞應力，孕震或構造應力主要引起深層介質電阻率的變化，因而大柏舍臺深井觀測應以第三層和底層介質為主要觀測對象。從圖5中可以看出，水平測道供電極距AB/2取250～350 m、裝置埋深H>150 m能滿足需求。考慮到地下潛水位在50 m左右變化，垂直測道的頂端電極埋深應在地下水位以下一段距離為宜，以減少地下潛水位變動對地電阻率觀測的影響。

圖5 大柏舍臺各層介質影響系數隨觀測極距和埋深的變化Fig.5 Changes of the influence coefficient with the each layer about the observation distance and depth in Dabaishe seismic station

5 討論與結論

1)深井觀測通過地層屏蔽效應和增加觀測裝置與地表干擾源的距離來實現對地表淺層干擾的抑制，因此，只要觀測裝置埋設足夠深，總是可以有效地抑制由非構造因素引起的地表淺層介質電阻率變化對觀測產生的干擾。另一方面，深井觀測也減小了觀測裝置與下伏地層之間的距離，有利于獲取深部介質變化信息。

2)根據目前深井觀測工程投入和對觀測系統穩定性的要求，觀測裝置埋深能到達250 m左右；對于多數臺站而言，在這一深度采用適當的極距已經可以有效地對基巖電阻率變化進行觀測。

3)大柏舍臺現有的深井觀測水平測道雖能有效地抑制來自地表的干擾信息，但是由于觀測極距較小，對深部信息的反映能力不夠。垂直測道年變動態明顯，第一個電極埋深為40 m，位于水位動態變化層位，受地下水位變動影響較大，掩埋了深部地層信息。由于垂直測道頂端供電電極靠近地表，因而垂直測道不能較好地抑制來自地表的干擾，建議取消垂直測道深井觀測。為盡可能記錄到與孕震有關的信息，大柏舍臺深井地電阻率應適當增加觀測極距和裝置埋深，水平測道極距AB/2取250～350 m、埋深H>150 m能滿足需求；在條件允許時，可進一步增加裝置埋深以更好地獲取深部信息。考慮到地電阻率變化與主應力方向有關，建議增加NS向水平測道深井觀測。

[1] 錢家棟，陳有發，金安忠. 地電阻率法在地震預報中的應用[M]. 北京：地震出版社，1985：83-103．

Qian Jiadong, Chen Youfa, Jin Anzhong. The Apparent Resistivity Method Used in Earthquake Prediction[M]. Beijing: Seismological Press, 1985: 83-103.

[2] 杜學彬，葉青，馬占虎，等. 強地震附近電阻率對稱四極觀測的探測深度[J]. 地球物理學報，2008，51(6)：1943-1949．

Du Xuebin, Ye Qing, Ma Zhanhu, et al. The Detection Depth of Symmetric Four-Electrode Resistivity Observation Innear the Epicentral Region of Strong Earthquakes[J]. Chinese Journal of Geophysics, 2008, 51(6): 1943-1949.

[3] 錢復業，趙玉林，于謀明，等. 地震前地電阻率異常變化[J]. 中國科學：B輯，1982(9)：831-839．

Qian Fuye, Zhao Yulin, Yu Mouming, et al. Geo-Resistivity Anomaly Before Earthquake[J]. Science in China: Series B, 1982 (9): 831-839.

[4] 錢復業，趙玉林，劉婕，等. 唐山7.8級地震地電阻率臨震功率譜異常[J]. 地震，1990(3)：33-38．

Qian Fuye, Zhao Yulin, Liu Jie, et al. Power Spectrum Anomaly of Earth Resistivity Immediately Before TangshanMs7.8 Earthquake[J]. Earthquake, 1990 (3): 33-38

[5] 桂燮泰，關華平，戴經安. 唐山、松潘地震前視電阻率短臨異常圖像重現性[J]. 西北地震學報，1989，11(4)：71-75.

Gui Xietai, Guan Huaping, Dai Jing’an. The Short -Term and Immediate Anomalous Pattern Recurrences of the Apparent Resistivity Before the Tangshan and Song Pan Earthquakes of 1976[J]. Northwestern Seismological Journal, 1989, 11(4): 71-75.

[6] Lu J, Qian F Y, Zhao Y L. Sensitivity Analysis of the Schlumberger Monitoring Array: Application to Changes of Resistivity Proior to the 1976 Rathquake in Tangshan, China[J]. Tectonophysics, 1999, 307(3): 397-405.

[7] Lu J, Xue S Z, Qian F Y, et al. Unexpected Changes in Resistivity Monitoring for Earthquakes of the Longmenshan in Sichuan China with a Fixed Schlumberger Sounding Array[J]. Pepi, 2004, 145(1/2/3/4): 87-97.

[8] 張學民，李美，關華平. 汶川8.0級地震前的地電阻率異常分析[J]. 地震，2009，29(1)：108-115.

Zhang Xuemin, Li Mei, Guan Huaping. Anomaly Analysis of Earth Resistivity Observations Before the Wenchuan Earthquake[J]. Earthquake, 2009, 29 (1): 108- 115.

[9] 張淑亮，李艷，王霞，等. 亞失穩態野外觀測證據[J]. 吉林大學學報(地球科學版)，2016，46(2)：603-609.

Zhang Shuliang, Li Yan, Wang Xia, et al. Outdoor Observation Evidence of the Meta-Instability State[J]. Journal of Jilin University (Earth Science Edition), 2016, 46(2): 603-609.

[10] 葉青，杜學彬，陳軍營，等. 2003年大姚和民樂--山丹地震1年尺度預測[J]. 地震研究，2005，28(3)：226-230.

Ye Qing, Du Xuebin, Chen Junying, et al. One-Year Prediction for the Dayao and Minle-Shandan Earthquakes in 2003[J]. Journal of Seismological Research, 2005, 28 (3): 226-230.

[11] 杜學彬. 在地震預報中的兩類視電阻率變化[J]. 中國科學：地球科學，2010，40(10)：1321-1330.

Du Xuebin. Two Types of Changes in Apparent Resistivity in Earthquake Prediction[J]. Scientia Sinica: Terrae, 2010, 40(10): 1321-1330.

[12] 衛定軍，李春貴. 搭建蔬菜大棚對固原地電臺視電阻率觀測影響的分析[J]. 地震地磁觀測與研究，2009，30(2)：77-82.

Wei Dingjun, Li Chungui. Analysis about Influence of Green House on Apparent Resistivity Observation at Guyuan Seismic Station[J]. Seismological and Geomagnetic Observation and Research, 2009, 30(2): 77-82.

[13] 解滔，盧軍，李美，等. 地埋鋼纜對寶昌臺地電阻率干擾的定量分析[J]. 地球物理學進展，2013，28(2)：727-734.

Xie Tao, Lu Jun, Li Mei, et al. Quantitative Analysis of Disturbance Caused by Burned Wirerope in Baochang Geoelectric Resistivity Station[J]. Progress in Geophys, 2013, 28(2): 727-734.

[14] 石富強，張國強，方煒，等. 陜西周至地電臺地電阻率年變特征分析[J]. 地震學報，2014，36(6)：1113-1123.

Shi Fuqiang, Zhang Guoqiang, Fang Wei, et al. Annual Variation Characteristics of Georesistivity at Zhouzhi Geoelectric Station Shaanxi Province[J]. Acta Seismologica Sinica, 2014, 36 (6) :1113-1123.

[15] 王幫本，劉永蘭，李驗軒. 深井地電阻率與地震預報[J].地震研究，1981，4(1)：378-403.

Wang Bangben, Liu Yonglan, Li Yanxuan. The Elecrtic Resisitivity in the Boreholes and the Earthquake Prediction[J]. Journal of Seismological Research, 1981, 4(1): 378-403.

[16] 蘇鸞聲，王邦本，夏良苗，等. 井下電極觀測地電阻率排除地面干擾的實驗[J]. 地震學報，1982，4(3)：274-276.

Su Luansheng, Wang Bangben, Xia Liangmiao, et al. Elimination of Surface Disturbances in Earth-Resistivity Measurement by Lowering the Electrodes in Shallow Wells[J]. Acta Seismologica Sinica, 1982, 4(3): 274-276.

[17] 劉允秀，吳國有，王蕃樹，等. 深埋電極地電阻率觀測的實驗結果[C]// 地震預測：地電方法論文集. 福州：福建科學技術出版社，1985：206-216.

Liu Yunxiu, Wu Guoyou, Wang Fanshu, et a1. Test Result for Buried Electrode Resistivity Observation[C]// Collect Works of Earthquake Prediction: Geoelectrical Methods. Fuzhou: Science and Technology Press of Fujian, 1985: 206-216.

[18] 劉昌謀，桂燮泰，柴劍勇，等. 河源地電臺全空間地電阻率試驗[J]. 華南地震，1994，14(3)：40-45.

Liu Changmou, Gui Xietai, Chai Jianyong, et al. The Observation Test in Deep-Hole Electrodes (Whole Space) Resistivity at the Heyuan Geoelectrical Observatory[J]. South China Journal of Seismology, 1994, 14(3): 40-45.

[19] 田山，劉允秀，聶永安，等.地震地電阻率觀測改進方法研究:電測井技術的移植應用與數值模型分析[J].地震學報，2009，31(3)：272-281.

Tian Shan, Liu Yunxiu, Nie Yong’an, et a1. Improved Method of Earth Resistivity Observation: Transplant Application and Numerical Model Analysis of Electric Well Logging Technique[J]. Acta Seismologica Sinica, 2009, 31(3) : 272-281.

[20] 聶永安，姚蘭予. 成層半空間深埋電極產生的電位分布[J].中國地震，2009，25(3)：246-225.

Nie Yong’an, Yao Lanyu. Study on Electrical Potential by Buried Source Electrode Within Horizontally Layered Half-Space Model[J]. Earthquake Reasearch in China, 2009, 25(3): 246-255.

[21] 聶永安，巴振寧，聶瑤. 深埋電極的地電阻率觀測研究[J]. 地震學報，2010，32(1)：33-40.

Nie Yong’an, Ba Zhenning, Nie Yao. Study on Buried Electrode Resistivity Monitoring System[J]. Acta Seismologica Sinica, 2010 , 32(1) :33 -40.

[22] 解滔，杜學彬，陳軍營，等. 井下地電阻率觀測中地表電流干擾影響計算[J]. 地球物理學進展，2012，27(1)：112-121.

Xie Tao, Du Xuebin, Chen Junying, et a1. Calculation for the Influence from the Surface Disturbance Current in the Deep-Well Geoelectrical Resistivity Observation[J]. Progress in Geophys, 2012, 27(1) : 112-121.

[23] 解滔，杜學彬，盧軍. 井下視電阻率觀測影響系數的分析[J]. 中國地震，2016，32(1)：40-53.

Xie Tao, Du Xuebin, Lu Jun. Sensitivity Coefficients Analysis of Deep-Well Apparent Resistivity Measurement[J]. Earthquake Reasearch in China, 2016, 32(1): 40-53.

[24] 解滔，盧軍. 地電阻率三維影響系數及其應用[J]. 地震地質，2015，37(4): 1125-1135.

Xie Tao, Lu Jun. Three-Dimensional Sensitivity Coefficients of Apparent Resistivity and Preliminary Application[J]. Seismology and Geology, 2015, 37(4): 1125-1135.

[25] 張磊，喬子云，羅娜，等. 河北大柏舍臺深井與淺層地電阻率觀測對比分析[J]. 華北地震科學，2015，33(4)：49-53.

Zhang Lei, Qiao Ziyun, Luo Na, et al. Contrastive Analysis of Georesistivity in Deep-Well and on Ground at Dabaishe Station [J]. North China Earthquake Sciences, 2015, 33(4): 49-53.

[26] 趙玉林，盧軍，張洪魁，等. 電測量在中國地震預報中的應用[J]. 地震地質，2001，23(2)：277-285.

Zhao Yulin, Lu Jun, Zhang Hongkui, et al. The Application of Electrical Measurements to Earthquake Prediction in China[J]. Seismology and Geology, 2001, 23(2): 277-285.

[27] 張國苓，喬子云，賈立峰，等. 隆堯地電阻率與地下水位關系分析[J]. 地震地磁觀測與研究，2013，34(5/6)：141-143.

Zhang Guoling, Qiao Ziyun, Jia Lifeng, et al. The Study of Relationship Between Earth Resistivity and Ground Water Level at Longyao Seismic Station[J]. Seismological and Geomagnetic Observation and Research, 2013, 34(5/6): 141-143.

[28] 錢家棟，曹愛民. 1976年唐山7.8級地震地電阻率和地下水前兆綜合物理機制研究[J]. 地震，1998，18(增刊1)：1-9.

Qian Jiadong, Cao Aimin. Research on Comprehensive Mechanism of Precursors in Apparant Resistivity and Water Table Associated with 1976 Tangshan Earthquake (M=7.8)[J]. Earthquake, 1998, 18(Sup. 1): 1-9.

[29] Park S K, Van G P. Inversion of Pole-Pole Data for 3-D Resistivity Structure Beneath Arrays of Electrodes[J]. Geophysics, 1991, 56(7): 951-960.

[30]Seigel H O. Mathematical Formulation and Type Curves for Induced Polarization[J]. Geophysics, 1959, 24(3): 547-565.