Three-Dimensional Modeling of the Casing Effect in Onshore Controlled-Source Electromagnetic Surveys

The presence of steel-cased wells and other infrastructure causes a significant change in the electromagnetic fields that has to be taken into consideration in modeling and interpretation of field data. A realistic and accurate simulation requires the borehole casing to be incorporated into the modeling scheme, which is numerically challenging. Due to the huge conductivity contrast between the casing and surrounding media, a spatial discretization that provides accurate results at different spatial scales ranging from millimeters to hundreds of meters is required. In this paper, we present a full 3D frequency-domain electromagnetic modeling based on a parallel finite-difference algorithm considering the casing effect and investigate its applicability on the borehole-to-surface configuration of the Hontomín CO₂ storage site. To guarantee a robust solution of linear systems with highly ill-conditioned matrices caused by huge conductivity contrasts and multiple spatial scales in the model, we employ direct sparse solvers. Different scenarios are simulated in order to study the influence of the source position, conductivity model, and the effect of the steel casing on the measured data. Several approximations of the real hollow casing that allow for a large reduction in the number of elements in the resulting meshes are studied. A good agreement between the modeled responses and the real field data demonstrates the feasibility of simulating casing effects in complex geological areas. The steel casing of the well greatly increases the amplitude of the surface electromagnetic fields and thus improves the signal-to-noise ratio and the sensitivity to deep targets.     

A land controlled‐source electromagnetic experiment using a deep vertical electric dipole: experimental settings,processing, and first data interpretation

A multichannel borehole‐to‐surface controlled‐source electromagnetic experiment was carried out at the onshore CO₂ storage site of Hontomín (Spain). The electromagnetic source consisted of a vertical electric dipole located 1.5 km deep, and the electric field was measured at the surface. The subsurface response has been obtained by calculating the transfer function between the transmitted signal and the electric field at the receiver positions. The dataset has been processed using a fast processing methodology, appropriate to be applied on controlled‐source electromagnetics (CSEM) data with a large signal‐to‐noise ratio. The dataset has been analysed in terms of data quality and repeatability errors, showing data with low experimental errors and good repeatability. We evaluate if the induction of current along the casing of the injection well can reproduce the behaviour of the experimental data.     

The influence of basement structure on the evolution of the Bicorb-Quesa Diapir (eastern Betics, Iberian Peninsula): contractive thin-skinned deformation above a pre-existing extensional basement fault

On the basis of an analysis of structural surface data in addition to previously available and new magnetotelluric data, it was possible to reconstruct the evolution of Bicorb-Quesa Diapir. This was initiated as a reactive diapir in relation to a basement fault. The reactive diapir was rejuvenated by a thin-skinned compression during the Paleogene and rose during an extensional early-middle Miocene phase. Later, in the middle Miocene, the diapir was squeezed and then, in the late Miocene, was extensionally reactivated. The current reconstruction allows us to analyse a diapir affected by a thin-skinned contractional deformation located on top of a pre-existing basement fault. Our study highlights the role played by the geometric relationship between the propagation direction of the cover deformation and the basement fault.     

Effects of galvanic distortion on magnetotelluric data over a three-dimensional regional structure

Although the galvanic distortion due to local, near-surface inhomogeneities is frequency-independent, its effect on the magnetotelluric data becomes, in a 3-D structure, frequency-dependent. Therefore, both the apparent resistivity and the phase responses are disturbed, and a correction should be carried out prior to the 3-D interpretation in order to retrieve the 3-D regional impedance tensor. In many cases, the structure is 2-D for depths corresponding to a first range of periods and 3-D for longer periods (called 2-D/3-D). For these cases, a simple method which allows us to retrieve the 3-D regional impedance tensor (except the static shift) is presented. The method proposed uses the Groom & Bailey decomposition of the distortion matrix for the short periods. Three examples are presented: two using synthetic data and one employing real data. These examples show the effect of the galvanic distortion over a regional 2-D/3-D model and the retrieval of the regional transfer functions from the distorted ones.     

Fuzzy Logic Determination of Lithologies from Well Log Data: Application to the KTB Project Data set (Germany)

Fuzzy logic has been used for lithology prediction with remarkable success. Several techniques such as fuzzy clustering or linguistic reasoning have proven to be useful for lithofacies determination. In this paper, a fuzzy inference methodology has been implemented as a MATLAB routine and applied for the first time to well log data from the German Continental Deep Drilling Program (KTB). The training of the fuzzy inference system is based on the analysis of the multi-class Matthews correlation coefficient computed for the classification matrix. For this particular data set, we have found that the best suited membership function type is the piecewise linear interpolation of the normalized histograms; that the best combination operator for obtaining the final lithology degrees of membership is the fuzzy gamma operator; and that all the available properties are relevant in the classification process. Results show that this fuzzy logic-based method is suited for rapidly and reasonably suggesting a lithology column from well log data, neatly identifying the main units and in some cases refining the classification, which can lead to a better interpretation. We have tested the trained system with synthetic data generated from property value distributions of the training data set to find that the differences in data distributions between both wells are significant enough to misdirect the inference process. However, a cross-validation analysis has revealed that, even with differences between data distributions and missing lithologies in the training data set, this fuzzy logic inference system is able to output a coherent classification.     

Forest biodiversity assessment in Peruvian Andean Montane cloud forest

Cloud forests are unusual and fragile habitats, being one of the least studied and least understood ecosystems. The tropical Andean dominion is considered one of the most significant places in the world as regards biological diversity, with a very high level of endemism. The biodiversity was analysed in an isolated remnant area of a tropical montane cloud forest known as the "Bosque de Neblina de Cuyas", in the North of the Peruvian Andean range. Composition, structure and dead wood were measured or estimated. The values obtained were compared with other cloud forests. The study revealed a high level of forest biodiversity, although the level of biodiversity differs from one area to another: in the inner areas, where human pressure is almost inexistent, the biodiversity values increase. The high species richness and the low dominance among species bear testimony to this montane cloud forest as a real enclave of biodiversity.     

Imaging the Seismic Crustal Structure of the Western Mexican Margin between 19°N and 21°N

Three thousand kilometres of multichannel (MCS) and wide-angle seismic profiles, gravity and magnetic, multibeam bathymetry and backscatter data were recorded in the offshore area of the west coast of Mexico and the Gulf of California during the spring 1996 (CORTES survey). The seismic images obtained off Puerto Vallarta, Mexico, in the Jalisco subduction zone extend from the oceanic domain up to the continental shelf, and significantly improve the knowledge of the internal crustal structure of the subduction zone between the Rivera and North American (NA) Plates. Analyzing the crustal images, we differentiate: (1) An oceanic domain with an important variation in sediment thickness ranging from 2.5 to 1 km southwards; (2) an accretionary prism comprised of highly deformed sediments, extending for a maximum width of 15 km; (3) a deformed forearc basin domain which is 25 km wide in the northern section, and is not seen towards the south where the continental slope connects directly with the accretionary prism and trench, thus suggesting a different deformational process; and (4) a continental domain consisting of a continental slope and a mid slope terrace, with a bottom simulating reflector (BSR) identified in the first second of the MCS profiles. The existence of a developed accretionary prism suggests a subduction–accretion type tectonic regime. Detailed analysis of the seismic reflection data in the oceanic domain reveals high amplitude reflections at around 6 s [two way travel time (twtt)] that clearly define the subduction plane. At 2 s (twtt) depth we identify a strong reflection which we interpret as the Moho discontinuity. We have measured a mean dip angle of 7° ± 1° at the subduction zone where the Rivera Plate begins to subduct, with the dip angle gently increasing towards the south. The oceanic crust has a mean crustal thickness of 6.0–6.5 km. We also find evidence indicating that the Rivera Plate possibly subducts at very low angles beneath the Tres Marias Islands.