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

1）中國西安710068陜西省地震局

2）中國石家莊050021河北省地震局

引言

我國自1966年邢臺地震后開展地電阻率觀測以來，已經(jīng)形成了由近百個(gè)定點(diǎn)連續(xù)觀測臺站組成的地電阻率觀測網(wǎng)，并且在一些大震之前已記錄到顯著的前兆異常（錢復(fù)業(yè)，趙玉林，1980；錢家棟，曹愛民，1998；Lu et al，1999；杜學(xué)彬，2010）.除一部分深井電阻率觀測外，絕大多數(shù)都記錄到明顯的地電阻率年變形態(tài)，有的波動形態(tài)規(guī)律整齊，有的復(fù)雜多變.但多數(shù)都表現(xiàn)為“夏高冬低”或“夏低冬高”，年變幅度差異較大，有的可達(dá)20%—30% （錢復(fù)業(yè)等，1987），對前兆異常判斷造成一定的困難.因此，探索地電阻率觀測的正常年變形態(tài)，對于提高異常判定和前兆識別是非常必要的.

研究表明，地電阻率年變一般與地下水位、降雨和氣溫等季節(jié)性變化因素有關(guān).王志賢（1981）分析了河北隆堯縣大柏舍臺不同極距下電阻率觀測和潛水位資料，認(rèn)為大柏舍電阻率觀測與極距長短、探測深度有關(guān)；徐世浙（1985）通過理論計(jì)算分析比較了地溫和地下水位對山東省臨沂臺地電阻率觀測的影響，并指出基巖埋深對年變大小有強(qiáng)烈的控制作用；錢復(fù)業(yè)等（1987）通過理論計(jì)算、現(xiàn)場觀測及實(shí)驗(yàn)結(jié)果分析指出，地電阻率年變是探測深度偏淺的情況下表層電阻率變化引起的一種干擾，同時(shí)分析了不同電性斷面對年變形態(tài)及其幅度的影響；張學(xué)民等（1996）提出一種褶積濾波消除降雨干擾的地電阻率觀測方法；毛先進(jìn)等（2003）計(jì)算分析了表面薄層干擾對地電阻率觀測的影響，指出表層電阻率增大對地電阻率影響很小，但表層電阻率減小對地電阻率影響很明顯；Lu等（2004）通過計(jì)算影響系數(shù)發(fā)現(xiàn)地電阻率“年變”和“反年變”與介質(zhì)表層電阻率的影響系數(shù)有關(guān)；宋曉磊和黃清華（2006）采用一維非飽和滲透模型分析了降雨及其滲透過程對地電阻率觀測的影響；嚴(yán)玲琴等（2013）分析對比了地下水與降水對甘肅省臨夏臺地電阻率觀測的影響；解滔等（2013a）對江蘇新沂臺地電阻率反向年變特征的分析表明，測區(qū)地下介質(zhì)的橫向不均勻也可以引起地電阻率年變特征的變化.

大量研究表明，地電阻率觀測的年變機(jī)制復(fù)雜.由于2005年以來周至地下淺水位變化穩(wěn)定，沒有明顯年變特征.本文分析認(rèn)為，周至地電臺地電阻率年變主要由地表薄層受降雨、溫度等季節(jié)因素干擾引起.首先利用薄層干擾模型分析周至地電臺地電阻率觀測的年變特征，然后利用張學(xué)民等（1996）提出的褶積算法分析周至地電臺地電阻率NS和NE向地電阻率年變的細(xì)節(jié)特征，最后通過改變干擾薄層厚度和薄層真電阻率大小模擬降雨及其滲透過程，定性分析所得結(jié)果的合理性.

1 周至地電臺介紹

周至地電臺位于陜西省周至縣廣濟(jì)鎮(zhèn)西部，測區(qū)平坦開闊，多為農(nóng)田，無大型灌溉水渠.覆蓋層基本為黃土層，有夾砂和礫石，厚度約為200—600m.成井記錄顯示，40m以上為土層，50—80m為土層夾砂.測區(qū)位于隴縣—岐山—馬召斷裂與周至—余下斷裂交匯處，距秦嶺北緣斷裂約4km.測深反演顯示，測區(qū)地下介質(zhì)為3層K型斷面，地下15—160m存在高阻層，如圖1b所示.周至地電臺周邊地震活動較為活躍，最大地震為1988年周至ML4.3地震，近期的顯著地震為2013年戶縣ML3.0小震群活動，圖1a所示為測區(qū)內(nèi)地震、觀測臺站及斷裂的分布①引自陜西省地震局“咸陽市活斷層探測項(xiàng)目”成果..

圖1 測區(qū)內(nèi)地震、觀測臺站及斷裂分布（a）和周至地電臺垂向電測深曲線解釋（b）Fig.1 （a）Distribution of earthquakes，seismic stations and faults in the surveying area；（b）Interpretation of vertical electric soundings data for the Zhouzhi geoelectric station

周至地電臺從2007年建成開始觀測至今，使用的觀測儀器為ZD8BI地電儀.地電阻率觀測有NS、NE和NW 3個(gè)測道，測量極距為400m，供電極距為1 200m.圖2給出了周至地電臺三測道地電阻率以及流體觀測臺深井動水位、降雨及淺層地下水位等觀測數(shù)據(jù)②引自西安市水利水土保持工作總站數(shù)據(jù)..從2007年5月27日開始觀測以來，除2008年4—5月更換電極造成數(shù)據(jù)中斷外，觀測數(shù)據(jù)年變形態(tài)穩(wěn)定清晰，三測道觀測方均差基本都保持在0.02Ω·m內(nèi).自2011年底開始，三測道地電阻率觀測曲線均下降一個(gè)臺階，但年變形態(tài)穩(wěn)定.地電測區(qū)地下淺水位埋深約47m，變化平穩(wěn)但年變形態(tài)不清晰，波動基本保持在±1m的范圍內(nèi)，遺憾的是該數(shù)據(jù)只收集到2011年底.另悉周至縣縣城供水2009年前一直采用地下水水源供水，有抽水井5眼，設(shè)計(jì)日供水能力為6 000m3；2010年2月設(shè)計(jì)日供水能力為2萬m3的黑河口水廠一期主體工程已建成，2010年4月開始從黑河金盆水庫引水向周至縣城供水③2011年西安市 城區(qū)及縣城 集中 供水調(diào)度計(jì) 劃.http：∥cds.nlc.gov.cn／sxsxafz／zfgb1／20117＿1／201108／t20110818＿948619.shtml？classid＝416..分析認(rèn)為2012年至今周至地下淺水位穩(wěn)定，不會出現(xiàn)大幅度下降變化.

2 數(shù)值模擬

2.1 物理模型

圖2 周至地電臺地電阻率、降雨、深井動水位及地下淺水位觀測數(shù)據(jù)（a）—（c）NS，NE，NW向地電阻率整點(diǎn)觀測值；（d）降雨觀測值；（e）深井動水位日均值；（f）測區(qū)表層地下淺水位月觀測值① 引自西安市水利水土保持工作總站數(shù)據(jù).Fig.2 Observations of georesistivity，rainfall and undergrounded water levels at Zhouzhi geoelectric station（a）-（c）Integral point observations of NS-，NE-，NW-direction georesistivity；（d）Ten-day means of rainfall；（e）Daily means of the dynamic water level in deep well（10km far from the geoelectric survey area）；（f）Month observations of shallow water level in the geoelectric surveying area

近年來數(shù)值模擬方法已經(jīng)廣泛應(yīng)用于地震電信號解釋及影響因素分析中（黃清華，林玉峰，2010；Huang，Lin，2010；解滔等，2013b）.本文旨在分析周至地電臺地電阻率年變形態(tài)以及地下水位變化對地電阻率觀測的影響，因而未考慮自然電場及其產(chǎn)生機(jī)制.根據(jù)圖1b電測深反演結(jié)果，地電測區(qū)等效為水平3層導(dǎo)電介質(zhì)：第一層厚度為h1＝15.3m，地電阻率ρ1＝29.87Ω·m；第二層厚度為h2＝151m，地電阻率ρ2＝161Ω·m；第三層電阻率ρ3＝36.6Ω·m.A和B為供電電極，M和N 為測量電極，對稱分布于x軸上，4個(gè)電極埋深均為h，如圖3a，b所示.其中ρ11為地表薄層受干擾后的地電阻率值，干擾厚度為h′.

圖3 （a）水平3層介質(zhì)中對稱四極裝置；（b）周至地電臺地下水位升降模型Fig.3 （a）Schlumberger sounding of georesistivity observation in horizontal 3-layer medium；（b）Model for the variation of groundwater level at Zhouzhi geoelectric station

利用ANSYS電磁場耦合實(shí)體單元（SOLID5）建立水平層狀結(jié)構(gòu)的三維有限元計(jì)算模型.為了有效地減小截?cái)噙吔鐚Φ叵码妶鼍€分布的影響（湯井田，公勁喆，2010），本文經(jīng)過反復(fù)試驗(yàn)，選取長、寬、高分別為10，10，5km的模型模擬無限半空間的地下介質(zhì)（解滔等（2013b）已經(jīng)論述了尺寸選擇的敏感性，本文試驗(yàn)過程與之類似）.該模型縱向分為3層，分別表示分層介質(zhì)不同層的電性屬性.周至地電臺地電觀測系統(tǒng)為供電極距AB＝1 200m的溫納裝置（MN＝400m），電極埋深h＝2m.A，B，M 和N 的坐標(biāo)分別為（－600m，0，2m），（600m，0，2m），（－200m，0，2m），（200m，0，2m）.供電采用偶極直流電流源，在供電電極A與B兩點(diǎn)分別施加荷載電壓UA＝100V及UB＝－100V.求解穩(wěn)恒電流源給出M與N之間電勢差ΔUMN和電流I的正常理論值，并給出地電阻率觀測值為ρs＝KΔUMN／I.當(dāng)存在干擾時(shí)，觀測電壓和電流受電性結(jié)構(gòu)變化的影響而改變.設(shè)干擾后的觀測電壓和電流分別為ΔU′MN和I′，且裝置系數(shù)不變，則此時(shí)地電阻率觀測值為?ρs＝KΔU′MN／I′.定義η＝?ρs／ρs為干擾系數(shù)，則

2.2 地下水位對周至地電臺地電阻率影響的定性分析

地下水位影響地電阻率觀測的機(jī)制比較復(fù)雜，本文基于錢復(fù)業(yè)等（1987）給出的地下水位影響地電阻率觀測模型b（類）和c（類），從定性角度分析周至臺地下水位對地電阻率觀測的影響.根據(jù)錢復(fù)業(yè)等（1987）分析可知：當(dāng)?shù)叵聹\水位降低時(shí)，介質(zhì)表層低阻層厚度h1減??；當(dāng)?shù)叵聹\水位上升時(shí)，介質(zhì)表層低阻層厚度h1增加；當(dāng)?shù)叵律畈亢畬铀唤档蜁r(shí)，介質(zhì)第二層高阻層厚度h2增加；當(dāng)?shù)叵律畈亢畬铀簧邥r(shí)，介質(zhì)第二層高阻層厚度h2減小.如圖3b所示.

令參數(shù)δ＝（η－1）×100%表征水位變化影響地電阻率觀測幅度的百分比.通過改變圖3b中分層介質(zhì)第一層厚度h1和第二層厚度h2，分別模擬地下淺水位和深井水位變化對周至地電臺地電阻率觀測的影響，模擬結(jié)果如圖4a所示.當(dāng)?shù)叵聹\水位下降時(shí)，表層介質(zhì)含水量減小，真電阻率值增大，在錢復(fù)業(yè)等（1987）的模型b中表現(xiàn)為表面低阻層厚度減小，模擬結(jié)果顯示此時(shí)地電阻率觀測值降低；當(dāng)?shù)叵律罹幌陆禃r(shí)，地下第三層介質(zhì)含水量和礦物活化度降低，真電阻率值增大，在錢復(fù)業(yè)等（1987）的模型c中表現(xiàn)為地下第二層高阻層厚度增大，第一層厚度不變，模擬結(jié)果顯示此時(shí)地電阻率觀測值增加.

周至地電臺地電阻率三測道均有明顯年變，而地下淺水位則沒有明顯年變；2009年地下淺水位波動幅度較大，最大變幅達(dá)5m，而地電阻率三測道卻無明顯變化.前面介紹2010年后周至地電臺地下水開采量大幅減小，分析認(rèn)為地下淺水位不會出現(xiàn)大幅下降.模擬結(jié)果顯示，地下水位上升，地電阻率觀測值增加.而2012年開始周至地電臺地電阻率三測道均出現(xiàn)大幅度下降（不排除環(huán)境干擾因素）.因此作者認(rèn)為周至地電臺地電阻率與地下淺水位關(guān)系不明顯.距離周至地電臺東北方向10km處有深井水位觀測（井深3 200m），該深井水位在2008年汶川MS8.0地震時(shí)產(chǎn)生上升的同震變化，之后呈趨勢下降，至2013年底降幅約1.2m.周至地電臺地電阻率從2007年開始觀測以來該深井水位呈趨勢下降變化.模擬結(jié)果顯示深井水位下降，地電阻率觀測值增加.因此作者認(rèn)為周至地電臺地電阻率觀測趨勢變化與深部水位變化無關(guān).

圖4 （a）地下水位對周至地電臺地電阻率影響的模擬結(jié)果；（b）模擬給出的地電阻率年變形態(tài)Fig.4 （a）Simulation results of the effects of groundwater level on georesistivity at Zhouzhi geoelectric station；（b）Simulation results of the annual variation of georesistivity at Zhouzhi geoelectric station

2.3 地電阻率年變形態(tài)模擬

地電阻率年變一般與地下水位或表層介質(zhì)季節(jié)性變化有關(guān)，根據(jù)前面的分析本文認(rèn)為周至地電臺地電阻率年變與地下水位變化關(guān)系不大，而與地表介質(zhì)電阻率季節(jié)性變化有關(guān).假定測區(qū)表層介質(zhì)電阻率受季節(jié)影響的變化規(guī)律服從正弦函數(shù)，即夏季溫度高，表層土壤中導(dǎo)電離子活躍，電阻率?。欢颈韺与娮杪蚀?又因?yàn)橹苤恋仉娕_建臺電測深試驗(yàn)在夏季進(jìn)行，故假定表層范圍內(nèi)介質(zhì)電阻率ρ11滿足如下關(guān)系：

季節(jié)性凍結(jié)、融化對薄層的影響一般僅限于表層1—2m（錢復(fù)業(yè)等，1987），徐世浙（1985）指出地溫對地電阻率觀測的影響主要在地表5m內(nèi).本文假定季節(jié)因素的干擾厚度為5m，根據(jù)式（2）改變前面建立的物理模型中薄層h′內(nèi)的電阻率ρ11，并利用式（1）計(jì)算季節(jié)干擾系數(shù)η，如圖4b所示.

圖4b中計(jì)算給出的地電阻率年變呈“夏高冬低”形態(tài)，與周至地電臺地電阻率觀測NS向和NE向年變形態(tài)一致，較好地解釋了周至地電臺地電阻率觀測NS向和NE向年變形態(tài)，但模擬并不能給出NW向地電阻率“夏低冬高”年變形態(tài).由于周至地電臺地電測區(qū)位于隴縣—岐山—馬召斷裂與周至—余下斷裂交匯處，構(gòu)造復(fù)雜.杜學(xué)彬等（2007）建立了真、視電阻率各向異性變化與介質(zhì)裂隙率、骨架電阻率及飽水電阻率的本構(gòu)關(guān)系，分析表明介質(zhì)變形是視電阻率變化及其各向異性變化的根本原因（杜學(xué)彬，2010）.NW向地電阻率“夏低冬高”年變形態(tài)可能與構(gòu)造作用下裂隙排布及裂隙水變化有關(guān).根據(jù)解滔等（2013a）分析新沂臺反向年變特征認(rèn)為，上述地電阻率年變形態(tài)也可能與周至臺地電測區(qū)的復(fù)雜構(gòu)造有關(guān).但目前缺少該測區(qū)地質(zhì)剖面圖和沿NW向測線的電測深剖面資料，關(guān)于NW向地電阻率“夏低冬高”這一年變形態(tài)還需深入調(diào)查研究.

3 褶積計(jì)算分析降雨對周至地電臺地電阻率的影響

圖2a，b顯示，2011年前周至地電臺地電阻率NS向與NE向觀測值每年夏季的“峰值”基本保持一致，而冬季的“低值”卻不盡一致，從2008年開始地電阻率年變幅度逐年增大.同時(shí)距離地電測區(qū)10km處地下流體觀測點(diǎn)觀測的降雨數(shù)據(jù)顯示，2011年前周至地電臺年均降雨量逐年增大（圖2d）.為此本文利用張學(xué)民等（1996）給出的降雨褶積分析方法分析降雨對周至臺地電阻率觀測的可能影響.褶積算法為

式中，Δρ（t）為降雨對地電阻率的影響量，CR（t）為降雨量時(shí)間序列，R（t）為褶積算法的系統(tǒng)函數(shù)，

聯(lián)合式（3）與（4）可得

根據(jù)張學(xué)民等（1996）的結(jié)果，選擇瞬時(shí)影響時(shí)間M為1天，滯后影響時(shí)間N為6個(gè)月.選取2009年5月8—31日、2009年6月7—18日、2009年8月1—15日、2010年5月12—31日、2010年7月1—31日和2010年8月13—31日等6個(gè)時(shí)段共計(jì)116天的降雨數(shù)據(jù)和NS向地電阻率整點(diǎn)觀測數(shù)據(jù)，分別給出系數(shù)矩陣K和降雨影響地電阻率變化量Δρ（t），并考慮t＝M時(shí)分段函數(shù)R（t）的連續(xù)性，求解超定方程組（式（5））得參數(shù)矩陣X為X＝ ［1.833150e－3，－2.800274e－5，－3.303554e－4，5.857361e－7，4.826080e－2，－2.673098e－1］T.

將參數(shù)X代入式（5），并將其作用于周至流體觀測臺歷年降雨觀測數(shù)據(jù)整點(diǎn)值，給出周至地電臺降雨對地電阻率觀測影響量Δρ（t）的短期變化和長趨勢變化曲線，如圖5所示.由于NS向觀測曲線與NE向形態(tài)基本一致，圖中只給出NS向觀測曲線.

圖5 周至地電臺降雨對地電阻率觀測（NS向）的短期影響（a）及長趨勢影響（b）Fig.5 Short-term （a）and long-term impacts（b）of the rainfall on georesistivity in NS-direction at Zhouzhi georesistivity station

通過分析上述褶積模型計(jì)算的周至地電臺降雨對地電阻率觀測的影響可見，降雨短期內(nèi)會造成地電阻率觀測值增大，且褶積擬合結(jié)果與實(shí)際觀測結(jié)果基本一致，如圖5a所示.但對于長趨勢，褶積擬合結(jié)果顯示為降雨對地電阻率呈現(xiàn)觀測“夏低冬高”的年變形態(tài)，且夏季“低值”基本保持恒定，而冬季“峰值”卻呈現(xiàn)出逐年增大的變化趨勢，如圖5b所示.這與2007—2011年底周至地電臺地電阻率觀測NS向與NE向?qū)嶋H觀測值呈負(fù)相關(guān)關(guān)系，也與圖5a所示短期相關(guān)性相反，如圖5b所示.如何解釋“短期影響呈正相關(guān)”，而“長趨勢影響呈負(fù)相關(guān)”這一現(xiàn)象？周至地電臺地電阻率觀測NS向和NE向年變幅度增大是否由降雨入滲引起？由于缺乏水文地質(zhì)資料以及土壤滲流相關(guān)的物理參數(shù)，為了解釋上述問題，本文僅在前面建立的物理模型基礎(chǔ)上改變薄層的厚度h′和薄層電阻率ρ11，定性模擬降雨入滲過程對地電阻率觀測的影響.假定降雨后表面薄層電阻率減小，以ρ11＝0.5，2，5，25Ω·m定性表示薄層土壤含水后電阻率的變化.通過增加薄層厚度h′定性模擬降水滲透過程對地電阻率觀測的影響，模擬結(jié)果如圖6所示.可以看出，隨著深度h′的增加，地電阻率呈先增大后轉(zhuǎn)折為減小的趨勢.表層電阻率越小，轉(zhuǎn)折出現(xiàn)得越早.

降雨之后，短時(shí)間內(nèi)地表受干擾的薄層厚度很小，但薄層介質(zhì)電阻率卻大幅度降低；隨著時(shí)間推移，雨水滲透深度逐漸增大，干擾薄層厚度也逐漸增大，并且隨著含水率的減小，薄層內(nèi)介質(zhì)電阻率也慢慢增大.模擬結(jié)果表明，在降雨后的短時(shí)間內(nèi)，地電阻率隨干擾薄層厚度的增大而增加，達(dá)到一個(gè)臨界厚度后，地電阻率會隨著干擾薄層厚度的增大而減小，如圖6所示.該定性模擬結(jié)果與褶積擬合分析結(jié)果相吻合，因此本文認(rèn)為周至地電臺地電阻率NS向和NE向年變幅度逐年變化與當(dāng)?shù)啬杲涤炅康牟煌嘘P(guān).

圖6 周至地電臺降雨入滲對地電阻率觀測影響的定性分析Fig.6 Qualitative analysis on the effect of rainfall infiltration on georesistivity at Zhouzhi geoelectric station

2011年底至今周至地電臺地電阻率NS向和NE向夏季“峰值”較2012年前出現(xiàn)了一個(gè)平移式的下降，并且在其曲線形態(tài)上與褶積擬合給出的長趨勢影響量基本平行（圖5b）.這表明基于褶積計(jì)算給出的曲線動態(tài)特征在2011年后與實(shí)際觀測結(jié)果仍然一致.而上述降雨影響分析及地下水位變化分析并不能對二者之間的平移式下降臺階給出合理的解釋.本文認(rèn)為2011年底—2012年3月周至地電臺NS和NE測道的測值下降可能來自測區(qū)某種固定干擾源，該干擾從出現(xiàn)至穩(wěn)定的過程與NS和NE測道測值逐漸下降直至出現(xiàn)穩(wěn)定年變形態(tài)相對應(yīng).有關(guān)這一問題還需進(jìn)一步的調(diào)查落實(shí)和分析.

4 討論與結(jié)論

本文基于測深反演結(jié)果，通過數(shù)值模擬定性分析了地下水位和降雨對周至地電臺地電阻率觀測的影響，并利用褶積算法定量分析了降雨對周至地電臺NS向和NE向地電阻率觀測的短期影響和長趨勢影響.得到以下結(jié)論：

1）周至地電臺地電阻率NS和NE向“夏高冬低”年變形態(tài)主要由地表薄層介質(zhì)電阻率季節(jié)變化引起，受地下水位變化影響不明顯.

2）周至地電臺地電阻率NS向和NE向夏季“高值”穩(wěn)定、冬季“低值”逐年變化與同期降雨量逐年變化有關(guān).這可能與降雨引起的表層介質(zhì)電阻率變化有關(guān).盡管每年夏季降雨量和降雨持續(xù)時(shí)間不同，但降雨對地電阻率的影響需要半年左右的時(shí)間（張學(xué)民等，1996；宋曉磊，黃清華，2006）才能穩(wěn)定；此外夏季溫度較高，降水多，表層介質(zhì)內(nèi)離子活躍度穩(wěn)定，介質(zhì)真電阻率不易出現(xiàn)較大差異變化，而冬季降水少，一年中降水的時(shí)間滯后效應(yīng)開始逐漸顯現(xiàn)，由于年降雨量不同，其對地電阻率影響的滯后效應(yīng)累加也不盡一致.

3）周至地電臺地電阻率NS和NE測道從2011年底—2012年3月出現(xiàn)大幅下降，之后轉(zhuǎn)為出現(xiàn)年變形態(tài)，且對降雨的敏感性較往年更為強(qiáng)烈.同時(shí)，褶積擬合結(jié)果給出的地電阻率年變曲線在2012年3月后與實(shí)際觀測曲線基本平行.這表明2011年底—2012年3月NS和NE測道的測值下降可能來自測區(qū)某種固定干擾源，該干擾從出現(xiàn)至穩(wěn)定的過程與NS和NE測道測值逐漸下降直至出現(xiàn)穩(wěn)定年變形態(tài)相對應(yīng).

4）由于缺乏測區(qū)地質(zhì)剖面資料，本文采取的計(jì)算模型簡單.模擬結(jié)果沒有給出NW向地電阻率“夏低冬高”的年變形態(tài).測區(qū)附近NW走向的隴縣—岐山—馬召斷裂與測區(qū)NW向測線近似平行.與之巧合的是解滔等（2013a）給出的新沂臺測區(qū)也存在與NS測線近似平行的大賀山—橋北鎮(zhèn)斷裂以及斷裂f5.二者之間是否存在必然的聯(lián)系，還需要進(jìn)一步分析.另外，周至地電臺位于岐山—馬召斷裂、周至—余下斷裂及秦嶺北緣斷裂附近，其構(gòu)造復(fù)雜，因此NW向地電阻率反向年變也有可能與構(gòu)造作用下裂隙排布及裂隙水變化有關(guān).

周至地電臺彭玉柱臺長、張世民高工在地電阻率異常落實(shí)中給予協(xié)助，陜西省地震預(yù)報(bào)中心趙小茂高工在獲取地下淺水位資料中給予幫助，監(jiān)測中心李守廣高工提供電測深數(shù)據(jù)，與中國地震臺網(wǎng)中心解滔同志就本文進(jìn)行了有益討論和交流，以及評審老師對本文提出了寶貴修改意見.作者在此一并表示感謝！

杜學(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é)彬.2010.在地震預(yù)報(bào)中的兩類視電阻率變化［J］.中國科學(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.

黃清華，林玉峰.2010.地震電信號選擇性數(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 influence［J］.Chinese Journal of Geophysics，53（3）：535-543（in Chinese）.

毛先進(jìn)，李志祥，楊玲英.2003.臺址下表層介質(zhì)電阻率變化對視電阻率觀測值影響特征的初步研究［J］.地震研究，26（4）：338-342.

Mao X J，Li Z X，Yang L Y.2003.Primary research of interference on observed apparent resistivity value due to resistivity variation in surface layer under seismic station［J］.Journal of Seismological Research，26（4）：338-342（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）.

錢復(fù)業(yè)，趙玉林，許同春.1987.地電阻率季節(jié)干擾變化分析［J］.地震學(xué)報(bào)，9（3）：289-302.

Qian F Y，Zhao Y L，Xu T C.1987.Analyses of the seasonal variation of disturbance in georesistivity［J］.Acta Seismologica Sinica，9（3）：289-302（in Chinese）.

錢家棟，曹愛民.1998.1976年唐山7.8級地震地電阻率和地下水前兆綜合物理機(jī)制研究［J］.地震，18（增刊）：1-9.

Qian J D，Cao A M.1998.Research on comprehensive mechanism of precursors in apparent resistivity and water table associated with 1976Tangshan earthquake（M＝7.8）［J］.Earthquake，18（Suppl.）：1-9（in Chinese）.

宋曉磊，黃清華.2006.降水滲透過程分析及其對視電阻率的影響［J］.北京大學(xué)學(xué)報(bào)：自然科學(xué)版，42（4）：470-477.

Song X L，Huang Q H.2006.Rainfall infiltration process and its effect on apparent resistivity［J］.Acta Scientiarum Naturalium Universitatis Pekinensis，42（4）：470-477（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）.

王志賢.1981.地電阻率年變化的一種物理解釋：大柏舍臺不同極距地電阻率觀測結(jié)果［J］.西北地震學(xué)報(bào)，3（3）：43-47.

Wang Z X.1981.A physical explanation of annual variation of earth resistivity：Observations at Daboshe station of earth resistivity［J］.Northwestern Seismological Journal，3（3）：43-47（in Chinese）.

解滔，李飛，沈紅會，盧軍.2013a.新沂地震臺地電阻率反向年變分析［J］.地震學(xué)報(bào)，35（6）：856-864.

Xie T，Li F，Shen H H，Lu J.2013a.Reverse annual variation of georesistivity at Xinyi seismic station，Jiangsu Province［J］.Acta Seismologica Sinica，35（6）：856-864（in Chinese）.

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

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

徐世浙.1985.視電阻率年變的定量計(jì)算［J］.地震學(xué)報(bào)，7（4）：422-427.

Xu S Z.1985.Quantitative calculation of the annual variation of apparent resistivity［J］.Acta Seismologica Sinica，7（4）：422-427（in Chinese）.

嚴(yán)玲琴，鄭衛(wèi)平，張輝，張昱.2013.臨夏臺地電阻率變化與震兆現(xiàn)象分析［J］.中國地震，29（1）：168-176.

Yan L Q，Zheng W P，Zhang H，Zhang Y.2013.Analysis on the variation of the earth resistivity and seismic precursory phenomena at Linxia seismic station，Gansu［J］.Earthquake Research in China，29（1）：168-176（in Chinese）.

張學(xué)民，王志賢，臧明珍，鄭慰坤.1996.降雨對地電阻率干擾的分析［J］.華北地震科學(xué)，14（4）：71-75.

Zhang X M，Wang Z X，Zang M Z，Zheng W K.1996.Analysis of interference of rainfall on earth resistivity［J］.North China Earthquake Sciences，14（4）：71-75（in Chinese）.

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）：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］.Physics Earth Planet Inter，145（1／2／3／4）：87-97.