还剩83页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
GMS使用说明第1章概述该文档包含了:GMS及各组成模块的使用指南。第2章用TINS建立地面模型
2.1导入顶点文件
(1)__图标然后选择“FILE”菜单下的“OPEN”菜单项打开“GMS
3.1\tutorial\tins\verts.tins”文件;
(2)然后选择“BuildTIN”下的“Triangulate”菜单项;建立三角网格并自动利用Delauney进行优化避免病态矩阵;
(3)选择“Display”菜单中的“DisplayOptions”命令;选中“Contours”和“TINBoundary”复选框绘制地面等值线图;
(4)选择“Display”菜单中的“DisplayOptions”命令;取消“Contours”和“TINBoundary”复选框选中“________s”复选框按“OK”;然后选择“View”菜单下的“ObliqueView”命令再选择“Display”菜单中的“Shade”命令绘制带阴影的曲面图;编辑顶点:1选择“Display”菜单中的“DisplayOptions”命令;选中“Contours”和“Verti__s”复选框并单击“Contours”右边的的“Options”按钮;更改等值线间隔;2拖动顶点:从“ToolPalette”上选择“SelectVerti__s”工具选择“View”菜单下的“PlanView“命令然后用鼠标选中顶点拖动。选择“View”菜单下的“ObliqueView”命令也可以进行拖动。3用键盘进行精确__:选中顶点在工具条中的X、Y、Z坐标文本框中输入精确数值。锁定顶点:按下SHIFT键选择要锁定的顶点然后选择“ModifyTIN”菜单下的“Lock/UnlockVerti__s”命令。锁定了的顶点为另一种颜色。添加顶点:在“PlanView”模式下选择“CreateVertex”下面一组工具中的然后用鼠标单击来添加顶点。删除顶点:择“Edit”菜单下的“Delete”命令。然后利用“Create________”工具建立新的三角形。在“ModifyTIN”菜单下的“VertexOptions”命令在弹出的对话框中选择“Retriangulateafterdeleting”可以在删除顶点后自动连接新的三角形。平滑TIN文件:选择“Edit”菜单下的“DeleteAll”命令删除所有TIN然后打开比较稀疏的新的TIN文件然后选择“BuildTIN”菜单下的“TIN-ScatterPoints”命令在弹出的消息框中选择NO;选择“ModifyTIN”菜单下的“UniformlySubdivideTIN”命令并选择合适的系数;然后单击“2DScatterPoint”选择“Interpolation”菜单下的“toActiveTIN”命令输入新的文件名选中“__pelevations”复选框即可。可以通过选择“Edit”菜单下的“DeleteAll”命令删除所有TIN。在GMS中可以同时建立几个TIN通过__选择TIN图标可以选种TIN;双击TIN则在其上出现一个A说明该TIN激活;可以通过“Display”菜单下的“Hide”或“Show”命令隐藏或显示TIN。第三章立体地层模型TIN模块Borehole模块Solid模块__p模块建立立体模型需要三步:
1、读入钻孔数据:可以根据参考手册说明的文件格式输入钻孔数据;或者通过“Boreholes”菜单下的“BoreholeEditor”命令。读入数据文件:单击“”进入Borehole模块选择“File”菜单下的“Open”命令选择“tutorial\solids”下的“holes.gpr”文件。改变Z轴刻度:选择“View”菜单下的“Z__gnification”命令输入比例系数。显示钻孔名称:选择“Display”菜单中的“DisplayOptions”命令或按钮;选中“HoleNames”复选框。2、创建扩展多边形:1选择“PlanView”图标选择“Zoom”工具2选择“Display”菜单下的“DrawingGridOptions”命令将“Spacing”设置为
50.0选择“Snaptogrid”和“Displaygridlines”复选框将“Gridlinespacingincrement”设置为1“OK”。3定义边界:扩展多边形用__p模块中的特征对象工具输入。首先确定一个包围研究区的单个弧段:进入__p模块选择“CreateArc”工具然后用鼠标按顺时针或逆时针方向依次单击边界角点绘制研究区的边界在最后双击起始点封闭区域。4创建多边形:选择“FeatureO__ects”菜单下的“BuildPolygons”命令5选择“Display”菜单下的“DrawingGridOptions”命令取消“Snaptogrid”和“Displaygridlines”复选框将“Gridlinespacingincrement”设置为1“OK”。6选择“View”菜单下的“ObliqueView”命令。3、创建地表TIN
3.1选择接触点:单击“”进入Borehole模块选择“SelectContacts”工具在按下SHIFT的同时单击每个孔的顶部。
3.2创建TIN选择“Boreholes”菜单下的“Contacts-TIN”命令然后选择“Auto-extrapolateonly”选项OK。
3.3隐藏TIN转到TIN模式选择SelectTINs工具并单击来选择TIN选择Hide命令返回Borehole模块。
3.4创建绿色层TIN单击“”进入Borehole模块选择“SelectContacts”工具单击3G孔的绿/蓝界面。
3.
4.1自动选择接触面:选择“Boreholes”菜单下的“AutoSelect”命令
3.
4.2创建TIN选择“Boreholes”菜单下的“Contacts-TIN”命令然后选择“Auto-extrapolatewithtrimmingtoselectedboreholes”选项OK。;输入名字“topgreen”。
3.
4.2藏TIN转到TIN模式选择SelectTINs工具并单击来选择TIN选择Hide命令返回Borehole模块。
3.5创建红色层TIN单击“”进入Borehole模块选择“SelectContacts”工具单击4G孔的红/蓝界面;选择“Boreholes”菜单下的“AutoSelect”命令;关闭“__tch__terialabove”选项;选择“Boreholes”菜单下的“Contacts-TIN”命令然后选择“Auto-extrapolateonly”选项OK;输入名字“topred”。转到TIN模式选择SelectTINs工具并单击来选择TIN选择Hide命令返回Borehole模块。
3.6创建蓝色层TIN单击“”进入Borehole模块选择“SelectContacts”工具单击8G孔的红/蓝界面上红下蓝;选择“Boreholes”菜单下的“AutoSelect”命令;选择“Boreholes”菜单下的“Contacts-TIN”命令然后选择“Auto-extrapolatewithtrimmingtoselectedboreholes”选项OK;输入名字“topblue2”。转到TIN模式选择SelectTINs工具并单击来选择TIN选择Hide命令返回Borehole模块。
3.7创建底层TIN单击“”进入Borehole模块选择“SelectContacts”工具按下SHIFT键依次选择8G孔的蓝/红下6G的底面7G的蓝/红5G的蓝/红下2G的底面;选择“Boreholes”菜单下的“Contacts-TIN”命令然后选择“Auto-extrapolatewithtrimmingtoselectedboreholes”选项OK;输入名字“botblue2”。转到TIN模式选择SelectTINs工具并单击来选择TIN选择Hide命令返回Borehole模块。
3.12创建红色立体
3.
12.1创建立体选择名字为topred的TIN选择“BuildTIN”菜单下的“ExtrudeTIN-Solid”命令输入“red”在Elevation中输入-
45.0。
3.
12.2遮蔽立体选择Shade命令创建阴影;选择Refresh命令则返回线框图。
3.13创建兰色层选择名称为topblue2的TIN可以用“Edit”菜单下的“SelectfromList”命令来选择按下SHIFT键选择名称为“botblue2”的TIN选择“BuildTIN”菜单下的“FillBetweenTINs-Solid”命令输入“blue2”OK。
3.14去掉兰色层此时兰色层和红色层占据相同的空间要创建准确的地层模型必须在红色层中的兰色层位置做一个孔这可用在红色层中减去兰色层来实现。进入“Solid”模式选择SelectSolids工具选择名称为“red”的立体按下SHIFT键选择名为“blue2”的立体选择“Solids”菜单下的“SetOperations”命令在“Operand#1”框架中选则Delete选项选择减号输入名称。选择Shade命令创建阴影;选择Refresh命令则返回线框图。
3.15创建绿色立体转到TIN模式选择SelectTINs工具选择名为“topgreen”的TIN按下SHIFT键选择名为“topred”的TINS选择“BuildTIN”菜单下的“FillBetweenTINs-Solid”命令输入“green”OK。
3.16创建顶部兰色立体转到TIN模式选择SelectTINs工具选择名为“topblue1”的TIN按下SHIFT键选择名为“topgreen”的TINS及名为“topred”的TIN;选择“BuildTIN”菜单下的“FillBetweenTINs-Solid”命令输入“blue1”OK。
3.17显示立体图选择Shade命令
3.18切面
3.
18.1创建切面进入Solid模式选择CreateCrossSection工具选择PlanView
3.
18.2隐藏立体选择SelectSolids工具从“Edit”菜单中选择“SelectAll”选择Hide命令进入“ObliqueView”模式。
3.
18.3删除所有钻孔进入Boreholes模式从“Edit”菜单中选择“DeleteAll”。
3.
18.4删除多边形进入__p模式从“Edit”菜单中选择“DeleteAll”。
3.
18.5显示立体剖面图选择Shade命令
3.19层边界绿和兰色层边界在垂向上有突变而不是渐变。边界的明显尖灭可以通过在创建立体之前用TIN的编辑工具改变TIN顶点坐标来实现。
3.20删除立体和TIN从“File”菜单中选择“New”命令选择“No”。第4章2维地质统计在GMS中可以用2DScatterPoint模块进行2维插值。
4.3输入离散点集文件格式参见GMSFileFor__ts文件中的TabularScatterPointFile格式。进入2DScatterPoint模块从“File”菜单下选择“Import”命令选择“tutorial\geos2\plumedat.sp2”文件。
4.4改变显示模式进入DisplayOptions选中Datacolors选择Scatterpointsymbols左侧的按钮选择一种图形OKOK。颜色表示了污染物的浓度。从“Data”菜单下选择“ColorRampOptions”命令选中ShowColorLegend选项来显示图例。
4.5创建限定网格从“ScatterPoints”菜单下选择“BoundingGrid”命令输入X60、Y40方向上的数目OKOK。
4.6选择插值方式
4.7线性插值选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Linear”选项OK。选择“Interpolation”菜单下的“to2DGrid”命令输入数据名“Linear”OK。
4.8绘制等值线进入“2DGrid”模式进入DisplayOptions选择“Contours”和“Fringes”关掉“Nodes”OK。
4.9映射高程进入“ObliqueView”模式。选择“Data”菜单下的“__pElevations”命令选择名称为“linear”的数据集。
4.10阴影化网格选择Shade命令
4.11深谷插值进入2DScatterPoint模块选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Clough-Tocher”选项OK。选择“Interpolation”菜单下的“to2DGrid”命令输入数据名“ct”选中“__pElevations”可以将Z方向的数值直自动赋上OK。选择Shade命令
4.12简单IDW插值选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Inversedistan__weighted”选项选择“Inversedistan__weighted”选项右侧的“Options”的按钮在此对话框顶部的“Nodalfunction”框架中选择“Constant”在“Computationofinterpolationweights”框架中选择“Useallpoints”OKCLOSE。选择“Interpolation”菜单下的“to2DGrid”命令输入数据名“idw_const”OK。选择Shade命令
4.14带梯度面的IDW插值选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Inversedistan__weighted”选项选择“Inversedistan__weighted”选项右侧的“Options”的按钮在此对话框顶部的“Nodalfunction”框架中选择“Quadratic”OKCLOSE。选择“Interpolation”菜单下的“to2DGrid”命令输入数据名“idw_quad”OK。选择Shade命令
4.15截切由于插值的问题出现复数这是不合实际的。可以利用截切的方法进行修正。选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Truncatevalues”选项选择“Truncatetospecifiedrange”选项输入最小和最大值OK。选择“Interpolation”菜单下的“to2DGrid”命令输入数据名“idw_quad_trunc”OK。选择Shade命令
4.16普通邻近插值利用邻近点的拓扑关系进行插值当离散点聚集时效果好。选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Naturalnei___or”选项选则右侧的“Options”的按钮在此对话框顶部的“Boundingpseudo-points”框架中选择“Quadratic”在“Boundingpseudo-points”框架中取消“Extrapolatebeyondconvexhull”选项保证区外的浓度为0OKCLOSE。选择“Interpolation”菜单下的“to2DGrid”命令输入数据名“nn”OK。选择Shade命令
4.17Kriging插值选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Kriging”选项选则其右侧的“Options”的按钮。
4.
17.1创建实验方差图选择“EditVariograms”按钮在“VariogramEditor”对话框右上部的“Experimentalvariogrameditor”区的“New”按钮OK。
4.
17.2创建模拟方差图创建实验方差图后则需要创建模拟方差图。在“VariogramEditor”对话框右下部的“Nestedstructurespecification”区的“New”按钮在对话框左边名为“Modelfunctions”区中选择“Gaussian”选项在名为“Modelparameters”区中为“Contribution”输入
1965.0为“Range”输入
63.0OKCLOSECLOSE。选择“Interpolation”菜单下的“to2DGrid”命令输入数据名“krig”OK。选择Shade命令
4.18转换数据集进入“2DGrid”模式在“DataSet”组合框中选择“linear”项等值线的网格变了但是高程值仍然是前一数据集中的数据可以通过“__pElevations”命令修改。选择“PlanView”。用颜色填充等值线:进入DisplayOptions关掉“__lls”选择“Gridboundary”选择“Contours”右侧的“Options”按钮在对话框的左下角的“Contourmethod”区中选择“Colorfillbetweencontours”选项OK。
4.19利用数据计算器有时利用数据计算器DataCalculator比较插值所得的两组数据。从“Data”菜单中选择“DataCalculator”命令选择数据集输入表达式:在“Expression”字段中输入“absh-i”在“Result”字段中输入“Differen__”选择“Compute”按钮。选择“DataSetInfo”按钮查看新的数据集的统计信息。OKDONE。
4.20删除所有数据从“File”菜单中选择“New”命令选择“NO”。第5章3维地质统计
5.2所需模块和界面3DScatterPoint3DGrid
5.3输入离散点集文件格式参见GMSFileFor__ts文件中的TabularScatterPointFile格式。进入3DScatterPoint模块从“File”菜单下选择“Import”命令选择“tutorial\geos3d\tank.sp3”文件选择ObliqueView模式。
5.4显示数据颜色选择“Display”菜单下的“DisplayOptions”命令选中“Datacolors”选项OK。选择“Data”菜单下的“ColorRampOptions”命令选中“Showcolorlegend”选项OK。
5.5Z轴放大从“View”菜单中选择“Z__gnification”命令输入放大倍数。
5.6创建限界网格从“ScatterPoints”菜单下选择“BoundingGrid”命令输入X30、Y10、Z10方向上的数目并确认选中“Mesh__ntered”。OK。
5.7简单IDW插值选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Inversedistan__weighted”选项选择“Inversedistan__weighted”选项右侧的“Options”的按钮在此对话框顶部的“Nodalfunction”框架中选择“Constant”在“Computationofinterpolationweights”框架中选择“Usesubsetofpoints”在“Computationofinterpolationweights”区选择“Subset”按钮选择“Usenearest___points”并为点数输入64。OKCLOSE。选择“Interpolation”菜单下的“to3DGrid”命令输入数据名“idw_const”OK。
5.8显示等值面进入“3DGrid”模式选择“DisplayOptions”模式关掉“__lls”选项选中“Gridshell”、“Fringes”和“Isosu_____s”选项;选择“Isosu_____s”右边的“Options”按钮在第一个编辑区中输入
3000.0选中值1和值2之间的“Cap”选项。OKOK。选择“Display”菜单下的“ShadingOptions”命令选择“Hiddensu_____”选项选择“Uselightsour__”“__oothfeatures”和“Overlayedges”OK。选择Shade命令
5.9删除内部边从“Data”菜单中选择“Iso-Su_____Options”命令在对话框底部选中“Interioredgeremoval”。选择Shade命令
5.10边界指定范围改变显示选项以更清楚的显示颜色差异。选择“Data”菜单下的“FringeOptions”命令选择“Colorfringespecifiedrange”选项在“Minimumfringevalue”中输入3000在“__ximumfringevalue”中输入9000。选择Shade命令
5.11用Z轴选项进入3DScatterPoint模块选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令将“zscale”的值改为8。选择“Interpolation”菜单下的“to3DGrid”命令输入数据名“idw_const2”OK。选择Shade命令
5.12带梯度面的IDW插值选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Inversedistan__weighted”选项选择“Inversedistan__weighted”选项右侧的“Options”的按钮在此对话框顶部的“Nodalfunction”框架中选择“Gra___ntplane”OKCLOSE。选择“Interpolation”菜单下的“to3DGrid”命令输入数据名“idw_gp”OK。选择Shade命令
5.13带2次函数的IDW插值选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Inversedistan__weighted”选项选择“Inversedistan__weighted”选项右侧的“Options”的按钮在此对话框顶部的“Nodalfunction”框架中选择“Quadratic”在“Computationofnodalfunctioncoefficients”区中选择“Useallpoints”OKCLOSE。选择“Interpolation”菜单下的“to3DGrid”命令输入数据名“idw_quad”OK。选择Shade命令
5.14其他插值方式
5.15显示污染羽的截面进入“3DGrid”模式选择SideView命令选择CreateCrossSection工具切一截面然后进入ObliqueView模式选择DisplayOptions关掉“Iso-su_____s”OK。选择“Data”菜单下的“CrossSectionOptions”命令选择“Interioredgeremoval”和“Fringes”OK。选择“Data”菜单下的“FringeOptions”命令取消“Colorfringespecifiedrange”OK。选择Shade命令
4.16截切Truncation由于插值的问题出现复数这是不合实际的。可以利用截切的方法进行修正。进入3DScatterPoint选择“Interpolation”菜单下的“Interp.Optionscom__nd”命令选择“Truncatevalues”选项选择“Truncatetomin/__xofactivedataset”选项OK。选择“Interpolation”菜单下的“to3DGrid”命令输入数据名“idw_quad_trunc”OK。选择Shade命令
5.17建立__切面影像循环用FilmLoop建立污染羽的动画显示。
5.
17.1显示模式进入“3DGrid”模式选择SelectCrossSections工具选择切面选择“Edit”菜单下的“Delete”命令OK。选择“Data”菜单下的“ColorRampOptions”命令关掉“Showcolorlegend”选项OK。选择“Data”菜单下的“FringeOptions”命令选择“Colorfringespecifiedrange”选项为“Minimumfringevalue”输入1000为“__ximumfringevalue”输入
15000.0。OK。
5.
17.2设置影像循环选择“Data”菜单下的“FilmLoop”命令选择“Setup”按钮选择“Ani__tecuttingplane”选项选中“Z”选项OK。
5.
17.3回放影像循环选择“Play”按钮选择Stop停止动画。
5.18建立__等值线影像循环在“FilmLoop”对话框中选择“Setup”按钮关掉“Ani__tecuttingplane”选项选中“Ani__teiso-su_____”选项为“Beginvalue”输入1000为“Endvalue”输入15000。选择“Capabove”选项选择“Displayvalues”选项。OK。选择“Play”按钮选择Stop停止动画。选择“FilmLoop”对话框中的“Done”按钮退出。第6章MODFLOW—网格化在GMS中有两种方法可以构建MODFLOW模拟:网格化和概念模型法。前者直接利用3D网格在单个单元上使用源/汇项和其他模型参数;后者利用GMS中__p模块中的GIS工具建立概念模型模型中的数据拷贝到网格中。实例:__75000ft×75000ft;参数见图。
6.4创建网格进入“3DGrid”模式选择“Grid”菜单中的“CreateGrid”命令在名为“X-dimension”区域中的“Length”中输入75000在“Number__lls”中输入15;在名为“Y-dimension”区域中的“Length”中输入75000在“Number__lls”中输入15;在名为“Z-dimension”区域中的“Length”中输入15000在“Number__lls”中输入3;OK。
6.5MODFLOW模拟的初始化从“MODFLOW”菜单中选择“NewSimulation”命令。
6.6基本模块从“MODFLOW”菜单中选择“BasicPackage”命令。
6.
6.1标题首先输入标题。在头文件的第一行中输入:GMSMODFLOWTutorial–GridApproach;在第二行中输入你的名字和日期。
6.
6.2程序包单击“Packages”按钮在对话框中选择“Drain”、“Well”、“Recharge”模块在“Solver”区中选择“StronglyImplicitPro__dure”OK。
6.
6.3单位单击“Units”按钮为长度单位输入“ft”时间单位输入“d”OK。
6.
6.4IBOUND数组单击“IBOUND”按钮。IBOUND0:有效单元IBOUND=0:无效单元IBOUND0:定水头单元
6.
6.5初始水头单击“StartingHeads”按钮。
6.
6.6退出对话框CLOSE。
6.7直接将IBOUND值赋给单元利用“__llAttributes”命令可以直接将IBOUND值赋给单元首先选中单元最左侧两层的单元
6.
7.1显示最左列选择ViewtheJAxis
6.
7.2选择单元选择“Select__lls”工具拖拉鼠标选择。
6.
7.3改变IBOUND值从“MODFLOW”菜单中选择“__llAttributes”命令为“IBOUND”输入-1;OK。在“GraphicsWindow”外单击取消单元选择。选择“ViewtheKAxis”命令
6.
7.4检查值从“MODFLOW”菜单中选择“BasicPackage”命令。单击“IBOUND”按钮。
6.8BCF模块从“MODFLOW”菜单中选择“BCFPackage”命令。
6.
8.1层的类型
6.
8.2层的参数
6.
8.3顶层选择“TopElevation”按钮选择“Constant-Layer”按钮输入200OKOK。然后以同种方法输入:BottomElevation-150ft、HorizontalHydraulicConductivity50ft/d、VerticalHydraulicConductivity5ft/d。
6.
8.4中间层然后以同种方法输入:BottomElevation-400ft、HorizontalHydraulicConductivity3ft/d、VerticalHydraulicConductivity
0.5ft/d。
6.
8.5底层然后以同种方法输入:BottomElevation-700ft、HorizontalHydraulicConductivity7ft/d、VerticalHydraulicConductivity
1.5ft/d。Close。
6.9Recharge模块从“MODFLOW”菜单中的“Sour__/Sink”子菜单中选择“RechargePackage”命令选择“Constant-Array”按钮输入值
0.003OK。OK。
6.10Drain模块
6.
10.1选择单元选择“Select__lls”工具安下SHIFT键用鼠标选择。
6.
10.2指定Drain从“MODFLOW”菜单中的“PointSour__s/Sinks”命令在“Drain”区选择“Add”选项选择在“conductan__”中输入80000OK。
6.
10.3指定Drain高程从“MODFLOW”菜单中的“Sour__/Sink”子菜单中选择“DrainPackage”命令在“elev.”栏中输入高程值OK。
6.11井模块
6.
11.1顶层井选择“Select__lls”工具安下SHIFT键用鼠标选择。从“MODFLOW”菜单中的“PointSour__s/Sinks”命令在“Well”区选择“Add”选项选择在“Flow”中输入-432000OK。
6.
11.2中间层井单击左下侧的“”进入下一层。选择“Select__lls”工具安下SHIFT键用鼠标选择。从“MODFLOW”菜单中的“PointSour__s/Sinks”命令在“Well”区选择“Add”选项选择在“Flow”中输入-432000OK。
6.
11.3底层井选择“Select__lls”工具安下SHIFT键用鼠标选择。从“MODFLOW”菜单中的“PointSour__s/Sinks”命令在“Well”区选择“Add”选项选择在“Flow”中输入-432000OK。
6.12保存模拟从“MODFLOW”菜单中的“S__eAs”命令输入文件名。
6.13运行MODFLOW从“MODFLOW”菜单中的“RunMODFLOW”命令OK。
6.14查看结果从“MODFLOW”菜单中的“ReadSolution”命令OK。
6.
14.1改变层单击左下侧的“”进入下一层。单击左下侧的“”进入上一层。
6.
14.2颜色填充等值线从“Data”菜单中选择“ContourOptions”命令在对话框左下角许选择“Colorfillbetweencontours”。
6.
14.3显示颜色图例从“Data”菜单中选择“ColorRampOptions”命令在对话框左下角许选择“Showcolorlegend”。第7章利用概念模型建立MODFLOW
7.1问题描述
7.2启动GMS
7.3所需模块和界面3DGrid模块2DScatterPoint模块__p模块MODFLOW界面
7.4导入背景图
7.
4.1读入图像文件选择“File”菜单下的“Open”命令打开文件“tutorial\modf__p\start.gpr”单击Open按钮。
7.
4.2I__GE文件VSTIFF文件装入TIFF文件则要进行注册也就是选择I__GE上的三个点并且输入其相应的实际坐标。
7.5特征目标对象点、弧段、结点弧段两端的点、顶点弧段中间的点、面
7.6创建局部源/汇图层Coverage进入__p模块选择“FeatureO__ects”菜单下的“Coverages”命令将“Defaultelevation”改为700将默认名字改为“Sour__s/Sinks”将“Coveragetype”改为“MODF/MT3D/MODP”选择“Options”按钮检查图层类型是否是MODF/MT3Dlocalsour__s/sinksOK。
7.
6.1定义单位单击“Units”按钮为长度单位输入“ft”时间单位输入“d”OK。
7.
6.2定义边界选择“CreateArc”工具沿边界造弧退一点时按Backspa__键取消造弧时按ESC键最后在始点双击。
7.
6.3拷贝边界首先为每一层创建水力传导度图层:选择“FeatureO__ects”菜单下的“Coverages”命令选择“Copy”命令将新图层名改为“Layer1”选择“Options”按钮将“Coveragetype”改为“MODF/MT3D/MODPlayerattributes”OK。当图层Layer1coverage高亮度显示时再次按“Copy”按钮将新图层名改为“Layer2”选择“Options”按钮将“Coveragetype”改为“MODF/MT3D/MODPlayerattributes”在中间对话框中将“Assignfromlayer”和“tolaye”均改为2OK。然后创建补给多边形图层:选择“FeatureO__ects”菜单下的“Coverages”命令选择“Copy”命令将新图层名改为“Recharge”选择“Options”按钮将“Coveragetype”改为“MODF/MT3D/MODParealattributes”OK。最后:选择“Sour__s/Sinks”图层单击OK退出Coverages对话框。
7.
6.4定义指定水头弧段首先选择“SelectVerti__s”工具选择两个顶点选择“FeatureO__ects”菜单下的“Verti__s-Nod”命令;选择“SelectArcs”工具选择弧段选择“FeatureO__ects”菜单下的“Attributes”选择“SpecifiedHead”在底部对话框的第二个文本框内输入2单击OK。然后指定水头值为弧段两端点指定水头值假设其间的水头线性变化:选择“SelectPoints/Nodes”工具双击南部边界左侧的端点输入值
697.0OK。输入其他2个交点的水头值。
7.
6.5定义排水沟弧段选择“CreateArc”工具沿排水沟造弧退一点时按Backspa__键取消造弧时按ESC键最后在终点双击。选择“SelectArcs”工具选择弧段选择“FeatureO__ects”菜单下的“Attributes”选择“Drain”在文本框内输入
6000.0单击OK。然后指定排水沟的高程值为弧段两端点指定高程值假设其间的水头线性变化:选择“SelectPoints/Nodes”工具双击弧段的端点输入值OK。
7.
6.6创建多边形先选择弧段然后选择“FeatureO__ects”菜单下的“BuildPolygons”命令。然后选择“FeatureO__ects”菜单下的“DisplayOptions”命令选择PolygonsFill选项可以显示区。
7.
6.7创建井选择“CreatePoint”工具在井位处单击要准确定位则可以输入坐标值;选择“FeatureO__ects”菜单下的“Attributes”选择“Well”在文本框内输入开采量并确定开采层位单击OK。网格细化:选择“SelectPoints/Nodes”工具选择所有井选择“FeatureO__ects”菜单下的“Attributes”单击“RefinePoint”按钮选择“RefinegridinXdirection”选项为Base__llsize输入75为__x__llsize输入500;选择“RefinegridinYdirection”选项为Base__llsize输入75为__x__llsize输入500;OK。
7.7刻画补给区
7.
7.1选择Recharge图层从Coverage组合框中选择Recharge图层
7.
7.2创建垃圾填埋场边界选择“CreateArc”工具选择“SelectPoints/Nodes”工具选择端点输入准确的坐标值选择“SelectVerti__s”工具依次选择4个端点。
7.
7.3创建多边形先选择弧段然后选择“FeatureO__ects”菜单下的“BuildPolygons”命令。然后选择“FeatureO__ects”菜单下的“DisplayOptions”命令选择PolygonsFill选项可以显示区。
7.
7.4指定补给值选择SelectPolygons工具双击补给区输入补给强度;双击__区输入__区的补给强度。
7.8定义水力传导系数
7.
8.1顶层从Coverage组合框中选择Layer1图层。后选择“FeatureO__ects”菜单下的“BuildPolygons”命令。双击区输入水平和垂直渗透系数。
7.
8.2底层从Coverage组合框中选择Layer2图层。后选择“FeatureO__ects”菜单下的“BuildPolygons”命令。双击区输入水平和垂直渗透系数。从Coverage组合框中选择“Sour__s/Sinks”图层。
7.9装入网格边框选择“FeatureO__ects”菜单下的“GridFrame”命令选择“NewFrame”按钮然后用拖拉或输入数值确定边界。
7.10创建网格选择“FeatureO__ects”菜单下的“__p-3DGrid”命令输入数值。
7.11指定有效单元选择“FeatureO__ects”菜单下的“Activate__llsinCoverages”命令。
7.12初始化MODFLOW的值进入3DGrid模式选择“MODFLOW”菜单下的“NewSimulation”命令。
7.13转换概念模型转到__p模块选择“FeatureO__ects”菜单下的“__p-MODFLOW”命令;确定选中“Allapplicablecoverages”选项OK。
7.14层的高程插值
7.
14.1导入地面高程离散数据文件转到2DScatterPoint模式选择“File”菜单下的“Import”命令打开“tutorial\modf__p\terrain.sp2”文件。
7.
14.2计算初始水头数据集通过将地面高程减少10ft作为初始水头值。选择“Data”菜单下的“DataCalculator”命令在“Expression”字段中输入计算公式“a-10”在“Result”字段输入结果名“startinghead”单击“Compute”按钮单击“Done”按钮退出。
7.
14.3高程和水头插值选择“Interpolation”菜单下的“toMODFLOWLayers”命令对话框上部左侧为数据文件名右侧为其相应的用途MODFLOW中的数据名。选择数据集名及其相应的用途名单击“__p”按钮单击“OK”进行插值。
7.
14.4导入层底高程离散数据文件选择“File”菜单下的“Import”命令打开“tutorial\modf__p\elevs.sp2”文件。
7.
14.5层底高程插值选择“Interpolation”菜单下的“toMODFLOWLayers”命令对应关系已经默认。
7.
14.6调整显示选择SelectScatterPointSets工具在离散点符号周围拖拉一个矩形框选择Hide命令。关掉网格外框模拟边界的显示:转到__p模块选择“Display”菜单下的“DisplayOptions”命令关掉“Gridframe”选项。
7.
14.7显示模型剖面转到3DGrid模块选择靠近模型中心的一个单元选择ViewJAxis命令要更好的显示可以调整垂向比例选择“View”菜单下的“Z__gnification”命令输入放大值OK。通过由__的两个水平方向的箭头可以看不同列的剖面。
7.
14.8修正高程值选择“MODFLOW”菜单下的“CheckSimulation”命令单击“RunCheck”按钮选择“FixLayerErrors”按钮。则出现错误信息并提供了修正方法。选择“Truncatetobedrock”按钮单击“FixAffectedLayers”按钮OKDone。同样通过可以看水平切面。
7.15检查模拟在运行前应该运行MODFLOWModelChecker多次选择“MODFLOW”菜单下的“CheckSimulation”命令单击“RunCheck”按钮
7.16保存工程
7.17运行MODFLOW选择“MODFLOW”菜单下的“RunMODFLOW”命令。
7.18显示水头等值线选择“MODFLOW”菜单下的“ReadSolution”命令。改变等值线颜色:选择“Data”菜单下的“ContourOptions”命令单击“Contourcolor”按钮单击“Color”窗口选择颜色OK。单击右下部的按钮可以看第二层的水头等值线。
7.19显示侧面上的水位线在图中任选一点单击。
7.20查看流量均衡进入__p模块选择SelectArcs工具单击最右侧的排水沟弧段则通过该排水沟的总流量显示在底部倒数第二条带上。通过按下SHIFT键可以选择多条弧段显示单元的流量:转到3DGrid模块选择一些单元选择“Data”菜单下的“FlowBudget”命令。
7.21添加注释进入__p模块选择CreateText工具在添加注释的地方单击输入文本选择是否填充和字体颜色等。利用SelectDrawingO__ects工具可以选择注释并进行编辑修改。通过创建CreateLine线工具绘制线然后利用利用SelectDrawingO__ects工具可以选择该线双击在Arrowheads区内选择适当的箭头。第8章MODPATH
8.1问题描述
8.2开始
8.3所需模块
8.4打开工程选择“File”菜单下的“Open”命令打开“tutorial\modf__p\sample\sample.gpr”文件。
8.5初始化MODPATH模拟进入3DGrid模块选择“MODPATH”菜单下的“NewSimulation”命令。
8.6指定参数转到__p模块在GMS窗口中的Coverages组合框中选择“Layer1”选择SelectPolygons工具双击该层多边形选中“AquiferPorosity”选项并输入值OK。在GMS窗口中的Coverages组合框中选择“Layer2”选择SelectPolygons工具双击该层多边形选中“AquiferPorosity”选项并输入值OK。
8.
6.1为单元指定参数选择“FeatureO__ects”菜单下的“__p-MODPATH”命令OK。
8.7指定起始位置
8.
7.1选择单元进入3DGrid模块选择进入第二层选择Zoom工具选择Select__lls工具选择单元。
8.
7.2创建起始位置选择“MODPATH”菜单下的“GenerateParticles”命令选择“Distributestartingpointson__llfa__s”选项输入点数OK。保存显示:选择Frame工具;单击返回上一层。
8.8指定跟踪方向选择“MODPATH”菜单下的“GeneralOptions”命令在对话框的顶部选择“Backward”选项OK。
8.9保存模拟选择“MODPATH”菜单下的“S__eAs”命令输入保存的文件名“path
1.rsp”单击“OK”按钮。
8.10运行MODPATH选择“MODPATH”菜单下的“RunMODPATH”命令单击“S__e”按钮。
8.11查看结果选择“MODPATH”菜单下的“ReadSolution”命令。在图中任选一点单击可从侧面看。同样通过可以看水平切面。
8.12从垃圾场追踪颗粒
8.
12.1改变跟踪方向选择“MODPATH”菜单下的“GeneralOptions”命令在对话框的顶部选择“Forward”选项OK。
8.
12.2删除原来的颗粒选择SelectStartingLocations工具用鼠标拖拉选择示踪源按Backspa__键或选择“Edit”菜单下的“Delete”命令删除OKOK。
8.
12.3指定新的起始位置在GMS窗口中的Coverages组合框中选择“Recharge”选择Select__lls工具拖拉鼠标选择单元。选择“MODPATH”菜单下的“GenerateParticles”命令选择“Distributestartingpointsonw.t.su_____”选项OK。
8.13保存模拟选择“MODPATH”菜单下的“S__eAs”命令输入保存的文件名“path
2.rsp”单击“OK”按钮。
8.14运行MODPATH选择“MODPATH”菜单下的“RunMODPATH”命令单击“S__e”按钮。
8.15查看结果选择“MODPATH”菜单下的“ReadSolution”命令。第9章MT3DMS进行MT3DMS模拟之前必须先建立水流模型根据前面所介绍的MODFLOW使用说明。
9.5建立溶质运移模型
9.
5.1初始化模拟选择“MT3D”菜单下的“NewSimulation”命令。
9.
5.2基本运移模块选择“MT3D”菜单下的“BasicTransportPackage”命令。标题:输入注释说明及__姓名、日期等物种:选择“DefineSpecies”按钮单击“New”按钮给物种命名OK。模块:单击“Packages”按钮选中“AdvectionPackage”选项选中“DispersionPackage”选项选中“Sour__/SinkMixingPackage”选项OK。应力期:单击“StressPeriods”按钮将“Length”字段改为365模拟时间长度;保证Transportstepsize的值为00表示MT3DMS自动计算适宜的时间步长;OK。输出控制:单击“OutputControl”按钮选中“Printors__eatspecifiedinterval”选项指定间隔大小OK。浓度边界数组ICBUND:初始浓度数组:HTOP和厚度数组:从MODFLOW模型中自动获取。孔隙度数组:
9.
5.3对流模块选择“MT3D”菜单下的“AdvectionPackage”命令在“Solutionschem”下拉框中选择“ThirdorderTVDschemeULTI__TE”OK。
9.
5.4弥散模块选择“MT3D”菜单下的“DispersionPackage”命令单击“LongitudinalDispersivity”按钮选择“Constant-Grid”选项输入纵向弥散系数值OK;分别为“Ratiooftrans.dispersivitytolong.dispersivity”和字段“Ratioofvert.dispersivitytolong.dispersivity”输入横/纵和垂/纵弥散度比例值。
9.
5.5源/汇项模块选择单元选择“MT3D”菜单下的“PointSour__s/Sinks”命令打开“Well”选项输入浓度OK。
9.
5.6保存模拟选择“MT3D”菜单下的“S__eAs”命令。
9.
5.7运行MT3DMS选择“MT3D”菜单下的“RunMT3DMS”命令。
9.
5.8读入结果选择“MT3D”菜单下的“ReadSolution”命令。
9.
5.9改变等值线选项选择“Data”菜单下的“ContourOptions”命令打开“Colorfillbetweencontours”选项OK。
9.
5.10设置影象循环选择“Data”菜单下的“FilmLoop”命令单击“Setup”按钮打开“Displayclock”选项选择“Useconstantinterval”选项并将“Timeinterval”的值改为所要的值OK。第10章MT3DMS——概念模型方法
10.4导入工程选择“File”菜单下的“Open”命令打开“tutorial\modf__p\sample\sample.gpr”文件。
10.5初始化MT3DMS模拟转到3DGrid模式选择“MT3D”菜单下的“NewSimulation”命令选择“MT3D”菜单下的“BasicTransportPackage”命令。
10.
5.1指定单位单击“”按钮选择适当的单位slugfor__ssppmforcon__ntration。
10.
5.2定义物种单击“DefineSpecies”按钮输入物种名。
10.
5.3定义应力期单击“StressPeriods”按钮;输入适当的值。
10.
5.5设定输出控制单击“OutputControl”按钮;选择“Printatspecifiedtime”选项单击“Times”按钮单击“InitializeValues”按钮输入适当的值;OK。
10.
5.6选择模块单击“Packages”按钮;选择“Advectionpackage”、“Dispersionpackage”和“Sour__/SinkMixingpackage”选项OK。
10.6指定含水层参数
10.
6.1按区赋参数转到__p模块在GMS窗口中的Coverages组合框中选择“Layer1”选择SelectPolygons工具双击该层多边形选中“AquiferPorosity”选项并输入值选择“Longitudinaldispersivity”选项并输入数值OK。在GMS窗口中的Coverages组合框中选择“Layer2”选择SelectPolygons工具双击该层多边形选中“AquiferPorosity”选项并输入值选择“Longitudinaldispersivity”选项并输入数值OK。
10.7指定补给浓度值在GMS窗口中的Coverages组合框中选择“Recharge”选择SelectPolygons工具双击该层多边形在“Recharge”区的“Contaminantcon__ntration”文本框内输入浓度值。OK。
10.8转换概念模型选择“FeatureO__ects”菜单下的“__p-MT3DMS”命令保证“Allapplicablecoverages”选中OK。
10.9层的厚度MT3DMS需要HTOP数组来指定最上部含水层的顶板高程;并为各层输入厚度值。由于此处采用“truelayer”技术系统自动赋值。
10.10对流模块转到3DGrid模式选择“MT3D”菜单下的“AdvectionPackage”命令采用“MethodofCharacteristicsMOC”求解方案OK。
10.11弥散模块选择“MT3D”菜单下的“DispersionPackage”命令。纵向弥散系数自动赋值其他三个参数则需要人工赋值:Ratiooftransversedisp.tolong.disp.:
0.2Ratioofvert.disp.tolong.disp.:
0.2Effectivemoleculardiff.Coefficient:0dispersivityratios:
0.2进入第2层输入同样的值。
10.12源/汇项模块此处已经由概念模型自动赋值。
10.13保存模拟选择“MT3D”菜单下的“S__eAs”命令输入文件名。
10.14运行选择“MT3D”菜单下的“RunMT3DMS”命令。
10.15查看结果选择“MT3D”菜单下的“ReadSolution”命令OK。在窗口顶部的TS组合框内选择最后一步。用颜色填充等值线:选择“Data”菜单下的“ContourOptions”命令选择“Contourspecifiedrange”选项。
10.16设置影象循环选择“Data”菜单下的“FilmLoop”命令单击“Setup”按钮打开“Displayclock”选项选择“Useconstantinterval”选项并将“Timeinterval”的值改为所要的值OK。
10.17模拟吸附和衰变
10.
17.1选中化学反应选项ChemicalReactionsPackage选择“MT3D”菜单下的“BasicTransportPackage”命令单击“Packages”按钮选择“Chemicalreactionpackage”选项OK。
10.
17.2输入吸附和生物降解数据选择“MT3D”菜单下的“ChemicalReactionPackage”命令在“Sorption”区中选择“Linearisotherm”选项在“Kineticratereaction”区中选择“First-orderirreversiblekineticreaction”选项分别为各层输入数据层的选择通过对话框下部的框架顶部带有上下箭头的列表框OK。
10.18保存模拟选择“MT3D”菜单下的“S__eAs”命令。
10.19运行MT3DMS选择“MT3D”菜单下的“RunMT3DMS”命令。
10.20查看结果选择“MT3D”菜单下的“ReadSolution”命令OK。在窗口顶部的TS组合框内选择最后一步。
10.21绘制时间历时曲线
10.
21.1创建观测图层转到__p模块选择“FeatureO__ects”菜单下的“Coverages”命令单击“New”按钮输入图层名将图层的类型改为“Observation”单击“Options”按钮单击“New”按钮输入名称OK。
10.
21.2创建观测点选择CreatePoint工具在窗口中单击添加点并通过编辑框输入准确的坐标保持该点处于被选中状态选择“FeatureO__ect”菜单下的“Attributes”命令选中“Assigntospecifiedlayer”选项并选择layer2OK。
10.
21.3创建时间序列图然后创建时间序列图:保持该点处于被选中状态选择“Display”菜单下的“Obs.PlotOptions”命令单击“New”按钮在“Plottype”列表框中选择“Timeseries”OK。最后选择“Display”菜单下的“ShowPlotWindow”命令选择Frame命令。
10.
21.4绘制多条曲线选择“Display”菜单下的“Obs.PlotOptions”命令单击“Options”按钮选择“Useselecteddataset”选项选择要显示的物种名称选中“Display”选项。
10.
21.5__观测点选择SelectPoints/Nodes工具拖动点。第11章SEAM3DSEAM3D是一个用于模拟涉及多酶化作用和多电子接受剂的复杂生物降解的反应溶质运移模型。它基于MT3DMS程序并添加了生物降解和NAPL溶解模块。
11.1问题描述模拟NAPL包括:苯和甲苯benzeneandtoluene的污染羽的分布。
11.2启动GMS
11.3所需模块和界面
11.4导入水流模型转到3DGrid模式选择“MODFLOW”菜单下的“ReadSimulation”命令打开MODFLOW模拟文件tutorial\seam3d\flowmod.mfs选择“MODFLOW”菜单下的“ReadSolution”命令查看水流模拟结果OK。
11.5初始化SEAM3D模拟选择“MT3D”菜单下的“NewSimulation”命令。SEAM3D、RT3D和MT3DMS都用MT3D菜单。
11.6基本运移模块选择“MT3D”菜单下的“BasicTransportPackage”命令在“Model”区内选择“SEAM3D”选项。标题:输入注释说明及__姓名、日期等
11.
6.1指定单位长度和时间单位在MODFLOW中已经定义在此只需定义质量和浓度单位:单击“Units”按钮选择适当的单位。
11.
6.2设定应力期单击“StressPeriods”按钮将“Length”字段改为2000模拟时间长度;保证Transportstepsize的值为00表示MT3DMS自动计算适宜的时间步长;OK。
11.
6.3选择模块单击“Packages”按钮选中如下模块:AdvectionPackageDispersionPackageSour__/SinkMixingPackageChemicalReactionPackageBiodegradationPackageNAPLDissolutionPackageOK。
11.
6.4指定物种选择“DefineSpecies”按钮将“NondegradableTra__rs”改为2将“HydrocarbonSubstrates”改为2在“MicrobialPro__sses”区内打开“SO4reduction”选项在“Productstotrack”区内打开“H2S”选项单击“New”按钮给物种命名OK。在“Names”列表框内选择“Tra__r1”并改变其名字为“Tra__rConservative”;选择“Tra__r2”并改变其名字为“Tra__rNon-Conservative”;选择“Substrate1”并改变其名字为“Benzene”;选择“Substrate2”并改变其名字为“Toluene”;OK。
11.
6.5输出控制:单击“OutputControl”按钮选中“Printatspecifiedtimes”选项单击“Times”按钮单击“InitializeValues”按钮为“Initialtimestepsize”输入适当的值100为“__ximumtimestepsize”输入适当的值100为“__ximumsimulationtime”输入适当的值2000OKOKOK。
11.
6.6输入孔隙度单击“Porosity”按钮单击“Constant-Grid”按钮输入适当的值
0.25OK。
11.
6.7初始浓度必须为每一物种赋初始浓度其中O
24.0SO
49.0和H2S
0.01要赋背景值:在列表框中选择“O2”单击“StartingCon__ntration”按钮单击“Constant-Grid”按钮输入适当的值OK;依次输入O
24.0SO
49.0和H2S
0.01。
11.7对流模块通常要输入适当的值。但对于该例子默认值“ThirdorderTVDschemeULTI__TE”即可。选择“MT3D”菜单下的“AdvectionPackage”命令在“Solutionschem”下拉框中选择“ThirdorderTVDschemeULTI__TE”OK。
11.8弥散模块选择“MT3D”菜单下的“DispersionPackage”命令单击“LongitudinalDispersivity”按钮选择“Constant-Grid”选项输入纵向弥散系数值OK;分别为“Ratiooftrans.dispersivitytolong.dispersivity”和字段“Ratioofvert.dispersivitytolong.dispersivity”输入横/纵和垂/纵弥散度比例值。
11.9源/汇项模块通过工具选择单元选择“MT3D”菜单下的“PointSour__s/Sinks”命令打开“Constanthead”选项依次在列表框内选择O
24.0SO
49.0和H2S
0.01并输入相应的背景浓度OK。
11.10化学反应模块ChemicalReactionsPackage选择“MT3D”菜单下的“ChemicalReactionPackage”命令在“Sorption”区中选择“Linearisotherm”选项输入比重Bulkdensity5e-7;从列表框中选择“Tra__rNon-Conservative”输入“1stsorptionconstant”值
5.0e-9同样给“Benzene”和“Toluene”赋值都为
5.0e-9OK。
11.11NAPL溶解模块对污染物从污染羽中逐渐溶解到地下水中在MT3DMS中可通过定浓度单元、注水井或补给来模拟但是此法并不理想;NAPL溶解模块则比较理想。
11.
11.1选择单元单击Select__lls工具选择“Grid”菜单下的“Find__ll”命令输入单元的IJK值然后拖拉选择所要的单元污染羽中的单元。
11.
11.2指定浓度选择“MT3D”菜单下的“PointSour__s/Sinks”命令在对话框的底部选择“NAPL”选项为“Initialcon__ntration”输入适当的值2000为“Dissolutionrate”输入适当的值
3.0OK。
11.
11.3输入NAPL数据选择“MT3D”菜单下的“NAPLDissolutionPackage”命令将“Numberoftra__rsinNAPL”的值改为适当的值2将“NumberofhydrocarbonsinNAPL”的值改为适当的值2在对话框的底部为“Inertfractionmolecularweight”输入适当的值150依次选择左侧列表框内各项分别输入相应的“Initial__ssfraction”、“Solubility”和“Molecularweight”值
0.14400110OK。
11.12生物降解模块Biodegradationpackage需要输入:给水度yieldcoefficients阻碍系数inhibitioncoefficients和控制苯和甲苯降解的参数。选择“MT3D”菜单下的“BiodegradationPackage”命令
11.
12.1最小浓度选择“El.Acc.”表单为“Numberofbiostepspertransportstep”输入值2为“Microcolonyminimumcon__ntration”输入值
0.001;从文本框内选择“Benzene”并输入其Minimumcon__ntration最小浓度值
0.01从文本框内选择“Toluene”并输入其Minimumcon__ntration最小浓度值
0.05从文本框内选择“O2”并输入其Minimumcon__ntration最小浓度值
0.1从文本框内选择“SO4”并输入其Minimumcon__ntration最小浓度值
0.
211.
12.2电子接受体系数选择“Min.Conc.”表单为“SO4-O2”的“Inihibitioncoefficient”输入值
0.5选择下部文本框内每一项并输入它们的“Yieldcoefficient”值O2-Benzene:
0.5;O2-Toluene:
0.4;SO4-Benzene:
0.15;SO4-Toluene:
0.
111.
12.3生成系数generationcoefficient输入H2S的生成系数:选择“Gen.Coeffs.”表单为“Inihibitioncoefficient”输入值
0.
1511.
12.4利用系数usecoefficients输入电子接受体的利用系数:选择“UseCoeffs.”表单选择左上部文本框内每一项并输入它们的“Elec.ac__ptorusecoeff.”值O2-Benzene:
3.0;O2-Toluene:
3.2;SO4-Benzene:
8.0;SO4-Toluene:
8.
211.
12.5饱和常数saturationconstants选择“Sat.Const.”表单选择第一个文本框内每一项并输入它们的“Hydro.halfsat.const.”值O2-Benzene:
1.0;O2-Toluene:
2.0;SO4-Benzene:
1.0;SO4-Toluene:
2.0;选择第二个文本框内每一项并输入它们的“Elec.acc.halfsat.const.”值O2:
0.1;SO4:
0.
511.
12.6速率选择“Rates”表单确保对话框的“DeathRate”部分的“Calculatedbymodel”选项被选中;选择第二个文本框内每一项并输入它们的“__x.specificrateofsubstrateutilization”值O2-Benzene:
0.5;O2-Toluene:
0.05;SO4-Benzene:
0.005;SO4-Toluene:
0.
011.
12.7初始浓度利用默认值0。
11.13保存和运行选择“MT3D”菜单下的“S__eAs”命令。选择“MT3D”菜单下的“RunMT3DMS”命令
11.14读入结果选择“MT3D”菜单下的“ReadSolution”命令。
11.15设置等值线选项选择“Data”菜单下的“ContourOptions”命令打开“Colorfillbetweencontours”选项单击“ColorRampOptions”按钮打开“Showcolorlegend”选项OK
11.16查看浓度等值线首先看最后时刻2000d后conservativetra__r的解:从GMS窗口的“DataSet”组合框内选择“Tra__rConservative”数据集从“TS”中选择时刻该图则显示了没有吸附没有化学反应的浓度等值线。它可以作为比较其他污染物的基准。通过“DataSet”组合框显示其他物种。non-conservativetra__r表示有吸附但无反应。
11.17绘制时间历时曲线首先:转到__p模块选择“FeatureO__ects”菜单下的“Coverages”命令单击“New”按钮输入图层名将图层的类型改为“Observation”单击“Options”按钮单击“New”按钮输入名称OK。然后:创建观测点选择CreatePoint工具在窗口中NAPL的下游单击添加点。然后创建时间序列图:保持该点处于被选中状态选择“Display”菜单下的“Obs.PlotOptions”命令单击“New”按钮在“Plottype”列表框中选择“Timeseries”OK。最后选择“Display”菜单下的“ShowPlotWindow”命令选择Frame命令。
11.
17.1__观测点选择SelectPoints/Nodes工具拖动点。
11.
17.2绘制多条数据线选择“Display”菜单下的“Obs.PlotOptions”命令单击“Options”按钮选择“Useselecteddataset”选项选择要显示的物种名称选中“Display”选项OK;打开“Curvelegend”选项OK。
11.18其他显示选择“Data”菜单下的“FilmLoop”命令可以制作影象循环放映。第12章定义层的数据
12.1开始选择“File”菜单下的“New”命令。
12.2所需模块/界面__p模块;2DScatterPoint模块;3DGrid模块。
12.3利用“TrueLayer”模式对于多层MODFLOW模型要根据每个层的类型来赋参数。GMS可自动计算某些参数。层的参数可以通过两种方式输入:直接输入数组或“truelayer”模式。对于“truelayer”模式用户需要输入顶底板高程及水平和垂直渗透系数。
12.4通过插值利用GMS中“2DScatterPoint”模块的“Interpolation”菜单下的“toMODFLOW”命令对二维离散数据进行插值并通过“ModelChecker”检查错误。
12.5实例
12.6创建网格转到3DGrid模式选择“Grid”菜单下的“CreateGrid”命令输入如下值:X-DimensionY-DimensionZ-DimensionOrigin5005000Length25001500600#__lls50303选择ViewIAxis工具查看剖面:选择ViewKAxis工具返回平面。
12.7例1——完整层如图
12.1所示该模型包含3层且3层都扩展到整个区。
12.
7.1导入离散点数据离散点数据文件可以通过MicrosoftEx__l输入保存成以TAB分隔的文本文件扩展名为.sp2如图
12.2。其中ID列可以不要上一层的底等于下一层的顶。然后导入离散点数据文件:选择“File”菜单下的“Import”打开“tutorial\layerdata\case
1.sp2”文件。
12.
7.2对高程值进行插值在插值前需要初始化MODFLOW数据:选择“MODFLOW”菜单下的“NewSimulation”命令。转到2DScatterPoint模块选择“Interpolation”菜单下的“toMODFLOWLayers”命令将对话框左上部列表内数据集与右侧的MODFLOW数组对应单击“__p”按钮进行映射结果在下部列表中若要修改映射可单击“Un__p”按钮取消映射然后重新映射。
12.
7.3查看结果转到3DGrid模式选择某一位置的单元然后选择ViewIAxis工具查看所选单元处的剖面:选择ViewKAxis工具返回平面。
12.8例2——嵌入层如图
12.
312.
8.1导入离散点数据离散点数据文件可以通过MicrosoftEx__l输入保存成以TAB分隔的文本文件扩展名为.sp2如图
12.2。其中ID列可以不要上一层的底等于下一层的顶但是右边的数据需要调整:使第2层的底高于第2层的顶插值后再修改。然后导入离散点数据文件:导入新的数据前先删除原有的数据。转到2DScatterPoint模块选择“Edit”菜单下的“DeleteAll”;选择“File”菜单下的“Import”打开“tutorial\layerdata\case
2.sp2”文件。
12.
8.2对高程值进行插值选择“Interpolation”菜单下的“toMODFLOWLayers”命令将对话框左上部列表内数据集与右侧的MODFLOW数组对应单击“__p”按钮进行映射结果在下部列表中若要修改映射可单击“Un__p”按钮取消映射然后重新映射。
12.
8.3修正层的数据转到3DGrid模式选择“MODFLOW”菜单下的“CheckSimulation”命令单击“RunCheck”按钮左上侧列表框内显示错误信息单击“FixLayerErrors”按钮;在右侧列表框内选择“Layer2”选择“__erage”选项单击“FixSelectedLayer”按钮OK。
12.9例3——含水层出露如图
12.
412.
9.1导入离散点数据离散点数据文件可以通过MicrosoftEx__l输入保存成以TAB分隔的文本文件扩展名为.sp2如图
12.2。其中ID列可以不要上一层的底等于下一层的顶但是右边的数据需要调整:使第1层的底高于第2层的顶插值后再修改。然后导入离散点数据文件:导入新的数据前先删除原有的数据。转到2DScatterPoint模块选择“Edit”菜单下的“DeleteAll”;选择“File”菜单下的“Import”打开“tutorial\layerdata\case
3.sp2”文件。
12.
9.2对高程值进行插值选择“Interpolation”菜单下的“toMODFLOWLayers”命令将对话框左上部列表内数据集与右侧的MODFLOW数组对应单击“__p”按钮进行映射结果在下部列表中若要修改映射可单击“Un__p”按钮取消映射然后重新映射。
12.
9.3修正层的数据转到3DGrid模式选择“MODFLOW”菜单下的“CheckSimulation”命令单击“RunCheck”按钮左上侧列表框内显示错误信息单击“FixLayerErrors”按钮;在右侧列表框内选择“Layer1”选择“Preservetop”选项单击“FixSelectedLayer”按钮OK。
12.
9.4显示结果选择ViewIAxis工具查看剖面:选择ViewKAxis工具返回平面。
12.10例4——基岩切割如图
12.
512.
10.1导入离散点数据离散点数据文件可以通过MicrosoftEx__l输入保存成以TAB分隔的文本文件扩展名为.sp2如图
12.2。其中ID列可以不要上一层的底等于下一层的顶。然后导入离散点数据文件:导入新的数据前先删除原有的数据。转到2DScatterPoint模块选择“Edit”菜单下的“DeleteAll”;选择“File”菜单下的“Import”打开“tutorial\layerdata\case
3.sp2”文件。
12.
10.2对高程值进行插值选择“Interpolation”菜单下的“toMODFLOWLayers”命令将对话框左上部列表内数据集与右侧的MODFLOW数组对应单击“__p”按钮进行映射结果在下部列表中若要修改映射可单击“Un__p”按钮取消映射然后重新映射。
12.
10.3修正层的数据转到3DGrid模式选择“MODFLOW”菜单下的“CheckSimulation”命令单击“RunCheck”按钮左上侧列表框内显示错误信息单击“FixLayerErrors”按钮;在右侧列表框内选择“Layer3”选择“Truncatetobedrock”选项单击“FixSelectedLayer”按钮OK。
12.
10.4显示结果选择ViewIAxis工具查看剖面:选择ViewKAxis工具返回平面。选择“View”菜单下的“GeneralMode”命令选择ObliqueView工具选择Shade工具选择Rotate工具。第13章非稳定流数据管理对于非稳定流需要大量的数据GMS提供了适当的工具对它们进行输入和管理。
13.1问题描述
13.2开始
13.3所需模块3DGrid模块__p模块MODFLOW界面
13.4读入工程选择“File”菜单下的“Open”命令打开“tutorial\trans\start.gpr”工程。
13.5非稳定流数据通过概念模型获得。
13.6在__p模式下输入非稳定流数据
13.
6.1指定非稳定补给速率转到__p模块选择SelectPolygon工具选择一个区区2选择“FeatureO__ects”菜单下的“Attributes”命令选择“RechargeRate”下面的“Transient”选项单击“Transient”下面的区域启动XYSeriesEditor.编辑器单击“New”按钮创建新的时间序列输入新的名字zone2并输入新的时间数据序列;也可以通过“Import”按钮将已经准备好的文本文件导入。
13.
6.2导入开采量数据开采井的数据文件包括两个文件:第一个文件包括井名、滤管几何形状、井的XY坐标。第二个文件包括开采时间安排其文件各格式为:导入开采井数据:从GMS窗口的coverage组合框中选择“Sour__sSinks”图层选择“File”菜单下的“Import”命令选择文件pumping.wpfOpenYes。双击井单元查看数据。
13.
6.3指定重力给水度在__p模块下选择HydraulicConductivity图层选择SelectPolygon工具选择一个区区2和1选择“FeatureO__ects”菜单下的“Attributes”命令选中“SpecificYield”并输入适当的值
0.2;选择一个区区3和4选择“FeatureO__ects”菜单下的“Attributes”命令选中“SpecificYield”并输入适当的值
0.15。
13.7初始化MODFLOW
13.
7.1将MODFLOW设置为非稳定流转到3DGrid模块选择“MODFLOW”菜单下的“BCFPackage”命令在“ModelType”中选择“Transient”。
13.
7.2指定模拟参考时间选择“MODFLOW”菜单下的“BCFPackage”命令单击“StressPeriods”按钮在“Referen__Time”区下选择“Edit”按钮输入初始时间年月日时分秒1985_____00000。注意其下的时间显示选项保证“Date/time”。
13.
7.3设置应力期单击“Initialize”按钮输入No.ofstressperiods7、length100和No.oftimesteps1;然后在stressperiods对话框中修改应力期及时间步数:通过对话框中左侧的带上下箭头的列表框改变应力期数在右侧的No.oftimesteps文本框内输入步数;依次输入各应力期的结束时间和期间的步数。
13.8转换概念模型转到__p模块选择“FeatureO__ects”菜单下的“__p-MODFLOW”命令OK。
13.9保存和运行MODFLOW转到3DGrid模块选择“MODFLOW”菜单下的“S__eAs”和“RunMODFLOW”命令。
13.10读入结果并显示选择“MODFLOW”菜单下的“ReadSolution”命令。
13.11创建影象循环选择“Data”菜单下的“FilmLoop”命令单击“Setup”按钮。
13.12非稳定观测数据
13.
12.1导入非稳定观测数据选择“File”菜单下的“Open”命令打开扩展名为“*.__p”的文件obswells.__p。其格式为:
13.
12.2导入文件:选择“File”菜单下的“Import”命令打开文件trans_obs.tof。出现ObservationCoverageOptions对话框在“MeasurementType”列表中选择“trans_head”保证该种测量模式被用于残生校正目标。
13.
12.3创建非稳定流观测曲线转到__p模块选择“Display”菜单下的“Obs.PlotOptions”命令单击“New”按钮通过下拉列表框将图形类型设置为ErrorvstimestepOK。选择“Display”菜单下的“ShowPlotWindow”命令选择“FrameI__ge”命令。Me为平均误差;__e为平均绝对误差;rms为误差平方的平均值的平方根。此图显示了测量值与模拟值的误差。时间序列图:转到__p模块选择两个观测井选择“Display”菜单下的“Obs.PlotOptions”命令单击“New”按钮通过下拉列表框将图形类型设置为Timeseries选中“Curvelegend”选项OK。选择“Display”菜单下的“ShowPlotWindow”命令选择“FrameI__ge”命令。第14章由区域模型转为局域模型RegionaltoLocalModelConversion对于许多模拟研究确定合适的边界条件是非常困难的。可通过先粗后细的方法分两步进行模拟。问题描述:如图
14.1所示:大部分边界为隔水边界河流的出入口为已知水头边界;两层;局域模型位于其中心位置。由区域模型转为局域模型的基本目标是从区域模型中创建一个包含水头和层参数的二维离散点集建立局域模型并将水头和层参数插值到局域模型。基本步骤为:建立区域模型并进行求解利用“MODFLOWLayers-2DScatterPoints”命令从区域模型中创建二维离散数据集为局域模型创建三维网格并为其进行插值。
14.1起始
14.2所需模块/界面2DGrid模块3DGrid模块2DScatterPoint模块__p模块;MODFLOW界面。
14.3读入区域模型选择“File”菜单下的“Open”命令打开“tutorial\reg2loc\regmod.gpr”文件。该模型利用概念模型创建。
14.4将层参数转换成离散点数据进入3DGrid模块选择“Grid”菜单下的“MODFLOWLayers-2DScatterPoints”命令输入文件名RegionalDataOK。
14.5创建局域模型的方法方法有很多一种通用的方法是将边界作为已知水头边界水头值采用区域模型的计算值这可以通过两种方法来实现:1构建一个简单的矩形整个网格的周边作为已知水头边界;2采用一个四边形其中的两个相对的边平行于区域模型的等水头线即各边上的水头相等另外两个边垂直于区域模型的等水头线为隔水边。我们采用方法2。
14.6创建局域概念模型简单的办法是在__p模块中建立概念模型要创建一个新的sour__/sink图层。首先删除区域模型中的sour__/sink图层。转到__p模块选择“FeatureO__ect”菜单下的“Coverages”命令激活要删除的图层ss单击“Delete”按钮OK。
14.
6.1创建新的图层单击“New”按钮输入新的图层名localss将图层类型设置为“MODF/MT3D/MODP”OK。
14.
6.2创建边界弧段选择Zoom工具拖拉鼠标放大局域模型部位。按如__法创建边界弧段:选择CreateArc工具创建4条弧段两条平行于等水头线图
14.2在终点出双击鼠标
14.
6.3创建多边形选择“FeatureO__ect”菜单下的“BuildPolygons”命令OK。要确保多边形被指定到所有层。选择SelectPolygons工具双击多边形将层指定为“fromlayer1to2”OK。
14.
6.4标记指定水头弧单击SelectArcs工具选择弧段东西两条选择“FeatureO__ect”菜单下的“Attributes”命令改为“SpecifiedHead”类型将层指定为“fromlayer1to2”OK。
14.7创建局域MODFLOW模型选择“FeatureO__ects”菜单下的“GridFrame”命令单击“DeleteFrame”删除老的网格边框单击“NewFrame”按钮创建新的网格边框新边框应比弧段稍大如果需要可以通过旋转工具使网格边框更靠近局域模型边界。
14.
7.1创建网格选择“FeatureO__ects”菜单下的“__p-3DGrid”命令OKOK。在对话框内输入参数。
14.
7.2激活单元选择“FeatureO__ects”菜单下的“Activate__llsinCoverage”命令OKOK。在对话框内输入参数。
14.
7.3映射属性将MODFLOW数据转入网格:进入3DGrid模块选择“MODFLOW”菜单下的“NewSimulation”命令转到__p模块选择“FeatureO__ect”菜单下的“__p-MODFLOW”命令OK。
14.8层参数插值转到2DScatterPoint模块选择“Interpolation”菜单下的“toMODFLOWLayers”命令在对话框的底部选择“Applystartingheadtoall__lls”选项OK。
14.9保存和运行局域模型进入3DGrid模块选择“MODFLOW”菜单下的“S__eAs”命令输入文件名locmod.mfs;选择“MODFLOW”菜单下的“RunMODFLOW”命令OK。
14.10查看结果选择“MODFLOW”菜单下的“ReadSolution”命令OK。第15章模型校准地下水模拟的一个重要部分是模型校准。
15.1问题描述如图
15.1所示降水量
1.25ft/yr__28平方英里地下水通过排水渠排入高程为1000ft的湖中。盆地周围为弱透水的基岩。补给和水力参数分区如图
15.2。
15.2开始
15.3所需模块3DGrid模块__p模块;MODFLOW界面。
15.4读入模型选择“Open”菜单下的“File”命令打开文件“tutorial\calib\bigval.gpr”。
15.5观测数据校准中可以用两种观测数据:观测井的水位和河流观测流量。由于位于较干旱区假设河流中的绝大部分水流来自地下水。
15.6输入观测点首先要输入观测点就要建立观测图层。
15.
6.1创建观测图层输入观测点之前要先创建观测图层:转到__p模块选择“FeatureO__ects”菜单下的“Coverages”命令单击“New”按钮输入名称ObservationWells在“thecoveragetype”下拉列表中选择“Observation”然后单击其下面的“Options”按钮;在弹出的对话框中选择“3DGrid”和“2Dwithinlayer”选项在MODFLOW中这是最佳选择它严格地把插值和观测点放在同一层在“MeasurementType”区内单击“New”按钮输入观测类型名HeadOKOK。
15.
6.2创建观测点第一个点按如下值创建:其中interval表示观测值的估计误差confiden__为误差估计的可靠度interval可被用作校准目标。当误差在观测值的估计误差范围此处为
1.5ft之内时则达到校准的目的也就是说计算水头位于
997.5-
1000.5之间时此观测井的校准目的已经达到。创建点:选择CreatePoint工具在模型中的任意处单击双击GMS顶部的X字段编辑框输入准确的坐标值1466132694指定该点的属性值:选择“FeatureO__ects”菜单下的“Attributes”命令输入点的名字Point#1打开“Observed”选项输入观测值Observed
999.0输入误差范围Interval
1.5输入可靠度Confiden__95OK。
15.
6.3校准目标校准目标显示在点的旁边组分如图
15.3所示:如果误差柱整个在目标范围之内则显示为绿色;如果误差柱整个在目标范围之外但小于200%显示为黄色;若大于200%显示为红色。
15.
6.4点统计单击SelectPoints/Nodes按钮选择观测点在GMS底部的帮助条倒数第二个上显示统计信息。
15.7读入观测点数据集用上述方法输入其他的观测点。为节省时间直接读入已经存在的图形文件。
15.
7.1删除当前图层读入观测点之前首先删除当前图层:选择“FeatureO__ects”菜单下的“Coverages”命令选择“ObservationWells”图层单击“Delete”按钮OKOK。
15.
7.2读入点集选择“File”菜单下的“Open”命令打开文件obswells.__pOpen。
15.8输入观测河流流量观测河流流量直接指定到概念模型的localsour__/sinkcoverage中的弧段和多边形。GMS提供了两种方法进行赋值:单个弧段或一组弧段。由于区内河流流量是在出口观测点观测流量代表了整个区的流量因此需要创建弧段组并将观测流量赋值到弧段组。GMS自动将所有河流弧段的计算流量求和并与观测流量对比。创建弧段组:在GMS顶部的Coverage组合框内选择Sour__sSinks图层单击SelectArcs工具按下SHIFT键选择河流弧段选择“FeatureO__ect”菜单下的“CreateArcGroup”命令创建一个新的对象我们可以为弧段组指定观测流量单击SelectArcGroups工具双击河流弧段组中的任一条弧段选中“Observedflux”选项单击“Options”按钮在观测值区域为“Steadystatefluxrate”输入值–164000在“Confiden__interval”区域为“Interval”输入值7500。
15.9生成误差图在GMS顶部的Coverage组合框内选择ObservationWells图层选择“Display”菜单下的“Obs.PlotOptions”命令单击“New”按钮选择“Computedvs.observed”类型再次单击“New”按钮选择“ErrorSum__ry”类型OK。选择“Display”菜单下的“ShowPlotWindow”命令单击“FrameI__ge”按钮。
15.
9.1计算—观测图在PlotWindow窗口内单击Computedvs.Observed标志点在GMS底部的帮助条倒数第二个上显示统计信息。
15.
9.2误差统计信息
15.10编辑水力传导系数在GMS顶部的Coverage组合框内选择HydraulicConductivity图层修改渗透系数:单击SelectPolygons工具按下SHIFT键选择区域12选择“FeatureO__ects”菜单下的“Attributes”命令输入HorizontalK
2.0OK;双击区域3输入HorizontalK
2.0OK。
15.11转换模型修改数据后将模型转成网格数值模型:选择“FeatureO__ects”菜单下的“__p-MODFLOW”命令确定“Allapplicablecoverages”选项被选中OK。
15.12运行模型
15.
12.1保存模拟转到3DGrid模式选择“MODFLOW”菜单下的“S__eAs”命令。
15.
12.2运行选择“MODFLOW”菜单下的“RunMODFLOW”命令。
15.13读入结果选择“MODFLOW”菜单下的“ReadSolution”命令。
15.14误差——模拟图在校准过程中要注意误差的变化趋势可通过该图来显示。转到__p模块选择“Display”菜单下的“Obs.PlotOptions”命令选择Plot#2单击“Delete”按钮;单击“New”按钮为Plottype选择“Errorvs.Simulation.”类型再次单击“New”按钮为Observationtype选择“Flux”类型为Plottype选择“Errorvs.Simulation.”类型OK。选择“Display”菜单下的“ShowPlotWindow”命令单击“FrameI__ge”按钮。
15.15继续校准
15.
15.1改变值与改变区除了改变渗透系数和补给外还可以改变区的大小。
15.
15.2查看响应要查看响应需导入一个包含符合校准目标的参数数值的__p文件。读入新的图层之前首先删除其他3个图层:转到__p模块选择“FeatureO__ects”菜单下的“Coverages”命令删除图层。选择“File”菜单下的“Open”命令打开文件“answer.__p”。然后再运行。第16章自动参数识别
16.1问题概述同上一章。
16.2启动GMS
16.3所需模块3DGrid模块__p模块;MODFLOW界面PEST/UCODE界面。
16.4读入模型选择“Open”菜单下的“File”命令打开文件“tutorial\inberse\bigval.gpr”。
16.5定义参数区建立反演模型时首先要定义参数区可用概念模型为每一个参数区指定关键值一般为负值。
16.
5.1设立渗透系数区如图
16.15个区定义成4个参数区输入关键值图上:在GMS顶部的Coverage组合框内选择HydraulicConductivity图层单击SelectPolygons工具双击各区并为HorizontalK赋图上的值。
16.
5.2设立补给区如图
16.25个区定义成4个参数区顶端的区很小它隔离于大部分观测井而且位于井的下游所以它不好用来进行参数估计我们将该区的补给设置为0为其他3个区的指定关键值图上:在GMS顶部的Coverage组合框内选择Recharge图层单击SelectPolygons工具双击各区并为Rechargerate赋图上的值。
16.
5.3将关键值赋给网格单元选择“FeatureO__ects”菜单下的“__p-MODFLOW”命令OK。
16.6打开PEST选项转到3DGrid模块选择“MODFLOW”菜单下的“BasicPackage”命令单击“Packages”按钮在ParameterEsti__tion区下拉列表内选择PEST选项OK。
16.7编辑参数然后我们创建一个包含初始值、最大值和最小值的参数列表。选择“MODFLOW”菜单下的“PEST/UCODE/Options”命令:Parameters区用来管理用于反演模型参数列表;NewParameter按钮可用来创建参数集;每个参数包含一个名称、关键值、初始值、最小值和最大值。然而在大多数情况下GMS可以用InitializefromModel按钮自动产生参数列表。当选择该按钮GMS搜索MODFLOW输入数据中的负值并创建参数列表。如果参数没有被GMS自动发现则需要利用NewParameter按钮来添加。单击InitializefromModel按钮。8个参数自动被找到并有一个默认名下一步则需要输入初始值、最大值和最小值。其中初始值越靠近真实值越好否则可能不收敛。
16.8编辑PEST选项在名为PESTOptions的选项卡propertysheet保留之中的缺省值。
16.9收敛选项PEST在迭代求解过程中不断调用MODFLOW因此需要保证MODFLOW收敛:选择“MODFLOW”菜单下的“SolverPackages/Convergen__Options”命令打开“Prevent__lldrying”选项为Minimumsaturatedthickness输入值
0.1OK。
16.10初始水头运行PEST之前要将计算水头拷贝到初始水头数组中这将保证每次PEST调用MODFLOW时初始水头尽可能靠近最终水头值MODFLOW收敛速度加快。选择“MODFLOW”菜单下的“BasicPackage”命令单击“StartingHead”按钮选择“3DDataSet-Grid”命令OK;单击OK确认水头只用于变水头单元选择“run1”数据集OK;Close。
16.11最大迭代选择“MODFLOW”菜单下的“SolverPackages/PCG2Package”命令将__ximumouteriterations改为100。
16.12保存模拟转到3DGrid模式选择“MODFLOW”菜单下的“S__eAs”命令pest
1.mfs。
16.13保存工程转到__p模块选择“File”菜单下的“S__eAs”命令输入文件名bigval
2.gpr。
16.14运行GMS2PEST运行PEST之前首先运行GMS2PEST。GMS2PEST为一个分析GMS输出文件和并为PEST准备最终数据的单元。选择“MODFLOW”菜单下的“PEST/UCODED/RunGMS2PEST”命令OK。
16.15运行PEST选择“MODFLOW”菜单下的“PEST/UCODED/RunPEST”命令OK。φ值为最小化的目标函数值。
16.16查看结果PEST的结果为一系列最优化的参数值。方法1:查看包含PEST输出结果的文本文件然后将最优参数导入GMS并装入到MODFLOW的正确位置最后保存包含最优参数的模拟运行MODFLOW读入结果。
16.
16.1查看PEST输出文件选择“File”菜单下的“EditFile”命令打开文件pest
1.recOK。文件的第一部分为输入总结;第二部分为迭代总结;下面是最优参数值的总结和一系列残差值。
16.
16.2读入最优参数值选择“MODFLOW”菜单下的“PEST/UCODE/Options”命令选择“ImportParameterFile”按钮选择文件pest
1.parOK。
16.
16.3将最优参数值拷贝到MODFLOW模型单击“LoadtoModel”按钮选择“LoadtoMODFLOWparametersonly”选项OK。
16.
16.4关闭PEST选项在保存模拟和运行MODFLOW之前必须关掉PEST选项。转到3DGrid模块选择“MODFLOW”菜单下的“BasicPackage”命令单击“Packages”按钮在ParameterEsti__tion区下拉列表内选择none选项OK。
16.
16.5保存和运行MODFLOW转到3DGrid模式选择“MODFLOW”菜单下的“S__eAs”命令输入文件名solution
1.mfs。选择“MODFLOW”菜单下的“RunMODFLOW”命令。选择“MODFLOW”菜单下的“ReadSolution”命令。转到__p模块选择“Display”菜单下的“Obs.PlotOptions”命令查看结果。
16.17用UCODE进行反演
16.
17.1映射关键值选择“FeatureO__ects”菜单下的“__p-MODFLOW”命令
16.
17.2打开UCODE选项转到3DGrid模块选择“MODFLOW”菜单下的“BasicPackage”命令单击“Packages”按钮在ParameterEsti__tion区下拉列表内选择UCODE选项OK。
16.
17.3编辑关键值修改关键值由于模型中已经是PEST的最优值需要修改。选择“MODFLOW”菜单下的“PEST/UCODE/Options”命令:在选项卡中输入原来的参数初始值同PEST。
16.
17.4编辑UCODE选项在UCODEOptions选项卡中编辑此处仍保留默认值。
16.
17.5保存模拟转到3DGrid模式选择“MODFLOW”菜单下的“S__eAs”命令ucode
1.mfs。
16.
17.6运行GMS2UCODE运行PEST之前首先运行GMS2UCODE。GMS2UCODE为一个分析GMS输出文件和并为PEST准备最终数据的单元。选择“MODFLOW”菜单下的“PEST/UCODED/RunGMS2UCODE”命令OK。
16.
17.7运行UCODE选择“MODFLOW”菜单下的“PEST/UCODED/RunUCODE”命令OK。
16.
17.8查看结果选择“File”菜单下的“EditFile”命令打开文件ucode
1._oOK。读入最优参数值选择“MODFLOW”菜单下的“PEST/UCODE/Options”命令选择“ImportParameterFile”按钮选择文件ucode
1.parOK。将最优参数值拷贝到MODFLOW模型单击“LoadtoModel”按钮选择“LoadtoMODFLOWparametersonly”选项OK。保存和运行MODFLOW转到3DGrid模式选择“MODFLOW”菜单下的“S__eAs”命令输入文件名solution
2.mfs。选择“MODFLOW”菜单下的“RunMODFLOW”命令。选择“MODFLOW”菜单下的“ReadSolution”命令。转到__p模块选择“Display”菜单下的“Obs.PlotOptions”命令查看结果。第17章FEMWATERFEMWATER是模拟饱和—非饱和三维水流和溶质运移的有限元模型。
17.1问题概述如图
17.1所示。三口开采井开采量为_____0ft^3/day上部含水层的K为10ft/day下部的为30ft/day。
17.3所需模块Mesh模块2DMesh和3DMesh__p模块Subsu_____Characterization模块Borehole、TIN和SolidGeostatistics模块2DScatterPoint和3DScatterPointFEMWATER界面。
17.4建立概念模型FEMWATER可以用概念模型的方法进行建模。
17.
4.1导入背景图选择“File”菜单下的“Open”命令打开“tutorial\femwater\start.gpr”文件。
17.
4.2初始化FENWATER图层转到__p模块选择“FeatureO__ects”菜单下的“Coverages”命令将Coverages对话框内的默认图层的名字设置为“femwater”将图层类型Coveragetype设置为“FEMWATER”。
17.
4.3指定单位在Coverages对话框内单击“Units”按钮社设置Length和Time的单位采用默认值质量__ss单位为slug。
17.
4.4创建边界弧段单击“CreateArc”工具按钮
17.
4.5修正弧段的顶点弧段的两端端点称为结点nodes中间点称为顶点verti__s。这些弧段用于生成可以转换成3D网格的2D网格。由顶点限定的线段间距控制了单元的大小和数目。因此我们需要沿每一条弧段对定顶点进行调整以保证间隔比较均匀而且长度合适。单击“SelectArcs”按钮拖拉鼠标选择所有弧段选择“FeatureO__ects”菜单下的“RedistributeVerti__s”命令将“Targetspacing”改为300。
17.
4.6指定边界条件用弧段指定边界条件。单击“SelectArcs”按钮按下SHIFT键用鼠标单击选择定水头弧段选择“FeatureO__ects”菜单下的“Attributes”命令打开“Head/fluidflux”选项OK。此时并没有赋水头值水头值赋在两端结点上再通过线性插值获得。沿海弧段的水头值为0可直接用默认值河流下游点也为0不用再赋。选择SelectPoints/Nodes工具单击河流弧段的左上结点输入水头值200。
17.
4.7创建多边形需要多边形的原因:1当产生网格时定义模型范围;2用于赋补给值。选择“FeatureO__ects”菜单下的“BuildPolygon”命令OK。
17.
4.8指定补给在FEMWATER中由两种方法:用指定流量边界或变化的边界条件。后者更精确但更耗时且不稳定。此处采用简单的指定流量法。选择“SelectPolygons”按钮在区内双击打开“Fluidflux”选项输入
0.003OK。
17.
4.9创建井单击CreatePoint工具按钮在模型的右上部单击在GMS顶部的坐标框内输入准确的坐标值5290420645选择“FeatureO__ects”菜单下的“Attributes”命令打开“Refinemesharoundpoint”选项井附近的单元细化为“Elementsize”输入值150控制井处的单元大小。然后将点指定为井点指定开采量、滤管的顶底高程流量赋予期间的单元。打开“Well”选项输入Topofscreen–185和Bottomofscreen–195及Fluxrate–_____0OK。同上输入其他两个井:
17.5创建三维网格此时概念模型已经完成要准备创建三维有限元网格。网格有两个区组成一个区为上部含水层另一个区为下部含水层。上区由两层单元组成;下区由三层单元组成。要创建网格我们首先利用概念模型的特征对象创建二维投影网格。然后创建三个不规则三角形网格TINs:一个为顶地面面一个为上部含水层的底面另一个为下部含水层的底面。然后通过延伸二维网格的单元形成三维单元的单元柱从而形成三维网格单元。
17.
5.1指定材料创建网格之前我们要为含水层指定材料。岩性指定给TINs最后到3D网格。选择“Edit”菜单下的“__terials”命令改变缺省材料名为UpperAquifer在Display区内颜色窗口内单击将材料改为绿色。单击“New”按钮创建另一材料改名为LowerAquifer改为红色。
17.
5.2创建2D投影网格选择“FeatureO__ects”菜单下的“__p-2DMesh”命令保留缺省值。
17.
5.3创建TINs要创建定义地层的TINs我们要创建3个TINs每个TINs为二维网格的拷贝。首先这三个TINs具有与二维网格相同的高程0。我们然后导入离散点数据集并插值到TINs。创建顶部TINs:转到2DMesh模块选择“BuildMesh”菜单下的“Mesh-TIN”命令输入TINs名terrain;从材料__terials下拉列表中选择“UpperAquifer”OK;在提示框内单击No按钮。创建第二个TIN:选择“BuildMesh”菜单下的“Mesh-TIN”命令输入TINs名bottomupperaquifer;从材料__terials下拉列表中选择“LowerAquifer”OK;在提示框内单击No按钮不删除二维网格。创建第三个TIN:选择“BuildMesh”菜单下的“Mesh-TIN”命令输入TINs名bottomloweraquifer;从材料__terials下拉列表中选择“LowerAquifer”OK;在提示框内单击No按钮不删除二维网格。
17.
5.4导入和插值地面数据然后我们导入指定地面高程的离散点和将其插值到顶部TIN。地形点通过数字化获得并存成表形式的离散点文件。选择“File”菜单下的“Import”命令打开文件terrian.sp2。在插值前要保证顶部TIN是激活的。转到TIN模式在屏幕的左上角的TIN组合框内选择terrain。插值前我们将对插值选项做一些调整:转到2DScatterPoint模块选择“Interpolation”菜单下的“Interp.Options”命令单击“Inversedistan__weighted”旁边的“Options”按钮在“Nodalfunction”区选择“Constant”选项OK。将离散点数据插值到TIN:选择“Interpolation”菜单下的“toActiveTIN”命令打开“__pelevations”选项OK。查看插值高程:单击ObliqueView命令按钮选择“View”菜单下的“Z__gnification”命令输入放大倍数
4.0OK。
17.
5.5导入和插值地面数据然后我们导入和插值顶部和底部含水层的高程。数据通过钻孔资料获得并存成表形式的离散点文件。选择“File”菜单下的“Import”命令打开文件elevs.sp2。在插值前要保证所需TIN是激活的。转到TIN模式在屏幕的左上角的TIN组合框内选择bottomupperaquifer。插值前我们将对插值选项做一些调整:转到2DScatterPoint模块选择“Interpolation”菜单下的“Interp.Options”命令单击“Inversedistan__weighted”旁边的“Options”按钮在“Nodalfunction”区选择“Gra___ntplane”选项OK。离散点集有两个数据集:一个为顶部含水层的底板高程另一个为底板高程。我们将首先为底部含水层高程插值:在GMS顶部的DataSet组合框的botoflayer1数据。将离散点数据插值到TIN:选择“Interpolation”菜单下的“toActiveTIN”命令OK。转到TIN模式在屏幕的左上角的TIN组合框内选择bottomloweraquifer转到2DScatterPoint模块在GMS顶部的DataSet组合框的botoflayer2数据选择“Interpolation”菜单下的“toActiveTIN”命令OK。查看插值高程:转到TIN模式选择“Display”菜单下的“DisplayOptions”命令关掉“Verti__s”选项OK。转到FrameI__ge命令选择Shade命令。单击ObliqueView命令按钮选择“View”菜单下的“Z__gnification”命令输入放大倍数
4.0OK。
17.
5.6创建三维网格三维网格通过选择一对TINs并填充挤压单元形成三维单元。选择SelectTINs工具按下SHIFT键选择顶部和中间的TINs选择“BuildTIN”菜单下的“FillBetweenTINs-3DMesh”命令为“Numberofinterpolatedmeshlayers”输入2选择“Specify__terial”选项选择材料UpperAquiferOK。然后创建下部含水层的单元:按下SHIFT键选择底部和中间的TINs选择“BuildTIN”菜单下的“FillBetweenTINs-3DMesh”命令为“Numberofinterpolatedmeshlayers”输入2选择“Specify__terial”选项选择材料LowerAquiferOK。
17.6隐藏对象隐藏TINs:拖拉鼠标选中TINs单击Hide命令。隐藏离散点:转到2DScatterPoint模式单击SelectScatterPointSets工具拖拉鼠标选中TINs单击Hide命令。隐藏2D网格:转到2DMesh模式选择“Display”菜单下的“DisplayOptions”命令关掉“Nodes”和“Elements”选项OK。整理三维网格的显示:转到模块选择“Display”菜单下的“DisplayOptions”命令关掉“Nodes”选项OK。选择Shade命令。
17.7转换概念模型选择“FEMWATER”菜单下的“NewSimulation”命令转到__p模式选择“FeatureO__ects”菜单下的“__p-FEMWATER”命令在提示框内单击“OK”按钮;
17.8选择分析选项
17.
8.1输入运行选项转到3DMesh模块选择“FEMWATER”菜单下的“RunOptions”命令为“Typeofsimulation”选择“Flowonly”在“SteadyStatevs.Transient”区内选择“Steadystatesolution”。我们研究的问题大多是部分饱和区主要在模型的左上部。非饱和区越大FEMWATER越难收敛。对于这种问题对求解quadrature选择“Nodal/Nodal”较好虽然不及默认选项“Gaussian/Gaussian”精确但更稳定。在“Quadratureselection”区内选择“Nodal/Nodal”选项。
17.
8.2设置迭代参数选择“FEMWATER”菜单下的“IterationParameters”命令设置“__xiterationsfornon-linearequation”为100设置“__xiterationsforlinearequation”为1000设置“Steadystateconvergen__criterion”为
0.01OK。
17.
8.3选择输出控制选择“FEMWATER”菜单下的“OutputControl”命令关掉“S__enodalmoisturecontentfile”和“S__evelocityfile”选项OK。
17.
8.4指定流体属性最后指定流体属性。选择“FEMWATER”菜单下的“FluidProperties”命令。对话框内的各项依赖于我们指定的单位。由于是稳定流解忽略流体的粘滞性和可压缩性。为“Densityofwater”输入值
1.94为“Ac__lerationofgr__ity”输入值
2.4037e11OK。
17.9指定初始条件由于FEMWATER模拟非饱和带的水流FEMWATER对初始条件比MODFLOW更为敏感。如果初始条件与最后的水头分布差异较大则FEMWATER收敛较慢或不能收敛。对于水流模拟FEMWATER需要压力水头分布作为初始条件。GMS中的FEMWATER界面中包含能够通过用户指定的水面来自动产生压力水头数据。地下水面是通过离散点获得的。总的水头数据通过插值赋予三维网格单元。最后压力水头数据通过从总数据中减去节点高程获得。
17.
9.1创建离散点定义初始条件我们将创建一些需要得到水面高程的点。一种方法是创建包含高程的表式离散点文件。另一种快速方法为互动式创建点。用后者。选择“FEMWATER”菜单下的“BCDisplayOptions”命令关掉“Flux”选项OK。在屏幕上数字化点我们将创建新的TIN和创建TIN顶点。当我们创建每一个点GMS提示顶点高程。一旦创建了点后将顶点转成二维离散点数据集。转到PlanView模式转到TIN模块选择“BuildTIN”菜单下的“NewTIN”命令输入名称NamestartingheadOK。然后创建点顶点按图
17.2创建点。单击CreateVertex按钮在模型内单击鼠标输入高程值。一旦创建了TIN顶点转成离散数据点:选择“BuildTIN”菜单下的“TIN-ScatterPoints”命令输入数据名StartingHeadsOK单击“Yes”删除TIN。
17.
9.2创建数据集最后要定义初始条件我们必须创建压力水头数据集并存成文件目录同FEMWATER的模拟文件目录。转到3DMesh模块选择“FEMWATER”菜单下的“InitialConditions”命令在对话框左上角的中部的PressureHead区域中选择“Readfromdatasetfile”选项单击“GenerateIC”按钮为“Active2Dscatterpointset”选择“Starting_Heads”为“Mininumpressurehead”输入值–200OK。输入文件名starthd.phdOK。
17.10指定材料属性建模的最后一步是指定材料的属性。选择“FEMWATER”菜单下的“__terialProperties”命令。我们需要输入渗透系数和每一层的非饱和区曲线。对话框左边的项为溶质运移模拟使用。指定上部含水层的材料属性:在对话框的顶部列表框中选择“Upperaquifer”项在“Conductivity”框架内的左上部编辑框内输入
10.0单击“Isotropic”按钮给其他方向赋同样的渗透系数。单击“GenerateCurves”按钮弹出vanGenuchtenCurveGenerator对话框为“__x.heightabovewatertable”输入值
6.0在“Presetparametervalues”区域选择“Silt”单击“ComputeCurves”按钮OK。指定上部含水层的材料属性:在对话框的顶部列表框中选择“Loweraquifer”项在“Conductivity”框架内的左上部编辑框内输入
30.0单击“Isotropic”按钮给其他方向赋同样的渗透系数。单击“GenerateCurves”按钮弹出vanGenuchtenCurveGenerator对话框单击“ComputeCurves”按钮OK。
17.11保存和运行模型选择“File”菜单下的“S__eAs”命令输入工程名称femmod.gprS__e。注意:通过保存工程我们保存了所有的数据。也保存了FEMWATER模拟的一系列文件femmod.*也可以通过“FEMWATER”菜单下的“S__e/S__eAs”命令来保存FEMWATER模拟而不保存其他数据。运行FEMWATER:选择“FEMWATER”菜单下的“RunFEMWATER”命令OK。
17.12查看结果选择“FEMWATER”菜单下的“ReadSolution”命令OK。
17.
12.1查看条纹图选择“Display”菜单下的“DisplayOptions”命令打开“Fringes”选项选择“FEMWATER”选项卡关掉“Head”和“Flux”选项OK。转到oblique视图选择Shade命令。
17.
12.2查看水面等值面一个更好的显示解的方法是在压力水头为0时的生成等值面。这创建一个符合计算水位的面。而且如果将等值面盖在大于0的边上我们将获得饱和区的压力变化色彩条纹图。选择“Display”菜单下的“DisplayOptions”命令关掉“Elements”选项打开“Iso-su_____s”选项单击紧邻“Iso-su_____s”的“Options”按钮为第一个等值面输入值
0.0打开第一个等值面紧__“Cap”选项OK。选择Shade命令。
17.
12.3在地面上覆盖TIFF图像最后我们将TIFF图像覆盖在地面上用来说明计算水面与地面之间的空间__。首先我们显示顶部TIN。转到TIN模块选择SelectTINs工具选择顶部TIN选择Show命令在GMS窗口左上角的TIN组合框内选择terrain。转到__p模块选择“I__ges”菜单下的“DisplayOptions”命令打开“Texture__ptosu_____whenshaded”选项在组合框内选择“__ptoactiveTIN”选项OK。选择Shade命令。选择“Rotate”工具从不同角度进行观看。最后我们将平滑阴影选项:选择“Display”菜单下的“ShadingOptions”命令选择“__oothfeatures”选项OK。选择Shade命令。第18章SEEP2D承压SEEP2D是二维有限元稳定水流模型。它通常用于剖面模拟。可用于承压或潜水相应的分为两部分。
18.1问题描述如图
18.1:
18.2启动GMS
18.3所需模块和界面2DMesh模块SEEP2D界面。
18.4创建网格SEEP2D用两种类型的单元:三角形和四边形。问题简单时用现在介绍的方法若问题复杂则用__P模块法。
18.
4.1制定坐标系统如图
18.2所示坐标原点采用板桩上游90英尺处的基岩顶面点。我们可以通过在关键处创建节点大方法创建网格创建节点前我们必须设置窗口边界以致在我们的图形窗口中能看到实世界坐标系统的适当区域。选择“View”菜单下的“SetWindBounds”命令选择“Xrangetobespecified”选项设置“Xatleft”为-10设置“Xatright”为160设置“Yatbottom”为-10OK。我们还需要通过设置网格线来帮助我们将节点放置到适当的位置:选择“Display”菜单下的“DrawingGridOptions”命令设置“Gridspacing”为2打开“Snaptogrid”选项OK。
18.
4.2创建角点转到2DMesh模块从工具板上选择“CreateNode”工具在图形窗口内单击添加节点根据窗口底部的坐标值创建以下坐标值的点:如果点位不准确可以通过SelectNodes工具选择节点然后在GMS窗口顶部的坐标编辑框修改坐标。通过选中节点并按Delete键或选择“Edit”菜单下的“Delete”命令删除节点。单击FrameI__ge工具使图形居中显示。关闭网格显示:选择“Display”菜单下的“DrawingGridOptions”命令关掉“Snaptogrid”选项OK。创建的点如图
18.3所示:
18.
4.3插值选择SelectNodes工具选择相邻的两个节点选择“ModifyMesh”菜单下的“NodeInterpOptions”命令选择“Numberofintervalsinstring”选项并输入两点之间的间隔数之间的节点数是间隔数-1OK。插值后的图
18.
418.
4.4创建额外的点添加点来刻画板桩。从工具板上选择“CreateNode”工具添加点
90.
2532.0。
18.
4.5创建单元通常四边形单元比三角形单元更有效。选择“CreateLinear________Element”工具选择“BuildMesh”菜单下的“Triangulate”命令选择“ModifyMesh”菜单下的“Merge________s”命令OK。
18.
4.6删除片状板桩单元下一步要保证在单元中适当的表示板桩板桩为隔水边界因此需要在板桩位置的单元中添加一个裂缝。可以通过删除板桩位置的窄单元来实现。选择Zoom工具在板桩位置拖拉鼠标放大。显示如图
18.6。单击“SelectElements”工具选中PQR窄单元体按Delete键或选择“Edit”菜单下的“Delete”命令删除窄单元体。OK。选则FrameI__ge命令。
18.
4.7细化网格选择“ModifyMesh”菜单下的“RefineElements”命令OK。重复几次。
18.5重新对单元编号要看节点半带宽:选择“File”菜单下的“GetInfo”命令记下最大板带宽366OK。现在计算每个节点在单元中的ID号:完成节点和单元的构建:选择“Display”菜单下的“DisplayOptions”命令选中“Nodenumbers”选项OK。选择Zoom工具和FrameI__ge命令查看节点编号。要对节点和单元进行重新编号必须首先选择边界上的节点串。该串为重新编号指定起始点。选择SelectNodeString工具按从上到下的顺序依次选择最左列上的各节点方法1:先单击顶部节点然后依次向下单击节点;方法2:先单击顶部节点然后按下CTRL键单击底部节点选择“ModifyMesh”菜单下的“Renumber”命令。沿网格的最小轴进行节点编号则最大节点半带宽最小。取消显示节点编号:选择“Display”菜单下的“DisplayOptions”命令取消“Nodenumbers”选项。
18.6初始化SEEP2D模拟选择“SEEP2D”菜单下的“NewSimulation”命令。
18.7设置分析选项选择“SEEP2D”菜单下的“____ysisOptions”命令为Title输入SampleConfinedFlowProblem在“ModelType”区选择“Confined”单击“Units”按纽设置长度Length单位ft时间Time单位yr质量__ss单位slug=
32.2磅OK。
18.8指定材料性质每个区有一组属性值包括k1、k2为两个主渗透系数方向和angleX轴与渗透系数主轴方向在逆时针方向上的夹角如图
18.7。大多数土壤堆积物渗透系数主轴为X此次轴为Y角度为0。每个单元具有一个材料编号材料编号为材料属性列表的索引。由于模型只有一种岩性因此所有的单元具有同一的材料编号1。输入属性值:选择“SEEP2D”菜单下的“__terialProperties”命令为k1和k2输入值100OK。
18.9指定边界条件
18.
9.1定水头边界如图
18.8所示:单击SelectNodes工具选择粘土覆盖层左侧的顶部节点选择“SEEP2D”菜单下的“NodeBC”命令为“Head”输入水头值43。同样为板桩右侧的顶部单元赋值:要注意从板桩的右边单元
90.2532开始。
18.10保存模拟选择“SEEP2D”菜单下的“S__eAs”命令输入文件名tutorial\s2con\blanket.sps。
18.11运行SEEP2D选择“SEEP2D”菜单下的“RunSEEP2D”命令。
18.12查看结果选择“SEEP2D”菜单下的“ReadSolution”命令。选择“Display”菜单下的“DisplayOptions”命令选择SEEP2D选项卡打开“Title”和“Totalflowrate”选项。第19章SEEP2D潜水SEEP2D是二维有限元稳定水流模型。它通常用于剖面模拟。可用于承压或潜水相应的分为两部分。SEEP2D有两种方法用于非承压含水层问题。1只计算饱和区单元被切成潜水面。新的计算解组成新的网格单元文件和水头解文件。2单元不变形计算饱和和非饱和区的水流。
19.1问题描述如图
18.1:各向异性土坝并具有防渗心墙。
19.2启动GMS
19.3所需模块和界面2DMesh模块__p模块SEEP2D界面。
19.4创建网格在__p模块模式下用“FeatureO__ects”菜单下的“__p-2DMesh”命令。
19.
4.1制定坐标系统如图
19.2所示坐标原点采用坝的左下角。我们可以通过在关键处创建节点的方法创建网格创建节点前我们必须设置窗口边界以致在我们的图形窗口中能看到实世界坐标系统的适当区域。选择“View”菜单下的“SetWindBounds”命令选择“Xrangetobespecified”选项设置“Xatleft”为-10设置“Xatright”为370设置“Yatbottom”为-50OK。我们还需要通过设置网格线来帮助我们将节点放置到适当的位置:选择“Display”菜单下的“DrawingGridOptions”命令设置“Gridspacing”为2打开“Snaptogrid”选项OK。
19.
4.2创建图层转到__p模块选择“FeatureO__ects”菜单下的“Coverages”命令将“CoverageType”改为“SEEP2D”输入名称SEEP2DMesh。
19.
4.3创建角点转到2DMesh模块从工具板上选择“CreateNode”工具在图形窗口内单击添加节点根据窗口底部的坐标值创建以下坐标值的点:如果点位不准确可以通过SelectNodes工具选择节点然后在GMS窗口顶部的坐标编辑框修改坐标。通过选中节点并按Delete键或选择“Edit”菜单下的“Delete”命令删除节点。单击FrameI__ge工具使图形居中显示。关闭网格显示:选择“Display”菜单下的“DrawingGridOptions”命令关掉“Snaptogrid”选项OK。创建的点如图
19.3所示:
19.
4.4创建弧段创建弧段来指定边界。单击CreateArc工具沿坝一圈单击鼠标A-B-C-D-I-J-H-G-A创建心墙弧段G-E-F-H。
19.
4.5顶点重新分布我们的弧段只有一条边而无顶点。利用__p-2DMesh命令来控制顶点间隔。单击CreateArc工具拖拉鼠标选择所有弧段选择“FeatureO__ects”菜单下的“RedistributeVerti__s”命令选择“Targetspacing”选项输入spacing间距
4.0OK。
19.
4.6创建多边形选择“FeatureO__ects”菜单下的“BuildPolygons”命令在提示框内单击OK按钮。
19.
4.7指定材料类型单击SelectPolygons工具在坝体外壳内双击鼠标在弹出的对话框中其默认材料名为__terial#1改为Shell。双击心墙区在弹出的对话框中单击New按钮输入名称Core。
19.
4.8创建网格选择“FeatureO__ects”菜单下的“__p-2DMesh”命令在提示框内单击OK按钮。
19.5重新对单元编号现在计算每个节点在单元中的ID号:要对节点和单元进行重新编号必须首先选择边界上的节点串。该串为重新编号指定起始点。对于简单的土坝从最左或最右边的节点开始编号一般最优。所以只要选择一个节点。转到2DMesh模块选择SelectNodeString工具选择最左的节点选择“ModifyMesh”菜单下的“Renumber”命令。取消显示节点编号:选择“Display”菜单下的“DisplayOptions”命令取消“Nodenumbers”选项。
19.6初始化SEEP2D模拟选择“SEEP2D”菜单下的“NewSimulation”命令。
19.7设置分析选项选择“SEEP2D”菜单下的“____ysisOptions”命令为Title输入SampleUnconfinedFlowProblem在“ModelType”区选择“Unconfined”选择“Saturatedwithdeformingmesh”选项单击“Units”按纽设置长度Length单位ft时间Time单位yr质量__ss单位slug=
32.2磅OK。
19.8指定材料性质每个区有一组属性值包括k1、k2为两个主渗透系数方向和angleX轴与渗透系数主轴方向在逆时针方向上的夹角如图
18.7。每个单元具有一个材料编号材料编号为材料属性列表的索引。输入属性值:选择“SEEP2D”菜单下的“__terialProperties”命令:材料1Shell:为k1和k2输入值15060材料2Core:为k1和k2输入值
15619.9指定边界条件此处有三种边界条件:隔水边界定水头边界渗出面。
19.
9.1指定水头边界如图
19.4所示:选择SelectNodeString工具选择坝上游与水面交点Y=60ft按下CTRL键单击左下角点y=0则水面以下坝上游面节点都被选中选择“SEEP2D”菜单下的“NodeBC”命令为“Head”输入水头值60;选择SelectNodeString工具选择单击右下角点y=0按下CTRL键选择坝下游与水面交点Y=6ft则水面以下坝上游面节点都被选中选择“SEEP2D”菜单下的“NodeBC”命令为“Head”输入水头值6。
19.
9.2渗出面边界条件由于采用改变网格的方法SEEP2D通过迭代求得潜水面的位置在迭代过程中网格的形状发生改变以使网格顶面与潜水面一致。为了引导迭代过程我们需要将潜水面存在的节点标记为渗出面节点。潜水面可能存在的区域如图
19.5所示:选择渗出面的两端点选择“SEEP2D”菜单下的“NodeBC”命令选择“ExitFa__”选项。
19.10保存模拟选择“SEEP2D”菜单下的“S__eAs”命令输入文件名tutorial\s2unc\dam
1.sps。
19.11运行SEEP2D选择“SEEP2D”菜单下的“RunSEEP2D”命令。
19.12查看结果选择“SEEP2D”菜单下的“ReadSolution”命令。在提示框内的单击“OK”按钮删除网格。选择“Display”菜单下的“DisplayOptions”命令关闭“Nodes”和“Elements”选项选中“Meshboundary”选项。选择“Display”菜单下的“DisplayOptions”命令选择SEEP2D选项卡在“Base__terial”下拉列表框内选择材料CoreOK。
19.13模拟非饱和区的水流
19.
13.1读入原始模型选择“SEEP2D”菜单下的“ReadSimulation”命令在提示框内的单击“No”选择文件dam
1.spsOpen。OK。
19.
13.2改变网格的显示选择“Display”菜单下的“DisplayOptions”命令打开“Nodes”和“Elements”选项。
19.
13.3改变分析选项SEEP2D中有两种方法计算非饱和区的相对渗透系数:线性锋面法和VanGenuchten法。线性锋面法相对渗透系数从饱和值线性下降到给定负压水头下的用户指定最小值。VanGenuchten法VanGenuchten参数用来定义非饱和区的相对渗透系数。详见SEEP2DPrimer。此处用线性锋面法。选择“SEEP2D”菜单下的“____ysisOptions”命令输入标题SampleUnconfinedProblemPartII在Unconfined区选择“Saturated/unsaturatedwithlinearfront”选项OK。
19.
13.4编辑材料属性需要输入相对渗透系数和压力水头阈值。选择“SEEP2D”菜单下的“__terialProperties”命令:材料1Shell:ho-
1.0kro
0.001为k1和k2输入值15060材料2Core:ho-
4.0kro
0.001为k1和k2输入值
15619.14保存模拟选择“SEEP2D”菜单下的“S__eAs”命令输入文件名tutorial\s2unc\dam
2.sps。
19.15运行SEEP2D选择“SEEP2D”菜单下的“RunSEEP2D”命令。
19.12查看结果由于网格不变形产生包含水头、流速等的数据文件。选择“SEEP2D”菜单下的“ReadSolution”命令。选择“Display”菜单下的“DisplayOptions”命令关闭“Nodes”和“Elements”选项。另一种用来查看饱和/非饱和解的有效方法是创建潜水面图。通过压力水头等值线和等于0的等值线获得。选择“Display”菜单下的“DisplayOptions”命令选择SEEP2D选项卡关闭“Flowlines”选项OK。在GMS窗口顶部的数据集组合框内选择pressure_head数据选择“Data”菜单下的“ContourOptions”命令在ContourInterval区内选择“Specifiedvalues”选项单击“Specifiedvalues”右边的“Values”按钮输入第一条等值线的值
0.0OKOK。显示饱和/非饱和带的分界线。第20章UTCHEMUTCHEM模拟地下水流和溶质运移的成熟软件还可以模拟多相状态。UTCHEM是UniversityofTexasatAustin大学的__nterforPetroleumandGeosystemsEngineering__的。UTCHEM是一个三维化学流体模拟器。特别适合于抽水冲洗抽水处理包括表面活性剂加强的冲洗。
20.1问题概述本教程要解决的问题如图
20.1所示。此问题是UTCHEM文档中所介绍的ex21和ex22中所涉及的问题的简化。该问题是一个中间有P__溢出的X-Z方向的切片。处理方法是:在溢出处开挖一条沟两边有两条冲洗沟。背景渗流场为由左向右的单方向流。此问题可分为两部分:污染状态和处理状态。污染状态持续30天包括中上部的P__溢出源。处理状态包括一个60天的等待过程和持续630天抽出处理过程。处理方法是直接在溢出源上部放置开采井和在污染区的左侧和右侧放置两个注水井。
20.2启动
20.3所需模块和界面3DGrid模块UTCHEM界面。
20.4创建网格3DGrid模块选择“Grid”菜单下的“CreateGrid”命令在“X-dimension”区内输入长度值
160.0和单元数49在“Y-dimension”区内输入长度值
50.0和单元数1在“Z-dimension”区内输入长度值
39.0和单元数24OK。选择ViewIAxis命令。
20.5污染状态首先我们模拟P__从表面溢出可用在模型中心的单个注水井中以稳定流量连续注入污染物30天。
20.
5.1初始化UTCHEM模拟选择“UTCHEM”菜单下的“NewSimulation”命令。
20.
5.2一般选项首先我们输入标题和基本信息:为“Runnumber”输入“EX21”为“Title1”输入“EX21:2dX-ZP__contaminationinthesaturatedzone”为“Title2”输入“30dayP__spill”为“Title3”输入“__姓名和日期”保证默认语言为英语。从组合框内为数值弥散控制选项numericaldispersioncontrolchoose选择两点上游方法Twopointupstreammethod。然后制定模拟的物种4种污染物:水、石油、氯化物和钙。单击“DefineSpecies…”按钮打开“DefineSpecies”对话框取消“su___ctant”、“alcohol1”和“alcohol2”选项在名称列表框内只留下以上四种物质water、oil、chloride和calcium选择“oil”并将名称改为“P__”OKOK。
20.
5.3时间选项选择“UTCHEM”菜单下的“TimeOptions”命令保证“Output”和“__xInject”的选项为“Days”为“timestepsize”选择“Auto__tic”。注入时段:UTCHEM的模拟时间域分为注入时段。出水量和边界条件值只能在一个注入时段的开始时改变。在此模拟中的污染阶段我们将设置为一个长为30天的注入时段。在对话框内单击“InjectionPeriods…”按钮为“Length”输入
30.0。时间步长:当用字段时间步长选项时UTCHEM允许最大和最小时间步长受Courantnumbers和days的控制详见UTCHEM用户手册的附录E。初始时间步长总是为天的单位。为“Initialtimestepsize”输入
0.03为“Tol.forconc.change”输入
0.003在Timesteps区选择“CourantNumbers”作为最大最小时间步数的单位为“__xtimestepsize”输入
0.2为“Mintimestepsize”输入
0.04。输出间隔控制:“OutputIntervals”决定保存什么值最重要的是“Profiles”的值这决定以什么样的间隔将浓度值输出到GMS的数据文件:为“Printdata”输入
15.0为“Prod.data”输入
15.0为“Prod.histories”输入
5.0为“Profiles”输入
5.0OKOK。输出控制:“OutputControl”允许我们选择从UTCHEM中保存什么数据。选择“UTCHEM”菜单下的“OutputControl”命令。在“Profiles”区包括可以保存在GMS数据表的形式默认值为选中。在“Print”区包括可以以文本编辑器浏览的UTCHEM的重要输出数据。在“Profiles”区我们将输出“Phasepressures”、“Phasesaturations”和“Componentcon__ntrations”;在“Print”区关掉所有选项即默认值。OK。含水层性质ReservoirProperties:在“ReservoirProperties”区我们指定含水层特性。一些特性如孔隙度和压力我们输入成数组的形式每个单元一个值。其他的属性则采用整个区一个值如氯化物和钙。对于本问题假定所有的值:选择“UTCHEM”菜单下的“ReservoirProperties”命令确保RockCompressibility值为
0.0和“Referen__pressure”值为
14.7。指定孔隙度:单击“Porosity”按钮单击“Constant-Grid”按钮并输入值
0.34OK。初始条件:单击“Pressure”按钮单击“Constant-Grid”按钮并输入值
14.7OK。单击“WaterSat…”按钮单击“Constant-Grid”按钮并输入值
1.0OK。将“Chloride”和“Calcium”的值改为
0.009。渗透性:为“Kx”输入值8300darcies。单击“Kx”按钮单击“Constant-Grid”按钮并输入值8300darciesOK。将“Ky/Kx”和“Kz/Kx”的值分别改为
1.0和
0.5OK。要编辑单元的属性值可以选中单元选择“UTCHEM”菜单下的“__llAttributes”命令。
20.
5.6物理属性现在输入物理属性值:相对渗透性粘性密度可压缩性弥散性和吸附参数。相态:相态行为主要与表面活性剂和助溶剂有关但在此处不用。然而必须输入如下数据:选择“UTCHEM”菜单的“PhysicalProperties”子菜单下的“PhaseBeh__ior”命令选择“Con__ntrations”选项卡将“Oilcon__ntrationatplaitpoint+”改为
0.0将“Oilcon__ntrationatplaitpoint-”改为
1.0将“Criticalmi__llecon__ntration”改为
0.0001OK。界面张力:选择“UTCHEM”菜单的“PhysicalProperties”子菜单下的“Inte___cialTension”命令选择“HealyReed’scorrelation”选项为“Water/micro-emulsion”输入参数:G
113.0G2–
14.8和G
30.007;“Oil/micro-emulsion”输入参数:G
113.0G2–
14.5和G
30.01;设置“XIFTW”为
1.65OK。有机溶质运移:编辑有机质量传输参数:选择“UTCHEM”菜单的“PhysicalProperties”子菜单下的Organic__ssTransfer”命令选择“Allowsolubilityofoilinwaterinabsen__ofsu___ctant”选项为“Equilibriumconc.ofoilinwater”输入
0.00015OK。相对渗透性:选择“UTCHEM”菜单的“PhysicalProperties”子菜单下的“RelativePermeability”命令选择“Firstdrainagecorey”选项在“Aqueous”栏中单击“ResidualSaturation”按钮单击“Constant-Grid”按钮并输入值
0.24OK;单击“KrEndpoint”按钮单击“Constant-Grid”按钮并输入值
1.0OK;单击“KrExponent”按钮单击“Constant-Grid”按钮并输入值
2.85OK;在“Oleic”栏中单击“KrEndpoint”按钮单击“Constant-Grid”按钮并输入值
0.65OK;单击“KrExponent”按钮单击“Constant-Grid”按钮并输入值
2.7OK;在“Microemulsion”栏中单击“ResidualSaturation”按钮单击“Constant-Grid”按钮并输入值
0.24OK;单击“KrEndpoint”按钮单击“Constant-Grid”按钮并输入值
1.0OK;单击“KrExponent”按钮单击“Constant-Grid”按钮并输入值
2.85OK。粘滞性:选择“UTCHEM”菜单的“PhysicalProperties”子菜单下的“Viscosity”命令将“Waterviscosity”改为
1.0将“Oilviscosity”改为
0.__确保“Referen__temp.”为
0.0。在“Microemulsionviscosity”区内将“Alpha1”改为
3.4将“Alpha2”改为
1.0将“Alpha3”改为
3.0将“Alpha4”改为
1.0将“Alpha5”改为
1.0OK。组分密度:选择“UTCHEM”菜单的“PhysicalProperties”子菜单下的“ComponentDensity”命令为“Water”输入
0.433为“Oil”输入
0.7036为“Coef.Ofoilinmicroemulsion”输入
0.7036跳过“Su___ctant”、“Alcohol1”和“Alcohol2”选择“Considergr__ityeffect”选项OK。流体压缩性:不可压缩将各种液体的“fluidcompressibility”设置为0。毛细压力:选择“UTCHEM”菜单的“PhysicalProperties”子菜单下的“CapillaryPressure”命令选择“Stronglywaterwet”选项单击“Cap.Press.Endpoint”按钮单击“Constant-Grid”按钮并输入值
1.4OK;单击“Cap.Press.Exponent”按钮单击“Constant-Grid”按钮并输入值–
0.52OKOK。扩散和弥散选择“UTCHEM”菜单的“PhysicalProperties”子菜单下的“Diffusion/Dispersion”命令在“Dispersivity”区中将“Aqueousphase”的“Longitudinal”和“Transverse”分别改为
0.1和
0.03重复最后两步为“Oleic”和“Microemulsion”赋值OK。吸附本次模拟忽略。
20.
5.7边界选择“UTCHEM”菜单下的“Boundaries”命令确保“Zonetobemodeled”区内的“Saturatedonly”被选中选中“Leftboundaryisspecified”和“Rightboundaryisspecified”:左边界:在左侧栏单击“Pressureattoplayer”标题右侧的带有数字的图形框弹出对话框在“pressure”标题下的编辑框中输入
15.028OK;单击“Conc.ofwater”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的编辑框中输入
1.0OK;单击“Conc.ofchloride”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的编辑框中输入
0.009OK;单击“Conc.ofcalcium”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的编辑框中输入
0.009OK。右侧边界在右侧栏单击“Pressureattoplayer”标题右侧的带有数字的图形框弹出对话框在“pressure”标题下的编辑框中输入
14.7OK;单击“Conc.ofwater”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的编辑框中输入
1.0OK;单击“Conc.ofchloride”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的编辑框中输入
0.009OK;单击“Conc.ofcalcium”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的编辑框中输入
0.009OK。
20.
5.8井选择“UTCHEM”菜单下的“Wells”命令单击“New”按钮输入名称P__spill为“well”选择“Rateconstrainedinjectionwell”类型将“Column”改为25将“Skinfactor”改为
0.0从“Phases”列表框中选择“oleic”单击“Q”标签右侧的图形框在“flow”标题下的编辑框中输入
33.33OK。从“Components”列表框中选择“P__”单击“Conc.”标签右侧的图形框在“con__ntration”标题下的编辑框中输入
1.0OK。
20.
5.9保存模型选择“UTCHEM”菜单下的“S__eAs”命令输入文件名contamination.uts。
20.
5.10运行选择“UTCHEM”菜单下的“RunUTCHEM”命令OK。
20.
5.11查看结果选择“UTCHEM”菜单下的“ReadSolution”命令OK。UTCHEM的解包含每种组分在每一状态下的的浓度值数据集及相应的压力和饱和度数据集。我们主要研究P__它是第二种组分。从“DataSet”组合框内选择“contamination_c002”从“TS”组合框内选择最后一个时间步。由于P__是DNPL所以它迅速下降到含水层的底部。
20.
5.12改变等值线选项选择“Data”菜单下的“ContourOptions”命令打开“Usecolorramp”选项OK。
20.6处理阶段添加井表示处理沟。一个排水沟在溢流处两个注水沟位于溢流处的两侧。在该阶段开始关掉P__源并在打开沟前允许溢出的P__稳定60天。沟将运营2年。为了确保获得适当的初始条件我们将在此模拟中包含污染阶段的30天。这也确保我们有一个完整的后处理过程。
20.
6.1时间选项首先我们创建额外两个注入阶段第一个模拟停滞阶段第二个模拟处理阶段。选择“UTCHEM”菜单下的“TimeOptions”命令在对话框内单击“InjectionPeriods…”按钮。60天的等待阶段:单击“InsertRight”按钮添加另一个注入阶段在“Injectionperiod”编辑区域单击向上或向下箭头选择时段为“Length”输入
60.0。在“Timesteps”编辑域中选择“Days”为最大和最小时间步数的单位。为“Initialtimestepsize”输入
0.001为“Tol.forconc.change”输入
0.003在Timesteps区选择“CourantNumbers”作为最大最小时间步数的单位为“__xtimestepsize”输入
0.28为“Mintimestepsize”输入
0.014。“OutputIntervals”决定保存什么值最重要的是“Profiles”的值这决定以什么样的间隔将浓度值输出到GMS的数据文件:为“Printdata”输入
60.0为“Prod.data”输入
60.0为“Prod.histories”输入
5.0为“Profiles”输入
5.0OKOK。处理阶段:单击“InsertRight”按钮添加另一个注入阶段在“Injectionperiod”编辑区域单击向上或向下箭头选择时段为“Length”输入
630.0。在“Timesteps”编辑域中选择“CourantNumbers”为最大和最小时间步数的单位。为“Initialtimestepsize”输入
0.001为“Tol.forconc.change”输入
0.003在Timesteps区选择“CourantNumbers”作为最大最小时间步数的单位为“__xtimestepsize”输入
0.15为“Mintimestepsize”输入
0.01。“OutputIntervals”决定保存什么值最重要的是“Profiles”的值这决定以什么样的间隔将浓度值输出到GMS的数据文件:为“Printdata”输入
90.0为“Prod.data”输入
90.0为“Prod.histories”输入
5.0为“Profiles”输入
30.0OKOK。
20.
6.2边界选择“UTCHEM”菜单下的“Boundaries”命令确保“Zonetobemodeled”区内的“Saturatedonly”被选中选中“Leftboundaryisspecified”和“Rightboundaryisspecified”:左边界:在左侧栏单击“Pressureattoplayer”标题右侧的带有数字的图形框弹出对话框在“pressure”标题下的所有编辑框中输入
15.028OK;单击“Conc.ofwater”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的所有编辑框中输入
1.0OK;单击“Conc.ofchloride”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的所有编辑框中输入
0.009OK;单击“Conc.ofcalcium”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的所有编辑框中输入
0.009OK。右侧边界在右侧栏单击“Pressureattoplayer”标题右侧的带有数字的图形框弹出对话框在“pressure”标题下的所有编辑框中输入
14.7OK;单击“Conc.ofwater”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的所有编辑框中输入
1.0OK;单击“Conc.ofchloride”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的所有编辑框中输入
0.009OK;单击“Conc.ofcalcium”标题右侧的带有数字的图形框弹出对话框在“con__ntration”标题下的所有编辑框中输入
0.009OK。
20.
6.3井选择“UTCHEM”菜单下的“Wells”命令左侧注水井:单击“New”按钮输入名称LeftInjectionWell为“well”选择“Rateconstrainedinjectionwell”类型将“Column”改为16“Row”保留1改“Endinglayer”为24将“Skinfactor”改为
0.0;从“Phases”列表框中选择“Aqueous”单击“Q”标签右侧的图形框在“flow”标题下的三个编辑框中输入
1125.0OK;从“Components”列表框中选择“Water”单击“Conc.”标签右侧的图形框在“con__ntration”标题下的第三个编辑框中输入
1.0OK;从“Components”列表框中选择“Chloride”单击“Conc.”标签右侧的图形框在“con__ntration”标题下的第三个编辑框中输入
0.009OK;从“Components”列表框中选择“Calcium”单击“Conc.”标签右侧的图形框在“con__ntration”标题下的第三个编辑框中输入
0.009OK。开采井:在排水渠处创建开采井单击“New”按钮输入名称__nterProductionWell为“well”选择“RateConstrainedProductionWell”类型将“Column”改为25“Row”保留1改“Endinglayer”为24将“Skinfactor”改为
0.0;在“Phases”列表框中单击“Q”标签右侧的图形框在“flow”标题下的三个编辑框中输入-
2250.0OK。右侧注水井:单击“New”按钮输入名称RightInjectionWell为“well”选择“Rateconstrainedinjectionwell”类型将“Column”改为41“Row”保留1改“Endinglayer”为24将“Skinfactor”改为
0.0;从“Phases”列表框中选择“Aqueous”单击“Q”标签右侧的图形框在“flow”标题下的三个编辑框中输入
1125.0OK;从“Components”列表框中选择“Water”单击“Conc.”标签右侧的图形框在“con__ntration”标题下的第三个编辑框中输入
1.0OK;从“Components”列表框中选择“Chloride”单击“Conc.”标签右侧的图形框在“con__ntration”标题下的第三个编辑框中输入
0.009OK;从“Components”列表框中选择“Calcium”单击“Conc.”标签右侧的图形框在“con__ntration”标题下的第三个编辑框中输入
0.009OK。
20.
6.4保存模型选择“UTCHEM”菜单下的“S__eAs”命令输入文件名treatment.uts。
20.
6.5运行运行需要一段时间可以通过预先计算的解查看结果。用预先计算的解:选择“UTCHEM”菜单下的“ReadSolution”命令单击“SelectSimulation”按钮打开文件tutorial/utchem/sample/utchem
2.utsOpen。计算新的解:选择“UTCHEM”菜单下的“RunUTCHEM”命令OK。读入解:选择“UTCHEM”菜单下的“ReadSolution”命令OK。UTCHEM的解包含每种组分在每一状态下的的浓度值数据集及相应的压力和饱和度数据集。我们主要研究P__它是第二种组分。从“DataSet”组合框内选择“treatment_c002”从“TS”组合框内选择最后一个时间步。设置影像循环:选择“Data”菜单下的“FilmLoop”命令单击“Setup”按钮打开“Displayclock”选项OK。单击“Run”按钮开始动画显示单击“Stop”按钮停止动画单击“Step”按钮则一次只__一幅最后单击“Done”按钮退出。PAGE1。