李遠強

（北京市地質(zhì)研究所，北京 100120）

孔間電磁波CT技術(shù)在北京普安店巖溶勘查中的應(yīng)用

李遠強

（北京市地質(zhì)研究所，北京 100120）

孔間電磁波CT技術(shù)是工程勘探中的一種常用方法, 利用孔間層析成像進行探測，其成果清晰、直觀,較好地展現(xiàn)了孔間介質(zhì)電磁波吸收情況,清楚地顯示出地下不同地質(zhì)體的空間分布。北京普安店地區(qū)存在巖溶塌陷，前期進行了大量的鉆孔勘探，利用鉆孔進行孔間電磁波CT掃描，探測鉆孔之間溶洞的分布情況，填補了鉆孔之間溶洞展布情況的空白，并對溶洞的充填情況進行解釋。該技術(shù)提供的數(shù)據(jù)準(zhǔn)確可靠，信息豐富，為后期巖溶地區(qū)的穩(wěn)定性評價及工程處理提供了科學(xué)依據(jù)和基礎(chǔ)數(shù)據(jù)。

孔間電磁波CT；巖溶勘查；數(shù)據(jù)采集處理

0 引言

孔間電磁波CT是用發(fā)射電磁波的方法對井間的地質(zhì)剖面進行層析技術(shù)處理，得到井間地質(zhì)剖面的內(nèi)部結(jié)構(gòu)圖像，該項技術(shù)分辨率高、野外施工方便，可以避開地面建筑。同時該方法充分利用已有鉆孔，填補鉆孔之間地質(zhì)信息空白。相對地面物探，其受到地面建筑及電磁干擾較小，目前廣泛應(yīng)用于斷層破碎帶、軟弱夾層、巖溶等地質(zhì)構(gòu)造方面的探測，如井間電磁波技術(shù)在碳酸鹽巖縫洞勘察中的應(yīng)用[1]、在居民區(qū)地下洞穴調(diào)查中的應(yīng)用[2], 在水庫堤壩隱患探測中的應(yīng)用[3].

1 研究區(qū)概況

北京香山普安店地區(qū)在近幾年發(fā)生了多起地面塌陷災(zāi)害，其中較嚴(yán)重的有兩次，分別發(fā)生于2009年7月25日和2010年7月3日。前者形成直徑約4m，深約5m的塌陷坑，造成一名人員掉入塌陷坑內(nèi)、一棵大樹陷入坑內(nèi)及一條污水管道受損折斷；后者形成塌陷坑大小約2.5m×2.5m，深3～4m，造成了2名人員受傷和部分財產(chǎn)損失。地面塌陷給當(dāng)?shù)鼐用裨斐煽只?，成為該地區(qū)的不安定因素和安全隱患，當(dāng)?shù)卣疄槊裆捕?，開展了該地區(qū)的地面塌陷災(zāi)害勘查工作。

研究區(qū)位于北京西山山前丘陵向平原過渡的區(qū)域，地勢平坦，地面標(biāo)高57～59m，地形坡度在1°～2°之間。屬于暖溫帶半濕潤大陸性季風(fēng)氣候，夏季極端最高氣溫35℃～40℃之間，冬季極端最低氣溫－11℃～－19℃；多年平均降水量613mm，年降水量最多達1406mm，夏季（6～8月）降水集中，占年降水量的75%，降水強度大，雨量往往集中在幾次降雨過程中，地面塌陷多發(fā)生于該時間段。研究區(qū)基巖被第四系（Q4）所覆蓋，為坡積、沖積、沖洪積、洪積物，以粘土、粉土、黃土夾砂礫石及砂礫石為主；基巖主要為奧陶系馬家溝組灰?guī)r（Om），均隱伏于第四系之下，基巖埋深最淺處小于10m。區(qū)內(nèi)發(fā)育有近NS向冷泉村平移斷裂（象鼻子溝斷裂），斷面傾向西，傾角80°～90°，東盤北移，西盤南移。

鉆孔揭示研究區(qū)灰?guī)r為灰色，強-中等風(fēng)化，中-厚層狀，層面傾角約30°～40°，節(jié)理裂隙發(fā)育，巖芯上可見溶孔及溶槽，溶蝕現(xiàn)象較發(fā)育，上部巖芯采取率低，漏漿；多個鉆孔揭示地下發(fā)育的溶洞，埋深17～50m，高度約0.1～3.7m不等，多數(shù)溶洞充填碎石、粘土等，少數(shù)無充填。鉆孔揭露的土體自上而下為雜填土、素填土、粘質(zhì)粉土、粉質(zhì)粘土等，下部含有碎石、卵石；鉆孔揭示地下發(fā)育的多個土洞，埋深0.6～37m，高度約0.8～5.3m不等，少數(shù)土洞充填碎石土，多數(shù)無充填。

2 孔間電磁波CT的基本原理

孔間電磁波CT的基本原理是在兩個鉆孔分別發(fā)射和接收無線電波(工作頻率0.5~32M Hz) , 根據(jù)不同位置上接收的場強大小,來確定地下不同介質(zhì)分布的一種地下地球物理勘查方法[4](圖1)。由于地下介質(zhì)具有不同物理性質(zhì)，在電磁波的傳播過程中主要表現(xiàn)在對電磁波能量的吸收, 使其能力衰減，這種吸收作用與地下巖土體的類別、裂隙分布、含水程度、充填物質(zhì)等因素有關(guān), 通過2個鉆孔之間電磁波掃描性觀測, 利用層析成像反演算法, 將不同巖性導(dǎo)致的電磁波能量上的差異分布轉(zhuǎn)變成二維分布圖像, 進而推斷出地下地層及構(gòu)造情況。該方法具有電磁波在地下有耗半空間的輻射、傳播和接收的特征，其正反演問題的理論基礎(chǔ)是電磁場理論。下式為地下電磁波法中的場強觀測值公式：

式中：E 為接收點的場強值；

E0′為初始輻射常數(shù)；

β為測區(qū)介質(zhì)吸收系數(shù), 即介質(zhì)中單位距離對電磁波的吸收值；

f(θ)為收發(fā)天線的方向因子函數(shù)；

r為發(fā)射與接收點之間的距離。

式(1) 表明通過e-βγf(θ) r–1因子，E0′衰減到E , 吸收系數(shù)（β）是一個與介質(zhì)電阻率（ρ）、介電常數(shù)（ε）、磁導(dǎo)率（μ）以及電磁波頻率數(shù)（ω）有關(guān)介質(zhì)的主要參數(shù)，它表示介質(zhì)對電磁波的吸收特性, 當(dāng)介電常數(shù)（ε）和磁導(dǎo)率（μ）一定時, 吸收系數(shù)（β）主要與介質(zhì)電阻率（ρ）有關(guān)。一般地,電阻率（ρ）越高，吸收系數(shù)（β）就越小,即介質(zhì)的性狀愈好；反之，即介質(zhì)的性狀越差, 電阻率（ρ）越低，吸收系數(shù)（β）就越大[2]。

由此可見,強度高、堅硬完整、較純的灰?guī)r中, 地球物理特征常表現(xiàn)為高速、高阻、低吸收的特征，而當(dāng)巖層中出現(xiàn)斷裂破碎帶、溶洞、泥沙充填時，則表現(xiàn)為：波速降低、電阻率降低, 吸收系數(shù)增大的現(xiàn)象，與完整灰?guī)r間存在較大的地球物理差異。介質(zhì)吸收系數(shù)的大小表示著不同的巖性和同一巖性的結(jié)構(gòu)完整性。

圖1 孔間電磁波CT原理示意圖

3 野外數(shù)據(jù)采集和處理

在野外鉆孔施工完成后需要安放PVC套管保護鉆孔，以保證探測順利進行，避免巖溶發(fā)育段塌孔現(xiàn)象發(fā)生，造成損失?，F(xiàn)場使用湖南岳陽奧成科技有限公司生產(chǎn)的HX-JDT-02孔間電磁波透視儀，該設(shè)備工作頻率范圍寬、動態(tài)范圍大、功耗極低、主控機極輕便、操作簡單快捷、穩(wěn)定可靠等特點。

3.1 探測實驗

探測之前，先進行儀器參數(shù)測試，首先在空氣中，由發(fā)射探頭發(fā)射，接受探頭接收，隨著距離的增加場強衰減如表1。從表中可以看出在空氣中，頻率為5MHz時，場強衰減最慢。

其次，選取鉆孔間距最小的ZK07和ZK08（間距23m）鉆孔進行土層吸收系數(shù)測試，取得接收機在無發(fā)射和有發(fā)射時電磁波場強值數(shù)據(jù)，結(jié)果見表2。從表中可以看出：在土層中，由于吸收系數(shù)過大，接收機接收到的信號與發(fā)射機不發(fā)射時接受的信號差異太小，無法進行后期解釋。所以，在土層中無法應(yīng)用電磁波CT進行探測。

表1 空氣中不同頻率、不同距離場信號場強值一覽表

最后，選取鉆孔間距33m的 ZK04和ZK05鉆孔進行巖層吸收系數(shù)實驗測試，測試結(jié)果表明：5～8MHz在巖石完整孔段觀測到的場強多為-45～-80dB，9～12MHz的場強多為-33～-100dB；在較完整灰?guī)r的孔段，將測試到的場強可視為正常場強值，用“兩孔測定法”確定初始場強E0′，并由此來計算較完整灰?guī)r吸收系數(shù)，經(jīng)式（2）。

計算得到的視吸收系數(shù)βs多在0.01～0.08db/m，該數(shù)值用于成果解釋中對異常進行劃分。

在以上詳細實驗的基礎(chǔ)上，確定了工作頻率、掃描間隔、點距、空氣中的背景場值和較完整灰?guī)r的吸收系數(shù)等探測參數(shù)，見表3。

表3 測試實驗結(jié)果參數(shù)一覽表

3.2 數(shù)據(jù)采集及處理

現(xiàn)場采用定點觀測法進行數(shù)據(jù)采集，將發(fā)射機和接收機分別放置在兩個鉆孔之中，把發(fā)射機固定在鉆孔某一深度上不動而移動接收機進行，完成一組數(shù)據(jù)采集后，移動發(fā)射機點位，采集下一組數(shù)據(jù)，如此循環(huán)，直至所有數(shù)據(jù)采集完成；如果將接收機固定在鉆孔某一深度上不動而移動發(fā)射機，其數(shù)據(jù)采集結(jié)果與前者相同。

電磁波 CT 資料的處理采用了阻尼代數(shù)重建法(ART)，其基本思想是先給被重建區(qū)域一個初始值，然后將所得到的投影殘差逐個沿其射線方向均勻地反投影回去，并不斷地對重建圖像進行修正，直到滿足要求為止。該層析處理方法采用加入阻尼技術(shù)來控制約束其結(jié)果，以克服迭代時容易發(fā)散的缺點，該方法占用的內(nèi)存少且計算速度快。

首先對原始數(shù)據(jù)進行平滑濾波及人機聯(lián)作刪除壞值，然后選用1m x 1m有限單元，使用 “井下電磁波分析系統(tǒng)”，對預(yù)處理后的數(shù)據(jù)進行反演計算，生成一對鉆孔的電磁波視吸收系數(shù)βs等值線剖面圖。為了克服迭加時常出現(xiàn)的發(fā)散現(xiàn)象，采用經(jīng)試驗選定的阻尼參數(shù)控制。

表2 土層中發(fā)射機有無信號時場強值（無發(fā)射/有發(fā)射）

4 CT剖面特征及地質(zhì)解譯

從電磁波CT圖可看到視吸收系數(shù)βs等值線值，等值線間區(qū)域顏色由色譜圖示方法加以表示。見圖2。數(shù)值越小，介質(zhì)對電磁波的吸收愈小，介質(zhì)的性狀越好；值越大，介質(zhì)對電磁波的吸收愈大，介質(zhì)的性狀愈差。依據(jù)電磁波CT視吸收系數(shù)βs等值線剖面圖，結(jié)合測區(qū)地質(zhì)、鉆探資料，分3步進行定性定量分析。

圖2 ZK04-ZK05孔間電磁波CT視衰減系數(shù)βs等值線剖面圖

⑴剔除假異常

βs等值線剖面圖的上部和底部，由于射線較稀，易出現(xiàn)假異常，必須剔除這些假異常，剩下的為有效異常。

⑵劃分正常區(qū)和異常區(qū)

正常場的確定：依據(jù)各鉆孔對的探測數(shù)據(jù)，選取較完整灰?guī)r孔段的曲線中間較穩(wěn)定的部分，來確定視吸收系數(shù)βs的正常場值，該值具備連續(xù)呈連片狀分布的特征，最終確定βs＜0.1db/m為正常區(qū)。高于該值的其它區(qū)域劃分為異常區(qū)。

⑶ 異常分析

依據(jù)有效異常的視吸收系數(shù)βs值大小和其等值線的形態(tài)、分布特征，結(jié)合測區(qū)地質(zhì)資料，對ZK04-ZK05孔間電磁波CT視衰減系數(shù)βs等值線剖面圖進行定性定量分析。

鉆孔ZK04數(shù)據(jù)從-16m開始，鉆孔ZK05從-22m開始。地下水位在ZK04鉆孔為-26m、ZK05鉆孔為-33m。分析結(jié)果如下：

① 正常區(qū)

剖面圖像上紫色－淺蘭色區(qū)域βs≤0.1db/m，呈連片分布，為完整灰?guī)r的反映。ZK04鉆孔資料表明在βs≤0.1db/m孔段，巖體完整性總體較好。

② 弱異常區(qū)

剖面圖上藍色－黃色區(qū)域βs=0.1～0.3db/m，異常形態(tài)多呈條帶狀或串珠狀。鉆孔資料顯示異常區(qū)范圍內(nèi)巖體完整性較差，有節(jié)理、裂隙分布，部分區(qū)域巖體較為破碎。

③ 強異常區(qū)

在地下水位面以上的溶洞，可能有兩種情況：

a)溶洞中填充物為泥、土，由于電磁波在泥、土中的吸收很強，所以，溶洞表現(xiàn)形式為βs≥0.3 db/m的紅色區(qū)域。

b)溶洞中沒有填充物，此時，溶洞中只有空氣，由于空氣對電磁波的吸收很弱，此種情況下，溶洞的表現(xiàn)形式為βs≤0.06db/m的紫色~淺藍色區(qū)域。

在地下水面以下的溶洞，由于洞中的充填物為水、泥土或它們的混合物，在此種情況下，溶洞的視吸收系數(shù)處于βs≥0.2db/m的綠色-紅色區(qū)域。

依據(jù)上述異常判定原則，在剖面圖中共圈定異常10處，詳細情況見表4。

在以上異常解釋資料的基礎(chǔ)上，繪制孔間地質(zhì)剖面圖（圖3）。溶洞呈斜向順層理發(fā)育，這個特點符合巖溶發(fā)育的特征，也證實了異常解釋的準(zhǔn)確性。

圖3 ZK04-ZK05孔間地質(zhì)剖面圖

表4 ZK04-ZK05孔間電磁波CT剖面異常一覽表

5 結(jié)論

孔間電磁波CT技術(shù)在普安店巖溶勘查中效果是良好的，能夠準(zhǔn)確反映出孔間巖溶的發(fā)育情況，可以圈定完整灰?guī)r、較破碎灰?guī)r、極破碎灰?guī)r和溶洞的界限，為后期的地基穩(wěn)定性評價和工程治理提供準(zhǔn)確數(shù)據(jù)和科學(xué)依據(jù)。

在工作中應(yīng)注意幾點：

（1）在結(jié)合地質(zhì)資料、地面物探成果、鉆探資料時要合理采用有效數(shù)據(jù)，摒棄不良數(shù)據(jù)引起的異常。

（2）異常解釋時要注意溶洞的充填情況，充填空氣、碎石土、含水碎石土、水等不同介質(zhì)，其引起異常的特征也會有所變化。

（3）孔距較遠時，接收機無法接收到發(fā)射信號，所以鉆孔布置間距不宜太大。

此外，在以后的工程處理工作中，電磁波CT技術(shù)在還可以監(jiān)測和檢測注漿效果，保證施工的質(zhì)量。

[1]吳 巖，等. 電磁波CT在碳酸鹽巖縫洞勘察中的應(yīng)用[J].工程地球物理學(xué)報，2009. 4.

[2]羅傳華，等.井間電磁波技術(shù)在居民區(qū)地下洞穴調(diào)查中的應(yīng)用[J]工程地球物理學(xué)報，2010.2.

[3]蔡加興.電磁波CT技術(shù)在水庫堤壩隱患探測中的應(yīng)用[J].人民長江，2001.33.

[4]吳以仁，邢鳳桐，易永森，等.鉆孔電磁波法[M].北京：地質(zhì)出版社，1982.

[5]《工程物探手冊》中國水利電力出版社，2011.

Application of Cross-hole Electromagnetic CT to Beijing Puandian Karst Prospecting

LI Yuanqiang
(Beijing Institute of Geology，Beijing 100120)

Cross-hole electromagnetic CT technology is a normal geophysical prospect method in engineering. The method is basis on tomographic technology. The result will gain a simple and visual map，which fully displays the cross-hole electromagnetic signal attenuation. The result shows different rocks and soil areas clearly. The karst caves are developed under the ground in Beijing Puandian. Many explored drills have been f nished, after the cross-hole electromagnetic CT scanning, the karst caves were located exactly, which data may fill the blanks between the drills, and infer the karst cave fillings. The technology gains true position data and other information of the underground karst cave. For the later stability evaluation and engineering treatment, it takes basal data and scientif c foundation.

Cross-hole electromagnetic CT；Karst prospecting；Data collecting and analyzing

P631.3+25

A

1007-1903（2014）03-0047-06