E-mail Alert Rss
 

物探与化探, 2024, 48(1): 281-286 doi: 10.11720/wtyht.2024.2602

工程勘察

基于三维电阻率法的水电工程隧道地质条件探查

黄瑶,

南宁理工学院 土木与工程学院, 广西 桂林 541006

Exploring geological conditions for tunnel construction in hydropower engineering using a 3D resistivity method

HUANG Yao,

Nanning College of Technology School of Civil Engineering, Guilin 541006, China

责任编辑: 叶佩, 沈效群

收稿日期: 2022-12-9   修回日期: 2023-03-27  

基金资助: 广西高校中青年教师科研基础能力提升项目(2019KY1052)

Received: 2022-12-9   Revised: 2023-03-27  

作者简介 About authors

黄瑶(1988-),女,副教授,主要从事结构工程方面的教学与研究工作。Email:huangyao12309@outlook.com

摘要

为开展水电工程隧道施工地质条件的探查与研究,基于三维电阻率技术构建了隧道地质条件的计算模型,采用数值模拟获得隧道模型的三维电性分布特征;将该技术应用于云南某水资源配置工程现场探查,经钻孔验证,获得了良好的探查效果。研究结果表明三维电阻率法对水电工程隧道的施工地质条件探查具有好的适用性,能够准确判断地层层厚,圈定岩溶空洞的大小、位置及充填物性质。通过探测研究获得的定量、定性信息,为隧道施工管理、信息化建设和防灾提供了可靠的依据。

关键词： 三维电阻率法; 隧道施工; 岩溶; 数值模拟

Abstract

To explore the geological conditions for the tunnel construction in hydropower engineering, this study built a calculation model for tunnel geological conditions using a 3D resistivity method. Through numerical simulations, this study determined the 3D resistivity distribution of the tunnel model. Then, the model was applied to the field exploration of a water resource allocation project in Yunnan, yielding satisfactory exploration results, as verified through drilling. The findings suggest that the 3D resistivity method can be effectively applied to the exploration of geological conditions for tunnel construction in hydropower engineering by accurately determining formation thicknesses, as well as the sizes, locations, and filling characteristics of karst cavities. The quantitative and qualitative data obtained from exploration in this study lay a reliable foundation for the management, informatization, and disaster prevention of tunnel construction.

Keywords： 3D resistivity method; tunnel construction; karst; numerical simulation

PDF (2923KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

黄瑶. 基于三维电阻率法的水电工程隧道地质条件探查[J]. 物探与化探, 2024, 48(1): 281-286 doi:10.11720/wtyht.2024.2602

HUANG Yao. Exploring geological conditions for tunnel construction in hydropower engineering using a 3D resistivity method[J]. Geophysical and Geochemical Exploration, 2024, 48(1): 281-286 doi:10.11720/wtyht.2024.2602

0 引言

随着工程隧道的建设量不断增加,所遇到的工程地质问题也越来越复杂。水电工程隧道由于用途的特殊性,隧道选址要求较高,隧道施工中对地质条件探查的要求也比较高。在隧道施工时,探测异常地质体,并确定其空间位置、形状、大小及埋深,查明有害地质体中填充物(如有害气体、地下水、沉积物或泥浆)的性质,探测不同岩性构造之间的分界面至关重要。为了有效获取隧道工作面前的地质构造,隧道地质条件探查的理论和技术一直是研究者们关注的焦点。近50年来,超前勘探从钻井和地质分析发展到非破坏性地球物理勘探,由于地球物理方法具有探测快速、破坏性较小以及施工高效等特点,在隧道地质条件的探查以及异常体的圈定方面得到了广泛应用[1⇓⇓⇓⇓-6]。由于电阻率法对隧道各地质结构的电性特征响应明显,对地层电导率差异敏感,识别地质结构及地层的能力较强,所以电阻率法在隧道地质条件探查方面发展迅速。采用高密度电阻率法可以准确探查出隧道前方溶洞的位置及判断出岩层的富水区[7-8],也可以探测出隧道隐伏断层的位置[9]。但上述电阻率法对隧道施工地质条件的探查与研究只是针对二维平面,对地质条件的空间分布特征研究较少,不能从整体空间上判断隧道地质条件。而水电工程隧道前方地质异常体的空间位置、形状、大小是预测隧道前方围岩质量,分析隧道施工过程中可能遇到的地质灾害的类型、风险和规模等亟待解决的关键难题,因此,对水电工程隧道的地质条件三维空间探查非常重要。

近年来,随着计算机技术的进步,三维电阻率法取得了长足的发展,在数值模拟方面[10⇓-12]和三维电阻率法反演方面[13⇓⇓-16]都取得了不少的成果。本文以三维电阻率法为基础,构建水电工程隧道中不同地质异常体的三维数值模型,模拟地质异常体的空间电性分布特征及响应特征,为解决水电工程隧道施工所遇到的地质问题提供技术支持。同时,将三维电阻率技术应用于现场实际探测中,结合地质资料,实现了对隧道施工地质条件的三维探查,保障了隧道的安全施工,证明该方法具有良好的推广实用价值。

1 三维电阻率技术

三维电阻率法勘探是电阻率法的三维可视化的发展,以岩层、地层及地质结构的电性差异为基础,通过空间范围内隧道及周围地质结构的电阻率值的不同来判断隧道相应的地质条件[17⇓-19]。电阻率法对水、岩层、空洞的电阻率差异响应明显,因此对水电隧道地质条件的探查具有良好的探测效果,三维电阻率勘探技术能够有效判断异常体的空间位置、形状及大小,完成对异常体的定位。电阻率法的探测深度与测线的长度和供电电流等因素有关,一般的水电工程隧道施工深度较浅,三维电阻率法的探测深度能满足浅层隧道的探测要求,对含水构造的响应特征明显,适用于水电工程隧道的建设。

本文采用的是用有限元法求解正演问题的全三维反演方法。正演计算构建三维地电断面模型,将三维区域划分为大量的小区块,每个小区块赋予相应的电阻率值,根据相应的插值公式,确定整个空间内的电阻率值[20⇓-22]。反演计算通过现场探测,获得测区内的电位数据。反演问题的一般形式可表达为:

Δd=GΔm

式中:G为Jacobi矩阵;Δd为观测数据d和正演理论值d0的残差向量;Δm为初始模型m的修改向量[23⇓⇓-26]。为了提高复杂解的稳定性,在传统的阻尼最小二乘法的基础上加入光滑约束,反演求得光滑模型,其中Δm的算法为:

(GTG+λCTC)Δm=GTΔd

式中:C是模型光滑矩阵;λ为阻尼因子。通过对Jacobi矩阵G及大型矩阵逆的计算,获得各三维网格相应的电性数据。通过比较正演计算与反演计算的结果,利用反演程序在反演过程中调节各子块的电阻率,以均方差[27-28]作为相似程度的判断标准,均方差小于标准后,正演构建的地电模型即可近似被认为是测区内的真实地电模型。根据测区内电性结构的分布特征,判断其相应的地质条件及圈定相应的地质异常体。

2 数值模拟

为了研究三维电阻率法对水电工程隧道地质异常体的响应程度,基于AGI软件平台,通过构建2个水电工程隧道地质条件的数值模型,模拟测区内有含水溶洞、断层破碎带及岩层含水时空间范围内电性分布特征。

模型如图1所示,在电阻率为1 000 Ω·m的均匀空间中有一个异常体,其长×宽×高为30 m×20 m×10 m,中心埋深为20 m,异常体分别设置为低阻(10 Ω·m)和高阻(10 000 Ω·m)。布设5条测线,电极间距为5 m,测线长145 m,测线之间的间距为20 m。数据采集方式为温纳三极装置。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   数值模型示意

Fig.1   Layout of model 1


本文研究的是水电工程隧道的地质条件,因为地层中含水量较大,因此模拟时围岩电阻率的取值较小,同时也有利于有效区分低阻异常体。

模拟结果如图2所示。图2a显示低阻异常体在深度25 m时表现最为明显,异常体的大小及分布位置与模型的位置对应准确,低阻异常体与周围均匀介质的分界线明显。因此,三维电阻率法对低阻异常体的探查具有很高的准确性,并且能够准确探查低阻异常体的位置。图2b中可以明显分辨出高阻异常体,在深度25 m时中心位置电阻率值由最高向四周扩散,最后趋于背景电阻率,且高阻异常体的空间分布与模型设置较为一致。高阻异常体与背景电阻率的响应特征都比较明显,模拟结果可以准确圈定高阻异常体的范围。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   模拟结果

Fig.2   Simulation results


通过数值模拟可以判断三维电阻率法能够准确探明水电工程隧道的地质条件,对隧道施工过程中遇到的溶洞问题能够作出明显的判别。因此,三维电阻率勘探技术对水电工程隧道地质条件的探查具有很好的适用性及可靠性。

3 现场实测

3.1 工区地质条件及测线布置

云南某水资源配置工程,需要利用隧道进行水资源的配置。该地区地貌以丘陵为主,且岩溶发育较多,隧道位于岩溶发育区域,场地经人工填平,地势较平坦。前期钻孔资料揭示工区主要分布第四系覆盖层,以粉质黏土为主,电阻率变化较大,一般为20~220 Ω·m;覆盖层下发育二叠系灰岩及泥盆系灰岩,电阻率一般为420~1 200 Ω·m。为探查该隧道的施工地质条件,在隧道前方布设7条测线(图3),线距均为20 m,每条测线均有90个电极,电极距1.5 m,测线方向与隧道开挖方向一致。采用温纳装置进行数据采集。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   工区测线布置

Fig.3   Layout of survey line


3.2 反演结果及综合解释

对7条测线的数据采集采用二维剖面采集模式,获得了质量较好的数据。合并数据后基于AGI软件平台进行三维反演。由于现场环境噪声较小,阻尼因子取值较小,初始阻尼因子为0.15,最小阻尼因子为0.01,结合隧道开挖的位置构建三维探测坐标。图4为反演得到的电阻率分布,可以看出剖面中电阻率具有明显的分层性:浅部低阻层约为40Ω·m,埋藏深度1~4 m,推测为工区的第四系覆盖层;其下方为明显的高阻层,电阻率为620~1 100 Ω·m,推测为二叠系灰岩。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   三维电阻率反演结果

Fig.4   Result of 3D resistivity inversion


同时,三维反演结果中存在2个明显的高阻异常体,电阻率值约为10 000 Ω·m。为了准确判断异常体的大小及空间分布,绘制了三维电阻率等值线图(图5)。图中左侧的高阻异常体较右侧的异常体大,结合现场的地质资料与数值模拟结果,可以判断这两个高阻异常体为溶洞,溶洞内含水极少,为溶洞空洞。左侧溶洞上部发育大,具有向左下延伸的趋势,埋藏深度较大,大小约为20 m×25 m×20 m;右侧溶洞较小,约为15 m×10 m×10 m,埋藏深度也较浅。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   高阻异常体的三维图示

Fig.5   3D diagram of high resistivity anomalous bodies


为了详细研究该隧道开挖过程中的隧道前方的地质条件,截取了三维反演结果中平行隧道测线的L4的电阻率剖面(图6)。可以看出:各地层岩性电性特征明显不同,地层之间具有较好的分层性,基岩面的起伏形态能够很好地分辨,溶洞的位置、大小及含水性均得到很好的显示。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   L4线的反演电阻率断面及解释成果

Fig.6   Inversion result and interpretation of survey line L4


3.3 钻孔验证

为验证三维电阻率法探测的准确性,在测线L4上方布设钻孔D1,图7给出了D1的钻孔柱状图及典型数字岩心图。钻孔揭示第四系覆盖层厚度为1.5 m,溶洞埋深6.4~10 m,溶洞内以空洞为主,未发现水及泥沙等破碎充填物。覆盖层厚度、基岩面埋深、溶洞的大小及充填物的性质与三维电阻率反演结果一致,取得了良好的探测效果。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   D1柱状图及典型数字岩心图

Fig.7   Bar chart and typical digital core chart of borehole D1


4 结论与建议

1)三维电阻率技术对隧道地质条件具有较好的响应程度,能够有效探测河床地质条件的三维特征,对断层、含水破碎带等异常体能够有效辨别,对判断溶洞的大小、空间分布以及充填物的性质具有良好的效果;三维电阻率技术对水电工程隧道的施工地质条件的空间探查具有很好的适用性。

2)利用三维电阻率技术对水电工程隧道进行现场探测,经钻孔验证,取得了良好的探测效果。三维电阻率技术能够准确判断其空间范围内的覆盖层厚度、基岩面分布特征。

3)三维电阻率技术能够有效探测隧道的空间范围内的地质条件,为隧道的安全开采提供保障,同时为类似工程的施工提供技术支持。

参考文献

张勇, 张子新, 华安增.

TSP超前地质预报在公路隧道中的应用

[J]. 西部探矿工程, 2001, 13(5):71-72.

[本文引用: 1]

Zhang Y, Zhang Z X, Hua A Z.

Application of TSP advanced geological prediction in highway tunnel

[J]. West-china Exploration Engineering, 2001, 13(5):71-72.

[本文引用: 1]

蔡盛.

张吉怀铁路隧道超前预报技术应用研究

[J]. 物探与化探, 2021, 45(5):1275-1280.

[本文引用: 1]

Cai S.

The research on the application of geological prediction technology to Zhangjihuai railway tunnel

[J]. Geophysical and Geochemical Exploration, 2021, 45(5):1275-1280.

[本文引用: 1]

钟世航, 王荣.

陆地声呐法探查中小溶洞、溶管效果好的原因

[C]// 中国地球物理学会第22届年会论文集. 北京: 中国地球物理学会, 2006:265.

[本文引用: 1]

Zhong S L, Wang R.

Reasons for good effect of small and medium-sized caves and dissolved tubes by terrestrial sonar

[C]// The Proceedings of 22th Annual Meeting of the Chinese Geophysical Society. Beijing: Chinese Geophysical Society, 2006:213.

[本文引用: 1]

薛国强, 李貅.

瞬变电磁隧道超前预报成像技术

[J]. 地球物理学报, 2008, 51(3):894-900.

[本文引用: 1]

Xue G Q, Li X.

The technology of TEM tunnel prediction imaging

[J]. Chinese Journal of Geophysics, 2008, 51(3):894-900.

[本文引用: 1]

邓国文, 王齐仁, 廖建平, 等.

隧道不良地质现象的探地雷达正演模拟与超前探测

[J]. 物探与化探, 2015, 39(3):651-656.

[本文引用: 1]

Deng G W, Wang Q R, Liao J P, et al.

Forward modeling and advanced detection of radar in adverse geological phenomena tunnel

[J]. Geophysical and Geochemical Exploration, 2015, 39(3):651-656.

[本文引用: 1]

李术才, 刘斌, 李树忱, 等.

基于激发极化法的隧道含水地质构造超前探测研究

[J]. 岩石力学与工程学报, 2011, 30(7):1297-1309.

[本文引用: 1]

Li S C, Liu B, Li S C, et al.

Study of advanced detection for tunnel water-bearing geological structures with induced polarization method

[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(7):1297-1309.

[本文引用: 1]

周文龙, 席超强, 胡泽安, 等.

利用并行电阻率法探测隧道施工地质条件

[J]. 工程地球物理学报, 2016, 13(2):166-169.

[本文引用: 1]

Zhou W L, Xi C Q, Hu Z A, et al.

The application of parallel resistivity method to geological conditions detection of tunnel construction

[J]. Chinese Journal of Engineering Geophysics, 2016, 13(2):166-169.

[本文引用: 1]

陆云祥, 陈华根, 陈建荣, 等.

高密度电阻率法在舟山市隧道勘察中的应用

[J]. 地球物理学进展, 2010, 25(6):2174-2183.

[本文引用: 1]

Lu Y X, Chen H G, Chen J R, et al.

Application of the high-density resistivity method to survey tunnels in Zhoushan city

[J]. Progress in Geophysics, 2010, 25(6):2174-2183.

[本文引用: 1]

蒋全科, 雷宛, 黄霄寒, 等.

应用综合电阻率法勘察隐伏断层

[J]. 成都理工大学学报:自然科学版, 2016, 43(3):378-384.

[本文引用: 1]

Jiang Q K, Lei W, Huang X H, et al.

An application of prospecting concealed fault with multi-resistivity method

[J]. Journal of Chengdu University of Technology:Science & Technology Edition, 2016, 43(3):378-384.

[本文引用: 1]

刘树才, 刘志新, 姜志海, 等.

矿井直流电法三维正演计算的若干问题

[J]. 物探与化探, 2004, 28(2):170-172,176.

[本文引用: 1]

Liu S C, Liu Z X, Jiang Z H, et al.

Some problems in 3d forward simulation of mine direct current method

[J]. Geophysical and Geochemical Exploration, 2004, 28(2):170-172,176.

[本文引用: 1]

Ruan B Y, Xiong B.

A finite element modeling of 3-D resistivity sounding with continuous conductivity

[J]. Chinese Journal of Geophysics, 2002, 45(1):124-131.

DOI:10.1002/cjg2.v45.1      URL     [本文引用: 1]

戴世坤, 凌嘉宣, 陈轻蕊, 等.

各向异性介质空间—波数混合域直流电阻率法三维数值模拟研究

[J]. 地球物理学报, 2022, 65(7):2729-2740.

[本文引用: 1]

Dai S K, Ling J X, Chen Q R, et al.

Three dimensional DC anisotropic resistivity modeling in a mixed space-wavenumber domain

[J]. Chinese Journal of Geophysics, 2022, 65(7):2729-2740.

[本文引用: 1]

罗登贵, 刘江平, 王京, 等.

活动断层高密度电法响应特征与应用研究

[J]. 地球物理学进展, 2014, 29(4):1920-1925.

[本文引用: 1]

Luo D G, Liu J P, Wang J, et al.

Study on the response characteristics and application of high density resistivity method of active fault

[J]. Progress in Geophysics, 2014, 29(4):1920-1925.

[本文引用: 1]

李忠平.

基于高密度电法温纳装置的三维电阻率反演应用

[J]. 地球物理学进展, 2020, 35(3):970-975.

[本文引用: 1]

Li Z P.

Application of 3-D resistivity inversion based on Winner device of high density electricity method

[J]. Progress in Geophysics, 2020, 35(3):970-975.

[本文引用: 1]

王志鹏, 刘江平, 易磊.2D、

3D高密度电法探测断层效果及其应用

[J]. 科学技术与工程, 2019, 19(25):75-82.

[本文引用: 1]

Wang Z P, Liu J P, Yi L.

Effect and application of 2D and 3D high density resistivity method for fault detection

[J]. Science Technology and Engineering, 2019, 19(25):75-82.

[本文引用: 1]

朱瑞, 李朝辉, 时向阳, 等.

三维高密度电法在隐伏断层探测中的应用

[J]. 人民黄河, 2019, 41(11):106-109,143.

[本文引用: 1]

Zhu R, Li Z H, Shi X Y, et al.

Application of 3D electrical resistivity tomography to buried fault detection

[J]. Yellow River, 2019, 41(11):106-109,143.

[本文引用: 1]

刘斌, 李术才, 聂利超, 等.

隧道含水构造直流电阻率法超前探测三维反演成像

[J]. 岩土工程学报, 2012, 34(10):1866-1876.

[本文引用: 1]

Liu B, Li S C, Nie L C, et al.

Advanced detection of water-bearing geological structures in tunnels using 3D DC resistivity inversion tomography method

[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(10):1866-1876.

[本文引用: 1]

吴荣新, 刘盛东, 张平松.

双巷并行三维电法探测煤层工作面底板富水区

[J]. 煤炭学报, 2010, 35(3):454-457.

[本文引用: 1]

Wu R X, Liu S D, Zhang P S.

The exploration of two-gateways parallel 3-D electrical technology for water-rich area within coal face floor

[J]. Journal of China Coal Society, 2010, 35(3):454-457.

[本文引用: 1]

谭磊, 张平松, 吴荣新.

隐伏洞体并行电法探查试验研究

[J]. 昆明理工大学学报:自然科学版, 2015, 40(2):38-43.

[本文引用: 1]

Tan L, Zhang P S, Wu R X.

Experimental study of parallel electrical imaging in concealed cave

[J]. Journal of Kunming University of Science and Technology:Natural Science Edition, 2015, 40(2):38-43.

[本文引用: 1]

张欣, 赵明阶, 汪魁, 等.

电法三维成像技术在隧道岩溶探测中的应用

[J]. 中国岩溶, 2016, 35(3):291-298.

[本文引用: 1]

Zhang X, Zhao M J, Wang K, et al.

Application of 3D electrical resistivity tomography to a tunnel in a Karst area

[J]. Carsologica Sinica, 2016, 35(3):291-298.

[本文引用: 1]

张平松, 刘盛东, 吴荣新, 等.

采煤面覆岩变形与破坏立体电法动态测试

[J]. 岩石力学与工程学报, 2009, 28(9):1870-1875.

[本文引用: 1]

Zhang P S, Liu S D, Wu R X, et al.

Dynamic detection of overburden deformation and failure in mining workface by 3d resistivity method

[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(9):1870-1875.

[本文引用: 1]

吴荣新, 刘盛东, 张平松, 等.

地面钻孔并行三维电法探测煤矿灰岩导水通道

[J]. 岩石力学与工程学报, 2010, 29(S2):3585-3589.

[本文引用: 1]

Wu R X, Liu S D, Zhang P S, et al.

Detection of limestone water-conducting channels in coal mine by parallel 3d electric method of surface boreholes

[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S2):3585-3589.

[本文引用: 1]

Loke M H, Dahlin T.

A comparison of the Gauss-Newton and quasi-Newton methods in resistivity imaging inversion

[J]. Journal of Applied Geophysics, 2002, 49(3):149-162.

DOI:10.1016/S0926-9851(01)00106-9      URL     [本文引用: 1]

Loke M H, Barker R D.

Practical techniques for 3D resistivity surveys and data inversion 1

[J]. Geophysical Prospecting, 1996, 44(3):499-523.

DOI:10.1111/gpr.1996.44.issue-3      URL     [本文引用: 1]

刘盛东, 吴荣新, 张平松, 等.

三维并行电法勘探技术与矿井水害探查

[J]. 煤炭学报, 2009, 34(7):927-932.

[本文引用: 1]

Liu S D, Wu R X, Zhang P S, et al.

Three-dimensional parallel electric surveying and its applications in water disaster exploration in coal mines

[J]. Journal of China Coal Society, 2009, 34(7):927-932.

[本文引用: 1]

任政勇, 汤井田.

基于局部加密非结构化网格的三维电阻率法有限元数值模拟

[J]. 地球物理学报, 2009, 52(10):2627-2634.

[本文引用: 1]

Ren Z Y, Tang J T.

Finite element modeling of 3-D DC resistivity using locally refined unstructured meshes

[J]. Chinese Journal of Geophysics, 2009, 52(10):2627-2634.

[本文引用: 1]

周文龙, 吴荣新, 肖玉林.

充水溶洞特征的高密度电阻率法反演分析研究

[J]. 中国岩溶, 2016, 35(6):699-705.

[本文引用: 1]

Zhou W L, Wu R X, Xiao Y L.

Back analysis of high density resistivity method in the water-bearing Karst cave

[J]. Carsologica Sinica, 2016, 35(6):699-705.

[本文引用: 1]

欧元超, 张平松, 李建宁, 等.

基于并行电法的重阳木根系空间分布探测试验研究

[J]. 科学技术与工程, 2017, 17(10):131-135.

[本文引用: 1]

Ou Y C, Zhang P S, Li J N, et al.

Bischofia polycarpa root spatial distribution detection test research based on parallel electrical method

[J]. Science Technology and Engineering, 2017, 17(10):131-135.

[本文引用: 1]

/

〈 〉

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址：北京市学院路29号航遥中心 邮编：100083
电话：010-62060192；62060193 E-mail：whtbjb@sina.com