数值差分法在大坝渗漏工程中的应用.doc

数值差分法在大坝渗漏工程中的应用 赵华新 ，李跃群 （ 河南省水利勘测总队， 河南郑州 450003） 摘要：结合数值差分法在高密度电法解决大坝渗漏工程中的应用，实践证明，数值差分法对资料处理，压制区域场，提高信噪比有一定的实用价值。 关键词：大坝渗漏；高密度电法；数值差分法； 用高密度电法进行大坝渗漏勘探时，所获得的地质信息，往往因为区域场强，信噪比低而难于提取有用信息进行解释推断。为提高高密度电法资料的解释推断水平，从数据处理、信息提取入手，将数值差分法应用于高密度电法曲线解释工作中，突出了有用局部异常，从而使解释工作上了一个新台阶，并在井子沟水库大坝渗漏探测实践中取得了明显的地质效果。 数值差分法在高密度电法中应用的基本思想 高密度电法为人工场源，当地下为均匀介质时，所得到的信息为稳定的区域场。当较弱的有用局部异常叠加在较强的区域场时，等值线仅仅表现为局部发生微弱的弯曲，并未形成封闭异常，难以达到满意的地质效果。如果对数据做二阶微分处理，则区域场受到压制，而局部异常得到突出，从而把有用的局部异常分离出来，这样有利于对异常体的解释推断。由于高密度电法测得的数据是离散性数据，相邻极距差均为1米，所以二阶微分运算可以用二阶差分运算来近似。 通过差商的极限来定义导数 ( 1 ) 用泰来基数展开可得 ( 2 ) 式中，△f（χi）为向前差分。当步长△χ充分小时，（2）式中右边第二项趋于零，则 当步长呈等距离变化时，，其中。 由导数的定义可知，当是离散数据时，差商代替微商。 前差商代替微商： ； 后差商代替微商： ； 中心差商代替微商： ； 同理，二阶导数则用以下定义的二阶差商近似 其中。 由以上可知，当步长为等距离变化时，二阶微分运算可以用二阶差分来近似。 在高密度电法勘测中，AB/2电极系列为2，3，4，5，6，7，8，9，10，…20，均按等距离变化，即步长为常数。设为，AB/2为x，并且，则由二阶差分得式中，，分别为第，，电极距对应的电阻率。 经计算求得的电阻率值可以作为解释推断所用数据。测深断面绘制在算术坐标纸上，地下若存在均匀介质时，曲线为单层曲线；当地下有不均匀体存在时，电测深曲线发生畸变，在对应不均匀体两侧出现拐点，且拐点处的二阶微分为零，微分值过拐点后正负发生变化。在用数值计算方法处理后的电阻率值绘制的等值线图中，出现负异常封闭圈，即对应大坝渗漏体的范围。数值差分法在大坝渗漏中的应用，使区域场受到压制而突出局部异常，且叠加异常得到分离，有利于解释推断。 数值差分法在大坝渗漏中的应用实例 井子沟水库位于豫西某地，该水库一到雨季，水库水位上升，大坝出现渗漏现象，给大坝安全造成事故隐患。大坝近东西走向，北侧为迎水面，大坝渗漏点在大坝两端南侧。经查阅有关资料和实地踏勘，区域内地层、岩性等地质情况如下：大坝一带出露地层为三叠系上统椿树腰组地层，出露主要岩性为灰绿色粉砂岩和灰绿色粉砂质页岩。这种软硬相间的地层结构为大坝渗漏灾害的发生、发展提供了先决地件。 由于大坝渗漏影响范围相对较小，常规电法难以有效解决地质问题，所以在渗漏点附近， 布置了2条剖面，选用高密度电阻率测深法进行勘测，点距2米 ， AB/2=2~20米。 图1a是井子沟大坝1号测线实测视电阻率数据绘制的断面等值线，可以看出18~28号点之间AB/2=4~12米处存在25~35欧姆.米高阻半封闭圈，下部电阻率均匀；58~68号点之间AB/2=6~14米处存在25~30欧姆.米高阻半封闭圈，难以进行范围划分。图1b是同一条测线经过数值差分法处理后绘制的等值线，其电阻率变化的内容更为丰富。我们感兴趣的是22~26、58~64号点之间，出现两个负异常带，对应的电极距深度为5~12m。这说明有地电不均匀低阻体存在。 图2a是井子沟大坝2号测线实测视电阻率数据绘制的断面等值线，同样在18~28号点之间，58~68号点之间出现高电阻率半封闭圈。图2b是同一条测线经过数值差分法处理后绘制的等值线，其电阻率变化明显。图2上20~26、58~64点号之间，出现与1剖面相对应的两个负异常带。结合地质特征和自然电场资料，认为西边22~26号点之间，异常较窄，埋深浅，是一些洞穴引起的大坝渗漏；东部异常在58~66号点之间，宽度8m，厚度2m，顶部埋深5~8m；是大坝回填土与坝肩基岩松软接触带的反映，也是构成大坝东部渗漏的主要原因。 以上由数值差分方法进行的解释推断，经西部工程开挖和东部钻孔验证，结果较吻合，说明数值差分方法是有效的。 结论 数值差分方法，在解决大坝渗漏中取得了一定的地质效果。尤其是城市物探或其它工程物探生产中，勘探的目的主要是浅