Simulation of multiphase flows in porous media

A 6M digital twin for modeling and simulation in subsurface reservoirs

The multiphase fluid flow in porous media is one of the most fundamental phenomena in various physical processes, such as oil/gas flow in reservoir, subsurface contamination dispersion, chemical separation, etc. Due to its importance, the efficient and accurate solution and prediction of multiphase flow in porous media is highly required in engineering applications and mechanism studies, which has been a research hot spot with increasing interest in recent years. However, the strong nonlinearity implicated in the multiphase flow model has brought great challenges for the computation and analysis. In addition, the permeability in Darcy-type pressure equation is always represented as a high contrast coefficient, which further complicates the simulation difficulty especially for high fidelity prediction and decision-making process. In this study, we attempt to apply one popular global model reduction method, namely Proper Orthogonal Decomposition (POD), to accelerate the simulation of incompressible multiphase flow in porous media. The cornerstone of POD is processing information from a sequence of snapshots (or instantaneous solutions) to identify a set of basis functions to construct a low-dimensional space, the most energetic part of these basis functions is then selected to represent the solution of the original system. In our work, one key point is how to establish the accurate POD reduced order model (ROM) for incompressible mul-tiphase porous flow problems. According to Li et al. [1], the traditional Galerkin projection method (directly project the original governing equation into the low-dimensional space spanned by basis functions) is incapable of introducing the fluid saturation and permeability located on the cell face in POD-ROM for incompressible multiphase flow whether in homogeneous or heterogeneous porous media, thus we adopt the matrix operation method (project the discrete governing equation instead of the original governing equation into the low-dimensional space spanned by basis functions) to build the POD-ROM in this study. Based on the proposed POD-ROM framework, we use the implicit method to solve the (non-) wetting-phase pressure and saturation equations. However, in the implicit solution of the pressure and saturation equations, the standard POD-ROM alone sometimes cannot yield an expected acceleration owing to two main reasons: (1) the nonlinearity between pressure and saturation is sensitive, thus the saturation would change largely even with slight pressure change. And the small error of pressure results would lead to a large deviation of saturation, which in turn exerts negative impacts on the computation of pressure; (2) the complexity for computing a projected nonlinear term in POD-ROM still depends on the dimension of the original full-order system, it is unavoidable in conventional POD-ROM. Therefore, more POD modes are needed to be used in order to achieve the prescribed numerical accuracy, which could worsen the speed up efficiency of proposed POD-ROM. To further reduce the computational complexity of POD-ROM to speed up the implicit solution of pressure and saturation equations, another model reduction method called Discrete Empirical Interpolation Method (DEIM) [2] is applied to treat the projected nonlinear term in the proposed incompressible multiphase porous flow POD-ROM. Different from previous studies, we develop an adaptive DEIM (ADEIM) in this work to achieve a more flexible selection of the interpolation points. Combined the POD-ROM with ADEIM, a POD-ADEIM-ROM for incompressible multiphase flow in porous media is established. The computational efficiency and numerical accuracy of the proposed model is validated through several representative numerical examples. Numerical results indicate that an outstanding acceleration (2~3 orders of magnitude faster) without sacrificing numerical accuracy significantly is obtained from the proposed model when comparing to the traditional solution method that without any acceleration technique. In addition, the effects of POD mode number, the distribution of permeability field, and the mesh density on the overall performances of the proposed model are investigated in detail.

Multiphase flow in porous media is involved in various natural and industrial applications, including water infiltration into soils, carbon geosequestration, and underground hydrogen storage. Understanding the invasion morphology at the pore scale is critical for better prediction of flow properties at the continuum scale in partially saturated permeable media. The deep learning method, as a promising technique to estimate the flow transport processes in porous media, has gained significant attention. However, existing works have mainly focused on single-phase flow, whereas the capability of data-driven techniques has yet to be applied to the pore-scale modeling of fluid–fluid displacement in porous media. Here, the conditional generative adversarial network is applied for pore-scale modeling of multiphase flow in two-dimensional porous media. The network is trained based on a data set of porous media generated using a particle-deposition method, with the corresponding invasion morphologies after the displacement processes calculated using a recently developed interface tracking algorithm. The results demonstrate the capability of data-driven techniques in predicting both fluid saturation and spatial distribution. It is also shown that the method can be generalized to estimate fluid distribution under different wetting conditions and particle shapes. This work represents the first effort at the application of the deep learning method for pore-scale modeling of immiscible fluid displacement and highlights the strength of data-driven techniques for surrogate modeling of multiphase flow in porous media.

The focus of the underlying research work is on the macroscopic modeling of unstable multiphase fluid flow in deformable porous media, where a lower‐viscous fluid is displaced by a more viscous fluid. This process leads to the formation of channel‐like networks, called viscous fingering. The instability effect is involved in a wide range of different fields in engineering. Some of the most common applications include carbon sequestration to store carbon dioxide (CO2) in underground reservoirs, contaminant transport in geostructures, in industrial processes, such as filtration, catalytic reactions, and in the operation of polymer electrolyte membrane fuel cells with multiphase flow in the gas diffusion layers. In this work, ideal miscible water–glycerin fluids are considered. It is assumed that the interacting fluids in the deformable porous media are incompressible. In addition, a dispersion–diffusion law is applied to capture the fluid–fluid interactions. The presence of a deformable porous material adds additional effects to the problem of the multiphase flow in porous media. The stresses in the porous solid matrix lead to changes in the porosity and influence the flow velocity of the fluids. To couple the deformable porous media with the multiphase flow, a macroscopic approach is used that relies on the theory of porous media. For the porous solid matrix, a linear elastic material model within small strains assumption is applied. The influence of the deformation‐dependent porosity on the instability is studied for 2D simulations, such as a multilayered geometry with different elastic parameters. The presented coupled nonlinear system of differential equations is simulated with the finite element method. Furthermore, a stabilization technique based on the quasi‐compressibility method is used.

A phase-field lattice Boltzmann model for simulating multiphase flows in porous media: Application and comparison to experiments of CO2 sequestration at pore scale

This article considered the numerical simulation of multicomponent multiphase flow in porous media. The resulting system of nonlinear equations linearized by the Newton-Raphson method and solved wi...

This study aims to bridge length scales in immiscible multiphase flow simulation by connecting two published governing equations at the pore-scale and continuum-scale through a novel validation framework. We employ Niessner and Hassnaizadeh's [“A model for two-phase flow in porous media including fluid-fluid interfacial area,” Water Resour. Res. 44(8), W08439 (2008)] continuum-scale model for multiphase flow in porous media, combined with the geometric equation of state of McClure et al. [“Modeling geometric state for fluids in porous media: Evolution of the Euler characteristic,” Transp. Porous Med. 133(2), 229–250 (2020)]. Pore-scale fluid configurations simulated with the lattice-Boltzmann method are used to validate the continuum-scale results. We propose a mapping from the continuum-scale to pore-scale utilizing a generalized additive model to predict non-wetting phase Euler characteristics during imbibition, effectively bridging the continuum-to-pore length scale gap. Continuum-scale simulated measures of specific interfacial area, saturation, and capillary pressure are directly compared to up-scaled pore-scale simulation results. This research develops a numerical framework capable of capturing multiscale flow equations establishing a connection between pore-scale and continuum-scale simulations.

Multiphase flow in porous media is relevant to amount of engineering processes, such as hydrocarbon extraction from reservoir rock, water contamination, CO2 geological storage and sequestration. Pore scale modeling, as an alternative approach to lab measurement, firstly serves as an effective bridge to link the pore scale properties (pore geometry and wettability) and displacement mechanisms to continuous scale multiphase flow in porous media; and secondly allows us to determine essential flow functions, such as capillary pressure and relative permeability curves, which are required for continuous scale modeling. In the literature, three methodologies, Bundle of Capillary Tube Modeling (BCTM), Direct Pore Scale Modeling (DPSM) and Pore Network Modeling (PNM), have appeared to be mostly widely adopted in the investigation of the pore-scale mechanics of fluid-fluid and fluid-solid interactions in porous media by numerical simulation. In this review article, a comprehensive review is provided to show their strengths and weaknesses and to highlight challenges that are faced in modelling of multiphase flow, key challenges include: are contact angle characterization, validation and upscale pore scale findings to core, or even field scale.

Numerical simulation of steam-chemical oil recovery processes require high temperature multiphase rock-fluid properties such as relative permeabilities and capillary numbers. These parameters are needed to predict and evaluate oil recovery performance and efficiency. Experiments in both Berea sandstone and sandpacks have been conducted to measure relative permeabilities of high and low tension oil-water mixtures at elevated temperatures. A new mathematical model consisting of sets of empirical correlations based on experimental data have been developed for relating endpoint relative permeabilities to both residual and mobile oil-water saturations and capillary numbers to residual oil-water saturations. The tests were run over a wide range of temperature, 72° to 400°F, and overburden pressure range of 500 to 1000 psia. It is our belief that this model is new for the range of temperature and pressure of steam oil recovery processes with chemical additives. Predictive results made with empirical correlations are in good agreement with experimental results.

This is the first year in the proposed three-year effort to develop high-resolution numerical methods for multi-phase flow in hierarchical porous media. The issues being addressed in this research are: Computational efficiency: Field-scale simulation of enhanced oil recovery, whether for energy production or aquifer remediation, is typically highly under-resolved. This is because rock transport properties vary on many scales, and because current numerical methods have low resolution. Effective media properties: Since porous media are formed through complex geologic processes, they involve significant uncertainty and scale-dependence. Given this uncertainty, knowledge of ensemble averages of flow in porous media can be preferable to knowledge of flow in specific realizations of the reservoir. However, current models of effective properties do not represent the observed behavior very well. Relative permeability models present a good example of this problem. In practice, these models seldom provide realistic representations of hysteresis, interfacial tension effects or three-phase flow; there are no models that represent well all three effects simultaneously. Wave propagation: It is common in the petroleum industry to assume that the models have the same well-posedness properties as the physical system. An example of this fallacy is given by the three-phase relative permeability models; they were widely assumed by the petroleum community to produce hyperbolic systems for the Buckley-Leverett equations, but later the mathematics community proved that these models inherently produce local elliptic regions. Since numerical methods must use the models for computations, oscillations that develop could erroneously be attributed to numerical error rather than modeling difficulties. During this year, we have made significant progress on several tasks aimed at addressing these issues.

Control volume function approximation methods and their applications to modeling porous media flow

Relative permeability is a critical parameter characterizing multiphase flow in porous media and it is strongly dependent on the wettability. In many situations, the porous media are nonuniformly wet. To investigate the effect of wettability heterogeneity on relative permeability of two‐phase flow in porous media, a multi‐relaxation‐time color‐gradient lattice Boltzmann model is adopted to simulate oil/water two‐phase flow in porous media with different oil‐wet solid fractions. For the water phase, when the water saturation is high, the relative permeability of water increases with the increase of oil‐wet solid fraction under a constant water saturation. However, as the water saturation decreases to an intermediate value (about 0.4–0.7), the relative permeability of water in fractionally wet porous media could be lower than that in purely water‐wet porous media, meaning additional flow resistance exists in the fractionally wet porous media. For the oil phase, similar phenomenon is observed. This phenomenon is mainly caused by the wettability‐related microscale fluid distribution. According to both our simulation results and theoretical analysis, it is found that the relative permeability of two‐phase flow in porous media is strongly related to three parameters: the fluid saturation, the specific interfacial length of fluid, and the fluid tortuosity in the flow direction. The relationship between the relative permeability and these parameters under different capillary numbers is explored in this paper.

Modeling multiphase flow in porous structures remains a challenge due to the complexity of handling multiple interfaces. This paper presents a one-domain pore-resolved simulation approach for immiscible two-phase flows in porous media, using a monolithic fluid–solid coupling framework to implicitly consider the existence of solid objects, with the fluid–fluid interfaces captured through solving an algebraic volume of fluid equation. Fluid interfacial tension is considered by integrating a continuum surface force, and the wall wettability condition is imposed by modifying the contact angle of the fluid interface at the embedded solid surface. The resulting equations are simple and stable, as there are no empirical models or parameters involved for the interface representation. This approach has been validated through performing a series of test-case simulations, including capillary-dominated flow, capillary rise with gravity, Taylor film formation, and finally two-phase flow in a heterogeneous porous structure. The numerical approach is demonstrated to be well suited for investigating pore-scale two-phase flows in realistic porous media.

Forced oil–water displacement is the crucial mechanisms of secondary oil recovery. The knowledge of relative permeability is required in the simulation of multiphase flow in porous media. Obvious dynamic effect of capillary pressure occurs in that the formation of ultra-low permeability reservoir (the permeability is <1 × 10−3 μm2) is tight and the pores and throats are very small. In addition, the significant capillary end effect causes serious errors when calculating relative permeabilities. For these reasons, the JBN method (neglecting capillary pressure) does not apply. Therefore, the dynamic capillary pressure and capillary end effects should be taken into account. This work focuses on calculating two-phase relative permeability of ultra-low permeability reservoir through considering the dynamic capillary pressure and eliminating the influence of capillary end effects. Firstly, laboratory core scale measurements of in situ water phase saturation history based on X-ray CT scanning technique were used to estimate relative permeability. Secondly, a mathematical model of two-phase relative permeability considering the dynamic capillary pressure was established. The basic problem formulations, as well as the more specific equations, were given, and the results of comparison using experimental data are presented and discussed. Results indicate that the dynamic capillary pressure measured at laboratory core scale in ultra-low permeability rocks has a significant influence on the estimation of unsteady-state relative permeability. The mathematical calculating method was compared with the history matching method and the results were close, suggesting reliability for ultra-low permeability reservoirs. Importantly, the proposed methods allow measurement of relative permeability from a single experiment. Potentially this represents a great time savings.

Mechanics and thermodynamics of multiphase flow in porous media including interphase boundaries