各种地球物理勘探方法在查找地下洞穴中的应用研究.doc

1各种地球物理勘探方法在查找地下洞穴中的应用研究孙跃军 （广东省地质物探工程勘察院 广州 510800）摘要：本文通过对各类地下洞穴探测的总结分析和研究，根据物性差异、勘查目的和要求及场地条件等方面的不同，设计出相应的地球物理勘探方法和手段，并通过工程实例进行了论证说明，进而得出地下洞穴探测的最佳方法和手段关键词 地球物理勘探 采空区 岩溶 路面脱空 防空洞 方法选择1. 引言在一些工程建设、矿产资源开采和环境评估中，常会遇到各种不同类型的地下洞穴，使得项目的安全和质量受到极大的影响因此，项目的设计施工前都必须查明场地内地下洞穴的规模及分布情况地下洞穴有天然生成的，比如溶洞、土洞等；也有人为因素而成的，比如采空区、防空洞等地下洞穴有的存在充填物（水、流塑状物质等） ，有的为空洞因此地下洞穴与周围岩土层存在明显的物性差异，可采用地球物理勘探方法进行探测但由于场地条件、勘察目的及方法本身的局限性等方面的原因，对不同类型的地下洞穴采用的物探方法也存在差异2. 采空区的探测2.1 采空区的地球物理特性存在矿产资源的地层由于成因不同而表现出不同的分布状态，多呈现为层状沉积，表现为层状展布。

当矿脉被开采后形成采空区，破坏了原有的应力状态，使原有地层错动、产生裂缝、发生塌陷等一般情况下，采空区表现为两种形式：当开采时间较短或支护措施较充分时，其矿脉层顶板为塑型岩石并保存完整，部分压力转移到地层上，未引起地层塌落、变形，采空区以不充水或充水的空洞形式保存下来；当采面较大或开采时间较长时，采空区在重力和地层应力作用下，顶板塌落， 形成冒落带、裂隙带和弯曲带等， 使得采空区被上层塌落的松散物充填而地球物理特征发生变化，导致塌落带上部地层的特性也发生了变化 2.2 探测方法选择对于采空区的探测方法有浅层地震法、电法、瞬变电磁法等一般来说，存在矿脉、煤田等矿产资源的地层地质构造和岩性都很复杂，开采前高低阻体就较多，采用2依靠电性差异进行勘探的电磁类物探方法对采空区的异常反映较为复杂，较难分辨真假异常，因此浅层地震反射波法才是最好的复杂地层采空区的勘探方法地震波中含有大量的地质信息，地层构造、岩性变化、矿脉厚度变化、分岔合并等均会引起地震波的变化，其变化主要包括密度、速度及其它弹性参量的差异，这些差异导致了地震波在传播时间、振幅、相位、频率等方面的变化当采空区未被充填时，矿脉反射波消失；当上覆地层塌落、回填，则采空区的反射波凌乱、一定范围内的反射波同相轴向下错动，呈现出松散、塌落的形态。

在解释采空区时，以波形变面积显示的时间剖面为主，结合彩色剖面显示识别采空区重点放在标准反射波同相轴的消失、变弱、 强相位转换及在采空区边界处绕射波振幅的强弱变化等多种因素，来确认采空区的存在及其范围2.3 实例分析2.3.1 工区概况广东省高要市长坑金矿在二十世纪九十年代进行过开采，为了评估矿产剩余储量，进行了本次物探因此本次物探主要目的和任务是查明采空区走向及其埋深，为评估矿产剩余储量提供有力的依据据地质勘测资料，场地岩土层主要包括：a. 第四系（Q）土层：主要为洪冲积、残坡积层和填土等组成，厚度 0～10m，平均约 3m；b. 石炭系灰岩（C 1Z）：分布于全区，工区西北角出露，其它地方被覆盖，该层岩溶发育，厚度 10～100m；c. 三迭系泥岩、砂岩、页岩（T 3xly 和 T3xf）：裂隙发育，且多为方解石脉充填；d. 断裂破碎蚀变带：主要由硅化岩、灰岩组成，裂隙发育，且充填，导水性差；e. 金矿体2.3.2 地球物理特征区内各主要矿脉埋藏较浅，厚度在 0~100m 之间，与围岩能形成较好的波阻抗界面， 形成较好的反射波由于矿脉采空区被空气或水填充，而地震波在空气或水中传播的速度比在矿脉中传播的速度要低得多( 地震波在空气中传播的速度约为340m／s 、在水中传播的速度约为 1400m／s、在矿脉中传播的速度在 3000m／s 以上)，故产生的反射波相位会有所延迟，所以把在连续的矿脉反射波发生显畸变的地方解释为采空区。

2.3.3 资料处理与成果推断资料的处理流程如下：解编→预处理→静校正→叠前滤波→干扰切除→道间平衡→抽 CDP 道集→速度分析→动校叠加→叠后滤波→修饰处理→剖面打印→初步解释→成果图输出3推断采空区的推断主要依据：①、存在弧形地震波组；②、地震波组突然出现错断或杂乱；③、波组振幅在一定走时以下突然减弱甚至歼灭如图 1 所示的 A 线 259-396 段时间剖面，T1 波组时间在 10ms~40ms 之间变化，反映了覆盖层的深浅变化，中间部分不连续表示地层的起伏过大、错断缺失或存在局部异常体等，并且受到采空区影响在 40ms 以后的地层信息反映较差在点号320～349、350～370、371～396 段分别出现弧形地震波组，尤其是在320～349、371～396 段弧形异常明显，原地层波组几乎缺失，289～396 段 120ms 以下的地震波振幅明显减弱甚至歼灭，因此结合地质资料可以推断 280～396 段存在一个大的采空区，深度范围在 20～80m 之间图 1 A 线纵波反射时间剖面图如图 2 所示的 L08-1 线时间剖面， T1 波组时间在 20ms~40ms 之间变化，反映了覆盖层的深浅变化，地层基本连续。

在点号 182～277、227~272 段出现明显的弧形波组，结合地质资料可以推断 182～272 段存在采空区，宽度约 90m，厚度在 30～70m之间4图 2 L08-1 线纵波反射时间剖面图2.3.4 结论综上所述并通过钻孔验证，说明浅层地震反射波法对采空区的探测是非常有效的一种物探方法根据不同情况可以采用综合物探的方法会更精确地查找出目的体3. 岩溶的探测3.1 岩溶的地球物理特性岩溶一般存在灰岩地区当灰岩遭受强烈溶蚀作用将发育溶洞，在地下水作用下，上覆土层中的砂土沿溶洞通道被带走，即在土层中发育土洞土洞和溶洞内一般存在充填物溶洞内充填物与周围岩体之间存在明显的物性差异，土洞内充填物与覆盖层之间也存在一定的物性差异3.2 探测方法选择对于岩溶探测的主要方法有电法、弹性波 CT 或电磁波 CT、管波探测和探地雷达等对各种探测方法的选择因勘察阶段不同而不同，一般来说在场地普查阶段，常采用电法中的复合联合剖面法，在场地初勘阶段常采用高密度电法或高密度单点测深法，在在场地详细勘察阶段常采用弹性波 CT 或电磁波 CT，在桩位超前钻或施工阶段常采用管波探测如果覆盖层较浅（5m 以内） ，还可以采用探地雷达直接对场地进行扫描。

其中管波探测是一种最新的针对桩位岩溶探测而发明的方法，可探测出孔旁的岩溶发育情况，精度相当高，可为桩基的设计提供最可靠的地质依据3.3 弹性波 CT 实例分析3.3.1 工区概况广州白云国际会议中心场地（原东方乐园）根据区域资料，探测地段位于广花盆地内，地层走向为北东，倾角约 60°，区域性活动断裂—广从断裂场地外围通过场地岩土分层主要有填土层、冲洪积层、残积层及下伏石炭系灰岩、炭质页岩、砂岩基岩主要为石灰岩，且不同程度地发育有土洞及溶洞3.3.2 地球物理特征场地发现溶洞、土洞，溶洞内一般有充填物鉴于溶洞位于地下水水位以下，溶洞内充填物一般饱和含水，溶洞内无充填物的部分一般充填有水溶洞内介质的纵波波速 Vp 在 1500～2000m/s 之间溶洞外介质一般为灰岩、炭质灰岩，介质的纵波波速 Vp 在 4000～6000m/s 之间溶洞内外介质存在极为明显的波速差异根据岩溶发育的一般规律，本区可能有土洞发育鉴于土洞位于地下水水位以下，5土洞内必然充填水土洞内介质为水，纵波波速 Vp=1500m/s如土洞内介质为流塑状软土，则纵波波速 Vp=1400~1600m/s土洞外介质为饱和含水的第四系土层，纵波波速 Vp 在 1800～2000m/s 之间。

土洞内外介质存在一定的波速差异3.3.3 资料处理与成果推断跨孔弹性波 CT 成像资料的处理流程如下：抽道集→检查激发点道集→初至时间拾取→初至时间检查→射线平均波速计算→初始速度模型预测→CT 反演→波速影像图绘制→ 综合地质解释剖面图绘制成果推断的主要依据见表1：以 L02 测线后半段（即WTZK14 至 WTZK18 剖面）为例如图 3，波速影像图中红色区域（高速区）与绿色区域（低速区）之间的分界线与钻孔中的灰岩顶面具有很好的对应性由此推断，波速影像图红色与绿色之间的分界线即为灰岩顶面，零星红色小块为孤石或风化岩块在波速影像图中，红色区域（高速区）中出现非红色区域与钻孔揭露的溶洞具有很好的对应性，如 WTZK17 孔-12.67～-13.67m 段和 WTZK18 孔-9.78 ～-11.38m、-11.98～ -14.28m 和-17.38～-20.88m 段由此推断，波速影像图红色区域（高速区）中出现非红色的区域即为溶洞，绿色区域（低速区）出现的蓝色区域即为土洞各岩土层反演波速范围一览表 表1层号 岩土层 波速范围（m/s）③2 粉质粘土、含炭粉质粘土 1600～1800③3 粉质粘土 1700～2000③4 粉质粘土、含砾粉质粘土、含 碎石粉质粘土 1800～2100⑤ 强风化岩 1900～2200⑦1 微风化岩 3500～55006图 3 L02 线跨孔弹性波 CT 反演波速影像及综合地质解释剖面图3.3.4 结论本次探测基本查清了场地的岩土分层及土、溶洞等不良地质现象的分布和发育特征，较好地反映了场地的实际地质特征，为场地土、溶洞勾画出一幅较为清晰的图像。

3.4 高密度电法实例分析3.4.1 工区概况广 州 市 花 都 区 赤 坭 镇 岩 溶 塌 陷 测区处于广花盆地碳酸盐岩地区根据区域地质资料，测区位于花县复向斜褶皱构造的次一级背斜部位，地层主要由第四系冲积层（Q） ；石炭系下统测水组（C 1dc） ；石磴子组（C 1ds） ；岩关阶孟公坳组（C 1ym） 地层主要有第四系冲积层（Q）：为冲积、冲洪积相灰黄色粘土，粉质粘土，灰黑色淤泥、淤泥质土，砂土，局部夹砂层土，厚度 8～35m 不等；下伏石炭系灰岩，岩溶发育溶洞、土洞受到地下水活动影响，而产生地面塌陷等地质灾害现象3.4.2 地球物理特征第四系覆盖层一般饱和含水，电阻率较低；完整的微风化灰岩电阻率极高如果灰岩遭受强烈溶蚀作用，即发育溶洞，溶洞内充填物电阻率将明显低于周围岩体，与覆盖层电阻率相当3.4.3 资料处理与成果推断高密度电法按下列步骤进行处理：参数数据输入→视电阻率测量→数据检查→数据转换格式→绘制 ρ s曲线类型图、绘制等 ρ s断面图、用最小二乘法反演→CAD 成图→物探地质解释→剖面输出→完成成果推断可通过以下三种手段进行：a.测深类型曲线解释：电测深 ρs 曲线上的低阻异常，一般情况下预示着有低阻地质体的存在，而岩溶（包括土洞、溶洞、溶蚀裂隙等）洞穴充填水或流塑状土等物质，其视电阻率（ρs）往往明显下降，从而形成低阻异常，对应 ρs 等值线图或反演断面图呈低阻闭合圈或舌状分布，一般为岩溶发育的反映。

当然，这些低阻异常也有可能由炭质或其它不明干扰引起b. 等 ρ s 断面图解释：作为岩土分界线的基岩面，在等 ρ s 断面图解释上一般表现为界面可以连续追综、相邻等值线上下间距发生蔬、密变化明显处；岩溶则表现为7在平缓的等值线段中突然发生畸变，呈陡坎、凹陷、低阻封闭圈等特征因此，若异常变化在岩土分界线以上深度，则解释推断为土洞或软弱土层；若异常变化在岩土分界线以下深度，则解释推断为溶洞或裂隙发育如图 4 所示图 4 D1 测线等 ρ s 断面图c. 视电阻率反演断面图解释：视电阻率反演经过软件多个参数合理设置、至少 3次以上拟合生成，均方误差控制合理视电阻率反演断面图的异常一般表现为带状，带内视电阻率为低阻值，即岩溶发育如图 5 所示图 5 D1 测线。