《武汉工程大学学报》 2015年10期
11-15
出版日期：2015-10-31
ISSN:1674-2869
CN:42-1779/TQ
高密度电阻率法在矿山岩溶探测中的应用
０ 引 言青山河位于湖北省大冶市铜录山矿区西侧，该河中下游段与该矿南露天坑的距离约５００ ｍ，河床流域内的岩溶发育带与矿区相连. 矿山的长期疏干排水导致青山河流域岩溶塌陷，洪水期河水倒灌流入矿坑. 矿山曾多次采用堤坝加固、回填封堵塌陷区等方法治理，但收效甚微. 地下水长期对河床底部冲刷和掏蚀，导致每逢汛期河床便再次复活塌陷，河水和大冶湖水通过复活的塌坑、岩溶裂隙流入矿坑，严重威胁矿山安全生产. 为了彻底解决矿山水患问题，采用高密度电阻率法探测了解河床底部的岩溶发育特征及矿区的地下水主径流通道，并根据探测成果对河床底部岩溶进行了注浆加固试验，有效地封堵了小清河河水及矿区地下水流入矿坑的主径流通道，消除了河床岩溶塌陷复活、洪水倒灌对矿山安全生产的威胁. １ 高密度电阻率法的基本原理及特征１．１ 高密度电阻率法的基本原理２０世纪７０年代末期英国学者Ｊｏｈａｎｓｓｏｎ提出了高密度电阻率法的最初模式（电测深偏置系统）［１］，８０年代中期日本地质计测株式会社借助电极转换板转换电极排列方式实现了野外高密度电阻率法的数据采集，８０年代后期我国原地质矿产部开始了高密度电阻率法及其应用技术研究［２］. 高密度电阻率法的基本原理与常规直流电阻率法相同，以探测目标体与周围物质导电性差异为基础，在人工电流场的作用下，通过观察和分析视电阻率的变化规律来分析地质问题［３］. 高密度电阻率法可以在探测断面上同时布置多个电极，通过电极转换装置控制电极的不同排列组合（如α、β、偶极、施伦贝谢和温纳－施伦贝谢等［４］）可以同时完成测剖面和测深工作，可以提供更加丰富的地断面信息，有利于分析和比较，见图１，其中ａ为极间距. 图１ 高密度电阻率法电极排列示意图Ｆｉｇ．１ Ｃｏｍｍｏｎ ａｒｒａｙｓ ｕｓｅｄ ｆｏｒ ｈｉｇｈ ｄｅｎｓｉｔｙ ｒｅｓｉｓｔｉｖｉｔｙ ｍｅｔｈｏｄ１．２ 高密度电阻率法的特点高密度电阻率法具有测点密度高，获取信息量大，对探测对象无损伤，成果精确、工作效率高等特点［５］. 与常规电阻率法相比具有如下优点：①电极布设一次性完成，有助于减少因电极设置引起的干扰和故障，减小测量误差；②通过转换电极排列方式可获得丰富的地电断面信息；③采用自动化的数据采集方式，不仅可以提高数据采集效率，还可以减小人工操作造成的误差；④可以对资料进行预处理并显示剖面曲线形态. ２ 测区概况本次高密度电法探测剖面布设在铜录山青山河两岸和露采坑西南侧半坡路面. 青山河河床及其两侧为第四系冲积亚粘土，细—中砂、砂砾石岩层. 结构松散，稳定性较差. 其下为三叠系碳酸岩地层，岩溶发育，富水性强. 近地表浅部发育有溶洞、开口型的洞穴及溶沟. 矿床开采过程中不断疏排地下水，引起区域地下水位下降，溶洞、洞穴上部第四系土层失去了地下水的“浮托力”，同时在地下水潜蚀、掏空、搬运第四系土层及溶洞、洞穴充填物的作用下，自然条件下的相对平衡遭到了破坏，导致河床出现下沉、开裂、塌陷现象. 为了解矿区地下水流入矿坑的主径流通道，采用了高密度电阻率法探测矿区地下岩溶分布特征及规律. ３ 工作方法与技术本次外业资料采集采用ＷＧＭＤ－９超级高密度电阻率法系统，该系统具有抗干扰性强、自动实现自然电位及电极极化补偿、可随时检查接地电阻等优点. 经现场初步试验后，决定在青山河岸Ⅰ、Ⅱ号剖面采用温纳装置，露采坑南侧 Ⅲ号剖面综合采用温纳装置、斯伦贝谢装置. 测线参数如表１所示，Ⅰ、Ⅱ、Ⅲ剖面位置如图２、３所示. 表1 测线参数表Table 1 Survey line parameters图２ 物探测试Ⅰ、Ⅱ号剖面位置图Ｆｉｇ．２ Ｐｏｓｉｔｉｏｎ ｏｆ ｐｒｏｆｉｌｅ １ ａｎｄ ２ ｂｙ ｇｅｏｐｈｙｓｉｃａｌ ｔｅｓｔ４ 测试成果解释外业采集的数据经过圆滑处理后，输入计算机进行二维、三维反演后获得地电断面信息，数据圆滑处理一般采用坏点剔除和滑动平均等方法［６］. 二维反演计算采用基于准牛顿最优化非线性最小二乘新算法，其计算速度可以达到常规最小二乘法的１０倍以上. 反演结果主要通过均方误差（ＲＭＳ）来衡量［７］，ＲＭＳ越大，反演结果的可靠度越低. 通过多次迭代可以减小均方误差，为保证反演结果的可靠度，本次探测成果采用五次迭代反演. ４．１ 青山河东岸Ⅰ号剖面Ⅰ号剖面位于青山河东岸，外业采集的两条温纳装置数据剖面反演成果一致性良好，实测数据可以客观反映测区的地电条件，数据成果可靠. 高、低阻区在断面上有清楚的反映，在断面上存在１＃、２＃两个低阻异常体（见图４），根据地质资料综合推测两处为强岩溶发育地段. １＃异常体平面位置位于测点号５３～９６ ｍ处，上顶埋深约７ ｍ，下底埋深约２５ ｍ，呈漏斗状，推断为强岩溶发育区. ２＃异常体平面位置位于测点号１１５～１５０ ｍ处，与现场塌方处一致，异常上顶埋深２８ ｍ，由于探测深度的原因，异常下部没有闭合，下底埋深大于４０ ｍ，推断为强岩溶发育区.图３ 物探测试Ⅲ号剖面位置图Ｆｉｇ．３ Ｐｏｓｉｔｉｏｎ ｏｆ ｐｒｏｆｉｌｅ ３ ｂｙ ｇｅｏｐｈｙｓｉｃａｌ ｔｅｓｔ图４ Ⅰ号剖面地质解释图Ｆｉｇ．４ Ｇｅｏｌｏｇｉｃａｌ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｏｆ ｐｒｏｆｉｌｅ １４．２ 青山河岸Ⅱ号剖面Ⅱ号剖面位于青山河西岸，平面位置与Ⅰ号剖面隔河对称，外业采集的两条温纳装置数据剖面反演成果一致性良好，实测数据可以客观反映测区的地电条件，数据成果可靠.高、低阻区在断面上有清楚的反映，在断面上存在３＃、４＃ 两个物探异常体（见图５），根据相关地质资料综合分析推断这两处为岩溶发育地带. ３＃异常体平面位置位于测点号５３～９８ ｍ处，上顶埋深约９ ｍ，下底埋深约２２ ｍ，推断为强岩溶发育区. ４＃异常体平面位置位于测点１３０～１７２ ｍ处，与现场塌方处一致，异常上顶埋深５ ｍ，下底埋深１７ ｍ，推断为强岩溶发育区.图５ Ⅱ号剖面地质解释图Ｆｉｇ．５ Ｇｅｏｌｏｇｉｃａｌ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｏｆ ｐｒｏｆｉｌｅ ２图６ Ⅲ号剖面地质解释图Ｆｉｇ．６ Ｇｅｏｌｏｇｉｃａｌ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｏｆ ｐｒｏｆｉｌｅ ３４．３ Ⅰ、Ⅱ号剖面对比分析Ⅰ、Ⅱ号剖面所反映的地层吻合度较好. 都存在两个高阻区即两个完整性较好的灰岩岩体，其中，Ⅰ号剖面完整性较好的灰岩岩体被埋深较大的岩溶隔开，Ⅱ号剖面因岩溶埋深较浅，探测成果反映两个完整性较好的灰岩岩体连成一片. １＃、３＃岩溶区的形态、位置、埋深基本一致，推断为同一岩溶发育区. ２＃、４＃ 岩溶发育区分别与现场塌方处一致，２＃岩溶发育区埋深大，４＃岩溶发育区埋深小，２＃岩溶发育区位置较４＃ 岩溶发育区偏北，推断为同一岩溶发育区. 该岩溶发育区走向北北东，沿南南西向下延伸，延伸深度大于４０ ｍ，形成的空洞较大、较深，这正是此处塌方的原因之所在.４．４ 露采坑南侧 Ⅲ号剖面Ⅲ号剖面位于露采坑南侧半坡路面，外业采集了两条温纳装置的数据剖面和一条斯伦贝尔装置的数据剖面，三个剖面反演成果一致性良好，实测数据客观地反映了测区的地电条件，数据成果可靠.Ⅲ号剖面存在５＃、６＃ 两个物探异常体（见图６）. 根据相关地质资料，剖面所处地层：测点０～４１０ ｍ为中下三叠统白色中粗粒大理岩，测点４１０～５９０ ｍ为石英正长闪长玢岩，测点４１０ ｍ处附近为岩性接触带，接触带大角度倾向东南. 因此，接触带与Ⅲ号剖面在测点４１０ ｍ处附近呈近似正交且大角度倾向大测点号方向.５＃异常体平面位置位于测点号２６０～３２５ ｍ，上顶埋深约１８ ｍ，下底埋深大于９０ ｍ，异常平面呈垂向延伸、陡立状，推断为金属矿化低阻体或裂隙发育带. ６＃异常体平面位置位于测点３６５～４５０ ｍ，上顶埋深１７ ｍ，下底埋深７０ ｍ，异常平面近似等轴状. 该异常为复合异常体，异常的主体为矿坑陡立金属排水管群，异常右边界为接触带右边界的反映，物探成果显示的岩性接触带位置与现场地表出露的岩性接触带走向基本吻合，接触带呈大角度倾向剖面大测点号方向. ５ 结 语高密度电法重复观测成果资料一致性良好，实测数据客观地反映了测区的地电条件，数据成果可靠. 后期对青山河沿岸物探异常体进行钻孔验证，发现实际地质情况与物探结果基本相符. 通过对河床岩溶塌陷进行注浆充填，封堵了小清河岩溶塌陷区与矿坑之间的主径流通道，基本消除了丰水期洪水通过河床塌陷灌入矿区矿坑的风险. 本案例说明先通过高密度电阻率法探测出岩溶分布特征及规律，然后采用帷幕注浆技术对矿区地下水流入矿坑的主径流通道进行重点封堵， 不仅可以大幅降低治水成本，而且可以达到较好的堵水效果. 该方法在类似岩溶矿山的地下水防治中可以推广使用.致 谢感谢武汉工程大学对本研究的资助！