Taking into account general boundary conditions in upscaling absolute permeability

Upscaling methods that need to solve local problems subject to boundary conditions are addressed in this article. We define a new upscaling method based on optimization problems, which can take into account general boundary conditions applied to local problems. The determination of upscaled permeability leads to minimizing the difference of dissipated energies (or averaged velocity) at fine and large scale. Using optimal control techniques, we obtain an effective computing algorithm that allows us to recover, with classical boundary conditions, the well-known results. The uniqueness issue is tackled for the optimization problems introduced in our approach. We show that the method is stable with respect to G-convergence, a property that establishes a link with homogenization theory, and finally, 2D numerical experiments are presented.     

Adaptive hierarchical upscaling of flow in heterogeneous reservoirs based on an a posteriori error estimate

This paper treats the upscaling of the absolute permeability in a heterogeneous reservoir. By replacing the fine scale permeability tensor with an upscaled, or effective permeability tensor, a modelling error is introduced. An a posteriori error estimate on this modelling error is formulated and tested. An implementation of the theory, based on domain decomposition coupled with a hierarchical representation of the absolute permeability field, is given. As hierarchical basis functions we have chosen the Haar system, which leads to a wavelet representation of the permeability. The wavelet representation offers a natural upscaling technique which resembles the highcut filters commonly used in signal analysis. This procedure represents an adaptive upscaling method. The numerical results show that this method conserves both the dissipation and the mean velocity in the problem fairly well. The a posteriori error estimate on the modelling error coupled with domain decomposition methods constitutes a powerful modelling tool.     

Unsteady source flow in weakly heterogeneous porous media

Average nonuniform flows in heterogeneous formations are modeled with the aid of the nonlocal effective Darcy's law. The mean head for flow toward source of instantaneous discharge in a heterogeneous medium of given statistics represents the fundamental solution of the average flow equation and is called the Mean Green Function (MGF). The general representation of the MGF is obtained for weakly heterogeneous formations as a functional of the logconductivity correlation function. For Gaussian logconductivity correlation, the MGF is derived in terms of one quadrature in time t and it is analyzed for isotropic media of any dimensionality d and for 3D axisymmetric formations. The MGF is further applied to determining the mean head distribution for flow driven by a continuous source of constant discharge. The large time asymptotic of the mean head is analyzed in details.     

多孔介质油水两相k~s~p关系数学模型的实验研究

数学模型是研究相对渗透率与饱和度关系曲线的重要方法。采用自行开发设计的人工平面多孔介质模型,测定了相对渗透率与饱和度的关系曲线。多孔介质选择粒径为0.5~1mm、1~2mm的标准砂,纯净的水为湿润相,用3号苏丹红染色的93#汽油为非湿润相,组成多孔介质油水两相流动系统。采用Van Genuchten and Mualeum(VGM)和Brooks-Corey-Burdine(BCB)两种数学模型计算相对渗透率与饱和度的关系曲线,通过比较两种数学模型计算结果之间和模型计算结果与实测结果的差异以及模型的应用、多相渗流系统自身特征,得出VGM、BCB两种数学模型计算结果符合实际情况,VGM模型应用过程更为简便,但VGM模型具有一定适用条件;在砂性多孔介质中,BCB模型计算相对渗透率与饱和度关系曲线更准确。     

Modelling of steady‐state fluid flow in 3D fractured isotropic porous media: application to effective permeability calculation

In this paper, 3D steady‐state fluid flow in a porous medium with a large number of intersecting fractures is derived numerically by using collocation method. Fluid flow in the matrix and fractures is described by Darcy's law and Poiseuille's law, respectively. The recent theoretical development presented a general potential solution to model the steady‐state flow in fractured porous media under a far‐field condition. This solution is a hypersingular integral equation with pressure field in the fracture surfaces as the main unknown. The numerical procedure can resolve the problem for any form of fractures and also takes into account the interactions and the intersection between fractures. Once the pressure field and then the flux field in fractures have been determined, the pressure field in the porous matrix is computed completely. The basic problem of a single fracture is investigated, and a semi‐analytical solution is presented. Using the solution obtained for a single fracture, Mori‐Tanaka and self‐consistent schemes are employed for upscaling the effective permeability of 3D fractured porous media. Copyright © 2012 John Wiley & Sons, Ltd.     

Calculating the internodal transmissibilities using finite analytic method and its application for multi‐phase flow in heterogeneous porous media

In the traditional numerical reservoir simulations, the internodal transmissibility is usually defined as the harmonic mean of the permeabilities of the adjacent grids. This definition underestimates the phase flux and the speed of the saturation front, especially for the strong heterogeneous case. In this article, the internodal transmissibility is recalculated according to the nodal analytic solution. The redefined internodal transmissibility can be used directly to calculate the multiphase flow in the numerical reservoir simulations. Numerical examples show that, compared to the traditional numerical methods, the proposed scheme makes the convergences much faster as the refinement parameter increases, and the accuracy is independent of the heterogeneity. Copyright © 2016 John Wiley & Sons, Ltd.     

Equivalent Permeability and Simulation of Two-Phase Flow in Heterogeneous Porous Media

The problem of calculating equivalent grid block permeability tensors for heterogeneous porous media is addressed. The homogenization method used involves solving Darcy's equation subject to linear boundary conditions with flux conservation in subregions of the reservoir and can be readily applied to unstructured grids. The resulting equivalent permeability tensor is stable as defined relative to G-convergence. It is proposed to use both conforming and mixed finite elements to solve the local problems and compute approximations from above and below of the equivalent permeability, respectively. Comparisons with results obtained using periodic, pressure and no-flux boundary conditions and the renormalization method are presented. A series of numerical examples demonstrates the effectiveness of the methodology for two-phase flow in heterogeneous reservoirs.     

Development of a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method

In the numerical modeling of fluid flow in heterogeneous geological media, large material contrasts associated with complexly intersected material interfaces are challenging, not only related to mesh discretization but also for the accurate realization of the corresponding boundary constraints. To address these challenges, we developed a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method (NMM) and the Lagrange multiplier method (LMM) for modeling boundary constraints. The advantages of NMM include meshing efficiency with fixed mathematical grids (covers), the convenience of increasing the approximation precision, and the high integration precision provided by simplex integration. In this discontinuous approach, the elements intersected by material interfaces are divided into different elements and linked together using the LMM. We derive and compare different forms of LMMs and arrive at a new LMM that is efficient in terms of not requiring additional Lagrange multiplier topology, yet stringently derived by physical principles, and accurate in numerical performance. To demonstrate the accuracy and efficiency of the NMM with the developed LMM for boundary constraints, we simulate a number of verification and demonstration examples, involving a Dirichlet boundary condition and dense and intersected material interfaces. Last, we applied the developed model for modeling fluid flow in heterogeneous media with several material zones containing a fault and an opening. We show that the developed discontinuous approach is very suitable for modeling fluid flow in strongly heterogeneous media with good accuracy for large material contrasts, complex Dirichlet boundary conditions, or complexly intersected material interfaces. Copyright © 2015 John Wiley & Sons, Ltd.     

A pore-scale numerical model for flow through porous media

A pore-scale numerical model based on Smoothed Particle Hydrodynamics (SPH) is described for modelling fluid flow phenomena in porous media. Originally developed for astrophysics applications, SPH is extended to model incompressible flows of low Reynolds number as encountered in groundwater flow systems. In this paper, an overview of SPH is provided and the required modifications for modelling flow through porous media are described, including treatment of viscosity, equation of state, and no-slip boundary conditions. The performance of the model is demonstrated for two-dimensional flow through idealized porous media composed of spatially periodic square and hexagonal arrays of cylinders. The results are in close agreement with solutions obtained using the finite element method and published solutions in the literature. Copyright © 1999 John Wiley & Sons, Ltd.     

A bridge between dual porosity and multiscale models of heterogeneous deformable porous media

In this paper, a multiscale homogenization approach is developed for fully coupled saturated porous media to represent the idealized sugar cube model, which is generally employed in fractured porous media on the basis of dual porosity models. In this manner, an extended version of the Hill-Mandel theory that incorporates the microdynamic effects into the multiscale analysis is presented, and the concept of the deformable dual porosity model is demonstrated. Numerical simulations are performed employing the multiscale analysis and dual porosity model, and the results are compared with the direct numerical simulation through 2 numerical examples. Finally, a combined multiscale-dual porosity technique is introduced by employing a bridge between these 2 techniques as an alternative approach that reduces the computational cost of numerical simulation in modeling of heterogeneous deformable porous media.     

A mixed solution for two-dimensional unsteady flow in fractured porous media

A mixed finite element–boundary element solution for the analysis of two-dimensional flow in porous media composed of rock blocks and discrete fractures is described. The rock blocks are modelled implicitly by using boundary elements whereas finite elements are adopted to model the discrete fractures. The computational procedure has been implemented in a hybrid code which has been validated first by comparing the numerical results with the closed-form solution for flow in a porous aquifer intercepted by a vertical fracture only. Then, a more complex problem has been solved where a pervious, homogeneous and isotropic matrix containing a net of fractures is considered. The results obtained are shown to describe satisfactorily the main features of the flow problem under study. © 1997 by John Wiley & Sons, Ltd.     

A finite element method for modeling coupled flow and deformation in porous fractured media

Modeling the flow in highly fractured porous media by finite element method (FEM) has met two difficulties: mesh generation for fractured domains and a rigorous formulation of the flow problem accounting for fracture/matrix, fracture/fracture, and fracture/boundary fluid mass exchanges. Based on the recent theoretical progress for mass balance conditions in multifractured porous bodies, the governing equations for coupled flow and deformation in these bodies are first established in this paper. A weak formulation for this problem is then established allowing to build a FEM. Taking benefit from recent development of mesh‐generating tools for fractured media, this weak formulation has been implemented in a numerical code and applied to some typical problems of hydromechanical coupling in fractured porous media. It is shown that in this way, the FEM that has proved its efficiency to model hydromechanical phenomena in porous media is extended with all its performances (calculation time, couplings, and nonlinearities) to fractured porous media. Copyright © 2015 John Wiley & Sons, Ltd.     

Dynamic capillary effects in heterogeneous porous media

In standard multi-phase flow models on porous media, a capillary pressure saturation relationship developed under static conditions is assumed. Recent experiments have shown that this static relationship cannot explain dynamic effects as seen for example in outflow experiments. In this paper, we use a static capillary pressure model and a dynamic capillary pressure model based on the concept of Hassanizadeh and Gray and examine the behavior with respect to material interfaces. We introduce a new numerical scheme for the one-dimensional case using a Lagrange multiplier approach and develop a suitable interface condition. The behavior at the interface is discussed and verified by various numerical simulations.     

Spatial averaging of hydraulic conductivity in three-dimensional heterogeneous porous media

Numerical models of groundwater flow require the assignment of hydraulic conductivities to large grid blocks discretizing the flow domain; however, conductivity data is usually available only at the much smaller scale of core samples. This paper describes a geostatistical model for hydraulic conductivity at both the core or point scale and that of grid blocks. Conductivity at the block scale is obtained empirically as a spatial power-average of point scale values. Assuming a multivariate Gaussian model for point log-conductivity, expressions are derived for the ensemble mean and variance of block conductivity. The expression for the ensemble mean of block scale conductivity is found to be similar to an expression for the ensemble effective conductivity of an infinite field derived analytically by earlier authors. Here, block conductivities obtained by power averaging are compared with effective conductivities obtained from a numerical flow model and are found to be in excellent agreement for a suitably chosen averaging exponent. This agreement deteriorates gradually as the log variance of conductivity increases beyond 2. For arbitrary flow field geometry and anisotropic conductivity covariances, the averaging exponent can be calibrated by recourse to numerical flow experiments. For cubic fields and an isotropic spatial covariance, the averaging exponent is found to be 1/3. In this particular case, it was found that flow field discretization at the block scale through local averaging of point conductivities gave similar results to those obtained directly using a point scale discretization of the flow field.     

非均匀多孔介质等效渗透率的普适表达式

多孔介质渗流是普遍的物理过程,涉及地下工程、地热开采、环境工程等各行各业,尤其是工程建设,常面临防渗问题。由于地质条件的复杂性,工程区域地层受到成岩、压实、风化、生物作用等各种影响,故渗流性质复杂,常需要对建设区域的渗流状况进行数值模拟,从而为工程的设计施工提供决策依据。数值仿真结果依赖于对地层介质关键参数的选取,但目前工程多将其视为均匀介质处理,对于介质的非均匀特性考虑较少。文章旨在研究非均质多孔介质渗透率空间分布与等效渗透率的关系。基于连续介质假定、达西定律以及非均匀多孔介质渗透率空间分布函数,建立一维到三维的达西渗流问题模型,通过求解偏微分方程和理论推导,得到基于渗透率空间分布函数的等效渗透率理论表达式,并与有限元计算的数值解进行对比分析,结果表明理论值和数值解误差很小,证明等效渗透率的表达式的合理性。利用该成果可通过多点局部渗透率的测定构建渗透率空间分布函数,从而对整体渗流区域的渗透性质进行快速计算和评估,从而简化异常复杂的工程地质模型以减少计算量需求,对于工程仿真的快速计算和结果评估有重要意义。     

Micromechanics of saturated and unsaturated porous media

The homogenization method is used to determine the formulation of the behaviour of both saturated and unsaturated porous media. This approach makes it possible to assess the validity of the effective stress concept as a function of the properties of the porous media at the microscopic scale. Furthermore, the influence of the morphologies of the solid and fluid phases on the macroscopic behaviour is studied. The strain induced by drying is examined as a function of the morphological properties. Copyright © 2002 John Wiley & Sons, Ltd.     

Multiscale finite-volume method for density-driven flow in porous media

The multiscale finite-volume (MSFV) method has been developed to solve multiphase flow problems on large and highly heterogeneous domains efficiently. It employs an auxiliary coarse grid, together with its dual, to define and solve a coarse-scale pressure problem. A set of basis functions, which are local solutions on dual cells, is used to interpolate the coarse-grid pressure and obtain an approximate fine-scale pressure distribution. However, if flow takes place in presence of gravity (or capillarity), the basis functions are not good interpolators. To treat this case correctly, a correction function is added to the basis function interpolated pressure. This function, which is similar to a supplementary basis function independent of the coarse-scale pressure, allows for a very accurate fine-scale approximation. In the coarse-scale pressure equation, it appears as an additional source term and can be regarded as a local correction to the coarse-scale operator: It modifies the fluxes across the coarse-cell interfaces defined by the basis functions. Given the closure assumption that localizes the pressure problem in a dual cell, the derivation of the local problem that defines the correction function is exact, and no additional hypothesis is needed. Therefore, as in the original MSFV method, the only closure approximation is the localization assumption. The numerical experiments performed for density-driven flow problems (counter-current flow and lock exchange) demonstrate excellent agreement between the MSFV solutions and the corresponding fine-scale reference solutions.     

模拟裂隙多孔介质中变饱和渗流的广义等效连续体方法

描述了一种计算裂隙多孔介质中变饱和渗流的广义等效连续体方法。这种方法忽略裂隙的毛细作用,设定一个与某孔隙饱和度相对应的综合饱和度极限值,并假定：（1）如果裂隙多孔介质的综合饱和度小于该极限值,水只在孔隙中存在并流动,而裂隙中则没有水的流动;（2）如果综合饱和度等于或大于该极限值,水将进入裂隙,并在裂隙内运动。分析比较了等效连续体模型的不同计算方法,并给出了一个模拟裂隙岩体中变饱和渗流与传热耦合问题的应用算例。结果表明,所述方法具有一般性,可以有效地模拟裂隙多孔介质中变饱和渗流的基本特征。     

多孔介质渗透率的分形描述

针对土壤、岩石等多孔介质结构的复杂性,从其结构形成的物理机制和达西定律出发,利用分形几何理论,将土壤等作为在统计意义上具有分形特征的多孔介质来研究其水力参数与结构之间的关系,建立了饱和多孔介质渗透率与其分维数之间的定量化的函数式。试验应用扫描电镜法研究了多孔介质断面微结构并算出分维数。试验结果表明:利用该模型预测的多孔介质渗透率与实测值基本吻合,能够比较精确地预测多孔介质水力参数。     

多孔介质中含均相/非均相反应传质过程数值模拟

结合化学反应方程式,并应用多孔多相介质溶混污染物输运过程的数值模型,对多孔多相介质中含均相/非均相化学反应传质过程进行了数值模拟。化学反应主要包含均相快速/慢速和非均相快速/慢速等5种化学反应过程,溶质输运行为的控制机制主要考虑对流、扩散及降解、吸附等。基于原有的隐式特征线Galerkin离散化的有限元方法,求解模型控制方程的边值初值问题,求解过程中把均相化学反应物质中按照反应物和生成物分开,非均相反应物质按照固相和液相分开,对均相反应物及非均相液相物质浓度耦合求解,而均相生成物和非均相固相物质独立求解。使方程组按照其不同类型进行分类,同时可减少未知数的个数。对于含有非线性内状态变量的右端项进行迭代求解。数值例题结果验证了所提出的数值方法的有效性、计算精度和稳定性。