地下水流数值模拟基础.pptVIP

地下水流数值模拟基础 中国地质大学环境学院 2009.6 本章主要内容 掌握有限差分法原理 导数的有限差分近似 建立有限差分方程 求解有限差分方程 收敛性和稳定性概念 截断误差 无压含水层水流问题 二类边界处理 抽水井计算水位与实际水位差别及影响因素分析 对于偏导数（偏微商），类似可以得到相应的差商，例如在地下水流微分方程中，常见的偏微商 3-2 一维承压水流问题有限差分法 例：河间地块承压水流模型 设两条河流平行、完全切割含水层，含水层等厚、均质各向同性。 步骤： （1）基础资料的分析 （2）概念模型 （3）数学模型 （4）数值方法及计算机程序 （5）参数 （6）结果分析 建立数学模型 （4）差分方程及其解法—显式格式 ①将（0—L）分成 N 等份， 3）显示差分方程的求解 计算各结点初始时刻水头值 利用差分方程计算t1时刻水头值 利用边界条件计算边界结点水头值 重复2、3步，直到计算出拟计算的各个时刻的水头值 4）差分方程的收敛性和稳定性 截断误差：用差商代替微商时，地下水流动方程产生的误差为截断误差。 收敛性：当空间步长和时间步长趋于0时，有限差分方程的精确解趋于地下水流动问题微分方程定解问题的精确解。则称该差分格式是收敛的。 稳定性：如果在求解差分方程过程中，某时间步引入某个误差，而在以后的各时段计算中，该误差不再扩大，则称该差分格式是稳定的。 5）算例 河间地块承压含水层地下水流动问题 设两条平行河流，完全切割承压含水层，含水层是均质、等厚、各向同性，隔水底板是水平。 ① 水文地质条件概化，建立概念模型 ② 建立坐标系，建立数学模型 ③ 确定参数 ④ 数值方法及计算 ⑤ 计算结果分析 5）算例（续1） 5）算例（续2） （5）差分方程及其解法—隐式格式 ①将（0—L）分成N 等份， 4）一维隐式差分方程的收敛性和稳定性条件 截断误差 隐式格式的收敛条件和稳定条件是： 5）算例 3）二维显示差分方程的求解 计算各结点初始时刻水头值 计算t1时刻水头值 计算边界结点水头值 重复2、3步，直到计算出拟计算的各个时刻的水头值 4）差分方程的收敛性和稳定性条件 例 一承压含水层，平面分布如图所示，地表高程为20m，隔水顶板底面高程为-20m，厚度 M 为20m，渗透系数K为8m/d，2008年9月18日测得的水位埋深为10m，南北边界为隔水边界，西部边界水位埋深保持为5米，东部边界水位埋深为10m，请（1）建立该含水层数学模型；（2）预测在含水层中部已2000方/天开采地下水，含水层水头分布动态。 （2）隐式格式 3）差分方程的求解 计算各结点初始时刻水头值 计算t1时刻水头值 计算边界结点水头值 重复2、3步，直到计算出拟计算的各个时刻的水头值 3-3 三维问题有限差分法 （1）显式格式 （1）显式格式（续1） 3）差分方程的求解 计算各结点初始时刻水头值 计算t1时刻水头值 计算边界结点水头值 重复2、3步，直到计算出拟计算的各个时刻的水头值 （1）显式格式（续3） （2）隐式格式 （2）隐式格式（续1） 3）差分方程的求解 计算各结点初始时刻水头值 计算t1时刻水头值 计算边界结点水头值 重复2、3步，直到计算出拟计算的各个时刻的水头值 （2）隐式格式（续3） 3-4 不等距矩形网格有限差分法 矩形网格有限差分法 例 如图所示，是一个承压含水层，含水层厚度为20米，渗透系数为8米/天，区域范围是一个矩形，x-方向长度为1000米，y-方向长度为1000米，储水系数为0.005，承压含水层顶板埋海拔高度-10米，与x-轴平行的边界为隔水边界，与y-轴平行的边界为水头已知的边界，并且左边界水头为海拔20米，右边界水头为海拔10米，试按图中网格系统建立地下水数值模型。 3-5 潜水（无压）水流有限差分法 潜水含水层平面二维流 3.5.1 显-隐式法 潜水含水层与承压含水层不同，随着潜水面的升降，会引起导水系数的变化，所以潜水含水层的导水系数是时间的函数。所谓凝固系数法就是将某时段的潜水含水层的导水系数利用该时段初的水头值计算，即令 3.5.2 全隐式法 3-6 源汇项的处理及井孔水头校正 3-6-1 越流、入渗和抽水井等问题的处理 如果考虑垂直渗流项（即源汇项），则二维承压流动微分方程可写成 建立差分方程时，在结点处应加上这一项，它可具体表示为 式中：为垂向渗流强度（包括越流、入渗或蒸发、井流等）；为相邻弱透水层垂向渗透系数；为相邻弱透水层厚度；为相邻含水层水头；为主含水层水头；为入渗强度或蒸发强度；为抽水或注水流量。 1） 网格剖分 将（0,b）分成Ny等份 ， (j=0,1,2，……Ny) ,记 取时间步长 （1）显式格式 (i=0,1,2，……Nx) （n=0，1，2，