_四小时学会Modflow.ppt

《Modflow 四小时速成》;主要参考书目： 1 地下水流的数学模型和数值法， 孙讷正，1981.5，地质出版社。 2 地下水动力学（第二版）， 薛禹群，1997.9，地质出版社。 3 地下水流动问题数值方法， 陈崇希等，1990.6，中国地质大学出版社。 4 地下水流数值模拟，1989.9,李俊亭，地质出版社 教学中涉及软件： 1 Processing Modflow 2 Surfer (V.6或更高) 3 Auto CAD 2000 ;数值解与数值方法;二维渗流问题矩形网格有限差剖分示意图;二维渗流问题矩形网格有限差剖分示意图;二维问题矩形网格有限差解之几何意义;三维问题的“简隐式”有限差分方程(Modflow用该格式) 1 不仅计算含水层水位，还要计算弱透水层中的水位 2 根据需要，含水层、相对隔水层还要进一步分成若干层;准三维问题的“简隐式”有限差分方程 1 仅求解各含水层（组）水位，弱透水层不计算 2 一般情况下，单一含水层（组）不再进一步分层(?);1 那些结（格）点需要求解 内格点 ==== 需要（列出差分方程式） 流量边界（第二类）格点 ====需要（列出差分方程式） 水位边界（第一类）格点==== 不需要（因为水位已知啦） 外格点 ==== 不需要 2 通用差分方程的特点 a 所有内格点的方程式都是相同的（给编写程序带来方便） b 当导水系数用调和平均值时（假设:计算区域之外的导水系数为零）， 隔水边界（零流量边界）的差分式与内结点相同。 c 流量边界“等价于”隔水边界 +“注(抽)水井”，注水井流量的大小 等于边界流入量。;1 概化水文地质条件－数学模型 2 离散化－数值模型 3 制备数据文件集合(BAS.dat,BCF.dat, Well.dat…) 4 运行模拟程序Modflow.exe 5 输出模拟结果文件集合 Heads.dat, Drawdown.dat… 6 解读模拟结果(均衡分析、绘图、表等) －－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－ 前处理软件－协助自动或半自动制作数据文件集合， 具有校对、图形显示、人机交互特征 后处理软件－具有自动整理结果数据文件为简洁易懂的图表 等功能;地下水模拟程序 美国地质调查局(UGS): MODFLOW 地下水流模拟程序 MT3D、MT3DMS…溶质运移模拟程序 集成前后处理程序 Visual ModFlow Processing ModFlow(PM) Groundwater Modeling System(GMS) Visual Modflow ……;;Modflow剖分特点与格元编号;Modflow格元均衡示意图;Modflow格元侧面流量计算表达式;垂向传导（越流）系数的计算;Modflow通用差分方程式;Modflow88模块化程序结构;递推求解模式;时间剖分问题 1 总模拟期分成若干“应力期” 2 每个“应力期”分解为若干时间段 3 每个应力期内各模块包参数为 “常数”，不同应力期间参数可以 变化，以实现随时间动态变化。;Processing Modflow集成前后处理系统介绍 File－文件菜单;Processing Modflow集成前后处理系统介绍 File/Preferences模拟程序设置;Processing Modflow集成前后处理系统介绍 Grid－剖分菜单;Processing Modflow集成前后处理系统介绍 Parameters－参数菜单;Processing Modflow集成前后处理系统介绍 Models－模型菜单;Processing Modflow集成前后处理系统介绍 Tools－前后处理工具菜单;;;基本模型－无任何外应力包;井模块包Well－源汇项强度与水位无关;“面状补给”模块包Recharge(RCH) －源汇项强度与水位无关;排水沟模块包Drain(DRN) －源汇项强度与水位“有关”！;;Drain Hydraulic Conductance 在该格元所控制的“排水沟”中， 单位水头差所溢出的水量 Cd;Equivalent Hydraulic Conducti