佟　鑫，郭建芳，周劍青，張思遠

(河北省地震局秦皇島中心臺，河北 北戴河　066100)

?

昌黎地電臺多極距觀測系統(tǒng)的設計

佟鑫，郭建芳*，周劍青，張思遠

(河北省地震局秦皇島中心臺，河北 北戴河066100)

摘要：基于對昌黎地電臺地下介質(zhì)的一維電性結(jié)構(gòu)模型，分析了不同深層介質(zhì)的響應系數(shù)、多極距布設方式及理論探測深度。結(jié)果表明：當MN分別取1/3 AB、1/4 AB、1/5 AB時，各層的響應系數(shù)變化趨勢非常一致；據(jù)初始電性結(jié)構(gòu)模型分析，4個不同深層中第3層的電阻率測值貢獻率最高，達到70%，因此，在設計多極距觀測布極時主要考慮第3層的貢獻，該層的供電極距對取300～500 m較合適，第3層最大貢獻率為60%，當供電極距對不超過1 000 m時，第4層貢獻率僅為20%；根據(jù)不同深層介質(zhì)響應系數(shù)分布情況和臺站實際觀測條件限制，多極距布設方式為5種，分別為1 000 m、500(或400)m、400(或300)m、125 m、30 m。5種布極的理論探測深度在14～96 m，此結(jié)果由昌黎臺介質(zhì)的電性結(jié)構(gòu)所決定。

關(guān)鍵詞：電性結(jié)構(gòu)模型；響應系數(shù)；多極距布極探測深度；昌黎臺

0引言

昌黎臺在1975年海城MS7.3地震、1976年唐山MS7.8大震等觀測到了突出的地電阻率異常[1]，唐山老震區(qū)及周邊自1976年大地震以來未發(fā)生過6級以上地震。華北東部及東北地區(qū)自2010年以來發(fā)生了數(shù)次MS4～5級地震，且震級逐年增大，其中以2012年5月唐山MS4.8，2013年5月、10月吉林5級震群活動最為顯著。而備受關(guān)注的唐山老震區(qū)在我國西、南部中強地震活躍、毗鄰各國強震多發(fā)的時期，維持了近3年的4級地震平靜、3級地震偏少的狀態(tài)。近期研究也證明：唐山斷裂處于孕震時期，其中灤縣—樂亭斷裂在NWW向主壓應變場中[2]，值得關(guān)注。

昌黎臺地電阻率多年來一直作為華北片區(qū)、首都圈及唐山地區(qū)中強震監(jiān)測的重點臺站，是河北省地震局和中國地震局分析預報工作的重點監(jiān)測手段之一，但主要以地電阻率時間序列上提取趨勢異常為重點。因多數(shù)前兆觀測受干擾影響，無法顯現(xiàn)真實的年變動態(tài)及異常特征，故而多數(shù)人采用不同數(shù)學方法提取異常。鑒于此，該臺于2015年開始實施地電阻率多極距觀測實驗，反演計算不同深度地電阻率動態(tài)變化的方法，研究地下分層介質(zhì)電阻率變化對地電阻率地震前兆信息的貢獻，并結(jié)合地下水觀測研究地電阻率變化的機理。本文重點對地電阻率多極距觀測系統(tǒng)的設計問題做分析和整理。

早期我國地電阻率觀測臺站的布極方式多為單測道(或正交雙測道)單極距觀測，采用對稱四極裝置，得到的視電阻率變化是觀測范圍內(nèi)位于不同深度層電阻率的綜合反映。年變特征主要受淺層水位及季節(jié)變化導致淺部地電阻率發(fā)生顯著變化的影響[3-4]，而深部地電阻率變化則是我們需要的與映震有關(guān)的信息。為了得到有價值的地震信息，早期有些臺站采用多極距觀測的方法獲得等效(真)電阻率[5-6]，采用深埋電極裝置進行實驗研究[7-8]：聶永安等通過計算水平層狀介質(zhì)電源激勵解析解的方式，從理論上研究了3層結(jié)構(gòu)在地表和基巖中對稱四級電阻率觀測結(jié)果與裝置深度、供電電極之間的關(guān)系[9-10]；解滔等對井下電阻率觀測中地表電流的干擾進行了計算分析，為如何在壓制地表電流干擾影響的同時選擇臺址電性結(jié)構(gòu)和電極埋深、確定干擾源避讓距離等方面提供參考[11]。另外，電阻率層析成像技術(shù)也在相關(guān)領(lǐng)域獲得新的研究[12]。

對于多極距觀測系統(tǒng)而言，了解不同層位的響應系數(shù)很關(guān)鍵，其能反映各深層的視電阻率對觀測值的貢獻大小。本文重點對昌黎臺的地下電性結(jié)構(gòu)層做響應系數(shù)分析，從而擬定多極距觀測的合理布設方式，同時對探測深度也做淺析。

1臺站概述

昌黎地電臺始建于1972年，該臺位于河北省東部沿海地區(qū)，NE走向的寧河—昌黎深斷裂位于測區(qū)東約5 km處，西側(cè)約20 km和65 km處為NW走向的灤縣—樂亭斷裂和NE走向的唐山斷裂。第四系覆蓋層厚70 m ，含水層發(fā)育，總厚度20～30 m，基巖風化厚度8 m，主要開采20 m 和30 m兩含水層，含水層巖性為粗砂和中砂。下伏基巖為燕山期花崗巖，電性條件屬2層電性斷面。上覆第四系為低阻層，ρ2大體在幾到幾百歐姆·米，下伏基巖為高阻層，ρ2大體在幾千至數(shù)萬歐姆·米。據(jù)歷史鉆孔資料及層析二維圖像，2個高阻層在10～40 m和130 m深度以下，分別對應于含水量甚少的膠泥層和花崗巖層；在40～130 m深度段的低阻層則是富含地下水的中粗砂層[13]。臺站地勢平坦，西側(cè)1.5 km處有深部承壓熱礦水，水位和水溫高，礦化度大，屬CaNa水型，近年來開采較多。西側(cè)、南側(cè)數(shù)公里范圍內(nèi)分布著較多的粉絲廠和淀粉廠。近年來，EW測向的西電極附近也分布多家粉絲廠，對測區(qū)環(huán)境有一定影響。昌黎臺所處地質(zhì)環(huán)境及斷層分布如圖1所示。

圖1　昌黎臺空間分布圖

1997年7月開始采用數(shù)字地電儀ZD8B觀測，背景場改造更換為ZD8M觀測，布極方式為2個相互垂直的測道NS、EW，AB/2=500 m，NS向測量極距MN=250 m，EW向測量極距MN=200 m。鉛板電極接地良好，外線路采用屏蔽絕緣銅線觀測。

2響應系數(shù)的計算

對于電性剖面參數(shù)已知的臺站，為確定多極距觀測的最佳布極方式，首先應計算各層介質(zhì)的響應系數(shù)Si隨觀測極距的變化情況。昌黎臺初始電性結(jié)構(gòu)模型為4層，第1層h1=6 m,ρ2=45 Ω·m；第2層h2=24 m,ρ2=90 Ω·m；第3層h3=85 m,ρ3=55 Ω·m；第4層h4=∞,ρ4=1 300 Ω·m[14-15](圖2)。電測深曲線為KH型。

圖2　昌黎臺初始電性結(jié)構(gòu)模型示意圖

以水平層狀介質(zhì)及其響應系數(shù)理論作為昌黎臺觀測資料反演的基本理論和模型，計算公式參見文獻[2]。對于理想的n層水平介質(zhì)，利用四極對稱裝置觀測時，得到的地電阻率相對變化量Δρs/ρs由公式：Δρs/ρs=s1×Δρ1/ρ1+s2×Δρ2/ρ2+…+sn×Δρn/ρn決定，s1，s2,…sn為響應系數(shù)，各層系數(shù)之和為1，每層介質(zhì)對Δρs/ρs的貢獻量取決于介質(zhì)本身的電阻率變化量和這一深度上的響應系數(shù)。響應系數(shù)大，介質(zhì)電阻率變化對地表觀測值影響也大；響應系數(shù)小，對地表觀測值的影響也越小。它實際上反映了四極對稱裝置在某種電性結(jié)構(gòu)介質(zhì)上的探測能力，其曲線就是水平層狀介質(zhì)不同深度對地表觀測電阻率變化的響應特征曲線[16]。

經(jīng)計算，得到昌黎臺不同極距情況下各層響應系數(shù)隨極距的變化情況(圖3)。

3裝置布設

主要包括：

1)極距個數(shù)。要反演計算n層電阻率值，觀測資料的數(shù)據(jù)應該等于或大于n[17]。一般來說，地電臺站實際電性結(jié)構(gòu)和理想的水平層狀模型有一定差異，因此，對于n層介質(zhì)來說，極距個數(shù)應為n+1或n+2。昌黎臺電性結(jié)構(gòu)模型為4層，結(jié)合實際觀測環(huán)境，極距個數(shù)設計為5個。

圖3　昌黎臺地下介質(zhì)響應系數(shù)分層曲線

2)極距長度分為供電極距AB和測量極距MN。電阻率隨供電極距逐漸增大，探測深度也越深，因此AB的取值范圍與響應系數(shù)密切相關(guān)。從圖2可以看出，對于昌黎臺來說，MN取AB/3、AB/4、AB/5的各層的響應系數(shù)的變化趨勢是非常一致的。對于多極距觀測，主要是捕捉或反映更深層的電阻率變化情況。當供電極距不大于1 000 m時，第4層的響應系數(shù)小于0.2，即對觀測到的視電阻率的貢獻不可能超過20%；而對于第三層來說，響應系數(shù)最高值在AB=300～400 m，當AB大于300 m時，則其貢獻可以達到70%，因此，在設計觀測極距時最主要是要顧及到第3層的貢獻；第2層響應系數(shù)最高值在AB=100 m左右，最高貢獻率60%。為了使觀測裝置也能反映淺部的變化，短極距選擇時可考慮1～2個。另外，實施中充分考慮了當?shù)丨h(huán)境的復雜性和多變性，既能利用原有裝置系統(tǒng)，又能合理埋設新的觀測裝置。當響應系數(shù)曲線極值點附近與直線偏離較大時，相應的極距應選在對應于響應系數(shù)曲線極值點附近的線性度較好的一段曲線上；當響應系數(shù)曲線極值點附近與直線偏離較小時，極距可選在響應系數(shù)曲線極值點上[17]。5個極距長度分別設計為：

NS向①AB=30 m，MN=6 m；②AB=125 m，MN=30 m；③AB=400 m，MN=80 m；④AB=500 m，MN=125 m；⑤AB=1 000 m，MN=250 m。

EW向①AB=30 m，MN=6 m；②AB=125 m，MN=30 m；③AB=300 m，MN=80 m；④AB=400 m，MN=125 m；⑤AB=1 000 m，MN=200 m。其中，2個測向1 000 m的極距布設利用原有裝置。

實際地電臺站進行多極距布設時，多受到觀測環(huán)境和建筑設施的影響，為了方便施工和維護工作，在節(jié)約費用的基礎上，充分考慮好布極方法是十分重要的，尤其提取貢獻層較大的電阻率時，對反演精度有一定影響。

4探測深度計算

關(guān)于地電阻率探測深度的定義，經(jīng)查文獻資料，認為其與物探電法中的探測深度含義有本質(zhì)區(qū)別。較早文獻有：1938年Evjen[18]認為物探中探測深度的目的，在于將地表的觀測值與某個具體深度聯(lián)系起來，使測量參數(shù)與地質(zhì)目標最相關(guān)，探測深度h與供電極距AB的關(guān)系最佳，為h=AB/9；1969年，AdelARZ和DallasBJ[19]在地下水勘探中，根據(jù)對多層電測深曲線的分析，認為介于AB/3～AB/4之間；而現(xiàn)在教學中的探測深度一般選擇AB/4 ～AB/6；1982年，趙和云等[20]定義視電阻率探測深度為：在一定的地電斷面條件下,用一定的裝置(包括裝置的類型和裝置系數(shù)2個方面)探測時,所具有探測深部介質(zhì)電阻率變化的能力；1986年，傅良魁[21]定義：在給定的各種主觀、客觀條件下，通過對電或電磁場信息的收錄、處理和分析研(正演，反演)，以查明探測目標的最大深度，主要影響因素包括主觀、客觀兩大方面等多種復雜情況。

該定義對現(xiàn)實觀測介質(zhì)多層分布的復雜性有較好的說服力，且當測量極距對與供電極距對互換時，公式的物理意義仍較清晰。公式推導過程中，采用的電流密度jMN和面電流分布IMN稱之為“等效”值，本身沒有物理意義。對于昌黎臺N=4層的復雜情況，由于電流在界面處可能出現(xiàn)很大的躍變，因此jMN應當分別在不同層內(nèi)計算，這樣推導過程十分復雜；加之實際各層介質(zhì)隨著深、淺層水位等復雜情況的變化，其實際探測深度會發(fā)生不同程度的變化，理論推導公式的物理意義已經(jīng)不大。本文為了簡化計算得到近似深度值，可將4層介質(zhì)大致看成2層，以基巖上層面為分界線，第1層h1=115m,ρ2=63.3Ω·m，h2=∞, ρ2=1 300Ω·m探測深度公式詳見文獻[20]。計算結(jié)果見表1。

表1　昌黎臺多極距對探測深度計算結(jié)果

5結(jié)論

利用昌黎臺初始電性結(jié)構(gòu)模型，對不同深層地下介質(zhì)的響應系數(shù)進行計算分析。當MN取AB/3、AB/4、AB/5時，各層響應系數(shù)的變化趨勢非常一致；當供電極距不大于1 000 m時，第4層的響應系數(shù)小于0.2，即對觀測到的視電阻率的貢獻不可能超過20%；而對于第3層來說，只要當AB大于300 m，則其貢獻可以達到70%，是多極距觀測裝置中主要考慮的布設方式；第2層貢獻率最高不會超過60%。

昌黎臺多極距觀測裝置的布極個數(shù)設計為5個。各供電極距長度選取時充分考慮了測區(qū)環(huán)境及建筑設施，為了方便施工和維護工作，在現(xiàn)有觀測裝置的前提下，依據(jù)不同深層介質(zhì)響應系數(shù)不同的特點，采取了較為合理的布極方式進行鋪設，盡最大可能捕捉到貢獻率較高的介質(zhì)層電阻率變化特征。5種布極方式計算出的探測深度在14～96 m左右，主要由介質(zhì)電性結(jié)構(gòu)決定，它取決于如下條件：①被探測地質(zhì)體本身的因素，如：大小、形狀、埋深及與圍巖的電阻率差別；②供電電極距(AB)的大??；③觀測精度；④地形和不均勻體的干擾；⑤外來電場的干擾。另外，不同的計算方法可能略有差異，對此問題的探討將另文敘述。

昌黎臺初始電性結(jié)構(gòu)模型的建立，為以后相關(guān)研究奠定了基礎，但不一定是最佳電性結(jié)構(gòu)模型，需要在以后的多極距觀測基礎上做大量實驗研究和反演計算。研究中，我們可以利用初始電性結(jié)構(gòu)模型計算獲得視電阻率，將其與相應的觀測極距對下獲得的實際觀測值相減，獲得電阻率殘差；然后進行反演，獲得模型參數(shù)的調(diào)整量；再與初始電性結(jié)構(gòu)模型的參數(shù)相加，即可獲得新的初始電性結(jié)構(gòu)模型參數(shù)。如此反復，直至獲得滿足條件的最佳模型參數(shù)[22]。

參考文獻：

[1]國家地震局《一九七六年唐山地震》編輯組. 一九七六年唐山地震[M]. 北京: 地震出版社, 1982:28-50.

[2]周海濤, 郭良遷, 張立成. 唐山斷裂現(xiàn)代活動性研究[J]. 華北地震科學, 2009, 27(3): 17-22.

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

[4]張國苓，喬子云，賈立峰，等.陽原地電阻率資料變化分析[J].震災防御技術(shù)，2015，10(2)：464-471.

[5]錢家棟, 陳有發(fā), 金安忠. 地電阻率法在地震預報中的應用[M]. 北京: 地震出版社, 1985: 82-90.

[6]趙和云, 錢家棟. 郫縣臺多極距觀測資料的反演和分析[J]. 中國地震, 1987, 3(增刊): 73-78.

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

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

[9]聶永安, 姚蘭予. 成層半空間深埋電極產(chǎn)生的電位分布[J]. 中國地震, 2009, 25(3): 246-255.

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

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

[12]孟凡月，曹均鋒，李光.電阻率層析成像技術(shù)在地震安全性評價中的應用[J].華北地震科學，2012,30(2)：52-55.

[13]張德信, 王文禮. 我國部分地電臺基本情況表[C]//錢復業(yè). 地震預報——地電方法論文集.福州: 福建科學技術(shù)出版社, 1985: 318-321.

[14]馮銳, 郝錦綺, 周建國. 地震監(jiān)測中的電阻率層析技術(shù)[J]. 地球物理學報, 2001, 44(6): 833-842.

[15]張學民, 翟彥忠, 郭建芳, 等. 2004年印尼8.7級強震前昌黎臺的遠場電性異常特征探討[J]. 地震, 2006, 26(4): 82-93.

[16]趙和云, 張文孝, 楊明芝. 水平層狀介質(zhì)中不同深度對地表ρs變化的響應特性[J]. 地震, 1992, 14(6): 3-10.

[17]李艷東, 趙家騮, 龐麗娜. 地電阻率多極距觀測中的最佳布極方式[J]. 地震, 2004, 24(2): 17-24.

[18]Evjen H M. Depth factors and resolving power of electrical measurements[J]. Geophysics, 1938, 3(2): 78-95.

[19]ZohdyAA R, Jackson D B. Application of deep electrical soundings for groundwater exploration in Hawaii[J]. Geophysics, 1969, 34(4): 584-600.

[20]趙和云, 錢家棟. 地電阻率法中勘探深度和探測范圍的理論討論和計算[J]. 西北地震學報, 1982, 4(1): 40-56.

[21]傅良魁. 電法勘探文集[M]. 北京: 地質(zhì)出版社, 1986: 30-72.

[22]朱濤, 王蘭煒, 胡哲, 等. 電阻率多極距觀測臺站的電性結(jié)構(gòu)模型及參數(shù)反演研究——以山丹臺為例[J]. 工程地球物理學報, 2011, 8(1): 16-23.

The Design of Multi-separation Array Observation System of Changli Geoelectricity Station

TONG Xin, GUO Jian-fang*, ZHOU Jan-qing, ZHANG Si-yuan

(Qinhuangdao Central Seismic Station of Earthquake Administration of Hebei Province, Qinhuangdao 066100, China)

Abstract:This study analyze the response coefficient, theoretical apparent resistivity of different deep medium and the depth of probe based on 1-dimension model of electrical structure in Changli station. The results show that when MN is 1/3AB, 1/4AB, 1/5AB respectively, the variation trends of response coefficient of each layer are very consistent. Because the contribution rate of third layer which achieved 70 percent is highest in four layers of different depth, the contribution rate of the third layer is the dominant factor in arranging the electrodes for multi-spacing system in apparent resistivity observation. When the scale of electrode distance of the third layer is from 300 to 500 meters, the results are best. The third layer contribution rate is no more than 60%. When the electrode distance is less than 1000 meters, the contribution rate of the forth layer is just 20 percent. According to the distribution of the response coefficient of the deep layer and the actual observation conditions of the station, the multi-separation array layout is 5 kinds: AB=1 000 m, 500(or 400)m, 400(or 300)m, 125 m, 30 m, the detecting depths is around 14-96 m, the results rely on the electrical structure of medium under the Changli station.

Key words:initial electrical structure; response coefficient; multi-separation array layout; Changli geoelectricity station

doi:10.3969/j.issn.1003-1375.2016.01.011

中圖分類號:P315.75

文獻標志碼:A

文章編號:1003-1375(2016)01-0065-05

作者簡介：佟鑫(1987—)，男，吉林白城人，助理工程師，從事電磁監(jiān)測、分析預報和儀器維修工作.E-mail:278507679@qq.com*通訊作者：郭建芳(1977—)，女，河北昌黎人，高級工程師，從事前兆臺網(wǎng)數(shù)據(jù)處理及分析預報工作.E-mail: gjfzyy76@sina.com

基金項目：中國地震局科技星火計劃“昌黎臺地電阻率多極距觀測實驗與地下水位關(guān)系研究”(XH15005Y)

收稿日期:2015-06-22

佟鑫，郭建芳，周劍青，等．昌黎地電臺多極距觀測系統(tǒng)的設計[J]．華北地震科學，2016，34(1):65-69.