综合电法在岩溶山区地下水勘探中的应用——以湖南怀化长塘村为例

邬健强，赵茹玥，甘伏平，张伟，刘永亮，朱超强

(1.中国地质科学院 岩溶地质研究所，广西 桂林 541004; 2.自然资源部广西岩溶动力学重点实验室，广西 桂林 541004; 3.武汉市测绘研究院，湖北 武汉 430022)

0 引言

岩溶山区碳酸盐岩溶蚀作用强烈，水文地质条件复杂，常形成孔隙、裂隙、断(层)裂破碎带等，而地表水能够通过这些通道渗漏于地下，导致地表干旱缺水，严重威胁了粮食的生产和居民的生活。岩溶山区地下水的勘探主要是以物探为常用手段[1-11]，焦彦杰等[12]综合利用电导率成像系统和高密度电法圈定了含水岩溶带并确定了井位，有效解决了当地饮水问题；陈玉玲等[13]通过不同物探方法的对比，分析了岩溶山区可控源音频大地电磁法、三级电测深法和充电法在探测地下暗河的优缺点；张银松等[14]利用高密度电法确定了岩溶区异常分布范围，接着利用联合剖面法进一步确定异常的具体位置。

文中以地处岩溶山区的湖南怀化长塘村为例，为了解决当地的灌溉水和饮用水问题，综合采用高密度电阻率法和音频大地电磁法进行了优化组合找水应用研究。通过对两种方法所得结果的综合讨论与分析，相互验证，确定井位并成功出水，最大涌水量大于146.88 m3/d。

1 测区概况

1.1 测区水文地质特征

工作区位于湖南怀化市长塘村，地形地貌上属于构造溶蚀低丘谷区，谷地开阔，宽3～7 km，残留缓坡孤丘高差10～30 m，地表均被黏土夹碎石覆盖，未见溶蚀漏斗，地下岩溶管道发育，是岩溶水的汇集、排泄区。测区出露主要地层为石炭系中统黄龙群组(C2hn)与震旦系下统南沱冰碛岩组(Zant)、黄龙群组(C2hn)岩性为灰、灰白色厚层灰岩，白云质灰岩及白云岩，地下水类型为碳酸盐岩裂隙岩溶水，为测区的含水层，南沱冰碛岩组(Zant)岩性为含砾粉砂岩、冰啧砾岩、泥砾岩，为相对隔水层。测区详细的地层岩性见表1。测区内发育一条近EW向断层(图1)。

表1 地层岩性一览

图1 测区水文地质略图

1.2 测区地球物理特征

通过野外踏勘可知测区内出露地层依次为：上部为第四系覆盖物(Qh)——褐灰色粉砂质黏土、粉质砂土及砂卵石层，厚度在8 m左右；底部为石炭系黄龙群(C2hn)白云质灰岩。土层及岩石的电阻率变化见表2。地下介质除卵石层和灰岩层电阻率较高外，其他介质电阻率均很小。但随着地层岩性变化、岩石风化程度等不同，介质电性的差异也会发生改变。根据地层岩性及水文地质条件分析，测区内的含水地层为碳酸盐岩岩溶含水层，为相对低阻区，而且在音频大地电磁法(EH-4)的电阻率等值线断面图上会有低阻的异常圈闭反映。综上所述，本测区存在地球物理勘探前提。

2 野外工作方法

根据测区内构造、地形地貌等特点，布置了两条平行的物探测线(见图1)。测线方位110°，两条测线相距15 m，与区域断层有一定的交角。先对两条测线均采用高密度电阻率法，然后根据高密度电阻率法的异常结果并结合场地条件，在对应的异常段采用音频大地电磁测深法进行综合对比研究，最后确定井位。

2.1 高密度电阻率法

高密度电阻率法是以地下岩(矿)石的电阻率差异为基础，研究在供电电场作用下地下传导电流的分布规律，从而达到勘探地质体的目的。野外测量时将全部电极(几十至上百跟)一次性布设于测点上，然后通过程控多路电极转换器选择不同的电极组合方式和不同的极距间隔，从而完成野外数据的快速采集。

主要对1线、2线采用伦贝尔(α2)装置进行对称四极电测深测量[15]。其中1线测量点距7 m，采集道数57道，采集层数为20层，长度为392 m；2线测量点距5 m，采集道数60道，采集层数为29层，长度为295 m。仪器采用重庆奔腾仪器厂生产的WGMD-3型多功能数字直流激电仪。

2.2 音频大地电磁法

音频大地电磁测深法是以岩土层(体)的导电性与导磁性差异为基础，研究电磁波在地下不同介质传播过程中在地面表现的不同特征，达到了解地下介质的电性变化之目的。对测线上的异常段进行音频大地电磁测深法测量[16]，不极化电极沿测线方向测量电场Ex，用磁棒垂直测线方向测量磁场Hy，测量方式为EH-4标量，采集频段为高频段(750～100 000 Hz)。仪器采用劳雷工业公司生产的EH-4连续电导率成像仪。

3 综合结果分析与讨论

对采集到的对称四极电测深数据采用RES2DINV软件进行反演，反演结果见图2、图3。

1线的地表地形起伏不大。纵向上整体存在分层现象，表层介质电阻率值很小，普遍在100 Ω·m以内，是第四系覆盖层的反映；随着深度的增大，介质电阻率值明显增大，普遍大于500 Ω·m，是基岩白云质灰岩的反映。横向上覆盖层厚度变化较大，有些地方很厚，有些地方则很薄，变化范围为1～10 m，而在基岩中存在明显的横向电阻率值突变和低阻异常带，主要在274～323 m、358～421 m和449～484 m，推测在358～421 m测线段可能存在一条局部断层，倾向往测线小号方向。

由图3可知，测线地表地形起伏不大，纵向上存在分层现象，表层介质电阻率值有大有小，部分测线段有基岩裸露的情况，但整体上介质电阻率值还是偏小，主要还是第四系覆盖层的反映，随着深度的增大，介质电阻率值明显增大，是基岩白云质灰岩的反映。横向上来看覆盖层厚度变化较大，有些地方很厚，有些地方则很薄，变化范围为0～10 m，而在基岩中存在明显的横向电阻率值突变和低阻异常圈闭，主要在测线段355～415 m和450～470，其中355～415 m测线段可能存在一条局部断层，倾向往测线小号方向。

对比图2、图3发现，1线和2线这两条平行线异常段有很好的对应关系，推测的断层倾向也一致，为同一断层，平行测线在消除旁侧影响的同时，也提高了异常解释的可靠性。为了进一步提高解释的可靠性，对两条测线的异常段进行了EH-4音频大地电磁测量。具体为：1线测线段365～463 m，点距7 m；2线测线段350～450 m，点距5 m。使用IMAGEM软件对探测结果进行二维反演，在地形比较陡的地段取平滑因子为0.1，较平缓的地段取 0.5[17]，最后采用Surfer软件成图。图4和图5分别为1线和2线音频大地电磁测深电阻率反演断面。

图2 1线对称四极电测深反演结果断面

图3 2线对称四极电测深反演结果断面

由图4可知，在小号段存在明显的低阻异常圈闭，主要在测线段365～400 m，电阻率值偏小，而大号段电阻率值偏大，横向上形成电阻率值突变现象，推测突变处存在局部断层，倾角较陡，倾向往小号方向。由图5可知，小号段相对于大号段整体电阻率值偏小，横向上有明显电阻率值突变现象，推测突变处存在局部断层，倾角较陡，倾向往小号方向。而在小号段还存在低阻异常圈闭，主要在测线段360～400 m。

综上所述，基本可以确定1线365～421 m异常段和2线360～415 m异常段存在断层，且为同一断层，而根据低阻异常圈闭的位置可以初步确定钻井位置为1线375 m和2线375 m处。为了最终确定钻井的具体位置，分别抽取了1线375 m和2线375 m附近的3条电阻率测深曲线(单支曲线)进行分析讨论，见图6、图7。

由图6可知，3条单支曲线形态上大体相似，最小极距从10.5 m开始，电阻率值从大约75 Ω·m开始增大至约200 Ω·m，期间增速有所变化；初始阶段电阻率增大得较快，表明覆盖层薄，接着快速进入高阻层，后续增速有所下降，但下降得不明显，表明存在些许破碎；之后从极距大约60 m处电阻率值开始减小，随后又增大，表明存在明显破碎或存在充水溶洞等现象。而由图7可知，最小极距从7.5 m开始，电阻率值在不同测点有增大有减小，表明部分地方有基岩出露；之后的形态与图6很相似，电阻率值增大进入高阻层，期间有增速变化表明存在些许破碎；随后有一段增速平缓然后又加快，表明存在明显破碎或充水溶洞等现象。

通过以上综合对比分析，可以确定钻井位置为1线375 m和2线375 m处，且异常为同一断层引起，故最终选择了钻机容易开展工作的1线375 m处。图8为1线地质推断解释剖面，其中375 m处的钻井信息显示：0～3.4 m为耕植土、粉质黏土；3.4～7.8 m为卵石；7.8～38 m为白云质灰岩，其中在12～13.6 m、14～15.6 m、19～20.5 m、21～22.7 m、26.3～26.7 m、28.5～29.3 m及31～31.7 m为全充填充水溶洞段。由此可见，推断结果与钻井结果基本一致，该处抽水试验成果显示最大涌水量大于146.88 m3/d。

图4 1线音频大地电磁法电阻率反演断面

图5 2线音频大地电磁法电阻率反演断面

图6 1线375 m处三条单支曲线

图7 2线375 m处三条单支曲线

图8 1线地质推断解释剖面

4 结论

高密度电阻率法中的对称四极测深法能较好地对地下横向电阻率的变化进行响应，识别断层破碎带及岩溶破碎带的位置及发育方向，且纵向上可以利用电阻率测深曲线(单支曲线)对地下电性结构进行有效分层和评价。音频大地电磁法探测深度明显较大，对地下倾角较大的构造响应明显，能够很好识别断层破碎带及岩溶破碎带的位置及发育方向。两种方法互相验证，基本可以确定异常特征和具体的异常位置。