Search Results
Benchmark study using a multi-scale, multi-methodological approach for the petrophysical characterization of reservoir sandstones
2021, Haruzi, Peleg, Katsman, Regina, Halisch, Matthias, Waldmann, Nicolas, Spiro, Baruch
This paper presents a detailed description and evaluation of a multi-methodological petrophysical approach for the comprehensive multi-scale characterization of reservoir sandstones. The suggested methodology enables the identification of links between Darcy-scale permeability and an extensive set of geometrical, textural and topological rock descriptors quantified at the pore scale. This approach is applied to the study of samples from three consecutive sandstone layers of Lower Cretaceous age in northern Israel. These layers differ in features observed at the outcrop, hand specimen, petrographic microscope and micro-CT scales. Specifically, laboratory porosity and permeability measurements of several centimetre-sized samples show low variability in the quartz arenite (top and bottom) layers but high variability in the quartz wacke (middle) layer. The magnitudes of this variability are also confirmed by representative volume sizes and by anisotropy evaluations conducted on micro-CT-imaged 3-D pore geometries. Two scales of directional porosity variability are revealed in quartz arenite sandstone of the top layer: the pore size scale of ∼0.1 mm in all directions and ∼3.5 mm scale related to the occurrence of high- and low-porosity horizontal bands occluded by Fe oxide cementation. This millimetre-scale variability controls the laboratory-measured macroscopic rock permeability. More heterogeneous pore structures were revealed in the quartz wacke sandstone of the intermediate layer, which shows high inverse correlation between porosity and clay matrix in the vertical direction attributed to depositional processes and comprises an internal spatial irregularity. Quartz arenite sandstone of the bottom layer is homogenous and isotropic in the investigated domain, revealing porosity variability at a ∼0.1 mm scale, which is associated with the average pore size. Good agreement between the permeability upscaled from the pore-scale modelling and the estimates based on laboratory measurements is shown for the quartz arenite layers. The proposed multi-methodological approach leads to an accurate petrophysical characterization of reservoir sandstones with broad ranges of textural, topological and mineralogical characteristics and is particularly applicable for describing anisotropy and heterogeneity of sandstones on various rock scales. The results of this study also contribute to the geological interpretation of the studied stratigraphic units.
Enhanced pore space analysis by use of μ-CT, MIP, NMR, and SIP
2018, Zhang, Zeyu, Kruschwitz, Sabine, Weller, Andreas, Halisch, Matthias
We investigate the pore space of rock samples with respect to different petrophysical parameters using various methods, which provide data on pore size distributions, including micro computed tomography (μ-CT), mercury intrusion porosimetry (MIP), nuclear magnetic resonance (NMR), and spectral-induced polarization (SIP). The resulting cumulative distributions of pore volume as a function of pore size are compared. Considering that the methods differ with regard to their limits of resolution, a multiple-length-scale characterization of the pore space is proposed, that is based on a combination of the results from all of these methods. The approach is demonstrated using samples of Bentheimer and Röttbacher sandstone. Additionally, we compare the potential of SIP to provide a pore size distribution with other commonly used methods (MIP, NMR). The limits of resolution of SIP depend on the usable frequency range (between 0.002 and 100 Hz). The methods with similar resolution show a similar behavior of the cumulative pore volume distribution in the overlapping pore size range. We assume that μ-CT and NMR provide the pore body size while MIP and SIP characterize the pore throat size. Our study shows that a good agreement between the pore radius distributions can only be achieved if the curves are adjusted considering the resolution and pore volume in the relevant range of pore radii. The MIP curve with the widest range in resolution should be used as reference.
Development of a numerical workflow based on μ-CT imaging for the determination of capillary pressure–saturation-specific interfacial area relationship in 2-phase flow pore-scale porous-media systems: a case study on Heletz sandstone
2016, Peche, Aaron, Halisch, Matthias, Bogdan Tatomir, Alexandru, Sauter, Martin
In this case study, we present the implementation of a finite element method (FEM)-based numerical pore-scale model that is able to track and quantify the propagating fluid–fluid interfacial area on highly complex micro-computed tomography (μ-CT)-obtained geometries. Special focus is drawn to the relationship between reservoir-specific capillary pressure (pc), wetting phase saturation (Sw) and interfacial area (awn). The basis of this approach is high-resolution μ-CT images representing the geometrical characteristics of a georeservoir sample. The successfully validated 2-phase flow model is based on the Navier–Stokes equations, including the surface tension force, in order to consider capillary effects for the computation of flow and the phase-field method for the emulation of a sharp fluid–fluid interface. In combination with specialized software packages, a complex high-resolution modelling domain can be obtained. A numerical workflow based on representative elementary volume (REV)-scale pore-size distributions is introduced. This workflow aims at the successive modification of model and model set-up for simulating, such as a type of 2-phase problem on asymmetric μ-CT-based model domains. The geometrical complexity is gradually increased, starting from idealized pore geometries until complex μ-CT-based pore network domains, whereas all domains represent geostatistics of the REV-scale core sample pore-size distribution. Finally, the model can be applied to a complex μ-CT-based model domain and the pc–Sw–awn relationship can be computed.