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

1各种地球物理勘探方法在查找地下洞穴中的应用研究孙跃军（广东省地质物探工程勘察院广州510800）摘要：本文通过对各类地下洞穴探测的总结分析和研究，根据物性差异、勘查目的和要求及场地条件等方面的不同，设计出相应的地球物理勘探方法和手段，并通过工程实例进行了论证说明，进而得出地下洞穴探测的最佳方法和手段。关键词地球物理勘探采空区岩溶路面脱空防空洞方法选择1.引言在一些工程建设、矿产资源开采和环境评估中，常会遇到各种不同类型的地下洞穴，使得项目的安全和质量受到极大的影响。因此，项目的设计施工前都必须查明场地内地下洞穴的规模及分布情况。地下洞穴有天然生成的，比如溶洞、土洞等；也有人为因素而成的，比如采空区、防空洞等。地下洞穴有的存在充填物（水、流塑状物质等），有的为空洞。因此地下洞穴与周围岩土层存在明显的物性差异，可采用地球物理勘探方法进行探测。但由于场地条件、勘察目的及方法本身的局限性等方面的原因，对不同类型的地下洞穴采用的物探方法也存在差异。2.采空区的探测2.1采空区的地球物理特性存在矿产资源的地层由于成因不同而表现出不同的分布状态，多呈现为层状沉积，表现为层状展布。当矿脉被开采后形成采空区，破坏了原有的应力状态，使原有地层错动、产生裂缝、发生塌陷等。一般情况下，采空区表现为两种形式：当开采时间较短或支护措施较充分时，其矿脉层顶板为塑型岩石并保存完整，部分压力转移到地层上，未引起地层塌落、变形，采空区以不充水或充水的空洞形式保存下来；当采面较大或开采时间较长时，采空区在重力和地层应力作用下，顶板塌落，形成冒落带、裂隙带和弯曲带等，使得采空区被上层塌落的松散物充填而地球物理特征发生变化，导致塌落带上部地层的特性也发生了变化。2.2探测方法选择对于采空区的探测方法有浅层地震法、电法、瞬变电磁法等。一般来说，存在矿脉、煤田等矿产资源的地层地质构造和岩性都很复杂，开采前高低阻体就较多，采用2依靠电性差异进行勘探的电磁类物探方法对采空区的异常反映较为复杂，较难分辨真假异常，因此浅层地震反射波法才是最好的复杂地层采空区的勘探方法。地震波中含有大量的地质信息，地层构造、岩性变化、矿脉厚度变化、分岔合并等均会引起地震波的变化，其变化主要包括密度、速度及其它弹性参量的差异，这些差异导致了地震波在传播时间、振幅、相位、频率等方面的变化。当采空区未被充填时，矿脉反射波消失；当上覆地层塌落、回填，则采空区的反射波凌乱、一定范围内的反射波同相轴向下错动，呈现出松散、塌落的形态。在解释采空区时，以波形变面积显示的时间剖面为主，结合彩色剖面显示识别采空区。重点放在标准反射波同相轴的消失、变弱、强相位转换及在采空区边界处绕射波振幅的强弱变化等多种因素，来确认采空区的存在及其范围。2.3实例分析2.3.1工区概况广东省高要市长坑金矿在二十世纪九十年代进行过开采，为了评估矿产剩余储量，进行了本次物探。因此本次物探主要目的和任务是查明采空区走向及其埋深，为评估矿产剩余储量提供有力的依据。据地质勘测资料，场地岩土层主要包括：a.第四系（Q）土层：主要为洪冲积、残坡积层和填土等组成，厚度0～10m，平均约3m；b.石炭系灰岩（C1Z）：分布于全区，工区西北角出露，其它地方被覆盖，该层岩溶发育，厚度10～100m；c.三迭系泥岩、砂岩、页岩（T3xly和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之间。图1A线纵波反射时间剖面图如图2所示的L08-1线时间剖面，T1波组时间在20ms~40ms之间变化，反映了覆盖层的深浅变化，地层基本连续。在点号182～277、227~272段出现明显的弧形波组，结合地质资料可以推断182～272段存在采空区，宽度约90m，厚度在30～70m之间。图2L08-1线纵波反射时间剖面图2.3.4结论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之间。溶洞内外介质存在极为明显的波速差异。根据岩溶发育的一般规律，本区可能有土洞发育。鉴于土洞位于地下水水位以下，土洞内必然充填水。土洞内介质为水，纵波波速Vp=1500m/s。如土洞内介质为流塑状软土，则纵波波速Vp=1400~1600m/s。土洞外介质为饱和含水的第四系土层，纵波波5速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段。由此推断，波速影像图红色区域（高速区）中出现非红色的区域即为溶洞，绿色区域（低速区）出现的蓝色区域即为土洞。图3L02线跨孔弹性波CT反演波速影像及综合地质解释剖面图各岩土层反演波速范围一览表表1层号岩土层波速范围（m/s）③2粉质粘土、含炭粉质粘土1600～1800③3粉质粘土1700～2000③4粉质粘土、含砾粉质粘土、含碎石粉质粘土1800～2100⑤强风化岩1900～2200⑦1微风化岩3500～550063.3.4结论本次探测基本查清了场地的岩土分层及土、溶洞等不良地质现象的分布和发育特征，较好地反映了场地的实际地质特征，为场地土、溶洞勾画出一幅较为清晰的图像。3.4高密度电法实例分析3.4.1工区概况广州市花都区赤坭镇岩溶塌陷测区处于广花盆地碳酸盐岩地区。根据区域地质资料，测区位于花县复向斜褶皱构造的次一级背斜部位，地层主要由第四系冲积层（Q）；石炭系下统测水组（C1dc）；石磴子组（C1ds）；岩关阶孟公坳组（C1ym）。地层主要有第四系冲积层（Q）：为冲积、冲洪积相灰黄色粘土，粉质粘土，灰黑色淤泥、淤泥质土，砂土，局部夹砂层土，厚度8～35m不等；下伏石炭系灰岩，岩溶发育。溶洞、土洞受到地下水活动影响，而产生地面塌陷等地质灾害现象。3.4.2地球物理特征第四系覆盖层一般饱和含水，电阻率较低；完整的微风化灰岩电阻率极高。如果灰岩遭受强烈溶蚀作用，即发育溶洞，溶洞内充填物电阻率将明显低于周围岩体，与覆盖层电阻率相当。3.4.3资料处理与成果推断高密度电法按下列步骤进行处理：参数数据输入→视电阻率测量→数据检查→数据转换格式→绘制ρs曲线类型图、绘制等ρs断面图、用最小二乘法反演→CAD成图→物探地质解释→剖面输出→完成。成果推断可通过以下三种手段进行：a.测深类型曲线解释：电测深ρs曲线上的低阻异常，一般情况下预示着有低阻地质体的存在，而岩溶（包括土洞、溶洞、溶蚀裂隙等）洞穴充填水或流塑状土等物质，其视电阻率（ρs）往往明显下降，从而形成低阻异常，对应ρs等值线图或反演断面图呈低阻闭合圈或舌状分布，一般为岩溶发育的反映。当然，这些低阻异常也有可能由炭质或其它不明干扰引起。b.等ρs断面图解释：作为岩土分界线的基岩面，在等ρs断面图解释上一般表现为界面可以连续追综、相邻等值线上下间距发生蔬、密变化明显处；岩溶则表现为在平缓的等值线段中突然发生畸变，呈陡坎、凹陷、低阻封闭圈等特征。因此，若异常变化在岩土分界线以上深度，则解释推断为土洞或软弱土层；若异常变化在岩土分界7线以下深度，则解释推断为溶洞或裂隙发育。如图4所示。图4D1测线等ρs断面图c.视电阻率反演断面图解释：视电阻率反演经过软件多个参数合理设置、至少3次以上拟合生成，均方误差控制合理。视电阻率反演断面图的异常一般表现为带状，带内视电阻率为低阻值，即岩溶发育。如图5所示。图5D1测线ρs反演断面图3.4.4结论本次工作遵循物探勘察从点至线到面的解释推断原则，逐点、逐线进行了解释推断，根据物探异常在平面上的分布特征，最终对整个测区综合归纳解释，推断划分了66条物探异常带，根据其分布情况及岩溶发育的强弱程度，划分强、中、弱岩溶发育区。圈定的物探异常带定义为存在土洞、溶洞、裂隙发育、塌陷的集合，属不良地质体较密集的区域。83.4管波探测法实例分析3.4.1管波探测法简介管波探测法由广东省地质物探工程勘察院工程技术人员在长期的工程实践中发明，2007年荣获广东省科技进步二等奖，被广东省建设厅列入“广东省建设