趙東東，陳基煒，宗全兵，張寶松，邸兵葉，朱紅兵，王佳龍

(1.中國(guó)地質(zhì)調(diào)查局 南京地質(zhì)調(diào)查中心，江蘇 南京 210016；2.福州地鐵集團(tuán)有限公司，福建 福州 350000)

1 引 言

花崗巖風(fēng)化層中的球狀風(fēng)化核，俗稱“孤石”，主要是不易風(fēng)化的石英礦物形成石英角礫質(zhì)殘留核組成，在我國(guó)南方沿海地區(qū)廣泛發(fā)育燕山期花崗巖地區(qū)普遍存在[1,2]。未探明的孤石會(huì)給地鐵盾構(gòu)施工帶來(lái)重大安全隱患,在花崗巖殘積層中鉆遇孤石時(shí)，盾構(gòu)掘進(jìn)非常困難，盾構(gòu)機(jī)姿態(tài)難以控制，刀盤頻繁被卡或嚴(yán)重變形甚至磨損，即使能通過(guò)地面土壤加固、排石或換刀等技術(shù)措施處理，也會(huì)極大地增加施工成本，對(duì)工期和投資控制產(chǎn)生重大不利影響,更嚴(yán)重的情況下時(shí)，甚至導(dǎo)致工作面噴涌、塌方，危及地面行車或建筑物安全。因此，必須針對(duì)隧道盾構(gòu)施工中的孤石進(jìn)行精準(zhǔn)探測(cè)，同時(shí)采取合理的措施對(duì)孤石進(jìn)行破碎處理，以確保盾構(gòu)施工作業(yè)的順利進(jìn)行。

孤石發(fā)育的主要影響因素有花崗巖的礦物組成、結(jié)構(gòu)、構(gòu)造、巖體節(jié)理發(fā)育情況、溫度、地形、水文條件等，導(dǎo)致其具有埋藏深度不定、粒徑不均等特征，但在空間上也具有以下分布規(guī)律：①孤石一般分布于花崗巖隆起區(qū)兩側(cè)的全-強(qiáng)風(fēng)化帶中，越靠近花崗巖隆起中心位置，粒徑越大，數(shù)量越多。②隨著與花崗巖風(fēng)化面的距離變小，孤石分布密度減小但體積增大，即有“上多下少，上小下大”的分布規(guī)律特性。③風(fēng)化程度增強(qiáng)時(shí)，孤石體積減小，數(shù)量會(huì)增多，故全風(fēng)化帶中一般存在體積較小的孤石，強(qiáng)風(fēng)化帶中一般存在體積較大的孤石。

目前，針對(duì)花崗巖風(fēng)化層中孤石探測(cè)主要采取鉆探和物探等方法。其中，鉆探可直觀地揭露地層，對(duì)孤石進(jìn)行采樣，是最為精準(zhǔn)的點(diǎn)位探測(cè)方法。然而，相對(duì)于地鐵工程詳勘階段30～50 m的鉆孔間距而言，孤石尺寸要小得多，故通過(guò)詳勘鉆探所揭露的孤石十分有限，即使通過(guò)加密鉆孔提高揭露孤石的機(jī)率，但受成本、工期、場(chǎng)地條件等限制，難以實(shí)現(xiàn)孤石的空間精準(zhǔn)定位[3]。由于孤石與圍巖的物性差異主要是電性差異和速度差異，故在實(shí)際應(yīng)用中開展孤石探測(cè)的物探方法主要有：常規(guī)瞬變電磁法[4]、高密度電法[5,6]、探地雷達(dá)法[7]、跨孔電阻率層析成像法[8-11]、淺層反射地震法[12,13]、地震波層析成像法[14]、多道瞬態(tài)面波法[15]、微動(dòng)探測(cè)法[16-22]等。但是，由于地鐵建設(shè)所處的特殊環(huán)境，常規(guī)的物探方法受狹小地面空間和復(fù)雜噪聲情況影響較大，會(huì)面臨以下問(wèn)題：①高密度電法受場(chǎng)地條件和接地條件的限制，難以施工；②傳統(tǒng)探地雷達(dá)法探測(cè)深度有限，且極易受干擾的影響；③淺層反射地震法抗干擾能力弱，分辨能力較差；④跨孔電阻率層析成像法和地震波層析成像法探測(cè)結(jié)果直觀，效果較好，但受孔距較小的制約從而增加了成本；⑤多道瞬態(tài)面波法在城市馬路上有噪聲干擾大的問(wèn)題，且傳感器與地面的耦合也不易克服；⑥微動(dòng)探測(cè)法抗干擾能力強(qiáng)，有效劃分地層結(jié)構(gòu)，對(duì)孤石體垂向定位效果明顯，但是對(duì)孤石體的橫向分辨能力較弱；⑦常規(guī)瞬變電磁法存在勘探盲區(qū)和收發(fā)互感等問(wèn)題，縱、橫向分辨率較低，體積效應(yīng)明顯。

為了消除瞬變電磁法(TEM)收發(fā)線圈互感，實(shí)現(xiàn)瞬變電磁早期純二次場(chǎng)響應(yīng)觀測(cè)，席振銖等(2016)[23]提出采用反向?qū)ε即旁?由兩個(gè)完全相同的平行共軸且相隔一定間距、負(fù)載電流等值反向的線圈組成)發(fā)射一次場(chǎng)，且在反向?qū)ε即旁吹恼虚g平面接收二次場(chǎng)信號(hào)的方法，簡(jiǎn)稱等值反磁通瞬變電磁法(opposing coils transient electromagnetic method,OCTEM)。該方法繼承了TEM觀測(cè)純二次場(chǎng)、不接地、近源測(cè)量以及具有多種耦合裝置等優(yōu)點(diǎn)；從物理上直接避免了發(fā)射線圈和接收線圈之間的互感現(xiàn)象，消除了一次場(chǎng)對(duì)接收線圈的影響，實(shí)現(xiàn)了純二次場(chǎng)觀測(cè)；采用反向?qū)ε即旁词挂淮螆?chǎng)能量主要聚焦在近地表某一小范圍內(nèi)，故其渦流向下擴(kuò)散的圓環(huán)截面減小，增強(qiáng)了渦流聚集，降低了體積效應(yīng)，從而有效提高了縱、橫向分辨率[24]；淺層信號(hào)強(qiáng)度提高3～4個(gè)數(shù)量級(jí)，大大提高了數(shù)據(jù)信噪比，增強(qiáng)了抗干擾能力[24，25],為開展城市地鐵盾構(gòu)區(qū)間孤石精準(zhǔn)探測(cè)任務(wù)提供了一種可靠的地球物理探測(cè)手段。本文主要將等值反磁通瞬變電磁法應(yīng)用于福州地鐵F1線孤石精準(zhǔn)探測(cè)試驗(yàn)中，研究該方法在地鐵盾構(gòu)區(qū)間孤石精準(zhǔn)探測(cè)的可行性和有效性，并通過(guò)福州地鐵F1線機(jī)場(chǎng)-大鶴段的等值反磁通瞬變電磁反演成像結(jié)果圈定孤石體的位置、大小和空間分布。

2 基本原理及技術(shù)特點(diǎn)

等值反磁通瞬變電磁法(OCTEM,Opposing Coils Transient Electromagnetic Method)是測(cè)量等值反磁通瞬態(tài)電磁場(chǎng)衰減擴(kuò)散的一種新的瞬變電磁法。該方法以相同兩組線圈通以反向電流時(shí)產(chǎn)生等值反向磁通的電磁場(chǎng)時(shí)空分布規(guī)律，采用上下平行共軸的兩組相同線圈為發(fā)射源，且在該雙線圈源合成的一次場(chǎng)零磁通平面上，測(cè)量對(duì)地中心耦合的純二次場(chǎng)，方法原理及裝置示意如圖1所示。

圖1 等值反磁通瞬變電磁方法原理示意圖(席振銖等,2014)Fig.1 Schematic Diagram of OCTEM (Xi et al., 2014)

基于以上基本原理，等值反磁通瞬變電磁法相對(duì)常規(guī)瞬變電磁法而言,主要在以下幾個(gè)方面有明顯優(yōu)勢(shì)：①零磁通面接收純二次場(chǎng)優(yōu)勢(shì)。相比于傳統(tǒng)的單線圈源瞬變電磁法，等值反磁通瞬變電磁法的接收線圈處于一次場(chǎng)零磁通平面上，當(dāng)發(fā)射關(guān)斷時(shí)，接收面的一次場(chǎng)磁通不變，不會(huì)產(chǎn)生一次場(chǎng)感應(yīng)電動(dòng)勢(shì)，可免受一次場(chǎng)干擾，接收到純二次場(chǎng)信號(hào)。②中心耦合優(yōu)勢(shì)。等值反磁通瞬變電磁法一般選擇半徑為0.5 m左右的雙發(fā)射線圈(簡(jiǎn)稱“微線圈”)。相同磁矩條件下，相對(duì)于常用的小線圈中心回線裝置，雙線圈正下方的淺地表磁場(chǎng)能量集中，初始渦流場(chǎng)明顯較強(qiáng)，中心耦合度越高；雙線圈源激勵(lì)的感應(yīng)電流密度分布范圍遠(yuǎn)小于單獨(dú)發(fā)射時(shí)的范圍，體積效應(yīng)較弱，橫向分辨率明顯增強(qiáng)。③高密度小間隔采樣優(yōu)勢(shì)。采用24 bit，625 ksps 采樣間隔1.6 μs，以多周期疊加采樣取平均值作為實(shí)測(cè)記錄信號(hào)。由于采樣間隔小且無(wú)時(shí)窗抽道，相對(duì)于傳統(tǒng)瞬變電磁系統(tǒng)的大間隔采樣和時(shí)窗抽道得到實(shí)測(cè)記錄信號(hào)，其數(shù)據(jù)準(zhǔn)確度更高，縱向分辨率更高。④其它優(yōu)勢(shì)。等值反磁通瞬變電磁的微線圈為收發(fā)一體天線，設(shè)備體積小，重量輕，施工方便；設(shè)備抗干擾能力強(qiáng)，能夠在城市復(fù)雜環(huán)境下作業(yè)；減弱瞬變電磁法早期過(guò)渡過(guò)程影響，分辨率高。

3 工程應(yīng)用案例

3.1 試驗(yàn)區(qū)地質(zhì)和物性特征

根據(jù)福州地鐵F1線工程詳勘報(bào)告[26]及工程鉆孔地質(zhì)剖面附圖等地質(zhì)資料，系統(tǒng)總結(jié)福州地鐵沿線自上而下地層、孤石波速和電阻率特征。從表1統(tǒng)計(jì)數(shù)據(jù)可以看出：①孤石電阻率是圍巖(殘積土層及全-強(qiáng)風(fēng)化花崗巖)電阻率的3～7倍，有比較明顯的電性差異，具備電磁法(高密度電法、等值反磁通瞬變電磁法和大深度探地雷達(dá)法等)進(jìn)行孤石探測(cè)的物性條件；②填砂、粉質(zhì)黏土、粉細(xì)砂可以劃為一個(gè)電性層，(泥質(zhì))粉細(xì)砂、(含泥)粉細(xì)砂、淤泥質(zhì)土、粉質(zhì)黏土、殘積砂質(zhì)黏性土、全風(fēng)化花崗巖可以劃為一個(gè)電性層，二者之間有1.5～9倍的電阻率差異，全風(fēng)化花崗巖與砂土狀強(qiáng)風(fēng)化花崗巖之間的電阻率差異近3倍，砂土狀強(qiáng)風(fēng)化花崗巖與中風(fēng)化花崗巖之間的電阻率差異近6倍，具備電磁法進(jìn)行地層巖性分層的物性條件。

表1 福州地鐵F1沿線機(jī)場(chǎng)-大鶴段鉆孔電阻率

3.2 野外數(shù)據(jù)采集

等值反磁通瞬變電磁野外數(shù)據(jù)使用的儀器設(shè)備為中南大學(xué)和湖南五維地質(zhì)科技有限公司聯(lián)合研發(fā)的HPTEM-18型高精度等值反磁通瞬變電磁儀(如圖2所示)，測(cè)量二次磁場(chǎng)的感應(yīng)電壓計(jì)算視電阻率。該儀器設(shè)備主要由四部分組成：①HPTEM-18儀器主機(jī)；②一體式發(fā)送接收線框；③專用傳輸纜；④操作筆記本一臺(tái)。本次試驗(yàn)采用收發(fā)一體反磁通天線定點(diǎn)測(cè)量方式，數(shù)據(jù)采集之前調(diào)節(jié)接收天線至零磁通面，數(shù)據(jù)采集過(guò)程保持天線水平放置，操作筆記本與天線、主機(jī)及傳輸線保持5 m以上距離，盡量避免對(duì)數(shù)據(jù)造成干擾。本次測(cè)量工作的主要儀器采集參數(shù)為：發(fā)射基頻為25 Hz，發(fā)送電壓為14 V，最大發(fā)送電流為14 A，疊加周期為500次，測(cè)量重復(fù)次數(shù)為2次。

圖2 HPTEM-18工作示意圖Fig.2 Work diagram of HPTEM-18

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

采用HPTEM-18型儀器配套的專用軟件對(duì)野外采集的數(shù)據(jù)資料進(jìn)行抽道、地形校正、去噪濾波、一維反演及二維反演等處理。圖3為等值反磁通瞬變電磁數(shù)據(jù)處理流程示意圖，結(jié)合圖3對(duì)數(shù)據(jù)處理的幾個(gè)主要環(huán)節(jié)步驟進(jìn)行說(shuō)明：

圖3 等值反磁通瞬變電磁數(shù)據(jù)處理流程示意圖Fig.3 Data processing flow diagram of OCTEM

1)采用任意時(shí)間抽道方式有利于根據(jù)野外實(shí)測(cè)數(shù)據(jù)特點(diǎn)選取合適的抽道方式，便于原始數(shù)據(jù)編輯。

2)在數(shù)據(jù)預(yù)處理過(guò)程中，根據(jù)對(duì)稱四極法測(cè)量的地表電阻率和等值反磁通瞬變電磁測(cè)量初始值的耦合相關(guān)性，建立設(shè)置地表初始電阻率與初始幅值的關(guān)系，提高參與反演成像數(shù)據(jù)的可靠性[24]。

3)參考已知地質(zhì)、巖性和鉆孔資料建立二維反演成像的初始模型，進(jìn)一步提升二維反演成像的可靠性。

4)將上述步驟3)建立的初始模型代入二維反演成像模塊。其中，二維反演成像算法中數(shù)值模擬采用自適應(yīng)有限單元法，有利于在模型突變區(qū)域自動(dòng)加密網(wǎng)格，二維反演成像算法采用共軛梯度法，保證迭代穩(wěn)定收斂。

5)反演完成后，判斷反演結(jié)果是否滿足探測(cè)孤石所需達(dá)到的精度，若滿足則直接對(duì)反演結(jié)果進(jìn)行分析解釋，若不滿足則返回到步驟1)中，直到滿足探測(cè)孤石所需要的精度，再進(jìn)行分析解釋。

3.4 資料解釋

試驗(yàn)區(qū)位于福州市長(zhǎng)樂(lè)區(qū)葫蘆山北側(cè)的農(nóng)田荒地，周邊電磁干擾相對(duì)較少，地表主要上覆填砂，中間第四系巖性地層依次為粉質(zhì)黏土、粉細(xì)砂、(含泥)粉細(xì)砂、淤泥質(zhì)土、殘積砂質(zhì)黏性土等，下伏基巖依次為全風(fēng)化花崗巖、砂土狀強(qiáng)風(fēng)化花崗巖、碎塊狀強(qiáng)風(fēng)化花崗巖、中風(fēng)化花崗巖。本次工作主要目的是查明盾構(gòu)區(qū)間內(nèi)孤石發(fā)育位置，分布范圍，為保障地鐵盾構(gòu)安全施工提供基礎(chǔ)資料。試驗(yàn)區(qū)隸屬于福州地鐵F1線機(jī)場(chǎng)-大鶴段鉆MJKZ3-25-163～MJKZ3-25-172盾構(gòu)區(qū)間，雙向盾構(gòu)中心線完成6條等值反磁通瞬變電磁測(cè)線，點(diǎn)距1 m，線距3 m，測(cè)線長(zhǎng)度均為140 m。8-2線為穿過(guò)詳勘鉆孔MJKZ3-25-165～MJKZ3-25-169的測(cè)線，故下面以8-2線的反演電阻率斷面圖為例進(jìn)行分析。

圖4為8-2線反演電阻率斷面圖。從圖中可以看出：淺部電阻率相對(duì)較低，橫向連續(xù)性較好，垂向分層比較明顯，海拔高10 m以淺的相對(duì)低阻層為填土層和含泥粉細(xì)砂。深部電阻率很高，且橫向連續(xù)性較好，海拔-15 m以下的高電阻層為中風(fēng)化花崗巖，局部表現(xiàn)為基巖凸起。盾構(gòu)區(qū)間深度(海拔-15～5 m)表現(xiàn)為中-高阻，且橫向連續(xù)性相對(duì)較差，主要為殘積砂質(zhì)黏性土、全-強(qiáng)風(fēng)化花崗巖發(fā)育層，局部高阻異常為推測(cè)孤石體引起的。根據(jù)已知詳勘鉆孔資料可以看出MJKZ3-25-165～MJKZ3-25-169鉆遇孤石與反演電阻率斷面圖顯示的高阻異常吻合度較高，尤其橫向位置分辨率較高。但是，當(dāng)鉆孔位置正下方存在多層孤石，且孤石體在垂向上間隔較小時(shí)，在垂向上難以區(qū)分鉆遇的多個(gè)孤石體，僅表現(xiàn)為一個(gè)整體高阻凸起異常。8-2線累計(jì)推斷孤石體12個(gè)，與三個(gè)鉆孔鉆遇孤石體吻合較好的有7個(gè)，表明采用該方法探測(cè)孤石體可行、有效，且可靠性相對(duì)較高。

圖4 8-2線反演電阻率斷面Fig.4 Inversion resistivity profile of line 8-2

根據(jù)葫蘆山北側(cè)6條等值反磁通瞬變電磁測(cè)線的反演結(jié)果建立三維電阻率模型，如圖5所示。從圖5中可以看出明顯的電性界面和高電阻率異常，根據(jù)高電阻率異常和低電阻率異常等值線圈出比較明顯的孤石體和破碎富水帶。

圖5 葫蘆山北三維電阻率模型Fig.5 Three-dimensional resistivity model of the north of Calabash mountain

結(jié)合已知地質(zhì)、鉆孔資料和微動(dòng)橫波視速度模型建立了該地區(qū)的三維地質(zhì)結(jié)構(gòu)模型，從圖6可以看出：相對(duì)于二維地質(zhì)解譯剖面，三維地質(zhì)結(jié)構(gòu)模型刻畫了各地層的空間展布特征、雕刻出盾構(gòu)區(qū)間孤石的空間分布形態(tài)，表達(dá)方式更直觀、清晰。此段共計(jì)推斷孤石體39個(gè)，破碎帶8個(gè)，得到已知鉆孔驗(yàn)證的孤石體共計(jì)13個(gè)。

圖6 葫蘆山北三維地質(zhì)結(jié)構(gòu)模型Fig.6 Three-dimensional geological structure model of the north of Calabash Mountain

4 結(jié) 語(yǔ)

等值反磁通瞬變電磁法采用更小的發(fā)射和接收一體式線圈，設(shè)備體積小，重量輕，野外數(shù)據(jù)采集施工相對(duì)方便，改善了傳統(tǒng)瞬變電磁法的勘探效率和精度。該方法在地鐵盾構(gòu)區(qū)間孤石探測(cè)中的應(yīng)用試驗(yàn),表明其獲得的反演電阻率剖面分層明顯，淺部信息豐富，分辨率相對(duì)較高，推測(cè)孤石與鉆孔驗(yàn)證吻合度較好。但在資料綜合解釋中，需要充分結(jié)合已知鉆孔、地質(zhì)和其他物探方法用于克服瞬變電磁成果的多解性，提高地質(zhì)解譯精度。