数值模拟( 二 ) _数值

有限差分方程实际上是基本微分方程的近似表达式，其近似程度可用泰勒级数进行分析。通过微分方程的差分表达式，可以看出在利用差分格式代替微分式时，是存在误差的，即用有限差分方程组模拟地下水流系统会产生误差。
（三）对于边界条件和垂向水量交换的处理
不论是已知水头的一类边界或已知流量的二类边界，计算点落在边界上，该点就不需要列入均衡方程。垂向水量交换的处理也是如此，若点与抽水井重合，该点已列入均衡离散方程时，抽水量就直接参与该点所在均衡区的水均衡。
（四）均衡离散方程的解算
显然，在含水层参数和边界条件都给定的条件下，只要知道某时刻流场中所有点的水头值，就可计算出下个时间步长的所有点的水头。即在已知初始条件的基础上，可以计算不同时刻各点水头值、不同时刻的流场。对于这类问题的求解方法，从广泛使用微机处理的角度来看，超松弛迭代为许多研究者所采用。
（五）应用
综上所述，在已知初始条件、边界条件、垂向水量交换以及给定含水层参数的情况下，可计算渗流区内不同时刻、不同节点的水头值。当前，不论是在地下水资源评价的水量计算中，还是在矿山开采地下水的疏干计算或在因大面积地下水位下降引起的地质灾害防治中，数值法都得到了广泛应用。
目前有许多地下水数值法计算软件，适应性强、有较高的仿真性，广为采用，例如，MOP-FLOW（孔隙水三维有限差分法数值计算软件），GWMS-3D（二维或三维地下水流和污染物质运移数值模拟软件）等。
（六）实例
通过实例的学习，使同学们对用数值法求解过程有所了解。这个过程包括：①水文地质条件概化，建立概念模型；②根据水文地质概念模型，建立数值模型；③剖分计算区，整理计算资料；④校正数值模型；⑤验证数值模型；⑥运用模型进行预报。
实例位于太行山东麓冲洪积扇的交界处。含水层为第四纪松散层，上部为细砂和粉砂层，下部为砂卵砾石、粗砂砾石加土层、含粘土砾石层等。上部含水层地下水已被疏干，当前开采层埋深为40～80m，水位埋深多在10m以下，漏斗中心区已达30m 。边沿部分地区水位埋深为2～10m 。
1.水文地质概念模型
①含水层底板为隔水粘土层；②含水层主要为非均质各向同性的潜水含水层；③计算区的边界三面为已知水头的一类边界，另一面为不同程度的弱透水层，计算区面积近600km2；④区内有开采井；⑤地下水流为非稳定平面流，水流符合达西流。
2.数值模型
1）微分方程：
地下水动力学
地下水动力学
2）初始条件：H（x，y，0）=H0（x，y）
3）一类边界条件：H（x，y，t）｜Γ1=H1（x，y，t）
4）二类边界条件：
地下水动力学
地下水动力学
式中：W为汇源项，由降水入渗量和井的开采量代数和求出；n为内法线；其他符号同前。
3.剖分计算区并整理计算资料
将计算区剖分为506个小区、230个节点，其中第一类边界点40个，二类14个，取旱季为模型校正时段，给出10个分区参数并经过试验给出参数初值。
4.校正数值模型
【数值模拟】校正结果表明，微分方程和边界条件吻合。
5.验证数值模型
取雨季水位资料，分7个时段进行水位验证。根据验证资料绘制高低水位拟合图以及其他所需拟合图件，证明拟合程度良好，符合规范要求。
6.模型使用
利用验证过的符合实际的模拟模型，根据设计水位预计开采量，或根据设计的开采量预计不同时段的水位降低，尤其是漏斗中心的水位降低。

数值模拟( 二 )

猜你喜欢