王 薇，鄧小虎，金 聰，周紅偉，林 松,*

(1.中國(guó)地震局地震研究所,地震預(yù)警湖北省重點(diǎn)實(shí)驗(yàn)室，武漢 430071； 2.武漢地震工程研究院有限公司，武漢 430071)

地質(zhì)勘察工作是重大工程建設(shè)前期的重要基礎(chǔ)工作之一，灰?guī)r地區(qū)地質(zhì)勘察中，查明地下空間的巖溶發(fā)育情況顯得尤為重要。前期基礎(chǔ)地質(zhì)勘察中，巖溶發(fā)育規(guī)模、埋深、空間延伸狀態(tài)直接影響工程勘察進(jìn)度以及施工人員和設(shè)備安全，后期重大工程建設(shè)中及建成后，巖溶的存在直接關(guān)系到建筑物和人民生命財(cái)產(chǎn)安全。巖溶作為灰?guī)r地區(qū)最典型的不良地質(zhì)現(xiàn)象[1]，發(fā)育狀態(tài)千奇百怪，常規(guī)鉆探雖然直觀，但由于數(shù)量的局限性，不能完整提供地下地質(zhì)體巖溶發(fā)育具體規(guī)模以及分布情況，且成本高、噪聲污染大[2]。為快速、準(zhǔn)確查明巖溶具體位置、埋深以及發(fā)育規(guī)模，地球物理手段在眾多場(chǎng)地寬闊，受干擾較小地區(qū)的重大工程建設(shè)中發(fā)揮著重要作用：眾多學(xué)者利用高密度電法[3-5]、地震反射波法[6-8]以及探地雷達(dá)[9-11]進(jìn)行巖溶探測(cè)和巖溶路基注漿質(zhì)量檢測(cè)。然而，城市巖溶探測(cè)相比其他巖溶探測(cè)具有場(chǎng)地狹窄、干擾因素多。協(xié)調(diào)難度大等特點(diǎn)。高密度電法和地震反射波法測(cè)線無(wú)法展布，電極和檢波器的插拔亦會(huì)存在一定困難，對(duì)于城市路面會(huì)造成一定破壞；探地雷達(dá)則主要在探測(cè)深度方面受到局限。因此，重大工程巖溶探測(cè)必須根據(jù)場(chǎng)地條件選取合適的地球物理探測(cè)方法，隨著地球物理探測(cè)方法基礎(chǔ)理論研究方面的不斷深入以及勘探設(shè)備的更新，無(wú)論是外業(yè)數(shù)據(jù)采集還是室內(nèi)資料分析和處理，都取得了質(zhì)的飛躍[12]，井間電磁波CT技術(shù)基本不受場(chǎng)地范圍限制，尤其是在潛水面較低的地區(qū)，不像地震CT因?yàn)槿彼疅o(wú)法耦合而不能開展工作。因此，井間電磁波CT技術(shù)有其自身優(yōu)點(diǎn)，逐漸在城市重大工程巖溶探測(cè)中發(fā)揮出優(yōu)勢(shì)，進(jìn)而成為研究熱點(diǎn)，學(xué)者們?cè)谡莺头囱莘矫婢隽舜罅抗ぷ?，取得不少成果：陳湘華等[13]用數(shù)值模擬手段進(jìn)行電磁波正演模擬；孫茂銳[14]對(duì)電磁波CT二維與三維成像進(jìn)行研究；彭永良等[15]利用鉆孔電磁波 CT進(jìn)行復(fù)雜采空區(qū)注漿效果檢測(cè)；黃生根等[16]通過(guò)電磁波CT技術(shù)對(duì)溶洞進(jìn)行模擬分析與應(yīng)用。然而，通過(guò)電磁波CT技術(shù)進(jìn)行巖溶探測(cè)，分辨率卻有待進(jìn)一步提高。

現(xiàn)以某地鐵巖溶勘察為例，采用電磁波CT對(duì)灰?guī)r地區(qū)進(jìn)行巖溶精細(xì)探測(cè)，應(yīng)用聯(lián)合迭代重建技術(shù)(simultaneous iterative reconstruction technique, SIRT)進(jìn)行反演，獲取電磁波CT吸收系數(shù)等值線圖，并通過(guò)吸收系數(shù)判斷巖溶發(fā)育狀態(tài)。在城市重大工程巖溶勘察中，若采用電磁波CT技術(shù)可對(duì)地下空間巖溶發(fā)育狀況進(jìn)行較為精細(xì)的探測(cè)，較好查明巖溶發(fā)育規(guī)模與特征；則電磁波CT成果剖面可彌補(bǔ)單純依靠鉆探資料僅提供點(diǎn)信息不足的缺陷，與此同時(shí)，勘察速度得到提高，經(jīng)濟(jì)成本節(jié)約，具有較高的應(yīng)用價(jià)值和科學(xué)意義。

1 地質(zhì)背景與地球物理?xiàng)l件

1.1 地質(zhì)背景

新近紀(jì)以來(lái)，某地鐵線路經(jīng)過(guò)區(qū)內(nèi)地層結(jié)構(gòu)復(fù)雜，并有褶皺出現(xiàn)，多呈閉合狀態(tài)，且各地層厚度不均，區(qū)間內(nèi)第四系分布廣泛，巖性主要為河湖相、沖積相和坡積碎屑物，主要有雜填土、粉質(zhì)黏土夾粉土；下伏基巖主要為三疊系大冶組灰?guī)r(T1d)，成傾斜狀態(tài)，角度較陡，局部直立甚至出現(xiàn)倒轉(zhuǎn)[17-18]；中部為泥質(zhì)灰?guī)r夾黏土巖條帶、泥灰?guī)r夾灰?guī)r條帶、泥灰?guī)r、泥質(zhì)條帶灰?guī)r、灰?guī)r夾黏土巖(頁(yè)巖)；下部為頁(yè)巖及鈣質(zhì)頁(yè)巖、黏土巖夾灰?guī)r條紋[19]。

1.2 地球物理?xiàng)l件

不同的地下介質(zhì)對(duì)電磁波吸收程度不同，介質(zhì)的均勻性、礦物質(zhì)含量、裂隙發(fā)育程度、巖性以及含水量等因素都會(huì)導(dǎo)致電磁波的吸收差異。根據(jù)地質(zhì)背景條件、鉆孔資料揭露以及原位測(cè)試，某地鐵線路區(qū)間覆蓋層、下伏基巖以及基巖中的裂隙、溶蝕現(xiàn)象對(duì)電磁波的吸收系數(shù)存在較大差異(表1)，第四系覆蓋層吸收系數(shù)較大，完整基巖吸收系數(shù)較小，溶洞、溶蝕對(duì)電磁波的吸收系數(shù)與巖溶發(fā)育規(guī)模以及填充物質(zhì)有關(guān)，在碳酸鹽巖中，大部分溶洞為無(wú)填充或半填充，因此吸收系數(shù)范圍較為固定。基于上述條件，在某地鐵巖溶探測(cè)中開展電磁波CT具備良好的地球物理前提。

表1 不同介質(zhì)吸收系數(shù)Table 1 Absorption coefficient of different media

2 井間電磁波CT層析成像技術(shù)

2.1 基本原理

電磁波CT層析成像是基于電磁場(chǎng)理論和天線理論的一種工程物探方法，根本出發(fā)點(diǎn)是基于惠更斯原理，將對(duì)稱偶極天線分別置于2個(gè)相鄰鉆井中，如圖1所示，進(jìn)行同步后，通過(guò)一個(gè)鉆井中放置發(fā)射天線發(fā)射場(chǎng)強(qiáng)為E0的電磁波，另一鉆孔放置接收天線接收?qǐng)鰪?qiáng)為E的電磁波，二者關(guān)系為

(1)

式(1)中：β為介質(zhì)對(duì)電磁波的吸收系數(shù)；r為接受孔與發(fā)射孔間的距離；f(θ)為與接收和發(fā)射有關(guān)的方向因子；θ為接收點(diǎn)位置與發(fā)射電場(chǎng)方向的夾角[20-21]。

圖1 電磁波CT測(cè)試與成像原理圖Fig.1 Electromagnetic wave computed tomography test and imaging schematic diagram

數(shù)據(jù)采集時(shí)，先在某一深度固定發(fā)射天線，接收天線在另一孔中由上往下按照一定間距移動(dòng)，直至達(dá)到設(shè)計(jì)孔深，可得到若干條射線，每條射線記錄了該射線方向地質(zhì)體對(duì)電磁波的吸收狀況。當(dāng)測(cè)區(qū)中某點(diǎn)附近有數(shù)條射線通過(guò)時(shí)，便可以通過(guò)式(1)求得該點(diǎn)的物性參數(shù)。

數(shù)據(jù)處理時(shí)，可將相鄰兩孔的測(cè)量區(qū)間均勻劃分為小方格，以便于實(shí)現(xiàn)空間離散化。為方便計(jì)算，可設(shè)置m個(gè)橫向方格，n個(gè)縱向方格，方格總數(shù)為m×n，并對(duì)每個(gè)方格進(jìn)行編號(hào)1,2,…，k，方格大小可以根據(jù)分辨率要求及數(shù)據(jù)采集情況進(jìn)行約束調(diào)整。若Xj代表第j個(gè)方格地下介質(zhì)對(duì)電磁波的吸收系數(shù)，則對(duì)于任意一條射線i，到接收點(diǎn)時(shí)，其射線的電磁波損失總量Pi可表示為

Pi=A1X1+A2X2+…+AkXk

(2)

式(2)中：Ak為射線經(jīng)過(guò)第k個(gè)方格的長(zhǎng)度。因此，若在井間不同深度進(jìn)行N次數(shù)據(jù)發(fā)射與接收，則可列出如下矩陣：

(3)

利用矩陣方程可解出每個(gè)方格的電磁波吸收系數(shù)，降系數(shù)網(wǎng)格化后便可獲取電磁波吸收系數(shù)等值線圖[21]。

2.2 采集參數(shù)設(shè)置

圖2 不同頻率電磁波電場(chǎng)強(qiáng)度能量-深度變化曲線Fig.2 Electric field intensity energy depth curve of electromagnetic wave with different frequency

選用HX-JDT-02電磁波系統(tǒng)進(jìn)行數(shù)據(jù)采集，鉆孔間距為20～25 m，天線頻率選擇4 MHz-8 MHz-12 MHz三種頻率同時(shí)進(jìn)行掃描，采用2.1節(jié)所提及的定發(fā)方式進(jìn)行數(shù)據(jù)采集，激發(fā)與接收角度控制在45°以內(nèi)，發(fā)射步距為1 m，接收步距為0.5 m，探測(cè)深度根據(jù)孔深和PVC管深度共同決定。經(jīng)過(guò)測(cè)試與比較(圖2)，3種不同頻率的天線對(duì)于某地鐵巖溶探測(cè)電磁波衰減趨勢(shì)差別不大，選擇8 MHz進(jìn)行數(shù)據(jù)反演成圖。

2.3 數(shù)據(jù)處理

電磁波CT處理方法有代數(shù)重建、聯(lián)合迭代、反投影法、共軛梯度。最大熵法以及波前射線追蹤等[22]。根據(jù)不同的原始數(shù)據(jù)，數(shù)據(jù)處理方法也各有優(yōu)缺點(diǎn)，針對(duì)某地鐵數(shù)據(jù)采集過(guò)程中，高壓和民用電等干擾現(xiàn)象嚴(yán)重，為保證收斂，使得反演圖像準(zhǔn)確反映實(shí)際地質(zhì)情況，采用了聯(lián)合迭代重建技術(shù)(SIRT)進(jìn)行數(shù)據(jù)處理，處理流程按照?qǐng)D3所示。

圖3 數(shù)據(jù)處理流程Fig.3 Data processing flow

3 研究區(qū)探測(cè)成果及地質(zhì)解釋

在某地鐵軌道巖溶探測(cè)中，按照設(shè)計(jì)和規(guī)范要求，沿線路布設(shè)排鉆進(jìn)行電磁波CT探測(cè)，鉆孔間距20～25 m，相鄰兩鉆孔之間進(jìn)行一一配對(duì)，進(jìn)行井間電磁波CT掃描(圖4)，并按照上述采集參數(shù)和數(shù)據(jù)處理流程獲取CT吸收系數(shù)等值線剖面。選取了7個(gè)剖面進(jìn)行探討，剖面顯示的覆蓋層與灰?guī)r、灰?guī)r與巖溶、溶蝕裂隙區(qū)的吸收系數(shù)差異明顯，巖溶發(fā)育規(guī)模和特征較為清晰，具備較高的分辨率。地質(zhì)成果解釋以CT吸收系數(shù)分布特征為基礎(chǔ)，結(jié)合鉆孔資料對(duì)比解釋，最終形成地質(zhì)解釋剖面。

圖5為1號(hào)CT剖面、2號(hào)CT剖面、3號(hào)CT剖面、4號(hào)CT剖面聯(lián)合反演成像與地質(zhì)解釋圖。1號(hào)剖面以YMK1為發(fā)射孔，YMK2為接受孔；2號(hào)剖面、3號(hào)剖面以YMK3為發(fā)射孔，分別以YMK2和YMK4為接收孔；4號(hào)剖面以YMK5為發(fā)射孔，YMK4為接收孔。將原始數(shù)據(jù)整理后，4個(gè)剖面按照所提及的數(shù)據(jù)處理方法進(jìn)行聯(lián)合反演，從成果圖中可以明顯看出，基巖面與上覆地層界限明顯，巖溶、溶蝕與灰?guī)r之間的視吸收系數(shù)差異較大，覆蓋層視吸收系數(shù)大于0.45 Nper/m，在基巖面之下，巖溶視吸收系數(shù)大于0.72 Nper/m，溶蝕現(xiàn)象視吸收系數(shù)介于0.48～0.72 Nper/m。結(jié)合鉆孔資料對(duì)比CT剖面聯(lián)合反演成像：2號(hào)剖面與3號(hào)剖面之間存在1處規(guī)模較大的視吸收系數(shù)異常，深度范圍為32～38 m，巖溶橫向發(fā)育規(guī)模較大，寬度約18 m，形態(tài)呈不規(guī)則狀態(tài)，溶洞外圍呈不規(guī)則溶蝕狀態(tài)；4號(hào)剖面中存在1處形態(tài)不規(guī)則的視吸收系數(shù)異常，深度范圍為33～36 m，巖溶橫向規(guī)模約7 m。對(duì)比2號(hào)剖面與3號(hào)剖面，兩個(gè)剖面的視吸收系數(shù)異常在同一位置呈現(xiàn)高度吻合狀態(tài)，且縱向發(fā)育規(guī)模與特征相符。為驗(yàn)證異常是否存在，在3號(hào)剖面和4號(hào)剖面中進(jìn)行了鉆孔驗(yàn)證，驗(yàn)證孔分別為YZK1、YZK2。兩個(gè)驗(yàn)證孔揭露的覆蓋層厚度與視吸收系數(shù)所反映的覆蓋層深度一致；驗(yàn)證孔YZK1揭露3號(hào)剖面在深度33～35 m范圍內(nèi)，灰?guī)r呈破碎狀態(tài)，伴隨溶蝕現(xiàn)象，深度35～38.5 m范圍為半填充溶洞；14.3～22.3 m處存在巖溶。驗(yàn)證孔YZK2揭露4號(hào)剖面在深度32～34 m范圍內(nèi)，灰?guī)r呈破碎狀態(tài)，伴隨溶蝕現(xiàn)象，深度34～36 m范圍為半填充溶洞；通過(guò)對(duì)比發(fā)現(xiàn)，驗(yàn)證孔揭露的異常位置與CT是視收系數(shù)異常位置相符。

圖4 CT剖面及巖溶平面投影Fig.4 CT profile and karst plane projection

灰?guī)r地區(qū)地下巖溶發(fā)育狀態(tài)千奇百怪，為更加準(zhǔn)確揭露巖溶發(fā)育的空間形態(tài)，以YMK3號(hào)鉆孔為中發(fā)射孔，分別與鉆孔YRK3、YRK4、YLK3、YLK4配對(duì)進(jìn)行了相同的數(shù)據(jù)采集工作，形成圖4所示的其他剖面。圖6為聯(lián)合5號(hào)CT剖面數(shù)據(jù)與6號(hào)CT剖面數(shù)據(jù)獲取的反演成果圖。視吸收系數(shù)顯示，覆蓋層深度為26～28 m，5號(hào)剖面、6號(hào)剖面在同一深度位置32～36 m均出現(xiàn)視吸收系數(shù)異常，視吸收系數(shù)值為0.72 Nper/m。YMK3號(hào)鉆孔揭露，33～38 m范圍為半填充溶洞，由此可說(shuō)明，2、3號(hào)剖面揭示的溶洞在5、6號(hào)剖面中視吸收系數(shù)也同樣出現(xiàn)異常，圖5結(jié)合圖6，將使YMK3為中心的溶洞空間發(fā)育規(guī)模和形態(tài)更加明確。

圖5 YMK1-YMK2-YMK3-YMK4電磁波CT聯(lián)合反演成果與地質(zhì)解釋圖Fig.5 Joint inversion results and geological interpretation map of YMK1-YMK2-YMK3-YMK4 electromagnetic wave CT

圖7中出現(xiàn)兩處視吸收系數(shù)異常，第一處異常以YMK3為中心，在7號(hào)剖面與8號(hào)剖面中均有分布，異常深度分布在33～37 m，視吸收系數(shù)為0.76 Nper/m，高于3號(hào)剖面中驗(yàn)證孔YZK1揭露的半填充溶洞的視吸收系數(shù)(該處為0.72 Nper·m-1)，結(jié)合YMK3號(hào)鉆孔資料，可知以YMK3號(hào)為中心，深度32～38 m范圍內(nèi)的巖溶呈立體分布狀態(tài)，橫向最大規(guī)模為18 m。第二處異常在8號(hào)剖面中，31～37 m深度范圍內(nèi)，橫向?qū)挾?0～38 m，視吸收系數(shù)大于0.76 Nper/m，結(jié)合驗(yàn)證鉆孔YZK3，證實(shí)該處為半填充溶洞，填充物質(zhì)多為破碎灰?guī)r夾黏土。

圖6 YRK3-YMK3-YLK4電磁波CT聯(lián)合反演成果與地質(zhì)解釋圖Fig.6 Joint inversion results of YRK3-YMK3-YLK4 electromagnetic wave CT and geological interpretation map

圖7 YRK4-YMK3-YLK3電磁波CT聯(lián)合反演成果與地質(zhì)解釋圖Fig.7 Joint inversion results of YRK4-YMK3-YLK3 electromagnetic wave CT and geological interpretation map

4 結(jié)論

充分利用井間電磁波CT不受場(chǎng)地限制、分辨率較高等優(yōu)點(diǎn)對(duì)某地鐵進(jìn)行巖溶探測(cè)，采用SIRT進(jìn)行反演，獲取了電磁波CT吸收系數(shù)等值線圖，通過(guò)視吸收系數(shù)異常與鉆孔資料進(jìn)行對(duì)比解釋，查明地下空間巖溶發(fā)育狀態(tài)，取得良好效果。

研究實(shí)例中，以鉆孔YMK3為中心，整合三個(gè)方向的8個(gè)電磁波CT剖面，最終形成三個(gè)聯(lián)合剖面(圖5～圖7)，三個(gè)聯(lián)合剖面共揭露5處吸收系數(shù)異常，通過(guò)鉆孔資料發(fā)現(xiàn)，覆蓋層吸收系數(shù)大于0.45 Nper/m，溶蝕與半填充溶洞視吸收系數(shù)下界介于0.48～0.76 Nper/m，完整基巖視吸收系數(shù)介于0.35～0.45 Nper/m。異常除利用鉆孔進(jìn)行驗(yàn)證之外，三個(gè)聯(lián)合剖面呈梅花狀交叉，交叉位置處的異常也得到了相互印證。

隨著城市發(fā)展，地下空間的安全越來(lái)越引起人們的重視，重大工程巖溶探測(cè)僅依靠鉆探手段提供的點(diǎn)信息來(lái)推測(cè)巖溶發(fā)育狀態(tài)，甚至推進(jìn)后期工程建設(shè)將會(huì)導(dǎo)致很多隱患?；?guī)r地區(qū)上覆地層界線變化劇烈，巖溶發(fā)育情況復(fù)雜，不能簡(jiǎn)單地依靠工程地質(zhì)剖面來(lái)確定地下地質(zhì)情況。井間電磁波CT技術(shù)可充分利用優(yōu)勢(shì)，為巖溶探測(cè)提供二維面信息甚至是三維信息，以彌補(bǔ)鉆探的局限性，從而為重大工程建設(shè)提供地質(zhì)及地球物理參考依據(jù)，具有較好的科學(xué)意義和研究?jī)r(jià)值。