Multi-scale imaging process for computations of porosity and permeability from carbonate rocks. 9th Middle East Geosciences Conference, GEO 2010.

Abstract

Reservoir rock material collected during drilling is one of the main sources used to derive reservoir fluid transport and rock mechanics properties. Carbonate reservoirs may have heterogeneities that create multi porosity/permeability systems that are very difficult to describe, and to determine their flow properties. Conventional methods use laboratory procedures to perform experiments that directly or indirectly yield required rock properties. Some of these procedures, such as the determination of relative permeabilities, may take several months to perform. Also, in some cases, it is very difficult, or impractical to perform the experiments in the first place. Yet, as reservoir characterization is becoming ever more important for oil and gas production, a much larger portion of reservoir rocks, from cuttings to full cores, will need to be analyzed than what are currently evaluated. This study offers an example of the use of digital rock physics to determine porosity, permeability, and relative permeabilities for a carbonate sample using multi-scale imaging. Digital rock physics using the Lattice Boltzmann (LBM) for fluid dynamic calculations is at a point where for a proper digital pore space the resulting calculated flow properties are reasonably correct. The main issue facing digital rock physics is the need to up scale the computed properties to the scale of the core. The process presented in this study includes sample preparation, imaging, image processing, property computations, and property integration to the core scale. The sample is subjected to a descending scale of X-ray CT imaging, along with physical sub-sampling of the core. The descending size of scanning leads to increased resolution of the three-dimensional digital core, keeping the sample volumes registered in place. The resulting digital rocks are segmented and the pore structure is determined on the X-ray CT grid system. The resulting three-dimensional pore structure, that is the same as the actual pore structure subjected to resolution limits, is used as the input grid system for direct fluid dynamic computations that are second-order accurate representation of the Navier-Stokes fluid flow equations. These computations yield porosity, absolute permeability, relative permeabilities, and capillary pressure. In this study we focus only on porosity and permeabilities. Multiple scale imaging permits the estimation of permeability at the core scale.