Capability of cross-hole electrical configurations for monitoring rapid plume migration experiments

Cross-hole electrical resistivity tomography is a useful tool in geotechnical, hydrogeological or fluid/gas plume migration studies. It allows better characterization of deep subsurface structures and monitoring of the involved processes. However, due to the large amount of possible four-electrode combinations between boreholes, the choice of the most efficient ones for rapid plume migration experiments (real-time monitoring), becomes a challenge. In this work, a numerical simulation to assess the capabilities and constraints of the most common cross-hole configurations for real-time monitoring is presented. Four-electrode configurations, sensitivity, dependence on the body location and amount of data were taken into account. The analysis of anomaly detection and the symmetry of the sensitivity pattern of cross-hole configurations allowed significant reduction of the amount of data and maintaining the maximum potential resolution of each configuration for real-time monitoring. The obtained results also highlighted the benefit of using the cross-hole AB-MN configuration (with both current - or potential - electrodes located in the same borehole) combined with other configurations with complementary sensitivity pattern

Land CSEM Simulations and Experimental Test Using Metallic Casing in a Geothermal Exploration Context: Vallès Basin (NE Spain) Case Study

Controlled-source electromagnetic (CSEM) measurements are complementary data for magnetotelluric (MT) characterization although its methodology on land is not sufficiently developed and tested as in marine environments. Acquiring expertise in CSEM is crucial for surveys in places where MT cannot be performed due to high levels of cultural noise. To acquire that expertise, we perform CSEM experiments in the Vallès fault [Northeast (NE), Spain], where MT results have been satisfactory and allow us to verify the CSEM results. The Vallès basin is relevant for potential heat generation because of the presence of several geothermal anomalies and its nearby location in urban areas. In this article, we present the experimental setup for that region, a 2-D joint MT+CSEM inverse model, several 3-D CSEM simulations in the presence of metallic casing, and its comparison with real data measurements. We employ a parallel and high-order vector finite element algorithm to discretize the governing equations. By using an adapted meshing strategy, different scenarios are simulated to study the influence of the source position/direction and the conductivity model in a metallic casing presence. An excellent agreement between the simulated data and analytical/real field data demonstrates the feasibility of study metallic structures in realistic configurations. Our numerical results confirm that metallic casing strongly influences electromagnetic (EM) responses, making surface measurements more sensitive to resistivity variations near the metallic structure. It could be beneficial getting higher signal-to-noise ratios and sensitivity to deep targets. However, such a casing effect depends on the input model (e.g., conductivity contrasts, frequency, and geometry)

Comparison of stresses in 3D v. 2D geomechanical modelling of salt structures in the Tarfaya Basin, West African coast

We predict stresses and strains in the Tarfaya salt basin on the West African coast using a 3D static geomechanical model and compare the results against a simplified 2D plane-strain model. Both models are based on present-day basin geometries, are drained, and use a poroelastic description for the sediments and visco-plastic description for salt. We focus on a salt diapir, where an exploratory well has been drilled crossing a major fault. The 3D model shows a significant horizontal stress reduction in sediments at the top of the diapir, validated with measured data later obtained from the well. The 2D model predicts comparable stress reduction in sediments at the crest of the diapir. However, it shows a broader area affected by the stress reduction, overestimating its magnitude by as much as 1.5 MPa. Both models predict a similar pattern of differential displacement in sediments along both sides of the major fault, above the diapir. These displacements are the main cause of horizontal stress reduction detected at the crest of the diapir. Sensitivity analysis in both models shows that the elastic parameters of the sediments have a minimal effect on the stress-strain behaviour. In addition, the 2D sensitivity analysis concludes that the main factors controlling stress and strain changes are the geometry of the salt and the difference in rock properties between encasing sediments and salt. Overall, our study demonstrates that carefully built 2D models at the exploration stage can provide stress information and useful insights comparable to those from more complex 3D geometrie

On the detectability of Teide volcano magma chambers (Tenerife, Canary Islands) with magnetotelluric data

Tenerife has been the subject of numerous studies covering a wide range of fields. Many studies have been focused on characterising the magmatic plumbing system. Even so, a controversy still exists regarding the location and size of the current magma chambers. Several magnetotelluric (MT) surveys have been carried out in the island, but no conductivity anomalies associated with the chambers have been detected. We report the results of a set of tests conducted against the 3-D resistivity model of the island, to determine the characteristics of the detectable chambers with the MT data. The most remarkable results indicate that the MT dataset is incompatible with a large-scale mafic reservoir located at shallower depths than 8 km b.s.l. However, shallower phonolitic chambers smaller than 3 x 3 x 1 km(3) could be undetected by the existing MT sites and new data should be acquired to confirm or not their existence. This new information is essential in volcanic islands like Tenerife, since many volcanic hazards are related to the size and depth of the sources of magma. Additionally, a joint interpretation of the obtained results together with other information is summarised in a hypothetical model, allowing us to better understand the internal structure of the island

A layer stripping approach for monitoring resistivity variations using surface magnetotelluric responses

The resolution of surface-acquired magnetotelluric data is typically not sufficiently high enough in monitoring surveys to detect and quantify small resistivity variations produced within an anomalous structure at a given depth within the subsurface. To address this deficiency we present an approach, called "layer stripping", based on the analytical solution of the one-dimensional magnetotelluric problem to enhance the sensitivity of surface magnetotelluric responses to such subtle subsurface temporal variations in resistivity within e.g. reservoirs. Given a well-known geoelectrical baseline model of a reservoir site, the layer stripping approach aims to remove the effect of the upper, unchanging structures in order to simulate the time-varying magnetotelluric responses at depth. This methodology is suggested for monitoring all kinds of reservoirs, e.g. hydrocarbons, gas, geothermal, compress air storage, etc., but here we focus on CO2 geological storage. We study one-dimensional and three-dimensional resistivity variations in the reservoir layer and the feasibility of the method is appraised by evaluating the error of the approach and defining different detectability parameters. The geoelectrical baseline model of the Hontomín site (Spain) for CO2 geological storage in a deep saline aquifer is taken as our exemplar for studying the validity of the 1D assumption in a real scenario. We conclude that layer stripping could help detect resistivity variations and locate them in the space, showing potential to also sense unforeseen resistivity variations at all depths. The proposed approach constitutes an innovative contribution to take greater advantage of surface magnetotelluric data and to use the method as a cost-effective permanent monitoring technique in suitable geoelectrical scenarios

El potencial espontáneo como herramienta para la caracterización de la afección por hidrocarburos

El objeto del presente documento es poner en conocimiento de la sociedad la disponibilidad de una técnica de prospección geofísica no destructiva que permite detectar y localizar la afección por hidrocarburos en el subsuelo de instalaciones industriales.La memoria recoge la descripción del fenómeno físico y los antecedentes de dicha técnica, citándose también novedades respecto a su aplicación en instalaciones industriales.El documento puede ser de utilidad para los agentes que intervienen en actividades potencialmente contaminantes del suelo, así como de cualquier entidad que pueda acometer trabajos de caracterización de suelos con técnicas que complementen la ejecución de estas tareas

Improving the modeling of geomagnetically induced currents in Spain

Vulnerability assessments of the risk posed by geomagnetically induced currents (GICs) to power transmission grids benefit from accurate knowledge of the geomagnetic field variations at each node of the grid, the Earth's geoelectrical structures beneath them, and the topology and relative resistances of the grid elements in the precise instant of a storm. The results of previous analyses on the threat posed by GICs to the Spanish 400 kV grid are improved in this study by resorting to different strategies to progress in the three aspects identified above. First, although at midlatitude regions the source fields are rather uniform, we have investigated the effect of their spatial changes by interpolating the field from the records of several close observatories with different techniques. Second, we have performed a magnetotelluric (MT) sounding in the vicinity of one of the transformers where GICs are measured to determine the geoelectrical structure of the Earth, and we have identified the importance of estimating the MT impedance tensor when predicting GIC, especially where the effect of lateral heterogeneities is important. Finally, a sensitivity analysis to network changes has allowed us to assess the reliability of both the information about the network topology and resistances, and the assumptions made when all the details or the network status are not available. In our case, the most essential issue to improve the coincidence between model predictions and actual observations came from the use of realistic geoelectric information involving local MT measurements

Validating GIC modeling in the Spanish power grid by differential magnetometry

series of experiences and recommendations are presented concerning the derivation of geomagnetically induced currents (GIC) by use of the differential magnetometry method (DMM) under power lines. This indirect technique, intended to obtain observations to validate GIC models, is an alternative to measuring the current flow in the transformer neutrals. It is a non-intrusive and autonomous technique, as the procedure does not depend on the grid operator. In contrast, the selection of suitable sites devoid of human interferences, the need for power to supply the magnetometer, the data acquisition and transmission system, along with the choice of the appropriate instrumentation are difficulties that make not just any site suitable for installation and often require costly solutions. We focus on the methodology followed to estimate the GIC flowing in several transmission lines of the Spanish power grid with the aim of validating our GIC models, and we share our experience on the installation of the measuring points. Uncertainty inherent in the DMM is assessed, showing that noise is the main handicap, although it can be minimized with appropriate filtering. According to such experience, on some occasions only total DC currents above a significant fraction of 1 A give magnetic signatures well above the noise level, so this figure can roughly be considered as the threshold limit for detection. The low solar activity, combined with the mid-latitude condition of Spain, limited the significance of available recorded data, but we can already report and analyze the results for several minor geomagnetic storms

Integrated interpretation of geophysical data from Zagros mountain belt (Iran)

Fluid composition and distribution, the key factors determining geoelectric structure in a seismically active region, are controlled by local and regional stresses and rheological contrasts. In the central Zagros collision zone, one of the world's most seismically active mountain belt, almost coincident magnetotelluric and seismic velocity profiles are jointly interpreted to recover more accurately structural boundaries and fluid distribution within the crust. A multi-site and multi-frequency approach was used for the strike analysis of regional structure and decomposition of distortion effects on magnetotelluric data. Distortion corrected magnetotelluric data were then used for two- dimensional inversion modeling. The results image a thick conductive overburden in the southwest of the profile, high conductivities attributed to the fault zone conductors (FZCs) and an almost concave conductor extending from middle to lower crust in the central- eastern portion of the mountain belt, beneath the High Zagros (HZ). Comparison with the already available S- velocity structure, obtained by joint inversion of P-wave receiver functions and surface wave dispersion data, shows that these main conductive features are spatially correlated with a low-velocity layer representative of the sedimentary cover overlying the Arabian platform and a velocity contrast bounded by the main Zagros thrust (MZT) fault, indicating the presence of fault zone fluids. The joint interpretation of magnetotelluric inverse modeling and seismicity data also shed light on fluid generation influencing rock deformation and seismicity in this region. It suggests that beneath the HZ, deep crustal fluids generated through metamorphism may promote aseismic deformations before high stresses are buildup and cause the north- eastern part of the Zagros Fold and Thrust Belt (ZFTB) to be seismically inactive compared to its south- western part

Mesozoic structural inheritance in the Cenozoic evolution of the central Catalan Coastal Ranges (western Mediterranean): Structural and magnetotelluric analysis in the Gaià-Montmell High

The control exerted by the Mesozoic basin configuration on the Cenozoic tectonic evolution of the Catalan Coastal Ranges has been frequently recognized as a key factor to explain its present-day structure. However, details of this structural inheritance and its evolution through geological time is still under discussion. In this work we present two structural cross-sections based on fieldwork, well and magnetotelluric data in order to illustrate the structural styles and tectonic evolution of the Gaià-Montmell High. Here, the Montmell Fault not only constitutes the SW segment of one of the major Neogene faults in the Catalan Coastal Ranges (the Montmell-Vallès Fault System), but also the NW limit of a Late Jurassic-Early Cretaceous extensional basin(the Montmell-Garraf Basin), facts that denote a major role of this fault in the tectonic evolution of the area. The present-day structure of the Gaià-Montmell High resulted, therefore, from two successive episodes of inversion during the Cenozoic. The first one reactivated the Montmell Fault as compressional during the Paleogene. As a result, and among other inversion-related structures, the Gaià-El Camp Thrust developed sa major NW-directed basement footwall shortcut. Later on, the previously formed compressional structure during the Paleogene became reactivated as extensional during the Neogene. During this phase, the reactivation of the Montmell Fault looks limited and, hence, the extension is transmitted to the Baix Penedès Fault. The reactivation of the Gaià-El Camp Thrust is also manifest in the development of an array of extensional faults in the backlimb of the Carme-Cabra Anticline that corresponds to the NE-end of El Camp Fault. This episode of negative inversion developed accommodation zones between the four major faults present in the area ( Vallès-Penedès, Montmell, El Camp and Baix Penedès faults) that are characterized by the presence of relay ramps with breaching faults