互易定理在水力层析法求解水文地质参数中的应用.pdfVIP

互易定理在水力层析法求解水文地质参数中的应用 １ １ 1 王文梅 ， 孙蓉琳 ，赵贝 １．中国地质大学 环境学院 武汉４３００７ 摘要：基于水力层析法（hydraulic tomography ）与互易定理理论，研究抽水实验组数中互易 定理存在，及对参数反演的影响。得到以下结论：（1）通过抽水实验，运用SLE方法可以有效 的估计参数的空间分布特征；（2 ）对于给定数量的井，通过增加抽水实验组数可以有效的提高 参数反演的效果；（3 ）HT方法抽水实验互易定理的存在，降低了抽水组数增加对参数反演的提 升速率；（4 ）适当的抽水组数可以有效的提高参数反演的效果，同时可以减少计算机的运行量， 对于实际工程而言，可以降低抽水实验的花费，减少工作量，取得经济实惠、有效的参数反演 方法。 关键词：水力层析法 互易定理 抽水组数 1 引言 linear estimator ）算法。SLE方法源于Yeh等[1]进 含水层水力传导系数K与储水系数S是地 行 似 协 克 里 金 迭 代 算 法 (Iterative 下水渗流、溶质运移模拟预测中的重要参数。 cokriging-like technique) 。即通过已知的K 、S 对于非均质含水层，为了有效刻画K与S的空 和观测点水头值H ，通过似协克里金迭代算法 * * 间分布，T-C.J.Yeh等人 [1-3]研究出了可以刻画 调整参数，得到估计的K 、S 和估计水头值 非均质含水层的新方法:水力层析法（hydraulic H* ，通过反复的调参迭代，最后使得估计值与 tomography ,HT ），即有序的单井抽水，多井 观测值达到精度要求或达到最大迭代步长而 观测抽水实验。通过各种的数值模拟、室内实 停止运算。从HT法的原理与核心算法分析知： 验、野外试验场实验[4-9]表明：HT法可以有效 抽水组数越多，获取的观测水头值H越多，似 的刻画含水层的非均质性。 协克里金迭代算法中计算量增大，参数反演速 HT法的原理为：通过止水器把完全设有 率减慢。 过滤器的水井分隔为许多垂直井段，在其中一 Bruggeman[10-13] 等人通过格林公式验证 个井段抽水，在其他井段中观测水位降深值， 了孔隙含水层中存在互易定理，即交换抽水井 从而得到一组抽水量下，其他井段的观测降深 与观测井位置，交换前后观测井的降深值H将 数据。通过顺序地交换抽水井位置，就可以获 不变。由于HT方法基于顺序交换抽水井与观 得更多的水位降深数据集。即抽水组数越多， 测井位置，存在互易现象，所以增加抽水组数， 获取的水位降深数据集越大，揭露含水层信息 获得新的观测井数据将越来越少，参数反演提 越多。HT 法的核心算法是SLE （successive 高效果将越来越小。本文基于互易定理研究有 基金项目：水力层析法刻画含水层非均质性的试验优化技术研究 作者简介：王文梅（1986－），女，重庆市人，硕士研究生，主要从事地下水量预报，地下水溶质运移和渗流的 研究，（电话 （电子信箱） wwmwjs @126.com 通讯作者：孙蓉琳（1979－），女， 湖北武汉人， 副教授，主要从事地下水资源管理和裂隙岩体地下水渗流模 拟研究，（电话 （电子信箱）likeoasis@163.com 效的抽水组数，减少参数反演的计算量。