堤坝集中渗漏温度场探测模型及数值试验

1、第25卷 增2岩石力学与工程学报 V ol.25 Supp.22006年10月 Chinese Journal of Rock Mechanics and Engineering Oct. ，2006收稿日期：20051214；修回日期：20060421基金项目：国家自然科学基金资助项目(50179009；国家自然科学基金重点项目(50139030堤坝集中渗漏温度场探测模型及数值试验王新建，陈建生(河海大学 岩土工程研究所，江苏 南京 210098摘要：为确定堤坝集中渗漏通道的位置和性质，在反分析理论的基础上，建立齐次边界条件下的各向同性土体热传导模型和复杂边界条件下的渗流传热模型；求解简化边

2、界条件后的解析解，对方程解进行变量转换，把渗漏通道位置和其他参数作为未知变量，建立反分析目标方程，指出实现回归优化的多条途径。为了验证其有效性，利用有限元软件设计数值试验，参照温度场特征，对两种主要模型不同探测线、最低温水平面及综合数据进行回归分析对比。工程实例优化结果表明，利用温度场进行集中渗漏探测是有效的方法。 关键词：水利工程；温度场；数值试验；渗流；回归分析；优化中图分类号：TV 31 文献标识码：A 文章编号：10006915(2006增2379408MODEL FOR DETECTING OF CONCENTRATED LEAKAGE IN DAM ANDDYKE AND ITS

3、NUMERICAL EXPERIMENTWANG Xinjian，CHEN Jiansheng(Institute of Geotechanical Engineering，Hohai University，Naijing ，Jiangsu 210098，China Abstract ：Concentration leakage is important cause of dam and dyke safety，so leakage intensity，scope and position of leakage passage need to be determined in order to

4、 intercept leakage in dyke and dam. In this paper，heat conduction model with homogeneous boundary condition in isotopic soil and heat conduction model including seepage are set up firstly. Then the characteristics of temperature field of those models are presented. Analytic solutions are given and e

5、quations are transferred with different variables. The position of the concentrated leakage passage and other parameters are treated as unknown parameters. Optimization target equations are established and different regression methods are referred to. At last，numerical experiments are carried out，an

6、d various regressions are done. The result indicates that this detecting method is very effective.Key words：hydraulic engineering；temperature field；numerical experiment；seepage ；regression analysis；optimization1 引 言堤防和土石坝破坏实例调查结果表明，由于渗透稳定性造成的破坏与失事，占整类事故的40%以上1。因而渗漏是危机堤坝安全的最主要隐患之一。确定渗漏发生的确切位置、集中渗漏通道的

7、大小和强度是对其有效治理的前提。因而开展渗漏定位的监(检 测技术研究很有必要。普遍存在于地质体中的温度场是一种天然物理场，渗流渗漏特征和热传导规律决定土体中温度场分布，是进行温度场示踪研究的基础。运用温度场第25卷 增2 王新建等. 堤坝集中渗漏温度场探测模型及数值试验 3795 对渗漏通道示踪具有成本低、效率高、环保等优点。堤坝温度获取的方法基本有4种：(1 采取勘探孔的水温代替堤坝温度，该方法往往应用在综合示踪法中，温度受孔内水自然对流换热影响，会改变地层原始温度分布，但是一般来讲由于钻孔直径不大，渗流和热交换不断进行，认为影响不大。测量工具是带有电缆的温度仪，每隔一段测量一个数据。(2

8、用一种管子楔入土体，在管子内部插入一连串的温度测试探头，待温度稳定后可读取数据。(3 在测点上用直径20 mm、长150 cm的钢钎打入并形成1 m深小孔，插入测温探头，封闭孔口，约过10 min待孔内温度平衡后读数。这种方法在饱和土的区域，同样是由孔内水体的温度代表土体。(4 应用光纤进行堤坝温度测量设置和数据处理方法，在温度数据收集上可以参考。目前堤坝渗漏温度场示踪法基本上是基于热源法2的定性分析，或主要作为其他示踪法的一种辅助手段3，4，或模型严重依赖其他方法的测试数据，对勘测孔有严格要求，参数难以确定3。尽管其原理较为明确，但真正、独立的温度场渗漏示踪模型还没有建立。国外相关研究是利用

9、测井水温通过作图法估计局部裂隙处地下水流速5，6以及地热开发中温度场数值方法正分析7，8。本文在建立经典热传导理论温度场模型和采用渗流模型的基础上，应用温度场反分析的方法确定渗漏位置。2 温度场反分析原理一般热传导物理问题的数学模型可以表达为方程的形式，即0( (=，x t x p A (1式中：A 为方程的算子；p (x 为方程参数向量结构体，可以是热传导系数、密度热、热容等，或者是边界上的温度、热源位置及其的组合等；x 为一维或者多维空间坐标。正问题的求解过程为正演，正演就是寻找一种映射关系M ，使得(，x t x p M = (2 在某些特殊工程实践过程中，需要通过测量土体内部某些温度分

10、布值t (x ， 来计算 (x p 中的全部或部分值。很显然，反演同样会出现正演中要考虑的适定问题。此外观测值一般不可避免地存在人为偶然误差和测量仪器的系统误差。这样观测值可以表示为( ( (，x x t x t +=(3式中： (，x t 为温度观测值； (，x t 为温度真实值； (，x 为温度误差。此问题通过解方程往往不能得出精确解，数值分析手段中回归分析或优化是求其最优解的有效方法。优化的目标是使某种残余量达到最小，非线性最优化目标9，即d | ( ( | (min ，x t x p M x p R = (4|可以根据不同要求取1，2L ，等范数，一般范数取2。通过Tikhonv 正则

11、化求解，正则解为 =(inf (x p R x p x p 202(| ( (|d | ( (|inf x p x p x t x p M x p +，(5式中：为正则化参数，20| ( (|x p x p 为正则化范函， (0x p 为参数初始值。将 (x p M 在初始点进行一阶Taylor 展开，可得 d ( ( (000p p pM x t x p M p p +=， (6将式(6代入式(5，使目标范函变分为0，并用0p Fr 表示算子M 在0p 处的Frechet 导数，于是有d d ( 00p p FrFr p pd ( ( 0，x t x t Fr p (7将式(7离散化，可得梯

12、度正则迭代公式为=+=+(d (d 1( 1( ( 1( (T ( (T (n n n n n n n n p M t p p p t t p K I K K (8 式中： (n K 为第n 次迭代步时离散后的导数矩阵， (n T 为第n 次迭代步时的计算值矢量，I 为单位矩阵。3 建立数学模型堤坝集中渗漏通道的样式主要有两大类：一类是面状或似面状，另一类是线柱状。本文主要研究线柱状渗漏通道。模型基本假设如下： 3796 岩石力学与工程学报 2006年 (1 渗漏通道充满流体，土体为饱和、各向同性导热介质。(2 在渗漏通道横截面的径向，水流自身对流换热快速，流体本身不存在温度差异。(3 渗漏对

13、温度场的影响范围和堤坝尺寸相比很小，不考虑源汇效应。 3.1 无渗流的稳定探测模型(1 在渗漏发生一段时间后，渗流和温度场达到稳定状态。(2 初始温度场均匀(一般情况下，由于热传导方程是线性的，通过正常地温梯度校正后，可以认为初始温度场均匀，或者采用最低温度平面法10计算，校正前后相差不大 ，根据初始等温均匀介质热传导理论，可得热传导方程11为0d d 1d d 22=+xtx x t (9 其一般解为21ln C x C t += (10任意温度测点T (x 0，y i 的温度、距离符合式(10，即得22201ln C y x C t i i += (i = 1，2L ，(11 若渗漏通道的

14、埋深为d ，则式(11变为22201ln C y d x C t i i +=）（ (12 式中：C 1，C 2为与边界条件和介质物性有关的参数11。当利用单线温度曲线做反分析时，把x 作为常量，目的是确定温度同竖直距离y 关系式中的参量。根据式(11可知，在某一竖直探测线上，即x 0固定时，当y = d 时温度最小，从曲线上可以判断温度异常最大的点就是渗漏所在的高程(不考虑垂向渗流和地温 。若先确定渗漏通道的埋深，则可有效地减少计算量。3.2 考虑渗流的探测模型假设传热过程由土颗粒热传导和渗流换热共同完成的，忽略流体的动能和耗散。建立有渗流时多孔介质单位土体的能量守衡方程为+y z y x

15、t k y x z y x t k x y x ( (， += ( (，z y x t cz z y x t k z z ztc v y t c v x t c v z y x +w w w w ww (13 式中：x k ，y k ，z k 分别为x ，y ，z 方向上的热传导系数；w w c 为水的体积比热。c为饱水多孔介质的体积比热。在热传导系数为常数k (不随位置和温度而改变 的情况下，具有非线性多孔介质渗流的方程变为=+=z z y y x x z y x z y x J k J k J k c c v v v c c u u u c ka z y x t z y x t z y x

16、 t a s s s w w w w 2( ( (U U ， (14式中：k 为热传导系数均匀饱和土体热传导系数；v x ，v y ，v z 分别为x ，y ，z 方向上渗流速度分量；J x ，J y ，J z 分别为x ，y ，z 方向上的水力梯度；k s x ，k s y ，k s z 分别为x ，y ，z 方向上的渗透系数；堤坝温度场边界条件及初始条件如下： (1 勘探数据显示，只有地温作用下，土体沿深度方向产生的温度梯度基本为常数。其实质是在某一深度处热流密度以常数传入几何模型，并达到稳定后形成温度场；可以利用勘探资料中不受渗漏影响的温度梯度 /(n t k q =来估计其值。由于土体

17、的热传导系数较低，在瞬态分析时，还要考虑因地热引起的初始温度分布，例如y y t0| (=，图1为渗流及地温作用下的渗漏温度等值线，显然温度场的变形可导致不同距离的探测线的最低温度不再在同一个水平面上。图1 渗流及地温作用下的渗漏温度等值线(单位： Fig.1 Leakage isothermal affected by geothermal andseepage(unit：水平距离x /m深度/m第25卷 增2 王新建等. 堤坝集中渗漏温度场探测模型及数值试验 3797 (2 土坝上边界为第一类边界条件：=t const 或对流换热条件。(3 横向上可认为土坝是无限宽，并保持恒温，即0( 0

18、( 0( 00 0( 00 000( ( 000( (21121120000M z M D y D x t y t M M z ht ht y tk D y t y t q y tk D y t y t t ht y tk y x z y x q y tk y x z y x q xt k i ，；，；，或，；，；，+=+±=+=±=± (15以上模型仍然是一种简化模型，事实上大气气温、水库水温是周期变化的，同时有日变化和年变化等周期变化，对流换热系数也是随堤坝表面特征而各异的。目前，对如此复杂边值条件和方程很难求出解析解，下面在简化边值情况下给出特殊情况下的解析

19、解。3.3 简化边值后的解析解参照数学物理方程有关资料，建立可解模型方程和边界条件分别如下：T a t tT2=+U (16 00 0=+=t r ；，r ，01t t =； r ，0002 0q rtrk t t =； 为了便于求解，方程作如下处理：对于具有热传导和扩散作用的二微偏微分方程一般形式：021=+ct t b t b t t y x yy xx (18令V t zy x e+=，并代入式(18，再令其中的x ，y 的一阶偏导项和零次项的系数为0，可得=121，221=， (412221b b c +=，原方程化减为04(2221=+V b b V V yy xx ）。对于稳态问题

20、含有V b b (412221+项，其边界条件或初始条件应作同样转化。为了便于计算机运算，采用过余温度0t t =，有渗流时的稳态传热传质问题解的表达式12为a yx U k a x U k q y x x x 22exp 2 (220l +=， (19 式中：K 0(x 为零阶第二类修正贝塞尔函数。结合数据采集实际情况有ak l x U k a l x U k q k l x x 2 (2 (exp2 (22000l+±±=， (20 模型1的温度场特征与探测状况见图2。一般探测线的距离是不均匀的，而探测线竖向测试距离是均匀的，k 为竖向测试步长。当下一个探测线在第一个参

21、考线(最接近渗漏通道 的同一侧时用“+”号；位于渗漏通道的另一侧时，用“”号。图2 模型1的温度场特征与探测状况 Fig.2 Characteric of temperature field and detectingconditions4 实现反分析的主要方法(1 回归分析。由正分析所求解，经变换变量，把热源位置、发热量等作为未知量，探测温度数据作为已知量。根据测温曲线数据特征，把探测间距l i (见图2 和相对于渗漏点深度k i 作为已知条件，最大异常温度点作为渗漏通道的埋深。根据上节热传导模型传热特点分析及假设，渗流只沿发生在顺x 方向，热源在y 方向上不会发生偏移，为了减少参数回归或优

22、化过程的复杂性，增强稳定性和准确性，在数据校正后，深度可作为已知条件(一般取校正后最大异常温度所在深度作为渗漏通道位置 ，以减少目标函数回归变量数。渗漏探测非线性回归是基于一般非线性回归模(17x 0x i(x i ，y j (x i ，y j +1 l il i +1o 3798 岩石力学与工程学报 2006年 型的，可以表示为+= (，X f t (21式中：X 为可观察的独立测试位置随机向量；为待估的渗漏通道位置等参数向量；t 为独立观测变量温度，其均值依赖于X 与；为随机误差。函数形式f ( 为已知热传导方程的解，迭代过程中算法取Levenberg-Marquart 法。(2 转化为最

23、优问题求解。用模拟退火算法实现参数最优化。随着研究问题的复杂性增加，利用一般的拟线性回归分析，要得到全局最优解就很困难，模拟退火算法可以增加现实全局最优化的可能性。最小化的目标函数为=ni ni i i iX f t S 1122( (，=ni i i iX X t12221(3 数值最优化方法。无论是拟线性回归还是模拟退火算法，都是建立在正分析问题解析解的基础之上。非稳态温度场绝大多数对含有级数，其中又有超越方程的根的情况，要实现这样的回归分析比较困难，此时正分析和反分析都采用数值方法是较好选择，借助于ANASYS 对结构的优化可以实现。数值最优化是正、反分析中都采用数值方法，当变量较多和搜

24、索范围较大时，常见问题是全局最优化的实现、是否收敛及收敛时间长短。5 数值试验单线回归定方位与综合回归定位置相结合计算渗漏通道的确切位置。在反分析中，为了确定渗漏通道的大致位置(即确定回归模型中的正、负号 ，先对每条探测线数据进行回归，这样每个探测线距离渗漏通道的位置大致确定，通过距离和温度曲线特征就可以判断是否位于渗漏通道的同一侧。假如所有曲线都是正态分布，则渗流方向上不同侧的曲线1与2，3，4在方差和最值特征上明显不同，因而可以确定综合回归方程中的符号。模型中的最低温平面法和综合回归(所有数据参与回归 都是建立在单线分析的基础上。(1 数值试验模型1取有关数据10，由于土体的热性参数随成分

25、、孔隙度、饱和度而变化，工程应用中需要对研究对象进行针对性热参数试验。设计试验参数为：饱和岩土体热传导系数，热容等数值可见相关研究13。取砂土样孔隙率n = 0.396，k = 2.87 W/(m· ，比热容加权平均值为2.381 J/(g·K(土的比热容为1.050 J/(g·K ，水的比热容为4.200 J/(g·K ，密度为2.6×103 kg/m3。用ANSYS 做试验模型，设计外半径为100 m，外边界温度为26 ，渗漏通道0.1 m近似为线热源，在内径上作用热源强度为5.06 W/m2的热流密度并保持在22 ，其他模型参数见表1。布

26、置勘探线3条，距离渗漏通道分别为10，30和50 m。提取数据，3条探测线的试验数据图形见图3，用MA TLAB 做最小二乘(p = 2非线性回归分析，其初始值参见相关研究10的最低温平面应用方程解。根据曲线判断渗漏通道位于3个曲线的最低温度决定的平面，从而渗漏通道的高程(埋深 可以决定。针对在50 m处及最低温平面出现局部最优解，运用模拟退火算法可得全局最优解。表1 模型参数 Table1 Parameters of models介质比热容c /(J·(g·K 1热传导系数k /(W·(m· 1密度/(kg·m 3导温系数a /(m2

27、3;s 1渗流速度v /(m·s 1饱和土体2.381 2.870 2.6×103 4.636×107 1.7×107渗漏水4.200 0.597图3 数值试验模型1的试验数据图形Fig.3 Temperature curves of numerical experiment model 1模拟退火算法是使目标函数分别为 = (2101C C x S ，=+ni i iC y x C Y 12220150 40ln ）（= (002q x S ，=+±ni x x i a k l x U k a l x U k q Y 12200002 (2

28、(exp 2 y /m温度/第 25 卷 增2 王新建等. 堤坝集中渗漏温度场探测模型及数值试验 3799 取得最小值 n 为测量数据的次数，其基本原理 见相关研究14，其参数设置见表 2。 表2 Table 2 模拟退火算法的参数表 y/m Parameters used simulated annealing at algorithm 衰减参数 0.97 步长因子 0.02 初始温度 100 容差 1×10 8 马可夫链长度 10 000 模型 1 的回归分析结果见表 3(C1，C2 的真值分 别为 0.36 和 22.6；最低温平面法和综合回归法中的 x0 的真值为 10 m。

29、 表3 Table 3 部位(或方法 10 m 30 m 50 m 50 m(模拟退火 最低温平面法 综合回归法 x/m 模型 1 的回归分析结果 残差平方和 7.585 2×10 1.197 7×10 3.140 2×10 7.018 1×10 4.149 9×10 8.346 5×10 6 图4 Fig.4 x 方向有渗流的渗漏管道横截面温度等值线图 (单位： Regression analysis results of model 1 C1 0.361 6 0.365 8 0.412 5 0.361 5 0.368 8 0.36

30、2 3 C2 22.669 2 22.633 6 22.192 3 22.669 7 22.604 5 22.642 1 x0/m 9.984 0 30.284 7 55.148 1 49.992 8 10.442 9 10.041 2 Temperature isothermal at cross-section of leakage with seepage in x direction (unit： 5 4 6 5 3 由表 3 可知，温度探测线距离渗漏通道越近， 参数的回归越准确，精度越高，残差越小；同时， 参数的初始值对拟合精度有影响，非线性回归中， 最低温平面法参数回归精度最低，综

31、合回归法的残 差最大。 (2 数值试验模型 2 利用 ANSYS 平面单元 PLANE55 中的 MASS TRANSPORT 模拟稳态热传导中的渗流问题。在 x 负轴大部分区域(温度变化只集中在靠近渗漏通道 的很小范围内和 y 轴方向上温度变化不大，特别是 远离渗漏通道的地方，以半径为 0.1 个单位的圆等 效线热源；热源强度为5.06 W/m2。远场单元采用 PLANE110 并赋值为 26.0 。取圆域半径为 200 个 单位，一是突出受渗流影响的温度区域，二是尽量 减少边界对渗流温度场产生的影响；取 x 轴方向上 的渗流速度为 1.7×10 7 y/m 距离/m 110 m；

32、210 m；330 m；450 m 图5 Fig.5 有渗流时不同探测线温度曲线示意图 Temperature curves of detecting lines with seepage 等的前提下，在渗流方向上的数据受渗漏通道影响 大，预测渗漏通道的位置更为精确些。总的回归分 析结果见表 4(q1 的真值为0.28 W/m2；最低温平面 法和综合回归方法中 x0 的真值为 10 m。 表4 部位(或方法 10 m 10 m 30 m 50 m 最低温平面法 综合回归法 模型 2 计算结果 Table 4 Computation results of model 2 q1/(W·m

33、 2 0.2776 1 0.280 08 0.279 95 0.279 88 0.208 80 0.280 03 x0/m 10.132 1 9.994 3 29.994 0 49.913 0 11.126 4 9.997 1 残差平方和 7.610 4×10 9.210 3×10 8.337 0×10 5.697 7×10 3.171 7×10 3 4 4 4 4 m/s，其他参数取值见表 1， 同时取 x =10， 10， 和 50 m 时的数据进行分析。 30 x 方向有渗流的渗漏管道横截面等值线见图 4。图 5 给出了有渗流时不同探测线

34、温度曲线示意图。 渗流影响下造成 x 轴方向温度场不对称，与渗 漏通道距离相同的 x =10 m 和 x =10 m 的温度曲线 明显有差异，结合回归分析的结果可得出：距离相 由表 4 可知，在渗漏通道的同一侧，温度探测 线距离渗漏通道越近，单线参数的回归越准确，残 差越小；由于渗流引起的温度场不对称，同样距离 3800 岩石力学与工程学报 2006 年 的探测线参数回归精度和残差都不同，最低温平面 法精度不是很高，而且残差异常。 离，因而通过垂直距离还可以计算渗漏通道与勘探 线的夹角。4#6#钻孔和 1#，2#钻孔受次要渗漏通 道的影响很大，并且通过同位素测试发现其中的垂 向渗流也较大，因而

35、进行通道计算误差较大，这里 只计算该区域的主要渗漏通道两边的 4 个钻孔。计 算前利用温度曲线的直线部分对 7#，8#钻孔(2#，3# 钻孔的温度曲线没有揭示地温梯度的温度曲线进 行地温矫正，矫正温度斜率为 0.07 /m。8#钻孔 剔除高程为 68 m 处的异常点。利用最小二乘法建 立目标函数，采用无约束最优化方法。工程实例的 计算结果见表 5。 表5 钻孔号 7# 8# 3# 2# C1 0.119 4 0.109 2 2.845 8 0.095 4 6 工程应用 濛浬灰场 4 段坝体垂向上是由三期子坝构成 # 的，在上一级坝形成的库区中放满灰渣，待固结稳 定后，再在上面修筑子坝，填土前先

36、在灰渣表面干 彻 0.50.6 m 厚的块石，库区离第一级子坝坝踵 18.0 m 处沿坝轴线纵向布置了水平排渗滤管直 径为 250 mm 的塑料管。每隔 50 m 设集渗井一个， 再用直径为 400 mm 的混凝土管把集渗井的积水排 到后坡坝面排水沟。经二级子坝加高后，目前 4#坝 顶高程为 83.0 m，坝高为 35.6 m。 利用水力输送灰渣到该灰场深槽堆放后。采用 原有的坝体过滤固结排水系统排除积水，40 d 后发 现 4#坝体后坡混凝土排水管正上方 71.0 m 高程处漏 水，漏水孔直径为 30 cm，并出现塌方、沉陷，漏 水呈黄色， 带出大量黄泥， 后漏水孔直径不断扩大。 为查明渗漏

37、的原因及确定渗漏通道的位置，以 漏水点垂直坝轴线的剖面为中心布置了 2 排探测 孔，第 1 排位于一级平台以下 1.0 m，即 76.0 m 高 程处；第 2 排位于坝顶以下 2.0 m，即 81.0 m 高程 处。第 1 排探测孔有 3 个，孔距为 15 m；第 2 排探 测孔有 5 个，孔距为 20 和 30 m。2 排共 8 个探测 孔，排距为 18 m。最低温平面工程布置及温度等值 线如图 6 所示。 工程实例的计算结果 x0/m 7.024 2 14.000 0 2.868 6 7.042 4 残差平方和 0.311 5 0.312 2 1.047 1 1.488 1 夹角/(

38、76; 44.49 Table 5 Results of the case C2 19.943 8 20.027 0 21.840 5 23.218 5 41.36 参考以上计算分析，对该工程防渗治理，防渗 修复加固措施采用技术上成熟的充填灌浆技术处 理。考虑到实际渗漏通道的复杂性，在平面上，上 排灌浆孔布置在上排桩轴线上，范围以计算出的渗 漏通道位置为中心向两侧横向分别延伸 25 m；在下 排桩轴线上，以计算的渗漏通道为中心向左延至 2# 钻孔，向右延伸 15 m，充填灌浆深度为 6374 m； 在以上范围内，灌浆孔排距为 3 m，孔距为 2 m，共 布置 154 个灌浆孔。在灌浆竣工 12

39、 个月后灰坝放 灰运行，漏水点没有出现渗水现象，进行钻孔注水 试验，测试的渗透系数为 4.5×10 明显。 4 cm/s，防渗效果 7 结论与讨论 论 7.1 结 本文建立的模型和数值试验结果，充分说明了 图6 Fig.6 最低温平面工程布置及温度等值线图 利用温度场进行堤坝集中渗漏探测是合理有效的。 模型及反分析结果表明： (1 渗流微弱或者只沿纵向发生时，温度测试 数据距离渗漏通道越近，反分析结果的精度越高； 考虑渗流时，这一结论在渗漏通道的同一侧才成立。 (2 由于一般温度场的渗漏模型(方程和边界 的解是非线性的，以及初始值的取值不同，用不同 的回归分析方法和手段可能得到不同的

40、结果，甚至 Isothermal curves and engineering arrangement at the lowest temperature plane 运用最低温平面法判断位于渗漏通道同一侧的 钻孔，结果表明，4 7 钻孔和 1 ，2 钻孔位于同 一侧。由于没有该坝段的渗流资料，不考虑渗流影 响，计算得到的 x0 应该是钻孔到渗漏通道的垂直距 # # # # 第 25 卷 增2 王新建等. 堤坝集中渗漏温度场探测模型及数值试验 3801 出现局部极小值。 (3 地热强度严重影响有渗漏通道发生时形成 的温度场形状。 (4 总体上看，渗漏通道位置的回归精度要比 其他参数精确，简单模

41、型回归精度比复杂模型回归 精度高。 (5 工程实例计算结果表明了该模型是合理有 效的。 7.2 讨 论 7 6 5 段 堤 防 的 渗 漏 J. 岩 土 工 程 界 ， 2002 ， 5(12 ： 37 39.(Chen Jiansheng， Bo， Yu Chen Liang. Study on seepage in dike of Gaoshuihe shipyard utilizing temperature field of underwaterJ. Geotechnical Engineering World，2002，5(12：3739.(in Chinese Beckera M

42、W， Georgianb T， Ambrosea H. Estimating flow and flux of ground water discharge using water temperature and velocityJ. Journal of Hydrology，2004，296：221233. Shemin G. Estimation of groundwater velocity in localized fracture zones from well temperature profilesJ. Journal of Volcanology and Geothermal

43、Research，1998，84：93101. Cenk Y，Jeffrey D S，Simon J R. A transient two-dimensional finite volume model for the simulation of vertical U-tube ground heat exchangesJ. Ashrae，1985，24(2：465744. 8 Mei V C，Emersond C J. New approach for analysis of ground-coil design for applied heat ump systemJ. Ashrae，19

44、85，24(4： 1 2171 224. 9 肖庭延，于慎根，王彦飞. 反问题的数值法M. 北京：科学出版 社，2003.(Xiao Tingyan，Yu Shengen，Wang Yanfei. Numerical Method of Inverse ProblemM. Beijing：Science Press，2003.(in Chinese 10 王新建，陈建生，陈 亮. 天然示踪法在灰土坝渗漏探测中的 (1 至今为止，没有保证取得全局最优的统一 有效工具，开发研制面向温度探漏的回归优化工具 是有必要的。 (2 温度测试精度能否达到本文需要的精度， 目前温度测试最高精度是 0.01，而

45、本文模型的温度 数据在远离渗漏通道 50 m 的地方，大部分温度只有 千分度的差别，因而在导温系数不大的情况下温度 探测渗漏通道可能会受到取样距离的限制。 (3 由于采用的是数值分析的方法，因而不存 在人为偶然误差和仪器的系统误差，实际数据回归 或者优化的精度会偏低。 参考文献(References： 1 太沙基，泼 克. 工程实用土力学M. 蒋彭年译. 北京：水利电力 出 版 社 ， 1960.(Terzaghi K ， Peck R. Engineering Practical Soil MechanicsM. Translated by Jiang Pengnian. Beijing：Water Resources and Electric Power Press，1960.(in Chinese 2 董海洲， 陈建生. 利用孔中温度场分布确定堤坝渗透流速的热源法