Digital rock physics: numerical prediction of pressure-dependent ultrasonic velocities using micro-CT imaging

Digital rock physics combines modern microscopic imaging with advanced numerical simulations to analyse the physical properties of rocks. Elastic-wave propagation modelling based on the microstructure images is used to estimate the effective elastic properties of the rock. The goal of this paper is to describe and understand how laboratory experiments compare with digital rock physics results using Berea sandstone. We experimentally measure pressure-dependent ultrasonic velocities and the pore size distribution. The effective elastic properties resulting from numerical simulations are based on microcomputed tomography (micro-CT) images, which are systematically stiffer than the laboratory measures. Because the tomographic images do not resolve the small-scale pore and crack network of the sample, we hypothesize that the numerical overprediction is attributable to the smallest pores and grain-to-grain contacts that are missing in the images. To reconcile the difference between numerical and experimental data, we suggest to use a grain boundary reconstruction algorithm. This allows to implement and approximate so far unresolved features in the virtual rock model. As a result, we can predict pressure-dependent effective velocity using micro-CT image

Syn-emplacement fracturing in the sandfell laccolith, eastern iceland—implications for rhyolite intrusion growth and volcanic hazards

Field work was funded by Uddeholms travel stipend (Värmlands nation, Uppsala, Sweden), Otterborgs travel stipend and the Swedish Royal Academy of Science (KVA). The research was funded by the Swedish Research Council (VR) grant 2015-03931_VR. ER is funded by the Center of Natural Hazards and Disaster Science (CNDS).Felsic magma commonly pools within shallow mushroom-shaped magmatic intrusions, so-called laccoliths or cryptodomes, which can cause both explosive eruptions and collapse of the volcanic edifice. Deformation during laccolith emplacement is primarily considered to occur in the host rock. However, shallowly emplaced laccoliths (cryptodomes) show extensive internal deformation. While deformation of magma in volcanic conduits is an important process for regulating eruptive behavior, the effects of magma deformation on intrusion emplacement remain largely unexplored. In this study, we investigate the emplacement of the 0.57km3 rhyolitic Sandfell laccolith, Iceland, which formed at a depth of 500m in a single intrusive event. By combining field measurements, 3D modeling, anisotropy of magnetic susceptibility (AMS), microstructural analysis, and FEM modeling we examine deformation in the magma to constrain its influence on intrusion emplacement. Concentric flow bands and S-C fabrics reveal contact-parallel magma flow during the initial stages of laccolith inflation. The magma flow fabric is overprinted by strain-localization bands (SLBs) and more than one third of the volume of the Sandfell laccolith displays concentric intensely fractured layers. A dominantly oblate magmatic fabric in the fractured areas and conjugate geometry of SLBs, and fractures in the fracture layers demonstrate that the magma was deformed by intrusive stresses. This implies that a large volume of magma became viscously stalled and was unable to flow during intrusion. Fine-grained groundmass and vesicle-poor rock adjacent to the fracture layers point to that the interaction between the SLBs and the flow bands at sub-solidus state caused the brittle-failure and triggered decompression degassing and crystallization, which led to rapid viscosity increase in the magma. The extent of syn-emplacement fracturing in the Sandfell laccolith further shows that strain-induced degassing limited the amount of eruptible magma by essentially solidifying the rim of the magma body. Our observations indicate that syn-emplacement changes in rheology, and the associated fracturing of intruding magma not only occur in volcanic conduits, but also play a major role in the emplacement of viscous magma intrusions in the upper kilometer of the crust.Publisher PDFPeer reviewe

Decoupling of paramagnetic and ferrimagnetic AMS development during the experimental chemical compaction of illite shale powder

Inclination shallowing of detrital remanent magnetization in sedimentary strata has solely been constrained for the mechanical processes associated with mud deposition and shallow compaction of clay-rich sediment, even though a significant part of mud diagenesis involves chemical compaction. Here we report, for the first time, on the laboratory simulation of magnetic assemblage development in a chemically compacting illite shale powder of natural origin. The experimental procedure comprised three compaction stages that, when combined, simulate the diagenesis and low-grade metamorphism of illite mud. First, the full extent of load-sensitive mechanical compaction is simulated by room temperature dry axial compression. Subsequently, temperature controlled chemical compaction is initiated by exposing the sample in two stages to amphibolite or granulite facies conditions (temperature is 490 to 750°C and confining pressure is 170 or 300 MPa) both in the absence (confining pressure only) and presence of a deformation stress field (axial compression or confined torsion). Thermodynamic equilibrium in the last two compaction stages was not reached, but illite and mica dehydroxylation initiated, thus providing a wet environment. Magnetic properties were characterized by magnetic susceptibility and its anisotropy (AMS) in both high- and low-applied field. Acquisition of isothermal remanent magnetization (IRM), stepwise three-component thermal de-magnetization of IRM and first-order reversal curves were used to characterize the remanence-bearing minerals. During the chemical compaction experiments ferrimagnetic iron-sulphides formed after reduction of magnetite and detrital pyrite in a low sulphur fugacity environment. The degree of low-field AMS is unaffected by porosity reduction from 15 to ∼1 per cent, regardless of operating conditions and compaction history. High-field paramagnetic AMS increases with compaction for all employed stress regimes and conditions, and is attributed to illite transformation to iron-bearing mica. AMS of authigenic iron-sulphide minerals remained constant during compaction indicating an independence of ferrimagnetic fabric development to chemical compaction in illite shale powder. The decoupling of paramagnetic and ferrimagnetic AMS development during chemical compaction of pelite contrasts with findings from mechanical compaction studie

Unraveling magnetic fabrics

ISSN:1437-3254ISSN:1437-326

Digital rock physics: Numerical prediction of pressure-dependent ultrasonic velocities using micro-CT imaging

ISSN:0956-540XISSN:1365-246

Magnetic anisotropy reveals Neogene tectonic overprint in highly strained carbonate mylonites from the Morcles nappe, Switzerland

The anisotropy of magnetic susceptibility (AMS) has been measured with low- and high-field methods, in deformed carbonate rocks along the Morcles nappe shear zone (Helvetic Alps). High-field measurements at room temperature and 77 K enable the separation of the ferrimagnetic, paramagnetic and diamagnetic anisotropy. The ferrimagnetic sub-fabric is generally insignificant in these rocks, contributing less than 10% to the total AMS. AMS results for both the separated diamagnetic and paramagnetic subfabrics are consistent with the regional shear movement in the late-stage formation of the Helvetic nappes, as seen in the Morcles nappe, whose inverted limb indicate shear displacement towards the northwest. The diamagnetic anisotropy correlates well quantitatively with the calculated magnetic anisotropy based on the calcite texture. There is a gradational change in the degree of anisotropy related to the strain gradient along the shear zone. A more complex magnetic fabric, resulting from partial overprinting due to displacement along the Simplon–Rhône fault, is evident at one site near the root zone of the nappe. Partial overprinting of the magnetic fabric appears to have taken place in two locations farther up the shear zone as well. This late phase deformation is associated with recent exhumation of the Mont Blanc and Belledonne external massifs and orogen parallel extension, and is reflected by the AMS. Rocks with bulk susceptibility ∼0 SI, and simple mineral compositions are ideal for low temperature high-field torque, as this method helps to enhance the paramagnetic susceptibility and anisotropy, which may otherwise be masked by the mixed magnetic contributions of the composite magnetic fabric

Decoupling of paramagnetic and ferrimagnetic AMS development during the experimental chemical compaction of illite shale powder

ISSN:0956-540XISSN:1365-246

Internal flow structures in columnar jointed basalt from Hrepphlar, Iceland II: Magnetic anisotropy and rock magnetic properties

ISSN:0258-8900ISSN:1432-081