Cinar, Yildiray (The University of New South Wales) | Arns, Christoph (The University of New South Wales) | Dehghan Khalili, Ahmad (The University of New South Wales) | Yanici, Sefer (The University of New South Wales)
Resistivity measurements play a key role in hydrocarbon in place calculations for oil and gas reservoirs. They are a direct indi-cator of fluid saturation and connected pore space available in the formation. Carbonate rocks, which host around half of the world's hydrocarbons, exhibit a wide range of porosities with scales spanning from nanometres to centimetres. The often sig-nificant amount of microporosity displayed by Carbonate rocks emphasizes the necessity of an adequate characterization of their micro-features and their contribution to hydrocarbon in place. In this paper we examine upscaling methods to probe for-mation factor of a fully saturated carbonate sample using an X-ray CT based numerical approach and compare to experimental measurements.
Three-dimensional high-resolution X-ray CT enables the numerical calculation of petrophysical properties of interest at the pore scale with resolutions down to a few microns per voxel. For more complex and heterogeneous samples however, a direct calculation of petrophysical properties is not feasible, since the required resolution and a sufficient field of view cannot be obtained simultaneously. Thus an integration of measurements at different scale is required. In this study a carbonate sample of 38mm in diameter is first scanned using the X-ray CT method with a resolution of 26µm. After accompanying experimental measurements on the full plug, four 5mm plugs were drilled vertically from this sample and X-ray CT images of these plugs acquired at resolutions down to 2.74 µm. We calculate the porosity of the sample (macro- and micro-porosities) using the phase separation methods and then predict the formation factor of the sample at several scales using a Laplace solver. The formation factor is calculated by using a general value of m=2 as cementation factor for intermediate porosity voxels. We compare to experimental measurements of formation factor and porosity both at the small plug and full plug scale and find good agreement.
To assess the degree of uncertainty of the numerical estimate, we probe the extent of heterogeneity by investigating the size of a representative elementary volume (REV) for formation factor. We find that for the considered heterogeneous carbonate sam-ple, formation factor varies considerably over intervals less than a centimetre. Our results show that this variation could be explained by different cementation exponents applied at the micro-voxel scale, with the exemption of one plug, for which the cementation exponent would have to be unreasonably low. These cementation factors are derived by direct comparison be-tween numerical simulation and experiment. We conclude that for one plug an error in experimental measurement might have occurred. The numerical approach presented here therefore aids in quality control. Excluding this plug in the upscaling proce-dure improves the agreement with the experimental result for the whole core while still underestimating formation factor. Al-lowing for a constant m=2 in the simulation at the small scale and using directly the resulting relationship between porosity and formation factor in the upscaling process leads to an overestimation of formation factor.
Feali, Mostafa (University of New South Wales) | Pinczewski, Wolf (University of New South Wales) | Cinar, Yildiray (University of New South Wales) | Arns, Christoph H. (University of New South Wales) | Arns, Ji-Youn (University of New South Wales) | Francois, Nicolas (Australian National University) | Turner, Michael L. (Australian National University) | Senden, Tim (Australian National University) | Knackstedt, Mark A. (Digitalcore)
It is now widely acknowledged that continuous oil-spreading films observed in 2D glass-micromodel studies for strongly water-wet three-phase oil, water, and gas systems are also present in real porous media, and they result in lower tertiary-gasflood residual oil saturations than for corresponding negative spreading systems that do not display oil-spreading behavior. However, it has not yet been possible to directly confirm the presence of continuous spreading films in real porous media in three dimensions, and little is understood of the distribution of the phases within the complex geometry and topology of actual porous media for different spreading conditions. This paper describes a study with high-resolution X-ray microtomography to image the distribution of oil, water, and gas after tertiary gasflooding to recover waterflood residual oil for two sets of fluids, one positive spreading and the other negative spreading, in strongly water-wet Bentheimer sandstone. We show that, for the positive spreading system, oil-spreading films maintain the connectivity of the oil phase down to low oil saturation. At similar oil saturation, no oil films are observed for the negative spreading system, and the oil phase is disconnected. The spatial continuity of the oil-spreading films over the imaged volume is confirmed by the computed Euler characteristic for the oil phase.
Pevzner, Roman (Curtin University) | Galvin, Robert J. (Curtin University) | Madadi, Mahyar (Curtin University) | Urosevic, Milovan (Curtin University) | Caspari, Eva (Curtin University) | Gurevich, Boris (Curtin University) | Lumley, David (University of Western Australia) | Shulakova, Valeriya (CSIRO) | Cinar, Yildiray (University of New South Wales) | Tcheverda, Vladimir (Trofimuk Institute of Geology)
Feali, Mostafa (U Of New South Wales) | Pinczewski, Wolf Val (U. of New South Wales) | Cinar, Yildiray (U. of New South Wales) | Arn, Christoph (U. of New South Wales) | Arns, Ji-Youn (Australian National Univ.) | Francois, Nicolas (Australian National Univ.) | Turner, Michael (Digital Core Laboratories Pty Ltd) | Senden, Tim | Knackstedt, Mark
It is now widely acknowledged that continuous oil spreading films observed in two-dimensional glass micro-model studies for strongly water wet three-phase oil, water and gas systems are also present in real porous media and result in lower tertiary gas flood residual oil saturations than for corresponding negative spreading systems which do not display oil spreading behavior. However, it has not been possible to directly confirm the presence of spreading films in real porous media in threedimensions and little is understood of the distribution of the phases within the complex geometry and topology of actual porous
media for different spreading conditions. This paper describes a preliminary study using high resolution X-ray microtomography to image the distribution of oil, water and gas after tertiary gas flooding to recover waterflood residual oil for two set of fluids, one positive spreading and the other negative spreading, for strongly water wet conditions in Bentheimer sandstone.
We show that for strongly water-wet conditions and a positive spreading system the oil phase remains connected throughout the pore space and results in a low tertiary gas flood residual oil saturation. The residual oil saturation for the corresponding negative spreading system is significantly higher and this is shown to be related to the absence of oil films in this system. The presence of films for positive spreading systems and the absence of such films for negative spreading systems is further confirmed by the computation of the Eurler characteristic for each phase.
Dehghan Khalili, Ahmad (U Of New South Wales) | Arns, Christoph Hermann (University of New South Wales) | Arns, Jiyoun (U. of New South Wales) | Hussain, Furqan (U. of New South Wales) | Cinar, Yildiray (U. of New South Wales) | Pinczewski, Wolf Val (Australian National University) | Latham, Shane (Saudi Aramco) | Funk, James Joseph
High-resolution Xray-CT images are increasingly used to numerically derive petrophysical properties of interest at the pore scale, in particular effective permeability. Current micro Xray-CT facilities typically offer a resolution of a few microns per voxel resulting in a field of view of about 5 mm3 for a 20482 CCD. At this scale the resolution is normally sufficient to resolve pore space connectivity and transport properties. For samples exhibiting heterogeneity above the field of view of such a single high resolution tomogram with resolved pore space, a second low resolution tomogram can provide a larger scale porosity
map. The problem then reduces to rock-typing the low resolution Xray-CT image, deriving viable porosity-permeability transforms from the high resolution Xray-CT image(s) for all rock types present, and upscaling of the permeability field to derive a plug-scale permeability.
In this study we characterize spatially heterogeneity using overlapping registered Xray-CT images derived at different resolutions spanning orders of magnitude in length scales. A 38mm diameter carbonate core is studied in detail and imaged at low resolution - and at high resolution by taking four 5mm diameter subsets, one of which is imaged using full length helical scanning. Fine-scale permeability transforms are derived using direct porosity-permeability relationships, random sampling of the porosity-permeability scatter-plot as function of porosity, and structural correlations combined with stochastic simulation. A range of these methods are applied at the coarse scale. We compare various upscaling methods including renormalization theory with direct solutions using a Laplace solver and report error bounds.
We find that for the heterogeneous samples permeability typically increases with scale. Conventional methods using basic averaging techniques fail to provide truthful vertical permeability due to large permeability contrasts. The most accurate upscaling technique is employing Darcy's law. A key part of the study is the establishment of porosity transforms between highresolution and low-resolution images to arrive at a calibrated porosity map to constraint permeability estimates for the whole core.