石富強(qiáng) 邵輝成 張國(guó)強(qiáng) 方 煒

（中國(guó)西安710068陜西省地震局）

引言

地電阻率觀測(cè)作為一項(xiàng)重要的地震前兆觀測(cè)手段，受到國(guó)內(nèi)外大量學(xué)者關(guān)注（Scholz et al，1973；Park et al，1993；Lu et al，1999，2004；Yang et al，2002；Balasco et al，2004）.在實(shí)驗(yàn)研究中也發(fā)現(xiàn)了被認(rèn)為是地震前兆的地電阻率異常信息（Brace，Orange，1968；Wang et al，1975，1978；陳峰等，2000；郝錦綺等，2002）.我國(guó)自1966年邢臺(tái)地震后開展地電阻率連續(xù)觀測(cè)以來，已經(jīng)形成了由近百個(gè)定點(diǎn)連續(xù)觀測(cè)臺(tái)站組成的地電阻率觀測(cè)網(wǎng).在觀測(cè)網(wǎng)內(nèi)也監(jiān)測(cè)到過中強(qiáng)地震的可能前兆異常信息，并開展了一年尺度的預(yù)測(cè)實(shí)踐（葉青等，2005）.長(zhǎng)期的監(jiān)測(cè)預(yù)報(bào)工作積累了大量震例總結(jié)（錢復(fù)業(yè)，趙玉林，1980；錢家棟，1993；Du et al，2000；杜學(xué)彬，2010），也認(rèn)識(shí)到了地震前兆機(jī)制的復(fù)雜性.結(jié)合地下介質(zhì)各向異性以及震源特征，國(guó)內(nèi)外學(xué)者開展了大量關(guān)于地電阻率地震前兆異常的探索研究工作（Qian et al，1996；Busby，2000；杜學(xué)彬等，2007，2008；）.

近年來，隨著我國(guó)農(nóng)村經(jīng)濟(jì)結(jié)構(gòu)的快速轉(zhuǎn)型以及城市化進(jìn)程的高速發(fā)展，過去多年定點(diǎn)連續(xù)觀測(cè)的很多地電臺(tái)都出現(xiàn)了各種各樣的干擾，歸納起來主要分為8類（汪志亮等，2002）：① 儀器系統(tǒng)干擾；② 觀測(cè)系統(tǒng)漏電干擾；③工業(yè)游散電流隨機(jī)干擾；④ 高壓線路干擾；⑤風(fēng)擾或磁擾；⑥地埋金屬管線干擾；⑦大型水庫(kù)和蓄水干擾；⑧年變干擾.針對(duì)這些類型的干擾源，人們采取井下地電觀測(cè)方法（楊興悅等，2012）以及發(fā)展無人值守臺(tái)站和可控源觀測(cè)系統(tǒng)等途徑抑制干擾影響.田山等（2009）提出將三向石油測(cè)井技術(shù)移植到地震地電阻率觀測(cè)來提高地電阻率的抗干擾能力；解滔等（2012）研究了布極區(qū)地表、井下地電觀測(cè)中抑制地表雜散電流的干擾影響；方煒等（2010）結(jié)合實(shí)際觀測(cè)理論計(jì)算了高壓直流輸電對(duì)地電場(chǎng)觀測(cè)的影響；張世中等（2013）通過測(cè)試分析，研究了城市軌道交通對(duì)地電阻率觀測(cè)的影響；安張輝等（2010，2011）闡述了希爾伯特－黃變換（HHT）方法在地電觀測(cè)數(shù)據(jù)處理中的可靠性和應(yīng)用前景，并利用該方法分析了城市軌道交通對(duì)地電觀測(cè)的干擾；解滔等（2013）利用數(shù)值模擬方法定量分析了地埋鋼纜線對(duì)寶昌地電阻率觀測(cè)的影響；張秀霞等（2009）以及衛(wèi)定軍和李春貴（2009）分別介紹了新沂臺(tái)和固原臺(tái)地電阻率觀測(cè)受蔬菜大棚搭建的干擾情況.本文以偶極接地線為例，重點(diǎn)研究布極區(qū)內(nèi)以蔬菜大棚等為代表的偶極接地線對(duì)地表、井下地電阻率觀測(cè)的固定干擾影響.

偶極接地線纜是金屬管網(wǎng)類干擾的一種，主要表現(xiàn)為布極區(qū)內(nèi)以蔬菜大棚等為代表的兩端與大地連接，中間架空的一類固定干擾源.其金屬線一般較長(zhǎng)（從幾十米到幾百米不等），橫貫地表，兩端直接與大地相連，中間輔以混凝土柱、木棒等支撐，如圖1所示.這類干擾使得地電阻率連續(xù)觀測(cè)出現(xiàn)下降、上升或年變形態(tài)畸變并隨著干擾的持續(xù)呈現(xiàn)趨勢(shì)性變化，表現(xiàn)出來的形態(tài)與地震前兆異常形態(tài)（杜學(xué)彬，2010）類似.同時(shí)在異常調(diào)查落實(shí)中也存在著一些矛盾的現(xiàn)象.以偶極接地線干擾為例，在江蘇新沂地電臺(tái)（張秀霞等，2009）表現(xiàn)為“下降”，而在寧夏固原地電臺(tái)（衛(wèi)定軍，李春貴，2009）卻表現(xiàn)為“上升”，并且同一觀測(cè)點(diǎn)不同測(cè)道的變化形態(tài)和大小也不盡相同.此外這類干擾還在北京延慶臺(tái)、天津?qū)氎媾_(tái)、河北昌黎臺(tái)、遼寧新城子臺(tái)以及陜西周至臺(tái)等均有出現(xiàn)①據(jù)2013年全國(guó)電磁會(huì)商會(huì)相關(guān)單位資料報(bào)告..現(xiàn)場(chǎng)實(shí)驗(yàn)結(jié)果表明，平行測(cè)線埋設(shè)的地下管線對(duì)地電阻率觀測(cè)的影響最為顯著，而斜交或垂直敷設(shè)的次之或不明顯（汪志亮等，2002）.目前針對(duì)偶極接地線纜這類特殊干擾影響地電阻率觀測(cè)的定量分析研究的報(bào)道尚不多見.

圖1 測(cè)區(qū)內(nèi)金屬線纜布設(shè)情況（a）及其簡(jiǎn)化模型（b）Fig.1 The layout of the metallic cables（a）and its simplified model（b）in surveying area

本文將日常觀測(cè)中存在的布極區(qū)內(nèi)以蔬菜大棚金屬構(gòu)架為代表的固定干擾源簡(jiǎn)化為偶極接地裝置，建立了偶極接地線干擾地電阻率觀測(cè)的耦合物理模型，利用數(shù)值模擬方法系統(tǒng)地分析計(jì)算了偶極接地線干擾地電阻率測(cè)量的干擾機(jī)制，并與現(xiàn)有觀測(cè)資料進(jìn)行對(duì)比驗(yàn)證分析.該分析結(jié)果對(duì)今后認(rèn)識(shí)類似干擾提供了可靠的分析處理方法.

1 物理模型構(gòu)建

基于數(shù)值模擬技術(shù)的地震電信號(hào)研究已成為解釋地震電信號(hào)資料，探索地震電信號(hào)產(chǎn)生機(jī)制，以及影響因素分析的新方法（Huang，Lin，2010；黃清華，林玉峰，2010）.本文旨在分析偶極接地線這種特定干擾對(duì)地電阻率觀測(cè)的影響，為此在建模時(shí)沒有考慮自然電場(chǎng)的產(chǎn)生機(jī)制，而是將地電測(cè)區(qū)等效為內(nèi)部無初始電場(chǎng)分布的水平3層導(dǎo)電介質(zhì)：第一層厚度為h1，電阻率ρ1；第二層厚度為h2，電阻率ρ2；第三層電阻率ρ3.A和B為供電電極，M和N為測(cè)量電極，對(duì)稱分布于x軸上，且4個(gè)電極埋深均為h，如圖2a所示.設(shè)PQ為地電測(cè)區(qū)一組偶極接地線，跨度為L(zhǎng)；P和Q為兩端接地點(diǎn)，P點(diǎn)的位置參數(shù)為x，到測(cè)線的垂直距離為s.接地線PQ走向的延長(zhǎng)線與x軸的正方向夾角為α（0°≤α≤90°）.令線纜兩端接地電阻分別為R1和R2，線纜自身電阻為Rf（圖2b）.因此一組接地線干擾系統(tǒng)可以簡(jiǎn)化為兩端接地的等效電阻體R0.定義R0為地電阻率觀測(cè)的干擾電阻，則

圖2 水平3層介質(zhì)中對(duì)稱四極裝置（剖面）（a）和偶極接地線模型示意圖（俯視）（b）圖（b）中左上角（b1）小圖為其右下角（b2）小圖中接地線PQ的細(xì)節(jié)結(jié)構(gòu)示意圖Fig.2 （a）Schlumberger array of apparent resistivity observation in horizontally inhomogeneous three-layered medium；（b）The physical model of two-end grounded cable.And the part（b1）in the upper-left corner is the details of grounded cable PQshown in part（b2）in the lower-right corner

模擬過程采用直流恒壓供電方式，在地電阻率正常觀測(cè)時(shí)，通過測(cè)量M與N 之間的電勢(shì)差UMN以及供電電流I，地電阻率由公式ρs＝K（UMN／I）計(jì)算可得，其中K為裝置系數(shù).當(dāng)接地線干擾發(fā)生后，測(cè)量電極M與N 之間的電勢(shì)差為U′MN，回路中的供電電流為I′.由于極距固定，電極埋深固定，所以裝置系數(shù)K不變.則干擾后實(shí)際測(cè)量的地電阻率為?ρs＝K（U′MN／I′）.令η＝?ρs／ρs為接地線干擾地電阻率觀測(cè)的影響系數(shù)，則

當(dāng)η＞1時(shí)，偶極接地線對(duì)地電阻率的干擾表現(xiàn)為增加；當(dāng)η＜1時(shí)，偶極接地線對(duì)地電阻率的干擾表現(xiàn)為減??；當(dāng)η＝1時(shí)，偶極接地線對(duì)地電阻率的影響可以忽略不計(jì).

2 模擬計(jì)算及實(shí)例分析

ANSYS軟件是美國(guó)ANSYS公司研制的大型的集結(jié)構(gòu)、流體、電場(chǎng)、磁場(chǎng)、聲場(chǎng)分析于一體的通用有限元模擬軟件.在航空航天、機(jī)械土木、輕工日化、能源化工、地礦水利等領(lǐng)域有著廣泛的應(yīng)用.根據(jù)所研究問題的情況，建模時(shí)可以在軟件內(nèi)選擇對(duì)應(yīng)的單元類型，如模擬桁架結(jié)構(gòu)的LINK10單元、模擬殼體變形的SHELL43單元、模擬力－電－磁-熱的SOLID5單元、模擬電路的CIRCUIT124單元等百余種單元類型.本文根據(jù)圖2a，利用ANSYS電磁場(chǎng)耦合實(shí)體單元（SOLID5）建立了水平層狀結(jié)構(gòu)的三維有限元計(jì)算模型.在全國(guó)100多個(gè)地電臺(tái)站中，大多數(shù)臺(tái)站的供電極極距約為1.2km，為了有效地減小截?cái)噙吔鐚?duì)地下電場(chǎng)線分布的影響（湯井田，公勁喆，2010），本文經(jīng)過反復(fù)試驗(yàn)，選取大?。ㄩL(zhǎng)×寬×高）為10km×10km×5km的模型模擬無限半空間的地下介質(zhì)（關(guān)于幾何尺寸對(duì)模擬結(jié)果穩(wěn)定性的影響，解滔等（2013）作了詳細(xì)的論證，本文模擬試驗(yàn)過程與之類似，不再贅述）.該模型在縱向分為3層，分別表示分層介質(zhì)不同層的電性屬性.以常見的4種3層電性結(jié)構(gòu)為例，本文采用AB＝1 200m的對(duì)稱四極溫納裝置（MN＝400m）.其中A，B，M，N 點(diǎn)的坐標(biāo)分別為（－600m，0，h）、（600m，0，h）、（－200m，0，h）和（200m，0，h），供電采用偶極直流電源，在供電電極A和B兩點(diǎn)分別施加荷載電壓UA＝100V、UB＝－100V.求解給出電勢(shì)差ΔUMN和電流I的正常理論值；然后利用ANSYS電路單元（CIRCUIT124）創(chuàng)建干擾電阻R0，將R0兩端耦合于已建立的電性結(jié)構(gòu)模型表面上的P點(diǎn)和Q點(diǎn)來模擬偶極接地線，耦合自由度為電壓，同時(shí)選取干擾電阻跨度PQ＝200m；通過改變接地線位置參數(shù)s，α，x及電極埋深h，求解給出干擾后測(cè)量電極M與N之間的電位差ΔU′MN和流過供電電極A 和B 兩點(diǎn)的電流I′；進(jìn)而根據(jù)式（2）計(jì)算出偶極接地線對(duì)地電阻率觀測(cè)的影響.本文計(jì)算所采用模型的電性結(jié)構(gòu)參數(shù)見表1.

2.1 影響系數(shù)η隨位置x的變化

以地表觀測(cè)為例：取電極埋深h＝2m；根據(jù)對(duì)稱性簡(jiǎn)化取方位角α＝0°，45°和90°；假定接地線分布于測(cè)線附近（s＝0m）且干擾電阻R0＝20Ω·m.通過平移改變位置參數(shù)x模擬偶極接地線對(duì)地電阻率的影響.模擬計(jì)算結(jié)果如圖3所示.

表1 模型的電性結(jié)構(gòu)參數(shù)Table 1 Physical parameters of the electric structures for the three-layered model

圖3 影響系數(shù)η隨位置參數(shù)x的變化關(guān)系.（a）α＝0°；（b）α＝45°；（c）α＝90°Fig.3 Relationship between the influence coefficientηand positional parameter x（a）α＝0°；（b）α＝45°；（c）α＝90°

由圖3可見，偶極接地線對(duì)地電阻率的影響主要源自于測(cè)區(qū)內(nèi)的干擾，當(dāng)接地線位于測(cè)線兩端的測(cè)區(qū)外時(shí)（接地線PQ兩端點(diǎn)遠(yuǎn)離供電極約200m），這種影響可以忽略不計(jì).當(dāng)接地線平行于測(cè)線時(shí)（α＝0°），對(duì)地電阻率的干擾主要表現(xiàn)為增加，且增加幅度較大，最大可達(dá)到10%的水平；當(dāng)接地線垂直于測(cè)線時(shí)（α＝90°），對(duì)地電阻率的干擾主要表現(xiàn)為減小，且減小幅度較小，最大不超過1%的水平；當(dāng)接地線與測(cè)線斜交時(shí)，接地線對(duì)地電阻率的干擾會(huì)表現(xiàn)出隨電極分布增大、減小交替出現(xiàn)的情況，變化幅度約為5%的水平.這與汪志亮等（2002）現(xiàn)場(chǎng)實(shí)驗(yàn)結(jié)果在定性趨勢(shì)上吻合.對(duì)比不同電性結(jié)構(gòu)的模擬結(jié)果發(fā)現(xiàn)，偶極接地線對(duì)地電阻率觀測(cè)的影響是增加還是減小，主要取決于其相對(duì)測(cè)線的位置x以及方位角α，與測(cè)區(qū)電性結(jié)構(gòu)關(guān)系不大.這是因?yàn)榻拥鼐€搭建后，線纜會(huì)將地表P點(diǎn)與Q點(diǎn)導(dǎo)通，人工電場(chǎng)的分布發(fā)生變化，進(jìn)而影響地電阻率ρs在原裝置下的測(cè)量值.

電性結(jié)構(gòu)的差異性決定干擾幅度的大小，對(duì)于H型和Q型這兩類表層電阻較大的電性結(jié)構(gòu)而言，電流主要分布于地下低阻層，表層電流密度較小.而測(cè)量又在地表進(jìn)行，因此當(dāng)偶極接地線將地表的P點(diǎn)與Q點(diǎn)導(dǎo)通后，地表微弱的對(duì)稱性電場(chǎng)分布將對(duì)干擾非常敏感，發(fā)生大幅度調(diào)整，使得地電阻率觀測(cè)值發(fā)生較大幅度的變化；反之，對(duì)于A型和K型這兩類下伏高阻層的電性結(jié)構(gòu)而言，表層電流密度較大，對(duì)干擾出現(xiàn)的敏感度很小，故地電阻率變化幅度較小.同時(shí)從另一個(gè)角度也可看出，A型和K型電性斷面相對(duì)于H型和Q型能夠較好地抑制來自地表的干擾，這與解滔等（2012）對(duì)地表電流干擾的分析結(jié)果相同.

2.2 影響系數(shù)η隨電極埋深h的變化

同樣以表1中列舉的4種電性結(jié)構(gòu)為例，取接地線方位角α＝0°，45°和90°，干擾電阻R0＝20Ω·m，跨度為200m.通過改變測(cè)量電極M，N和供電電極A，B的埋深h來模擬偶極接地線對(duì)地震地電阻率觀測(cè)的影響.由圖3可知，接地線對(duì)地電阻率觀測(cè)的影響的最大值集中在供電電極和測(cè)量電極附近，因此本文僅取位置參數(shù)x＝200m（P點(diǎn)在測(cè)量電極N處），且線纜位于測(cè)線附近（s＝0m）進(jìn)行討論分析.模擬計(jì)算結(jié)果如圖4所示.

圖4顯示了數(shù)值模擬給出的影響系數(shù)隨電極埋深的變化關(guān)系.總體而言，增加電極埋深可以顯著減小地表接地線搭建對(duì)地電阻率觀測(cè)的影響.合理的電極埋深選擇取決于電性結(jié)構(gòu)的差異性.對(duì)于H型和Q型這兩類下伏低阻層的電性斷面而言，在第一層介質(zhì)中（h＜15m）增大電極埋深便可顯著地減小線纜搭建對(duì)地電阻率干擾的影響系數(shù)，當(dāng)埋深穿透第一層介質(zhì)時(shí)，降幅可達(dá)10%左右；繼續(xù)增大電極埋深到第二層介質(zhì)，地表接地線對(duì)地電阻率觀測(cè)的影響將繼續(xù)降低，但影響系數(shù)的衰減速率急劇減小.因此，對(duì)于H型和Q型電性斷面而言，盡管在第一層介質(zhì)中增大電極埋深可以顯著減小接地線纜的影響，但是由于其引起的干擾幅度較大，當(dāng)電極埋設(shè)于第一層介質(zhì)的底面還沒有達(dá)到預(yù)期效果的時(shí)候，需要在第二層介質(zhì)中大幅度增大電極深埋，甚至需要埋到第三層介質(zhì)中.而對(duì)于A型和K型兩類表層電阻較小的電性斷面而言，盡管在第一層介質(zhì)（h＜15m）中增大電極埋深時(shí)，影響系數(shù)變化速度緩慢不明顯，降幅約為0.5%，但是在第二層介質(zhì)中增加電極埋深，影響系數(shù)衰減速率將迅速加大；加之，接地線對(duì)A型和K型結(jié)構(gòu)的電性斷面地電阻率觀測(cè)影響幅度較小，因此在第二層介質(zhì)中增加電極埋深便可顯著提高地電阻率連續(xù)觀測(cè)的質(zhì)量和抗干擾能力.

圖4 影響系數(shù)η隨電極埋深h的變化關(guān)系.（a）H型；（b）Q型；（c）A型；（d）K型Fig.4 Relationship between the influence coefficientηand buried depth h（a）H-type；（b）Q-type；（c）A-type；（d）K-type

2.3 實(shí)例分析

衛(wèi)定軍和李春貴（2009）以及張秀霞等（2009）分別根據(jù)相關(guān)臺(tái)站的異常落實(shí)工作介紹了接地金屬線纜干擾地電阻率觀測(cè)的情況.本文利用已建立的耦合物理模型，通過數(shù)值模擬分析偶極接地線對(duì)固原臺(tái)和新沂臺(tái)地電阻率觀測(cè)的影響.電性結(jié)構(gòu)參數(shù)由擬合核函數(shù)法反演測(cè)深曲線給出.反演結(jié)果表明，新沂臺(tái)臺(tái)址下方為3層K型斷面介質(zhì)，而固原臺(tái)臺(tái)址下方為4層HK型斷面介質(zhì)（圖5）.

圖5 垂直向電測(cè)深數(shù)據(jù)解釋.（a）新沂臺(tái)；（b）固原臺(tái)Fig.5 Interpretation of vertical electric soundings data for the stations Xinyi（a）and Guyuan（b）

新沂臺(tái)電極埋深h＝3m，極距AB＝1 000m，MN＝240m，臺(tái)址下方介質(zhì)為表層電阻較小的K型斷面.地電測(cè)區(qū)內(nèi)于2007年9—11月開始搭建偶極接地線，搭建區(qū)域如圖6a所示.每根線長(zhǎng)約100m，主要分布于NE測(cè)道的供電極A2與測(cè)量極M2之間，呈EW走向.張秀霞等（2009）指出，2007年8月后新沂臺(tái)N45°E和EW向受到線間漏電影響，觀測(cè)數(shù)據(jù)無法反映真實(shí)變化.2007年12月21日更換3個(gè)測(cè)向的埋地電纜，更換前后只有NS向數(shù)據(jù)未發(fā)生變化.所以他們并沒有給出這兩道干擾的影響量，本文沒有對(duì)此作進(jìn)一步分析，僅以NS向?qū)嶋H測(cè)值為對(duì)比研究對(duì)象.經(jīng)實(shí)測(cè)新沂臺(tái)線纜接地電阻約為1.7—6.0Ω（張秀霞等，2009），本文近似取接地電阻為R1＝R2＝1.5Ω.利用圖5a反演給出的參數(shù)建立耦合物理模型，模擬給出接地線對(duì)新沂臺(tái)NS向地電阻率的影響系數(shù)，模擬結(jié)果與實(shí)際觀測(cè)結(jié)果的對(duì)比情況見表2.

圖6 測(cè)區(qū)電極和金屬線分布示意圖.（a）新沂臺(tái)；（b）固原臺(tái)Fig.6 Schematic diagram of the arrangement of electrodes and metallic cables for the stations Xinyi（a）and Guyuan（b）in the surveying area

表2 新沂臺(tái)和固原臺(tái)偶極接地線干擾地電阻率觀測(cè)的模擬結(jié)果與實(shí)測(cè)結(jié)果對(duì)比Table 2 Comparison of the simulated and observed georesistivities at the stations Xinyi and Guyuan under the influence of two-end grounded cables

固原臺(tái)位于山坡前，臺(tái)址下方為表層電阻較大的4層HK型斷面，極距為AB＝1 000m，MN＝200m.地電測(cè)區(qū)內(nèi)2008年3月22日開始搭建接地線纜，搭建區(qū)域如圖6b所示.總體上線一端固定于供電電極A2附近，另一端固定于測(cè)量電極M2附近，與NE向地電阻率測(cè)線近乎平行，線纜長(zhǎng)約350m.由于固原地表干燥且表層電阻大，本文近似取干擾電阻R0＝25Ω，由圖5b反演參數(shù)建立耦合物理模型，模擬給出接地線對(duì)固原地電阻率觀測(cè)的影響系數(shù)，模擬結(jié)果與實(shí)際觀測(cè)結(jié)果對(duì)比情況見表2.

經(jīng)對(duì)比發(fā)現(xiàn)，對(duì)于新沂臺(tái)NS測(cè)項(xiàng)以及固原臺(tái)NS測(cè)項(xiàng)和EW測(cè)項(xiàng)，模擬計(jì)算結(jié)果與實(shí)測(cè)結(jié)果吻合良好，能夠?yàn)槭卟舜笈?、鐵絲網(wǎng)等偶極接地裝置干擾地電阻率觀測(cè)的定量計(jì)算以及干擾排查等提供有力依據(jù).但模擬結(jié)果對(duì)固原臺(tái)NE向預(yù)測(cè)偏大，這主要是由于地下電性結(jié)構(gòu)的各項(xiàng)異性差異造成的.本文計(jì)算采用的電性結(jié)構(gòu)是由EW向測(cè)深反演給出.由圖5b可以看出，NE向底層存在高阻層，而根據(jù)EW向反演給出的電性結(jié)構(gòu)并不能反映NE向底層的高阻結(jié)構(gòu).因此，對(duì)于NE向地電阻率的實(shí)際觀測(cè)，由于下伏高阻層的存在，人工電流將趨于地表，使得地表電流密度大于本文模型計(jì)算結(jié)果；進(jìn)而降低了其對(duì)地表干擾感應(yīng)的靈敏度，抑制了干擾的大幅度增加.故本文模型預(yù)測(cè)值偏大.為此，我們以NE向測(cè)深反演結(jié)果建立干擾物理模型，其它參數(shù)和條件不變.模擬結(jié)果表明，蔬菜大棚對(duì)固原臺(tái)NE向地電阻率觀測(cè)的影響約為8%.綜上，本文建立的耦合物理模型不但解釋了“同為蔬菜大棚干擾，新沂臺(tái)表現(xiàn)為‘下降’，而固原臺(tái)表現(xiàn)為‘上升’”，并且在定量上也吻合良好.

3 討論與結(jié)論

農(nóng)村經(jīng)濟(jì)結(jié)構(gòu)的轉(zhuǎn)型和城市化進(jìn)程的發(fā)展給地震地電阻率觀測(cè)帶來的影響已不容忽視，探索新的、基于物理模型的定點(diǎn)連續(xù)觀測(cè)地電阻率數(shù)據(jù)處理方法對(duì)于異常判別和地震監(jiān)測(cè)預(yù)報(bào)實(shí)踐顯得更加重要.本文將臺(tái)站下方介質(zhì)簡(jiǎn)化為水平3層介質(zhì)，將偶極接地線簡(jiǎn)化為電阻體，建立耦合物理模型，就偶極接地線這類特定干擾作了比較詳細(xì)的模擬研究，并結(jié)合實(shí)際觀測(cè)資料作了對(duì)比分析研究，從理論上分析給出了偶極接地線對(duì)定點(diǎn)連續(xù)觀測(cè)單極距地電阻率觀測(cè)的影響機(jī)制，主要包括：

1）本文根據(jù)實(shí)際調(diào)查和理論分析建立了偶極接地線干擾地電阻率觀測(cè)的耦合物理模型，并結(jié)合常見電性結(jié)構(gòu)綜合分析了線纜分布、線纜方位角以及電極埋深等對(duì)地電阻率的干擾.結(jié)果表明，偶極接地線對(duì)地電阻率的干擾是增大還是減小主要由地表電勢(shì)的重分布造成，與線纜分布的位置、方位角等有關(guān)，電性結(jié)構(gòu)的差異決定干擾幅度的大小.該結(jié)果對(duì)于相關(guān)地震前兆觀測(cè)臺(tái)站的異常判別分析以及臺(tái)站改造具有一定的指導(dǎo)意義.

2）根據(jù)圖4模擬結(jié)果可知，不同電性結(jié)構(gòu)下地表干擾隨電極埋深的變化關(guān)系也不盡一致.在第一層介質(zhì)增加電極埋深，當(dāng)表層電阻較大時(shí)，地表接地線干擾的影響系數(shù)快速衰減；而當(dāng)表層電阻較小時(shí)，影響系數(shù)的衰減速率很慢，甚至不明顯.因此在采用深埋電極觀測(cè)時(shí)，避免盲目參考借鑒，需要根據(jù)當(dāng)?shù)貙?shí)際的電性斷面選擇合適的電極埋深和避讓距離可以有效地減小接地線搭建對(duì)地電阻率的干擾.

3）在實(shí)例分析中的物理模型參數(shù)由擬合核函數(shù)法反演給出，盡管忽略了實(shí)際水平向不均勻性，但從模擬結(jié)果來看，本文基于測(cè)深反演結(jié)果，由電磁耦合單元的塊導(dǎo)體和電路單元電阻體耦合建立的水平均勻?qū)訝罱橘|(zhì)耦合物理模型，能夠反映臺(tái)站下方介質(zhì)的電性性質(zhì).將該模型應(yīng)用于已確定干擾的新沂臺(tái)和固原臺(tái)地電阻率異常，通過對(duì)比發(fā)現(xiàn)，本文建立的耦合物理模型能夠從定量的角度解釋分析受接地線干擾，新沂臺(tái)地電阻率觀測(cè)值趨勢(shì)下降而固原臺(tái)地電阻率觀測(cè)值趨勢(shì)上升的原因.同時(shí)也體現(xiàn)出基于物理模型的數(shù)值模擬方法在處理地震地電阻率觀測(cè)中存在的“虛假異?！狈矫娴膬?yōu)勢(shì).

寧夏回族自治區(qū)地震局衛(wèi)定軍、江蘇省地震局李飛和河北省地震局張國(guó)苓為本文提供了相關(guān)資料，在工作中與中國(guó)地震臺(tái)網(wǎng)中心解滔進(jìn)行了有益的討論和交流，審稿專家為本文的修改提出了寶貴的建議.作者在此一并表示誠(chéng)摯的謝意.

安張輝，元麗華，李寧，范瑩瑩，譚大誠(chéng)，陳軍營(yíng)，鄭國(guó)磊.2010.HHT方法在地電場(chǎng)數(shù)據(jù)處理中的應(yīng)用［J］.地球物理學(xué)進(jìn)展，25（2）：525-532.

An Z H，Yuan L H，Li N，F(xiàn)an Y Y，Tan D C，Chen J Y，Zheng G L.2010.Discussion on the application of HHT method to geoelectric field data processing［J］.Progress in Geophysics，25（2）：525-532（in Chinese）.

安張輝，杜學(xué)彬，元麗華，譚大誠(chéng)，陳軍營(yíng).2011.HHT方法在受城市軌道交通干擾地電場(chǎng)觀測(cè)數(shù)據(jù)中的應(yīng)用［J］.地震學(xué)報(bào)，33（2）：243-251.

An Z H，Du X B，Yuan L H，Tan D C，Chen J Y.2011.Application of Hilbert-Huang transform to geo-electric data disturbed by urban track traffic［J］.Acta Seismologica Sinica，33（2）：243-251（in Chinese）.

陳峰，修濟(jì)剛，安金珍，廖椿庭，陳大元.2000.巖石電阻率變化各向異性與微裂隙擴(kuò)展方位實(shí)驗(yàn)研究［J］.地震學(xué)報(bào)，22（3）：310-318.

Chen F，Xiu J G，An J Z，Liao C T，Chen D Y.2000.Research on dependence of resistivity changing anisotropy on microcracks extending in rock with experiment［J］.Acta Seismologica Sinica，22（3）：310-318（in Chinese）.

杜學(xué)彬，李寧，葉青，馬占虎，閆睿.2007.強(qiáng)地震附近視電阻率各向異性變化的原因［J］.地球物理學(xué)報(bào)，50（6）：1802-1810.

Du X B，Li N，Ye Q，Ma Z H，Yan R.2007.A possible reason for the anisotropic changes in apparent resistivity near the focal region of strong earthquake［J］.Chinese Journal of Geophysics，50（6）：1802-1810（in Chinese）.

杜學(xué)彬，葉青，馬占虎，李寧，陳軍營(yíng)，譚大誠(chéng).2008.強(qiáng)地震附近電阻率對(duì)稱四極觀測(cè)的探測(cè)深度［J］.地球物理學(xué)報(bào)，51（6）：1943-1949.

Du X B，Ye Q，Ma Z H，Li N，Chen J Y，Tan D C.2008.The detection depth of symmetric four-electrode resistivity observation in／near the epicentral region of strong earthquakes［J］.Chinese Journal of Geophysics，51（6）：1943-1949（in Chinese）.

杜學(xué)彬.2010.在地震預(yù)報(bào)中的兩類視電阻率變化［J］.中國(guó)科學(xué)：地球科學(xué)，40（10）：1321-1330.

Du X B.2011.Two types of changes in apparent resistivity in earthquake prediction［J］.Science China：Earth Sciences，54（1）：145-156.

方煒，張國(guó)強(qiáng)，邵輝成.2010.高壓直流輸電對(duì)地電場(chǎng)觀測(cè)的影響［J］.地震地質(zhì)，32（3）：434-441.

Fang W，Zhang G Q，Shao H C.2010.Study on the impact of HVDC to geoelectric field observation［J］.Seismology and Geology，32（3）：434-441（in Chinese）.

郝錦綺，馮銳，周建國(guó)，錢書清，高金田.2002.巖石破裂過程中電阻率變化機(jī)理的探討［J］.地球物理學(xué)報(bào)，45（3）：426-434.

Hao J Q，F(xiàn)eng R，Zhou J G，Qian S Q，Gao J T.2002.Study on the mechanism of resistivity changes during rock cracking［J］.Chinese Journal of Geophysics，45（3）：426-434（in Chinese）.

黃清華，林玉峰.2010.地震電信號(hào)選擇性數(shù)值模擬及可能影響因素［J］.地球物理學(xué)報(bào)，53（3）：535-543.

Huang Q H，Lin Y F.2010.Numerical simulation of selectivity of seismic electric signal and its possible influences［J］.Chinese Journal of Geophysics，53（3）：535-543（in Chinese）.

錢復(fù)業(yè)，趙玉林.1980.地震前地電阻率變化十例［J］.地震學(xué)報(bào)，2（2）：186-197.

Qian F Y，Zhao Y L.1980.The examples of changes in earth-resistivity prior to strong earthquakes［J］.Acta Seismologica Sinica，2（2）：186-197（in Chinese）.

錢家棟.1993.與大震孕育過程有關(guān)的地電阻率變化研究［J］.中國(guó)地震，9（4）：341-350.

Qian J D.1993.A study on the changes in geoelectrical resistivity associated with preparatory process of great earthquakes in China［J］.Earthquake Research in China，9（4）：341-350（in Chinese）.

湯井田，公勁喆.2010.三維直流電阻率有限元－無限元耦合數(shù)值模擬［J］.地球物理學(xué)報(bào)，53（3）：717-728.

Tang J T，Gong J Z.2010.3DDC resistivity forward modeling by finite-infinite element coupling method［J］.Chinese Journal of Geophysics，53（3）：717-728（in Chinese）.

田山，劉允秀，聶永安，馬驥，龐群英，杜學(xué)彬.2009.地震地電阻率觀測(cè)改進(jìn)方法研究：電測(cè)井技術(shù)的移植應(yīng)用與數(shù)值模型分析［J］.地震學(xué)報(bào)，31（3）：272-281.

Tian S，Liu Y X，Nie Y A，Ma J，Pang Q Y，Du X B.2009.Improved method of earth resistivity observation：Transplant application and numerical model analysis of electric well logging technique［J］.Acta Seismologica Sinica，31（3）：272-281（in Chinese）.

汪志亮，鄭大林，余索榮.2002.地震地電阻率前兆異常現(xiàn)象［M］.北京：地震出版社：51-79.

Wang Z L，Zheng D L，Yu S R.2002.Geoelectric Resistivity Precursor Anomalies of Earthquake［M］.Beijing：Seismological Press：51-79（in Chinese）.

衛(wèi)定軍，李春貴.2009.搭建蔬菜大棚對(duì)固原地電臺(tái)視電阻率觀測(cè)影響的分析［J］.地震地磁觀測(cè)與研究，30（2）：77-82.

Wei D J，Li C G.2009.Analysis about influence of green house on apparent resistivity observation at Guyuan seismic station［J］.Seismological and Geomagnetic Observation and Research，30（2）：77-82（in Chinese）.

解滔，杜學(xué)彬，陳軍營(yíng)，安張輝，譚大成，范瑩瑩，劉君.2012.井下地電阻率觀測(cè)中地表電流干擾影響計(jì)算［J］.地球物理學(xué)進(jìn)展，27（1）：112-121.

Xie T，Du X B，Chen J Y，An Z H，Tan D C，F(xiàn)an Y Y，Liu J.2012.Calculation for the influence from the surface disturbance current in the deep-well geoelectrical resistivity observation［J］.Progress in Geophysics，27（1）：112-121（in Chinese）.

解滔，盧軍，李美，王中平，王亞麗.2013.地埋鋼纜對(duì)寶昌臺(tái)地電阻率干擾的定量分析［J］.地球物理學(xué)進(jìn)展，28（2）：724-734.

Xie T，Lu J，Li M，Wang Z P，Wang Y L.2013.Quantitative analysis of disturbance caused by burned wirerope in Baochang geoelectric resistivity station［J］.Progress in Geophysics，28（2）：724-734（in Chinese）.

楊興悅，楊立明，康云生，安海靜，郝臻，武銀，王燕.2012.天水地震臺(tái)井下地電觀測(cè)系統(tǒng)介紹及其分析［J］.地震研究，35（1）：92-97.

Yang X Y，Yang L M，Kang Y S，An H J，Hao Z，Wu Y，Wang Y.2012.Introduction and analysis of geo-electrical observation system in underground well at Tianshui station［J］.Journal of Seismological Research，35（1）：92-97（in Chinese）.

葉青，杜學(xué)彬，陳軍營(yíng)，譚大成，馬占虎.2005.2003年大姚和民樂－山丹地震1年尺度預(yù)測(cè)［J］.地震研究，28（3）：226-230.

Ye Q，Du X B，Chen J Y，Tan D C，Ma Z H.2005.One-year prediction for the Dayao and Minle-Shandan earthquakes in 2003［J］.Journal of Seismological Research，28（3）：226-230（in Chinese）.

張世中，石航，王蘭煒，胡哲，劉大鵬，魏連生，鞠永.2013.地電臺(tái)站受城市軌道交通干擾的測(cè)試分析與抗干擾措施研究［J］.地震學(xué)報(bào)，35（1）：117-124.

Zhang S Z，Shi H，Wang L W，Hu Z，Liu D P，Wei L S，Ju Y.2013.Test analysis on disturbances caused by urban rail transit at geoelectric stations and measures to reduce its influence［J］.Acta Seismologica Sinica，35（1）：117-124（in Chinese）.

張秀霞，李飛，楊馮威，孫春仙，肖飛.2009.蔬菜大棚對(duì)新沂地震臺(tái)地電阻率的影響［J］.地震地磁觀測(cè)與研究，30（增刊）：32-36.

Zhang X X，Li F，Yang F W，Sun C X，Xiao F.2009.About the effect of the vegetable shelter for the earth resistivity at Xinyi seismic station［J］.Seismological and Geomagnetic Observation and Research，30（Suppl.）：32-36（in Chinese）.

Balasco M，Colangelo G，Lapenna V，Loddo M，Siniscalchi A，Telesca L.2004.Measuring apparent resistivity in a seismically active area of southern Italy［J］.Phys Chem Earth，Parts A／B／C，29（4／5／6／7／8／9）：329-337.

Brace W F，Orange A S.1968.Electrical resistivity changes in saturated rocks during fracture and frictional sliding［J］.J Geophys Res，73（4）：1433-1445.

Busby J P.2000.The effectiveness of azimuthal apparent-resistivity measurements as a method for determining fracture strike orientations［J］.Geophys Prospect，48（4）：677-695.

Du X B，Xue S Z，Hao Z，Zhang S Z.2000.On the relation of moderate-short term anomaly of earth resistivity to earthquake［J］.Acta Seismological Sinica，13（4）：393-403.

Huang Q H，Lin Y F.2010.Selectivity of seismic electric signal（SES）of the 2000Izu earthquake swarm：A 3DFEM numerical simulation model［J］.Proc Japan Acad：Series B，86（3）：257-264.

Lu J，Qian F Y，Zhao Y L.1999.Sensitivity analysis of the Schlumberger monitoring array：Application to changes of resistivity prior to the 1976earthquake in Tangshan，China［J］.Tectonophysics，307（3／4）：397-405.

Lu J，Xue S Z，Qian F Y，Zhao Y L，Guan H P，Mao X J，Ruan A G，Yu S R，Xiao W J.2004.Unexpected changes in resistivity monitoring for earthquakes of the Longmen Shan in Sichuan，China，with a fixed Schlumberger sounding array［J］.Phys Earth Planet Inter，145（1／2／3／4）：87-97.

Park S K，Johnston M J S，Madden T R，Morgan F D，Morrison H F.1993.Electromagnetic precursors to earthquakes in the ULF band：A review of observations and mechanisms［J］.Rev Geophys，31（2）：117-132.

Qian F Y，Zhao Y L，Huang Y N.1996.Calculation of the parameters of georesistivity anisotropy and case history of earthquake precursors［J］.Acta Seismologica Sinica，9（4）：617-627.

Scholz C H，Sykes L R，Aggarwal Y P.1973.Earthquake prediction：A physical basis［J］.Science，181（4102）：803-810.

Wang C Y，Goodman R E，Sundaram P N，Morrison H F.1975.Electrical resistivity of granite in frictional sliding：Application to earthquake prediction［J］.Geophys Res Lett，2（12）：525-528.

Wang C Y，Sundaram P N，Goodman R E.1978.Electrical resistivity changes in rocks during frictional sliding and fracture［J］.Pure Appl Geophys，116（4／5）：717-731.

Yang C H，Cheng P H，You J I，Tsai L L.2002.Significant resistivity changes in the fault zone associated with the 1999 Chi-Chi earthquake，west-central Taiwan［J］.Tectonophysics，350（4）：299-313.