综合物探方法在隐伏岩溶区的应用分析——以长宁页岩气区块宁2019H18钻井平台为例
蔡从德1, 谢林成1, 何智浩2, 朱霞飞3, 黄腾4
1.四川长宁天然气开发有限责任公司,四川 宜宾 644000
2.成都理工大学环境与土木工程学院, 四川 成都 610059
3.中冶成都勘察研究总院有限公司,四川 成都 610023
4.西南科技大学环境与资源学院,四川 绵阳 621010

第一作者简介: 蔡从德(1970—),男,高级工程师,主要从事地质勘察研究工作。Email: 315628878@qq.com

摘要

岩溶塌陷具有隐蔽性和突发性,是碳酸盐岩地区工程建设过程中常见的不良地质问题。长宁地区岩溶发育广泛,其页岩气开发利用受困于岩溶塌陷危害,综合运用微动、高精度重力、瞬态面波和等值反磁通多项物探方法,调查长宁页岩气区块宁2019H18钻井平台下伏地质体特征,查明土洞和溶洞的发育状况,结合钻孔数据进行验证,主要成果包括: ①微重力探测呈现重力值南东高、北西低,重力场响应特征表明该场地南东放向基岩密度较大且完整,而场地NE向基岩密度低; ②微动剖面响应结果与微重力成果对应较好,场地速度结构南东高、北西低,岩石的波速与密度呈正相关; ③进一步查明了场地下方软弱岩层的分布范围,场地SN向基岩完整性好,场地NW向存在着一条裂隙带,裂隙带内充填较为破碎松软的全风化砂岩,其承载力较低。通过物探三维可视化成果,可以为下一步场地建设指明方向,为类似工程案例提供科学依据。

关键词: 隐伏岩溶; 综合物探方法; 塌陷; 长宁页岩气区块
中图分类号:P681.7 文献标志码:A 文章编号:2095-8706(2023)04-0018-11
Application of comprehensive geophysical prospecting method in concealed karst area: A case study from Ning 2019H18 drilling platform in Changning shale gas area
CAI Congde1, XIE Lincheng1, HE Zhihao2, ZHU Xiafei3, HUANG Teng4
1. Sichuan Changning Natural Gas Development Co., Ltd, Sichuan Yibin 644000, China
2. College of Environment and Civil engineering Chengdu University of Technology, Sichuan Chengdu 610059, China
3. Chengdu Surveying Geotechnical Research Institute Co., Ltd of MCC, Sichuan Chengdu 610023, China
4. School of Environment and Resource, Southwest University of Science and Technology, Sichuan Mianyang, 621010, China
Abstract

Karst collapse is a common geological problem in the process of engineering construction in carbonate rock areas because of its concealment and sudden occurrence. The karst collapse is widely developed in Changning area, therefore, the exploitation and utilization of shale gas is hampered by extensive karst collapse. In this paper, the authors used micromotion, high-precision gravity, transient surface wave and equivalent anti-flux geophysical prospecting methods to investigate the subsurface geological condition of Ning 2019H18 drilling platform, and ascertained the development of earth cave and karst cave, which were tested by the drilling data. The main achievements are as follous. ①The microgravity is high in southeast and low in northwest, and the gravity response indicates that the southeast bedrock has large density and complete structure, and that the northwest bedrock has low density. ②The micromotion response profile is well corresponded with the microgravity, and the court page is high in southeast and low in northwest. And the wave velocity is positively related to its density. ③The distribution of underground soft strata was also verified, and the southeast bedrocks were completely preserved. A fractured zone exists in the northwest of the area, filled with cracked and soft completely weathered sandstane, which has low bearing capacity. The 3D visualized physical prospecting achievements could provide directions for the next establishment of the area, and some scientific references for similar engineering cases.

Keyword: concealed karst; comprehensive geophysical prospecting method; collapse; Changning shale gas area
0 引言

岩溶主要发育在碳酸盐岩地区, 属于化学溶蚀, 岩溶作用的条件可归纳为以下4点: 岩石可溶蚀性、岩石中裂隙发育、地下水的溶蚀力以及流动性, 地下水将上覆土体通过裂隙带走, 在土体掏空形成土洞[1, 2, 3]。岩溶塌陷具有隐蔽性和突发性的特点, 对上部工程的顺利建设和安全运营构成了潜在威胁, 是工程建设中亟待解决的工程地质问题[4, 5, 6, 7]。查明溶洞走向、大小、埋深, 揭示岩溶发育规律, 是治理岩溶危害的前提。由于传统的钻孔检测方法盲目性强、施工量大、易受人为因素影响, 利用地球物理勘探的方法可以避免对地基损害。同时, 开展多种物探方法可以针对同一工区地质条件进行综合解译, 所得成果相比单一物探方法会更为精确[8]。刘伟等[9]利用高密度点发和微动技术组合对广西来宾吉利塌陷进行调查, 圈出两个易塌陷区, 利用钻探验证了物探解译成果。McCormack等[10]在爱尔兰西部开展高密度电法和三维建模研究, 准确还原了下伏岩溶管道的展布形态。徐联泽等[8]运用地质雷达、瞬态面波和高密度电法等综合物探方法在场地岩溶调查中, 探明了下伏复杂岩溶的发育状况。

长宁公司页岩气区块隶属于四川省川东南宜宾— 自贡— 泸州一带属页岩气资源勘探优质区块, 地处四川盆地西南缘, 南临云贵高原, 地表水文地质、地形复杂, 碳酸盐岩广泛分布, 属于岩溶现象发育的区域[11, 12, 13]。页岩气钻井勘探过程中, 在前期1 000 m 以浅常受一些不良地质现象(岩溶、裂隙、滑坡、破碎堆积体等)影响, 造成卡钻、掉钻等孔内事故, 甚至引发页岩气钻井平台地面基础塌陷事故; 钻遇岩溶通道时, 会引起钻井液严重漏失。以上事故不仅造成了较大经济损失, 还严重危害区块的自然生态环境, 为长宁地区页岩气的勘探开发带来了困难[14, 15]

宁2019H18钻井平台位于宜宾市珙县曹营镇龙胜村4组, 距曹营镇约8.6 km。本次研究工作仅针对设在井场井架处附近以及平台硬化场地之上的两个钻井, 其中靠内测井架处于挖方区的钻井已完工, 而靠外侧井架处于填方区的钻井在钻遇过程中钻遇浅部砂泥岩堆积层, 造成事故发生。为了进一步查明钻井平台下的隐伏土洞和溶洞位置和洞隙分布情况, 本文利用地质雷达和瞬态面波组合方法探测隐伏土洞和高密度电法探测溶洞洞隙, 为该岩溶地区下一步的勘察和治理提供科学依据。

1 宁2019H18钻井平台地质条件
1.1 地层岩性

宁2019H18钻井平台地层主要有上三叠统— 下侏罗统香溪群(T3-J1xx)及中三叠统雷口坡组(T2l)(图1), 岩性主要为素填土(Q4ml)、第四系残坡积层(Q4(el+dl))(图2)。素填土(0.3~4.6 m)为黄灰— 黄褐色, 松散、稍湿的粉土夹碎石组成, 碎石含量15%~35%, 一般大小2~5 cm, 次棱角状, 成分主要为灰岩, 其余为黄灰— 黄褐色粉土(图2(a))。第四系残坡积层(Q4(el+dl))(4.6~21.8 m)为黄灰色、黄褐色, 可塑、稍湿的粉质黏土, 韧性中等, 干强度中等, 稍有光泽, 局部偶夹有碎石角砾(图2(b))。上三叠统— 下侏罗统香溪群(T3-J1xx)(20~25 m)为黄灰— 黄褐色的全风化砂岩组成, 结构构造特征已破坏, 岩石已风化呈砂土状, 手触即碎(图2(c))。中三叠统雷口坡组(T2l)(顶部界面 25 m)为灰色— 深灰色的强风化灰岩组成, 微晶结构, 中层状构造, 节理较为发育, 局部表面可见有溶蚀孔洞和裂隙。岩心破碎, 多呈碎块状。其中, 中风化灰岩为灰色— 深灰色, 微晶结构, 中层状构造, 节理较为发育。局部表面可见有溶蚀孔洞和裂隙。岩心较为破碎, 多呈碎块状、短柱状, 节长为10~25 cm。勘探深度范围内未揭穿。溶洞为石灰岩地层中溶蚀形成的溶洞或裂隙, 发育及分布无规律, 纵向高度0.3~5.8 m。

图1 研究区区域地质简图
1.第四系; 2.白垩系夹关组; 3.白垩系天马山组; 4.侏罗系沙溪庙组; 5.侏罗系自流井组; 6.三叠系须家河组; 7.三叠系雷口坡组; 8.上二叠统; 9.下二叠统; 10.志留系; 11.奥陶系桐梓组; 12.寒武系; 13.钻井平台位置; 14.地名
Fig.1 Regional geological sketch of the study area

图2 2019H18钻井平台松散沉积物和灰岩岩心Fig.2 Loose sediment and limestone core of 2019H18 drilling platform

1.2 水文地质

根据地下水赋存条件、含水介质的岩石特征等划分, 本区地下水主要有松散堆积物孔隙水和岩溶水, 主要特征如下。

(1)松散堆积物孔隙水赋存于第四系残坡积层中, 以潜水为主, 水位埋深变化大, 含水层较薄, 分布面积小; 受季节性影响明显, 含水条件和水力联系均较好, 渗透性较强, 与地表水有密切的水力联系, 在不同的地段表现为互补关系。这类地下水主要接受地表水和大气降水补给, 具有地下径流较短、甚至就近补给排泄、水位受季节(丰水期、枯水期)及地表水体影响变化幅度较大等动态特征。

(2)碳酸盐岩岩溶水主要赋存于三叠系中统雷口坡组(T2l)的灰岩中, 由于地下水的赋存条件和裂隙岩溶发育程度的差异, 各地层中的富水性很不均一。以大气降水补给为主, 由于基岩裂隙分布不均匀, 大多不具有统一水力联系, 在接受大气降水和地表水补给后, 表现出渗流明显受裂隙岩溶控制的特点, 运移带有局限性, 多以下降泉形式排向河谷、沟谷。此类型水的存在可降低岩体的抗剪强度, 加快岩石风化, 降低岩体的稳定性。勘察时未发现碳酸盐岩岩溶水。

(3)岩石的电性特征。根据岩土体的物性经验资料, 不良地质体与围岩的物性分析如表1, 不良地质体与围岩物性差异明显, 具备地球物理工作的前提条件。

表1 岩土介质的电性特征 Tab.1 Electrical characteristics of geotechnical media
2 综合物探工作方法及测线布置
2.1 综合物探工作方法

本次研究综合采用微动、高精度重力、瞬态面波和等值反磁通多项物探方法, 对场地下伏地质体情况进行了调查(表2)。

表2 宁2019H18钻井平台工作量统计 Tab.2 Workload statistics of Ning 2019H18 drilling platform

2.1.1 微动

微动法勘探深度变化幅度大, 浅层勘探可以在数米之内, 深部勘探现今已有在页岩气勘探中开展微动法并取得一定效果的先例[21]。此外微动台阵的不同勘探深度也有所不同, 环形台阵勘探深度大约为最大半径的3~5倍, 线形台阵的勘探深度大约为排列长度的0.7~1.0倍。试验中发现在投入相同数量检波器的情况下, 线形台阵反演成果不同深度的异常分布较均匀, 环形台阵深部异常的细节则不够突出, 推测可能与线形台阵检波器分布均匀有关。勘探深度大时, 需要增加采集时间, 还需要加大排列长度或台阵半径, 检波器数量不足会带来细节上的损失, 增加检波器数量又会带来成本上的增加。微动法主要依据是分析计算列阵中心地下的横波速度随深度的变化情况推断异常值。当工作区为单一土层时, 三分量记录曲线的周期彼此相同或接近, 且卓越周期只有一个; 当工作区为多层土且厚度较大时会出现多个峰值, 每个峰值代表该层与上覆层具有某种联系[16]。本次研究采用的微动点距5 m, 直线形台阵布设。

2.1.2 高精度重力

地面微重力方法不受电磁场等人为干扰和接地条件的影响及工作场地大小等因素的限制, 对于埋藏浅的微小探测对象及其低缓微弱异常具有较高的分辨能力, 并且野外工作方法简单, 成本低, 效率高, 从而弥补了某些物探方法的不足。高精度重力法依据不同岩性、不同年代地层和局部地质体的密度值具有明显的差异来推断异常。但岩溶溶洞引起的重力异常幅值通常在几十微伽级, 为区分局部异常和背景异常, 需参照钻孔剖面正演得到的模型[17]。高精度重力采用重力点距2.5 m, 读数时间为45 s。

2.1.3 瞬态面波

混合源面波中的被动源勘探深度与微动类似, 主动源面波的勘探深度受震源能量、排列长度等限制, 一般约为20~30 m。被动源面波勘探深度较大, 视投入检波器数量以及所选台阵而不同。该方法能够较为准确地划分第四系, 通过与钻孔分层对比可知该方法划分的第四系黏土层和砂卵砾石层与钻孔的深度误差在5%以内, 该方法划分的基覆界面与钻孔的深度误差在5%左右, 反映该方法探测精度较高。瞬态面波可以通过锤击落重等震源, 产生一定频率范围的瑞利波, 并通过振幅谱和相位谱分析把不同频率的瑞利波分离出来, 得出频谱曲线[18]。依据不同岩土面波波速绘制频谱曲线, 并结合不同深度地层的瑞利面波对地下介质结构差异进行分析, 从而得出地层结构和特点的分析预测[19]。本次采用的是多道瞬态地震采集系统sinoX, 运用落锤震源, 偏移距6 m, 道间距2.5 m, 点距5 m, 24道检波器, 时间窗口为300 ms。在被动源面波采集数据过程, 我们改进的增加人工震源这一措施, 只适用于浅层勘探, 此时并非为了增加勘探深度, 而是增加高频信号, 提高浅表的分辨率, 但仍然会存在0.5~1.0 m左右的盲区。

2.1.4 等值反磁通

勘探深度受瞬变电磁法的供电电流、磁矩大小影响, 一般有效的深度在200 m以浅, 采用小发射线框的瞬变仪勘探深度大致在50 m以浅。该方法对于浅部的黏土层识别效果较好, 与钻孔的深度误差较小, 但其识别浅部砂卵砾石层的效果一般; 在其探测深度范围内能够识别出基覆界面, 分层误差在10%以内。等值反磁通瞬变电磁法通过收集地下到地表的磁场变化推断异常值。大地的电导率和局部良导体的埋深影响涡流扩散速度和衰减幅度, 均匀大地介质的电导率一般为常数, 局部地质体出现的电导率、埋深、形态等参数变化会影响到磁场的极大值的扩散速度和衰减幅度改变, 从而识别出异常值[20]。等值反磁通瞬变电磁点距5 m, 采集频率6.25 Hz。

2.2 测线布置

宁2019H18钻井平台工区共设计微动测线4条, 重力测线4条, 等值反磁通测线4条, 瞬态面波测线1条, 4条测线平行布置, 测线方位角为137° 。测线沿钻井平台基础布设, 宁2019H18钻井平台研究区测线布置工作如图3表2所示。

图3 宁2019H18钻井平台测线布置Fig.3 Layout of Ning 2019H18 drilling platform survey line

(1)微动、重力和等值反磁通方法布设4条测线: 10线、11线、12线和13线, 4条测线均为NE— SW向, 作为综合物探剖面, 研究不同方法在岩溶勘察中的作用。

(2)瞬态面波方法布设1条测线, 即12测线, 测线方向为NE— SW向。

3 综合物探方法解译结果
3.1 L10测线物探推断成果

测线L10微重、微动和反磁通量的反演结果如图4所示。在重力异常曲线图中(图4(a)), 26号重力点位的局部重力异常中呈现一处较为明显的局部重力低异常, 在水平位置上与微动反演剖面比较接近, 疑似溶洞或者类似松散体的异常特征。基覆界面介于17~20 m, 从小号点到大号点逐渐增加, 与重力背景场分布比较接近(图4(b))。图4(b)中红色虚线圈所示位置为一低速异常区域, 中心值为500~600 m/s之间, 与该区域溶洞波速接近, 疑似溶洞充填泥沙。等值反磁通结果如图4(c), 从测线75 m 处开始出现两处电阻局部低异常(红色圈闭处), 水平位置和垂直位置均与微动反演结果相似。但并没有在测线40 m附近勘测到低电阻异常区。

ZK1钻孔资料显示, 深度为20.4 m和40.7 m处为2个明显的分界点。0~20.4 m为黄灰色、黄褐色的第四系粉质黏土、粉质砂土和碎石土; 20.4~40.7 m段为全风化砂岩, 且岩石结构特征已被破环, 呈现风化砂土状, 手摸即碎的特征; 40.7 m 以下为中风化灰岩层为主, 勘探深度范围内未揭穿。在深度59.8 m处出现全风化砂岩异常(图2), 推测与测线位置65~70 m、深度50~60 m处的溶洞松散体有关(图4(b))。

图4 L10物探结果解译图Fig.4 Interpretation diagram of L10 geophysical exploration results

可见, 反磁通的探测结果在水平位置上稍有误差, 且未探出测线40 m深度20~40 m的全风化砂岩异常。微动和微重的探测结果与ZK1工程钻孔分层结果基本一致, 探测结果较为可靠。

3.2 L11测线物探推断成果

L11测线微动、微重和反磁通量反演剖面如图5所示: 在重力点位23附近的重力局部场有一重力局部低值, 由于该溶洞规模较小, 所以在重力异常上体现不明显, 但仍有响应(图5(a))。基岩面在15~22 m之间, 其趋势和重力区域场相近(图5(b))。在图5(b)中红色虚线圈所示的异常区域, 其波速中心值为250~300 m/s, 低于该区域泥沙充填溶洞的波速值为疑似溶洞, 充填物未知。等值反磁通结果如图5(c), 从测线60 m开始出现电阻局部低异常(红色圈闭处), 较浅的两个电阻异常位置与微动反演一致。但深度在60 m处的电阻异常与微动反演不符合, 这可能是由于等值反磁通法对于深部的探测效果较差而产生的误差。

图5 L11物探结果解译Fig.5 Interpretation diagram of L11 geophysical exploration results

ZK3钻孔资料显示, 钻孔总深为60.2 m, 0~21.8 m为黄灰色、黄褐色的第四系松散粉质黏土, 局部偶夹有碎石角砾; 21.8~25.1 m段为全风化砂岩层; 25.1~31.4 m为强风化石灰岩, 岩心破碎, 多层碎块状; 31.4~60.2 m为中风化灰岩层, 节理较为发育, 岩心较为破碎, 在43~60.2 m段, 局部表面可见有溶蚀孔洞和裂隙。

可见, 除了反磁通在测线40~50 m、深度40~ 60 m处的探测结果与ZK03的分层结果不符之外, 微动、微重与反磁通整体的探测结果与ZK3工程钻孔分层结果基本一致, 探测结果较为可靠。

3.3 L12测线物探推断成果

测线L12微动、微重、面波和反磁通量的反演结果如图6所示。ZK5号钻孔位于测线L12上65 m 左侧, 反演结果显示表层覆盖层约20~25 m, 其趋势重力区域场一致(图6(a))。等值反磁通反演结果如图6(d)所示, 该钻孔进钻30 m见溶洞, 溶洞顶板约30 m, 底板约47 m, 在图6(b) 面波反演结果和图6(c)微动反演剖面中表现为相对低速异常区域, 溶洞充填物为泥沙, 其中心值介于400~500 m/s 之间。从测线65 m处开始, 有两处电阻率局部低异常分别位于图中圈闭位置, 在水平位置与垂向位置上与面波和微动反演结果比较接近。

图6 L12物探结果解译Fig.6 Interpretation diagram of L12 geophysical exploration results

ZK5钻孔揭露的地层与ZK01相似, 0~24.1 m 为黄灰色、黄褐色的第四系松散粉质黏土, 可塑、稍湿, 局部偶夹有碎石角砾; 24.1~30 m 段为灰色— 深灰色的强风化石灰, 岩心破碎, 多层碎块状; 31.4~60.2 m为灰色— 深灰色的中风化灰岩层, 中层状构造, 节理较为发育, 但在31.4~40.7 m段, 岩心更破碎, 局部表面可见有溶蚀孔洞和裂隙。60.2 m处出现了较为紧密的粉质黏土层(图2), 这可能是溶洞塌陷堆积的松散体形成。

可见, 微动、微重和瞬态面波法的探测结果与ZK5工程钻孔分层结果基本一致, 探测结果较为可靠。

3.4 L13测线物探推断成果

测线L13微动、微重和反磁通量的反演结果如图7所示。基岩面在20~25 m之间, 表层疑似有风化破碎层, 与重力背景场分布接近(图7(a)), 比较明显的速度分界面如红线所示(图7(b)), 图中两处圈闭位置所示的局部波速低异常, 其中心值在400~500 m/s之间, 疑似岩溶。在该测线的等值反磁通反演剖面上共有4处电阻率局部低异常区域, 其中浅层40 m深两个低阻异常在空间位置上与微动剖面低速异常位置接近, 为疑似溶洞区域(图7(c))。

图7 L13物探结果解译Fig.7 Interpretation diagram of L13 geophysical exploration results

ZK7钻孔资料显示, 深度17.2 m和23.2 m为两个明显的分层点。0~17.2 m为黄灰色、黄褐色的第四系松散粉质黏土、粉质砂土以及一些碎石角砾; 17.2~23.2 m段为黄灰色、黄褐色的全风化砂岩层, 结构特征已被破坏; 23.2~27.8 m 段为灰色— 深灰色的强风化石灰, 岩心破碎, 多层碎块状; 27.8~68.5 m为灰色— 深灰色的中风化灰岩层, 中层状构造, 节理较为发育, 但在27.8~52.5 m段, 岩心更破碎, 局部表面可见有溶蚀孔洞和裂隙。

可见, 反磁通的探测结果在0~30 m段探测效果不如微动, 未能探出较为浅层的局部异常, 而微动和微重的探测结果与ZK7工程钻孔分层结果基本相符, 探测结果较为可靠。

4 宁2019H18钻井平台岩溶勘察三维可视化分析

将宁2019H18钻井平台微动视横波速度反演数据进行三维网格化成图, 可以清楚直观地发现在基覆界面的位置及埋深, 强风化基岩层的厚度以及隐伏岩溶在三维空间内的展布、形态及其规模: 基覆界面埋深15~25 m, 基岩地层内存在多个大小不等的隐伏岩溶, 最大的岩溶位于平台中部L12测线70 m处, 岩溶直径可达10 m以上, 其余岩溶直径0.5~5 m 不等, 为进一步摸清钻井平台场地下部地质空间结构提供了有效的指导依据。

(1)在宁2019H18钻井平台上进行了微动、瞬态面波、等值反磁通和高精度重力四种物探方法的场地试验, 并着重的对已进钻打孔的L12测线进行了探测, 得到了较好的试验结果。

(2)通过在该平台上开展的微重力测量结果表明: 微重力法能够有效地应用于隐伏岩溶勘察、基覆界面形态刻画等方面, 其横向分辨率高, 当重力测点下方存在软弱土体或全风化基岩时, 重力值呈明显的下降趋势, 当测点下方基岩较为完整且覆盖层较薄时, 重力值会呈现上升趋势, 利用该重力特征可以十分有效地在平面上圈定出软弱地层、含隐伏岩溶地层的平面位置, 经钻孔约束后还能反演出隐伏岩溶的规模大小及其埋深情况。

(3)利用在该平台上开展的微动和混合源瞬态面波工作成果, 可以有效地在垂向上识别出具有传播速度差异的岩层分界线, 通常基岩与第四系岩性差异大, 速度变化大, 可以利用该特点识别基覆界面埋深及其起伏形态; 而隐伏岩溶更是与围岩形成了明显的传播速度差异, 因此利用该方法识别隐伏岩溶效果显著。

(4)等值反磁通瞬变电磁法对场地条件要求较高, 在施工过程中场地大量积水, 该方法探测效果不明显。

(5)综合利用微重力和地震弹性波法(微动、面波法)能有效地识别隐伏岩溶的空间位置, 共圈定了隐伏岩溶27个(图8), 经钻孔验证探测结果可靠。该套方法组合可在长宁区块页岩气钻井平台钻前工程建设中推广应用。

图8 宁2019H18平台异常平面示意图Fig.8 Plan diagram of Ning 2019H18 platform anomalies

5 结论

基于已取得良好应用效果的“ 微重力+微动” 探测技术组合模式, 将其应用于宁2019H18钻井平台钻前工程场地地下地质探测, 取得4点研究成果。

(1)宁2019H18钻井平台第四系覆盖层较薄, 厚度约10 m, 基覆界面埋深10 m以浅。

(2)微重力探测结果呈现出重力值南东高、北西低, 重力场响应特征表明了该场地SE方向基岩密度较大, 基岩较为完整, 而场地NE向基岩密度低, 推测由松散软弱岩层和裂隙带引起。

(3)微动剖面响应结果与微重力成果对应较好, 场地速度结构南东高、北西低, 岩石的波速与密度呈正相关, 结果表明两种方法探测结果一致, 经过剖面工程钻孔验证, 精度较高。

(4)根据收集的该场地地勘成果报告, 场地下伏基岩为三叠系须家河组砂岩地层为主, 局部夹杂泥土层, 岩性风化程度较高, 承载力较低。物探探测成果也表明了该特征, 又进一步查明了场地下方软弱岩层的分布范围, 场地SE向基岩较为隆起, 岩石密度大, 速度高, 完整性好, 场地NW向存在着一条裂隙带, 裂隙带内充填较为破碎松软的全风化砂岩, 其承载力较低, 通过物探三维可视化成果, 可以为下一步场地建设指明方向, 规避井架沉降等风险隐患。

(责任编辑: 常艳)

参考文献
[1] 罗小杰, 罗程. 岩溶地面塌陷三机理理论及其应用[J]. 中国岩溶, 2021, 40(2): 171-188.
Luo X J, Luo C. Three-Mechanism Theory(TMT) of karst ground collapse and its application[J]. Carsologica Sinica, 2021, 40(2): 171-188. [本文引用:1]
[2] 刘洋, 郑俊杰, 曾彦. 基于极限分析下限法的溶洞顶板承载力研究[J]. 岩土工程学报, 2019, 41(S2): 181-184.
Liu Y, Zheng J J, Zeng Y. Bearing capacity of karst roof based on lower bound method[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(S2): 181-184. [本文引用:1]
[3] 钟宇, 张明堃, 潘玲, . 基于GIS的武汉市武昌区岩溶塌陷危险性评价[J]. 天津师范大学学报: 自然科学版, 2015, 35(1): 48-53.
Zhong Y, Zhang M K, Pan L, et al. Risk assessment for urban karst collapse in Wuchang District of Wuhan based on GIS[J]. Journal of Tianjin Normal University: Natural Science Edition, 2015, 35(1): 48-53. [本文引用:1]
[4] 潘健, 周森, 林培源, . 广州市白云区岩溶塌陷风险初探[J]. 岩土力学, 2013, 34(9): 2589-2600.
Pan J, Zhou S, Lin P Y, et al. Preliminary study of risks of karst collapse in Guangzhou Baiyun district[J]. Rock and Soil Mechanics, 2013, 34(9): 2589-2600. [本文引用:1]
[5] 肖先煊. 覆盖型岩溶区水气相互驱动盖层变形演化及塌陷机理研究[D]. 成都: 成都理工大学, 2018.
Xiao X X. Deformation Behavior Evolution and Collapse Mechanism of Karst Covers Under Water-Air Interaction in Karst Area[D]. Chengdu: Chengdu University of Technology, 2018. [本文引用:1]
[6] 张丽芬, 曾夏生, 姚运生, . 我国岩溶塌陷研究综述[J]. 中国地质灾害与防治学报, 2007, 18(3): 126-130.
Zhang L F, Zeng X S, Yao Y S, et al. Review on karst collapse in China[J]. The Chinese Journal of Geological Hazard and Control, 2007, 18(3): 126-130. [本文引用:1]
[7] 徐联泽. 综合物探方法在复杂岩溶勘察中的应用[J]. 科技经济导刊, 2020, 28(33): 51-52.
Xu L Z. The application of comprehensive geophysical methods in complex karst exploration[J]. Technology and Economic Guide, 2020, 28(33): 51-52. [本文引用:1]
[8] 刘伟, 甘伏平, 赵伟, . 高密度电法与微动技术组合在岩溶塌陷分区中的应用分析——以广西来宾吉利塌陷为例[J]. 中国岩溶, 2014, 33(1): 118-122.
Liu W, Gan F P, Zhao W, et al. Application analysis of combining high density resistivity and microtremor survey methods in areas of karst collapse: A case study of the collapse in Jili village, Laibin, Guangxi[J]. Carsologica Sinica, 2014, 33(1): 118-122. [本文引用:2]
[9] 王鹏, 谢丹, 罗乃菲. 长宁地区页岩气地质特征与勘探开发前景[J]. 宜宾学院学报, 2020, 20(12): 10-15, 99.
Wang P, Xie D, Luo N F. Geological characteristics of shale gas and prospects for exploration and development in Changning area[J]. Journal of Yibin University, 2020, 20(12): 10-15, 99. [本文引用:1]
[10] 张坦. 川南地区断裂特征与构造样式分析[D]. 北京: 中国地质大学, 2020.
Zhang T. Fracture Characteristics and Construal Styles in Southern Sichuan Area[D]. Beijing: China University of Geosciences, 2020. [本文引用:1]
[11] 赵瑞. 四川盆地南缘地形梯度带区域岩溶水系统研究[D]. 成都: 成都理工大学, 2016.
Zhao R. Research on the System of Regional Karst Water in the Topographic Gradient Zone of the South Edge of Sichuan Basin[D]. Chengdu: Chengdu University of Technology, 2016. [本文引用:1]
[12] 余长恒, 周昊, 邹忠平, . 长宁地区页岩气钻井平台不同开孔层位不良地质体勘查[J]. 工程地球物理学报, 2019, 16(1): 61-69.
Yu C H, Zhou H, Zou Z P, et al. The exploration of unfavorable geological body of shale gas drilling platform at different openings in Changning Area[J]. Chinese Journal of Engineering Geophysics, 2019, 16(1): 61-69. [本文引用:1]
[13] 余长恒, 张旭林, 王强, . 高密度电法在岩溶勘查中的参数试验——以四川南部宜宾市长宁页岩气开发区为例[J]. 矿产勘查, 2020, 11(9): 1986-1992.
Yu C H, Zhang X L, Wang Q, et al. Parameter test of high density electrical method in Karst exploration: Take the Changning shale gas development area in Yibin, southern Sichuan Province as an example[J]. Mineral Exploration, 2020, 11(9): 1986-1992. [本文引用:1]
[14] Foudili D, Bouzid A, Berguig M C, et al. Investigating karst collapse geohazards using magnetotellurics: A case study of M’rara basin, Algerian Sahara[J]. Journal of Applied Geophysics, 2019, 160: 144-156. [本文引用:1]
[15] McCormack T, O’Connell Y, Daly E, et al. Characterisation of karst hydrogeology in Western Ireland using geophysical and hydraulic modelling techniques[J]. Journal of Hydrology: Regional Studies, 2017, 10: 1-17. [本文引用:1]
[16] 张伟, 甘伏平, 梁东辉, . 多道瞬态面波法与微动法联合推断岩溶塌陷区覆盖层的结构和厚度分布[J]. 中国岩溶, 2018, 37(2): 272-279.
Zhang W, Gan F P, Liang D H, et al. Joint inference of multi-channel transient surface wave method and microtremor method to investigate the structure and thickness of the overburden layer in the typical area of karst collapse[J]. Carsologica Sinica, 2018, 37(2): 272-279. [本文引用:1]
[17] 徐燕君, 陈明, 伍卓鹤, . 高精度重力在城市地下空间探测研究中的应用[J]. 地质学刊, 2021, 45(4): 425-432.
Xu Y J, Chen M, Wu Z H, et al. Application of high precision gravity in the research of urban underground space exploration[J]. Journal of Geology, 2021, 45(4): 425-432. [本文引用:1]
[18] 李巧灵, 张辉, 雷晓东, . 综合利用多道瞬态面波和微动探测分析斜坡内部结构[J]. 物探与化探, 2022, 46(1): 258-267.
Li Q L, Zhang H, Lei X D, et al. Analyses of internal structure of slopes using multi-channel transient surface wave exploration and microtremor survey[J]. Geophysical and Geochemical Exploration, 2022, 46(1): 258-267. [本文引用:1]
[19] 王鑫文. 瞬态面波勘探技术在岩土工程勘察中的应用[J]. 中国水运(下半月), 2017, 17(11): 233-234.
Wang X W. Application of transient surface wave exploration technology in geotechnical investigation[J]. China Water Transport, 2017, 17(11): 233-234. [本文引用:1]
[20] 席振铢, 龙霞, 周胜, . 基于等值反磁通原理的浅层瞬变电磁法[J]. 地球物理学报, 2016, 59(9): 3428-3435.
Xi Z Z, Long X, Zhou S, et al. Opposing coils transient electromagnetic method for shallow subsurface detection[J]. Chinese Journal of Geophysics, 2016, 59(9): 3428-3435. [本文引用:1]
[21] 郑福龙, 余舟, 蔡从德, . 微重力和微动法相结合的隐伏岩溶探测效果分析——以长宁某页岩气钻井平台为例[J]. 地球物理学进展, 2023, 38(2): 967-976.
Zheng F L, Yu Z, Cai C D, et al. Analysis of the effect of combining microgravity and microdynamic methods for the detection of hidden karst: Taking a shale gas drilling platform in Changning as an example[J]. Progress in Geophysics, 2023, 38(2): 967-976. [本文引用:1]