綜合物探方法在地鐵隧道過山段勘察中的應用

摘"要：隨著城市軌道交通的發展，城市山嶺隧道勘察已成巖土工程勘察中熱門課題之一。文章以南沙至珠海（中山）城際工程（萬頃沙—興中段）地鐵工程過山段勘察為例，選用微動探測法、高密度電法和土壤氡濃度測量法進行綜合探測，劃分了該區域巖土分界面，推斷出2條隱伏斷裂的位置及影響范圍。探測結果表明，微動探測法在山嶺地區干擾小的條件下，能夠采集到高質量的數據，劃分巖土分界面和推斷斷裂的效果較好；高密度電法對于斷裂破碎帶較為敏感，但對于地層結構劃分的能力較弱；土壤氡濃度測量能大致圈定斷裂的位置，但無法確定斷裂的影響深度和傾向。

關鍵詞：山嶺隧道；微動探測法；高密度電法；土壤氡濃度測量

Application of integrated geophysical prospecting methods in investigation of metro tunnel mountain-crossing sections

SUN Siliang

（Guangzhou Metro Design amp; Research Institute Co.， Ltd.， Guangzhou 510010， Guangdong， China）

Abstract： With the development of urban rail transit， the survey of urban mountain tunnels has become one of the hot topics in geotechnical engineering survey. This paper takes the survey of the mountain crossing section of Nansha to Zhuhai intercity subway project as an example. Three methods are selected for comprehensive detection including microtremor survey method， high-density electrical method， and soil radon measurement method. The rock soil interface in the area is divided and the locations and impact ranges of two hidden faults are inferred with good detection results. The research indicates that the microtremor survey method can obtain high-quality data in mountainous areas with minimal interference， and is effective in dividing rock soil interfaces and inferring faults. The high-density electrical method is more sensitive to fracture zones， but less ideal for classifying geological structures. The soil radon concentration measurement can roughly delineate the location of faults， but cannot determine the fault tendency and depth of influence.

Keywords： mountain tunnel; microtremor survey method; high-density electrical method; soil radon measurement

目前城市軌道交通工程逐步由城市中心走向近郊區，隨著城市山嶺隧道的不斷增多，巖土工程勘察越來越受到重點關注。在地質環境復雜、地形起伏劇烈和隧道埋深大的區域，傳統單一的物探勘察方法各有優劣，無法取得理想的探測效果，而綜合物探技術優勢互補，可快速、高效獲取多源、多尺度的地下空間數據（Mohamed Metwaly et al.，2014；Martin et al.，2015；賴劉保等，2016；胡程亮，2022；王軍成等，2023；徐琳等，2023；劉衡秋，2016），已被廣泛地應用到山嶺隧道勘察中。閆強剛等（2019）采用可控源電磁法、淺層折射波法、高密度電法和地質雷達法對山嶺隧道進行工程勘察，為設計、施工的順利開展奠定基礎；李軍等（2009）利用高密度電法與高頻大地電磁法對湘桂鐵路石山嶺隧道進行探測研究，查明了斷層構造的發育情況；袁偉等（2022）結合音頻大地電磁法、綜合測井法和高密度電法對鐵山坪隧道進行勘察，探測效果較好。每種物探方法都有其特定的適用范圍，綜合物探方法可以相互驗證，減少物探多解性，提高探測的準確性，可以有效解決工程勘察中遇到的地質難題。

本文以南沙至珠海（中山）城際工程（萬頃沙—興中段）里程段為WAK47+000 m—WAK47+540 m（明角山）的過山段巖土工程勘察為例，根據場地的環境條件和地球物理特征，選用微動勘探法、高密度電法、土壤氡濃度測量法3種物探方法對該區段隧道沿線進行勘察，綜合3種物探結果來確定覆蓋層的厚度、基巖面的起伏形態和埋深，圈定斷裂的位置和影響范圍，以驗證綜合3種物探方法在城市地鐵隧道過山段勘察中的有效性和適用性。

1 "測區地質及地球物理特征

測區位于南沙至珠海（中山）城際工程（萬頃沙—興中段）里程段為WAK47+000 m——WAK47+540 m（明角山）的過山段，長度約為540 m，地面高程為5～72 m，該段設計隧道埋深為18～85 m。測區地貌單元為低山，地形起伏相對較小，局部為山間溝谷地貌單元，植被茂盛。

探測范圍在大地構造上位于華南褶皺系（一級構造單元），粵北、粵東北-粵中坳陷帶（二級構造單元），粵中坳陷（三級構造單元），增城-臺山隆斷束（四級構造單元）的西南段，該區域整體上發育有3組斷裂，走向分別為NE向、近EW向和NW向，這些斷裂多形成于白堊紀的燕山運動，規模較大，見圖1。距本次測區最近的是五桂山北斷裂，該斷裂位于中山市張家邊、石岐、大涌西延至新會市睦洲，再往西可能伸展到臺山，長度超過50 km。

根據廣東省地質調查院（2000）《1∶25萬中華人民共和國地質圖（江門市幅）》可知，測區第四系主要由人工填土層和殘積層組成，基巖為早奧陶紀侵入的片麻狀細粒、細粒斑狀（含斑）黑云母二長花崗巖。根據收集到的資料，結合以往工作經驗，歸納得到測區巖土層的物理特性參數（表1）。

在地質解釋中，通常將中風化巖層上界面定義為巖土分界面。從表1可明顯得到覆蓋土層與下伏基巖之間和斷裂破碎巖與中風化基巖之間存在明顯的橫波波速和電阻率差異。因此，根據這些物性差異可區分巖土層分界面和圈定斷裂破碎帶。

2 "方法原理

2.1 "微動探測法

微動探測法是以平穩隨機過程理論為基礎（Aki，1957），采用數據處理技術從地表的微動信號中提取瑞雷波相速度頻散曲線（楊振濤等，2019；Capon，1970；Ling et al.，1993），經過反演以獲取地下介質的橫波速度結構的一種地球物理探測方法。

獲取面波頻散曲線后，用式（1）計算視S波速度V_x（Xu Peifen et al.，"2013），可將相速度頻散曲線（V_"r"－f曲線）轉換成視S波速度V_x隨深度的變化曲線（V_x－H曲線），再通過插值、光滑計算，最終得到視S波速度彩色剖面。

（1）

式中，V_r為瑞雷波速度，t _i為周期。

視S波速度V_x為既不同于相速度（V_r）也不同于S波速度（V_s）的面波物性參數，具有速度量綱。由于避免了反演過程中設置初始模型、選取反演結果等人為因素的影響，微動剖面結果能更客觀、直觀地反映地層巖性及構造變化。

2.2 "高密度電法

高密度電法以巖、土導電性的差異為理論基礎，地下電場的分布規律與地下巖土介質的電阻率分布密切相關，采用電極向地下供電形成人工電場，通過測量這種人工電場即可得到地下介質的視電阻率分布規律，并推斷出地層巖性與地質結構。

2.3 "土壤氡濃度測量

土壤氡濃度測量是指通過測量氡及其衰變子體產生的粒子的數量來圈定斷裂發育范圍的一種放射性方法（陳希泉等，2011；馮軍等，2011；李巧靈等，2024）。由于含有放射性元素的地下水從深部以地質構造為通道向地表運動，或者放射性氣體沿斷裂和裂隙向地表擴散遷移，從而在斷裂帶附近地表土壤中相對富集，斷裂破碎帶往往成為地下水或氣體運移的良好通道和儲存場所（劉偉等，2021）。放射性氣體主要為氡氣及其子體，都能逸散，其中以氡的半衰期最長，因此在地質構造及不同巖層的接觸帶上土壤中形成放射性（氡）的異常。斷裂帶中的氡氣向上遷移，常常會在上方土壤中形成線狀氡氣異常，通過這種氡氣異常推斷斷裂的位置及大致范圍。

3 "數據采集及處理

3.1 "測線的布置及采集參數的選取

結合測區地形地貌及地鐵隧道工程埋深等情況，本次物探工作主要采用微動法、土壤氡濃度測量和高密度電法進行探測，沿兩隧道內側中間線布置測線，微動法和土壤氡濃度測量法實施過程中測點布置相同，間距為20 m；高密度電法電極距為5 m。

1）微動探測法：采用深圳面元智能科技有限公司的SmartSolo IGU-16HR 3C型三分量智能地震傳感器進行微動數據采集。探測點距20 m，異常段適當加密。數據采集觀測方式采用直線型臺陣，臺陣半徑為30 m，臺陣中地震傳感器的個數為13個，間距為5 m，臺陣記錄時間為20 min。微動探測法數據采集現場工作照見圖2。

2）高密度電法：采用測控主機為中地裝（重慶）地質儀器有限公司生產的DUK-4高密度電法測量系統進行數據采集；采用溫納α排列和單極-偶極排列（AMN排列和MNB）進行探測，電極距為5 m，93道電極測量，供電時間0.7 s，斷電時間0.3 s，記錄2個波形數據，單極-偶極裝置中MN的間隔為2。高密度電法數據采集現場工作照見圖3。

3）土壤氡濃度測量法：采用HS-01測氡儀進行土壤氡濃度測量。探測點距20 m，異常段適當加密。土壤氡濃度測量法數據采集現場工作照見圖4。

3.2 "數據處理

1）微動數據處理方法：首先將原始數據按記錄時間進行切割，格式轉換成SG2，并導入驕佳面波軟件Surface plus；然后設置合適的參數計算頻譜，提取頻散曲線，將每個測點的頻散曲線通過公式（1）轉化為視V_s－深度關系；最后利用Surfer軟件將上述斷面數據轉換成Surfer網格化文件，繪制成視橫波速度等值線剖面圖，進行成果解釋。數據處理流程見圖5。

2）高密度電法數據處理方法：野外采集完數據后進行數據傳輸，將原始數據轉換成Surfer格式和RES2DINV格式。利用Surfer制作探測斷面的電阻率擬斷面圖，并進行初步解譯，然后對RES2DINV格式數據進行處理，主要有剔除畸變點、平滑處理、地形校正等，再進行反演，得到視電阻率成果剖面圖。

3）氡氣數據處理方法：對采集到的氡氣數據進行統計分析，將氡氣異常下限值設為測區測量背景值的1.5倍，測區氡氣背景值為去掉樣本中的一個最大值和一個最小值之后得到的平均值。

4 "結果與分析

1）微動法探測結果

圖6為二維微動法探測結果圖（微動測線里程段為WAK47+010 m—WAK47+530 m）。根據測區地層視橫波速度特征，可將視橫波速度結構剖面劃分為2個地質界面：覆蓋層（土層及強全風化花崗巖地層）和基巖層（中微風化花崗巖層）。根據表1測區巖土層物理特性參數表，將視橫波速度600 m·s^-1的界面定義為測區巖土分界面，視橫波速度在200～600 m·s^-1的地層為土層及強全風化花崗巖，視橫波速度在600～2 000 m·s^-1的地層為中微風化花崗巖。從微動成果剖面圖可以發現：①里程段WAK47+010 m—WAK47+444 m覆蓋層的平均厚度約為16 m，橫向分布均勻；里程段WAK47+444 m—WAK47+520 m，由于中微風化巖石出露地面，視橫波速度整體表現為高速，無覆蓋層。②里程段WAK47+189 m—WAK47+248 m（微動點WD010—WD013）存在1個向下貫穿低速異常區（編號1-1），與前后兩側的視橫波速度存在明顯的差異，速度值為430～780"m·s^-1，影響深度約65 m，影響深度和范圍較廣，設計隧道穿過該異常區。③里程段WAK47+270 m—WAK47+351 m（微動點WD015—WD018）存在1個低速異常帶（編號1-2），視橫波速度值為500～800 m·s^-1，影響深度約為18 m，影響范圍較小。初步推斷上述2個異常區為裂隙發育帶，巖石較為破碎。

2）土壤氡濃度探測結果

表2為明角山土壤氡濃度測量統計結果，該場地土壤氡濃度背景值為6 000 Bq·m^-3，異常閾值為8 000 Bq·m^-3，將氡氣濃度值大于8 000 Bq·m^-3的區域定義為氡氣異常段。

從圖7可以看到，測區內存在2個氡氣異常區域為：WAK47+211 m—WAK47+371 m（微動點在WD011—WD014之間）和WAK47+312 m—WAK47+371 m（微動點在WD016—WD019之間）。氡氣異常峰值之間存在1"個低于異常閾值的極小值段，推斷可能是由于復雜的地表環境和斷層破碎帶的非均勻性等因素造成的。氡氣異常段的分布范圍和微動法的低速異常范圍存在重疊，推斷氡氣異常段為斷裂的氡氣響應帶。

3）高密度電法探測結果

圖8為高密度電法探測結果圖（高密度電法測線里程段為WAK47+010 m—WAK47+425 m）。從圖8可見：①電阻率整體上表現為高阻，局部出現低阻異常。②電阻率在縱向上表現為自上而下先減小再增大，分層性較差，難以進行地質分層，出現該現象的原因為：a）地表為殘＼坡積土、全風化花崗巖，地下水位較深，所以淺地表電阻率值較高；b）受地下水位及強風化花崗巖地層影響，電阻率逐漸降低；c）在向下進入中—微風化花崗巖，電阻率逐漸增大。③從橫向上看，低阻異常區主要分布在50～100 m（里程段：WAK47+055 m—WAK47+105 m）、150～200 m（里程段：WAK47+155—WAK47+205）、230～270 m（里程段：WAK47+235 m—WAK47+275 m）、330～360 m（里程段：WAK47+335 m—WAK47+365 m），編號分別為2-1、2-2、2-3、2-4，異常區的電阻率值與周圍存在明顯差異。

4）綜合解譯結果

將微動法、土壤氡濃度測量和高密度電法探測結果進行綜合分析。在地層劃分方面，微動探測法可以劃分巖土分界面，確定覆蓋層的厚度和基巖面的起伏形態和埋深，而高密度電法難以劃分巖土分界面和確定覆蓋層的厚度，故地質分層以微動為主。在斷裂構造方面：低阻異常區（編號2-2、編號2-4）在平面位置上與微動成果低速異常區（編號1-1、編號2-2）和氡氣異常區位置重合，影響深度相符，推斷該區域受斷裂破碎帶影響，為斷裂帶，分別定名為F1和F2；低阻異常區（編號2-1、編號2-3）影響范圍和深度小，結合微動法和土壤氡濃度測量成果，推斷為裂隙發育帶。3種物探方法都能夠識別斷裂，確定斷裂的位置，但土壤氡濃度測量法無法確定斷裂的影響深度。

綜合解譯剖面圖如圖9所示。綜合解譯結果：測區覆蓋層平均厚度約16 m。區內分布2條斷裂，斷裂F1影響深度約為60 m，影響寬度約為50 m；斷裂F2影響深度約為40 m，影響寬度約為40 m；走向都為北東向，傾向南偏東，設計隧道穿過2個斷裂帶區域，可能對隧道產生影響。

依據上述解釋成果和前人的研究成果（李廣才等，2023；易隆科，2022），全面總結了3種物探方法的適用范圍及其優缺點，見表3。

5 "結論

1）采用微動勘探法、高密度電法、土壤氡濃度測量法3種物探方法對南沙至珠海（中山）城際工程（萬頃沙—興中段）地鐵工程過山段隧道沿線進行綜合探測。探測結果顯示，測區內100 m以淺的地層主要分為2層，即殘積土、坡積土及強全風化花崗巖的覆蓋層和中微風化花崗巖的基巖層，測區覆蓋層平均厚度約16 m。

2）本次物探成果推斷測區存在2處斷裂。F1：位于WAK47+175 m—WAK47+225 m里程段，斷裂走向約北東向，傾向為南偏東，影響寬度約50 m；F2：位于WAK47+324 m—WAK47+364 m里程段，斷裂走向約北東向，傾向為南偏東，影響寬度約40 m。

3）微動探測法在山嶺地區干擾小的條件下，能夠獲取得到較高質量的數據，劃分巖土分界面和推斷斷裂的效果較好；高密度電法對于斷裂破碎帶較為敏感，但對于地層結構劃分的能力較弱；土壤氡濃度測量能大致確定斷裂的位置，但無法確定斷裂影響深度、傾向等。

參考文獻

陳希泉，陳頡，羅先熔，歐陽菲，文美蘭2011.地氣（氡氣）測量方法尋找隱伏含礦斷裂試驗［J］.物探與化探，35（6）：817-820.

馮軍，李紅光，吳濤，玄月，李偉華，2011.河北隆堯隱伏斷裂地球化學探測［J］.物探與化探，35（5）：597-599.

廣東省地質調查研究院，2000．中華人民共和國區域地質調查報告廣州市幅（1：250 000）"［R］.

胡程亮，2022.綜合物探方法在常德沅江過江隧道工程勘察中的應用［J］.工程建設與設計（6）：177-179.

賴劉保，陳昌彥，張輝，王鐵領，張健，2016.綜合物探技術在長大公路隧道勘察中的應用研究［J］.工程勘察，44（12）：70-74.

李廣才，李培，姜春香，張鵬輝，王興宇2023.我國城市地球物理勘探方法應用進展［J］.地球物理學進展，38（4）：1 799-1 814.

李軍，薛軍平，2009.綜合物探在湘桂鐵路石山嶺隧道勘察中的應用［J］.工程地球物理學報，6（S1）：130-133.

李開朗，劉紅兵，2023.廣州市軌道交通穿越房屋密集區物探方法的選擇與應用［J］. 城市地質，"18 （2）："225-231.

李巧靈，雷曉東，李晨，雷坤超，韓宇達，趙旭辰，羅勇，2024.高密度電法和氡氣測量在地裂縫調查中的應用分析［J］.城市地質，19（1）：1-8 .

劉衡秋，2016.山嶺隧道工程綜合地質勘察方法探討［J］.城市地質，11（4）：34-38

劉偉，黃韜，王庭勇，劉怡，張繼，劉文濤，張琦斌，李強，2021.綜合物探方法在城市隱伏斷裂探測中的應用［J］.物探與化探，45（4）：1 077-1 087.

王軍成，"趙振國，"高士銀，"羅傳根，"李琳，"徐明鉆，"李勇，"袁國境，"2023. 綜合物探方法在濱海縣月亮灣地熱資源勘查中的應用［J］. 物探與化探，"47 （2）："321-330.

徐琳，陳永東，楊榮，官云彬，李海，龐仁俊，熊昌利，秦宇龍，2023.四川綿陽科技城新區玉皇溝斷層探測與活動性評價［J］.城市地質，18（2）：239-247

閆強剛，何壽迎，何松，2019.綜合物探技術在城市山嶺隧道勘察中的應用研究［J］.城市勘測（6）：204-208.

楊振濤，陳曉非，潘磊，王建楠，徐劍俠，張大洲，2019.基于矢量波數變換法（VWTM）的多道Rayleigh波分析方法［J］.地球物理學報，62（1）：298-305.

姚衍桃，詹文歡，劉再峰，張志強，詹美珍，2008.珠江三角洲的新構造運動及其與三角洲演化的關系［J］.華南地震（1）：29-40.

易隆科，2022.交叉梯度聯合反演等多物探方法在城市地質調查中的應用［J］.城市地質，17（2）：200-209.

袁偉，朱永珠，幸大軍，鄒維，2022.綜合地球物理勘察技術在鐵山坪隧道勘察中的應用［J］.工程地球物理學報，19（3）：334-340.

AKI K，"1957. Space and time spectra of stationary stochastic waves，"with special reference to microtremors ［J］. Bulletin，"Earthquake Research Institute，"35："415-456.

CAPON J，1970. Applications of detection and estimation theory to large array seismology ［J］. Proceedings of the IEEЕ，"58（5）：760-770.

LING S，"OKADA H，"1993. An extended use of the spatial autocorrelation method for the estimation of geological structure using microtremors ［C］. Proceedings of the 89th SEGJ Conference，"44-48.

MARTIN B，"ERWAN G，"2015. Comprehensive geophysical data integration and stratigraphic contacts delineation in a regional hydrogeological characterization study ［J］. Journal of Environmental amp; Engineering Geophysics，"20（2）：183-193.

Mohamed Metwaly，"Eslam Elawadi，"Sayed S R Moustafa，"Nasser Al-Arifi，"2014. Combined inversion of electrical resistivity and transient electromagnetic soundings for mapping groundwater contamination plumes in Al Quwy’"yia area，"Saudi Arabia ［J］. Journal of Environmental and Engineering Geophysics，"19（1）：45-52.

XU Peifen，"LING Suqun，"RAN Weiyan，"LIU Qingxiao，"LIU Jinguang，"2013. Estimating Cenozoic thickness in the Beijing Plain Area using array microtremor data ［J］. Seismological Research Letters，"84（6）："1 039-1 047.