Predicting elastic properties and permeability of rocks from 2D thin sections

2018 ◽  
Vol 37 (6) ◽  
pp. 421-427 ◽  
Author(s):  
Nattavadee Srisutthiyakorn ◽  
Sander Hunter ◽  
Rituparna Sarker ◽  
Ronny Hofmann ◽  
Irene Espejo

Predicting rock elastic properties and permeability from high-resolution 2D thin sections has been a challenging problem in rock physics because the 2D thin sections reveal very little about how the microstructure connects in the third dimension. However, 2D thin sections are widely available and inexpensive because they are often produced as a part of the reservoir-quality workflow. Furthermore, they have much higher resolution and greater field of view than micro X-ray computed tomography images, which are commonly used for rock properties estimation. The 2D thin sections we studied are from various hydrocarbon-bearing clastic formations with a variety of provenances, depositional environments, and burial histories. The high-resolution 2D images were scanned from these physical 2D thin sections. K-means segmentation was then employed to identify different minerals and pores for creating 2D binary images. The focus of this study is to simulate 2D elastic properties and permeability from 2D thin sections and then to employ various empirical relations to transform these 2D simulation results to 3D intrinsic rock properties. We compared the rock properties from this process to those from core measurements and measured wireline logs and found that these 2D to 3D rock property transformations yield promising results, especially for elastic properties. The results show that 2D thin section images have high enough resolution to resolve grain contacts very well. Predicting the permeability from 2D thin sections is still challenging since the process requires fitting the physical equation in order to retrieve the fitting coefficient for prediction due to our lack of understanding of the difference between 2D and 3D pore size distribution.

2016 ◽  
Author(s):  
Steven Henkel ◽  
Dieter Pudlo ◽  
Frieder Enzmann ◽  
Viktor Reitenbach ◽  
Daniel Albrecht ◽  
...  

Abstract. An essential part of the collaborative research project H2STORE ("hydrogen to store"), which is founded by the German government, was a comparison of various analytical methods to characterize reservoir sandstones from different stratigraphic units. In this context Permian, Triassic and Tertiary reservoir sandstones were analysed. Rock core materials, provided by RWE Gasspeicher GmbH (Dortmund), GFD Suez E&P Deutschland GmbH (Lingen), E.ON Gas Storage GmbH (Essen) and RAG Rohöl-Aufsuchungs Aktiengesellschaft (Wien), was processed by different laboratory techniques; thin sections were prepared, rock fragments were crushed, cubes of 1 cm edge length and plugs of 5 cm in length were sawn from macroscopic homogenous cores. With this prepared sample material, polarized light microscopy and scanning electron microscopy – coupled with image analyses, specific surface area measurements (BET), He-porosity and N2-permeability measurements and high resolution micro-computer-tomography (µ-CT), which were used for numerical simulations were conducted. All these methods were applied to most of the same sample material, before and after static CO2 experiments under reservoir conditions. A major concern in comparing the results of these methods is an appraisal of the reliability of the given porosity, permeability and mineral specific reactive (inner) surface areas data. The CO2 experiments are modifying the petrophysical as well the mineralogical/geochemical rock properties. These changes are detectable by all applied analytical methods. Nevertheless, a major outcome of the high resolution µ-CT analyses and proceeded numerical data simulations results in quite similar data sets and data interpretations maintained by the different standard methods; even regarding only CT-single scan of the rock samples. Moreover, this technique is not only time saving, but also none destructive. This is an important point, if only minor sample material is available and a detailed comparison before and after the experimental tests on micro meter, pore scale of specific rock features is envisaged.


Geophysics ◽  
2021 ◽  
pp. 1-76
Author(s):  
Jin Hao ◽  
Guoliang Li ◽  
Jiao Su ◽  
Yuan Yuan ◽  
Zhongming Du ◽  
...  

Digital rock physics (DRP) is an emerging technique that has rapidly become an indispensable tool to estimate elastic properties. The success of DRP mainly depends on three factors: acquiring a 3D rock structure image, accurately identifying 3D minerals, and using a proper numerical simulation method. Shales present a substantial challenge for DRP owing to their heterogeneous structure, composition, and properties from micron to centimeter scale. To obtain a sufficiently large field-of-view (FOV) image of a sample that reflects the detailed and representative internal structure and composition, we have developed a new DRP workflow to obtain large-FOV high-resolution digital rocks with 3D mineralogical information. Using the “divide-and-stitch” technique, a long shale sample is divided into several subunits, imaged separately by high-resolution X-ray microscopy (XRM), and then stitched to obtain a large-FOV 3D digital rock. An FOV of a rock cylinder (736 μm in diameter, 2358 μm in height, and 1 μm resolution) is used as an example. By correlating XRM and automated mineralogy, a large-FOV 3D mineral digital rock is obtained from a shale sample. Six mineral phases are identified and verified by automated mineralogy, and four laminae are detected according to the grain size, which offer a new perspective to study sedimentary processes and heterogeneities at the millimeter scale. The finite-difference method is used to compute the elastic properties of the large-FOV 3D mineral digital rock, and the results of Young’s modulus are within the limit of the Voigt/Reuss bounds. It also reveals that there is a difference in simulated elastic properties in the four laminae. The large-FOV 3D mineral digital rock offers the potential to explore the relationship between elastic properties and mineral phases, as well as the heterogeneities of elastic properties at the millimeter scale.


2020 ◽  
Author(s):  
Martin Balcewicz ◽  
Erik H. Saenger

<p>Digital rock physics (DRP) became a complementary part in reservoir characterization during the last two decades. Deriving transport, thermal, or effective elastic rock properties from a digital twin requires a three-step workflow: (1) Preparation of a high-resolution X-ray computed tomography image, (2) segmentation of pore and grain phases, respectively, and (3) solving equations due to the demanded properties. Despite the high resolution µ-CT images, the numerical predictions of rock properties have their specific uncertainties compared to laboratory measurements. Missing unresolved features in the µ-CT image might be the key issue. These findings indicate the importance of a full understanding of the rocks microfabrics. Most digital models used in DRP treat the rock as a heterogeneous, isotropic, intact medium which neglect unresolved features. However, we expect features within the microfabrics like micro-cracks, small-scale fluid inclusions, or stressed grains which may influence the elastic rock properties but have not been taken into account in DRP, yet. Former studies have shown resolution-issues in grain-to-grain contacts within sandstones or inaccuracies due to micro-porosity in carbonates, this means the micritic phase. Within the scope of this abstract, we image two different sandstone samples, Bentheim and Ruhrsandstone, as well as one carbonate sample. Here, we compare the mentioned difficulties of X-ray visualization with further analytical methods, this means thin section and focused ion beam measurements. This results into a better understanding of the rocks microstructures and allows us to segment unresolved features in the X-ray computed tomography image. Those features can be described with effective properties at the µ-scale in the DRP workflow to reduce the uncertainty of the predicted rock properties at the meso- and fieldscale.</p>


2017 ◽  
Vol 5 (2) ◽  
pp. B17-B27 ◽  
Author(s):  
Mark Sams ◽  
David Carter

Predicting the low-frequency component to be used for seismic inversion to absolute elastic rock properties is often problematic. The most common technique is to interpolate well data within a structural framework. This workflow is very often not appropriate because it is too dependent on the number and distribution of wells and the interpolation algorithm chosen. The inclusion of seismic velocity information can reduce prediction error, but it more often introduces additional uncertainties because seismic velocities are often unreliable and require conditioning, calibration to wells, and conversion to S-velocity and density. Alternative techniques exist that rely on the information from within the seismic bandwidth to predict the variations below the seismic bandwidth; for example, using an interpretation of relative properties to update the low-frequency model. Such methods can provide improved predictions, especially when constrained by a conceptual geologic model and known rock-physics relationships, but they clearly have limitations. On the other hand, interpretation of relative elastic properties can be equally challenging and therefore interpreters may find themselves stuck — unsure how to interpret relative properties and seemingly unable to construct a useful low-frequency model. There is no immediate solution to this dilemma; however, it is clear that low-frequency models should not be a fixed input to seismic inversion, but low-frequency model building should be considered as a means to interpret relative elastic properties from inversion.


2021 ◽  
Author(s):  
Martin Balcewicz ◽  
Mirko Siegert ◽  
Marcel Gurris ◽  
David Krach ◽  
Matthias Ruf ◽  
...  

<p>Over the last two decades, Digital Rock Physics (DRP) has become a complementary part of the characterization of reservoir rocks due to, among other things, the non-destructive testing character of this technique. The use of high-resolution X-ray Computed Tomography (XRCT) has become widely accepted to create a digital twin of the material under investigation. Compared to other imaging techniques, XRCT technology allows a location-dependent resolution of the individual material particles in volume. However, there are still challenges in assigning physical properties to a particular voxel within the digital twin, due to standard histogram analysis or sub-resolution features in the rock. For this reason, high-resolution image-based data from XRCT, transmitted-light microscope, Scanning Electron Microscope (SEM) as well as inherent material properties like porosity are combined to obtain an optimal spatial image of the studied Ruhr sandstone by a geologically driven segmentation workflow. On the basis of a homogeneity test, which corresponds to the evaluation of the grayscale image histogram, the preferred scan sample sizes in terms of transport, thermal, and effective elastic rock properties are determined. In addition, the advanced numerical simulation results are compared with laboratory tests to provide possible upper limits for sample size, segmentation accuracy, and a calibrated digital twin of the Ruhr sandstone. The comparison of representative grayscale image histograms as a function of sample sizes with the corresponding advanced numerical simulations, provides a unique workflow for reservoir characterization of the Ruhr sandstone.</p>


Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. D53-D64 ◽  
Author(s):  
Claudio Madonna ◽  
Beatriz Quintal ◽  
Marcel Frehner ◽  
Bjarne S. G. Almqvist ◽  
Nicola Tisato ◽  
...  

Synchrotron radiation X-ray tomographic microscopy is a nondestructive method providing ultra-high-resolution 3D digital images of rock microstructures. We describe this method and, to demonstrate its wide applicability, we present 3D images of very different rock types: Berea sandstone, Fontainebleau sandstone, dolomite, calcitic dolomite, and three-phase magmatic glasses. For some samples, full and partial saturation scenarios are considered using oil, water, and air. The rock images precisely reveal the 3D rock microstructure, the pore space morphology, and the interfaces between fluids saturating the same pore. We provide the raw image data sets as online supplementary material, along with laboratory data describing the rock properties. By making these data sets available to other research groups, we aim to stimulate work based on digital rock images of high quality and high resolution. We also discuss and suggest possible applications and research directions that can be pursued on the basis of our data.


Solid Earth ◽  
2016 ◽  
Vol 7 (3) ◽  
pp. 917-927 ◽  
Author(s):  
Steven Henkel ◽  
Dieter Pudlo ◽  
Frieder Enzmann ◽  
Viktor Reitenbach ◽  
Daniel Albrecht ◽  
...  

Abstract. An essential part of the collaborative research project H2STORE (hydrogen to store), which is funded by the German government, was a comparison of various analytical methods for characterizing reservoir sandstones from different stratigraphic units. In this context Permian, Triassic and Tertiary reservoir sandstones were analysed. Rock core materials, provided by RWE Gasspeicher GmbH (Dortmund, Germany), GDF Suez E&P Deutschland GmbH (Lingen, Germany), E.ON Gas Storage GmbH (Essen, Germany) and RAG Rohöl-Aufsuchungs Aktiengesellschaft (Vienna, Austria), were processed by different laboratory techniques; thin sections were prepared, rock fragments were crushed and cubes of 1 cm edge length and plugs 3 to 5 cm in length with a diameter of about 2.5 cm were sawn from macroscopic homogeneous cores. With this prepared sample material, polarized light microscopy and scanning electron microscopy, coupled with image analyses, specific surface area measurements (after Brunauer, Emmet and Teller, 1938; BET), He-porosity and N2-permeability measurements and high-resolution microcomputer tomography (μ-CT), which were used for numerical simulations, were applied. All these methods were practised on most of the same sample material, before and on selected Permian sandstones also after static CO2 experiments under reservoir conditions. A major concern in comparing the results of these methods is an appraisal of the reliability of the given porosity, permeability and mineral-specific reactive (inner) surface area data. The CO2 experiments modified the petrophysical as well as the mineralogical/geochemical rock properties. These changes are detectable by all applied analytical methods. Nevertheless, a major outcome of the high-resolution μ-CT analyses and following numerical data simulations was that quite similar data sets and data interpretations were maintained by the different petrophysical standard methods. Moreover, the μ-CT analyses are not only time saving, but also non-destructive. This is an important point if only minor sample material is available and a detailed comparison before and after the experimental tests on micrometre pore scale of specific rock features is envisaged.


2021 ◽  
Vol 11 (6) ◽  
pp. 2497-2518
Author(s):  
Syed Haroon Ali ◽  
Osman M. Abdullatif ◽  
Lamidi O. Babalola ◽  
Fawwaz M. Alkhaldi ◽  
Yasir Bashir ◽  
...  

AbstractThis paper presents the facies and depositional environment of the early Miocene Dam Formation, Eastern Arabian platform, Saudi Arabia. Deposition of Dam Formation (Fm.) was considered as a restricted shallow marine deposition. Few studies suggest the role of sea-level change in its deposition but were without decisive substantiation. Here, we describe the facies and high-resolution model of Dam Fm. under varying depositional conditions. The depositional conditions were subjected to changing relative sea level and tectonics. High-resolution outcrop photographs, sedimentological logs, and thin sections present that the mixed carbonate–siliciclastic sequence was affected by a regional tectonics. The lower part of Dam Fm. presents the development of carbonate ramp conditions that are represented by limestones and marl. The depositional conditions fluctuated with the fall of sea level, and uplift in the region pushed the siliciclastic down-dip and covered the whole platform. The subsequent rise in sea level was not as pronounced and thus allowed the deposition of microbial laminites and stromatolitic facies. The southeast outcrops, down-dip, are more carbonate prone as compared to the northwest outcrop, which allowed the deposition of siliciclastic-prone sedimentation up-dip. All facies, architecture, heterogeneity, and deposition were controlled by tectonic events including uplift, subsidence, tilting, and syn-sedimentary faulting, consequently affecting relative sea level. The resulting conceptual outcrop model would help to improve our understanding of mixed carbonate–siliciclastic systems and serve as an analogue for other stratigraphic units in the Arabian plate and region. Our results show that Dam Fm. can be a good target for exploration in the Northern Arabian Gulf.


2015 ◽  
Vol 18 (03) ◽  
pp. 432-440 ◽  
Author(s):  
C.. Germay ◽  
T.. Richard ◽  
E.. Mappanyompa ◽  
C.. Lindsay ◽  
D.. Kitching ◽  
...  

Summary Knowledge of rock properties is essential to predict and optimize the performance of oil and gas reservoirs by means of the reduction of the uncertainty pertaining to standard subsurface issues such as the mechanical integrity of the borehole (Tiab and Donaldson 1996; Moos et al. 2003), the risk of sanding (Tronvoll et al. 2004), and the geometry and efficiency of hydraulic fractures. These properties are evaluated by combining different field-measurement techniques (wireline logs, results of well tests, seismic surveys) and laboratory-test results (Archie 1942, 1950; Serra 1986; Bassiouni 1994). When cores are available, empirical models are built from correlations derived between well logs and laboratory measurements to estimate rock properties in noncored wells. The validity of these empirical models is often limited to specific litho-facies (see reviews by Chang et al. 2006; Blasingame 2008; Khaksar et al. 2009), which makes the identification of lithofacies a necessity before applying the model for predictions in uncored wells (Massonnat 1999). Because of the heterogeneity of rocks (Haldorsen 1996), with characteristic length scales commonly smaller than the resolution of wireline logs or even the core-plug size, the robustness of correlations is determined by how plug samples capture the dispersion in rock properties over the lithofacies under consideration. The correlation between a very localized core-plug measurement and a low-resolution wireline log with inherent low-pass filtering properties raises issues related to the upscaling of a property from one length scale (few centimeters for core plugs) to another (up to 1 m for wireline log). As an illustration, consider the high-resolution, continuous profile X, where the variations of the measured property are quantified for length scales smaller than typical plug sizes. We filter this data to produce the profiles X5 and X50 (the subscript stands for the length scale in centimeters at which the signal is averaged out) with lower spatial resolutions similar to the plug and the well-log resolutions, respectively (Fig. 1). The resulting crossplot, shown in Fig. 2, of X5 vs. X50 exhibits a cloud of points in which the dispersion is governed by the properties of the signal (the degree of heterogeneity or the frequency content) and the difference between the two resolution length scales. Two linear-fit optimizations were carried out with the low-resolution-data X50 and the high-resolution-data X5 as the dependent variables, respectively. It is interesting to note that these linear fits yield different results, with a slope of 0.96 in the first case and 0.69 in the second case. This is a mathematical artifact caused by the minimization process inherent in the search for the best linear fit, which is most commonly a minimization of the vertical distance between the representative data points and the best-fit line. On the basis of this result, it should always be advisable to select the high-resolution data (plug) as the dependent variable. Discrete sampling (i.e., plugging) and the dispersion caused by the difference in resolution scales of two measurements are two important root causes of the errors often seen in correlations between two variables. The examples shown in Fig. 2 illustrate how the correlations derived from several sampling schemes can deviate from the expected one-to-one relation between the two variables. To circumvent these issues, petrophysicists usually select large quantities of plugs to build representative statistical data sets, with the hope that they are large enough to attenuate the effects listed previously. However, extensive plugging strategies imply longer lead times and higher costs, and are therefore not always viable (e.g., in the cases of rock-mechanics testing or special-core-analysis programs). As an illustration, consider the modeling of the variations of rock strength, one of the key geomechanical properties along a well trajectory. Such an exercise relies heavily on correlations derived between well logs and laboratory tests (uniaxial or triaxial compressive tests), because there is no wireline log providing a direct measure of a mechanical property related to strength. In their comprehensive review of existing literature, Khaksar et al. (2009) listed approximately 40 models designed to derive strength properties from wireline logs. The authors showed that the relevance of these as empirical is limited to specific rock types. A broader application of these models would require the considerations of additional complexity such as the coexistence of several facies within the same data set or the impact of diagenesis on petrophysical variability within one facies. The elements of reflection introduced previously all suggest that a continuous measure of a physical property such as the strength profiles generated from the scratch test, which provides some useful elements for the mapping of rock heterogeneity, could partially fill the gap between measurements on plugs and well logs and help with the optimization of the selection of plug samples. In the main sections of this paper, we first describe briefly the scratch test and outline the key intrinsic benefits of the test. We then discuss how standard and special core analysis could benefit most from all the features of the scratch test when introduced at a very early stage of the work flows. In particular, we illustrate with some examples how rock-strength profiles averaged to the relevant length scale can be correlated with other petrophysical properties either measured on core plugs or inferred from well logs.


Solid Earth ◽  
2016 ◽  
Vol 7 (4) ◽  
pp. 1185-1197 ◽  
Author(s):  
Erik H. Saenger ◽  
Stephanie Vialle ◽  
Maxim Lebedev ◽  
David Uribe ◽  
Maria Osorno ◽  
...  

Abstract. Modern estimation of rock properties combines imaging with advanced numerical simulations, an approach known as digital rock physics (DRP). In this paper we suggest a specific segmentation procedure of X-ray micro-computed tomography data with two different resolutions in the µm range for two sets of carbonate rock samples. These carbonates were already characterized in detail in a previous laboratory study which we complement with nanoindentation experiments (for local elastic properties). In a first step a non-local mean filter is applied to the raw image data. We then apply different thresholds to identify pores and solid phases. Because of a non-neglectable amount of unresolved microporosity (micritic phase) we also define intermediate threshold values for distinct phases. Based on this segmentation we determine porosity-dependent values for effective P- and S-wave velocities as well as for the intrinsic permeability. For effective velocities we confirm an observed two-phase trend reported in another study using a different carbonate data set. As an upscaling approach we use this two-phase trend as an effective medium approach to estimate the porosity-dependent elastic properties of the micritic phase for the low-resolution images. The porosity measured in the laboratory is then used to predict the effective rock properties from the observed trends for a comparison with experimental data. The two-phase trend can be regarded as an upper bound for elastic properties; the use of the two-phase trend for low-resolution images led to a good estimate for a lower bound of effective elastic properties. Anisotropy is observed for some of the considered subvolumes, but seems to be insignificant for the analysed rocks at the DRP scale. Because of the complexity of carbonates we suggest using DRP as a complementary tool for rock characterization in addition to classical experimental methods.


Sign in / Sign up

Export Citation Format

Share Document