Modeling of Two-Phase Flow in Fractured Porous Media on Unstructured Non-Uniformly Coarsened Grids

In this article, we investigate two strategies for coarsening fractured geological models. The first approach, which generates grids that resolve the fractures, is referred to as explicit fracture-matrix separation (EFMS). The second approach is based on a non-uniform coarsening strategy introduced in Aarnes et al. (Adv Water Resour 30(11):2177–2193, 2007a). A series of two-phase flow simulations where the saturation is modeled on the respective coarse grids are performed. The accuracy of the resulting solutions is examined, and the robustness of the two strategies is assessed with respect to number of fractures, degree of coarsening, well locations, phase viscosities, and fracture permeability. The numerical results show that saturation solutions obtained on the non-uniform coarse grids are consistently more accurate than the corresponding saturation solutions obtained on the EFMS grids. The numerical results also reveal that it is much easier to tune the upscaling factor with the non-uniform coarsening approach.     

A Least-Squares Mixed Scheme for the Simulation of Two-Phase Flow in Porous Media on Unstructured Grids

A least-squares mixed formulation is developed for simulation of two-phase flow in porous media. Such problems arise in petroleum applications and ground-water flow. An adaptive strategy based on the element residual as an error indicator is developed in conjunction with unstructured remeshing and tested for the two-phase flow of oil and water. An element-by-element conjugate-gradient scheme (EBE-CG) is compared to a band solution algorithm.     

Finite-Volume Discretisations for Flow in Fractured Porous Media

Over the last decades, finite-volume discretisations for flow in porous media have been extended to handle situations where fractures dominate the flow. Successful discretisations have been based on the discrete fracture-matrix models to yield mass conservative methods capable of explicitly incorporating the impact of fractures and their geometry. When combined with a hybrid-dimensional formulation, two central concerns are the restrictions arising from small cell sizes at fracture intersections and the coupling between fractures and matrix. Focusing on these aspects, we demonstrate how finite-volume methods can be efficiently extended to handle fractures, providing generalisations of previous work. We address the finite-volume methods applying a general hierarchical formulation, facilitating implementation with extensive code reuse and providing a natural framework for coupling of different subdomains. Furthermore, we demonstrate how a Schur complement technique may be used to obtain a robust and versatile method for fracture intersection cell elimination. We investigate the accuracy of the proposed elimination method through a series of numerical simulations in 3D and 2D. The simulations, performed on fractured domains containing permeability heterogeneity and anisotropy, also demonstrate the flexibility of the hierarchical framework.     

Effective Behavior Near Clogging in Upscaled Equations for Non-isothermal Reactive Porous Media Flow

For a non-isothermal reactive flow process, effective properties such as permeability and heat conductivity change as the underlying pore structure evolves. We investigate changes of the effective properties for a two-dimensional periodic porous medium as the grain geometry changes. We consider specific grain shapes and study the evolution by solving the cell problems numerically for an upscaled model derived in Bringedal et al. (Transp Porous Media 114(2):371–393, 2016. doi: 10.1007/s11242-015-0530-9). In particular, we focus on the limit behavior near clogging. The effective heat conductivities are compared to common porosity-weighted volume averaging approximations, and we find that geometric averages perform better than arithmetic and harmonic for isotropic media, while the optimal choice for anisotropic media depends on the degree and direction of the anisotropy. An approximate analytical expression is found to perform well for the isotropic effective heat conductivity. The permeability is compared to some commonly used approaches focusing on the limiting behavior near clogging, where a fitted power law is found to behave reasonably well. The resulting macroscale equations are tested on a case where the geochemical reactions cause pore clogging and a corresponding change in the flow and transport behavior at Darcy scale. As pores clog the flow paths shift away, while heat conduction increases in regions with lower porosity.     

Comparison of Three FE-FV Numerical Schemes for Single- and Two-Phase Flow Simulation of Fractured Porous Media

We benchmark a family of hybrid finite element–node-centered finite volume discretization methods (FEFV) for single- and two-phase flow/transport through porous media with discrete fracture representations. Special emphasis is placed on a new method we call DFEFVM in which the mesh is split along fracture–matrix interfaces so that discontinuities in concentration or saturation can evolve rather than being suppressed by nodal averaging of these variables. The main objective is to illustrate differences among three discretization schemes suitable for discrete fracture modeling: (a) FEFVM with volumetric finite elements for both fractures and porous rock matrix, (b) FEFVM with lower dimensional finite elements for fractures and volumetric finite elements for the matrix, and (c) DFEFVM with a mesh that is split along material discontinuities. Fracture discontinuities strongly influence single- and multi-phase fluid flow. Continuum methods, when used to model transport across such interfaces, smear out concentration/saturation. We show that the new DFEFVM addresses this problem producing significantly more accurate results. Sealed and open single fractures as well as a realistic fracture geometry are used to conduct tracer and water-flooding numerical experiments. The benchmarking results also reveal the limitations/mesh refinement requirements of FE node-centered FV hybrid methods. We show that the DFEFVM method produces more accurate results even for much coarser meshes.     

A Numerical Study of Micro-Heterogeneity Effects on Upscaled Properties of Two-Phase Flow in Porous Media

Commonly, capillary pressure–saturation–relative permeability (P ^c–S–K _r) relationships are obtained by means of laboratory experiments carried out on soil samples that are up to 10–12 cm long. In obtaining these relationships, it is implicitly assumed that the soil sample is homogeneous. However, it is well known that even at such scales, some micro-heterogeneities may exist. These heterogeneous regions will have distinct multiphase flow properties and will affect saturation and distribution of wetting and non-wetting phases within the soil sample. This, in turn, may affect the measured two-phase flow relationships. In the present work, numerical simulations have been carried out to investigate how the variations in nature, amount, and distribution of sub-sample scale heterogeneities affect P ^c–S–K _r relationships for dense non-aqueous phase liquid (DNAPL) and water flow. Fourteen combinations of sand types and heterogeneous patterns have been defined. These include binary combinations of coarse sand imbedded in fine sand and vice versa. The domains size is chosen so that it represents typical laboratory samples used in the measurements of P ^c–S–K _r curves. Upscaled drainage and imbibition P ^c–S–K _r relationships for various heterogeneity patterns have been obtained and compared in order to determine the relative significance of the heterogeneity patterns. Our results show that for micro-heterogeneities of the type shown here, the upscaled P ^c–S curve mainly follows the corresponding curve for the background sand. Only irreducible water saturation (in drainage) and residual DNAPL saturation (in imbibition) are affected by the presence and intensity of heterogeneities.     

裂缝介质渗流多尺度混合有限元数值模拟研究

针对裂缝介质具有多尺度特点,建立了Darcy/Stokes-Brinkman多尺度耦合模型,采用多尺度混合有限元方法,对裂缝介质渗流问题进行了研究.阐述了多尺度混合有限元方法的基本原理,并推导得到Darcy/Stokes-Brinkman方程的多尺度混合有限元计算格式.数值计算结果表明,大尺度Darcy模型能够捕捉到小尺度上裂缝网络渗流特征;与网格粗化、传统有限元方法相比,多尺度混合有限元方法的基函数具有能反映单元内参数变化的优点,在保证计算精度的同时能够减少计算量,对于裂缝油藏具有良好的适用性.     

Stokesian Dynamics Simulation of Suspension Flow in Porous Media

Transport in Porous Media - Suspension flow through porous medium was studied using the Stokesian dynamics simulation method. Stokesian dynamics is an efficient tool to carry out numerical...     

Modeling Fluid Flow in Fractured Porous Media with the Interfacial Conditions Between Porous Medium and Fracture

One of the most popular models that has been applied to predict the fluid velocity inside the fracture with impermeable walls is the cubic law. It highlights that the mean flux along the fracture is proportional to the cubic of fracture aperture. However, for a fractured porous medium, the normal and tangential interface conditions between the fracture and porous matrix can change the velocity profile inside the fracture. In this paper, a correction factor is introduced for flow equation along the fracture by imposing the continuity of normal and tangential components of velocity at the interface between the fracture and porous matrix. As a result, the mean velocity inside the fracture depends not only on the fracture aperture, but also on a set of non-dimensional numbers, including the matrix porosity, the ratio of intrinsic permeability of fracture to that of matrix, the wall Reynolds number, and the ratio of normal velocity on one wall to the other. Finally, the introduced correction factor is employed within the extended finite element method, which is widely used for numerical simulation of fluid flow within the fractured porous media. Several numerical results are presented for the fluid flow through a specimen containing single fracture, in order to investigate the deviation from the cubic law in different case studies.

     

High-Resolution Temporo-Ensemble PIV to Resolve Pore-Scale Flow in 3D-Printed Fractured Porous Media

Transport in Porous Media - Fractures are conduits that can enable fast advective transfer of (fluid, solute, reactant, particle, etc.) mass and energy. Such fast transfer can significantly affect...     

Flow in Fractured Porous Media: A Review of Conceptual Models and Discretization Approaches

Transport in Porous Media - The last decade has seen a strong increase of research into flows in fractured porous media, mainly related to subsurface processes but also in materials science and...     

Poroelastic Response of Orthotropic Fractured Porous Media

An algorithm is presented for inverting either laboratory or field poroelastic data for all the drained constants of an anisotropic (specifically orthotropic) fractured poroelastic system. While fractures normally weaken the system by increasing the mechanical compliance, any liquids present in these fractures are expected to increase the stiffness somewhat, thus negating to some extent the mechanical weakening influence of the fractures themselves. The analysis presented in this article quantifies these effects and shows that the key physical variable needed to account for the pore-fluid effects is a factor of (1 − B), where B is Skempton’s second coefficient and satisfies 0 ≤ B < 1. This scalar factor uniformly reduces the increase in compliance due to the presence of communicating fractures, thereby stiffening the fractured composite medium by a predictable amount. One further aim of the discussion is to determine the number of the poroelastic constants that needs to be known by other means to determine the rest from remote measurements, such as seismic wave propagation data in the field. Quantitative examples arising in the analysis show that, if the fracture aspect ratio a_f @ 0.1{a_f \simeq 0.1} and the pore fluid is liquid water, then for several cases considered, Skempton’s B @ 0.9{B \simeq 0.9}, and so the stiffening effect of the pore-liquid reduces the change in compliance due to the fractures by a factor 1 - B @ 0.1{1 - B \simeq 0.1}, in these examples. The results do, however, depend on the actual moduli of the unfractured elastic material, as well as on the pore-liquid bulk modulus, so these quantitative predictions are just examples, and should not be treated as universal results. Attention is also given to two previously unremarked poroelastic identities, both being useful variants of Gassmann’s equations for homogeneous—but anisotropic—poroelasticity. Relationships to Skempton’s analysis of saturated soils are also noted. The article concludes with a discussion of alternative methods of analyzing and quantifying fluid-substitution behavior in poroelastic systems, especially for those systems having heterogeneous constitution.     

Numerical Simulation of Two-Phase Inertial Flow in Heterogeneous Porous Media

In this study, non-Darcy inertial two-phase incompressible and non-stationary flow in heterogeneous porous media is analyzed using numerical simulations. For the purpose, a 3D numerical tool was fully developed using a finite volume formulation, although for clarity, results are presented in 1D and 2D configurations only. Since a formalized theoretical model confirmed by experimental data is still lacking, our study is based on the widely used generalized Darcy–Forchheimer model. First, a validation is performed by comparing numerical results of the saturation front kinetics with a semi-analytical solution inspired from the Buckley–Leverett model extended to take into account inertia. Second, we highlight the importance of inertial terms on the evolution of saturation fronts as a function of a suitable Reynolds number. Saturation fields are shown to have a structure markedly different from the classical case without inertia, especially for heterogeneous media, thereby, emphasizing the necessity of a more complete model than the classical generalized Darcy’s one when inertial effects are not negligible.     

Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs

A numerical method as well as a theoretical study of non-Darcy fluid flow through porous and fractured reservoirs is described. The non-Darcy behavior is handled in a three-dimensional, multiphase flow reservoir simulator, while the model formulation incorporates the Forchheimer equation for describing single-phase or multiphase non-Darcy flow and displacement. The non-Darcy flow through a fractured reservoir is handled using a general dual-continuum approach. The numerical scheme has been verified by comparing its results against those of analytical methods. Numerical solutions are used to obtain some insight into the physics of non-Darcy flow and displacement in reservoirs. In addition, several type curves are provided for well-test analyses of non-Darcy flow to demonstrate a methodology for modeling this type of flow in porous and fractured rocks, including flow in petroleum and geothermal reservoirs.     

Mesoscopic Simulation of Rarefied Flow in Narrow Channels and Porous Media

Rarefied gas flow in channels and computer-aided reconstructions of porous media is simulated using the direct simulation Monte Carlo (DSMC) method and a modified lattice-Boltzmann (LB) method that can account for rarefaction effects. An increase of two orders of magnitude was noted for the gas permeability as the Knudsen number increased from 0.1 to 10. It was found that incorporation of a Bosanquet expression for the viscosity in the dusty gas flux equations leads to the recovery of the well known Klinkenberg expression for the gas permeability, revealing an explicit relation of the, thus far empirical, permeability correction factor to the fluid and structure properties. The expression for the effective gas viscosity in the transition flow regime is also incorporated in the LB method, which is then validated against the DSMC method by comparing predictions for the velocity profiles in straight channels over the entire transition flow regime. The new LB method offers the additional advantages of simplicity in the code implementation and great savings in computational time and memory compared to the DSMC method. It is shown that the rough adjustment of a single parameter suffices to make the LB method suitable for the reliable prediction of the gas permeability in porous media over the entire transition flow regime.     

Flow in Random Porous Media

Flow in a porous medium with a random hydraulic conductivity tensor K(x) is analyzed when the mean conductivity tensor (x) is a non-constant function of position x. The results are a non-local expression for the mean flux vector (x) in terms of the gradient of the mean hydraulic head (x), an integrodifferential equation for (x), and expressions for the two point covariance functions of q(x) and (x). When K(x) is a Gaussian random function, the joint probability distribution of the functions q(x) and (x) is determined.     

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．