femwater记录.docVIP

第17章FEMWATERFEMWATER是模拟饱和—非饱和三维水流和溶质运移的有限元模型。17.1 问题概述如图17.1所示。三口开采井，开采量为100000 ft^3/day，上部含水层的K为10 ft/day，下部的为30 ft/day。17.3 所需模块Mesh模块（2D Mesh 和3D Mesh ），Map模块，Subsurface Characterization模块37GMS 使用说明刘延锋译（Borehole 、TIN 和Solid ），Geostatistics模块（2D Scatter Point 和3D Scatter Point ），FEMWATER界面。17.4 建立概念模型FEMWATER可以用概念模型的方法进行建模。17.4.1 导入背景图选择“File”菜单下的“Open”命令，打开“tutorial\femwater \ start.gpr”文件。17.4.2 初始化FENWATER图层转到Map模块，选择“Feature Objects”菜单下的“Coverages”命令，将Coverages对话框内的默认图层的名字设置为“femwater”，将图层类型Coverage type设置为“FEMWATER”。17.4.3 指定单位在Coverages对话框内，单击“Units”按钮，社设置Length和Time的单位（采用默认值），质量Mass单位为（slug）。17.4.4 创建边界弧段单击“Create Arc”工具按钮，17.4.5 修正弧段的顶点弧段的两端端点称为结点nodes，中间点称为顶点vertices。这些弧段用于生成可以转换成3D网格的2D网格。由顶点限定的线段间距控制了单元的大小和数目。因此，我们需要沿每一条弧段对定顶点进行调整，以保证间隔比较均匀而且长度合适。单击“Select Arcs”按钮，拖拉鼠标选择所有弧段，选择“Feature Objects”菜单下的“RedistributeVertices”命令，将“Target spacing”改为（300）。17.4.6 指定边界条件用弧段指定边界条件。单击“Select Arcs”按钮，按下SHIFT键，用鼠标单击选择定水头弧段，选择“Feature Objects”菜单下的“Attributes”命令，打开“Head/fluid flux”选项，OK。38GMS 使用说明刘延锋译此时并没有赋水头值，水头值赋在两端结点上，再通过线性插值获得。沿海弧段的水头值为0，可直接用默认值，河流下游点也为0，不用再赋。选择Select Points/Nodes工具，单击河流弧段的左上结点，输入水头值（200）。17.4.7 创建多边形需要多边形的原因：1）当产生网格时定义模型范围；2）用于赋补给值。选择“Feature Objects”菜单下的“Build Polygon”命令，OK。17.4.8 指定补给在FEMWATER中由两种方法：用指定流量边界或变化的边界条件。后者更精确，但更耗时，且不稳定。此处采用简单的指定流量法。选择“Select Polygons”按钮，在区内双击，打开“Fluid flux”选项，输入（0.003），OK。17.4.9 创建井单击Create Point工具按钮，在模型的右上部单击，在GMS顶部的坐标框内输入准确的坐标值（5290,4206, 45），选择“Feature Objects”菜单下的“Attributes”命令，打开“Refine mesh around point”选项（井附近的单元细化），为“Element size”输入值（150）（控制井处的单元大小）。然后将点指定为井点，指定开采量、滤管的顶底高程（流量赋予期间的单元）。打开“Well”选项，输入Top of screen（–185）和Bottom of screen（–195）及Flux rate（–100000），OK。同上，输入其他两个井：17.5 创建三维网格此时概念模型已经完成，要准备创建三维有限元网格。网格有两个区组成，一个区为上部含水层，另一个区为下部含水层。上区由两层单元组成；下区由三层单元组成。要创建网格，我们首先利用概念模型的特征对象创建二维投影网格。然后创建三个不规则三角形网格（TINs）：一个为顶（地面）面，一个为上部含水层的底面，另一个为下部含水层的底面。然后通过延伸二维网格的单元形成三维单元的单元柱，从而形成三维网格单元。17.5.1 指定材料创建网格之前，我们要为含水层指定材料。岩性指定给TINs，最后到3D网格。选择“Edit”菜单下的“Materials”命令，改变缺省材料名