張吉濤

(山東省第四地質(zhì)礦產(chǎn)勘查院)

高密度電阻率法是一種以地下巖土介質(zhì)的電性差異為基礎(chǔ)，通過觀測人工構(gòu)建的地下穩(wěn)定電流場的分布規(guī)律以解決地下地質(zhì)問題的地球物理勘探方法[1-4]。與常規(guī)電阻率法相比，該方法具有點距小、數(shù)據(jù)采集密度大、成本低、效率高、信息豐富、解譯方便等優(yōu)點[5-8]，近年來，該方法在工程勘查領(lǐng)域得到了廣泛應(yīng)用，取得了顯著的經(jīng)濟和社會效益[9-14]。本研究通過詳細梳理該方法的技術(shù)原理、特點、裝置形式及數(shù)據(jù)處理方法，并結(jié)合工程實例，進一步探討該方法在地下水探測、采空區(qū)勘查中的應(yīng)用。

1 高密度電阻率法概述

1.1 技術(shù)原理

高密度電阻率法是一種陣列勘探方法，實際上是一種集電阻率剖面法和電阻率測深法于一體的多裝置、多極距組合方法，其原理與普通電阻率法一致，區(qū)別在于觀測中設(shè)置了高密度的觀測點。該方法是一類以巖土體導電性差異為基礎(chǔ)的電法勘查方法，當人工向地下施加直流電流時，在地表利用儀器觀測其電場分布，將所測電流和電壓值換算為電阻率值。但由于受到地形起伏和地下介質(zhì)不均勻性等因素的影響，測量所得的電阻率值并非巖石的電阻率真值，而為視電阻率值。通過進一步研究視電阻率值的差異和變化規(guī)律，有助于發(fā)現(xiàn)和勘查地下目標體，從而達到解決工程地質(zhì)問題的目的。

1.2 特 點

采用高密度電阻率法進行野外測量時可將全部電極置于觀測剖面的各測點上，利用程控電極轉(zhuǎn)換裝置和微機工程電測儀便可實現(xiàn)不同極距、不同電極排列方式的數(shù)據(jù)快速自動采集。與常規(guī)電阻率法相比，該方法的特點為：①電極布設(shè)一次性完成，減少了因電極設(shè)置而引起的故障和干擾，大大提高了工作效率；②可進行多種電極排列方式測量，獲取的數(shù)據(jù)量大、信息豐富，能夠獲得較豐富的有關(guān)地電斷面的地質(zhì)信息；③野外數(shù)據(jù)采集實現(xiàn)了自動化或半自動化，不僅采集速度快，而且減輕了勞動強度；④觀測精度和分辨率顯著提高，探測深度更加靈活；⑤可以對探測數(shù)據(jù)進行現(xiàn)場實時處理并顯示剖面曲線形態(tài)，數(shù)據(jù)解譯更為方便。

1.3 裝置形式

高密度電阻率法的裝置形式由最初的α、β、γ3種方式發(fā)展至目前的十余種[6,15]，主要包括二極裝置(Pole-Pole)、單邊三極裝置(Pole-Dpole，Dpole-pole)、溫納裝置(Wenner)、偶極裝置(Dpole-Dpole)和斯隆貝格裝置(Schlumberger)等。該類裝置具有各自的優(yōu)缺點及相應(yīng)的應(yīng)用條件，因此在實際工作中應(yīng)根據(jù)實際情況進行合理選擇。根據(jù)相關(guān)實踐經(jīng)驗，四極裝置為公認的最穩(wěn)妥的裝置，比較常用的有溫納裝置和斯隆貝格裝置。

溫納裝置方式(Wenner)又稱為對稱四極裝置方式，電極排列方式為A、M、N、B等間距排列(A、B為供電電極，M、N為測量電極)，裝置系數(shù)K=2πna(a為電極間距，n為隔離系數(shù))，AM=MN=NB=na，隨著隔離系數(shù)按照由小到大的順序等間隔增大，4個電極之間的間距也被均勻拉開。該測量方式屬剖面測量方式，所得斷面呈倒梯形[16](圖1)。

圖1 溫納裝置跑極示意

斯隆貝格裝置方式(Schlumberger)屬測深測量方式，裝置系數(shù)K=πn(n+1)a(a為電極間距，n為隔離系數(shù))，MN=a，AM=NB=na。測量時，測量電極M、N保持不動，供電電極A、B同時逐點分別向左、向右移動，得到1條滾動掃描測量線，而后A、M、N、B同時向右移動1個電極，再按照同樣方式跑極，得到另1條滾動掃描測量線。按照該測量方式，可得到如圖2所示的矩形斷面。

圖2 斯隆貝格裝置跑極示意

1.4 數(shù)據(jù)處理方法

高密度電法儀器可以實現(xiàn)野外數(shù)據(jù)自動采集并存入主機，主機可通過通訊軟件將原始數(shù)據(jù)傳輸至計算機。野外數(shù)據(jù)采集過程中難免存在各類干擾，為獲取真實的探測信息，一般需對采集的原始數(shù)據(jù)進行預(yù)處理，常用方法為剔除突變點、數(shù)據(jù)平滑、地形校正和數(shù)據(jù)拼接等[17-19]。剔除突變點是將與相鄰測點的電阻率值有數(shù)十倍差異的測點數(shù)據(jù)進行剔除，而后進行插值，以壓制虛假異常。數(shù)據(jù)平滑是為了消除測量過程中的隨機干擾，可采用滑動平均等方法進行處理。地形校正是為了消除因地形起伏而形成的假異常。數(shù)據(jù)拼接主要是對長剖面的分段測量數(shù)據(jù)進行合并處理，形成1條完整剖面，便于對數(shù)據(jù)進行整體解譯。經(jīng)過預(yù)處理后的數(shù)據(jù)可利用Sufer軟件繪制出視電阻率等值線擬斷面圖。當生成的斷面圖分辨率不高、解譯較困難時，還可利用相應(yīng)的軟件進行反演計算成圖，便于對探測目標進行定性或半定量解譯。

2 實例分析

2.1 地下水探測

2.1.1 工程概況

山東省濰坊市經(jīng)濟發(fā)展迅速，工農(nóng)業(yè)生產(chǎn)中對水的需求量越來越大。自2013年8月開始，濰坊市經(jīng)歷了一段有水文記載以來最嚴重的持續(xù)干旱時期。為應(yīng)對突發(fā)飲用水供水危機，濰坊市政府決定開展應(yīng)急供水水源地勘查評價工作。由水文地質(zhì)調(diào)查資料知，勘查區(qū)內(nèi)主要富水地段的地下水類型為松散巖類孔隙水和基巖裂隙水。松散巖類孔隙水位于第四系覆蓋區(qū)，主要含水層位為古河道沉積形成的砂層，與周圍黏性土壤相比，砂層孔隙度大，導電性差，視電阻率值明顯高于周圍黏土的視電阻率值?；鶐r裂隙水主要賦存于玄武巖和砂巖的風化裂隙中，裂隙巖層通常由部分裂隙在巖層中局部范圍內(nèi)連通構(gòu)成帶狀或脈狀裂隙含水系統(tǒng)，其視電阻率值較圍巖明顯下降。因此，利用高密度電阻率法尋找地下水具備地球物理前提。

2.1.2 野外數(shù)據(jù)采集與處理

本研究地下水探測工作采用吉林大學工程技術(shù)研究所生產(chǎn)的E60M型高密度電法儀，鑒于野外工作場地開闊且地形較平坦，故優(yōu)先考慮四極裝置形式，以便獲得最大的測量電位，從而有效壓制干擾，增強有效信號。但經(jīng)現(xiàn)場試驗，溫納裝置信噪比最高，具有較好的垂向分辨率，對于解決垂向變化(例如水平層狀結(jié)構(gòu))問題比較有利，因此本研究最終采用溫納裝置方式進行測量。本研究所用的高密度電法儀中，電纜上的電極間距呈10 m均勻分布，經(jīng)現(xiàn)場試驗?zāi)軌驖M足野外工作需要。輸入電壓取12 V，供電脈沖寬度設(shè)置為2 s。為保證野外數(shù)據(jù)采集精度，每天工作前需對儀器、電瓶進行充分充電。在正式測量前，應(yīng)對每個電極開關(guān)和電極接地電阻進行檢測，當接地電阻過大時，應(yīng)采取措施(如加水或有機活性物質(zhì)等)改善電極的接地條件，確保探測資料準確可靠。

將野外采集的數(shù)據(jù)傳輸至計算機后，需對數(shù)據(jù)進行預(yù)處理和反演處理。數(shù)據(jù)預(yù)處理包括剔除異常點、濾波去噪等，為反演處理提供穩(wěn)定數(shù)據(jù)體。采用瑞典RES2DINV電法反演軟件進行反演處理，網(wǎng)格采用“有限元”方法，最優(yōu)化阻尼系數(shù)，利用線性搜索的最小二乘反演方法，確保反演結(jié)果的穩(wěn)定性，迭代次數(shù)為2～5次，設(shè)置均方根誤差門限值為2.5%。處理后的數(shù)據(jù)可采用Surfer軟件繪制出視電阻率斷面圖。

2.1.3 數(shù)據(jù)解譯與驗證

圖線剖面高密度電法綜合斷面

分析圖3可知：視電阻率斷面總體呈水平層狀展布，橫向連續(xù)性較好，縱向上視電阻率變化較大，表現(xiàn)為高低阻相間特征。淺部(標高小于10 m)視電阻率值整體較小，一般為10～20 Ω·m，推斷為表層粉土的電性反映；淺部視電阻率值沿縱向延深不均勻，剖面兩端向下延深較深，中部向下延深較淺，反映了粉土厚度的不均勻性。中部(標高為10～35 m)視電阻率值為20～40 Ω·m，整體呈高阻特征，推斷由砂層引起，因砂層孔隙度較大，致使其導電性較差，視電阻率值相對較大，該層為主要的富水層位。深部(標高大于35 m)視電阻率值隨著深度增加逐漸變小，其值小于20 Ω·m，呈低阻特征，推斷為粉質(zhì)黏土的電性反映。

為驗證高密度電法探測成果的可靠性，本研究在該剖面附近布置了2個水文地質(zhì)鉆孔。鉆探揭露結(jié)果為：標高0～10.4 m區(qū)段為粉土；標高10.4～32 m區(qū)段為砂層與粉質(zhì)黏土互層，其中，10.4～14.4 m 區(qū)段為粉砂層，25.6～27.2 m區(qū)段為礫砂層，30.4～32 m區(qū)段為粗砂層，該3層砂層為主要的含水層，粉質(zhì)黏土中含少量砂質(zhì)，為弱含水層；標高32～82 m區(qū)段為粉質(zhì)黏土，僅在標高53，69 m附近存在的薄層粗砂中含水。由此可見，鉆孔資料與高密度電法解譯成果基本吻合，后經(jīng)抽水試驗，單井涌水量均在24.5 m3/h左右。

2.2 采空區(qū)勘查

2.2.1 工程概況

山東省乳山市金碃嶺金礦開采歷史悠久，目前礦區(qū)已由早先的城郊農(nóng)村擴充至城市規(guī)劃核心區(qū)，礦山閉坑后遺留的采空區(qū)已成為較大的災(zāi)害源，影響了城市發(fā)展。該區(qū)采礦形成的采空區(qū)可大致分為2個不同時期：一是早期民采形成的采空區(qū)；二是20世紀70年代開采硫鐵礦時形成的采空區(qū)；三是20世紀90年代開采金礦時形成的采空區(qū)。民采和硫鐵礦開采期間的開采深度較淺，一般不超過30 m，但開采隨意性大，采空區(qū)分布無規(guī)律，地表采坑和塌陷的采空區(qū)已被回填，僅局部遺留有采空區(qū)。金礦開采期間采用地下巷道方式，有正規(guī)的巷道開拓圖，巷道深度為地表以下30～70 m，該時期的巷道及采空區(qū)勘查以巷道及采空區(qū)水平投影圖為基礎(chǔ)，輔以物探勘查和鉆探驗證。

勘查區(qū)內(nèi)的地層主要為新生界第四系臨沂組，巖性為灰黃色黏土質(zhì)粉砂、含礫中細砂。巖漿巖大面積分布，巖性主要為中生代燕山早期垛崮山序列大孤山單元斑狀中細粒含黑云花崗閃長巖，此外，還發(fā)育有燕山晚期花崗斑巖脈和煌斑巖脈。區(qū)內(nèi)斷裂構(gòu)造發(fā)育，有4條大致平行分布的斷裂構(gòu)造蝕變帶，走向2°～10°，傾向NW，傾角80°～85°，斷裂帶沿走向呈舒緩波狀延伸，局部膨大、分枝復(fù)合現(xiàn)象明顯。理論上講，當?shù)V體被采空后形成一定的充氣空間，破壞了巖體的完整性和連續(xù)性，造成采空區(qū)與圍巖有顯著的電性差異，表現(xiàn)出明顯的高阻特征。當經(jīng)過一段時間后，采空區(qū)上覆巖石在重力作用下會發(fā)生坍塌變形，致使巖石破碎并出現(xiàn)裂縫，地下水會沿著破碎巖石和裂縫向采空區(qū)匯集，采空區(qū)充水后，其電阻率會呈低阻反映。該類特性成為應(yīng)用高密度電阻率法探測地下采空區(qū)良好的地球物理前提。

2.2.2 野外數(shù)據(jù)采集與處理

仍采用E60M型高密度電法儀進行野外數(shù)據(jù)采集，經(jīng)現(xiàn)場試驗，以斯隆貝格裝置方式進行測量，點距5 m，個別剖面又選擇2，3 m點距進行了驗證測量，輸入供電電壓為12 V，供電脈沖寬度為2 s，控制深度40～80 m不等。將采集的數(shù)據(jù)進行預(yù)處理后，采用瑞典高密度電法處理軟件RES2DINV進行反演計算并成圖。

2.2.3 數(shù)據(jù)解譯與驗證

圖線視電阻率等值線斷面

3 結(jié) 語

高密度電阻率法以地下介質(zhì)的電性差異為基礎(chǔ)，具有自動化程度高、工作效率高、信息豐富、解譯方便等特點。結(jié)合工程實例，通過數(shù)據(jù)解譯并輔以鉆探驗證的方式，詳細討論了該方法在地下水探測、采空區(qū)勘查等方面的適用性，對于類似工程有一定的借鑒價值。

圖5 ZK27鉆孔巖芯

[1] 劉國興.電法勘探原理與方法[M].北京：地質(zhì)出版社，2005.

[2] 劉 斌，張光保.高密度電法在隧道涌水通道勘查中的應(yīng)用[J].工程地球物理學報，2012，9(6)：750-754.

[3] 馬法成.某鐵礦山隱伏采空區(qū)精準探測方法[J].金屬礦山，2016(5)：196-199.

[4] 顏秉超.某礦山采空區(qū)探測的高密度電阻率法層析成像技術(shù)[J].現(xiàn)代礦業(yè)，2017(8)：250-252.

[5] 鄧文超，周孝宇.高密度電法的原理及工程應(yīng)用[J].西部探礦工程，2006(S)：278-279.

[6] 董浩斌，王傳雷.高密度電法的發(fā)展與應(yīng)用[J].地學前緣，2003，10(1)：171-176.

[7] 呂惠進，劉少華，劉伯根.高密度電阻率法在地面塌陷調(diào)查中的應(yīng)用[J].地球物理學進展，2005，20(2)：381-386.

[8] 張先林，許 強，彭大雷，等.高密度電法在黑方臺地下水探測中的應(yīng)用[J].地球物理學進展，2017，32(4)：1862-1867.

[9] 朱紫祥，胡俊杰.高密度電法在巖溶地區(qū)溶洞勘查中的應(yīng)用[J].工程地球物理學報，2017，14(3)：290-293.

[10] 李洪嘉，閆紹波，張 超.綜合物探技術(shù)在煤礦采空區(qū)探測中的應(yīng)用研究[J].工程地球物理學報，2014，11(5)：714-720.

[11] 王士鵬.高密度電法在水文地質(zhì)和工程地質(zhì)中的應(yīng)用[J].水文地質(zhì)工程地質(zhì)，2000(1)：52-56.

[12] 劉曉東，張虎生，朱偉忠.高密度電法在工程物探中的應(yīng)用[J].工程勘察，2001(4)：64-66.

[13] 劉宏岳.直流高密度電法在淺海工程勘察中的應(yīng)用[J].物探與化探，2013，37(4)：756-760.

[14] 蘇兆鋒，陳昌彥，肖 敏，等.精細高密度電阻率法在白云巖礦采空區(qū)中的應(yīng)用[J].物探與化探，2012(S)：45-47.

[15] 趙光輝.高密度電法勘探技術(shù)及其應(yīng)用[J].礦產(chǎn)與地質(zhì)，2006，20(2)：166-168.

[16] 馬志飛，劉鴻福，葉 章，等.高密度電法不同跑極方式的對比及效果分析[J].工程地質(zhì)計算機應(yīng)用，2008(3)：11-15.

[17] 嚴加永，孟貴祥，呂慶田，等.高密度電法的進展與展望[J].物探與化探，2012，36(4)：576-584.

[18] 李銀真.高密度電阻率法物探技術(shù)及其應(yīng)用研究[D].阜新：遼寧工程技術(shù)大學，2007.

[19] 李學敏，方玉滿，李 爭，等.高密度電阻率法在滑坡體勘探中的應(yīng)用[J].現(xiàn)代礦業(yè)，2012(11)：37-38.