高密度電法在巖溶地區(qū)溶洞探測中的應用

別發(fā)禮 朱蒙蒙 何云飛

(1.宜昌星興藍天科技有限公司，湖北 宜昌 443311； 2.三峽大學，湖北 宜昌 443300)

1 工程概況

擬建項目位于湖北省神農(nóng)架國際滑雪場，場區(qū)內(nèi)原有造雪用水水庫由一土壩攔水而成。位于滑雪場接待中心南部600 m，為粘土心墻壩。由于雪場造雪用水量逐年增加，須進一步擴大水庫蓄水量。而原有水庫不能滿足用水量需求。因此規(guī)劃建設神農(nóng)架國際滑雪場土壩擴建工程，該工程主要是加高大壩8.6 m，加高后壩頂高程為2 221.5 m，同時興建輸水設備，改建溢洪道等，增加有效蓄水。工程前期開展了地質(zhì)勘查工作，在左岸壩肩位置的鉆孔ZK5處發(fā)現(xiàn)溶洞。而地下溶洞為工程建設帶來了較大的安全隱患，故需要探測地下溶洞的規(guī)模及分布。溶洞發(fā)育的不規(guī)律性，為勘察工作帶來了一定的難度。為了調(diào)查溶洞大小規(guī)模以及分布情況，采用高密度電法結合地質(zhì)資料和鉆探資料的物探手段對該工程范圍的地質(zhì)情況進行勘查。以此勘測結果為土壩工程設計、施工提供數(shù)據(jù)資料。

2 場區(qū)地質(zhì)條件概述

在地質(zhì)方面，場區(qū)位置在構造格架上均屬揚子準地臺。地處神農(nóng)架穹窿中部，由一系列北西向壓性為主的斷裂及褶曲組成，構造形跡以木魚坪向斜為主，向斜較開闊，區(qū)域地殼穩(wěn)定性好。場區(qū)出露的地層巖性為元古界神農(nóng)架群中亞群石槽河組(Pts)白云巖、砂巖，巖層走向在庫區(qū)為N20°W，產(chǎn)狀近于垂直。由地層巖性組合結合地下水賦存條件綜合分析可知，場區(qū)主要含水巖組為第四系孔隙含水巖組，基巖中基本不含水，構成相對隔水層。局部裂隙較發(fā)育帶含水，水量較少。

3 高密度電法工作原理

高密度電法是通過電極向地下供電，從而形成人工電場以測出地層的電阻率，以地層的導電性差異來分析地層結構的方法。它是一種陣列勘測的方法，是在直流電法的基礎上發(fā)展出來的。其功能相當于四極測深與電剖面法的結合，故此提高了電法勘測的效率，同時自動化程度及勘測精度有較大提升。

當巖體完整性較差，且?guī)r石裂隙中含水時，其視電阻率明顯較低。而完整性好且堅硬的巖體，其視電阻率明顯高于斷層破碎帶及富水帶圍巖的視電阻率。這種對圍巖的含水情況特別敏感探測方法，能夠較好的反映出巖溶地區(qū)溶洞的大小及分布情況。高密度電法反演出的圖像直觀、原理清晰、分辨率高等優(yōu)點，使它成為一種較好的物探方法。

高密度電法現(xiàn)場勘測時，先是沿需要觀測的剖面一次性等間距布好電極,電極與多芯電纜連接、多芯電纜與轉換裝置連接、轉換裝置與電測儀相連。通電測量時，轉換開關可按照特定的裝置方式，通過特定的電極轉換裝置，依次測量該測線上不同位置、不同深度處的視電阻率值。再將采集的數(shù)據(jù)傳送至計算機中，采用專用的高密度電法處理軟件系統(tǒng)計算、分析、成圖，從而得到一條完整的二維視電阻率剖面。其工作系統(tǒng)如圖1所示。

3.1 裝置選擇

白云質(zhì)灰?guī)r等巖石的視電阻率較高，而黏土等物質(zhì)的視電阻率較低。這種不同物質(zhì)間的電性差異，是使用高密度電法研究溶洞和溶蝕的先決條件。高密度電法常用的測量裝置主要是溫納和施倫貝爾裝置。本項目通過在測區(qū)內(nèi)試測，經(jīng)比較分析，選擇溫納裝置為主要測量裝置。

3.2 溫納裝置原理

溫納裝置又稱為對稱四級裝置方式。A，M，N，B等間距排列，其中A，B是供電電極，M，N是測量電極。其視電阻率表達式為：

其中，k為裝置系數(shù)。

測量時，AM=MN=NB=AB/3為一個電極間距，探測深度為AB/3，A，B，M，N逐點同時向右移動，得到第一層剖面；接著AM，MN，NB增大一個電極間距，A，B，M，N逐點同時向右移動，得到另一層剖面；這樣不斷掃描測量下去，得到倒梯形斷面(跑極方式見圖2)。

4 場區(qū)布置測線

在場區(qū)內(nèi)總共布設三條測線，編號為A,B,C。A—A′測線平行于壩軸線，經(jīng)過鉆孔ZK5位置，電極數(shù)30，電極間距5 m，測線長150 m；B—B′測線垂直于壩軸線，沿著左岸邊坡縱向布置，經(jīng)過鉆孔ZK5位置且保持鉆孔ZK5孔口高程向下游延伸，電極數(shù)60，電極間距5 m，測線長300 m；C—C′測線：垂直于壩軸線，沿著左岸邊坡縱向布置，一半測線(上游)與B—B′測線重疊且經(jīng)過鉆孔ZK5位置，另一半測線沿左岸邊坡坡腳向下游延伸，電極數(shù)60，電極間距5 m，測線長300 m。

5 勘探數(shù)據(jù)的解算

5.1 高密度電法剖面的解譯方法

本次數(shù)據(jù)反演采用RES2DINV高密度電阻率反演程序，它是基于圓滑約束最小二乘法的計算機反演程序。它使用了基于準牛頓最優(yōu)化非線性最小二乘新算法，在大數(shù)據(jù)下的計算速度較常規(guī)最小二乘法快很多。在原始資料的處理的基礎上，通過建立適宜的計算機反演模型，輸入處理好的數(shù)據(jù)文件后，計算機自動進行反演，其反演次數(shù)由模型的均方誤差來決定，一般是選取迭加后均方誤差不再明顯改變的模型，一般在反演3次～5次可達到目的。經(jīng)解算得到三條測線的電阻率分布圖如圖3～圖5所示。

5.2 解譯成果分析

1)從測線A—A′所解譯的電阻率剖面圖中可以發(fā)現(xiàn)在新的壩軸線中間部分是大片的低阻區(qū)，電阻率小于100 Ω·m，這片低阻區(qū)之所以表現(xiàn)出低阻特性，應當是現(xiàn)有蓄水水庫向下游滲透形成高含水區(qū)，其表現(xiàn)出兩個分布特征：一是深層電阻率比表層電阻率更低，第二是壩軸線中部電阻率比兩邊壩肩位置低。

2)在測線A—A′電阻率分布圖壩軸線靠左岸壩肩位置，出現(xiàn)了一片高阻區(qū)域，電阻率超過1 600 Ω·m，這片高阻區(qū)表明該區(qū)域的含水量應顯著低于周邊地層，而此處正是鉆孔ZK5發(fā)現(xiàn)溶洞的位置，由于溶洞空腔電阻率高于周邊區(qū)域，因此可以判定該高阻區(qū)域應是溶洞范圍。溶洞在A—A′剖面上的分布范圍應自鉆孔ZK5向主河道延伸20 m～25 m，溶洞高5 m左右。溶洞深度信息受局部地形影響(局部突然變陡),不能準確判讀，大致應在主河道地面線以下10 m附近。受兩邊地形的影響，溶洞向左岸邊坡內(nèi)部延伸情況難以評估。

3)觀察測線B—B′和C—C′電阻率分布圖，可以發(fā)現(xiàn)左岸土壩上游部分為連續(xù)的低阻區(qū)，下游部分是連續(xù)的高阻區(qū)，在兩個區(qū)域交界的位置出現(xiàn)溶洞。B—B′測線剖面顯示這一區(qū)域可能存在多個溶洞。溶洞大小寬2 m～10 m不等，高2 m～5 m。從圖上量測的溶洞埋深距鉆孔ZK5孔口20 m左右，與鉆孔情況相符(如圖6所示)。

4)綜合在三個勘探剖面中發(fā)現(xiàn)的溶洞跡象，判斷溶洞走向與壩軸線平行的可能性較大。

5)根據(jù)溶洞區(qū)域的相對高阻狀態(tài)，可判定溶洞填充狀態(tài)至少是欠填充。

6 結語

隨著計算機數(shù)據(jù)采集技術的發(fā)展，高密度電法勘探效率也得到很大提高。同時，它也向著更大的剖面覆蓋面積、增大探測深度、提高信噪比等方向發(fā)展。通過本次高密度電法在巖溶地區(qū)溶洞探測，并結合地質(zhì)、鉆探資料綜合分析，取得了顯著的效果。為類似的工程勘測提供了一個經(jīng)典的案例。