Coupling upscaling finite element method for consolidation analysis of heterogeneous saturated porous media

The coupling upscaling finite element method is developed for solving the coupling problems of deformation and consolidation of heterogeneous saturated porous media under external loading conditions. The method couples two kinds of fully developed methodologies together, i.e., the numerical techniques developed for calculating the apparent and effective physical properties of the heterogeneous media and the upscaling techniques developed for simulating the fluid flow and mass transport properties in heterogeneous porous media. Equivalent permeability tensors and equivalent elastic modulus tensors are calculated for every coarse grid block in the coarse-scale model of the heterogeneous saturated porous media. Moreover, an oversampling technique is introduced to improve the calculation accuracy of the equivalent elastic modulus tensors. A numerical integration process is performed over the fine mesh within every coarse grid element to capture the small scale information induced by non-uniform scalar field properties such as density, compressibility, etc. Numerical experiments are carried out to examine the accuracy of the developed method. It shows that the numerical results obtained by the coupling upscaling finite element method on the coarse-scale models fit fairly well with the reference solutions obtained by traditional finite element method on the fine-scale models. Moreover, this method gets more accurate coarse-scale results than the previously developed coupling multiscale finite element method for solving this kind of coupling problems though it cannot recover the fine-scale solutions. At the same time, the method developed reduces dramatically the computing effort in both CPU time and memory for solving the transient problems, and therefore more large and computational-demanding coupling problems can be solved by computers.     

A new upscaling method for fractured porous media

We present a method to determine equivalent permeability of fractured porous media. Inspired by the previous flow-based upscaling methods, we use a multi-boundary integration approach to compute flow rates within fractures. We apply a recently developed multi-point flux approximation Finite Volume method for discrete fracture model simulation. The method is verified by upscaling an arbitrarily oriented fracture which is crossing a Cartesian grid. We demonstrate the method by applying it to a long fracture, a fracture network and the fracture network with different matrix permeabilities. The equivalent permeability tensors of a long fracture crossing Cartesian grids are symmetric, and have identical values. The application to the fracture network case with increasing matrix permeabilities shows that the matrix permeability influences more the diagonal terms of the equivalent permeability tensor than the off-diagonal terms, but the off-diagonal terms remain important to correctly assess the flow field.     

Time-of-flight (TOF)-based two-phase upscaling for subsurface flow and transport

Subsurface formations are characterized by heterogeneity over multiple length scales, which can have a strong impact on flow and transport. In this paper, we present a new upscaling approach, based on time-of-flight (TOF), to generate upscaled two-phase flow functions. The method focuses on more accurate representations of local saturation boundary conditions, which are found to have a dominant impact (in comparison to the pressure boundary conditions) on the upscaled two-phase flow models. The TOF-based upscaling approach effectively incorporates single-phase flow and transport information into local upscaling calculations, accounting for the global flow effects on saturation, as well as the local variations due to subgrid heterogeneity. The method can be categorized into quasi-global upscaling techniques, as the global single-phase flow and transport information is incorporated in the local boundary conditions. The TOF-based two-phase upscaling can be readily integrated into any existing local two-phase upscaling framework, thus more flexible than local–global two-phase upscaling approaches developed recently. The method was applied to permeability fields with different correlation lengths and various fluid-mobility ratios. It was shown that the new method consistently outperforms existing local two-phase upscaling techniques, including recently developed methods with improved local boundary conditions (such as effective flux boundary conditions), and provides accurate coarse-scale models for both flow and transport.     

Transport upscaling using multi-rate mass transfer in three-dimensional highly heterogeneous porous media

A methodology for transport upscaling of three-dimensional highly heterogeneous formations is developed and demonstrated. The overall approach requires a prior hydraulic conductivity upscaling using an interblock-centered full-tensor Laplacian-with-skin method followed by transport upscaling. The coarse scale transport equation includes a multi-rate mass transfer term to compensate for the loss of heterogeneity inherent to all upscaling processes. The upscaling procedures for flow and transport are described in detail and then applied to a three-dimensional highly heterogeneous synthetic example. The proposed approach not only reproduces flow and transport at the coarse scale, but it also reproduces the uncertainty associated with the predictions as measured by the ensemble variability of the breakthrough curves.     

Biofilm growth in porous media: Experiments,computational modeling at the porescale,and upscaling

Biofilm growth changes many physical properties of porous media such as porosity, permeability and mass transport parameters. The growth depends on various environmental conditions, and in particular, on flow rates. Modeling the evolution of such properties is difficult both at the porescale where the phase morphology can be distinguished, as well as during upscaling to the corescale effective properties. Experimental data on biofilm growth is also limited because its collection can interfere with the growth, while imaging itself presents challenges.In this paper we combine insight from imaging, experiments, and numerical simulations and visualization. The experimental dataset is based on glass beads domain inoculated by biomass which is subjected to various flow conditions promoting the growth of biomass and the appearance of a biofilm phase. The domain is imaged and the imaging data is used directly by a computational model for flow and transport. The results of the computational flow model are upscaled to produce conductivities which compare well with the experimentally obtained hydraulic properties of the medium. The flow model is also coupled to a newly developed biomass–nutrient growth model, and the model reproduces morphologies qualitatively similar to those observed in the experiment.     

Renormalization group analysis of permeability upscaling

The heterogeneity of the subsurface permeability is considered as the most influential factor in determining groundwater flow and the transport of toxic contaminants. Numerical simulators cannot handle the large grids required to represent the small-scale variability of permeability, and thus explicit estimates of the large-scale behavior in terms of coarse-grained parameters are often required. Perturbation formulations of the effective permeability are based on simplifying assumptions that are valid only for certain probability distributions and weak heterogeneity. A generalized perturbation ansatz that involves higher orders has been proposed (Gelhar and Axness, 1983), but to our knowledge its validity has not been rigorously proved before in three dimensions. In this work we propose a general upscaling formulation valid for strong heterogeneity, general permeability distributions, and media with impermeable zones. We show that the effective permeability is determined by the self-energy series of the permeability fluctuations at zero frequency. Using the diagrammatic representation, we obtain a Dyson equation that involves only irreducible diagrams of the proper self-energy series. We develop a renormalization group (RG) analysis for isotropic lognormal media that proves the generalized perturbation ansatz to all orders. We show that the RG result accurately estimates laboratory permeability measurements in limestone (strong heterogeneity) and sandstone (weak heterogeneity). We also propose an explicit RG estimate for the preasymptotic effective permeability. We compare our results with an approach based on a leading order Green's function expansion (Paleologos et?al., 1996), which, however, requires intensive numerical computations. Finally, we investigate the relation between the RG expression and the algebraic means used in numerical upscaling.     

渗流方程的三维粗化算法

本文提出了地下流体渗流问题的三维解粗化算法,在粗网格内流体压强分布用直接解法求解三维渗流方程,用这些解计算粗网格的等效渗透率,在流体流速大的区域仍采用精细网格的计算方法.用所得等效渗透率计算了粗化网格的渗流场的压强分布,结果表明渗流方程的三维粗化解非常逼近采用精细网格的解,但计算的速度比采用精细网格提高了100多倍.     

Direct simulations of two-phase flow on micro-CT images of porous media and upscaling of pore-scale forces

Pore-scale forces have a significant effect on the macroscopic behaviour of multiphase flow through porous media. This paper studies the effect of these forces using a new volume-of-fluid based finite volume method developed for simulating two-phase flow directly on micro-CT images of porous media. An analytical analysis of the relationship between the pore-scale forces and the Darcy-scale pressure drops is presented. We use this analysis to propose unambiguous definitions of Darcy-scale viscous pressure drops as the rate of energy dissipation per unit flow rate of each phase, and then use them to obtain the relative permeability curves. We show that this definition is consistent with conventional laboratory/field measurements by comparing our predictions with experimental relative permeability. We present single and two-phase flow simulations for primary oil injection followed by water injection on a sandpack and a Berea sandstone. The two-phase flow simulations are presented at different capillary numbers which cover the transition from capillary fingering at low capillary numbers to a more viscous fingering displacement pattern at higher capillary numbers, and the effect of capillary number on the relative permeability curves is investigated. Overall, this paper presents a new finite volume-based methodology for the detailed analysis of two-phase flow directly on micro-CT images of porous media and upscaling of the results to the Darcy scale.     

Automatic Generation of Locally Refined Composite Grids for Efficient Solute Transport Modeling

Widely used numerical models of solute transport processes in subsurface aquifers are limited to nonlocally refined rectangular, or logically rectangular, structured grids. This presents an unsuitable option to efficient numerical simulations maintaining an acceptable level of accuracy. Optimal selection of locally refined cells for efficient solute transport models is challenging to the current generation of numerical models. We present a novel and relatively simple to implement algorithm addressing these shortcomings. This method operates in four steps involving travel times simulations, a grid coarsening stage followed by a selective local grid refinement based on a cell-wise indicator, and a final postprocessing step. The refinement index is the sum of weighted logarithmic distributions of scaled forward and backward travel times. We calculate representative flow and transport properties at the two scales of the composite grid with a flow-based upscaling technique. We present two test problems to demonstrate the performances of this new gridding algorithm. We obtain the most important speedups for composite grids generated with the highest indicator thresholds. When hydrodynamic dispersion effects increase, we obtain less important speedups. An important outcome of this work is that grid design depends on nature and strength of the underlying flow and solute transport processes. Therefore, we suggest developing solute transport workflows integrating this grid generation algorithm as an integral component to build comprehensive and efficient groundwater models.     

Scaling of flow and transport behavior in heterogeneous groundwater systems

Three-dimensional numerical simulations using a detailed synthetic hydraulic conductivity field developed from geological considerations provide insight into the scaling of subsurface flow and transport processes. Flow and advective transport in the highly resolved heterogeneous field were modeled using massively parallel computers, providing a realistic baseline for evaluation of the impacts of parameter scaling. Upscaling of hydraulic conductivity was performed at a variety of scales using a flexible power law averaging technique. A series of tests were performed to determine the effects of varying the scaling exponent on a number of metrics of flow and transport behavior. Flow and transport simulation on high-performance computers and three-dimensional scientific visualization combine to form a powerful tool for gaining insight into the behavior of complex heterogeneous systems.Many quantitative groundwater models utilize upscaled hydraulic conductivity parameters, either implicitly or explicitly. These parameters are designed to reproduce the bulk flow characteristics at the grid or field scale while not requiring detailed quantification of local-scale conductivity variations. An example from applied groundwater modeling is the common practice of calibrating grid-scale model hydraulic conductivity or transmissivity parameters so as to approximate observed hydraulic head and boundary flux values. Such parameterizations, perhaps with a bulk dispersivity imposed, are then sometimes used to predict transport of reactive or non-reactive solutes. However, this work demonstrates that those parameters that lead to the best upscaling for hydraulic conductivity and head do not necessarily correspond to the best upscaling for prediction of a variety of transport behaviors. This result reflects the fact that transport is strongly impacted by the existence and connectedness of extreme-valued hydraulic conductivities, in contrast to bulk flow which depends more strongly on mean values. It provides motivation for continued research into upscaling methods for transport that directly address advection in heterogeneous porous media.An electronic version of this article is available online at the journal's homepage at http://www.elsevier.nl/locate/advwatres or http://www.elsevier.com/locate/advwatres (see “Special section on vizualization”. The online version contains additional supporting information, graphics, and a 3D animation of simulated particle movement.©1998 Elsevier Science Limited. All rights reserved     

Upscaling of hydraulic conductivity and telescopic mesh refinement

Performance assessments of repositories for the underground disposal of nuclear fuel and waste include models of ground water flow and transport in the host rocks. Estimates of hydraulic conductivity, K, based on field measurements may require adjustment (upscaling) for use in numerical models, but the choice of upscaling approach can be complicated by the use of nested modeling, large-scale fracture zones, and a high degree of heterogeneity. Four approaches to upscaling K are examined using a reference case based on exhaustive site data and an application of nested modeling to evaluate performance assessment of a waste repository. The upscaling approaches are evaluated for their effects on the flow balance between nested modeling domains and on simple measures of repository performance. Of the upscaling approaches examined in this study, the greatest consistency of boundary flows was achieved using the observed scale dependence for the rock domains, measured values from the large-scale interference test for the conductor domain, and a semivariogram regularization based on the Moye model for packer test interpretation. Making the assumption that large fracture zones are two-dimensional media results in the greatest changes to the median of travel time and improves the flow balance between the nested models. The uncertainty of upscaling methods apparently has a small impact on median performance measures, but a significant impact on the variances and earliest arrival times.     

Application of a data assimilation method via an ensemble Kalman filter to reactive urea hydrolysis transport modeling

With growing importance of water resources in the world, remediations of anthropogenic contaminations due to reactive solute transport become even more important. A good understanding of reactive rate parameters such as kinetic parameters is the key to accurately predicting reactive solute transport processes and designing corresponding remediation schemes. For modeling reactive solute transport, it is very difficult to estimate chemical reaction rate parameters due to complex processes of chemical reactions and limited available data. To find a method to get the reactive rate parameters for the reactive urea hydrolysis transport modeling and obtain more accurate prediction for the chemical concentrations, we developed a data assimilation method based on an ensemble Kalman filter (EnKF) method to calibrate reactive rate parameters for modeling urea hydrolysis transport in a synthetic one-dimensional column at laboratory scale and to update modeling prediction. We applied a constrained EnKF method to pose constraints to the updated reactive rate parameters and the predicted solute concentrations based on their physical meanings after the data assimilation calibration. From the study results we concluded that we could efficiently improve the chemical reactive rate parameters with the data assimilation method via the EnKF, and at the same time we could improve solute concentration prediction. The more data we assimilated, the more accurate the reactive rate parameters and concentration prediction. The filter divergence problem was also solved in this study.     

Another look at the conceptual fundamentals of porous media upscaling

We suggest a critical look at the epistemic foundations of the porous media upscaling problem that focuses on conceptual processes at work and not merely on form manipulations. We explore the way in which critical aspects of scientific methodology make their appearance in the upscaling context, thus generating useful effective parameters in practice. The fons et origo of our approach is a conceptual blending of knowledge states that requires the revision of the traditional method of scientific argument underlying most upscaling techniques. By contrast to previous techniques, the scientific reasoning of the proposed upscaling approach is based on a stochastic model that involves teleologic solutions and stochastic logic integration principles. The syllogistic form of the approach has important advantages over the traditional reasoning scheme of porous media upscaling, such as: it allows the rigorous derivation of the joint probability distributions of hydraulic gradients and conductivities across space; it imposes no restriction on the functional form of the effective parameters or the shape of the probability laws governing the random media (non-Gaussian distributions, multiple-point statistics and non-linear models are automatically incorporated); it relies on sound methodological principles rather than being ad hoc; and it offers the rational means for integrating the multifarious core knowledge bases and uncertain site-specific information sources about the subsurface system. Previous upscaling results are derived as special cases of the proposed upscaling approach under limited conditions of porous media flow, a fact that further demonstrates the generalization power of the approach. Our hope is that looking at the upscaling problem in this novel way will direct further attention to the methodological exploration of the problem at the length and the detail that it deserves.I would like to thank Drs. A. Kolovos and D.T. Hristopulos for their valuable comments. The work was supported by grants from the Army Research Office (Grant no. DAAG55–98–1-0289) and the National Institute of Environmental Health Sciences (P42-ES05948 & P30-ES10126).     

Uncertainty propagation of arbitrary probability density functions applied to upscaling of transmissivities

In many fields of study, and certainly in hydrogeology, uncertainty propagation is a recurring subject. Usually, parametrized probability density functions (PDFs) are used to represent data uncertainty, which limits their use to particular distributions. Often, this problem is solved by Monte Carlo simulation, with the disadvantage that one needs a large number of calculations to achieve reliable results. In this paper, a method is proposed based on a piecewise linear approximation of PDFs. The uncertainty propagation with these discretized PDFs is distribution independent. The method is applied to the upscaling of transmissivity data, and carried out in two steps: the vertical upscaling of conductivity values from borehole data to aquifer scale, and the spatial interpolation of the transmissivities. The results of this first step are complete PDFs of the transmissivities at borehole locations reflecting the uncertainties of the conductivities and the layer thicknesses. The second step results in a spatially distributed transmissivity field with a complete PDF at every grid cell. We argue that the proposed method is applicable to a wide range of uncertainty propagation problems.     

Real gas flow through heterogeneous porous media: theoretical aspects of upscaling

We consider the problem of upscaling transient real gas flow through heterogeneous bounded reservoirs. One of the commonly used methods for deriving effective permeabilities is based on stochastic averaging of nonlinear flow equations. Such an approach, however, would require rather restrictive assumptions about pressure-dependent coefficients. Instead, we use Kirchhoff transformation to linearize the governing stochastic equations prior to their averaging. The linearized problem is similar to that used in stochastic analysis of groundwater flow. We discuss the effects of temporal localization of the nonlocal averaged Darcy's law, as well as boundary effects, on the upscaled gas permeability. Extension of the results obtained by means of small perturbation analysis to highly heterogeneous porous formations is also discussed.     

Inverse upscaling of hydraulic parameters during constant flux infiltration using borehole radar

Modeling unsaturated flow in porous media requires constitutive relations that describe the soil water retention and soil hydraulic conductivity as a function of either potential or water content. Often, the hydraulic parameters that describe these relations are directly measured on small soil cores, and many cores are needed to upscale to the entire heterogeneous flow field. An alternative to the forward upscaling method using small samples are inverse upscaling methods that incorporate soft data from geophysical measurements observed directly on the larger flow field. In this paper, we demonstrate that the hydraulic parameters can be obtained from cross borehole ground penetrating radar by measuring the first arrival travel time of electromagnetic waves (represented by raypaths) from stationary antennae during a constant flux infiltration experiment. The formulation and coupling of the hydrological and geophysical models rely on a constant velocity wetting front that causes critical refraction at the edge of the front as it passes by the antennae. During this critical refraction period, the slope of the first arrival data can be used to calculate (1) the wetting velocity and (2) the hydraulic conductivity of the wet (or saturated) soil. If the soil is undersaturated during infiltration, then an estimate of the saturated water content is needed before calculating the saturated hydraulic conductivity. The hydraulic conductivity value is then used in a nonlinear global optimization scheme to estimate the remaining two parameters of a Broadbridge and White soil.     

Mixed multiscale finite elements and streamline methods for reservoir simulation of large geomodels

We present a new approach to reservoir simulation that gives accurate resolution of both large-scale and fine-scale flow patterns. The method uses a mixed multiscale finite-element method (MMsFEM) to solve the pressure equation on a coarse grid and a streamline-based technique to solve the fluid transport on a fine-scale subgrid. The MMsFEM is based on the construction of special approximation velocity spaces that are adaptive to the local properties of the differential operator. As such, MMsFEM produces a detailed subgrid velocity field that reflects the impact of the fine-scale heterogeneous structures. By combining MMsFEM with rapid streamline simulation of the fluid transport, we aim towards a numerical scheme that facilitates routine reservoir simulation of large heterogeneous geomodels without upscaling. The new method is applied to two different test cases. The first test case consists of two (strongly) heterogeneous quarter five-spot problems in 2D. The second test case is a 3D upscaling benchmark taken from the 10th SPE Comparative Solution Project, a project whose purpose is to compare and validate upscaling techniques. The test cases demonstrate that the combination of multiscale methods and streamlines is a robust and viable alternative to traditional upscaling-based reservoir simulation.     

Influence of permeability on hydrothermal circulation in the sediment-buried oceanic crust

1 Introduction Hydrothermal circulation is the key process of hydrothermal activity. Modern seafloor hydrothermal circulation can be divided into three parts: convective cells in the oceanic curst, interface between seafloor and ocean and hydrothermal plume. Hydrothermal convection in the crust is the dominant part of the whole seafloor hydrothermal circulation. The distribu-tion and nature of hydrothermal system in the oceanic crust are controlled by crust thermal structures and permeability …     

A coupled local–global upscaling approach for simulating flow in highly heterogeneous formations

A new technique for generating coarse scale models of highly heterogeneous subsurface formations is developed and applied. The method uses generic global coarse scale simulations to determine the boundary conditions for the local calculation of upscaled properties (permeability or transmissibility). An iteration procedure assures consistency between the local and global calculations. Transport processes are simulated using a subgrid velocity reconstruction technique applied in conjunction with the local–global upscaling procedure. For highly heterogeneous (e.g., channelized) systems, the new method is shown to provide considerably more accurate coarse scale models for flow and transport, relative to reference fine scale results, than do existing local (and extended local) upscaling techniques. The applicability of the upscaled models for different global boundary conditions is also considered.