Solution Construction to a Class of Riemann Problems of Multiphase Flow in Porous Media

Transport in Porous Media - Fluid displacement in porous media can usually be formulated as a Riemann problem. Finding the solution to such a problem helps shed light on the dynamics of flow and...     

A Computationally Efficient Approach to Applying the SAFT Equation for CO2 + H2O Phase Equilibrium Calculations

The statistical associating fluid theory equation of state (EoS) is employed in a time efficient way for the correlation and prediction of vapor–liquid equilibrium of the CO₂ + H₂O binary system for the temperature (10–100 °C) and pressure (1–600 bar) ranges suitable for simulation of CO₂ geologic sequestration. The effective number of segments and energy parameter are correlated with the reduced temperature. Simple mixing rules are applied to obtain binary interaction parameters. Assigning a fixed H₂O composition in the mixing rule makes the phase equilibrium calculations relatively fast compared to other EoS’s. The results obtained by the model used were found to be in satisfactory agreement with the literature data.     

Displacement Theory of Low-Tension Gas Flooding

Low-tension gas (LTG) flooding is a promising chemical enhanced oil recovery technique in tight sandstone and carbonate reservoirs where polymer may not be used because of plugging and degradation issues. This process has been the subject of many experimental studies. However, theoretical investigation of the LTG process is scarce in the literature. Hence, in this study, we lay out a displacement theory for LTG flooding, with a constant mobility reduction factor, which lays the groundwork for further theoretical studies. The proposed model is based on the three-phase flow of water, oil, and gas in the presence of a water-soluble surfactant component. Under the developed model, we study the effect of MRF and oil viscosity on the flow dynamics and oil recovery. Moreover, we explain experimental observations on early gas breakthrough that occurs during LTG core floods even in the presence of a stable foam drive.

     

New Trapping Mechanism in Carbon Sequestration

The modes of geologic storage of CO₂ are usually categorized as structural, dissolution, residual, and mineral trapping. Here we argue that the heterogeneity intrinsic to sedimentary rocks gives rise to a fifth category of storage, which we call local capillary trapping. Local capillary trapping occurs during buoyancy-driven migration of bulk phase CO₂ within a saline aquifer. When the rising CO₂ plume encounters a region (10⁻² to 10⁺¹m) where capillary entry pressure is locally larger than average, CO₂ accumulates beneath the region. This form of storage differs from structural trapping in that much of the accumulated saturation will not escape, should the integrity of the seal overlying the aquifer be compromised. Local capillary trapping differs from residual trapping in that the accumulated saturation can be much larger than the residual saturation for the rock. We examine local capillary trapping in a series of numerical simulations. The essential feature is that the drainage curves (capillary pressure versus saturation for CO₂ displacing brine) are required to be consistent with permeabilities in a heterogeneous domain. In this work, we accomplish this with the Leverett J-function, so that each grid block has its own drainage curve, scaled from a reference curve to the permeability and porosity in that block. We find that capillary heterogeneity controls the path taken by rising CO₂. The displacement front is much more ramified than in a homogeneous domain, or in a heterogeneous domain with a single drainage curve. Consequently, residual trapping is overestimated in simulations that ignore capillary heterogeneity. In the cases studied here, the reduction in residual trapping is compensated by local capillary trapping, which yields larger saturations held in a smaller volume of pore space. Moreover, the amount of CO₂ phase remaining mobile after a leak develops in the caprock is smaller. Therefore, the extent of immobilization in a heterogeneous formation exceeds that reported in previous studies of buoyancy-driven plume movement.     

Application of higher-order flux-limited methods in compositional simulation

A higher-order flux-limited finite-difference scheme has been implemented into a compositional simulator to discretize the convection terms of the component conservation equations and the relative permeability terms of the phase fluxes. Harten's total variation diminishing criteria are imposed directly to the finite-difference equations and the bounds of the flux limiters which are suitable for larger Courant numbers and nonuniform grid systems are obtained. A time-correction technique is applied to increase the time accuracy and improve the stability condition.The scheme has been tested for miscible and immiscible flow problems in one and two dimensions, and the results were compared with those using a third-order method without flux limiting and some available analytical solutions.     

Finite-Difference Approximation for Fluid-Flow Simulation and Calculation of Permeability in Porous Media

We introduce a finite-difference method to simulate pore scale steady-state creeping fluid flow in porous media. First, a geometrical approximation is invoked to describe the interstitial space of grid-based images of porous media. Subsequently, a generalized Laplace equation is derived and solved to calculate fluid pressure and velocity distributions in the interstitial space domain. We use a previously validated lattice-Boltzmann method (LBM) as ground truth for modeling comparison purposes. Our method requires on average 17 % of the CPU time used by LBM to calculate permeability in the same pore-scale distributions. After grid refinement, calculations of permeability performed from velocity distributions converge with both methods, and our modeling results differ within 6 % from those yielded by LBM. However, without grid refinement, permeability calculations differ within 20 % from those yielded by LBM for the case of high-porosity rocks and by as much as 100 % in low-porosity and highly tortuous porous media. We confirm that grid refinement is essential to secure reliable results when modeling fluid flow in porous media. Without grid refinement, permeability results obtained with our modeling method are closer to converged results than those yielded by LBM in low-porosity and highly tortuous media. However, the accuracy of the presented model decreases in pores with elongated cross sections.