MOVER - 三相有限元模型

MOVER是三相（水，石油和天然气）的有限元模型。MOVE可用于水、石油和天然气的流动模型，并通过不饱和区域中的NAPL截留来LNAPL和谁的回收率。MOVER的子模型可用于模拟水和LNAPL与静态大气气相的耦合流动。MOVER模拟非均质各向异性多孔介质或压裂介质。它允许使用等参元素来表示材料和物理/液压边界。MOVER可用于设计复杂水文地质条件下的NAPL回收好水力围堵系统。MOVER可用于模拟真空增强的恢复方案。

在常规的自由相回收系统（具有静态气相）中，当空气-油位下降时，油被捕集在不饱和中；如果油-水位升高，则油被捕集在饱和中。在计划不周的回收系统中，如果不自由项羽流，则由于凹陷锥内或超出影响半径的区域中的液位波动，可能会捕获70%的NAPL。真空增强的采纳率增加了水和油势的梯度，而流体表的波动较小。因此，真空增强的回收或生物制浆有助于较少残留产品的体积并提高自有产品的回收率，从而降低清洁成本。

关键功能

• 初始条件和游离油量是从监测井液位数据内部估算的

• 矩形或二维等参四边形元素，可以的对不规则的材料边界，水力和物理边界建模

• 模拟了NAPL和水相的真空增强回收率（生物浆化）

• 计算没有真空的油和水的回收率

• 剩余碳氢化合物的地域分布

• 交互式有限元网格生成器，矩形/等参四边形网格

• 水的补给，注入或LNAPL泄露在空间上是可变的

• 模拟多个抽水井和/或注入井

• 对指定的扬程和通量条件建模

• 模拟了破裂的介质或颗粒状多孔介质

• 模拟了NAPL和水相的生物浆化

移动输入

• 网格离散化数据

• 初始条件：水和油压力分布

• 流动的边界条件：指定的顶部边界和通量边界

• 源/汇边界：土壤水力van Genuchten参数，水力传导率分布和孔隙率

移动用户界面

Windows 3.x，Windows 95和Windows NT中的图形用户界面网格生成器/图形编辑器，它允许DXF站点地图导入，高度不规则的几何图形以及16位和32位版本。

移动输出

• 流体压力随时间的空间分布

• 流体饱和度随时间的空间分布

• 流体速度随时间的分布

• 流体泵送/诸如速率和体积与时间的关系

• MOVER的输出可以地用于使用BIOF& T 2-D /  3-D模拟多组分水相传输

MOVER 要求：（1500个节点）Windows 95 / 98 / 2000和Windows NT，带有16 MB RAM。

【英文介绍】

MOVER is an areal three-phase (water, oil and gas) finite-element model. Currently, the most advanced model of its kind, MOVER can be used to model flow of water, oil and gas, and optimize the recovery of LNAPL and water by minimizing NAPL entrapment in the saturated/unsaturated zones. A submodel of MOVER can be used to simulate coupled flow of water and LNAPL with a static atmospheric gas phase. MOVER simulates heterogeneous, anisotropic porous media or fractured media. It allows use of isoparametric elements to accurately represent material and physical/hydraulic boundaries. MOVER can be used to design NAPL recovery and hydraulic containment systems under complex hydrogeological conditions.

In a conventional free-phase recovery system (with static gas phase), oil is trapped in the unsaturated zone when the air-oil table falls and oil is trapped in the saturated zone if the oil-water table rises. In a poorly planned recovery system, 70 percent of the NAPL may be trapped due to fluid tables fluctuations within the cone of depression or in areas beyond the radius of influence if the free phase plume is not contained. Vacuum enhanced recovery increases gradients in water and oil potentials with minimal fluctuations in the fluid tables. Thus vacuum enhanced recovery or bioslurping helps to reduce volume of residual product and enhanced free product recovery, thus reducing cleanup costs.

MOVER KEY FEATURES

• Initial conditions and free oil volume are estimated internally from the monitoring well fluid level data.

• Rectangular or 2-D isoparametric quadrilateral elements to accurately model irregular material boundaries, hydraulic, and physical boundaries.

• Vacuum enhanced recovery (bioslurping) of NAPL and water phases are simulated.

• Oil and water recovery rates without vacuum are computed.

• Areal distribution of residual hydrocarbon.

• Interactive finite-element mesh generator, rectangular/isoparametric quadrilateral mesh.

• Spatially-variable water recharge, injection or LNAPL leakage.

• Model multiple pumping and/or injection wells.

• Model specified head and flux conditions.

• Simulates fractured media or granular porous media.

• Bioslurping of NAPL and water phases are simulated.

MOVER INPUT

• Mesh discretization data

• Initial conditions: water and oil pressure distribution

• Boundary conditions for flow: specified head boundaries and flux boundaries

• Source/sink boundary: Soil hydraulic properties include van Genuchten parameters, hydraulic conductivity distribution, and porosity

MOVER USER INTERFACE

Graphical user interface in Windows 3.x, Windows 95 and Windows NT includes a mesh generator/graphical editor which allows DXF site map import, highly irregular geometry and comes with a 16-bit and 32-bit version.

MOVER OUTPUT

• Spatial distribution of fluid pressure with time

• Spatial distribution of fluid saturation with time

• Fluid velocity distribution with time

• Fluid pumping/injection rates and volume vs. time

• Output from MOVER can be used seamlessly to simulate multicomponent aqueous phase transport using BIOF&T 2-D/3-D

MOVER TECHNICAL INFORMATION

MOVER (Multiphase Areal Organic Vacuum Enhanced Recovery Simulator) can be used to model recovery and migration of light nonaqueous phase liquids with vacuum enhanced recovery in unconfined heterogeneous, anisotropic aquifers. BIOF&T can be used in conjunction with MOVER to simulate multispecies dissolved phase transport in heterogeneous, anisotropic, fractured media, or unfractured granular porous media.

Groundwater contamination from hydrocarbon spills/leaks is a serious environmental problem. Nonaqueous phase liquids (NAPL) are immiscible fluids that have insignificant solubility in water. NAPLs in subsurface migrate under the influence of capillary, gravity, and buoyancy forces as a separate phase. Light NAPLs (LNAPLs) float and migrate on  of the water table posing a continuous source of contamination to the groundwater. Due to water table fluctuations, some of the NAPL gets trapped in the unsaturated and saturated zones. NAPL trapped in the soil and groundwater acts as a continuing source of groundwater contamination resulting in expensive restoration of these aquifers. MOVER optimizes recovery of LNAPL by vacuum enhanced recovery by increasing gradient in the oil potentials with minimum fluctuations in the fluid tables.

MOVER consists of:

• MOVER flow module simulates vacuum enhanced recovery and migration of water and LNAPL in unconfined aquifers following an LNAPL spill or leakage at a facility. It can also simulate NAPL recovery with skimmers and trenches with static gas phase. MOVER can be used to optimize the number, location, and recovery rates for water, oil, and soil vapor extraction wells.

• MOVER writes input files for BIOF&T, a transport model that simulates decoupled 2-D or 3-D multispecies aqueous phase transport from the free and residual NAPLs.

The MOVER flow module invokes an assumption of near-equilibrium conditions in the vertical direction. This reduces the nonlinearity in the constitutive model and transforms a 3-D problem into a 2-D areal problem, thereby drastically reducing computational time for the simulation.

MOVER gives the initial distribution of NAPL specific volume in the domain for BIOF&T which models the aqueous phase transport, computes and updates the temporal and spatial variation in the source during the simulation.

This software is accompanied by a user-friendly pre-processor, Mesh Editor and post-processor. The pre-processor and Mesh Editor can be used to create input data files for MOVER. The pre-processor and Mesh Editor include modules for: mesh generation; allocating heterogeneous and anisotropic soil properties; defining fixed head, flux, source/sink boundary conditions for water, oil, and gas phases; and allocating spatially-variable recharge in the domain. Two-dimensional rectangular or isoparametric quadrilateral elements are permissible to accurately model irregular domain and material boundaries.

Required input for flow analyses consists of initial air-oil table, air-water table distribution, soil hydraulic properties, fluid properties, time integration parameters, boundary conditions and mesh parameters. The van Genuchten constitutive model, along with fluid scaling parameters, is used to compute water and oil phase volumes.

MOVER output includes a list of the input parameters, initial and boundary conditions, and the mesh connectivity. It also includes simulated water, oil, and gas phase pressures, water, oil, and gas phase velocities at each node, total volume of water and oil versus time, and water and oil recovery/injection rates for each sink/source location versus time. Volume of free oil, and residual oil, and their spatial distributions are also printed versus time. Flow simulations can be performed in stages. MOVER creates an auxiliary file at the end of the current stage that can be used to define initial conditions for the next stage.

MOVER INPUT PARAMETERS

Estimation of Soil Properties

Soil properties needed for a MOVER flow simulation are: saturated hydraulic conductivity in principal flow directions, anisotropy angle of the main principal flow direction in the areal plane with the x-direction of the model domain, soil porosity , irreducible water saturation, and van Genuchten retention parameters. SOILPARA 1995, a proprietary computer model, provides an easy-to-use tool for estimating soil hydraulic parameters from soil texture based on: 1) the public domain model RETC developed by M. Th. van Genuchten et al., 1991, 2) the work of Shirazi and Boersma, 1984 and Campbell, 1985, and 3) a selection of USDA-recommended typical parameter values for various texture classes available in the SOILPARA database are included in the MOVER document.

Fluid Properties

Fluid properties required by MOVER are specific gravity, oil to water dynamic viscosity ratio, and fluid scaling parameters. Methods to estimate these parameters are included in the MOVER document.

Creating Input Data Files

The sequence of the input parameters and their definitions have been furnished in Appendix D of the document. This section explains the procedure for spatial discretization and mesh generation, defining initial conditions, boundary conditions, and the maximum permissible array dimensions.

Spatial Discretization and Mesh Generation

The MOVER modules allow use of rectangular and isoparametric elements. The element size and shape can be changed to obtain mesh refinement that are necessary to obtain accurate results.

Initial Conditions for Flow

Initial head distribution in the domain for water and oil can be specified by:

1) A bilinear interpolation with heads defined on the left and right boundaries

2) A non-uniform head distribution defined by fluid levels in the monitoring wells

Boundary Conditions

Specified pressure head (type-1) boundary conditions can be defined at selected nodes versus time.

Type-2 (specified flux) and source/sink boundary conditions can be defined by specifying the volumetric rate [L3 T -1] versus time for respective nodes. For a type-2 boundary condition, when flux [L T -1] is known at a node, the user should multiply flux with the area represented by the node in a plane perpendicular to the flux.

MOVER WINDOWS INTERFACE

What is the MOVER pre-processor?

The Windows pre-processor for MOVER is designed to help users create and edit input files for the MOVER numerical model. The pre-processor works in concert with the Mesh Editor to allow users to assign boundary condition schedules, soil types, recharge zones, etc., to the finite element mesh used in the MOVER numerical model. The pre-processor contains all control parameters that determine model run options, initial conditions, monitoring well information, fluid properties, boundary schedule data and soil type definitions as well as serving as a binder for Mesh Editor files. The pre-processor also contains a module for writing input files for the MOVER numerical model and for actually running the numerical model.

Using the MOVER pre-processor

The MOVER pre-processor runs under Windows 3.X, Windows 95 and Windows NT. The pre-processor uses the familiar tabbed notebook interface to allow quick editing of input files. The main program has two sets of tabs, one along the bottom which separates major sections of the interface, and, on some of the large notebook pages, tabs along the  that separate subsections to make the most use of available screen space. For example, clicking on the bottom tab "Boundary Schedules" takes the user to the boundary schedule notebook. Here there is a tabbed notebook for editing type 1 and type 2 boundary condition schedules.

The Tools selection on the main menu opens a tabbed notebook which includes Cue Cards, Files used in the pre-processor, a numerical model Runner, and a place to determine the location of the Mesh Editor. The files listed in the pre-processor are used to store variables and retrieve data and are generated automatically by the pre-processor and the Mesh Editor.

MOVER data files

All the MOVER variables for the numerical file are stored in ASCII text files that resemble Windows .ini files. These files are read and written by the MOVER pre-processor and Mesh Editor. They can be interchanged in the pre-processor setup window. For example, a material property file used in an earlier project can be assigned to a new project and all those soils will be available in the new project. A mesh file and all its associated files can be imported in the same manner. The data files can be edited with any ASCII text editor although this is not necessary. This open architecture was designed for future expansions of all DAEM models, or for third-party development of graphical interface tools.

What is the Mesh Editor?

The Mesh Editor was designed to work with these numerical models to create and edit finite element meshes. The Mesh Editor allows designing irregular quadrilateral meshes in two and three dimensions. Working with a numerical model pre-processor, the Mesh Editor provides a graphical interface for assigning properties to a mesh such as initial concentrations of contaminants, soil properties, boundary conditions, etc.

Speed Buttons

Pan

Panning allows the mesh to be moved in the Mesh Editor window. This puts the Mesh Editor into its "pan" state. You can then hold down the left mouse button and move the mesh on the screen. Letting go of the left mouse button "drops" the mesh in place. Nodes cannot be selected individually in this state.

Rotate

The rotate button is the leftmost red button on the tool bar. It puts the Mesh Editor into a "rotate" state where the mouse can be used to rotate the mesh in three dimensions. The mouse will rotate the mesh on the XYZ axis intersection. Rotating the Mesh Editor takes a little getting used to but if you ever get lost, the handy X-Y, Y-Z, and Z-X buttons will snap the mesh back into place.

Editing

The rightmost red button called Node Control puts the Mesh Editor in its "editing" state. Now a node or group of nodes can be selected to have values assigned to them, or to be moved in the X, Y and/or Z direction. To edit associations (i.e., soil type, recharge zones, type-1 boundary conditions, etc.) for a node, select the node, then right click the mouse. This will bring up a list of available associations. Clicking on an association will bring up a secondary list of associations that can be assigned to the node.

Moving Nodes

Nodes can be moved by holding down the Ctrl key (control) and the left mouse button, then moving the cursor on the screen. Nodes move according to the dimension displayed on the screen, so that two dimensional meshes should be in the default X-Y view for node movement. Nodes can only be moved when the Mesh Editor is in its "editing" state.

DXF Import

Version 1.1 of the DAEM Mesh Editor introduced DXF import. This tool allows for .dxf files to be placed on a mesh. This way, site files in CAD programs can be exported to the Mesh Editor, then used to aid mesh refinement and adjustment.

Post-processor

The post-processor is a data parsing tool, graphing package and contour export tool for these numerical models. The post-processor is designed to be a user-friendly tool for quickly discerning model results. Users of these models can also review model text output files for a more detailed view of model results.

MOVER VERIFICATION

Example: NAPL leak in an unconfined aquifer

This example is to test the accuracy of MOVER to simulate an LNAPL leak into an unconfined aquifer. The domain is 40 m x 40 m, discretized with uniformly spaced 21 rows and 21 columns (x = 2 m, y = 2 m). A leak occurred at the center of the domain (x = 20 m, y = 20 m) at a rate of 1 m 3/day for 20 days. Simulation was performed with an initial time step of 0.001 days, and the time incremental factor and maximum time step size were 1.03 and 0.2 days, respectively. The initial uniform piezometric head (Pao = Paw ) was 100 m throughout the domain. All boundaries were no-flow boundaries for water and oil phase. The water head at the out boundary was fixed at 100m throughout the simulation.

The simulated specific oil volume distribution is shown in the following figure along with the corresponding results from MOFAT. There are some differences in these solutions due to the difference in the formulations of the MOVER and MOFAT models. MOVER is a vertically-integrated areal model while MOFAT is a 2-D vertical slice (Planar or radially-symmetric vertical section) model. Nevertheless, it can be seen in the following figure that both solutions agree reasonably well over the range of the simulation.