毛先進(jìn) 段煒 莊儒新 楊玲英 趙晉民

摘要：為壓制地電阻率淺層干擾、突出深部以地震預(yù)報(bào)為目的的有用信息，選取了小江斷裂中段一個(gè)場(chǎng)點(diǎn)作為實(shí)例，研究了井下地電阻率前兆觀測(cè)中確定電極埋深、電極間距等布極參數(shù)的方法。結(jié)果表明：選擇井下電阻率觀測(cè)布極參數(shù)時(shí)要考慮影響系數(shù)和探測(cè)深度2個(gè)主要因素，即地下潛水位面以上各地層的影響系數(shù)應(yīng)遠(yuǎn)小于深部各層，探測(cè)深度范圍內(nèi)最底層（受孕震影響最大的層）的影響系數(shù)應(yīng)遠(yuǎn)大于其它各層，觀測(cè)系統(tǒng)的探測(cè)深度最好不小于已有震例指示的地電阻率觀測(cè)系統(tǒng)的探測(cè)深度。按照該方法選擇了所給場(chǎng)點(diǎn)的井下電阻率觀測(cè)布極參數(shù)，即電極埋深200m、供電極距1050m、測(cè)量極距350m，按對(duì)稱四極布置，可獲得最佳觀測(cè)效果。

關(guān)鍵詞：地電阻率觀測(cè);布極參數(shù)確定方法;影響系數(shù);探測(cè)深度

中圖分類號(hào)：P315.3文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2019）01-0096-06

0引言

在國(guó)內(nèi)50余年的地電阻率前兆觀測(cè)中，1976年唐山7.8級(jí)、1988年瀾滄—耿馬7.6和7.2級(jí)、2008年汶川8.0級(jí)、2013年蘆山7.0級(jí)等多次大地震以及一些中強(qiáng)地震之前，震中區(qū)周圍300km范圍內(nèi)的地電阻率資料都出現(xiàn)過中期異常變化，這表明地電阻率是一種比較可靠的地震前兆（趙玉林，錢復(fù)業(yè)，1978;錢家棟等，1985，2013;桂燮泰等，1989;葉青等，2005;張學(xué)明等，2009;Luetal，1999;錢復(fù)業(yè)等，1980，1987;顧功敘，1981;張國(guó)明等，2001;杜學(xué)彬，2010）。但由于地下淺層介質(zhì)電阻率受季節(jié)變化的影響，數(shù)據(jù)變化具有多源性，實(shí)際工作中難以確定數(shù)據(jù)變化的真實(shí)原因，這可能是這些大震發(fā)生前未能做出短期或臨震預(yù)報(bào)的原因之一。因此尋求壓制與地震孕育無關(guān)的淺層變化（本文稱之為干擾）的觀測(cè)方法成為必然，井下地電阻率觀測(cè)是較為有效的方法之一。

近年來，我國(guó)陸續(xù)開展了井下地電阻率觀測(cè)的理論與試驗(yàn)研究（王蘭煒等，2015），已有研究結(jié)果表明井下觀測(cè)優(yōu)于傳統(tǒng)的地表觀測(cè)（張磊等，2015;髙曙德，2016），是一種有前景的觀測(cè)方法。目前國(guó)內(nèi)開展了井下地電阻率觀測(cè)的臺(tái)站中，不同臺(tái)站的電極埋深和電極間距（簡(jiǎn)稱為布極參數(shù)，用C表示，下同）有很大差別，觀測(cè)結(jié)果表明各臺(tái)站壓制淺層干擾的能力參差不齊。因此，如何確定布極參數(shù)有效壓制淺層干擾，從而有效獲取深部變化信息，是一個(gè)十分重要的問題。本文以小江斷裂中段附近一個(gè)場(chǎng)點(diǎn)為例，以有效壓制淺層干擾變化、獲取深部有用變化信息為目的，研究井下地電阻率觀測(cè)時(shí)最優(yōu)布極參數(shù)的確定方法。

1最優(yōu)布極參數(shù)確定原則

1.1影響系數(shù)

對(duì)于一維情形，假定地電阻率觀測(cè)區(qū)是水平層狀介質(zhì)，共n層，第i層電阻率為ρi，層厚為hi（i=1，2，…，n），且假定底層的層厚無限大，即hi→∞。在地表或地下四極對(duì)稱裝置觀測(cè)到的視電阻率ρs既是地下各層電阻率和層厚的函數(shù)，同時(shí)也是布極參數(shù)的函數(shù)，即：

式中：C代表布極參數(shù)。

假定各層只是地電阻率發(fā)生變化，層厚保持不變。由式（1）可以得到：

式中：Si=ρiρsfρi，稱為影響系數(shù)，它是在給定的布極參數(shù)下，第i層地電阻率變化率對(duì)觀測(cè)值變化率的貢獻(xiàn)權(quán)系數(shù)，由介質(zhì)電性參數(shù)和布極參數(shù)共同決定。

Seigel（1959）已證明影響系數(shù)之和等于1，某些層的影響系數(shù)之和增大，其余層的影響系數(shù)之和就會(huì)減小。由式（2）可見，通過控制影響系數(shù)在各層的分配就能壓制或突出某些層變化信息。近年來，關(guān)于水平層狀介質(zhì)、點(diǎn)源在地下或地表時(shí)各層介質(zhì)影響系數(shù)的計(jì)算方法已開展了逐步深入的研究（聶永安等，2009，2010;毛先進(jìn)等，2014;解滔等，2016），為影響系數(shù)的計(jì)算打下了良好的基礎(chǔ)。

由于地下潛水位受大氣降雨、灌溉、抽水等因素影響而隨時(shí)間波動(dòng)，因此潛水位面以上各層介質(zhì)含水量隨之波動(dòng)，這些地層的地電阻率也發(fā)生波動(dòng)，使觀測(cè)到的視電阻率發(fā)生變化，因這種變化與地震孕育無關(guān)而稱之為干擾;數(shù)十年的地電阻率觀測(cè)資料表明這種波動(dòng)大致具有一年周期，但波形并不規(guī)則（杜學(xué)彬等，2017），僅用數(shù)學(xué)手段并不能很好地消除其影響;從地電阻率觀測(cè)目的看，這些層位的影響系數(shù)越小越好。有關(guān)研究（毛先進(jìn)等，2014;解滔等，2016）表明，隨著電極埋深（布極參數(shù)之一）的變化，各層影響系數(shù)并非單調(diào)變化，而是呈現(xiàn)起伏變化。在一定的電極間距下，電極埋設(shè)在一定深度范圍內(nèi)時(shí)對(duì)淺層干擾具有放大作用，換言之，井下地電阻率觀測(cè)時(shí)布極參數(shù)須有所選擇。筆者認(rèn)為，要達(dá)到有效壓制淺層干擾、有效獲取深部有用信息的目的，應(yīng)當(dāng)使“深部地層”的影響系數(shù)明顯大于易受干擾的淺部地層的影響系數(shù)。

1.2探測(cè)深度

直流電阻率觀測(cè)的依據(jù)是電流的體積效應(yīng)，在井下觀測(cè)時(shí)，前述“深部地層”實(shí)際上是電流分布區(qū)域中靠下的部分，這一部分的深度多少合適，應(yīng)該由其電阻率的變化能夠反映地震的孕育和發(fā)展過程并且能被觀測(cè)到來確定，因此涉及探測(cè)深度問題。

關(guān)于點(diǎn)源在地下或地表時(shí)的探測(cè)深度，我們采用趙和云和錢家棟（1982）給出的探測(cè)深度的定義，利用二維邊界積分方程法對(duì)幾種典型電阻率結(jié)構(gòu)做初步研究，得到了關(guān)于探測(cè)深度的一些有益的認(rèn)識(shí)（毛先進(jìn)等，2017）。本文研究針對(duì)的是一維電性結(jié)構(gòu)，因此采用相同的探測(cè)深度的定義，但點(diǎn)源電場(chǎng)的計(jì)算采用一維方法（葛為中，1994）。

我國(guó)（地表）直流地電阻率前兆觀測(cè)的供電極距一般為1000m，50余年的觀測(cè)經(jīng)驗(yàn)表明，這樣的觀測(cè)裝置其探測(cè)深度范圍內(nèi)電阻率的變化是能夠反映大地震的孕育和發(fā)展過程，且能被觀測(cè)到的。對(duì)于這樣的裝置和較為常見的電阻率水平分層均勻結(jié)構(gòu)，下伏低阻與下伏高阻結(jié)構(gòu)探測(cè)深度分別約為550m與300m（毛先進(jìn)等，2017），井下地電阻率觀測(cè)時(shí)最好不小于這樣的深度。同時(shí)，在電極埋深方面還應(yīng)考慮建設(shè)成本和技術(shù)可行性。

2布極參數(shù)確定實(shí)例

2.1場(chǎng)點(diǎn)選擇

本文研究的場(chǎng)點(diǎn)位于小江斷裂帶中段西支附近（圖1）。小江斷裂帶北起巧家以北，南至建水東南，全長(zhǎng)超過400km，其形成時(shí)代較早、發(fā)育歷史長(zhǎng)，經(jīng)歷過多期構(gòu)造活動(dòng)。新構(gòu)造時(shí)期以來，該斷裂活動(dòng)表現(xiàn)為強(qiáng)烈的左旋走滑運(yùn)動(dòng)以及兩側(cè)斷塊垂直差異活動(dòng)，沿?cái)嗔训卣鸹顒?dòng)強(qiáng)烈。這些活動(dòng)表明斷裂附近有利于構(gòu)造應(yīng)力積累，而應(yīng)力積累過程中地下介質(zhì)的變形會(huì)導(dǎo)致其電阻率發(fā)生變化，因此選擇在小江斷裂帶中段開展壓制淺層非震變化、有效獲取深部構(gòu)造運(yùn)動(dòng)變化信息的井下地電阻率觀測(cè)方法研究是合適的。

2.2布極參數(shù)確定

通過電測(cè)深獲得場(chǎng)點(diǎn)的一維電阻率結(jié)構(gòu)，見表1。該場(chǎng)點(diǎn)按電阻率分為6層，為下伏低阻結(jié)構(gòu)，電測(cè)深測(cè)線方向?yàn)槟媳毕颉?/p>

該場(chǎng)點(diǎn)位于嵩明盆地，1～3層介質(zhì)厚度之和為28.1m。根據(jù)參考場(chǎng)點(diǎn)附近多個(gè)巖土工程勘察報(bào)告，該場(chǎng)點(diǎn)地下潛水位隨季節(jié)變化，雨季0.2m，旱季不超過10m，因此干擾（潛水位變化對(duì)介質(zhì)電阻率的影響）主要發(fā)生在1～3層，4～6層則不受影響。按照前述原則，確定布極參數(shù)時(shí)應(yīng)使1～3層的影響系數(shù)之和盡量小，而4～6層影響系數(shù)之和則盡量大。

由于影響系數(shù)隨電極埋深的變化很復(fù)雜（毛先進(jìn)等，2014;解滔等，2016），本文僅針對(duì)具體的電極埋深進(jìn)行討論。

表2是電極埋深2m（我國(guó)大多數(shù)地面地電阻率觀測(cè)臺(tái)站的電極埋深）、不同電極間距條件下各層介質(zhì)的影響系數(shù)和探測(cè)深度計(jì)算結(jié)果，其中AB，MN分別為供電極距與測(cè)量極距，Si（i=1，2，…，6）為1～6層介質(zhì)影響系數(shù)。從表2可見，在8種大小不同的電極距中，1～3層介質(zhì)（干擾層）的影響系數(shù)隨AB增大而減小;干擾層之下的4～6層介質(zhì)的影響系數(shù)隨AB增大而增大，同時(shí)探測(cè)深度也隨AB增大而增大。當(dāng)AB=150m時(shí)，1～3層介質(zhì)的影響系數(shù)之和為0.9031，4～6層的為0.0967，前者是后者的9.3倍，觀測(cè)數(shù)據(jù)的變化主要反應(yīng)1～3層介質(zhì)的干擾;當(dāng)AB=1200m時(shí)，1～3層介質(zhì)的影響系數(shù)之和為0.0704，4～6層的為0.9293，后者是前者的13.2倍，對(duì)壓制1～3層的干擾變化很有利，而且探測(cè)深度達(dá)到579.6m。但由于4～6層的影響系數(shù)分別為0.5727，0.2424，0.1142，觀測(cè)系統(tǒng)對(duì)底界面埋深只有66.1m的第4層最敏感，而對(duì)其探測(cè)深度范圍內(nèi)受孕震影響最大的第6層則不夠敏感;另外工業(yè)電磁干擾對(duì)地表觀測(cè)的影響是最大的。綜合來看，對(duì)于該場(chǎng)點(diǎn)，即使供電極距達(dá)到1200m，地表淺埋電極的觀測(cè)方式也難以達(dá)到地震前兆觀測(cè)的目的。

表3是電極埋深150m、不同電極間距條件下各層介質(zhì)的影響系數(shù)和探測(cè)深度計(jì)算結(jié)果。從表中可見，在8種大小不同的電極距中，第1～3層介質(zhì)（干擾層）的影響系數(shù)之和隨AB增大而增大，當(dāng)AB=1200m時(shí)為0.0881，第4～6層介質(zhì)的影響系數(shù)之和最小，為0.9117，后者是前者的10.3倍，對(duì)1～3層的干擾變化有較好的壓制能力;當(dāng)AB=150m時(shí)，第1～3層介質(zhì)的影響系數(shù)之和僅為0.002，第4～6層的為0.998，后者是前者的499.0倍，對(duì)壓制1～3層的干擾變化十分有利。

僅從影響系數(shù)看，表3中所有電極間距方案都可以考慮，然而不同方案的探測(cè)深度差異很大，當(dāng)AB=150m時(shí)探測(cè)深度只有194.9m，遠(yuǎn)達(dá)不到現(xiàn)有震例所指示的探測(cè)深度要求;只有最后2種方案達(dá)到探測(cè)深度不小于550m的要求。

表4是電極埋深200m、不同電極間距條件下各層介質(zhì)的影響系數(shù)和探測(cè)深度計(jì)算結(jié)果。與表2類似，8種大小不同的電極距下1～3層介質(zhì)的影響系數(shù)之和隨AB增大而增大，當(dāng)AB=1200m時(shí)為0.0597，第4～6層介質(zhì)的影響系數(shù)之和最大，為0.9401，后者是前者的15.7倍，對(duì)1～3層的干擾變化有良好的壓制能力;當(dāng)AB=150m時(shí)，第1～3層介質(zhì)的影響系數(shù)之和僅為0.0006，第4～6層的為0.9993，后者是前者的1665.5倍，對(duì)壓制來自于1～3層的干擾變化非常有利。

從影響系數(shù)看，表4中所有電極距方案都符合要求，然而不同方案的探測(cè)深度差異很大，當(dāng)AB=150m時(shí)，第5層介質(zhì)的影響系數(shù)達(dá)到0.8528，采用這種布極方案，其觀測(cè)數(shù)據(jù)主要反映第5層的變化，相應(yīng)的探測(cè)深度只有240.6m，盡管電極埋深200m，依然達(dá)不到已有震例指示的探測(cè)深度要求。只有最后3種方案符合探測(cè)深度不小于550m的要求。

綜合表3，4來看，從達(dá)到有效壓制淺層非震干擾、獲取深部有用信息的觀測(cè)目的，建設(shè)成本以及技術(shù)可行性等方面綜合考慮，該場(chǎng)點(diǎn)可選擇埋深200m、供電極距1050m、測(cè)量極距350m的對(duì)稱四極布置方案（表4），此時(shí)電阻率受地下水位影響最大的1～3層介質(zhì)影響系數(shù)之和為0.0566，4～6層影響系數(shù)之和為0.9432，后者為前者的16.7倍，有良好的壓制淺層干擾的能力，其探測(cè)深度為619.8m，完全滿足現(xiàn)有震例指示的地電阻率前兆觀測(cè)系統(tǒng)的探測(cè)深度要求，而且受孕震影響最大的第6層的影響系數(shù)達(dá)到0.7862，觀測(cè)數(shù)據(jù)的變化將能包含更多的有用信息。

3結(jié)論

壓制淺層干擾、有效獲取深部構(gòu)造運(yùn)動(dòng)信息是包括地電阻率的所有地震前兆觀測(cè)的目的。本文以小江斷裂北段嵩明盆地內(nèi)一個(gè)場(chǎng)點(diǎn)為例，計(jì)算并分析了不同布極參數(shù)條件下各層介質(zhì)的影響系數(shù)及觀測(cè)系統(tǒng)探測(cè)深度變化情況，對(duì)采用不同布極參數(shù)時(shí)壓制淺層干擾、獲取深部介質(zhì)電阻率變化信息的能力進(jìn)行了分析，研究表明：

（1）對(duì)于給定的一維電阻率結(jié)構(gòu)和不同布極參數(shù)的觀測(cè)系統(tǒng)，可以首先用影響系數(shù)和探測(cè)深度這2個(gè)指標(biāo)預(yù)判觀測(cè)效果，然后通過選擇合適的布極參數(shù)達(dá)到最佳觀測(cè)效果。

（2）分析表明，在確定井下電阻率地震前兆觀測(cè)系統(tǒng)的布極參數(shù)時(shí)，首先應(yīng)考慮影響系數(shù)分布，即地下潛水位面以上各地層的影響系數(shù)之和應(yīng)遠(yuǎn)小于深部各層影響系數(shù)之和，同時(shí)還應(yīng)使最下層（受孕震影響最大的層）的影響系數(shù)遠(yuǎn)大于其它各層的影響系數(shù)，其次要考慮觀測(cè)系統(tǒng)的探測(cè)深度，最好不小于過去50余年已有震例指示的探測(cè)深度要求。