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.

In geochemically perturbed systems where porewater and mineral assemblages are unequilibrated the processes of mineral precipitation and dissolution may change important transport properties such as porosity and pore diffusion coefficients. These reactions might alter the sealing capabilities of the rock by complete pore-scale precipitation (cementation) of the system or by opening new migration pathways through mineral dissolution. In actual 1D continuum reactive transport codes the coupling of transport and porosity is generally accomplished through the empirical Archie’s law. There is very little reported data on systems with changing porosity under well controlled conditions to constrain model input parameters. In this study celestite (SrSO4) was precipitated in the pore space of a compacted sand column under diffusion controlled conditions and the effect on the fluid migration properties was investigated by means of three complementary experimental approaches: (1) tritiated water (HTO) tracer through diffusion, (2) computed micro-tomography (µ-CT) imaging and (3) post-mortem analysis of the precipitate (selective dissolution, SEM/EDX). The through-diffusion experiments reached steady state after 15 days, at which point celestite precipitation ceased and the non-reactive HTO flux became constant. The pore space in the precipitation zone remained fully connected using a 6 µm µ-CT spatial resolution with 25 % porosity reduction in the approx. 0.35 mm thick dense precipitation zone. The porosity and transport parameters prior to pore-scale precipitation were in good agreement with a porosity of 0.42 ± 0.09 (HTO) and 0.40 ± 0.03 (µ-CT), as was the mass of SrSO4 precipitate estimated by µ-CT at 25 ± 5 mg and selective dissolution 21.7 ± 0.4 mg, respectively. However, using this data as input parameters the 1D single continuum reactive transport model was not able to accurately reproduce both the celestite precipitation front and the remaining connected porosity. The model assumed there was a direct linkage of porosity to the effective diffusivity using only one cementation value over the whole porosity range of the system investigated. The 1D single continuous model either underestimated the remaining connected porosity in the precipitation zone, or overestimated the amount of precipitate. These findings support the need to implement a modified, extended Archie’s law to the reactive transport model and show that pore-scale precipitation transforms a system (following Archie’s simple power law with only micropores present) towards a system similar to clays with micro- and nanoporosity.

Violence seems deeply rooted in human nature and an endemic potential for such is today frequently associated with differing ethnic, religious or socio-economic backgrounds. Ethnic nepotism is believed to be one of the main causes of inter-group violence in multi-ethnic societies. At the site of Els Trocs in the Spanish Pyrenees, rivalling groups of either migrating early farmers or farmers and indigenous hunter-gatherers collided violently around 5300 BCE. This clash apparently resulted in a massacre of the Els Trocs farmers. The overkill reaction was possibly triggered by xenophobia or massive disputes over resources or privileges. In the present, violence and xenophobia are controlled and sanctioned through social codes of conduct and institutions. So that, rather than representing an insurmountable evolutionary inheritance, violence and ethnic nepotism can be overcome and a sustainable future achieved through mutual respect, tolerance and openness to multi-ethnic societies.

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 R₀ 0.53 % and HAD R₀ 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⁻²¹ m² 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.

Computer X-ray microtomography (µXCT) represents a powerful tool for investigating the physical properties of porous rocks. While calculated porosities determined by this method typically match experimental measurements, computed permeabilities are often overestimated by more than 1 order of magnitude. This effect increases towards smaller pore sizes, as shown in this study, in which nanostructural features related to clay minerals reduce the permeability of tight reservoir sandstone samples. Focussed ion beam scanning electron microscopy (FIB-SEM) tomography was applied to determine the permeability effects of illites at the nanometre scale, and Navier–Stokes equations were applied to calculate the permeability of these domains. With these data, microporous domains (porous voxels) were defined using microtomography images of a tight reservoir sample. The distribution of these domains could be extrapolated by calibration against size distributions measured in FIB-SEM images. For this, we assumed a mean permeability for the dominant clay mineral (illite) in the rock and assigned it to the microporous domains within the structure. The results prove the applicability of our novel approach by combining FIB-SEM with X-ray tomographic rock core scans to achieve a good correspondence between measured and simulated permeabilities. This methodology results in a more accurate representation of reservoir rock permeability in comparison to that estimated purely based on µXCT images.