EMRP1.15

The continuously rising threat of global warming caused by human activities related to CO₂ emission is facilitating the development of greenhouse gas control technologies. Subsurface CO₂ injection and sequestration is one of the promising techniques to store CO₂ in the subsurface.  However, during CO₂ injection, the mechanisms of processes like injectant immobilizations and trapping and pore-scale geochemical reactions such as mineral dissolution/precipitation are not well understood. Consequently, the multi-physics modeling approach is essential to elucidate the impact of all potential factors during CO₂ injection, thus to facilitate the optimization of this engineered application. 

Here, we propose a coupled framework to fully utilize the capabilities of the geochemical reaction solver PHREEQC while preserving the Lattice-Boltzmann Method (LBM) high-resolution pore-scale fluid flow integrated with diffusion processes. The model can simulate the dynamic fluid-solid interactions with equilibrium, kinetics, and surface reactions under the reactive-transport scheme.  In a simplified 2D spherical pack, we focused on examining the impact of pore sizes, grain size distributions, porosity, and permeability on the calcite dissolution/precipitation rate. Our simulation results show that the higher permeability, injection rate, and more local pore connectivity would significantly increase the reaction rate, then accelerate the pore-scale geometrical evolutions. Meanwhile, model accuracy is not sacrificed by reducing the number of reactants/species within the system.

Our modeling framework provides high-resolution details of the pore-scale fluid-solid interaction dynamics. To gain more insights into the mineral-fluid interfacial properties during CO₂ sequestration, our next step is to combine the electrodynamic forces into the model. Potentially, the proposed framework can be used for model upscaling and adaptive subsurface management in the future.  

How to cite: Zhang, C., Liu, S., and Barati, R.: Pore-scale hydrodynamic evolution within carbonate rock during CO2 injection and sequestration, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9578, https://doi.org/10.5194/egusphere-egu21-9578, 2021.

Monitoring of fluid movements in the crust is one of the most discussed topics in oil & gas industry as well as in geothermal systems and CO₂ storage, but still remains a challenge. The seismic method is one of the most common ways to detect the fluid migration. However, the use of ultrasonic monitoring at the sample scale in laboratory experiments persists as the most effective way to highlight large scale observations in which the boundary conditions are not well constrained.

To unravel the fluid effect on P-wave and S-wave velocity, we performed mechanical experiments coupled with ultrasonic monitoring on Obourg chalk from Mons basin (Belgium). Water injection tests under critical loading, imbibition tests and evaporation tests provided a full spectrum of observations of fluid-induced wave alteration in term of propagation time and attenuation.

The analysis of these experimental results showed that significant velocity dispersion and attenuation developed through variations in water saturation, and that these processes are linked to the presence of patches of water and air in the pore space.

We used the White’s formulation to model the relaxation effects due to spherical pockets of air homogeneously distributed in a water-saturated medium. In this framework, the pressure induced by the passing wave, produces a fluid flow across the water-air boundary with consequent energy loss.

This model reproduces both qualitatively and quantitatively the experimental results observed on the water injection tests. Indeed, it is shown that the progressive water saturation or desaturation of this chalk, generates a shift of the critical frequency (from the undrained relaxed towards unrelaxed regimes) which at some point matches the resonance frequency of the piezoelectric transducers used in the experimental setup (0.5 MHz). This phenomenon allowed us to get a continuous recording of the relaxation processes induced by saturation variations.

The outcomes of this work can significantly improve the actual knowledge on coupled effects of waves and fluids which is a crucial aspect of fluid monitoring in the context of reservoir evaluation and production.

How to cite: Geremia, D., David, C., Barnes, C., Menéndez, B., Dautriat, J., Esteban, L., Sarout, J., Vandycke, S., and Descamps, F.: Laboratory experiments of water injection coupled with ultrasonic monitoring reveal wave-induced fluid flow in microporous carbonate rock, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9230, https://doi.org/10.5194/egusphere-egu21-9230, 2021.

Metamorphic and deformed rocks in thrust zones show particularly high seismic anisotropy causing challenges for seismic imaging and interpretation. A good example is the Seve Nappe Complex in Jämtland, Sweden, an exhumed orogenic thrust zone characterized by a strong but incoherent seismic reflectivity and considerable seismic anisotropy. However, only little is known about the origin of the anisotropy in relation to composition, structural influences, and implications for measurements at different seismic scales. We present an integrative study of the seismic anisotropy at different scales combining mineralogical composition, microstructural analyses and seismic laboratory experiments from samples of the 2.5 km-deep COSC-1 borehole. While there is a pronounced crystallographic preferred orientation in most of the core samples, variations in anisotropy correlate strongly with bulk mineral composition and dominant core lithology. Based on three major lithologic different facies (felsic gneiss, amphibole-rich rocks, and mica schists), we propose an anisotropy model for the full length of the borehole, which indicates two prevailing anisotropic units. Comparison of laboratory seismic measurements and electron-backscatter diffraction (EBSD) data reveals a strong scale-dependence, which is more pronounced in the highly deformed, heterogeneous samples. This highlights the need for comprehensive cross-validation of microscale anisotropy analyses with additional lithological data when integrating seismic anisotropy through seismic scales.

How to cite: Kästner, F., Pierdominici, S., Zappone, A., Morales, L. F. G., Schleicher, A. M., Wilke, F. D. H., and Berndt, C.: Seismic anisotropy in metamorphic rocks from the COSC-1 borehole, Sweden: A cross-scale investigation from thin section analysis to seismic scales, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-12483, https://doi.org/10.5194/egusphere-egu21-12483, 2021.

Fractures, which are ubiquitous in the Earth’s upper crust, have significant impacts on a wide range of human activities, and, hence, their adequate characterization is of wide interest and importance. Seismic methods have significant potential for effectively addressing this objective. When a seismic wave propagates across a fluid-filled fracture, its amplitude is diminished and its travel time is increased. Based on the linear slip theory, the associated amplitude decays and phase delays can be used to estimate the mechanical compliance of fractures. Full-waveform sonic (FWS) log data are particularly well-suited for this purpose. While the amplitudes of FWS data acquired during standard continuous logging runs (tool being moved uphole at a constant logging speed) can be somewhat unstable, the associated first-arrival travel times are generally quite robust. In this work, we exploit the relation between the time delay that seismic waves experience across fractures and relate them to the associated compliances. Specifically, we estimate fracture compliance from the differences in group time delay of the refracted P-wave between fractured and non-fractured sections along a borehole. Numerical simulations indicate that the proposed method provides reliable compliance estimates not only for individual fractures, but also for sets of multiple discrete fractures. This finding is corroborated by applying our approach to FWS log data acquired in the course of standard logging runs in the Bedretto Underground Laboratory (www.bedrettolab.ethz.ch). Our estimates are comparable to previously inferred compliance values in a closely comparable geological environment (Grimsel test site, www.grimsel.com). The latter were inferred under rather ideal conditions, involving the quasi-static acquisition of the FWS data as well as the combination of amplitude and travel time information for their interpretation. An interesting and important open question, which we plan to address in the following, concerns the influence of the heterogeneity of the host rock embedding the fractures on compliance estimation in general and on the proposed method in particular.

How to cite: Zhou, Z., Caspari, E., D. Barbosa, N., Greenwood, A., and Holliger, K.: Estimation of fracture compliance based on group delays inferred from full-waveform sonic log data, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-10820, https://doi.org/10.5194/egusphere-egu21-10820, 2021.

We are using 3D magnetic vector inversion (MVI) of scanning magnetic microscopy (SMM) data to investigate the fine‐scale magnetization of rock samples, and particularly of their remanence carriers, which can record geologically meaningful information. Previous investigations of magnetite grains suggest variable remanence intensities and directions coherent with multidomain behaviour. This research seeks to improve our understanding of the contribution of different microstructures on remanence acquisition.

SMM offers a spatial resolution down to tens of micrometers, allowing detailed investigation of discrete magnetic mineral grains, or magnetic textures and structures. However, all magnetic measurements are, at some scale, bulk measurements. Further analysis of the data is required to extract information about the magnetization within the samples: for this, we employ state-of-the-art MVI methods. The MVI problem suffers from a high degree of nonuniqueness. Additional constraints are required to obtain accurate, reliable and interpretable results. Such constraints are readily available for this application.

SMM instruments use magnetic shields or Helmholtz coils to allow collection of data in controlled magnetic fields, enabling the removal of induced magnetization effects. Measurements can be taken both above and below the sample. Individual magnetized mineral grains are easily outlined through optical and electron microscopy. The internal geometry of the oxide mineral phases and compositions can also be constrained. Physical property information constrains the range of magnetization intensity. As such, there is a tremendous amount of constraining information invaluable for reducing the nonuniqueness of the inverse problem. We use a highly flexible and functional inversion software package, MAGNUM, developed jointly at Mount Allison University and Memorial University of Newfoundland, that allows incorporation of all available constraints.

We take a multitiered approach for investigating specific magnetized grains. First, coarse regional inversions are performed to assess and remove any effects of other magnetized grains in the vicinity. The entire grain is then modelled with a homogeneous magnetization to obtain an approximate but representative bulk magnetization. The grain is then modelled as a collection of independent subdomains, each with a different homogeneous magnetization direction. Subsequently, more heterogeneous scenarios are considered by relaxing inversion constraints until the data can be fit to the desired degree.

Obtaining reliable information about the magnetic mineralogy of rock samples is vital for an understanding of the origin of rock bulk behaviour in both the laboratory and larger scale magnetic surveys. This work is among the first to simultaneously invert SMM data collected above and below a thin sample, which is critical for improving depth resolution on thicker samples. It is also the first time we have been able to incorporate all available constraints into inverse modelling to improve results.

How to cite: Lelièvre, P., Pastore, Z., Church, N., Lee, M., Oda, H., and McEnroe, S.: Constrained Magnetic Vector Inversion of Scanning Magnetic Microscopy Data for Modelling Magnetization of Multidomain Mineral Grains, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3466, https://doi.org/10.5194/egusphere-egu21-3466, 2021.

Seismic waves are affected by rock properties such as porosity, permeability, grain material and by their heterogeneities as well as by the fluid properties saturating the rocks. Consequently, seismic methods are a valuable tool for the indirect characterization of rocks. For example, at the microscale, the presence of compliant pores (cracks or grain contacts) in fluid-saturated rocks can cause strong seismic attenuation and velocity dispersion. In this case, the deformation caused by a passing wave induces a fluid pressure gradient between compressed compliant pores and much less compressed pores (stiff isometric pores or cracks having a different orientation than the most compressed ones) if they are hydraulically connected. The consequent fluid pressure diffusion (FPD) dissipates seismic energy due to viscous friction in the fluid.

Digital rock physics (DRP) aims to reproduce experimental measurements using numerical simulation in models derived from high resolution rock images. We developed a DRP workflow to calculate the frequency dependent seismic moduli dispersion and attenuation in fluid-saturated models derived from micro X-Ray Computed Tomography (µXRCT) images. Filtering, segmentation and meshing procedures are applied on sub-volumes of different rock images to create 3D numerical models. We apply our workflow to calculate seismic moduli attenuation due to FPD at the microscale (squirt flow). We consider a µXRCT image of a cracked (through thermal treatment) Carrara marble sample. A detailed visualization of the fluid pressure as well as of the energy dissipation rate in the 3D model helps to understand the squirt flow attenuation process at different frequencies.

How to cite: Lissa, S., Ruf, M., Steeb, H., and Quintal, B.: Seismic moduli dispersion and attenuation obtained using Digital Rock Physics, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-5887, https://doi.org/10.5194/egusphere-egu21-5887, 2021.

An understanding of micro- and macrofracture behavior in low porosity rocks is pertinent to several areas of energy and environmental science such as petroleum production, carbon sequestration, and enhancement of technologies based on geothermal energy, etc. For example, the carbonate reservoirs in dolomitic or micritic formations with matrix porosities below 6% suggest the importance of fracture-augmented permeability in production. Similarly, hydrocarbons have been found on nearly every continent in tight basement rocks, all of which have little matrix porosity and their permeability therefore rely solely on hydraulic connectivity from fractures. For geothermal energy, various igneous and sedimentary rocks (granites, basalts, and limestones) are being exploited across the globe, with some of the lowest porosity and permeability. In all these cases, fractures are necessary to improve rock permeability and thermal exchange between the rock and working fluid, which can be enabled by hydraulic stimulation, as well as by secondary cracking due to extreme temperature gradients from the injection of cold water. The fracture geometry, density, and distribution all together control not only fluid and thermal transport in the rocks, but also their seismic attributes that can be used to extract information about the fractures. 
In order to accurately interpret the seismo-acoustic data (usually, the velocities of compression and shear waves) reliable rock physics models are required. Here, we report the results of interpretation of such experimental data for both as-cored rock samples and those subjected to thermo-hydro-chemo-mechanical damage (THCMD) in the laboratory. For interpretation, we use a convenient model of fractured rock in which fractures are represented as planar defects with decoupled shear and normal compliances. The application of such an approach makes it possible to assess and compare the elastic properties of fractures in the rocks before and after application of THCMD procedures. For the analyzed samples of granites, basalts, and limestones it has been found that for a significant portion of rocks, the ratio of normal-to-shear compliances of cracks significantly differ from the value typical of conventionally assumed penny-shape cracks. Furthermore, for some samples, this ratio appears to be noticeably different for fractures existing in the as-cored rock and arising in the same rock after THCMD procedures. These results indicate that damage to a rock typically changes its compliance ratio since the old and new cracks are likely to have different elastic properties. Our results are also consistent with the notion that a specific damage process occurring for a given microstructure will consistently create cracks with a particular set of elastic properties. The proposed methodology for assessment of elastic properties of cracks in rock samples subjected to thermo-hydro-chemo-mechanical damage has given previously inaccessible useful information about the elastic properties of fractures and can be extended to interpretation of seismic attributes of rocks for a broad range of other applications.

How to cite: Clark, A., Vanorio, T., Radostin, A., and Zaitsev, V.: Assessing elastic properties of cracks in rock samples subjected to thermo-hydro-chemo-mechanical damage, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-6492, https://doi.org/10.5194/egusphere-egu21-6492, 2021.

Within the scope of the present work, the pressure-dependent effective thermal conductivity of rock samples is simulated. Our workflow can be assigned to the field of digital rock physics. In a first step, a 3D micro-CT scan of a rock sample is taken. Subsequently, the resulting greyscale images are analysed and segmented depending on the occurring phases. Based on this data set, a computational mesh is created and the corresponding thermal conductivities are assigned to each phase. Finally the numerical simulations can be carried out.
For the representation of the pressure dependency we use the approach proposed by Saenger [1]. By making use of the watershed algorithm, boundaries between the individual grains of the rock sample are detected and assigned to an artificial contact phase. In the course of several simulations, the thermal conductivity of the contact phase is continuously increased. Starting with the thermal conductivity of the pore phase and ending with the thermal conductivity of the grain phase. A linear correlation is used to match the thermal conductivity of the contact phase with the pressure of a given experimental data set. This enables a direct comparison between simulation and measurement.
In a further step, the numerical model is calibrated to optimise the agreement between experimental data and simulation results. In particular, starting from two calibration points of the experimental data set, an adjustment of the thermal conductivities in the numerical model is carried out. While the thermal conductivity of the pore phase is held constant during the whole calibration process, thermal conductivities of the grain and contact phase are adjusted.

References
[1] Saenger et al. 2016. Analysis of high-resolution X-ray computed tomography images of Bentheim sandstone under elevated confining pressures. Geophysical Prospecting, 64(4), 848–859.

How to cite: Siegert, M., Gurris, M., and Saenger, E. H.: Numerical determination of pressure-dependent effective thermal conductivity in Berea sandstone, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-4691, https://doi.org/10.5194/egusphere-egu21-4691, 2021.

Geoelectrical methods are increasingly being used for non-invasive characterization and monitoring of permafrost sites, since the electrical properties are sensitive to the phase change of liquid to frozen water. Here, electrical resistivity tomography (ERT) is most commonly applied, using resistivity as a proxy for various quantities, such as temperature or ice content. However, it is still challenging to distinguish between air and ice in the pore space of the rock based on resistivity alone due to their similarly low electrical conductivity. Meanwhile, geoelectrical methods that utilize electrical polarization effects to characterize permafrost are also being explored. For example, the usage of the spectral induced polarization (SIP) method, in which the complex, frequency-dependent impedance is measured, can reduce ambiguities in the subsurface conduction properties, considering the SIP signature of ice. These measurements seem to be suitable for the quantification of ice content (and thus the differentiation of ice and air), and for the improved thermal characterization of alpine permafrost sites. However, to improve the interpretation of SIP measurements, it is necessary to understand in more detail the electrical conduction and polarization properties as a function of temperature, ice content, texture, and mineralogy under frozen and partially frozen conditions.

In the study presented here, electrical impedance was measured continuously using SIP in the frequency range of 10 mHz to 45 kHz on various water-saturated solid rock and loose sediment samples during controlled freeze-thaw cycles (+20°C to -40°C). These measurements were performed on rock samples from different alpine permafrost sites with different mineralogical compositions and textures. For all samples, the resistance (impedance magnitude) shows a similar temperature dependence, with increasing resistance for decreasing temperature. Also, hysteresis between freezing and thawing behavior is observed for all measurements. During freezing, a jump within the temperature-dependent resistance is observed, suggesting a lowering of the freezing point to a critical temperature where an abrupt transition from liquid water to ice occurs. During thawing, on the other hand, there is a continuous decrease in the measured resistance, suggesting a continuous thawing of the sample. The spectra of impedance phase, which is a measure for the polarization, exhibit the same qualitative, well-known temperature-dependent relaxation behaviour of ice at higher frequencies (1 kHz - 45 kHz), with variations in shape and strength for different rock texture and mineralogy. At lower frequencies (1 Hz - 1 kHz), a polarization with a weak frequency dependence is observed in the unfrozen state of the samples. We interpret this response as membrane polarization, which likewise depends on the texture as well as on the mineralogy of the respective sample. This polarization response partially vanishes during freezing. Overall, the investigated SIP spectra do not only show a dependence on texture and mineralogy, but mainly a dependence on the presence of ice in the sample as well as temperature. This indicates the possibility of a thermal characterization, as well as a determination of the ice content, of permafrost rocks using SIP.

How to cite: Limbrock, J. K., Weigand, M., and Kemna, A.: Textural and mineralogical controls on temperature dependent SIP behavior during freezing and thawing, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-14273, https://doi.org/10.5194/egusphere-egu21-14273, 2021.

Magnetic pore fabrics (MPF) are an indirect measure of the 3D pore structure. They are defined by measuring anisotropy of magnetic susceptibility after samples have been impregnated with ferrofluid. Previous studies proposed that MPFs target pores down to 10 nm. Therefore, the method complements X-ray computed tomography (XRCT) datasets, with resolution on the order of 1-10 µm. Empirical relationships exist between MPF and pore fabric, and between MPF and permeability anisotropy. This study investigates quantitative correlations between these three properties, and between measured quantities and digital-rock-model-simulations of permeability anisotropy and MPF. Samples used for this study include natural sedimentary rocks and synthetic samples. Sediments are Plio-Pleistocene calcarenite (Apulia, Italy) with ~50% porosity and complex pore structure, and Upper Marine Molasse sandstone (Belpberg, Switzerland) with 10-20% porosity and relatively homogeneous pore space properties. Synthetic samples were made from quartz sand and calcite powder in different proportions, to simulate sandstone and carbonate rocks. Samples were characterized by pycnometry, XRCT scans, MPF determination and directional permeability measurements to obtain porosity, digital rock models, MPFs and permeability anisotropy. Porosity, permeability anisotropy, and MPFs were also computed based on digital rock models derived from XRCT data, and compared to direct measurements. Permeability anisotropy and MPF are both second-order tensors, representing the average property of the entire sample. To directly relate the XRCT-derived individual pore properties to these second-order tensor quantities, a total shape ellipsoid was computed by adding the second-order tensors reflecting the best-fit ellipsoids of single pores. Once all properties were described by second-order tensors, they were correlated in terms of fabric orientation, degree and shape of anisotropy. The MPF and total shape ellipsoids are coaxial when the samples have sufficiently large pores to be resolved, and good impregnation efficiency, and as expected, total shape ellipsoids have larger anisotropy degree. Preliminary results further indicate that the permeability anisotropy is partly consistent with total shape ellipsoids and MPFs. The defined quantitative relationships facilitate the interpretation of MPF data, thus making the method more applicable to geological and fluid migration studies.

How to cite: Zhou, Y., Pugnetti, M., Foubert, A., Lanari, P., Neururer, C., and Biedermann, A.: Correlations of magnetic pore fabrics with pore fabrics derived from high-resolution X-ray computed tomography and with permeability anisotropy in sedimentary rocks and synthetic samples, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-13527, https://doi.org/10.5194/egusphere-egu21-13527, 2021.