高密度电法与瞬变电磁法在戈壁区找水的联合应用_唐塑.pdf

唐塑,武银婷,邢浩,等.高密度电法与瞬变电磁法在戈壁区找水的联合应用J.CT 理论与应用研究,2023,32(1):27-34.DOI:10.15953/j.ctta.2022.081.TANG S,WU Y T,XING H,et al.Combined Application of High-density Electrical Method and Transient Electromagnetic Method inGobi Desert AreaJ.CT Theory and Applications,2023,32(1):27-34.DOI:10.15953/j.ctta.2022.081.(in Chinese).高密度电法与瞬变电磁法在戈壁区找水的联合应用唐塑1,2,武银婷2,邢浩1,娄雪聪1,石海涛11.中国地质调查局乌鲁木齐自然资源综合调查中心,乌鲁木齐 8300572.长安大学地质工程与测绘学院,西安 710061摘要：为准确掌握地下水资源的动态变化，有效支撑生态环境的可持续发展，在新疆哈密某矿集区利用高密度电法与瞬变电磁法联合进行物探勘查，对地下含水层与隔水层进行预判。经验证，高密度电法与瞬变电磁法的组合，有效指导水文钻探施工，为同类地质条件下物探手段的选取和应用提供参考。关键词：高密度电法；瞬变电磁法；荒漠戈壁区；地下水DOI:10.15953/j.ctta.2022.081中图分类号：P 631.3文献标识码：A荒漠戈壁区分布着丰富的矿产资源，但其生态环境极为脆弱，水资源尤为珍贵。随着社会经济发展对矿产资源的需求加大，矿业活动日渐频繁。目前我国荒漠戈壁区及其地下水资源的破坏加剧，在矿山建设水文观测基础设施、掌握水文动态变化数据，是地下水资源科学管控和维护生态可持续发展的有效手段。布设水文钻孔前，地球物理勘查在预判找水中发挥着重要的作用，但每种地球物理方法都有其适用条件以及局限性，单一方法的使用存在一定的多解性，所以采用综合地球物理方法找水，不仅可以多角度分析研究地下岩层性质和含水性，又可以相互验证方法的有效性1。传统地下水勘查工作中采用最多的是直流电阻率法和激发极化法2。2010 年李国占等3采用音频大地电磁法和激发极化法相结合的方式，在河北保定寻找裂隙构造水工作中取得了满意的效果；2011 年王星明等4采用对称四极电测深法、瞬变电磁法、高密度电阻率法及激发极化法等技术在陕西榆林，开发寻找永久性大型或特大型供水水源地方面解释精度高、成果可靠；2012 年宋希利等5在山东沂蒙山侵入岩地区，采用可控源音频大地电磁法和电阻率测深寻找地下水，有效地提高了效率和定井成功率；2015 年苏永军等6在河南新乡水文地质条件复杂的松散层和基岩干旱地区，采用高密度电阻率法和激发极化法组合模式勘查地下水，提高了地质解释的准确性和可靠性；2017 年孙中任等7综合利用磁法勘探、瞬变电磁法、可控源音频大地电磁法和大地电磁法，在辽宁某地热带实现了深部地下水的勘查；2020 年康方平等8利用视电阻率联合剖面法、高密度电阻率法和激电测深法组合，在湖南邵东板岩地区有效查明了区内赋水相对有利的断层构造等部位；2021 年潘剑伟等9在湖北白水河地区，综合利用地面核磁共振法、自然电场法确定地下水的空间分布和运动规律，为滑坡地下水模型的建立提供了丰富的水文地质信息。本文采用高密度电法和瞬变电磁法组合手段，在典型的荒漠戈壁区，开展综合地球物理方法找水研究与应用，为相似地质条件水文勘查施工中物探方法的选取提供实践参考。1工作区地质概况及地球物理特征工作区位于新疆哈密市南部某矿集区。区内出露的主要地层有：侏罗系中统西山窑组上段，泥收稿日期：20220422。基金项目：新疆哈密市南湖一带矿山集中区生态修复支撑调查项目（DD20208081）。第 32 卷第 1 期CT 理论与应用研究Vol.32,No.12023 年 1 月（2734）CT Theory and ApplicationsJan.,2023岩、粉砂岩、砾岩互层，夹炭质泥岩及不稳定的薄煤层；新近系桃树园组，泥岩、砂质泥岩、泥质砂岩夹砾岩及灰岩；第四系，黄土、砂质黏土、砾石、细砂、砂砾层、风成砂土；局部发育少量石炭系侵入岩。本区断裂构造不发育，未发现断层存在，构造类型属简单构造类型，全区大部分地层倾角平缓，一般为 515。所在区域自然地理条件为极度干旱，大气降水奇缺，无地表径流及水体，地下水主要补给来源为北部山区的冰雪融水。本区主要岩性由第四系松散岩类、新近系及侏罗系沉积碎屑岩类组成，找水目的层主要是新近系碎屑岩类孔隙裂隙含水岩组。其中沉积碎屑岩的各类岩石，其单层厚度沿走向方向的变化较大，可由几厘米变化到数米，尤其以砂岩最为明显，因此，只能以较大的岩性段划分含（隔）水层（段），其中隔水层主要岩性以泥岩、泥质砂岩为主，而含水层岩性主要以粗砂岩、砾岩为主。据工作区大量钻孔资料：泥岩、泥质砂岩电阻率约为 515 m，粗砂岩、砾岩电阻率约为 1040 m。本区的隔水层主要以相对低阻的泥岩、泥质砂岩为主，含水层主要以相对高阻的粗砂岩、砾岩为主。利用两者之间电性存在差异，综合采用高密度电法和瞬变电磁法进行测量。一方面圈出中高阻的粗砂岩、砾岩地层及不同岩相显著变化区空间分布特征，然后在中高阻的粗砂或砾岩地层中找到相对低阻的裂隙发育位置，最终找到含水裂隙带，实现找水目标10。另一方面，通过分析研究区域周边，特别是底部，找到低阻的泥岩、泥质砂岩地层，利用寻找底部隔水层实现找水目的11。2联合探测方法及工作布置2.1理论方法高密度电阻率法的基本理论与传统的电阻率法完全相同，由供电极 A、B 往地下供电建立电场，测定 M、N 电极之间的电位差，进而计算出视电阻率值，所不同的是高密度电法在观测中设置了较高密度的测点，现场测量时，只需将全部电极布置在一定间隔的测点上，由主机自动控制供电电极和接收电极的变化，即沿着测线不停的变换电极间距，逐点测量电流、电位差，然后通过下列公式，计算视电阻率。高密度电法测量系统采用先进的自动控制理论和大规模集成电路，使用的电极数量多，而且电极之间可自由组合，这样就可以提取更多的地电信息12-15。高密度电法主要是通过地下介质间的导电性差异，查明含水构造、裂隙和破碎带等空间分布来实现找水的目的16。s=KVI，（1）K=21AM1AN1BM+1BN，（2）sVIK式中，为视电阻率，为观测的电位差，为供电电流，为装置系数。另一方面，瞬变电磁法是一种时间域的电磁探测方法。介质在一次电流脉冲场激励下会产生涡流，在脉冲间断期间涡流不会立即消失，在其周围空间形成随时间衰减的二次磁场。二次磁场随时间衰减的规律主要取决于异常体的导电性、体积规模和埋深，以及发射电流的形态和频率17-21。地质体导电性越好，涡流的热耗损越小，瞬变的过程就越长，因此，我们可以通过接收线圈测量的二次场空间分布形态，了解异常体的空间分布。瞬变电磁法野外采集的数据就是二次场的感应电动势，通过公式（3）计算出视电阻率22。其基本原理和常规电法大同小异，仍然是以地质异常体与周围岩体之间的电性差异为基础来研究含水构造、裂隙和破碎带等空间分布。=04t20STSR5t(V(t)/I)23，（3）00=4107H/mSTSRV(t)V(t)/I式中，为视电阻率，为真空磁导率：，为发送线圈面积，为接收线圈面积，为感应电动势，为瞬变值。28CT 理论与应用研究32 卷2.2工作布置为便于观测矿业活动对地下水的影响，同时降低高密度电法接地电阻，物探测线部署在工作区某煤电一体化企业附近，一条地势平坦的干涸河床内（图 1）。河床内地表土壤均已板结龟裂，有部分植被生长，将板结的土块刨开后插入电极，接地电阻相对较小。河床两岸均为无植被覆盖的荒漠戈壁，相较于河床内，接地电阻极高，不适合布设开展高密度电法工作。因此，在干涸河床内布置两条测线，WT1 线与 WT2 线，均进行高密度电法和瞬变电磁法测量。本次高密度电阻率法测量使用 AGI 公司生产的高密度电法仪 SuperSting R8，采用温纳装置，电极距为 20 m。瞬变电磁法测量使用 Phoenix Geophysics 公司生产的 V8 多功能电法仪，采用中心回线装置，发射线圈为边长 240 m 的正方形，接收装置的等效面积为 100 m2，点距为 40 m。点号布设由北向南、由西向东依次增大，点号大小与该测点的实际距离值相等。3数据处理及异常解释3.1数据处理不论是高密度电法还是瞬变电磁法，数据处理中最重要的都是反演，通常采用一个简单的模型近似模拟真实的地质状况，然后尝试从观测数据中确定构造的真实地电参数，二维反演同一维反演的数学性质相同，都是要寻找一个地电模型，使其对应的理论计算值与实测视电阻率值在一定规则下拟合最好23。阻尼最小二乘法是高密度电法与瞬变电磁法在反演中最常用的方法之一，即求解目标方程（4）的最小值问题。P(m,d)=?WdA(m)Wdd?2+?WmmWmmapr?2，（4）P(m,d)=min，（5）mn+1=mn+m，（6）Qp3X CJ2XQhplWT1线/000WT1线/500WT2线/000WT2线/700zk03Qp3X CJ2X图例全新统洪积物上更新统新疆群西山窑组砾岩、砂泥岩石炭系石英闪长岩WT1线/5zk03整合界线不整合界线物探点位 及编号钻孔点位 及编号Qhpl图 1工作布设示意图Fig.1Sketch map of work layout1 期唐塑等：高密度电法与瞬变电磁法在戈壁区找水的联合应用29P(m,d)mmaprdA(m)mWdWmm式中，表示目标方程，为需要求解的地电模型，为一先验模型，为野外获得的观测数据，为由地电模型 算出的观测数据，和 分别为数据与模型权矩阵，为阻尼因子，为每次迭代计算的改变值。最小二乘法反演的计算过程实际就是通过迭代计算方程（6）满足方程（5）的条件，目标方程（4）的物理意义就是计算出的观测值与实际观测值的差达到最小，同时计算出的地电模型与提前给出的先验模型的差达到最小。从中可以得出，反演结果的质量与两个方面息息相关，一是准确无误的实际观数据，二是符合实际的先验模型，所以数据处理的过程要以这两点为准则。高密度电法数据处理软件为仪器厂家 AGI 配套的 EarthImager 2D，可以对数据进行预处理，删除和修正数据中的负值、孤立的高值等，然后可以对数据进行二维反演。瞬变电磁法数据处理软件为 Steminv，首先将原始数据利用程序 TEMHub 转换格式，然后将数据导入软件 Steminv，在软件中可以查看是否需要进行噪声剔除和曲线圆滑，预处理后就开始对数据进行一维反演。得到高密度电法和瞬变电磁法反演数据后，利用软件 Surfer 生成图件，对比两种物探方法的反演断面图，寻找对应的异常特征。3.2高密度电法分析图 2（a）和图 3（a）分别是 WT1 线、WT2 线高密度电法反演断面图，断面横向长度分别为500 m 和 700 m，纵向深度范围为 0-150 m，剖面整体地电信息均呈现水平层状分布，电阻率由浅部至深部表现为“低阻-高阻-低阻”特征。WT1 线断面中，深度-40-110 m 范围内，电阻率呈相对高阻，电阻率值大于 10 m，推断为砂岩、砾岩层；在点号 200 至 500 范围内等值线波动较大，有相对低阻异常反映，等值线形态呈现类似“U”型，其电阻率低于 20 m，根据本区地层情况，推断其为基岩裂隙含水区。WT2 线断面中，深度-40-90 m 范围内，电阻率呈相对高阻，电阻率值大于 10 m，推断为砂岩、砾岩层；在点号300400 范围内等值线形态呈现类似“U”型，有相对低阻异常反映，其电阻率低于 20 m，根据本区地层情况，推断其为基岩裂隙含水区。而且 WT1 线与 WT2 线推测的基岩裂隙含水区基本重合，均在两条剖面交点处附近。视电阻率/（m）10002003004005001500100150322722171272（a）10101818（a）高密度电法反演50100150200视电阻