A Hierarchical Sampling for Capturing Permeability Trend in Rock Physics

Among all properties of reservoir rocks, permeability is an important parameter which is typically measured from core samples in the laboratory. Due to limitations of core drilling all over a reservoir, simulation of rock porous media is demanded to explore more scenarios not seen in the available data. One of the most accurate methods is cross correlation based simulation (CCSIM) which recently has broadly applied in geoscience and porous media. The purpose of this study is producing realizations with the same permeability trend to a real sample. Berea sandstone sample is selected for this aim. Permeability results, extracted from smaller sub-samples of the original sample, showed that classic Kozeny–Carman permeability trend is not suitable for this sample. One reason can be due to lack of including geometrical and fractal properties of pore-space distribution in this equation. Thus, a general trend based on fractal dimensions of pore-space and tortuosity of the Berea sample is applied in this paper. Results show that direct 3D stochastic modeling of porous media preserves porous structure and fractal behavior of rock. On the other hand, using only 2D images for constructing the 3D pore structures does not reproduce the measured experimental permeability. For this aim, a hierarchical sampling is implemented in two and three steps using both 2D and 3D stochastic modeling. Results showed that two-step sampling is not suitable enough, while the utilized three-step sampling occurs to be show excellent performance by which different models of porous media with the same permeability trend as the Berea sandstone sample can be generated.     

Fractals — An overview of potential applications to transport in porous media

We present an overview of the potential applicability of fractal concepts to various aspects of transport phenomena in heterogeneous porous media. Three examples of phenomena where a fractal approach should prove illuminating are presented. In the first example we consider pore level heterogeneities as typified by pore surface roughness. We suggest that roughness may be usefully modelled by fractal curves and surfaces and also cite experimental evidence for regarding pores as fractals. In the second example we consider a fractal network approach to modelling large-scale heterogeneities. The presence of features on all length scales in simple fractal models should capture the essential role played by the presence of heterogeneities on many scales in natural reservoirs. Studies of transport phenomena in such models may yield valuable insights into the problems of macroscopic dispersion. The final example concerns dispersion in multiphase flow. Here the fractal character is attributed to the distribution of the fluid phases rather than the porous medium itself. Again studies of transport phenomena in simple fractal models should help to clarify various problems associated with the corresponding phenomena in real reservoirs.     

High-Resolution 3D FIB-SEM Image Analysis and Validation of Numerical Simulations of Nanometre-Scale Porous Ceramic with Comparisons to Experimental Results

The development of focused ion beam-scanning electron microscopy (FIB-SEM) techniques has allowed high-resolution 3D imaging of nanometre-scale porous materials. These systems are of important interest to the oil and gas sector, as well as for the safe long-term storage of carbon and nuclear waste. This work focuses on validating the accurate representation of sample pore space in FIB-SEM-reconstructed volumes and the predicted permeability of these systems from subsequent single-phase flow simulations using a highly homogeneous nanometre-scale, mesoporous (2–50 nm) to macroporous (>50 nm), porous ceramic in initial developments for digital rock physics. The limited volume of investigation available from FIB-SEM has precluded direct quantitative validation of petrophysical parameters estimated from such studies on rock samples due to sample heterogeneity, large variations in recorded sample pore sizes and lack of pore connectivity. By using homogeneous synthetic ceramic samples we have shown that lattice-Boltzmann flow simulations using processed FIB-SEM images are capable of predicting the permeability of a homogeneous material dominated by 10–100 nanometre-scale pores (similar, albeit simpler, to those in natural samples) at the much larger scale where permeability measurements become practical. This result shows the LB flow simulations can be used with confidence in pores at this scale allowing future work to focus on sample preparation techniques for samples sensitive to drying and multiple FIB-SEM site selection for the population of larger-scale models for heterogeneous systems.     

An Experimental Study on the Influence of Sub-Core Scale Heterogeneities on CO<Subscript>2</Subscript> Distribution in Reservoir Rocks

This article presents the results of CO₂/brine two-phase flow experiments in rocks at reservoir conditions. X-ray CT scanning is used to determine CO₂ saturation at a fine scale with a resolution of a few pore volumes and provide 3D porosity and saturation maps that can be use to correlate CO₂ saturations and rock properties. The study highlights the strong influence of sub-core scale heterogeneities on the spatial distribution of CO₂ at steady state and provides useful relative permeability data on a sample originated from an actual storage site (CO2CRC-Otway project, Victoria, South-West Australia). Two different samples tested, although different in nature, present strong heterogeneities, but differ in the detail of the connectivity of high porosity layers. In both samples, the results of the investigations show that sub-core scale heterogeneities control the sweep efficiency and may cause channeling through the porous medium. In one of the samples, CO₂ saturation appears uncorrelated to porosity close to the outlet end of the core. This observation is understood as a result of the position and the orientation of high porosity layers with respect to the inlet face of the core. Finally, in the operating conditions of the two experiments, the saturation maps demonstrate that gravity does not play a major role since no detectable buoyancy driven flow is observed.     

The effective homogeneous behavior of heterogeneous porous media

The macroscopic equations that govern the processes of one- and two-phase flow through heterogeneous porous media are derived by using the method of multiple scales. The resulting equations are mathematically similar to the point equations, with the fundamental difference that the local permeabilities are replaced by effective parameters. The method allows the determination of these parameters from a knowledge of the geometrical structure of the medium and its heterogeneities. The technique is applied to determine the effective parameters for one- and two-phase flows through heterogeneous porous media made up of two homogeneous porous media.     

An Analytical Model of Porosity–Permeability for Porous and Fractured Media

The classic Kozeny–Carman equation (KC) uses parameters that are empirically based or not readily measureable for predicting the permeability of unfractured porous media. Numerous published KC modifications share this disadvantage, which potentially limits the range of conditions under which the equations are applicable. It is not straightforward to formulate non-empirical general approaches due to the challenges of representing complex pore and fracture networks. Fractal-based expressions are increasingly popular in this regard, but have not yet been applied accurately and without empirical constants to estimating rock permeability. This study introduces a general non-empirical analytical KC-type expression for predicting matrix and fracture permeability during single-phase flow. It uses fractal methods to characterize geometric factors such as pore connectivity, non-uniform grain or crystal size distribution, pore arrangement, and fracture distribution in relation to pore distribution. Advances include (i) modification of the fractal approach used by Yu and coworkers for industrial applications to formulate KC-type expressions that are consistent with pore size observations on rocks. (ii) Consideration of cross-flow between pores that adhere to a fractal size distribution. (iii) Extension of the classic KC equation to fractured media absent empirical constants, a particular contribution of the study. Predictions based on the novel expression correspond well to measured matrix and fracture permeability data from natural sandstone and carbonate rocks, although the currently available dataset for fractures is sparse. The correspondence between model calculation results and matrix data is better than for existing models.     

Flow behaviour in the presence of wettability heterogeneities

The physical processes occurring during fluid flow and displacement within porous media having wettability heterogeneities have been investigated in specially designed heterogeneous visual models. The models were packed with glass beads, areas of which were treated with a water repellent to create wettability variations. Immiscible displacement experiments show visually the effect of wettability heterogeneities on the formation of residual oil and recovery due to capillary trapping. This work demonstrates by experiment the importance of incorporating reservoir heterogeneity into pore displacement analysis, essential for the correct interpretation of core data and for directing the route for scale-up of the processes to reservoir scale.     

THEORETICAL MODEL OF EFFECTIVE STRESS COEFFICIENT FOR ROCK/SOIL-LIKE POROUS MATERIALS

Physical mechanisms and influencing factors on the effective stress coefficient for rock/soil-like porous materials are investigated, based on which equivalent connectivity index is proposed. The equivalent connectivity index, relying on the meso-scale structure of porous material and the property of liquid, denotes the connectivity of pores in Representative Element Area （REA）. If the conductivity of the porous material is anisotropic, the equivalent connectivity index is a second order tensor. Based on the basic theories of continuous mechanics and tensor analysis, relationship between area porosity and volumetric porosity of porous materials is deduced. Then a generalized expression, describing the relation between effective stress coefficient tensor and equivalent connectivity tensor of pores, is proposed, and the expression can be applied to isotropic media and also to anisotropic materials. Furthermore, evolution of porosity and equivalent connectivity index of the pore are studied in the strain space, and the method to determine the corresponding functions in expressions above is proposed using genetic algorithm and genetic programming. Two applications show that the results obtained by the method in this paper perfectly agree with the test data. This paper provides an important theoretical support to the coupled hydro-mechanical research.     

The Effect of Rock Matrix Heterogeneities Near Fracture Walls on the Residence Time Distribution (RTD) of Solutes

For interpreting solute transport experiments within rocks, the Residence Time Distributions (RTDs) of solutes in the rock matrix were either derived from the assumption of a homogeneous matrix, or were considered through analytical distributions of diffusion coefficients. A numerical approach based on a Lagrangian framework was developed in order to investigate the effect of spatial heterogeneities on RTDs. The matrix diffusion was simulated over two-dimensional computation grids, representing virtual or real digitized porosity maps. First, virtual porosity maps were used to mimic porous features linked to conductive fractures, such as (i) low-porosity coatings on fracture walls, and (ii) porosity gradients within the rock matrix. Furthermore, RTD was calculated for the real pore network of the Palmottu granite (Finland). It was shown that the arrangement of spatial heterogeneities located in the immediate vicinity of a conductive fracture modifies the RTD of solutes within the rock matrix. Porous zones located near fractures are of particular importance, because they generate anomalies on the RTD.     

Impact of Rock Heterogeneity on Interactions of Microbial-Enhanced Oil Recovery Processes

Residual oil saturation reduction and microbial plugging are two crucial factors in microbial-enhanced oil recovery (MEOR) processes. In our previous study, the residual saturation was defined as a nonlinear function of the trapping number, and an explicit relation between the residual oil saturation and the trapping number was incorporated into a fully coupled biological (B) and hydrological (H) finite element model. In this study, the BH model is extended to consider the impact of rock heterogeneity on microbial-enhanced oil recovery phenomena. Numerical simulations of core flooding experiments are performed to demonstrate the influences of different parameters controlling the onset of oil mobilization. X-ray CT core scans are used to construct numerical porosity-permeability distributions for input to the simulations. Results show clear fine-scale fingering processing, and that trapping phenomena have significant effects on residual oil saturation and oil recovery in heterogeneous porous media. Water contents and bacterial distributions for heterogeneous porous media are compared with those for homogenous porous media. The evolution of the trapping number distribution is directly simulated and visualized. It is shown that the oil recovery efficiency of EOR/MEOR will be lower in heterogeneous media than in homogeneous media, largely due to the difficulty in supplying surfactant to unswept low-permeability zones. However, MEOR also provides efficient plugging along high-permeability zones which acts to increase sweep efficiency in heterogeneous media. Thus, MEOR may potentially be more suited for highly heterogeneous media than conventional EOR.     

Dependency of Tortuosity and Permeability of Porous Media on Directional Distribution of Pore Voids

Single-phase fluid flow in porous media is usually direction dependent owing to the tortuosity associated with the internal structures of materials that exhibit inherent anisotropy. This article presents an approach to determine the tortuosity and permeability of porous materials using a structural measure quantifying the anisotropic distribution of pore voids. The approach uses a volume averaging method through which the macroscopic tortuosity tensor is related to both the average porosity and the directional distribution of pore spaces. The permeability tensor is derived from the macroscopic momentum balance equation of fluid in a porous medium and expressed as a function of the tortuosity tensor and the internal structure of the material. The analytical results generally agree with experimental data in the literature.     

Global behavior of compressible three-phase flow in heterogeneous porous media

The homogenization method is used to analyze the equivalent behavior of a compressible three-phase flow model in heterogeneous porous media with periodic microstructure, including capillary effects. Asymptotic expansions lead to the definition of a global or effective model of an equivalent homogeneous reservoir. The resulting equations are of the same type as the points equations, with effective coefficients. The method allows the determination of these effective coefficients from a knowledge of the geometrical structure of the basic cell and its heterogeneities. Numerical computations to obtain the homogenized coefficients of the entire reservoir have been carried out via a finite element method.     

A new stochastic method of reconstructing porous media

We present a new stochastic method of reconstructing porous medium from limited morphological information obtained from two-dimensional micro- images of real porous medium. The method is similar to simulated annealing method in the capability of reconstructing both isotropic and anisotropic structures of multi-phase but differs from the latter in that voxels for exchange are not selected completely randomly as their neighborhood will also be checked and this new method is much simpler to implement and program. We applied it to reconstruct real sandstone utilizing morphological information contained in porosity, two-point probability function and linear-path function. Good agreement of those references verifies our developed method’s powerful capability. The existing isolated regions of both pore phase and matrix phase do quite minor harm to their good connectivity. The lattice Boltzmann method (LBM) is used to compute the permeability of the reconstructed system and the results show its good isotropy and conductivity. However, due to the disadvantage of this method that the connectivity of the reconstructed system’s pore space will decrease when porosity becomes small, we suggest the porosity of the system to be reconstructed be no less than 0.2 to ensure its connectivity and conductivity.     

Acoustical characterization of fluid-saturated porous media with local heterogeneities: Theory and application

A linear dynamic model of fully saturated porous media with local (either microscopic or mesoscopic) heterogeneities is developed within the context of Biot’s theory of poroelasticity. Viscoporoelastic behavior associated with local fluid flow is characterized by the notion of the dynamic compatibility condition on the interface between the solid and the fluid. Complex, frequency-dependent material parameters characterizing the viscoporoelasticity are derived. The complex properties can be obtained through determining the quasi-static poroelastic parameters, the properties of individual constituents, and the relaxation time of the dynamic compatibility condition on the interface. Relationships among various quasi-static poroelastic parameters are developed. It is shown that local fluid flow mechanism is significant only in the porous media with local heterogeneities. The relaxation time of the compatibility condition on the interface depends upon the details of local structure of porous media that control local fluid pressure diffusion. The new model is used to describe the velocity dispersion and attenuation in fully saturated porous media. The proposed model provides a theoretical framework to simulate the acoustical behavior of fully saturated porous media over a wide range of frequencies without making any explicit assumption about the structure of local heterogeneities.     

Time-Dependent Diffusion and Surface-Enhanced Relaxation in Stochastic Replicas of Porous Rock

Understanding the connection between pore structure and NMR behavior of fluid-saturated porous rock is essential in interpreting the results of NMR measurements in the field or laboratory and in establishing correlations between NMR parameters and petrophysical properties. In this paper we use random-walk simulation to study NMR relaxation and time-dependent diffusion in 3D stochastic replicas of real porous media. The microstructures are generated using low-order statistical information (porosity, void–void autocorrelation function) obtained from 2D images of thepore space. Pore size distributions obtained directly by a 3D pore space partitioning method and indirectly by inversion of NMR relaxation data are compared for the first time. For surface relaxation conditions typical of reservoir rock, diffusional coupling between pores of different size is observed to cause considerable deviations between the two distributions. Nevertheless, the pore space correlation length and the size of surface asperity are mirrored in the NMR relaxation data for the media studied. This observation is used to explain the performance of NMR-based permeability correlations. Additionally, the early time behavior of the time-dependent diffusion coefficient is shown to reflect the average pore surface-to-volume ratio. For sufficiently high values of the self-diffusion coefficient, the tortuosity of the pore space is also recovered from the long-time behavior of the time-dependent diffusion coefficient, even in the presence of surface relaxation. Finally, the simulations expose key limitations of the stochastic reconstruction method, and allow suggestions for future development to be made.     

Two-phase flow in heterogeneous porous media: The method of large-scale averaging

The analysis of two-phase flow in porous media begins with the Stokes equations and an appropriate set of boundary conditions. Local volume averaging can then be used to produce the well known extension of Darcy's law for two-phase flow. In addition, a method of closure exists that can be used to predict the individual permeability tensors for each phase. For a heterogeneous porous medium, the local volume average closure problem becomes exceedingly complex and an alternate theoretical resolution of the problem is necessary. This is provided by the method of large-scale averaging which is used to average the Darcy-scale equations over a region that is large compared to the length scale of the heterogeneities. In this paper we present the derivation of the large-scale averaged continuity and momentum equations, and we develop a method of closure that can be used to predict the large-scale permeability tensors and the large-scale capillary pressure. The closure problem is limited by the principle of local mechanical equilibrium. This means that the local fluid distribution is determined by capillary pressure-saturation relations and is not constrained by the solution of an evolutionary transport equation. Special attention is given to the fact that both fluids can be trapped in regions where the saturation is equal to the irreducible saturation, in addition to being trapped in regions where the saturation is greater than the irreducible saturation. Theoretical results are given for stratified porous media and a two-dimensional model for a heterogeneous porous medium.     

Pore-Scale Network Modeling of Three-Phase Flow Based on Thermodynamically Consistent Threshold Capillary Pressures. I. Cusp Formation and Collapse

We present a pore-scale network model of two- and three-phase flow in disordered porous media. The model reads three-dimensional pore networks representing the pore space in different porous materials. It simulates wide range of two- and three-phase pore-scale displacements in porous media with mixed-wet wettability. The networks are composed of pores and throats with circular and angular cross sections. The model allows the presence of multiple phases in each angular pore. It uses Helmholtz free energy balance and Mayer–Stowe–Princen (MSP) method to compute threshold capillary pressures for two- and three-phase displacements (fluid configuration changes) based on pore wettability, pore geometry, interfacial tension, and initial pore fluid occupancy. In particular, it generates thermodynamically consistent threshold capillary pressures for wetting and spreading fluid layers resulting from different displacement events. Threshold capillary pressure equations are presented for various possible fluid configuration changes. By solving the equations for the most favorable displacements, we show how threshold capillary pressures and final fluid configurations may vary with wettability, shape factor, and the maximum capillary pressure reached during preceding displacement processes. A new cusp pore fluid configuration is introduced to handle the connectivity of the intermediate wetting phase at low saturations and to improve model’s predictive capabilities. Based on energy balance and geometric equations, we show that, for instance, a gas-to-oil piston-like displacement in an angular pore can result in a pore fluid configuration with no oil, with oil layers, or with oil cusps. Oil layers can then collapse to form cusps. Cusps can shrink and disappear leaving no oil behind. Different displacement mechanisms for layer and cusp formation and collapse based on the MSP analysis are implemented in the model. We introduce four different layer collapse rules. A selected collapse rule may generate different corner configuration depending on fluid occupancies of the neighboring elements and capillary pressures. A new methodology based on the MSP method is introduced to handle newly created gas/water interfaces that eliminates inconsistencies in relation between capillary pressures and pore fluid occupancies. Minimization of Helmholtz free energy for each relevant displacement enables the model to accurately determine the most favorable displacement, and hence, improve its predictive capabilities for relative permeabilities, capillary pressures, and residual saturations. The results indicate that absence of oil cusps and the previously used geometric criterion for the collapse of oil layers could yield lower residual oil saturations than the experimentally measured values in two- and three-phase systems.     

Dependence of Pore-to-Core Up-scaled Reaction Rate on Flow Rate in Porous Media

Due to inherent heterogeneities in structure, mineral placement and fluid velocity in rock, bulk reaction rates realized during reactive flow through porous media may differ significantly from that predicted by laboratory-measured rate laws. In particular, rate laws determined in batch reactor experiments do not capture any of the flow dependence that will be experienced in the porous medium. Based on network flow model simulations of anorthite and kaolinite reactions in two sandstone pore networks under acidic conditions commensurate with CO2 sequestration, we compute up-scaled reaction rates at the core scale and investigate the dependence of the observed reaction rates on flow rate. For the anorthite reaction which, under these acidic conditions is far from equilibrium and dominated by pH, we find a power law dependence of reaction rate on flow rate. For the kaolinite reaction, which is near equilibrium, a more complex dependence emerges, with the up-scaled rate tending to rapidly increasing net precipitation at low-flow rates, then reversing and tending toward net dissolution at high-flow rates.