ALBERT

Description: Image processing of X-ray-computed polychromatic cone-beam micro-tomography (μXCT) data of geological samples mainly involves artefact reduction and phase segmentation. For the former, the main beam-hardening (BH) artefact is removed by applying a best-fit quadratic surface algorithm to a given image data set (reconstructed slice), which minimizes the BH offsets of the attenuation data points from that surface. A Matlab code for this approach is provided in the Appendix. The final BH-corrected image is extracted from the residual data or from the difference between the surface elevation values and the original grey-scale values. For the segmentation, we propose a novel least-squares support vector machine (LS-SVM, an algorithm for pixel-based multi-phase classification) approach. A receiver operating characteristic (ROC) analysis was performed on BH-corrected and uncorrected samples to show that BH correction is in fact an important prerequisite for accurate multi-phase classification. The combination of the two approaches was thus used to classify successfully three different more or less complex multi-phase rock core samples.

Description: 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.

Description: To date, very little is known about the distribution of natural gas hydrates in sedimentary matrices and its influence on the seismic properties of the host rock, in particular at low hydrate concentration. Digital rock physics offers a unique approach to this issue yet requires good quality, high-resolution 3-D representations for the accurate modeling of petrophysical and transport properties. Although such models are readily available via in situ synchrotron radiation X-ray tomography, the analysis of such data asks for complex workflows and high computational power to maintain valuable results. Here, we present a best-practice procedure complementing data from Chaouachi et al. (2015) with data post-processing, including image enhancement and segmentation as well as exemplary numerical simulations of an acoustic wave propagation in 3-D using the derived results. A combination of the tomography and 3-D modeling opens a path to a more reliable deduction of properties of gas hydrate-bearing sediments without a reliance on idealized and frequently imprecise models.

Description: The goal of this study is to better understand the porosity and permeability in shales to improve modelling fluid and gas flow related to shale diagenesis. Two samples (WIC and HAD) were investigated, both mid-Jurassic organic-rich Posidonia shales from Hils area, central Germany of different maturity (WIC R0 0.53 % and HAD R0 1.45 %). The method for image collection was focused ion beam (FIB) microscopy coupled with scanning electron microscopy (SEM). For image and data analysis Avizo and GeoDict was used. Porosity was calculated from segmented 3-D FIB based images and permeability was simulated by a Navier Stokes–Brinkman solver in the segmented images. Results show that the quantity and distribution of pore clusters and pores (≥  40 nm) are similar. The largest pores are located within carbonates and clay minerals, whereas the smallest pores are within the matured organic matter. Orientation of the pores calculated as pore paths showed minor directional differences between the samples. Both samples have no continuous connectivity of pore clusters along the axes in the x, y, and z direction on the scale of 10 to 20 of micrometer, but do show connectivity on the micrometer scale. The volume of organic matter in the studied volume is representative of the total organic carbon (TOC) in the samples. Organic matter does show axis connectivity in the x, y, and z directions. With increasing maturity the porosity in organic matter increases from close to 0 to more than 5 %. These pores are small and in the large organic particles have little connection to the mineral matrix. Continuous pore size distributions are compared with mercury intrusion porosimetry (MIP) data. Differences between both methods are caused by resolution limits of the FIB-SEM and by the development of small pores during the maturation of the organic matter. Calculations show no permeability when only considering visible pores due to the lack of axis connectivity. Adding the organic matter with a background permeability of 1 × 10−21 m2 to the calculations, the total permeability increased by up to 1 order of magnitude for the low mature and decreases slightly for the overmature sample from the gas window. Anisotropy of permeability was observed. Permeability coefficients increase by 1 order of magnitude if simulations are performed parallel to the bedding. Our results compare well with experimental data from the literature suggesting that upscaling may be possible in the future as soon as maturity dependent organic matter permeability coefficients can be determined.

Description: Various geoscientific applications require a fast prediction of fracture permeability for an optimal workflow. Hence, the objective of the current study is to introduce and validate a practical method to characterize and approximate single flow in fractures under different stress conditions by using a core-flooding apparatus, in situ X-ray computed tomography (CT) scans and a finite-volume method solving the Navier–Stokes–Brinkman equations. The permeability of the fractured sandstone sample was measured stepwise during a loading–unloading cycle (0.7 to 22.1 MPa and back) to validate the numerical results. Simultaneously, the pressurized core sample was imaged with a medical X-ray CT scanner with a voxel dimension of 0.5  ×  0.5  ×  1.0 mm3. Fracture geometries were obtained by CT images based on a modification of the simplified missing attenuation (MSMA) approach. Simulation results revealed both qualitative plausibility and a quantitative approximation of the experimentally derived permeabilities. The qualitative results indicate flow channeling along several preferential flow paths with less pronounced tortuosity. Significant changes in permeability can be assigned to temporal and permanent changes within the fracture due to applied stresses. The deviations of the quantitative results appear to be mainly caused by both local underestimation of hydraulic properties due to compositional matrix heterogeneities and the low CT resolution affecting the accurate capturing of sub-grid-scale features. Both affect the proper reproduction of the actual connectivity and therefore also the depiction of the expected permeability hysteresis. Furthermore, the threshold value CTmat (1862.6 HU) depicting the matrix material represents the most sensitive input parameter of the simulations. Small variations of CTmat can cause enormous changes in simulated permeability by up to a factor of 2.6 ± 0.1 and, thus, have to be defined with caution. Nevertheless, comparison with further CT-based flow simulations indicates that the proposed method represents a valuable method to approximate actual permeabilities, particularly for smooth fractures (