_四小时学会modflow文档.ppt

矩阵数据文件编辑 外部矩阵数据文件输入 外部矩阵数据文件输入（覆盖） 加载矢量图形文件(如Autu_CAD 文件) 格元图例－易于看懂加载各种模块分布 PMPath主界面 阻滞因子及其含义 PMPath－设置“粒子”分布与数量 PMPath－ “粒子”迹线演示 排水沟模块包Drain(DRN) －源汇项强度与水位“有关”！ 1 可模拟排水沟、泉水溢出等功能（只出不进，类似二极管） 2 只能向“排水沟”排，不能由排水沟进入含水层 3 排水沟水力传导系数(Drain Hydraulic Conductance)描述阻力大小，近似模拟向 排沟溢出的通畅程度 4 排水沟排水高程 Elevation of the Drain： －－泉口高程； －－排水沟内水位; 4 包含信息 各应力期的：排水沟位置、 哪层、排水高程、 排水沟水力传导系数 Drain Hydraulic Conductance 在该格元所控制的“排水沟”中， 单位水头差所溢出的水量 Cd Elevation of the Drain 该格元所控制的“排水沟”， 平均排水高程Di,j,k 。 所谓水头差，是指格元地下水位 与该排水高程之差。 用逐格元输入法输入排水沟参数： Equivalent Hydraulic Conductivity 排水沟在该节点所控制 单位长度、单位水头差 所溢出的水量 K Elevation of the Drain 排水沟在该节点所控制 排水高程 Di,j,k 用折线输入法输入排水沟参数： 自动化程度高！ 在某地下水模型中，选用了“排水沟”模块包 问题：模拟程序运行后，发现排水沟出水太多 了，怎么办？ 解决方法 沟排水太通畅了，把系数Cd或K减小一点，把地下水憋一下。或者 沟的排水高程 D 太低了，高点试试 3 两种方法一起上，再试试！ 河流模块包River(RIV) _1 当地下水向河流排泄时 当河流向地下水渗漏时 河流模块包River(RIV) _2 1 可模拟河流渗漏与溢出 2 河底水力传导系数为等效阻力，近似模拟水流的通畅程度 3 河水位指河水水面高程 4 包含信息 各应力期的：河流位置、 河水高程、河底（等效）阻力层底面高程 河底阻力层等效水力传导系数 用逐格元输入法输入River模块参数： 当地下水向河流排泄时 当河流向地下水渗漏时 Criv－格元单位水头差 渗漏（溢出）量 用折线输入法输入排水沟参数： 自动化程度高！ 蒸腾模块包Evapotranspiration(EVT) 1 可近似模拟地下水浅埋带蒸发 与 植被根系吸收地下水蒸发功能 2 涉及参数： 蒸发面高程、极限蒸发深度、水面蒸发强度 一般用外部 场文件输入 通用水位边界模块包General_Head_Boundary(GHB) 1 可近似模拟变水位边界、与水头有关的变流量边界等 2 涉及参数： 边界外侧水位Hbi,j,k 、边界水力传导系数CHBi,j,k 应用举例 1 变水头边界 2 边界流量随降深变化 3 阻水构造（如阻水断层），边界流量 与内外水头差成正比 CHBijk－－单位“格元水位与外部水位”差“诱发”的流入（出）量 注意： 1 在定水头边界上加任何“应力包”都是徒劳的，是白增加的。相当于在大海中“抽水”或“注水”！ 2 在定水头边界附近抽水，可“袭夺”定水头边界“取之不尽的水源”，而边界水位永远不会下降，这种情况自然界少见。 总之，应慎用定水头边界！！ 4 推荐：若边界即能提出水位边界，又能出提流量边界，则应首选流量边界。否则当出现边界量失真时不易察觉。 四种模型层类型 * 《Modflow 四小时速成》 1 水流问题有限差分法 2 Modflow 3 Processing Modflow 4 PM Path 5 例题练习 主要参考书目： 1 地下水流的数学模型和数值法， 孙讷正，1981.5，地质出版社。 2 地下水动力学（第二版）， 薛禹群，1997.9，地质出版社。 3 地下水流动问题数值方法， 陈崇希等，1990.6，中国地质大学出版社。 4 地下水流数值模拟，1989.9,李俊亭，地质出版社 教学中涉及软件： 1 Processing Modflow 2 Surfer (V.6或更高) 3 Auto CAD 2000 《Modflow 4小时速成》 数值解与数值方法 数值解 用有限个离散点的函数（水头）值，近似代替函数（水头）的 连续分布，称数值解。数值解属于近似解，在地下水渗流计算 中，是最“精确”的（概化误差最小） 地下水计算误差 概化误差(要命的误差)、测量误差、截断误差等…… 数值方法