第二讲数值模拟基础及方法理论.pptVIP

3）显示差分方程的求解 计算各结点初始时刻水头值 利用差分方程计算各结点t1时刻水头值 利用边界条件计算边界结点水头值 重复2、3步，直到计算出拟计算的各个时刻的水头值 显式格式（续3） 第二十六页，共四十九页，2022年，8月28日 算例（续4） 在上述模型中，设L=1000米 取空间步长为200米，时间步长为0.25天，分别计算各节点各时刻的水头值。 第二十七页，共四十九页，2022年，8月28日 Time/day x=0 m x=200 m x=400 m x=600 m x=800 m x=1000 m 0 20 10 10 10 10 10 0.25 20 12.5 10 10 10 10 0.50 20 13.75 10.625 10 10 10 0.75 20 14.531 11.250 10.156 10 10 1.00 20 15.078 11.797 10.391 10.039 10 1.25 20 15.488 12.266 10.654 10.117 10 算例（续5） 第二十八页，共四十九页，2022年，8月28日 Time/day x=0 m x=200 m x=400 m x=600 m x=800 m x=1000 m 0 20 10 10 10 10 10 0.25 20 20 10 10 10 10 0.50 20 10 20 10 10 10 0.75 20 30 0 20 10 10 1.00 20 -10 10 1.25 20 10 算例（续6） 如果 ⊿t=1,则 第二十九页，共四十九页，2022年，8月28日 (6)算例：隐式格式 在上述模型中，设L=1000米 取空间步长为200米，时间步长为0.25天，用隐式差分格式计算各节点个时刻的水头值。 第三十页，共四十九页，2022年，8月28日 在这个例子中， 解：隐式格式一般方程为 于是有 根据初始条件得 根据边界条件得 第三十一页，共四十九页，2022年，8月28日 由初始条件和边界条件 由此解得t1时刻的水头值为 在上述方程中取 n=0,可以得到计算t1时刻水头值的方程 所以上述方程变成 同理，可计算t2时刻的水头值 第三十二页，共四十九页，2022年，8月28日 2.地下水数值模型建模步骤 第三十三页，共四十九页，2022年，8月28日 模拟步骤 建立概念模型 建立数学模型 数值方法及软件（编程） 参数 拟合模拟：模型校正与检验 参数敏感性分析 预测模拟 软件 第三十四页，共四十九页，2022年，8月28日 一、概念模型（模型概化） 根据详细的地形地貌、地质、水文地质、构造地质、水文地球化学、岩石矿物、水文、气象、工农业利用情况等 模拟的区域: 含水层类型: 潜水（无压）、承压、混合、多层 维数:一维、二维、三维 水流状态：稳定流/非稳定流、饱和流/非饱和流 介质状况: 均质和非均质/各向同性和各向异性 孔隙/裂隙/双重介质 流体的密度差 边界条件和初始条件 必要时需进行一系列的室内试验与野外试验, 以获取有关参数, 如渗透系数、弥散系数、分配系数、反应速率常数等。 第三十五页，共四十九页，2022年，8月28日 第一页，共四十九页，2022年，8月28日 地下水建模方法和步骤 1.求解地下水运动方程的数值方法 2.地下水数值模型建模步骤 3.建模所需要的基本资料 第二页，共四十九页，2022年，8月28日 绝大部分数学模型是无法用解析法求解的，数值化就是将数学模型转化为可解的数值模型。 1.数值方法 1.1有限差分法 1.2有限单元法 1.3分有限差分法 1.4半解析半数值法 1.5边界元法 第三页，共四十九页，2022年，8月28日 (1)有限差分法原理 (2)两种方法建立有限差分方程 (3)求解有限差分方程 (4)收敛性和稳定性概念 (5)算例 1.1有限差分法 第四页，共四十九页，2022年，8月28日 (1)有限差分法的基本原理 将连续的问题离散后求解： 方法一．以地下水流基本微分方程及其定解条件为基础， 在渗流区剖分基础上，用差商代替微商，将地下水流微分方程的求解转化为差分方程（代数方程）求解。 方法二．在渗流区剖分的基础上，直接由达西定律和水均衡原理，建立各个均衡区的水均衡方程，即差分方程。 矩形网格 多边形网格 第五页，共四十九页，2022年，8月28日 网格划分的基本类型 （1）先划格线，格点位于网格中心 均衡网格 节点网格 （2）先规定格点位置，再垂直平分两相邻结点的连线作格线，形成的网格即为水均衡区 第六页，共四十九页，2022年，8月28日 MODFLOW网格系统 第七页，共四十九页，2022年，8月28日 导数的有限差商近似 导数的定