ALBERT

Description: Coal seam fires are a worldwide disaster of both ecological and economic importance. Their remote detection from the ground surface or using airborne techniques is required for developing efficient strategies to extinguish them. We investigate here the use of time-domain-induced polarization to localize coal seam fires. For laboratory experiments, we first introduce a modified time-domain-induced polarization methodology to quickly acquire and invert the secondary voltage distribution mapped after the shutdown of the primary current. A set of sandbox experiments is conducted in which coal is embedded into humidified sand. Raw coal alone generates significant induced polarization anomalies, above those shown by the sand. Even higher induced polarization anomalies are detected in presence of a coal seam fire. We postulate that the higher chargeability is due to the pyrolysis, which may enhance electronic polarization or the polarization associated with the cation exchange capacity (CEC) of the material. The position of the coal seam fire is well recovered inside the tank by inverting the secondary voltages in term of a source current density distribution. We also collected field data over a recognized coal seam fire in Colorado, USA. A chargeability anomaly (~800 mV V –1 ) and a resistivity anomaly (~1 Ohm m) are observed at the position of the coal seam fire. We propose a normalized burning front index (a scaled normalized chargeability) to image and localize, without ambiguity, the position of the coal seam fire in the subsurface. The 3-D reconstructed target is located below a negative self-potential anomaly (similarly to what is observed in laboratory experiments) and a temperature anomaly recorded at a depth of 30 cm.

Description: Joint inversion that simultaneously inverts multiple geophysical data sets to recover a common Earth model is increasingly being applied to exploration problems. Petrophysical data can serve as an effective constraint to link different physical property models in such inversions. There are two challenges, among others, associated with the petrophysical approach to joint inversion. One is related to the multimodality of petrophysical data because there often exist more than one relationship between different physical properties in a region of study. The other challenge arises from the fact that petrophysical relationships have different characteristics and can exhibit point, linear, quadratic, or exponential forms in a crossplot. The fuzzy c-means (FCM) clustering technique is effective in tackling the first challenge and has been applied successfully. We focus on the second challenge in this paper and develop a joint inversion method based on variations of the FCM clustering technique. To account for the specific shapes of petrophysical relationships, we introduce several different fuzzy clustering algorithms that are capable of handling different shapes of petrophysical relationships. We present two synthetic and one field data examples and demonstrate that, by choosing appropriate distance measures for the clustering component in the joint inversion algorithm, the proposed joint inversion method provides an effective means of handling common petrophysical situations we encounter in practice. The jointly inverted models have both enhanced structural similarity and increased petrophysical correlation, and better represent the subsurface in the spatial domain and the parameter domain of physical properties.

Description: In general, the complex electrical resistivity in the subsurface is anisotropic. Despite this, algorithms for the tomographic inversion of complex resistivity data commonly assume isotropy, mainly due to the lack of anisotropic modelling and inversion schemes, potentially leading to artefacts in the inversion results in the presence of anisotropy. The development of an effective anisotropic complex resistivity inversion algorithm which utilizes the gradient information of some cost function benefits from understanding the characteristics of the problem's sensitivities, that is, the partial derivative of the impedance forward response with respect to the complex conductivities in the different spatial directions, as well as with respect to the different ratios of complex conductivities, that is, the different anisotropy ratios. We here derive expressions for these sensitivities and, based on a 2.5-D finite-element modelling algorithm, we compute and discuss sensitivity distributions as well as measurement response curves of typical surface and cross-borehole measurement configurations for 2-D subsurface anisotropic complex resistivity distributions. Depending on the electrode layout and measurement configuration, the sensitivity with respect to the conductivity in a particular direction shows a unique pattern, while for other directions sensitivity patterns are qualitatively similar. These sensitivity characteristics translate into important equivalences between impedance responses of local anisotropic and isotropic anomalies, for both magnitude and phase. Accordingly, with collinear surface arrays only the complex conductivity in the direction of the electrode layout can be unambiguously resolved, and with cross-borehole arrays only the conductivity in the vertical direction, provided an in-hole current injection is used. Nevertheless, anisotropy ratios involving these resolvable conductivity components are likewise detectable. The distinct shape of the measurement response curves, reflecting the distinct spatial patterns of the corresponding sensitivity distributions, suggest that optimized measurement configurations can be inferred for specific exploration questions involving electrical anisotropy and given electrode layouts. The gained insight into the characteristics of the sensitivity distributions of complex resistivity measurements in case of subsurface anisotropy should guide the implementation of effective anisotropic complex resistivity inversion schemes and lead to a routine use of such schemes in any resistivity and induced polarization surveys whenever subsurface electrical anisotropy could be encountered.

Description: The streaming potential phenomenon is an electrokinetic effect that occurs in porous media. It is characterized by an electrokinetic (EK) coefficient. The aim of this paper is to simulate the EK coefficient in unsaturated conditions using the Lattice Boltzmann method in a 2-D capillary channel. The multiphase flow is simulated with the model of Shan & Chen. The Poisson–Boltzmann equation is solved by implementing the model of Chai & Shi. The streaming potential response shows a non-monotonous behaviour due to the combination of the increase of charge density and decrease of flow velocity with decreasing water saturation. Using a potential of –20 mV at the air–water interface, an enhancement of a factor 5–30 of the EK coefficient, compared to the saturated state, can be observed due to the positive charge excess at this interface which is magnified by the fluid velocity away from the rock surface. This enhancement is correlated to the fractioning of the bubbles, and to the dynamic state of these bubbles, moving or entrapped in the crevices of the channel.

Description: In recent years, marine controlled source electromagnetics (CSEM) has found increasing use in hydrocarbon exploration due to its ability to detect thin resistive zones beneath the seafloor. It is the purpose of this paper to evaluate the physics of CSEM for an ocean whose electrical thickness is comparable to or much thinner than that of the overburden using the in-line configuration through examination of the elliptically polarized seafloor electric field, the time-averaged energy flow depicted by the real part of the complex Poynting vector, energy dissipation through Joule heating and the Fréchet derivatives of the seafloor field with respect to the subseafloor conductivity that is assumed to be isotropic. The deep water (ocean layer electrically much thicker than the overburden) seafloor EM response for a model containing a resistive reservoir layer has a greater amplitude and reduced phase as a function of offset compared to that for a half-space, or a stronger and faster response. For an ocean whose electrical thickness is comparable to or much smaller than that of the overburden, the electric field displays a greater amplitude and reduced phase at small offsets, shifting to a stronger amplitude and increased phase at intermediate offsets and a weaker amplitude and enhanced phase at long offsets, or a stronger and faster response that first changes to stronger and slower, and then transitions to weaker and slower. These transitions can be understood by visualizing the energy flow throughout the structure caused by the competing influences of the dipole source and guided energy flow in the reservoir layer, and the air interaction caused by coupling of the entire subseafloor resistivity structure with the sea surface. A stronger and faster response occurs when guided energy flow is dominant, while a weaker and slower response occurs when the air interaction is dominant. However, at intermediate offsets for some models, the air interaction can partially or fully reverse the direction of energy flux in the reservoir layer toward rather than away from the source, resulting in a stronger and slower response. The Fréchet derivatives are dominated by preferential sensitivity to the reservoir layer conductivity for all water depths except at high frequencies, but also display a shift with offset from the galvanic to the inductive mode in the underburden and overburden due to the interplay of guided energy flow and the air interaction. This means that the sensitivity to the horizontal conductivity is almost as strong as to the vertical component in the shallow parts of the subsurface, and in fact is stronger than the vertical sensitivity deeper down. However, the sensitivity to horizontal conductivity is still weak compared to the vertical component within thin resistive regions. The horizontal sensitivity is gradually decreased when the water becomes deep. These observations in part explain the success of shallow towed CSEM using only measurements of the in-line component of the electric field.

Description: In this work, we study seismoelectric conversions generated in the vadose zone, when this region is traversed by a pure SH wave. We assume that the soil is a 1-D partially saturated lossy porous medium and we use the van Genuchten's constitutive model to describe the water saturation profile. Correspondingly, we extend Pride's formulation to deal with partially saturated media. In order to evaluate the influence of different soil textures we perform a numerical analysis considering, among other relevant properties, the electrokinetic coupling, coseismic responses and interface responses (IRs). We propose new analytical transfer functions for the electric and magnetic field as a function of the water saturation, modifying those of Bordes et al. and Garambois & Dietrich, respectively. Further, we introduce two substantially different saturation-dependent functions into the electrokinetic (EK) coupling linking the poroelastic and the electromagnetic wave equations. The numerical results show that the electric field IRs markedly depend on the soil texture and the chosen EK coupling model, and are several orders of magnitude stronger than the electric field coseismic ones. We also found that the IRs of the water table for the silty and clayey soils are stronger than those for the sandy soils, assuming a non-monotonous saturation dependence of the EK coupling, which takes into account the charged air–water interface. These IRs have been interpreted as the result of the jump in the viscous electric current density at the water table. The amplitude of the IR is obtained using a plane SH wave, neglecting both the spherical spreading and the restriction of its origin to the first Fresnel zone, effects that could lower the predicted values. However, we made an estimation of the expected electric field IR amplitudes detectable in the field by means of the analytical transfer functions, accounting for spherical spreading of the SH seismic waves. This prediction yields a value of 15 μV m –1 , which is compatible with reported values.

Description: Surface nuclear magnetic resonance (NMR) is a unique geophysical method due to its direct sensitivity to water. A key limitation to overcome is the difficulty of making surface NMR measurements in environments with anthropogenic electromagnetic noise, particularly constant frequency sources such as powerlines. Here we present a method of removing harmonic noise by utilizing frequency domain symmetry of surface NMR signals to reconstruct portions of the spectrum corrupted by frequency-domain noise peaks. This method supplements the existing NMR processing workflow and is applicable after despiking, coherent noise cancellation, and stacking. The symmetry based correction is simple, grounded in mathematical theory describing NMR signals, does not introduce errors into the data set, and requires no prior knowledge about the harmonics. Modelling and field examples show that symmetry based noise removal reduces the effects of harmonics. In one modelling example, symmetry based noise removal improved signal-to-noise ratio in the data by 10 per cent. This improvement had noticeable effects on inversion parameters including water content and the decay constant T 2 *. Within water content profiles, aquifer boundaries and water content are more accurate after harmonics are removed. Fewer spurious water content spikes appear within aquifers, which is especially useful for resolving multilayered structures. Within T 2 * profiles, estimates are more accurate after harmonics are removed, especially in the lower half of profiles.

Description: Electromagnetic responses reflect the interaction between applied electromagnetic fields and heterogeneous geoelectrical structures. Quantifying the relationship between multiscale electrical properties and the observed electromagnetic response is therefore important for meaningful geologic interpretation. We present here examples of near-surface electromagnetic responses whose spatial fluctuations appear on all length scales, are repeatable and fractally distributed, supporting the notion of a ‘rough geology’ exhibiting multiscale hierarchical structure. Bounded by end member cases from homogenized isotropic and anisotropic media, we present numerical modelling results of the electromagnetic responses of textured and spatially correlated, stochastic geologic media, demonstrating that the electromagnetic response is a power law distribution, rather than a smooth response polluted with random, incoherent noise as commonly assumed. Our modelling results show that these electromagnetic responses due to spatially correlated geologic textures are examples of fractional Brownian motion. Furthermore, our results suggest that the fractal behaviour of the electromagnetic responses is correlated with degree of the spatial correlation, the contrasts in ground conductivity, and the preferred orientation of small-scale heterogeneity. In addition, the EM responses acquired across a fault zone comprising different lithological units and varying wavelengths of geologic heterogeneity also support our inferences from numerical modelling.

Description: Variations in self-potentials (SP) measured at surface during pumping of a heterogeneous confined fractured rock aquifer have been monitored and modelled in order to investigate capabilities and limitations of SP methods in estimating aquifer hydraulic properties. SP variations were recorded around a pumping well using an irregular grid of 31 non-polarizing Pb-PbCl 2 that were referenced to a remote electrode and connected to a commercial multiplexer and digitizer/data logger through a passive lowpass filter on each channel. The lowpass filter reduced noise by a factor of 10 compared to levels obtained using the data logger's integration-based sampling method for powerline noise suppression alone. SP signals showed a linear relationship with water levels observed in the pumping and monitoring wells over the pumping period, with an apparent electrokinetic coupling coefficient of –3.4 mV · m –1 . Following recent developments in SP methodology, variability of the SP response between different electrodes is taken as a proxy for lateral variations in hydraulic head within the aquifer and used to infer lateral variations in the aquifer's apparent transmissivity. In order to demonstrate the viability of this approach, SP is modelled numerically to determine its sensitivity to (i) lateral variations in the hydraulic conductivity of the confined aquifer and (ii) the electrical conductivity of the confining layer and conductive well casing. In all cases, SP simulated on the surface still varies linearly with hydraulic head modelled at the base on the confining layer although the apparent coupling coefficient changes to varying degrees. Using the linear relationship observed in the field, drawdown curves were inferred for each electrode location using SP variations observed over the duration of the pumping period. Transmissivity estimates, obtained by fitting the Theis model to inferred drawdown curves at all 31 electrodes, fell within a narrow range of (2.0–4.2) x 10 –3 m 2 · s –1 and were consistent with values measured in the pumping and monitoring wells. This approach will be of particular interest where monitoring wells are lacking for direct measurement, and SP on the surface can be used to quickly estimate hydraulic properties.

Description: We present a numerical algorithm for 3-D electromagnetic (EM) simulations in conducting media with general electric anisotropy. The algorithm is based on the finite-difference discretization of frequency-domain Maxwell's equations on a Lebedev grid, in which all components of the electric field are collocated but half a spatial step staggered with respect to the magnetic field components, which also are collocated. This leads to a system of linear equations that is solved using a stabilized biconjugate gradient method with a multigrid preconditioner. We validate the accuracy of the numerical results for layered and 3-D tilted transverse isotropic (TTI) earth models representing typical scenarios used in the marine controlled-source EM method. It is then demonstrated that not taking into account the full anisotropy of the conductivity tensor can lead to misleading inversion results. For synthetic data corresponding to a 3-D model with a TTI anticlinal structure, a standard vertical transverse isotropic (VTI) inversion is not able to image a resistor, while for a 3-D model with a TTI synclinal structure it produces a false resistive anomaly. However, if the VTI forward solver used in the inversion is replaced by the proposed TTI solver with perfect knowledge of the strike and dip of the dipping structures, the resulting resistivity images become consistent with the true models.