Modeling Transport in Transient Ground-Water Flow: An Unacknowledged Approximation

Abstract. During unsteady or transient ground-water flow, the fluid mass per unit volume of aquifer changes as the potentiometric head changes, and solute transport is affected by this change in fluid storage. Three widely applied numerical models of two-dimensional transport partially account for the effects of transient flow by removing terms corresponding to the fluid continuity equation from the transport equation, resulting in a simpler governing equation. However, fluid-storage terms remaining in the transport equation that change during transient flow are, in certain cases, held constant in time in these models. For the case of increasing heads, this approximation, which is unacknowledged in these models'documentation, leads to transport velocities that are too high, and increased concentration at fluid and solute sources. If heads are dropping in time, computed transport velocities are too low. Using parameters that somewhat exaggerate the effects of this approximation, an example numerical simulation indicates solute travel time error of about 14 percent but only minor errors due to incorrect dilution volume. For horizontal flow and transport models that assume fluid density is constant, the product of porosity and aquifer thickness changes in time: initial porosity times initial thickness plus the change in head times the storage coefficient. This formula reduces to the saturated thickness in unconfined aquifers if porosity is assumed to be constant and equal to specific yield. The computational cost of this more accurate representation is insignificant and is easily incorporated in numerical models of solute transport.     

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 5. Single-Fluid-Phase Transport

This work is the fifth in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach for modeling flow and transport phenomena in multiscale porous medium systems. The general TCAT framework and the mathematical foundation presented in previous works are used to develop models that describe species transport and single-fluid-phase flow through a porous medium system in varying physical regimes. Classical irreversible thermodynamics formulations for species in fluids, solids, and interfaces are developed. Two different approaches are presented, one that makes use of a momentum equation for each entity along with constitutive relations for species diffusion and dispersion, and a second approach that makes use of a momentum equation for each species in an entity. The alternative models are developed by relying upon different approaches to constrain an entropy inequality using mass, momentum, and energy conservation equations. The resultant constrained entropy inequality is simplified and used to guide the development of closed models. Specific instances of dilute and non-dilute systems are examined and compared to alternative formulation approaches.     

利用三维非均匀不稳定渗流场的数值解求解垂直井与水平井的联合开采

用差分的方法给出了三维非均匀不稳定渗流方程的离散方程组，将配斯曼（Peaceman）方程拓展成能够被用于同时求解垂直井与水平井联合开采问题的形式，利用迭代解法求解差分方程和配斯曼方程所组成的方程组，并以实际层状不均匀介质为例计算了其压强的空间分布。     

Detection of an Impermeable Boundary in an Anisotropic Formation: A Case Study

Direct mirror-well image theory to locate an impermeable boundary does not apply in anisotropic formations because the relationship between flow lines and equipotential lines is distorted due to directional transmissivity. A simple approach enables an impermeable boundary to be located in anisotropic formations, and it can be verified with a field example.     

Two-Dimensional Flow and Mass Transport through Variably Saturated Porous Media from a Surface Waste Storage

Water Resources - Numerical calculations in DHI FEFLOW software environment were used to assess the effect of variably saturated zone parameters on capillary barrier formation. The sensitivity of...     

Flow Characteristics of Microbubble Suspensions in Porous Media as an Oxygen Carrier

Microbubble suspensions were generated as an oxygen carrier for aerobic biodegradation, and their flow characteristics in porous media were investigated. Commercial surfactants including sodium dodecyl sulfate (SDS), and dodecylethyldimethylammonium bromide (DEDAB), saponin (a natural surfactant), and collagen (a protein hydrolysate) were examined as base materials for generating microbubble suspensions. Among them, 2×CMC (critical micellar concentration) of SDS, DEDAB, and saponin developed microbubble suspensions with the highest gas hold‐up and half‐drainage time. Visualization of the flow patterns in sand showed that the microbubble suspensions were separated into a liquid and gas phase directly after injection, showing much faster movement of liquid phase flow. The gas front of the microbubble suspensions flowed in a plug‐flow manner, particularly in cases of SDS and DEDAB. The experimental results from both homogeneous and heterogeneous cells confirmed that the microbubble flow could overcome the heterogeneity in porous media. However, the plug‐flow characteristics and flow propagation of the microbubble suspensions to the low‐permeability zone was accompanied by a large pressure drop, which needs to be considered for future field application.     

Experiment and Simulation of Non-Reactive Solute Transport in Porous Media

To more accurately predict the migration behavior of pollutants in porous media, we conduct laboratory scale experiments and model simulation. Aniline (AN) is used in one-dimensional soil column experiments designed under various media and hydrodynamic conditions. The advection-dispersion equation (ADE) and the continuous-time random walk (CTRW) were used to simulate the breakthrough curves (BTCs) of the solute transport. The results show that the media and hydrodynamic conditions are two important factors affecting solute transport and are related to the degree of non-Fickian transport. The simulation results show that CTRW can more effectively describe the non-Fickian phenomenon in the solute transport process than ADE. The sensitive parameter in the CTRW simulation process is , which can reflect the degree of non-Fickian diffusion in the solute transport. Understanding the relationship of with velocity and media particle size is conducive to improving the reactive solute transport model. The results of this study provide a theoretical basis for better prediction of pollutant transport in groundwater.     

Energy Characteristics of the Transient Flow Effect on Channel and Bank Deformations and the Formation of Flow Secondary Pollution from Bottom Sediments

The ratio of the energy of water release into the lower pool of a waterworks facility for the time of release to the energy generated by the long-wave component of the wave spectrum caused by the release is considered. The effect of long-wave components of the wave spectrum on the channel and shore deformations is evaluated. An empirical relationship is proposed to describe the dependence of flow velocity on the release characteristics. Analysis of exchange processes between the flow and bottom deposits along with the obtained relationships is used to assess possible pollutant release from bottom deposits into water.     

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 8. Interface and Common Curve Dynamics

This work is the eighth in a series that develops the fundamental aspects of the thermodynamically constrained averaging theory (TCAT) that allows for a systematic increase in the scale at which multiphase transport phenomena is modeled in porous medium systems. In these systems, the explicit locations of interfaces between phases and common curves, where three or more interfaces meet, are not considered at scales above the microscale. Rather, the densities of these quantities arise as areas per volume or length per volume. Modeling of the dynamics of these measures is an important challenge for robust models of flow and transport phenomena in porous medium systems, as the extent of these regions can have important implications for mass, momentum, and energy transport between and among phases, and formulation of a capillary pressure relation with minimal hysteresis. These densities do not exist at the microscale, where the interfaces and common curves correspond to particular locations. Therefore, it is necessary for a well-developed macroscale theory to provide evolution equations that describe the dynamics of interface and common curve densities. Here we point out the challenges and pitfalls in producing such evolution equations, develop a set of such equations based on averaging theorems, and identify the terms that require particular attention in experimental and computational efforts to parameterize the equations. We use the evolution equations developed to specify a closed two-fluid-phase flow model.     

Simulation of Tidal Effects on Contaminant Transport in Porous Media

A one-dimensional numerical model is developed with oscillating velocities and dispersions to simulate the migration process of a contaminant plume within tidally influenced aquifers. Model simulations demonstrate that a major effect the tidal fluctuation has on the migration process of a contaminant plume is the exit concentration discharging to the tidal estuary. Tidal fluctuation causes the exit concentration levels to be significantly diluted by the surface-water body of the estuary. Sensitivity analyses demonstrate that tidal fluctuation hastens the rate of plume migration near the bank of the estuary because of the relatively high advective and dispersive fluxes induced by tides. However, tides affect the migration process only over a short distance from the tidal-water interface (about 40 ft for the parameters used in this study). If the contaminant plume is located far beyond the interface, tidal fluctuations will not affect the rate of plume migration until an existing regional ground-water flow velocity brings the plume to the tidally active zone. With or without tides, the rate of contaminant migration increases with higher regional hydraulic gradient. Furthermore, the effects of tidal fluctuations on the transport process become insignificant with higher regional hydraulic gradients.     

Hydraulic-Head Formulation for Density-Dependent Flow and Transport

Density-dependent flow and transport solutions for coastal saltwater intrusion investigations, analyses of fluid injection into deep brines, and studies of convective fingering and instabilities of denser fluids moving through less dense fluids typically formulate the groundwater flow equation in terms of pressure or equivalent freshwater head. A formulation of the flow equation in terms of hydraulic head is presented here as an alternative. The hydraulic-head formulation can facilitate adaptation of existing constant-density groundwater flow codes to include density-driven flow by avoiding the need to convert between freshwater head and hydraulic head within the code and by incorporating density-dependent terms as a compartmentalized “correction” to constant-density calculations already performed by the code. The hydraulic-head formulation also accommodates complexities such as unconfined groundwater flow and Newton-Raphson solution schemes more readily than the freshwater-head formulation. Simulation results are presented for four example problems solved using an implementation of the hydraulic-head formulation in MODFLOW.     

Analysis of Transient Flow to an Extended Fully Penetrating Well at Constant-Head Pumping

Vertical wells with radial extension at the well bottom can improve the rate of water production. No study has yet investigated the effects of the transient state and anisotropy in directional hydraulic conductivities on the wellbore flux rate for this type of well. This study derives a semianalytical transient drawdown solution for constant-head pumping at a fully penetrating well radially extended at the bottom of a confined, anisotropic aquifer by applying Laplace transform and separation of variables as well as conducting a Fourier analysis. The results of this new solution indicate that transient and steady-state wellbore flux rates can be increased by a factor of two for greater radial extension of the well. Compared with an isotropic aquifer (a ratio of vertical and horizontal hydraulic conductivities equal to one), an anisotropic aquifer with the ratio less than one may produce a higher transient wellbore flux rate and lower steady-state wellbore flux rate. Moreover, the time required to achieve the steady-state wellbore flux rate can be substantially affected by anisotropy of the aquifer.     

Transient and Transition Factors in Modeling Permafrost Thaw and Groundwater Flow

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near-surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three-dimensional model of a permafrost plateau that includes the supra-permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug-in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin-up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady-state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.     

Analysis of Flow and Sediment Transport Characteristics in Sand Beds

Analysis of flow velocities in stable beds is used to recognize the zones of eroding, transferring, and nontransferring velocities. The profiles of cross-sections in a pool and a rift are constructed using the author's formula.     

各向异性双相介质中多极源声波测井理论研究

将各向异性地层模拟为横向各向同性双相介质 ,使其对称主轴与井孔中心轴线相垂直 .本文针对井孔模型下的多极源声波测井问题提出了一种解析的摄动方法 .将各向异性双相介质看作是由对应的各向同性双相介质加上摄动部分而组成 ,并引入 5个关于弹性模量偏离各向同性时的摄动量 ,对对称主轴与井轴垂直的各向异性双相介质的多极源声波测井问题进行了系统的理论研究和分析 ,严格求解了井孔内外的零级和一级近似声波场 .发现在n级多极源激发下 ,井孔中不仅出现n级多极声场 ,而且还能激发出与各向同性地层时不相同的高于和低于n级的多极声波场     

Lattice Boltzmann Models for Flow and Transport in Saturated Karst

Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law–based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions.     

Studying the Flow Dynamics of a Karst Aquifer System with an Equivalent Porous Medium Model

The modeling of groundwater flow in karst aquifers is a challenge due to the extreme heterogeneity of its hydraulic parameters and the duality in their discharge behavior, that is, rapid response of highly conductive karst conduits and delayed drainage of the low‐permeability fractured matrix after recharge events. There are a number of different modeling approaches for the simulation of the karst groundwater dynamics, applicable to different aquifer as well as modeling problem types, ranging from continuum models to double continuum models to discrete and hybrid models. This study presents the application of an equivalent porous model approach (EPM, single continuum model) to construct a steady‐state numerical flow model for an important karst aquifer, that is, the Western Mountain Aquifer Basin (WMAB), shared by Israel and the West‐Bank, using MODFLOW2000. The WMAB was used as a catchment since it is a well‐constrained catchment with well‐defined recharge and discharge components and therefore allows a control on the modeling approach, a very rare opportunity for karst aquifer modeling. The model demonstrates the applicability of equivalent porous medium models for the simulation of karst systems, despite their large contrast in hydraulic conductivities. As long as the simulated saturated volume is large enough to average out the local influence of karst conduits and as long as transport velocities are not an issue, EPM models excellently simulate the observed head distribution. The model serves as a starting basis that will be used as a reference for developing a long‐term dynamic model for the WMAB, starting from the pre‐development period (i.e., 1940s) up to date.