地下水数值模拟讲义.pptVIP

地下水数值模拟;;地下水的现状;地下水位降落漏斗状况;;;;各地区均出现了不同程度的地下水降落漏斗;对地下水位的数值模拟;地下水层;;1.3　模拟模型的数值解法;1.3.4　系统的输入与输出 1) 降雨量及其保证率 研究区地处平原区, 包气带由粉细砂, 亚砂土、亚粘土等细粒物质组成, 而且植被覆盖程度高, 降雨输入是最重要的地下水补给量. 降水量保证程度选用皮尔逊(Pearson) ? 型概率密度曲线来描述. 2) 其它参数值 浅层地下水系统的蒸发强度采用郑州水均衡试验场的实测值; 越流补给参数引用了水文地质条件相似的邻区资料; 地表水、灌溉水的入渗补给量和地下水开采量均为统计实测数据. 1.3.5　模拟模型的离散 通过越流项将模型é 和ê 耦合, 构成该区地下水水位分布的数学模型. 用不规则网络有限差分方法,以每个节点为中心, 由三角形外心连线组成每个节点的小均衡域. 在均衡域内对数学模型中泛定方程和定解条件进行离散, 用差商代替微商, 求得每个计算节点的差分方程, 形成一系列非线性方程组.;1.3.6　差分方程解法 方程(1) 和(2) 为大型非线性方程组, 不能直接求解, 同时为了节省计算内存和缩短计算时间, 本次计算采用超松驰迭代法解方程组(1) 和(2). 迭代时增加一松驰因子W , 构成一个收敛速度较快的近似解序列, 1.3.7　模拟模型的识别 模型识别以1994 年2 月25 日全区观测网统测水位作为初始流场, 以1994 年8 月31 和1995 年2 月25 日两次实测水位作为对比流场. 识别方法是正演拟合水位间接校正参数. 识别后确定13 个参数区, 其中浅层地下水5 个参数区, 弱透水层和中深层地下水各4 个参数区. 识别时段内地下水长期观测网117 个井点实测水位与计算水位值的绝对误差小于015m 的在第一时段为70% , 在第二时段为71% , 均符合国家标准(BG/T 14497- 93) 要求. 1.3.8　数学模型的验证 模型验证时段为1995 年2 月25 日至7 月30 日. 该时段的降雨量采用郑州市气象局气象公报数据, 地下水开采量采用郑州市节水供水动态监测报告的统计数据. 其它源汇项的选取同模型识别阶段. 验证结果表明水位拟合相对误差小于5% 的结点占观测水位结点总数的80% 以上, 符合国家标准(BG/T14497-93).;　模拟模型的运行预报 浅层地下水系统1996 年7 月30 日、1997 年7 月30 和1998 年7 月预报流场的水位分布与初始流场相比, 其展布趋势基本一致. 表明决定浅层地下水系统动态变化的各要素在预报时段内相对稳定. 中深层地下水系统的预报流场表明, 若按现有的开采布局和规划开采量, 1996 年7 月30 日降落漏斗中心将出现小于45m 的水位分布区(8. 96km2 ) , 1997 年7 月30 日降落漏斗中心将出现低于40m 的水位分布区(3. 02km 2) ; 1998 年7 月30 日水位降落漏斗中心将出现近8km 2 的40m 水位分布区. 漏斗中心区地下水位平均降速3. 5m/a 　流场宏观调控方案 根据模拟预报结果, 拟定研究区地下水流场调控方案. 1) 调整地下水开采量的时空分布, 改变中深层地下水系统开采井、井深和开采时间三集中的现状. 2) 增加浅层和深层地下水系统( 250 米) 的开采量,实行浅、中、深地下水系统统一调蓄、统一调控. 3) 改变地表水在城市供水中的份额, 逐步加大其供水比例.4) 实行分层分质供水方案, 大幅度提高优质水供水价格, 以利于优质水优用. 5) 统一规化, 统一调度和统一管理城市地表水、地下水和污水资源, 走可持续发展之路. ;;;地下水污染状况;模拟三维地下水溶质运移过程;地下水流和溶质的运移模拟;;;;;;海水向地下含水层入侵;　3 数值模拟在地下水污染研究中的应用　　　　数值模拟可应用于海水入侵、石油类污染物泄漏、垃圾填埋场污染物迁移、核素迁移等众多领域的模拟研究 。 3.1 海水入侵数值模拟　　在滨海和海岛地区，由于地下水的超量开采，以及持续性的气候干旱或海平面上升等自然因素而引起的海水入侵正成为沿海地区特别是滨海平原区普遍面临的重大地质灾害。海水入侵研究经历了从理想假定到合理概化，从室内实验模型、理想模型到数值模型这一发展过程。数值方法已成为模拟和求解海水入侵问题的最有力工具。概括起来，研究海水入侵的模型按研究对象可分为突变界面模型与基于海水-淡水以弥散带接触的过渡带模型。海水入侵的过渡带模型成为主要的研究方向。;;;海水入侵使灌溉地下水变咸;;;;;;;;;地下水将核素转移;;