Stephan Thiel | Search USGS

USGS2016

Conductivity response to intraplate deformation: Evidence for metamorphic devolatilization and crustal‐scale fluid focusing

We present two‐dimensional electrical resistivity models of two 40 km magnetotelluric (MT) profiles across the Frome Embayment to the east of the northern Flinders Ranges, South Australia. The lower crust shows low resistivity of 10 Ω m at around 30 km depth. The middle crust is dominated by resistive (>1000 Ω m) basement rocks underlying the Flinders Ranges. Adjacent to the ranges, conductive lower crust is connected to vertical zones of higher conductivity extending to just below the brittle‐ductile transition at ∼10 km depth. The conductive zones narrow in the brittle upper crust and dip at roughly 45° beneath the surface. Zones of enhanced conductivity coincide with higher strain due to topographic loading and sparse seismicity. We propose that fluids are generated through neotectonic metamorphic devolatilization. Low‐resistivity zones display areas of fluid pathways along either preexisting faults or an effect of crustal compression leading to metamorphic fluid generation. The lower crustal conductors are responding to long‐wavelength flexure‐induced strain, while the upper crustal conductors are responding to short wavelength faulting in the brittle regime. MT is a useful tool for imaging crustal strain in response to far‐field stresses in an intraplate setting and provides important constraints for geodynamic modeling and crustal rheology.

South Australia↗

USGS2015

Probabilistic 3-D time-lapse inversion of magnetotelluric data: Application to an enhanced geothermal system

Surface-based monitoring of mass transfer caused by injections and extractions in deep boreholes is crucial to maximize oil, gas and geothermal production. Inductive electromagnetic methods, such as magnetotellurics, are appealing for these applications due to their large penetration depths and sensitivity to changes in fluid conductivity and fracture connectivity. In this work, we propose a 3-D Markov chain Monte Carlo inversion of time-lapse magnetotelluric data to image mass transfer following a saline fluid injection. The inversion estimates the posterior probability density function of the resulting plume, and thereby quantifies model uncertainty. To decrease computation times, we base the parametrization on a reduced Legendre moment decomposition of the plume. A synthetic test shows that our methodology is effective when the electrical resistivity structure prior to the injection is well known. The centre of mass and spread of the plume are well retrieved.We then apply our inversion strategy to an injection experiment in an enhanced geothermal system at Paralana, South Australia, and compare it to a 3-D deterministic time-lapse inversion. The latter retrieves resistivity changes that are more shallow than the actual injection interval, whereas the probabilistic inversion retrieves plumes that are located at the correct depths and oriented in a preferential north-south direction. To explain the time-lapse data, the inversion requires unrealistically large resistivity changes with respect to the base model. We suggest that this is partly explained by unaccounted subsurface heterogeneities in the base model from which time-lapse changes are inferred.

Geophysical Journal International↗