电剖面法探测花岗岩裂隙水

电阻率剖面法又称电剖面法。其工作原理与电测深法相同，该方法在工作中是采用不 变的电极距，并使整个装置沿着观测剖面移动，逐点观测视电阻率ρS的变化。由于电极 距固定不变，勘探深度就保持在同一个范围内，因而ρS值沿剖面的变化可以把地下某一 深度以上具有不同电阻率的地质体沿剖面方向的分布情况反映出来。
根据电极排列方式的不同，电剖面法又有许多变种。目前常用的有联合剖面法、对称 四极剖面法、中间梯度法以及偶极剖面法等。图4-12所示是它们的电极排列方式。由于 电极排列方式的差异，各种电剖面法所解决的地质问题也不同。但总地来讲，电剖面法适 于探测陡倾的层状或脉状金属矿体和高阻岩脉，划分接触带，配合地质填图，也能为寻找 含水断裂破碎带等水文、工程地质工作服务。

图4-12 电剖面法常用的电极排列简图

（一）联合剖面法
联合剖面法是两个三极排列AMN∞和∞MNB联合进行探测的一种电剖面法。所谓三极排列是指一个供电电极置于无穷远的排列，如图4-12（a）所示。A，M，N，B四 个电极位于同一直线（即测线）上，以MN之间的中点为测点，且AO＝BO，MO＝NO。电极C是两个三极排列共同的无穷远极，一般敷设在测线的中垂线上，与测线之间的距 离大于AO的5倍（CO＞5AO）。
工作中将A，M，N，B四个电极沿测线一起移动，并保持各电极间的距离不变。在 每个测点上分别测出A与C极供电时的电位差 和电流强度I，B与C极供电的电位差 和电流强度I，然后按式（4-12）求得两个视电阻率 和 ，即

勘探地球物理教程

式中：KA和KB分别为AMN∞和∞MNB两种排列的装置系数，根据计算可得

勘探地球物理教程

因此，联合剖面法的剖面曲线图上有两条视电阻率曲线 和 ，具体如图4-13所示。

图4-13 直立良导低阻薄板上的联合剖面曲线与对称四极剖面曲线的对比图 —联合剖面曲线； —对称四极曲线

图4-13 是直立良导低阻薄板（或称薄脉）上的联合剖面曲线与对称四极剖面曲线的对比图。从联合剖面曲线看， 和 两支曲线在低阻板顶相交，板左侧 ，板右侧 ，且交点处视电阻率值小于围岩的电阻率，因此在良导薄板上联剖视电阻率曲线出 现低阻正交点。直立低阻板上两条ρS曲线对称，在倾斜薄板上两条ρS曲线不对称，倾向 一侧的视电阻率极小值小于反倾向一侧的视电阻率极小值，因此， 极小值与 极小值连 线的倾斜方向与低阻板倾向一致。
当薄板电阻率ρ0大于围岩电阻率ρ1时，可称为高阻薄板。模型实验表明，直立高阻 薄板顶端联合剖面曲线有一个高阻反交点（即板左侧 ，板右侧 ）。倾斜高阻 薄板的联合剖面曲线呈不对称状，反倾向一侧的视电阻率极大值大于倾向一侧的视电阻率 极大值，两条曲线极大值连线的倾斜方向与高阻板倾向一致。
联合剖面测量结果可绘成ρS剖面图（图4-13），也可绘成ρS剖面平面图，通过剖面平面图中相邻剖面上ρS曲线的对比，可以了解地质体沿走向的变化情况，如图4-14 所示。
实际工作中情况比较复杂。地表电性不均匀可以引起ρMN的变化；地形起伏可以影响 电流密度jMN的分布；有时单纯地形影响就会引起与矿体相似的异常；相邻导电体的干扰 也可以造成异常的畸变。这些因素都将使ρS曲线大大复杂化。因此在进行资料解释之前，必须结合实际情况对实测曲线进行分析，消除一些干扰因素，辨认出由矿体引起的异常，才能进行正确的地质解释。

图4-14 某工区联合剖面法视电阻率剖面平面图

（二）对称四极法
对称四极法的电极排列如图4-12（b）所示。该方法在野外工作中始终保持A，M，N，B四个电极间距离不变，同时沿测线对称于M与N的中点O（即测点）向前移动，一般 ，这种方法测得的 值为

勘探地球物理教程

将条件AM＝BN和AN＝BM代入装置系数公式，可得

勘探地球物理教程

对称四极剖面ρS曲线见图4-13中的 。通过某些模型实验表明，对称四极 曲线 的异常幅度和分辨能力都不及联合剖面ρS曲线优越。但对称四极装置不需要笨重的无穷 远极，比较轻便，工作效率高，故多用于普查以及面积性的电阻率测绘等。
如果在对称四极装置的基础上，对称于O点再增加两个电极A′和B′，并且AB＞A′B′，这种装置称为复合对称四极剖面装置，见图4-12（c）。采用复合对称四极剖面法时，先 由大极距的A和B供电，M与N测量，求出 ；然后由小极距的A′和B′供电，M与N 测量，求出 。所以沿每条测线可得到两条不同勘探深度的视电阻率曲线。复合对称四极剖面法多用于探测覆盖层下基岩的起伏。图4-15（a）是基岩为高阻 的向斜构造；图4-15（b）是基岩为低阻的背斜构造。在上述两种情况下， 曲线（大 极距）具有相同的特征，即都有一极小值出现，所以单凭一条 曲线难以判别基岩起伏 情况。若用复合对称四极法就有可能解决这个问题。因为 曲线可以确定浅部的电性情 况，故在基岩为高阻的向斜上， 曲线低于 ；而在低阻的背斜上， 曲线位于 曲 线的上方。根据勘探目的和要求，复合对称四极法的极距选择，应力求大极距反映深部情 况，一般是AB/2≈（3～5）H（H是覆盖层的平均厚度）；小极距反映浅部情况，一般 A′B′/2≈（1～2）H。测量结果多用ρS剖面图表示。

图4-15 复合对称四极剖面法ρS曲线

对称剖面法的野外工作比较简单，生产效率高，表土不均匀和地形影响小。故常用于普 查阶段，以探测基底构造、基岩面起伏以及不同电阻率的岩层接触面等；在面积性工作中，常用于圈定某一深度范围内的视电阻率异常区，如岩溶分布范围与古河道的追索等。
图4-16是利用对称四极法追索古河道的例子。该图根据高阻异常特征清晰地反映了某地古河道的分布概貌。古河道是由一条主流和一条支流组成。古河道的视宽度由 剖 面曲线变化梯度最大的位置近似确定。如果曲线对称，说明古河床两岸也对称，极大值位 置对应古河床最深的中心部位。如果曲线不对称，可根据两翼曲线的陡缓来推断古河床两 岸坡度的陡缓。
图4-17是利用复合四极剖面法了解基岩起伏情况的例子。已知该区基岩深度大约为20～40m，因此选择小极距AB为40m，大极距AB为180m。由图可见，小极距的 有 一个高阻异常，推断这是浅部覆盖层中的高阻卵石透镜体的反映。大极距 曲线在小极 距曲线的上面，说明基岩的电阻率高于覆盖层电阻率。从大极距 曲线减去小极距曲线 即得出基岩起伏情况。

图4-16 对称四极法追索古河道的 剖面平面图


图4-17 某工区复合对称四极剖面法探测基岩起伏图

（三）中间梯度法
中间梯度法的装置特点见图4-12（d）。这种装置的供电电极距AB很大，通常选取AB为覆盖层厚度的70～80倍。测量电极距MN相对于AB小得多，一般选用MN＝ 。工作中保持A和B固定不动，M与N在A与B之间的中部约 的范围内同时移动，逐点进行测量。测点O为MN的中点。因为AB距离很大，在AB中 部测量范围内的电场可以认为是均匀电场。所以这种装置能最大限度地克服其他电剖面法 由于供电电极附近电性不均匀对视电阻率测量的影响。而且ρS曲线反映的必然是MN电 极附近地层电阻率的变化情况，这是该方法的优点。此外，这种装置在AB连线两侧 AB/6范围内的测线中部位置上，电场也近似于均匀电场。因此，在某一固定的A，B情 况下，又可以在相邻的AB/6范围内的测线中部进行中间梯度法测量。这种称为“一线布 置、多线测量”的方式，可以大大提高工作效率。
中间梯度法工作中，按下式计算ρS和K值：

勘探地球物理教程

但必须注意，装置系数K不是恒定的，测量电极每移动一次要计算一次K值。
中间梯度法用于寻找陡倾的高阻岩脉（如石英脉、伟晶岩脉等）可以取得显著的效果。因为在有浮土的情况下，高阻岩脉的屏蔽作用比较明显，排斥电流线使之汇聚于浮土，使 JMN急剧增加而ρS曲线上升形成突出的高峰；至于低阻岩脉，效果不明显。图4-18是我国 东北某铅锌矿区使用中间梯度法所得的ρS剖面平面图。该区铅锌矿位于倾角接近70°的高 阻石英脉中。图中两条ρS高值带是由含矿石英脉引起的。右边1号矿脉是已知的，左边2 号矿脉是根据中间梯度法ρS曲线形态与1号矿脉ρS曲线对比而圈定的。

图4-18 用中间梯度法追索石英脉

（四）偶极剖面法
偶极剖面法的装置特点见图4-12（e）。这种装置供电电极AB和测量电极MN均采用偶极并分开有一定距离，由于四个电极都在一条直线（即测线）上，故又称轴向偶极。其ρS表达式为

勘探地球物理教程

其中：

勘探地球物理教程

如果取AB＝MN，则

勘探地球物理教程

当令AB＝MN＝a，BN＝na（取n为正整数）时，则

勘探地球物理教程

偶极装置常取OO′中点为记录点。其中O为AB的中点，O′为MN的中点。OO′＝（n＋1）a。
偶极剖面法（简称偶极法）由于供电和测量均采用偶极装置，故又称双偶极或偶极— 偶极法，它一般采用轴向排列。该方法也是目前金属矿及其矿产资源调查中的一种常用方 法，尤其是与频率域激发极化法配合测量时，其应用更为普遍。
偶极法在各种金属矿上的异常反映是相当明显的，加之由于它的供电电极AB和测量电极MN是分开的，且所需导线均很短，因此它在减弱游散电流或电磁感应作用引起的 干扰方面，相对其他装置有明显的优越性。偶极法的主要缺点是，当极距较大时，在一个 矿体上往往可出现两个异常。故当有多个矿体存在或围岩电性不均匀时，曲线变得复杂，给解释带来困难。尽管如此，偶极法在实际找矿和填图工作中仍取得了较好的地质效果。
偶极装置的电极距OO′与联合剖面装置的电极距AO是相当的。因此，按照联合剖面 确定AO的原则，偶极的OO′应大于或等于3倍浮土厚度，而对于脉状矿体，则应取最佳 OO′＝L＋l（L和l分别为脉状体的走向长度和下延长度之半）。a的大小应为 不能取得太大，否则将会使异常减弱；但也不能取得过小，否则对围岩 电性不均匀将反映过于灵敏，使曲线呈锯齿状跳跃，并使观测电位差产生困难。
图4-19是偶极法进行地质填图的实测。图中为了查明某地覆盖层下花岗岩中破碎的 分布情况及确定古近-新近系中砂砾岩、泥岩与花岗岩的分界线，进行了偶极法的测量，为了便于对资料的解释，可将AB等于MN的所谓单边偶极ρ 曲线平移OO′后得到另一条单边 曲线。对于组成所谓的联合偶极剖面，这种“双边偶极”的作用与联剖相同，即可用两条曲线 的分布特征和交点性质来判断地下岩石的导电能力，从而较 准确地进行地质填图与找矿。

图4-19 用偶极法进行地质填图的实测ρs曲线

该区所用的电极距为AB＝MN＝10m，OO′＝30m。由图可见，在花岗岩破碎带上ρS 剖面曲线反映为低阻，并有明显的“正交点”。在古近-新近系砂砾岩、泥岩地层上，ρS更 低，形成一个明显的凹槽状。根据凹槽两侧 曲线上升最陡的地方，可较准确地 划定古近-新近系与花岗岩的分界线。而夹在古近-新近系与破碎带之间的花岗岩体上，ρS 表现为典型的宽高阻带，且有“反交点”存在。该区用“双边偶极剖面法”并根据ρS曲 线的特征进行的面积性地质填图取得了明显的效果。

在这一节，我们来看一些2-D高密度电法勘探的例子。
4.3.1 地下水探测
该测线数据是由丹麦Aarhus大学地球科学系完成的Grundfor区近地表松散沉积物和地下含水层中的岩性探测结果（Christensen et al.，1994）。反演模型中的低电阻率区（图4.28）为第四纪冰川粘土沉积层，而高电阻区则是砂质沉积层，并沿测线布置了一些钻孔，探测结论得到了证实。
温纳（Wenner）γ装置（或称为差分装置，图1.6c）有一个相对不寻常的布设方式，供电电极和电位电极相互交错，相对于温纳（Wenner）α、β装置来说，在实际勘探中，温纳（Wenner）γ装置很少使用。但是，在某些情况下，该装置可能会有一些优势，探测深度明显高于温纳（Wenner）α装置（0.59a与0.52a相比，见表3.2），然而电位电极比α装置少33％，相比而言，温纳（Wenner）β装置的观测电压是γ装置的1/3，在噪声严重的环境下，该装置是比较有利的。图4.29a为温纳（Wenner）γ装置在尼日利亚（Acworth，1981）Bauchi地区地下水探测的拟断面图。在该地区，地下水常见于结晶基岩的上覆风化层中，风化层较厚地区的基岩裂隙是含地下水良好的目标。在该地区，采用温纳（Wenner）α，β和γ装置进行了探测，同时使用Geonics EM34-3系统的电磁剖面探测（Acworth，1987），此处，作为一个例子，仅给出了温纳（Wenner）γ装置数据结果。为了突出土壤层和基岩之间的边界，使用了稳健反演技术。反演模型如图4.29b所示，低电阻率的风化层厚度一般为10～20 m，在桩号为190的下方有一不足20 m宽的垂直狭窄的低电阻率区域，可能是基岩裂隙带，在基岩裂隙带边缘，桩号为175的地方布设了一口钻井，出水量低于预期结果（Acworth，1987），在这种情况下，2-D电阻率模型对裂隙区中心位置精确定位将是有有益的，可以改善井的出水量。

图4.28 实测视电阻率视剖面及反演模型


图4.29 Bauchi地区地下水探测

4.3.2 岩脉探测
在图4.30的上部靠近剖面中间部位，岩脉出现了一个较为突出的高电阻率区域（Dahlin，1996）。在这套数据中，相对于较低的岩脉剖面来说，高阻岩脉的宽度远小于埋藏深度，反演非常困难，所以，浅层部位没有得到很好的解决，如果使用较小的阻尼因子，将会导致模型剖面中的电阻率值振荡。要反演这套数据，垂向平滑滤波值应该设置为横向滤波值的两倍，因为岩脉主要在垂向上，由于火成岩岩脉和周边的沉积岩之间的电阻率差异明显，使用了强制模型反演技术，在模型剖面中，岩脉表现为一个较突出的垂向高阻异常体。

图4.30 Odarslov岩脉实测视电阻率视剖面及反演模型


图4.31 滑坡勘探实例

4.3.3 滑坡探测
马来西亚所面临的一个自然灾害问题是山体滑坡，由于积水范围内的部分斜坡变得十分脆弱，从而诱发滑坡，图4.31为在一低阻山体滑坡地段上进行探测的结果，基岩（花岗岩）风化后，生成砂土与粘土相混合的不均匀物质。从反演处理结果来看，测线中心下方的低阻区较为突出，这可能是由于积水造成了电阻率降低，小于600Ω·m。为了加固边坡，就必须从该区域抽出多余的水，因此，对地下积水区域精确成像显得非常重要。
4.3.4 废旧工厂用地探测
在工业国家，一个常见的环境问题就是废弃的工业用地污染治理，这些用地还可以恢复，因此，有必要对埋在地下的老工业材料（如金属和混凝土块）进行成图。这些地区另一个问题是工厂地下的化学污染物，由于场地的自然现状，地下往往是非常复杂的，对大多数地球物理方法来说是一个挑战。图4.32是一个废弃工业用地垃圾渗漏液探测例子，探井较少，地面的废弃物正通过污物潟湖流入到地下砂岩内（Barker，1996），并从附近的小溪看到了渗透出的垃圾渗滤液，但是，该地下污染程度还不清楚。

图4.32 工业污染实例

沿着潟湖和小溪之间的老铁路路基进行高密度电法勘探，金属铁轨已拆除，但是部分嵌入地下的装载沥青大金属托架依然存在。在视电阻率剖面中（图4.32a），受污染的地下水聚集区显示为低阻，在测线桩号140右侧；金属物托架为一 “V” 字形低阻异常，测线桩号90附近。从反演模型（图4.32b）来看，反演处理已经成功重建了金属物在近地表正确的形状；在剖面右半部分有一低阻异常区域，由受污染的地下水所引起，依据剖面可以清晰地圈定出羽状边界在桩号140处，污染区域似乎已扩大至30 m深的地方。
4.3.5 粘土层洞穴探测
该项探测的目的是为了查清粘土层下8～20英尺（ft）处砂岩层中的洞穴，随后的钻井证实了高密度电阻率法探测的结果。此次探测的一条测线视剖面如图4.33a所示，该剖面中的数据采用横向重叠测线探测。这次探测，一个有趣的特点是展示了视剖面的误导性，尤其是对偶极-偶极（dipole-dipole）排列装置，在反演模型中，在200英尺以下存在一高阻异常，可能是一个低阻粘土层中的洞穴（图4.33b），在剖面中有塌陷特征的数据区域为一条裂隙。
但反演模型中的单元灵敏度值图表明，在视剖面中，该区域高阻体的模型单元（即更可靠的模型电阻率值）比在同一深度有更多数据点的邻近地区有更高的灵敏度值（图4.33c），这种现象是由于偶极-偶极（dipole-dipole）排列装置灵敏度等值线的不同所造成的，这些有最高灵敏度值的区域下方是C1-C2和P1-P2偶极，而并不是在排列中心下面的成图点。如果模型单元仅放置在数据点的位置，反演模型中将失去高阻体，这一地下重要的特征就可能不会被发现。

图4.33 粘土层中洞穴探测

4.3.6 抽水试验
在英格兰东中部的Hoveringham地区进行了抽水试验高密度电阻率法探测，该含水层为砂子和砾石层上覆泥岩，图4.34为初始时的电阻率拟断面、模型剖面和抽水220 min后的模型剖面，图4.35为抽水试验开始后40，120，220 min后的电阻率相对变化情况。

图4.34 Hoveringham抽水试验（英国）

图4.36清晰显示随抽水时间的推移，高电阻率值区域在增加，通过阿尔奇定律，并假设水的电阻率值不随时间变化，可以估计含水饱和度的变化情况，图4.36为含水层中的含水饱和度的降低情况，由于阿奇定律假设导电性由含水量单独引起，如果地下存在粘土，含水饱和度很可能比真实值低。
4.3.7 水下移动电极勘探
与传统的方法相反，实际可以在水下进行电阻率勘探，即使是海洋环境，这个例子是一套不寻常的数据，并且对任何电阻率成像反演软件来说均值得挑战，它不仅物理长度和电极位置数最长，而且还使用一个不同寻常的高度不对称的非传统电极布设，水下移动观测系统采集数据。移动观测系统具有观测速度较快的优势，但在陆地上电极与地面接触不理想（直接接触式）或信号强度低（对于静电型）。水下环境提供了一个直接接触式移动系统的理想环境，因为几乎不存在电极接触不良的问题。
图4.37a为一条沿河进行勘探的测线，长18km，该项勘探由比利时Sage Engineering承担，其目的是摸清河床近地表岩性，进而计划铺设电缆。这套数据共7479个电极位置和6636个数据。在反演模型（图4.37b）中，河床大部分物质电阻率小于120Ω·m，其中近地表几个区域的物质的电阻率明显较高，超过150Ω·m，不幸的是，在这个区域的地质信息相当有限。在高电阻率地区，潜水员面临的问题是取得沉积物标本；低电阻率的物质可能是较多的相干沉积物（可能有淤泥/粘土砂），而高电阻率区域可能是粗糙，并且是不太相干的物质。

图4.35 Hoveringham地下抽水试验电阻值相对变化百分比a—抽水40 min后反演模型剖面；b—抽水120 min后反演模型剖面；c—抽水220 min后反演模型剖面为了突出地下电阻率的变化，在模型中电阻率的变化显示出来，请注意，井孔下的模型电阻率随时间增加

4.3.8 水面浮动电极勘探
这项探测是沿着经过康涅狄格州的泰晤士河（美国），通过船尾拖动水面漂浮电缆进行的。该电缆固定2根供电电极和9根电位电极，观测采用一套8通道电阻率仪系统，在探测过程中，水深和水的导电率也进行了观测，采用偶极-偶极（dipole-dipole）装置，但是，在电位偶极与供电偶极长度不相同的时候，采用非对称布设采集数据。图4.38a为一条测线的视电阻率拟断面，没有采用特殊约束获得的反演模型（图4.38b），由于数据中存在噪声，模型内部分地方出现水层落差的假象，当河底物质成高阻时，该模型是相当准确的；但是，在低阻地区，精确程度就受到一定的影响，这些问题不会出现在水层电阻率固定的反演模型中（图4.38c）。
图4.39为一次沿巴西圣弗朗西斯科河勘探的模型单元和水底边界的布设情况。此次勘探装置为偶极-偶极（dipole-dipole）装置，偶极距为5 m，但是，数据采集时，采用每个正常单位电极距（1 m）读一个数，为了减少模型单元的数量，同时考虑到这种勘探的分辨率不可能超过半偶极长度，因此，该模型的单元格宽度设置为3倍单位电极距，即整个剖面上的单元格宽度基本上为3 m。视电阻率剖面和反演模型如图4.40所示，反演时，水层的电阻率变化减小到最低限度，因此，水层电阻率一般不均匀，除了测线左侧，因为那里数据点不多，在190～220m之间显示，河床近地表沉积物成低阻。

图4.36 利用阿尔奇定律进行Hoveringham抽水试验法

4.3.9 考古探测
晋阳古城遗址是太原宋代以前旧城遗址，位于太原市晋源区，自春秋末赵简子家臣董安于始筑，已有2500余年的历史。由于古城遗址内容丰富，遗址特征复杂，如古城墙具有占地范围大、埋藏较深、夯土构建、土质杂乱且坚硬，护城河具有占地范围大、埋藏较深、所含物质成分主要是泥质粉砂、含水丰富、物质比较均匀、成层性较好，而对于浅层古文化遗迹，主要是古建筑基址、古代人类活动文化层等等，这些文化遗迹层与原始未扰动地层都有非常明显的物性差异，尤其是电性（电阻率、磁导率等）与非文化层存在相当大的差异；另外，晋阳古城坐落在山前（龙山与蒙山）冲洪积扇上，古城底基由于分选性差，而且物质不均匀，砂砾层空隙度较大，含水相对较少等特点，这就为应用高密度电法探测技术探测在晋阳古城考古提供了非常好的地球物理前提和基础。

图4.37 水下河床探测4次迭代反演模型（Sage Engineering，Belgium）


图4.38 泰晤士河（CT，USA）浮动电极调查

采用温纳（Wenner）α装置，在勘探区布设了一条点距为6 m的高密度电法测线进行探测，除担负摸清护城河与古城墙的电性规律外，还肩负探清该区域砂砾层厚度和潜水面深度。由于古城墙就地取土，夯土构建，土质杂乱且坚硬，在电性上主要表现为高阻异常，虽然现在护城河年代已久，被废弃了，甚至被土掩埋，但所含物质成分主要是泥质粉砂，含水丰富，因此，常表现为低阻异常。由于潜水面以下含水相对丰富，电阻率表现为低阻，从图4.41可清晰看出，该区域西边在36 m以下，东边大概在24 m以下，即有西深东浅的分布特点，成一斜坡状分布趋势，视电阻率均比较低，可断定为潜水面；由于晋阳古城坐落在山前冲洪积扇上，西边为扇根，东边为扇前，即砂砾层西厚东薄，从该图上也表现得非常明了，由于砂砾层空隙度较大、干燥、含水相对较少，而且物质不均匀，常表现为高阻异常，因此，从该图上可定性断定潜水面以上均为砂砾层，当然，古城墙跟前的高阻区域也不排除是古建筑遗迹区域，有待进一步做工作。

图4.39 模型单元和水底边界的布设


图4.40 巴西的圣弗朗西斯科河床探测


图4.41 古城西城护城河、古城墙及下伏砂砾厚度探测视电阻率剖面（点距6m）山前冲洪积扇成东薄西厚的分布趋势

4.3.10 复杂岩溶地区路基勘测
柳肇铁路CK51＋200～CK51 ＋900段岩溶较发育，地质结构较复杂，其中地表以第四系粘土、卵石土为主，下伏基岩主要为石炭系泥灰岩、灰岩和硅质岩。采用温纳（Wenner）α装置进行勘探，电极间距5m，探测深度80 m。

图4.42 柳肇铁路复杂岩溶地区路基勘测结果

高密度电法勘探反演结果如图4.42a所示，通过探测，发现表层覆盖层电阻率为200～1000Ω·m，推测深度为10 m左右，主要为粘土。同时发现表层下面有2个高阻体，其电阻率为800～4000Ω·m，推测为完整灰岩或硅质岩，由于高阻体较厚，且下面还有相对低阻体，为了保险起见，在CK51＋400附近布置一个80 m深的钻孔，钻探结果为:前60 m为卵石土含硅质岩，60～80 m之间为灰岩。结合钻孔资料，得出左边的高阻体不是完整灰岩的反映，而是卵石土含硅质岩的反映，依此推断，右边的高阻体也是卵石土含硅质岩的反映，其最大厚度可达60 m，2个高阻体在浅部地区相通。左边高阻体下面有一相对低阻体，电阻率为30～800Ω·m，钻孔资料已经证实为灰岩，从灰岩在该区的电阻率值来判断，岩溶较发育。溶蚀灰岩右边有一低阻体，电阻率为50～250Ω·m，推断为泥盆系灰岩或泥灰岩，在CK51 ＋650附近布置一个80 m钻孔，结果证明该解释正确，钻孔打出来的是泥灰岩，解释成果详见图4.42b。通过此次高密度电法探测，准确而翔实地划分了路基通过区域的地层结构和岩溶发育情况，为以后的岩溶灾害整治工作提供了可靠的依据。
除了这些实例以外，2-D高密度电法已经成功应用于许多其他探测领域，如垃圾填埋场、基岩起伏的石灰岩地区污染物泄漏检测、基岩覆盖层厚度成像（Ritz et al.，1999）、大坝漏水、海水入侵、（Dahlin et al.，1998）、地下水示踪剂监测（Nyquist et al.，1999）、松散沉积物探测（Christensen et al.，1994），高密度电阻率法也可以用来于湖泊和水库水下探测。
4.3.11 城市活断层探测
长春地区第四系覆盖层厚18 m左右，基岩为白垩系地层，以砂砾岩和泥岩为主，胶结微弱孔隙较发育。北北东断层有火山岩体侵入是富水断裂带，产状近直立垂直断距不大水平断距相对较大。根据地质情况，采用联合剖面法、高密度电法、瞬变电磁测深探地雷达、高分辨率地震反射和重磁方法等进行对比试验研究，在活断层上取得了良好的效果。
图4.43是在长春市区南用高密度电法探测的NE40°断裂带，中部向西倾的低阻异常反映了断层的形态，它与电阻率联合剖面的低阻正交点和地质推断均吻合较好。

图4.43 高密度电阻率二维反演拟断面

4.3.12 水库大坝塌陷勘测
随着水库服务年限的增加，水库大坝坝体局部渗漏、管涌现象出现较多，按规范采用常规钻探、坑探方法，无法对坝体渗漏、管涌现象全面了解。而高密度电法是进行水库大坝坝体勘测的首选方法，它具有采集数据密度大，数据反演成果反映剖面和深度双重性质，横向、纵向分辨率高的特点。
吉木萨尔县贡拜沟水库坝型为均质土坝，基础为基岩，根据实测，坝体土视电阻率值一般为2～40Ω·m，坝体土视电阻率随土体含水量而变化，当坝体含水量高时，视电阻率降低，当坝体含水量底时，视电阻率升高。坝体空洞表现为高电阻率，一般大于100Ω·m，基岩视电阻率值一般大于500Ω·m。场地具有进行高密度电法工作的物理前提，根据大坝坝体塌陷情况坝体布置4条高密度测线，垂直坝体布置2条。
从探测反演结果来看（图4.44至图4.49），在坝体上发现局部异常点3处，异常为高阻，电阻率在110～300Ω·m之间变化，异常形态位置详见图图4.50，其余方向均没有发现延伸异常，探测结果基本控制贡拜沟水库塌陷影响范围。后期对3号异常、2号异常打钻验证，验证部位，3线桩号94.0，钻井深25.0 m；2线桩号29.0，钻井深15.0m；结果在2号、3号异常位置发现塌陷区，与高密度电法勘探异常吻合。

图4.44 测线1高密度电阻率反演结果


图4.45 测线2高密度电阻率反演结果


图4.46 测线3高密度电阻率反演结果


图4.47 测线4高密度电阻率反演结果

4.3.13 管线探测
城市地下管线主要包括煤气、自来水、污水、雨水、通信、暖气管线等等，管线在城市发展过程中发挥了极其重要的作用。但是，地下管线在地面以下层层交错，错综复杂，形成了网状的地下管网，随着城市的改扩建和管线维护过程中，管线的探测显得非常重要。高密度电法由于采样间距和点距小、速度快，能有效满足地下管线勘测高分辨率的需求，是一种有效的管线探测方法，图4.51为管线高密度电法探测模型正反演结果，其中管线材质为钢材，效果明显。

图4.48 测线5高密度电阻率反演结果


图4.49 测线6高密度电阻率反演结果


图4.50 坝体老放水涵洞高密度勘探断面

待测管道实例为上海南汇区某新建公路路面下埋设的排污混凝土管道，根据施工方原始资料此管道走向东西，直径大约有2m，且在某个位置被长约2.4m铁皮包裹，施工单位要求准确找出包铁皮管道在地面上的精确位置，以便进行新修路面下水管道的埋设与对接。该测区地质情况较简单，地表被第四纪地层所覆盖，岩性为粘土层，所勘探的地下管道基本上位于该层地层内。一般说来，粘土电阻率为0.5～30Ω·m，水泥混凝土管道的电阻率大大的高于此值，铁皮的电阻率则低于此值，这为分辨地下管道电性差异提供了良好的基础。

图4.51 高密度电法管线探测模型


图4.52 高密度电法管线探测反演结果


图4.53 测线布置及管线解释图

图4.52为高密度电法管线探测成像结果，根据4条测线的反演结果，得出包铁皮管道的中心位置应在原先假设中心点的东南方，具体位置偏离原中心点东约0.5 m，南约0.25 m，如图4.53所示。管道中心深度约为3m，由于管道半径为2m，故施工中下挖2 m即可发现管道。后经开挖验证，与推测的平面分布位置和深度均相吻合。
4.3.14 采空区及塌陷区勘查
采空塌陷是常见的地质灾害之一，主要发生在以井巷开采的矿区，特别是开采矿体埋藏较浅、产状较平缓的矿区。若开采矿体埋藏相对较深，且不能及时回填，当采空面积达到一定规模后，也会产生大面积塌陷。
实例1:采空塌陷区为贵州省毕威高速公路某段地质灾害，岩层单斜，产状180°∠24°，节理裂隙较发育，未见区域性断裂面通过，地质构造相对简单。依据以上分析，观察电法剖面（图4.54a）可知，塌陷区底界面的水平投影区间由120 m处向大桩号方向延伸至230 m处，埋深由57 m缓慢升高至地表，视电阻率范围为30～150Ω·m；塌陷中心位于120 m处，埋深约57m，视电阻率值为30～50 Ω·m。钻孔资料（图4.54b）揭示:ZK5于30 m钻遇灰岩夹层，ZK6于37 m钻遇灰岩夹层。结合岩层产状推断，塌陷中心位于ZK5与ZK6之间（即90～130 m之间），最大塌陷深度不小于7m，与物探解释的塌陷中心位置（120 m）相吻合。

图4.54 高密度电法勘探成像剖面及地质解释

实例2:河北省沙河市某铁矿采空区位于河北省邢台市境内，属太行山脉与华北平原过渡丘陵区，地表为第四纪黄土夹坡积碎石层覆盖，基岩有奥陶纪灰岩，以低缓角度向东倾覆，与下部闪长岩体接触，其矿体属位于接触带上的矽卡岩型磁铁矿。多年来的乱采乱掘使矿区出现了近SN走向的大面积塌陷和裂缝，导致耕地被破坏，甚至威胁到人们的生命和财产安全。为对该矿山采空区沉降程度作出正确评估，提出合理的治理方案，需要查明地下150 m内采空区的分布范围，采空区 “三带” 分布状况。对此，在开展地形测量和地面地质调查的基础上，采用了高密度电阻率法，布置4条近EW走向的剖面，线距:200～300 m不等，两条横穿测区中心的近SN向剖面，以基本控制采空区的分布范围。对于剖面I（图4.55），依据倾斜分布的异常梯度带作为采空区及冒落裂隙带边缘的标志，推断采空区边界分别位于10和285号点，两处梯度带均为向内相向倾斜（即倒 “八” 形）分布。依据反演剖面内的电阻率异常分布特征，结合地面裂隙分布和地下水位，即属未含水采空区呈高阻性特征，并将规模单一高阻异常划定为沉降非稳定区的反映，局部高、低电阻异常组合分布区为沉降稳定区的反映；其上部水平带状分布的低阻带，是因裂隙带内存有一定的滞留水所致，从而成为圈定裂隙带的依据；地面已出现明显的断陷裂隙带和多处塌陷坑，表明该采空区以沉降塌陷至地表，不存在弯曲带。上述解释成果可以对该剖面的沉降发展现状有所了解，采空区沉降裂隙带位于地表层，其厚度20～25 m，下部采空区的冒落带中间75～225点属不稳定区，仍有冒落可能，其两侧属冒落稳定区，不再有明显的塌落沉降。相邻剖面也存在与其可对比的解释结果，将各剖面解释结果绘制成采空区平面解释图，将可作为本区矿山开采沉降治理的依据。

图4.55 测线I视电阻率及反演断面图

实例3:辽宁省某大型露天铁矿，储量巨大，开采历史悠久，但是由于历史原因，尤其是多年前周围众多集体和私营小矿点在开采过程中的不规范行为，使矿区内遗留下多处采空区。随着近年来矿石需求量的猛增，矿石开采量增长迅速，原来深埋于地下的采空区逐渐暴露于地表，给矿山生产带来严重的安全隐患，治理采空区迫在眉睫。然而，现存的地下采空区资料稀少，难以满足采空区治理的要求。因此，对地下采空区进行探测，确定地下采空区的位置、埋深、规模和区域分布，不仅为矿山进行地下采空区的治理提供详细可靠的资料，还是保证矿山人员设备安全和生产计划顺利执行的必要前提。工程探测区域位于该矿山一开采平台之上，地形平坦，基岩裸露，南北长约300m，东西宽约200m，四周为矿石开采而形成的陡崖。区内主要地层有太古宙鞍山群、元古宙辽河群，其中太古宙鞍山群是赋存条带状铁矿的地层。区域内的岩石分布，主要为黑云绿泥石片岩、磁铁矿、绿泥石片岩，并可见少量石英脉。矿体的走向近东西向，倾向北东，倾角大约60°～85°。矿体与围岩间存在明显的地质界线，富矿体与贫矿体为渐变关系。矿石中矿物种类较少，主要为石英、氧化铁矿物、闪石矿物。矿石品位一般为20％～45％，平均品位一般为25％～35％。本地区条带状铁矿石矿物学特征可以分为磁铁石英岩及假象赤铁石英岩两大类型。据如上测区地质和地球物理特征和高密度电法探测的理论依据，采空区的地球物理特征与围岩的差异显著，符合利用高密度电阻率法进行综合探测的地球物理条件，理论上可应用该方法进行探测，并根据阻值异常分布确定空区分布。根据探测结果（图4.56），视电阻率大于2000Ω ·m的区域划定为采空区，高密度电阻率法对地下采空区产生明显的视电阻率值异常反应，如此划定的采空区范围为X∈（128，160）、Y∈（26，44），与实际范围的最大边界误差不超过10％。划定采空区的范围略大于实际范围，尤其是垂直方向误差较大。根据资料及电阻率法勘探原理分析，可能是由于采空区的存在，破坏了原有的力学平衡条件，使围岩受到巨大的荷载压力，超过了岩石的承载能力，从而产生了大量的裂隙，尤其是在荷载较大的垂向，裂隙更发育，岩石原有电性特征发生改变，电阻率增大。

图4.56 采空区电阻率拟断面图

2004年5月上旬为河北省涞源县支家庄铁选厂找水定井，使用了电剖面法。电剖面法在已知井位（20线）和待定井位（12线）的视电阻率ρS曲线如图1－6－3和图1－6－4所示。

图1－6－3 涞源县支家庄20测线ρS曲线

图1－6－4 涞源县支家庄12测线ρS曲线

井位处于南北向谷地中，东西两侧出露燕山期花岗岩。沟东西向宽约100m，第四系坡洪积相碎石层厚12m左右。所测ρS曲线低阻形态清晰完整，采用线性方程式法和机器模式识别法计算了95/12号点井位涌水量后，推测其涌水量约730m³/d，同年7月钻探成井后涌水量696m³/d。现将定量计算程序叙述如下。

1.线性方程式法

由图1－6－3可知：

地球物理找水方法技术与仪器

已知87/20号点水井涌水量Q_x为720m³/d。

由图1－6－4可知：

地球物理找水方法技术与仪器

换算系数：

地球物理找水方法技术与仪器

地球物理找水方法技术与仪器

将上述换算系数代入（1－5－3）式，（1－5－4）式，计算95/12号点待定井位涌水量Q_y。

地球物理找水方法技术与仪器

上述二式计算结果算术平均值：

地球物理找水方法技术与仪器

2.机器模式识别法

将图1－6－3和图1－6－4中ρS曲线低阻异常带放在虚线所圈定的3×4小方格内，已知：F＝12，列出矩阵，并求其和。

地球物理找水方法技术与仪器

得不符值D＝2，小方格数F＝12，计算贴近度：

地球物理找水方法技术与仪器

地球物理找水方法技术与仪器

计算结果和实际值（696m³/d）相比较：绝对误差是77m³/d，相对误差是11.1％。综上所述用电阻率法计算井位涌水量方法实际应用效果良好。