孔间电磁波CT技术在北京普安店岩溶勘查中的应用

技术应用46472014第9卷第3期2014Vo1.9No.32014Vo1.9No.32014第9卷第3期参考文献[1]杜华光.高密度电法在工程地质勘察中的应用[J].路基工程,2008，139(4):110～112.[2]张献民.应用高密度电法探测煤田陷落柱[J].物探与化探,1994(5):363～370.[3]肖宏跃主编.地电学教程[M].地质出版社,2008.9.[4]刘海飞.高密度电阻率法数据处理方法研究[J].中南大学,2004.[5]马德锡,于爱军.高密度电法在金矿勘查中的应用[J].地质与勘探,2008,44(3):65～69.[6]秦正.高密度电阻率法在工程勘察中的应用[J].物探装备,2005,15(3):205～209.[7]王兴泰.高密度电法在煤矿采空区探测中的应用[J].能源技术,2006,12(5):011～013.TheApplicationoftheHigh-densityElectricalMethodinDetectingtheUndergroundFaultsCHIYuelong1.2,ZHANGTing1.2,WANGXingjie1.2(1.SchoolofEngineeringTechnology，ChengduUniversityofTechnology，Leshan614000；KeyLabEarthExploration&InformationTechniquesofMinistryofEducation,ChengduUniversityofTechnology,Chengdu610059)Abstract:Thispaperintroducestheworkingprincipleandmethodofthehigh-densityelectricalmethod,andappliesthemethodtosurveytheundergroundfaultsofLongmenshanfaultarea.Thesurveyresultsshowthathigh-densityelectricalmethodisbeneficialtodetecttheundergroundfaults.Keywords:High-densityelectricalmethod;undergroundfaults;Geodeticelectricaldetection本刊讯（孙宇浩段金平段海明）在日前召开的“西藏地热专题技术论坛”上，北京市地热协会、北京地质学会发出倡议，北京地区乃至全国广大地热科技工作者和企事业单位，要以饱满的热情、严谨的科学态度，积极关注或投身西藏地热勘查开发事业中。来自中国地调局、中国煤炭地质总局、中科院青藏高原研究所、北京市国土资源局、北京市地勘局、西藏地热（发电）工程研究中心，以及中国矿业联合会、北京大学等单位的50余名专家学者在倡议书上签名。西藏得天独厚的自然地理位置，孕育了丰富的地热资源，是我国中高温地热资源最为集中的地区。截至2012年，西藏共发现地热爆炸、间歇喷泉、热泉等各种水热活动显示区700多处，其中可供开发的地热显示区342处。在已调查的169个热田和水热区中，温度高于80℃的占22%，温度在60℃～80℃之间的占26%，40℃～60℃的占35%，低于40℃的占17%。《西藏自治区地热资源区划》将西藏地热资源进行了三级区划：5个资源区、19个资源带、75个资源地。北京市地热协会、北京地质学会发出倡议积极投身西藏地热勘查开发事业0引言孔间电磁波CT是用发射电磁波的方法对井间的地质剖面进行层析技术处理，得到井间地质剖面的内部结构图像，该项技术分辨率高、野外施工方便，可以避开地面建筑。同时该方法充分利用已有钻孔，填补钻孔之间地质信息空白。相对地面物探，其受到地面建筑及电磁干扰较小，目前广泛应用于断层破碎带、软弱夹层、岩溶等地质构造方面的探测，如井间电磁波技术在碳酸盐岩缝洞勘察中的应用[1]、在居民区地下洞穴调查中的应用[2],在水库堤坝隐患探测中的应用[3].1研究区概况北京香山普安店地区在近几年发生了多起地面塌陷灾害，其中较严重的有两次，分别发生于2009年7月25日和2010年7月3日。前者形成直径约4m，深约5m的塌陷坑，造成一名人员掉入塌陷坑内、一棵大树陷入坑内及一条污水管道受损折断；后者形成塌陷坑大小约2.5m×2.5m，深3～4m，造成了2名人员受伤和部分财产损失。地面塌陷给当地居民造成恐慌，成为该地区的孔间电磁波CT技术在北京普安店岩溶勘查中的应用李远强（北京市地质研究所，北京100120）摘要：孔间电磁波CT技术是工程勘探中的一种常用方法,利用孔间层析成像进行探测，其成果清晰、直观,较好地展现了孔间介质电磁波吸收情况,清楚地显示出地下不同地质体的空间分布。北京普安店地区存在岩溶塌陷，前期进行了大量的钻孔勘探，利用钻孔进行孔间电磁波CT扫描，探测钻孔之间溶洞的分布情况，填补了钻孔之间溶洞展布情况的空白，并对溶洞的充填情况进行解释。该技术提供的数据准确可靠，信息丰富，为后期岩溶地区的稳定性评价及工程处理提供了科学依据和基础数据。关键词：孔间电磁波CT；岩溶勘查；数据采集处理中图分类号：P631.3+25文献标识码：A文章编号：1007-1903（2014）03-0047-06不安定因素和安全隐患，当地政府为民生安定，开展了该地区的地面塌陷灾害勘查工作。研究区位于北京西山山前丘陵向平原过渡的区域，地势平坦，地面标高57～59m，地形坡度在1°～2°之间。属于暖温带半湿润大陆性季风气候，夏季极端最高气温35℃～40℃之间，冬季极端最低气温－11℃～－19℃；多年平均降水量613mm，年降水量最多达1406mm，夏季（6～8月）降水集中，占年降水量的75％，降水强度大，雨量往往集中在几次降雨过程中，地面塌陷多发生于该时间段。研究区基岩被第四系（Q4）所覆盖，为坡积、冲积、冲洪积、洪积物，以粘土、粉土、黄土夹砂砾石及砂砾石为主；基岩主要为奥陶系马家沟组灰岩（Om），均隐伏于第四系之下，基岩埋深最浅处小于10m。区内发育有近NS向冷泉村平移断裂（象鼻子沟断裂），断面倾向西，倾角80°～90°，东盘北移，西盘南移。钻孔揭示研究区灰岩为灰色，强-中等风化，中-厚层状，层面倾角约30°～40°，节理裂隙发育，岩芯上可见溶孔及溶槽，溶蚀现象较发育，上部岩芯采取率低，漏浆；多个钻孔揭示地下发育的溶洞，埋深17～50m，高度约0.1～3.7m不等，多数溶洞充填碎石、粘土等，李远强等：孔间电磁波CT技术在北京普安店岩溶勘查中的应用技术应用48492014第9卷第3期2014Vo1.9No.32014Vo1.9No.32014第9卷第3期少数无充填。钻孔揭露的土体自上而下为杂填土、素填土、粘质粉土、粉质粘土等，下部含有碎石、卵石；钻孔揭示地下发育的多个土洞，埋深0.6～37m，高度约0.8～5.3m不等，少数土洞充填碎石土，多数无充填。2孔间电磁波CT的基本原理孔间电磁波CT的基本原理是在两个钻孔分别发射和接收无线电波(工作频率0.5~32MHz),根据不同位置上接收的场强大小,来确定地下不同介质分布的一种地下地球物理勘查方法[4](图1)。由于地下介质具有不同物理性质，在电磁波的传播过程中主要表现在对电磁波能量的吸收,使其能力衰减，这种吸收作用与地下岩土体的类别、裂隙分布、含水程度、充填物质等因素有关,通过2个钻孔之间电磁波扫描性观测,利用层析成像反演算法,将不同岩性导致的电磁波能量上的差异分布转变成二维分布图像,进而推断出地下地层及构造情况。该方法具有电磁波在地下有耗半空间的辐射、传播和接收的特征，其正反演问题的理论基础是电磁场理论。下式为地下电磁波法中的场强观测值公式：E=E0′e-βγf(θ)r–1…………………………(1)式中：E为接收点的场强值；E0′为初始辐射常数；β为测区介质吸收系数,即介质中单位距离对电磁波的吸收值；f(θ)为收发天线的方向因子函数；r为发射与接收点之间的距离。式(1)表明通过e-βγf(θ)r–1因子，E0′衰减到E,吸收系数（β）是一个与介质电阻率（ρ）、介电常数（ε）、磁导率（μ）以及电磁波频率数（ω）有关介质的主要参数，它表示介质对电磁波的吸收特性,当介电常数（ε）和磁导率（μ）一定时,吸收系数（β）主要与介质电阻率（ρ）有关。一般地,电阻率（ρ）越高，吸收系数（β）就越小,即介质的性状愈好；反之，即介质的性状越差,电阻率（ρ）越低，吸收系数（β）就越大[2]。由此可见,强度高、坚硬完整、较纯的灰岩中,地球物理特征常表现为高速、高阻、低吸收的特征，而当岩层中出现断裂破碎带、溶洞、泥沙充填时，则表现为：波速降低、电阻率降低,吸收系数增大的现象，与完整灰岩间存在较大的地球物理差异。介质吸收系数的大小表示着不同的岩性和同一岩性的结构完整性。图1孔间电磁波CT原理示意图3野外数据采集和处理在野外钻孔施工完成后需要安放PVC套管保护钻孔，以保证探测顺利进行，避免岩溶发育段塌孔现象发生，造成损失。现场使用湖南岳阳奥成科技有限公司生产的HX-JDT-02孔间电磁波透视仪，该设备工作频率范围宽、动态范围大、功耗极低、主控机极轻便、操作简单快捷、稳定可靠等特点。3.1探测实验探测之前，先进行仪器参数测试，首先在空气中，由发射探头发射，接受探头接收，随着距离的增加场强衰减如表1。从表中可以看出在空气中，频率为5MHz时，场强衰减最慢。其次，选取钻孔间距最小的ZK07和ZK08（间距23m）钻孔进行土层吸收系数测试，取得接收机在无发射和有发射时电磁波场强值数据，结果见表2。从表中可以看出：在土层中，由于吸收系数过大，接收机接收到的信号与发射机不发射时接受的信号差异太小，无法进行后期解释。所以，在土层中无法应用电磁波CT进行探测。李远强等：孔间电磁波CT技术在北京普安店岩溶勘查中的应用1m5m10m20m35m5MHz-53.8-75.2-77.1-90.4-90.36MHz-47.8-66.3-79.7-92.8-95.97MHz-58.1-64.5-83.4-86.2-108.18MHz-51.7-60.2-76.7-81.8-99.69MHz-43.2-57.4-81.1-92.9-101.710MHz-46.1-58.2-72.2-83.4-90.811MHz-46.1-58.6-75.6-81.5-111.212MHz-36.4-55.2-59.8-71.7-97.6表1空气中不同频率、不同距离场信号场强值一览表距离场强(dB)最后，选取钻孔间距33m的ZK04和ZK05钻孔进行岩层吸收系数实验测试，测试结果表明：5～8MHz在岩石完整孔段观测到的场强多为-45～-80dB，9～12MHz的场强多为-33～-100dB；在较完整灰岩的孔段，将测试到的场强可视为正常场强值，用“两孔测定法”确定初始场强E0′，并由此来计算较完整灰岩吸收系数，经式（2）。…………………………（2）计算得到的视吸收系数βs多在0.01～0.08db/m，该数值用于成果解释中对异常进行划分。在以上详细实验的基础上，确定了工作频率、扫描间隔、点距、空气中的背景场值和较完整灰岩的吸收系数等探测参数，见表3。表3测试实验结果参数一览表3.2数据采集及处理现场采用定点观测法进行数据采集，将发射机和接收机分别放置在两个钻孔之中，把发射机固定在钻孔某一深度上不动而移动接收机进行，完成一组数据采集后，移动发射机点位，采集下一组数据，如此循环，直至所有数据采集完成；如果将接收机固定在钻孔某一深度上不动而移动发射机，其数据采集结果与前者相同。电磁波CT资料的处理采用了阻尼代数重建法(ART)，其基本思想是先给被重建区域一个初始值，然后将所得到的投影残差逐个沿其射线方向均匀地反投影回去，并不断地对重建图像进行修正，直到满足要求为止。该层析处理方法采用加入阻尼技术来控制约束其结果，以克服迭代时容易发散的缺点，该方法占用的内存少且计算速度快。首先对原始数据进行平滑滤波及人机联