gms6 0tutor vo2翻译modflow模块在GMS中有两种应用方式可以构建模拟栅格和概念

CHAPTERMODFLOW-GridGMS中，有两种应用方式可以构建modflow模拟：栅格和概念模型。栅格方式包括直接在3D栅格的逐个单元（cell-by-cell）添加源/汇项和其他模型参数。概念模型应用方式包括在map模块中使用gis工具开发一个被模拟区域的概念模型。概念模型中的数据将会被到栅格中。MODFLOWConceptualModelApproachtutorial.ThisiswhatyouwillCreatea3DSetupaMODFLOWCheckthesimulationandrunAssignzonebudgetsandviewtheRequiredYouwillneedthefollowingcomponentsenabledtocompletethisYoucanseeifthesecomponentsareenabledbyselectingtheFile Descriptionof 此指南中需要解决的问题在图2.1中描述。此问题是MODFLOWReference（图中没有标示）0，稳Getting 选择Edit|Units，Time点击okCreatingthe 3D在ProjectExplorer中，空白出，在下拉菜单中选择New|3DGrid命度对于MODFLOW的计算不起作用。选择okMini-GridPlot.CreatingtheMODFLOW 下一步是建立modflow1.在ProjectExplorer中，3DGridData ，在下拉菜单中选择NewTheGlobal 这个输入是为了细分modflow为块。一些块是备用的，一些是必需的。必需的 这个块框用于指定建立模型所需要的块基础的块是必需的不能被关闭。在Pointsources/sinks部分，打开DrainDRN1)Well(WEL1)在Arealsources/sinks部分，打开RechargeRCH1)选项在SolverStonglyImplProcedureSIP1)选择ok推出框TheIBOUND IBOUND框在一个电子表格中显示IBOUND数列的值，每次显示一层。左上角的0外，我们需要其他所有的数30个定水头单元格中输入-1即可。尽管如此，还有另外更简单的方法修改IBOUND数列值。这将会在后续的选择ok退出框。StartingHeads 下一步建立开始水头数列选择StartingHeads选择ok退出框TopandBottomElevations 选择TopElevation确定当前层为1（Layeris60选择ok退出TopElevations框选择BottomElevation确认当前层为1输入-452Constant→Layer输入-1203Constant→Layer输入-215选择ok退出BottomElevation框AssigningIBOUNDValuesDirectlyto 综上所述，IBOUND值可以通过IBOUND数列框进行输入。从某些情况下说，这比直接对每个单元进行输入更为容易。可以通过CellProperties命令来完成。ViewingtheLeft 栅格非常薄。这可以通过调整Z选择Disy选将Zmagnification修改为15选择okSelectingthe 选择SelectCells工具ChangingtheIBOUNDValue 编辑IBOUND值在所选择的其中一个单元上选择ok退出CellProperties框选 nView按 CheckingtheValues 确保输入的IBOUND值的准确度。选择MODFLOW|GlobalOptions在左上角的框中选择右侧的向上箭头，让图层循环显示注意上面两层最左边一列单元的值均为-1GMS中使用电子表格方式或者是通过选择单元并直接输入（选择一个最为便利的输入方式）MODFLOW数据选择ok按钮退出IBOUNDArray框选择ok按钮退出MODFLOWGlobalPackage框TheLPF （LPF1.选择MODFLOW|LPFPackageLayer 在框中LayerData部分的选项被用于定于图层类型以及每一个图层的水力传导系数。在我们的问题中，有三个图层。最上面的一层是的，下面两层是有限制的（承Layer TopLayer 选择HorizontalHydraulicConductivity15选择ok选择ok按钮退出HorizontalHydraulicConductivity框10Middle 在LayerData部分选择向上箭头转到正确的图层编辑区域，转到图层22Bottom 3选择ok退出MODFLOWLPFPackage框TheRecharge 补给（回灌）选择MODFLOW|Source/SinkPackages|RechargePackage0.0009并点击ok选择ok按钮退出RechargePackage框TheDrain 格，再选择PointSources/Sinks菜单命令。Selectingthe 182-10列单元格。选择单元格：选择SelectCells工具 框中显示出来，如图2-3.选择i=8,j=2,k=1.按下shift3－102-Assigningthe 在所选单元格的图形窗口中，在下拉菜单中选择Sources/Sinks菜单选择okTheWell 下面使用PointSources/SinksTopLayer 按下shift2-5well选择new为所有的井输入一个flow值为-12,230（负值表示抽水选择okMiddleLayer 在Mini-GridPlot.中选择Decrement按钮。按下shift2-6well选择new为所有的井输入一个flow值为-12,230（负值表示抽水选择okBottomLayerWell 在Mini-GridPlot.中选择Decrement按钮。2-7well选择new为所有的井输入一个flow值为-0.15（负值表示抽水选择ok至此，所有的已定义，返回上部图层。在Mini-GridPlot.选择两次Checkingthe MODFLOWModelChecker的目的是分析定义当前MODFLOW模拟所输入的数据并给运行Model选择MODFLOW|CheckSimulation选择RunCheck选择Done按钮退出ModelChecker.框Savingthe 选择File|SaveAs移动到 下Running 选择MODFLOW|RunMODFLOWViewingthe ChangingLayers 在Mini-GridPlot.中选择向下箭头。选择向下箭头。选择向上箭头两次。ColorFillContours 选择Data|ContourOptions改变ContourMethod为ColorFill.选择okColorLegend 显示一个颜例。选择Data|ColorRampOptions选择okZone 区域预算是地质局USGS开发的，用于计算水流动模型中子区域AssigningZoneBudgetIds 选择SelectCells工具如果需要，转换 nView显示确保当前层为3D栅格中的第二层图层。在Mini-GridPlot在图层2中所有单元一个方框输入2为ZonebudgetID并选择ok通过在Mini-GridPlot中选择向下箭头转向图层34-73ZonebudgetIDViewingtheZoneBudgetReport 选择Data|FlowBudget选择Zones（/显示区域3的报告。做完后选择ok按钮退出框 这是MODFLOWGridApproach参数。GMS不会自动做任何的单位转换。MODFLOW3D可以通过选择MODFLOW|GlobalOptions菜单命令并选择块功能来调用你所栅格单元并应用MODFLOW|CellProperties菜单命令。令或者在所选择的单元格的下拉菜单中选择Sources/Sinks菜单命令来建立可以使用ModelChecker可以通过使用Data|FlowBudget菜单命令来指定区域预测ID给栅格来产生一3DCHAPTERMODAEM是一个单图层是GMS中加强的恒态分析单元水流动模型本MODAEMGMS中的应用。本指南并未详细解释单元分析方法。需要详细了解单元分析模型和MODAEMGMShelp。在此之前你必需先完成GMSBasics运行不同条件下的MODAEMRequiredYouwillneedthefollowingcomponentsenabledtocompletethisYoucanseeifthesecomponentsareenabledbyselectingtheFileDescriptionof 本章描述GMS水流动模型在BigWalnutCreek(seeFigure1).泛滥平原的Brazil,IN.Brazilpopulation8188)地区小镇东部5SettingandData 井场位于BigWalnutCreek.24.2m600ft（183m60ft（18.3mGetting Let’sget1.Ifnecessary,launchGMS.IfGMSisalreadyrunning,selecttheFilecommandtoensurethattheprogramsettingsarerestoredtotheirdefaultReadingintheBackground 查找并打 将Filesoftype改变为(*.tif*.tiff*.jpg;*.jpeg;*.png;*.sid).在提示框中选择yes，产生信息图示（角锥此时，信息记录框出现。如果信息文件为地球（？？）相关文件或者文件包含有地质信息，信息记录框将不会出现。此信息并非为包含地质信息，所以在RegisterImage框中选择2PointRegistration选项8.选择ok8.选择okDefiningthe GMS界面中应用于选择Edit|Units选择okCreatingtheConceptual 以下输入模型数据。首先，建立一个MODAEM概念模型。建立一个覆盖并定在ProjectExplorer框中空白区域，在下拉菜单中选择New|ConceptualModel菜单命令。 在na概念模 上并选择NewCoverage菜单命令将Coveragename改变为Boundary并选择Usetodefinemodelboundaryok选择Zoom工具3-3选择CreateArc工具3-3CreatingtheSpecifiedHead MODAEM默认的边界是零流量边界。意思是弧段类型为指定流量并且流量为ConvertVerticesto 选择SelectObjects工具选择下图所示的四个点（顶点和/或节点shift键来选择超过如果需要增加输入顶点，可以通过CreateVertex工具来实现。选择工具并在选择FeatureObjects|VerticesNodesAssigning 选择SelectArcs工 3.选择3.选择Properties。在All一行将类型改变为spechead.选择ok退出框选择SelectPoints/Nodes工具选定北部定水头弧段的两个节点（通过按住shift或者是在两个节点周围拉动 10－12（6－9）178.6EnteringtheAquifer 选择MODAEM|GlobalOptions在Backgroundaquiferproperties选项中输入170.0作为Base，厚度Thickness为18.0，Hyd.cond.为60.0。选择ok退出框选择FeatureObjects|BuildPolygonsSavingthe 当你运行MODAEM，当前在内存中的数据将被写入临时文件，MODAEM将会读入并计算出解。因此，在运行MODAEM之前不需要保存修改。但是，定时保存有选择File|SaveAs在模型进行过程中定时的点击按钮Running 运行可以通过选择MODAEM|Solve菜单命令或者是通过按F5键运行MODAEM命令会弹出一个框，显示MODAEM模型的输出。点击F5当MODAEM运行完成后，点击close按钮。选择 yOptions按钮将Contourmethod改变为LinearandColorFill.在Fill选项下的框中将Transparency值改变为0.4在Filloptions选项中选择ColorRamp在PaletteMethod选项中的框，打开Legend重复选择ok退出框Creatingthe 在ProjectExplorer ，并选择NewCoverage按钮将名字改变为River4.选择ok4.选择ok选择CreateArc工 图3-6所示。注意，河流 选择SelectArcs工 选择Properties按钮选择SelectPoints/Nodes工具选择ok退出框HeadElev178.6175.6Running 运行再次运行MODAEM点击F5当MODAEM运行结束后，选择closeAdding 在CopyofBoundary在框的Sources/Sinks/BCs选项中，去掉SpecifiedHead和SpecifiedFlow.在ArealProperties选项中，打开Recharge6.选择ok6.选择ok选择SelectPolygons工 在多边形上双击并在Recharge区域输入值0.000420。点击ok退出框Running 运行再次运行MODAEM点击F5当MODAEM运行结束后，选择closeProduction 现在从一个指定的文本文档中输入生产井在ProjectExplorer ，并选择NewCoverage按钮选择ok退出框选择open按 查找并打 选中Headingrow点击Next>将GMSdatatype改变为Welldata.在Filepreview区域中，将第一，二，三列的type改变为：X，Y，FlowRate.选择finish按钮退出框 y符号，在其Dis Observation 在再次运行MODAEM在ProjectExplorer ，并选择NewCoverage按钮在ObservationPoints区域中打开Head选择ok退出框选择open按 查找并打 选中Headingrow点击Next>将GMSdatatype改变为Observationdata.在Filepreview区域中，将第一，二，三列的type改变为：X，Y，Obs.Head选择finish按钮退出框Running 运行再次运行MODAEM点击F5当MODAEM运行结束后，选择close 包括MODAEMMODAEM是GMSCHAPTERMODFLOW-ConceptualModel在GMS中，有两种应用方式可以构建modflow模拟：栅格和概念模型。栅格方应用的有关步骤已经在MODFLOWGridApproach.中进行描述。概念模型应用方式包参数力传导系数，模型边界和其他所有模拟需要的数据都可以在概念模型水平建议在完成此指南前先完成InterpolatingLayerElevationsandMODAEM RequiredYouwillneedthefollowingcomponentsenabledtocompletethisYoucanseeifthesecomponentsareenabledbyselectingtheFileDescriptionof MT3DMS指南中流动区域的粒子运动轨迹以及污染物运移模拟。模拟北部丘陵脚下沉积边界和南部两条河流交界之间区域的水流动。4-1b显示了一条南北的剖面此区域下伏基岩为石灰石在北区丘陵山丘上有露头。具有两个原生沉积层。上层作为含水层，下层作为有限含水层。对应假设首先进入系统的是降雨入渗区域内有一些季节性河流也会补给水。注意：虽然模拟区域为一个实际区域，填埋场和水文地质条件是假设的。Getting GMSImportingtheBackground Readingthe 选择open按 查找并打 GMSSavingthe 选择File|SaveAs往下可以随时点击按钮保存模型进展Definingthe GMS界面中应用于选择Edit|Units选择okDefiningthe 首先定义模型的边界。可以通过创建一个弧段并形成圈闭区域来完成Createthe 在ProjectExplorer中的空白位置，在下拉菜单中选择New|Conceptual点击ok在EastTexas概念模型中 ，并在下拉菜单中选择NewCoverage菜将名字Coveragename改变为Boundary。将Defaultelevation改变为213；将Defaultlayerrange改变为1-2。选择ok退出框Createthe 选择CreateArc工具在模型的左边（西）点击一次开始弧段，如图-22的基岩出露想一致。BuildingtheLocalSource/SinkCoverage 在创建特征目标前，可以在图层创建框中改变选项。 层中并在下拉菜单中选择Duplicate菜单命令。将新图层名字修改为Sources&Sinks.。在Sources&Sinks. 层中并在下拉菜单中选择CoverageSetup菜单LayerRefineSpecified确定Usetodefinemodelboundaryactivearea)选择okDefiningtheSpecifiedHead 选择SelectVertices工具4-312位于河流在所选的任何一个顶点上并选择Vertex->Node命令选择CreateArc工具选择模型南部和东部（右边和下边）在所选的任意一条弧段上并在下拉菜单中选择AttributeTable菜单找到all行和typespechead选择ok选择SelectPoints/Nodes工具选择ok同理，为河流交接处的节点输入值208DefiningtheDrainArcs 选择CreateArc工 4-41。从弧段的底部开始，通过沿河床点击点来创建，4-42、3选择SelectArcs工具通过按下shift在所选弧段上并在下拉菜单中选择AttributeTable菜单命令在all555作为传导系数。这指示了传导系数的单位长度值。当季节性河流被分配给栅格单元时，GMS自动计算恰当的单元传导系数值。将每一个弧段的Fromlayer和Tolayer改变为1选择ok选择SelectPoints/Nodes工具4-52为drain特性的Bot.elev.输入212。而不改变spechead其他的任何特性。点击ok确定仅改变drain特性，而不改变specheadBuildingthe 在区域源汇项图层中，模拟的全部区域都必需用不的多边形覆盖。定义栅1.选择FeatureObjects|BuildPolygons注意：多边形是被填充的。可以通过选择Disy|DisyOptions菜单命令来改Creatingthe 选择CreatePoint工具4-11（8351425选择Properties按钮为Flowrate,输入值-680将Fromlayer和Tolayer都改变为1选择oklayerTolayer2，意味着井只应用于栅格的第二层。Grid 井是水流动的点，并且在井附近引起水头陡坡（水跃，即水头损失。GMS中可以通过在概念模型中的井直接指派精确度来自动完成。选择SelectPoints/Nodes工具通过按下shift选择Properties按钮将两个点的Basesize改变为25，Bias为1.1，Maxsize为150点击okDelineatingtheRecharge 下一步是构建概念模型中的定义补给区域。假设除填埋场以外的地区补给都是均一的。在填埋场区域的补给因为衬垫而减小。CopyingtheBoundary 通过边界来创建补给图层。在 层中并在下拉菜单中选择Duplicate菜单命令在 层中并在下拉菜单中选择CoverageSetup菜单命令在ArealProperties列表中，打开Rechargerate选择okCreatingtheLandfill 构建填埋场边下一步，将创建弧段来描绘填埋场的边界。首先在填埋场大致区域创建一个方框。接下来边界节点和顶点让弧段于填埋场的边界相一致。选择CreateArc工 创建一个如图4-6所示的方框来描绘填埋场边界。不必担心点位的精（下面会进行精确定位选择SelectVertices工具在填埋场的矩形方框周围拉动一个方框，选择所有的顶点在其中一个顶点 ，并选择Vertex->Node菜单命令选择SelectPoints/Nodes工 节点被选择之后，在编辑窗口EditWindow.输入节点的精确坐标。输入每一个重复这个过程精确定 BuildingthePolygons 1.选择FeatureObjects|BuildPolygonsAssigningtheRecharge 补给区域已经构建，可以指派补给值。指派一个值给填埋场多边形单元，选择SelectPolygons工具将Rechargerate改为0.00006注意：补给率相对其他多边形单元要小很多。填埋场因为衬垫，压实等导致一个很小的补给值。此补给值本质上是填埋场很小的渗滤液渗漏。选择ok将Rechargerate改变为0.00695选择okDefiningtheHydraulic 数CopyingtheBoundary 通过边界来创建图层。在 层中并在下拉菜单中选择Duplicate菜单命令将新图层名字修改为Layer1在Layer 层中并在下拉菜单中选择CoverageSetup菜单命令在ArealProperties列表中，打开以下选项：HorizontalKVertical将Defaultlayerrange改变为从1到1（from1to选择ok在Layer1 层中并在下拉菜单中选择Duplicate菜单命令。将新图层名字修改为Layer2。在Layer2层 并在下拉菜单中选择CoverageSetup菜单命令将Defaultlayerrange改变为从2到2（from2to2选择okTopLayer 首先，为上部图层指派一个K值。在ProjectExplorer.中选择Layer1图 选择FeatureObjects|BuildPolygons选择SelectPolygons工具，并在多边形上双击将HorizonalK改变为5.5将Verticalanis.改变为4选择okBottomLayer 在ProjectExplorer.中选择Layer2图 选择FeatureObjects|BuildPolygons将HorizonalK改变为10将Verticalanis.改变为4选择ok1.在ProjectExplorer.中选择Sources&Sinks图 LocatingtheGrid 在ProjectExplorer中的空白区域，在下拉菜单中选择New|Grid在ProjectExplorer中 GridFrame，在下拉菜单中选择FittoActive在ProjectExplorer中栅格骨架中双击打开参数框将Originz和Dimensionz分别修改为170和60。这是为MODFLOW图层高程数列选择ok退出GridFrame框Creatingthe 选择FeatureObjects|Map→3DGrid5X，Y方向上的单元格数量不能改变。这是因在Z-Dimension中将Numcells改变为2选择okInitializingtheMODFLOWData MODFLOW数据。在ProjectExplorer.中选择3DGrid 选择MODFLOW|NewSimulation选择okDefiningtheActive/InactiveZones 在ProjectExplorer.中选择MapData选择SelectPolygons工具选择Properties按钮1－2，选择ok选择FeatureObjects|ActivateCellsinCoverage(s)InterpolatingLayer 插值（篡改）图层插值法在InterpolatingLayerDataImportingtheGroundSurfaceScatterPoints 通过在ProjectExplorer.中选择terrain让散点阵为激活状态。InterpolatingtheHeadsandElevations 下一步通过插值修改MODFLOW栅格的地表高程以及初始水头。在terrain散点阵中 并选择InterpolationTo|MODFLOWLayers菜单在InterpolatingLayerData指南中有详细的解释。突出显示ground_elev数据阵以及TopElevationsLayer1数列，并选择Map选择okInterpolatingtheLayer 在elevs散点阵上并选择InterpolationTo|MODFLOWLayers菜单命令GMS自动映射BottomElevationsLayer1和BottomElevationsLayer2数列基于其1.选择okAdjustingtheDisy 1.在ProjectExplorer.让散点阵处于未选择状态。1.在ProjectExplorer.ViewingtheModelCrossSections 在ProjectExplorer.中选择3DGridData选择SideView按钮为了更好的显示剖面图，需要增加Z选择 yOptions命令在Zmagnification中输入5选择ok 可以使用ToolPalette