改进BT法在平原区潜水系统水位预报中的应用研讨.pdfVIP

改进“BT’’法在平原区潜水系统水位 预报中的应用研究 张辉 王献坤左正金郭东兴 罗文金 (河南省地质调查院，郑州450007) 摘要：本文是针对潜水系统水位预报有限单元计算方法——“BT”法存在的问题，提出了改进方法——改进“BT” 法，该法是在第一类边界水位的修正中，考虑水位下1降与潜水蒸发减少的相互作用；该方法的关键是计算潜水蒸发 减少引起的第一类边界水位下降减少值Sc。本文给出了用数值离散及解析公式相结合计算S。值的方法、原理、过程， 并给出了“BT”法、改进“BT”法的计算实例及结果对比。 关键词：潜水系统第一类边界(r、，) “BT”法改进“BT”法。 1 前言 地下水不稳定流有限单元计算方法——“BT”法是应用“解析法来确定渗流计算区第一类边 界(11t)上各结点的变水头值(BT值)，然后再利用有限单元法来计算内结点及渗流计算区第二类 边界(F2)上的变水头值”的一种方法u‘。它既节省勘探费用又能解决复杂的计算问题。它在供水 水源地勘察及矿坑排水设计中的应用非常广泛。 但传统的“BT”法在平原区潜水系统中的应用存在较大误差，原因是在应用解析公式计算、修 正边界水位时没有考虑潜水蒸发减少(由于新增抽水量引起水位埋深增大)反过来对水位的影响。 简单的说，改进的“BT法”就是在用“BT”法进行第一类边界水位修正时考虑潜水蒸发减少 与边界水位下降的相互影响的一种数值计算计方法。 2改进“BT法的原理方法(无限潜水系统) 在豫中东部平原区(可概化为无限潜水系统)，数值法预报局部水源地开采后的水位变化情况时， 四周边界一般处理为变水头的第一类边界(r1)，但预报时必须先给出其水头值(BT值)。开采后 边界(r1)水头值H有两部分组成：初始水位(水源地没有开采时的水位)H0及水源地开采而引 起的水位下降值S，即： H=Ho—S ④ 在传统的“BT”法中，s值的计算一般由解释公式直接算出，一般不考虑由于水源地开采引起 的区域上潜水蒸发量的变化(当然在内结点的数值计算中，考虑了数值计算区内各量的变化)。在水 源地开采后，在一定区域内潜水位下降，从而引起潜水蒸发大量减少，使实际水位下降值远小于“BT” 法中解析公式的直接计算值，即： S=SB—S。 ② 式中： S——改进“BT”法计算水位修正值； st；——“BT”法计算的水位修正值(单井公式迭加)； 作者简介：张辉，1966年生，hnzmdlwj@vip．sina．com 485 s。——水源地开采后潜水蒸发减少而引起的水位下降减少值。 由上式知，要获得s，最关键也是最复杂的是计算Se。第一类边界水位修正以后，就可以用常 规的有限元法计算各有关结点的水位。 3 se的计算 6 3．1开采前后潜水蒸发强度的变化A 设水源地开采前区域天然平均水位埋深为hm0，开采后为hm，根据阿维杨诺夫公式有： 厂E0=D·￡oo(1-hmo／Lo)“ ③ Le=p·eoo(1一hm／Lo)“ ④ 由③、④可得： S ⑤ 她≈。12面￡。 式中： n——因土质及气候而异的指数，河南中东部平原区一般取n=l； ￡0、E——开采前后的潜水蒸发强度(平均)； E“厂一水面蒸发强度； p——作物影响系数； A￡——潜水蒸发强度减少量(开采引起)； s——水位降深(开采后)，S=hm—hm0 Lc广一潜水蒸发极限深度，dh=L0一hm0； 3．2离散 E越小，时间越长A￡越