2020, Hiller, Thomas, Dlugosch, Raphael, Müller-Petke, Mike

Surface nuclear magnetic resonance (SNMR) is a well-established technique for the hydrogeological characterization of the subsurface up to depths of about 150 m. Recently, SNMR has been adapted to investigate also the shallow unsaturated zone with small surface loop setups. Due to the decreased volume, a pre-polarization (PP) field prior to the classical spin excitation is applied to enhance the measured response signal. Depending on the strength and orientation of the applied PP-field, the enhancement can often reach several orders of magnitude in the vicinity of the PP-loop. The theoretically achievable enhancement depends on the assumption of an adiabatic, that is perfect, switch-off of the corresponding PP-field. To study the effect of imperfect switch-off, we incorporate full spin dynamics simulations into the SNMR forward modelling. The affected subsurface volume strongly depends on the chosen PP switch-off ramp and the geometry of the loop setup. Due to the imperfect switch-off, the resulting SNMR sounding curves can have significantly decreased signal amplitudes. For comparison, the signal amplitudes of either a 1 ms exponential or linear switch-off ramp are reduced by 17 and 65 per cent, respectively. Disregarding this effect would therefore yield an underestimation of the corresponding subsurface water content of similar magnitude. © 2020 The Author(s) 2020. Published by Oxford University Press on behalf of The Royal Astronomical Society.

2021, Hiller, Thomas, Costabel, Stephan, Radić, Tino, Dlugosch, Raphael, Müller‐Petke, Mike

In the past few years, small-scale (2 m) prepolarized surface nuclear magnetic resonance (SNMR) has gained increasing interest in the research community. As recent studies demonstrated, the application of a strong prepolarization field enhances the SNMR signal of coils with a footprint <1 m2 up to a level that even enables investigations in urban areas. In particular, it is expected that this noninvasive method provides the soil moisture distribution in the upper 2 m of the subsurface in the near future. However, until now all field experiments have been carried out on water reservoirs only, in an approach to test and implement this rather new technique into the field of SNMR applications. We present the first prepolarized SNMR measurement on a real soil and demonstrate the general feasibility of this technique to qualitatively and quantitatively detect soil moisture in the upper first 0.5 m. Our soil moisture measurements are validated by independent time domain reflectometry data. To complement the field experiments with numerical simulations, we adapted the underlying SNMR spin dynamics simulations and account for prepolarization switch-off effects in the forward modeling of the SNMR excitation.

2020, Smodej, Jörg, Lemmens, Laurent, Reuning, Lars, Hiller, Thomas, Klitzsch, Norbert, Claes, Steven, Kukla, Peter A.

Carbonate reservoirs form important exploration targets for the oil and gas industry in many parts of the world. This study aims to differentiate and quantify pore types and their relation to petrophysical properties in the Permo-Triassic Khuff Formation, a major carbonate reservoir in Oman. For that purpose, we have employed a number of laboratory techniques to test their applicability for the characterization of respective rock types. Consequently, a workflow has been established utilizing a combined analysis of petrographic and petrophysical methods which provide the best results for pore-system characterization. Micro-computed tomography (μCT) analysis allows a representative 3D assessment of total porosity, pore connectivity, and effective porosity of the ooid-shoal facies but it cannot resolve the full pore-size spectrum of the highly microporous mud-/wackestone facies. In order to resolve the smallest pores, combined mercury injection capillary pressure (MICP), nuclear magnetic resonance (NMR), and BIB (broad ion beam)-SEM analyses allow covering a large pore size range from millimeter to nanometer scale. Combining these techniques, three different rock types with clearly discernible pore networks can be defined. Moldic porosity in combination with intercrystalline porosity results in the highest effective porosities and permeabilities in shoal facies. In back-shoal facies, dolomitization leads to low total porosity but well-connected and heterogeneously distributed vuggy and intercrystalline pores which improves permeability. Micro- and nanopores are present in all analyzed samples but their contribution to effective porosity depends on the textural context. Our results confirm that each individual rock type requires the application of appropriate laboratory techniques. Additionally, we observe a strong correlation between the inverse formation resistivity factor and permeability suggesting that pore connectivity is the dominating factor for permeability but not pore size. In the future, this relationship should be further investigated as it could potentially be used to predict permeability from wireline resistivity measured in the flushed zone close to the borehole wall. © 2019 by the authors. Licensee MDPI, Basel, Switzerland.

2021, Costabel, Stephan, Hiller, Thomas

Geophysical nuclear magnetic resonance (NMR) applications are used to estimate pore size distributions (PSDs) of rocks and sediments. This is commonly realized by empirical calibration using information about the surface-to-volume ratio of the material. Recent research has developed joint inversion concepts for NMR relaxation data that provides the PSD with a minimum of information. The application requires the NMR signal of a sample at saturation and at least one at partial saturation and at known suction. The new inversion concept physically simulates the desaturation process as part of the forward operator. The cross-section of the model capillaries in the underlying bundle can be either circular or triangular. Our study investigates the performance of the NMR joint inversion to predict water retention function (WRF) and capillary-based hydraulic conductivity (Kcap) as functions of saturation for different sands. The angularity of the pores has no significant impact on the estimated WRF but affects the Kcap estimation significantly. Our study shows that the WRF is predicted reliably for sand samples under fast diffusion conditions. The Kcap estimations are also plausible but tend to systematic overestimation, for which we identified the tortuosity being the main reason. Because NMR relaxation data generally do not provide tortuosity information, a plausible tortuosity model remains an issue of classical calibration. Further development of the approach will thus consider tortuosity measurements (e.g., by electrical resistivity measurements and/or gradient NMR) and will consider the relaxation mechanisms outside fast diffusion conditions to enhance its applicability for coarse soils.