Application of Digital Core Analysis (DCA) and Pore Network Modeling(PNM) based on 3D Micro-CT Images for an EOR project in a mature oil field in East Malaysia.

Characterizing internal microscopic structures of porous media is of vital importance to simulate fluid and electric current flow. Compared to traditional rock mechanics and geophysical experiments, digital core and pore network modeling is attracting more interests as it can provide more details on rock microstructure with much less time needed. The axis extraction algorithm, which has been widely applied for pore network modeling, mainly consists of a reduction and burning algorithm. However, the commonly used methods in an axis extraction algorithm have the disadvantages of complex judgment conditions and relatively low operating efficiency, thus losing the practicality in application to large-scale pore structure simulation. In this paper, the updated algorithm proposed by Palágyi and Kuba was used to perform digital core and pore network modeling. Firstly, digital core was reconstructed by using the Markov Chain Monte Carlo (MCMC) method based on the binary images of a rock cutting plane taken from heavy oil reservoir sandstone. The digital core accuracy was verified by comparing porosity and autocorrelation function. Then, we extracted the central axis of the digital pore space and characterize structural parameters through geometric transformation technology and maximal sphere method. The obtained geometric parameters were further assigned to the corresponding nodes of pore and throat on the central axis of the constructed model. Moreover, the accuracy of the new developed pore network model was measured by comparing pore/throat parameters, curves of mercury injection, and oil-water relative permeability. The modeling results showed that the new developed method is generally effective for digital core and pore network simulation. Meanwhile, the more homogeneity of the rock, which means the stronger “representative” of binary map the rock cutting plane, the more accurate simulated results can be obtained.

Download Full-text

Using Direct Numerical Simulation of Pore-Level Events to Improve Pore-Network Models for Prediction of Residual Trapping of CO2

Frontiers in Water ◽

10.3389/frwa.2021.710160 ◽

2022 ◽

Vol 3 ◽

Author(s):

Amir H. Kohanpur ◽

Yu Chen ◽

Albert J. Valocchi

Keyword(s):

Numerical Simulation ◽

Direct Numerical Simulation ◽

Network Modeling ◽

Pore Network ◽

Micro Ct ◽

Two Phase ◽

Modified Model ◽

Pore Network Modeling ◽

Modified Equations ◽

Pore Level

Direct numerical simulation and pore-network modeling are common approaches to study the physics of two-phase flow through natural rocks. For assessment of the long-term performance of geological sequestration of CO2, it is important to model the full drainage-imbibition cycle to provide an accurate estimate of the trapped CO2. While direct numerical simulation using pore geometry from micro-CT rock images accurately models two-phase flow physics, it is computationally prohibitive for large rock volumes. On the other hand, pore-network modeling on networks extracted from micro-CT rock images is computationally efficient but utilizes simplified physics in idealized geometric pore elements. This study uses the lattice-Boltzmann method for direct numerical simulation of CO2-brine flow in idealized pore elements to develop a new set of pore-level flow models for the pore-body filling and snap-off events in pore-network modeling of imbibition. Lattice-Boltzmann simulations are conducted on typical idealized pore-network configurations, and the interface evolution and local capillary pressure are evaluated to develop modified equations of local threshold capillary pressure of pore elements as a function of shape factor and other geometrical parameters. The modified equations are then incorporated into a quasi-static pore-network flow solver. The modified model is applied on extracted pore-network of sandstone samples, and saturation of residual trapped CO2 is computed for a drainage-imbibition cycle. The modified model yields different statistics of pore-level events compared with the original model; in particular, the occurrence of snap-off in pore-throats is reduced resulting in a more frontal displacement pattern along the main injection direction. Compared to the original model, the modified model is in closer agreement with the residual trapped CO2 obtained from core flow experiments and direct numerical simulation.

Download Full-text