1 、引言

水库淤积是一个普遍存在的问题，泥沙的大量淤积，一方面降低了水库的防洪能力，一方面水库淤积会引起回水上延，吞蚀既有农田耕地，影响农业生产及社会和谐稳定，因此解决水库淤积测量存在的问题是十分必要的。水库淤积的地球物理探测手段一般有水上地震法、声呐探测技术、探地雷达技术和高密度电法。本文应用高密度电法对某库区泥沙淤积进行调查，连续、快速、有效的探测出了水深和淤积层厚度，探测资料与水上验证钻孔资料有很好的对应性，为后续水库清淤设计方案提供支撑资料。

2 、水上高密度电法原理

高密度电法原理上属于电阻率法的范畴，是电测深与电剖面的组合，其观测点密度大，获得信息量丰富，可以较详细探测水平和垂直方向上的电性变化的一种电法勘探。高密度电法是一种阵列布置的物探方法，也称自动电阻率系统，是直流电法的发展，其功能相当于四极测深与电剖面法的结合。通过电极向地下供电形成人工电场，该电场的分布与地下岩土介质的视电阻率ρs的分布密切相关。测量时，AM=MN=NB=AB/3为一个电极间距，A、M、N、B同步向右移动，得到第一层深度的剖面线；接着AM、MN、NB增大一个电极间距，A、M、N、B同步向右移动，得到第二层深度的剖面线；依此类推，这样通过对地表不同部位人工电场的扫描测量，得到视电阻率断面图像，由此来了解地下介质视电阻率ρs的分布，根据地下介质视电阻率的分布推断解释地下地质体的情况。在此工程中即可用于探测淤积层厚度、分布情况。原理示意图见1-1，水上工作照片见1-2。

图1-1 水上高密度电法探测原理示意图

图1-2 水上高密度电法施工照片

3、工程应用实例

云南某水库自1957年兴建，1958年竣工，至今已运行59年，库功能以灌溉为主，兼顾下游农田、乡镇、国道、工矿及企业的防洪。现今，水库的淤积严重，影响着水库的灌溉效益，一定程度制约受益区经济社会的发展，所以水库清淤工作势在必行。现利用高密度电法对水库水下淤积情况进行探测，需初步查明水库库区淤积层厚度、分布情况，为后续水库清淤设计方案提供支撑资料。

3.1工作布置

某水库淤积层探测共在库区布置29条物探剖面，共布置14个验证钻孔，实际物探测线布置图见图1-3。

图1-3 某水库水上高密度测线、验证钻孔布置示意图

3.2解释原则

①下图1-4为某水库WT16剖面视电阻率断面成果图，测线长度为300m，桩号为0m～300m。由视电阻率断面图可以看出：剖面电阻率在2Ω•m～140Ω•m之间，电阻率分布均匀，基本成层状分布，视电阻率水平方向上随基岩起伏变化。根据电阻率等级分布，以等值线最密集（电性差异变化最大）的地方为界面，定性的将该剖面大致分为两个电性层：第一层底界面深度为10m左右，电阻率在22Ω•m左右；第二个电性层底界面深度为14m左右，电阻率在28Ω•m左右。电性层分好后，结合地质和钻探资料，定量进行综合解释推断，见图1-5。该剖面桩号165m处通过钻孔WT1601，由地质资料、水深测试资料、钻孔资料知：WT1601钻孔终孔深度为4.5m，该桩号处水深为8.5m，淤积层厚度为4m。综合定量解释后见图4.2-6，将第一个电性层的底界面解释为水底界面，第二个电性层的底界面解释为淤积底界面。

图1-4 某水库WT16测线视电阻率断面成果图

图1-5 某水库ZK1601钻孔岩芯与现场工作照片

图1-6 某水库WT16测线视电阻率断面成果图

②下图1-7为某水库WT2测线视电阻率断面成果图，测线长度为360m，桩号为0m～360m。解释原则基本同上介绍，但是在剖面桩号80m～160m的地方，视电阻率整体较低，推测为淤积层以下长石石英砂岩含泥岩、泥质粉砂岩夹层影响或基岩起伏导致的体积效应，此时不应该盲目的将淤底往深处勾画解释，而应该结合地质情况，钻孔情况，水库淤积的基本规律（如：在库水位变化幅度不大，淤积处于自由发展情况下，水库淤积一般呈三角洲形态；多沙河流上的水库普遍有淤积一大片，冲刷一条带的特点，淤积一大片指泥沙在横断面上基本呈均匀分布，等等）结合视电阻率等级分布断面图来进行合理的推断与解释。

图1-7  某水库WT2测线视视电阻率断面成果图风化。

3.3淤积量计算与探测结果分析

3.3.1淤积分布

某水库淤泥厚度分布成果图见1-8，淤积层厚度的区域分布图是根据每单条测线测试的点数据整合插值形成的二维空间面数据。由成果图可以看出：在纵断面上，泥沙淤积主要发展到坝前，成锥体状淤积形态，坝前淤积多，沿库尾方向淤积逐渐减少，甚至未发现明显淤积；在横断面上，水库淤积基本是平淤，泥沙在横断面上基本呈均匀分布，库区横断面上不存在明显的滩槽。淤积层厚度最大的地方在测线WT11桩号115m附近，该处测试水深为10.1m，水下淤积层厚度约为8.6m；29条测线平均水深为7.4m，最大水深12.3m，淤积层平均厚度为4.5m，最大淤积厚度为8.6m。

图1-8  某水库淤积厚度分布示意图

3.3.2淤积量计算

计算原理：物探共布置29条测线，测试后得到每条测线的二维视电阻率断面成果图，根据断面图推测出相应点的淤积层厚度，最后插值成二维空间的厚度数据体，得到淤积层厚度等值线图。淤积层体积算法为：将该整个库区面积近似划分1m×1m的单元网格，每个单元网格根据插值的二维数据体给出一个淤积层平均厚度，由此计算出每个网格淤积层的体积，累加每一个网格的体积即得到水库淤积层的总体积。经计算的狮子口水库淤积体积约为108.3万m3，按照水库原设计总库容707.8万m3计算，淤积体积约占库容的15.3%。

4 结论

(1)该水库钻孔解释厚度有数值对比的共有7对数据，平均误差0.21，最大相对误差10%。说明使用水上高密度电法进行水库淤积层厚度探测是可行的。实际探测中从探测剖面图上能正确识别水底及淤泥层底界的界面，进而计算出水深和淤积层厚度。由于水库淤积为稀软泥沙，周边基础松软，建议使用挖掘机清理淤积。

(2)总结整个探测过程形成的工法可保证探测技术有效、规范实施。

(3)探测实例表明, 该方法具有连续、快速、有效的特点。

参考文献：

[1]雷宛，肖宏跃，邓一谦.工程与环境物探教程[M].北京:地质出版社,2006.

[2] 王俊，雷宛，李星,等.综合物探方法在工程基础勘探中的应用[J].物探化探计算技术,2013(3):35

[3]王兴泰．工程与环境物探新方法新技术［Ｍ］．北京：地质出版社，1996．