地下水水流模型建立过程.doc

地下水水流现状模型建立 1.模拟时间可长可短，不影响水流模型过程，一般用非稳定流，溶质运移考虑稳定流。单位里只需变更渗透系数单位m/s变为m/d，模型已运行完需要修改运移时间时，点主菜单F10,点设置到编辑引擎，可修改运移时间。 2.在给定模型底图时，先确定画好好模型的边界，埋深线，渗透系数等参数分区线等，以便后期人为好分区。 3.模型的边界零流量边界有：天然断裂带、天然基底隆起阻水带及人为流畅零流量变为。定水头边界有：泉沟河及水位变化很微弱的等水位线。给定水头边界（变水头边界）：按上下游等水位线给定一条弧线，或者根据补给边界断面给定一条直线。 4. 底图校正时，原点坐标输入左下角坐标。角点坐标输入左下角和右上角坐标，角度为0。 5.导入地表高程和基地高程时，采用模型坐标，单位为米。 6.生成网格后，将模型区外围采用无效水流区多边形概化，无效区不参与计算，流入流出量外也概化为无效区，给定水头后水头边界模型会给定水流量参与计算，其他边界为零流量边界。模型无效概化前全是有效网格，因此不能采用有效网格多边形，只能采用无效网格多边形进行无效区，无效区可用有效多边形修改，可按F9进行无效区可见进行视图可视化。 7.网格菜单下的绘制等值线可绘制出模型地表高程、厚度及基底等高线。 8.导入抽水井时要注意滤水管的顶底高程，开采时段及开采量，概化井的开采量和总量要一致。添加水位观测井时滤水管的高程为滤水管的中点高程，没有顶底高程。 9.给定渗透系数时电脑可按井渗透系数自动分区，比较分散，最好是人为划定多边形区域赋值，最后好调整参数。调整参数时只需要点数据库进行调整。 10.存贮参数一般不分区，给水度0.1—0.26之间，有效孔隙度0.25左右，总孔隙度0.3左右，后两参数对模型影响不大。调整参数时只需要点数据库进行调整。 11.依据统测水位导入水位标高，生成初始水头，最好有年初（模型开始期）水位作为初始水头，年末（模型结束期）水头作为与模型运行至365天时长流畅做对比验证。 12.定水头边界：一是泉沟河边界，水头为地面标高，可按线性梯度给值，起始点和结束点一条线的两个端点给定值；二是水位变化很小的等水位线曲线，给定一个水头值，三是最常用的给定水头边界（变水头），可按等水位线给一个值，也可按边界补给断面线边界性梯度给值，最好按月地下水位动态赋值，按回车键可给时间段。调整定水头时点击编辑，出现对话框后调整那个水头就用鼠标点击那个水头边界，一条线为组，线性梯度时始末点可点击个体进行修改。 13.河流赋值时主要是河流入渗量，按照河流整体在河床入渗考虑，河床厚度为包气带，渗透系数为垂向包气带入渗系数，按时间段河流来水量给定始末河流长宽，开始和结束可按比例给定，然后按照运行后的入渗量调整河流长宽，最好其他量给定，调整河宽一个量好调整。 14.渠系入渗量按照补给强度菜单给定，年mm数和时间段给，依据模型计算结果再进行调整。 15，降雨、凝结及灌溉入渗均按照补给强度菜单给定，mm数按照报告中均衡计算的量除以各多边形水位埋深分区的面积计算，再依据模型计算结果再进行调整。所有用补给强度给定的量计算结果为一个合并的补给量。调整量时点击编辑，可从1区到4或多区进行参数调整。 16.蒸发蒸腾按照mm数按照报告中均衡计算的量除以各多边形水位埋深分区的面积计算，再依据模型计算结果再进行调整。 17.所有参数给定后，进行模型运行，返回主菜单，点击运行，点击菜单MODFIOW2000,选择设置模拟层，潜水选择非承压层，选择干湿交替，点击有效网格单元湿润，若是预测模型，初始水头可给与预测期开始的水头。点击运行，选择MODFIOW2000和水均衡，点击编译和运行。 18.模型运行完点击输出，可根据不同时间看到水头和降深的运行结果，点击各菜单下的选项和时间可进行设置。查看水均衡菜单下的质量平衡的结束时间，有不同的时间的，看计算量和模型量的差别，进行参数调整。待调到基本一致时可进行观测孔和末期水位验证。 19.点击上部菜单中的图表下的时间序列下的水位，可看运行观测孔水头与实际水头拟合曲线，再选择直线相关系数图表，点击显示全部时刻，点击可查看相关系数及误差情况。点击上部菜单中的图菜单可绘制多种等值线，不想看到的或调整覆盖顺序的可按F9进行调整。 20.按F2进行3D模型演示，包括模型框架视图（最上面的英文）、输入、输出等。可分菜单点击设置，配合F9和调整出不同颜色和覆盖层等，主要包括坐标轴参数、底图透明性、网格、输入输出中各项3D图，潜水面及等值线等等。 21.预测模型是，把增加的开采量或者开采井在原有开采井表格中添加，并加入20年等预测时段和开采量，导入复制后的预测模型中，运行后查看降深和格时期的水头，降深色差值菜单可给出降深颜色效果。进行地下水水位影响预测。