ALBERT

Description: Rock permeability is an important parameter for the formation evaluation. In this paper, a new method with streaming current is proposed to determine the sample permeability based on the electrokinetic effects, and is proved by the experimental measurements. Corresponding to this method, we have designed an experimental setup and a test system, then performed the streaming current (potential) and electro-osmosis pressure experiments with 23 sandstone samples at 0.05 mol l –1 NaCl solution. The streaming current (potential) coefficient and electro-osmosis pressure coefficient are obtained, respectively, with the experimental data at low frequencies with AC lock-in technique. The electrokinetic permeabilities are further calculated with these coefficients. The results are consistent well with the gas permeability measured with Darcy's law, which verifies the current method for estimating rock permeability. Our measurements are also analysed and compared with previous measurements. The results indicate that our method can reflect the essence of electrokinetic effects better and simplify the electrokinetic measurements as well. In addition, we discuss the influences of experimental artefacts (core holder and confining pressure installation) on the electrokinetic data. The results show that the trough phenomenon, appeared in frequency curves of streaming current (potential) coefficients, is induced by the resonance of the core-holder/vibrator system. This is important for the design of electrokinetic setup and the analysis of low-frequency response of the electrokinetic coupling coefficients.

Description: Multiparameter full waveform inversion (FWI) applied to an elastic orthorhombic model description of the subsurface requires in theory a nine-parameter representation of each pixel of the model. Even with optimal acquisition on the Earth surface that includes large offsets, full azimuth, and multicomponent sensors, the potential for trade-off between the elastic orthorhombic parameters are large. The first step to understanding such trade-off is analysing the scattering potential of each parameter, and specifically, its scattering radiation patterns. We investigate such radiation patterns for diffraction and for scattering from a horizontal reflector considering a background isotropic model. The radiation patterns show considerable potential for trade-off between the parameters and the potentially limited resolution in their recovery. The radiation patterns of C 11 , C 22 , and C 33 are well separated so that we expect to recover these parameters with limited trade-offs. However, the resolution of their recovery represented by recovered range of model wavenumbers varies between these parameters. We can only invert for the short wavelength components (reflection) of C 33 while we can mainly invert for the long wavelength components (transmission) of the elastic coefficients C 11 and C 22 if we have large enough offsets. The elastic coefficients C 13 , C 23 , and C 12 suffer from strong trade-offs with C 55 , C 44 , and C 66 , respectively. The trade-offs between C 13 and C 55 , as well as C 23 and C 44 , can be partially mitigated if we acquire P – SV and SV – SV waves. However, to reduce the trade-offs between C 12 and C 66 , we require credible SH – SH waves. The analytical radiation patterns of the elastic constants are supported by numerical gradients of these parameters.

Description: Typically, seismic data are sparsely and irregularly sampled due to limitations in the survey environment and these cause problems for key seismic processing steps such as surface-related multiple elimination or wave-equation-based migration. Various interpolation techniques have been developed to alleviate the problems caused by sparse and irregular sampling. Among many interpolation techniques, matching pursuit interpolation is a robust tool to interpolate the regularly sampled data with large receiver separation such as crossline data in marine seismic acquisition when both pressure and particle velocity data are used. Multicomponent matching pursuit methods generally used the sinusoidal basis function, which have shown to be effective for interpolating multicomponent marine seismic data in the crossline direction. In this paper, we report the use of wavelet basis functions which further enhances the performance of matching pursuit methods for de-aliasing than sinusoidal basis functions. We also found that the range of the peak wavenumber of the wavelet is critical to the stability of the interpolation results and the de-aliasing performance and that the range should be determined based on Nyquist criteria. In addition, we reduced the computational cost by adopting the inner product of the wavelet and the input data to find the parameters of the wavelet basis function instead of using L-2 norm minimization. Using synthetic data, we illustrate that for aliased data, wavelet-based matching pursuit interpolation yields more stable results than sinusoidal function-based one when we use not only pressure data only but also both pressure and particle velocity together.

Description: Deep sea pockmarks underlain by chimney-like or pipe structures that contain methane hydrate are abundant along the Norwegian continental margin. In such hydrate provinces the interaction between hydrate formation and fluid flow has significance for benthic ecosystems and possibly climate change. The Nyegga region, situated on the western Norwegian continental slope, is characterized by an extensive pockmark field known to accommodate substantial methane gas hydrate deposits. The aim of this study is to detect and delineate both the gas hydrate and free gas reservoirs at one of Nyegga's pockmarks. In 2012, a marine controlled-source electromagnetic (CSEM) survey was performed at a pockmark in this region, where high-resolution 3-D seismic data were previously collected in 2006. 2-D CSEM inversions were computed using the data acquired by ocean bottom electrical field receivers. Our results, derived from unconstrained and seismically constrained CSEM inversions, suggest the presence of two distinctive resistivity anomalies beneath the pockmark: a shallow vertical anomaly at the underlying pipe structure, likely due to gas hydrate accumulation, and a laterally extensive anomaly attributed to a free gas zone below the base of the gas hydrate stability zone. This work contributes to a robust characterization of gas hydrate deposits within subseafloor fluid flow pipe structures.

Description: One of the main applications of magnetic field measurements in boreholes is the detection of unexploded ordnance or buried utility structures like pipes or tiebacks. Even though the advantage of fully oriented magnetic vector measurements have long been recognized and could significantly reduce costs and risks, the tools used for those purposes typically measure only the total magnetic field, the vertical and horizontal components or gradients thereof. The Göttingen Bohrloch Magnetometer uses three fibre optic gyros to record its orientation and thus enables us to compute high-quality three-component magnetic vector data regardless of borehole orientation. The measurements described in this paper were run in the scientific borehole Cuxhaven Lüdingworth 1/1A, which was drilled as a part of the ‘Coastal Aquifer Test field’ project to study the dynamics of the saltwater/freshwater interface. As the drill string got stuck during drilling of the first borehole, a second hole was drilled in the immediate vicinity. The drill string lies at a depth between 80 and 114 m at a distance of only 2.5 m southeast of the borehole used for the measurements, making it an ideal target to demonstrate the benefits of vector magnetic surveys. Although the theories to calculate magnetic fields of objects with different shapes is well established and do not need to be tested, they almost exclusively include approximations of the geometry. It is not obvious whether these approximations are suited to describe real data, or whether additional effects or refinements have to be considered. We use both a simplified monopole model and a cylinder model to fit the data and are able to determine the position of the drill string within a statistical error of approximately 10 cm. Additionally, we show that the location of the drill string could not have been determined by measurements of the total field or its horizontal and vertical component alone and that those methods would require the drilling of additional boreholes to obtain an unambiguous result.

Description: Although seismic attenuation measurements have great potential to enhance our knowledge of physical conditions and rock properties, their application is limited because robust methods for improving both the resolution and accuracy of attenuation estimates have not yet been established. We propose attenuation estimation methods for zero-offset vertical seismic profile (VSP) data by combining seismic interferometry (SI) and the modified median frequency shift (MMFS) method developed for attenuation estimation using sonic waveform data. The configuration of zero-offset VSP data is redatumed to that of the sonic logging measurement by adopting two types of SI: deconvolution interferometry and crosscorrelation interferometry (CCI). Then, we can apply the MMFS method to the redatumed VSP data. Although the amplitude information estimated from CCI is biased, we propose a correction method for this bias to correctly estimate attenuation. First, to investigate the performance both in resolution and accuracy, we apply different trace separations to synthetic data with random noise at different signal-to-noise ratio levels. Second, we estimate the influence of residual reflection events after wavefield separation on attenuation estimation. The proposed methods provide more stable attenuation estimates in comparison with the spectral ratio method because the mean-median procedure suppresses random events and characteristic features caused by residual reflection events in spectral domain. Our numerical experiments also imply that the proposed methods can estimate scattering attenuation values even if frequency components are not changed between the two receiver depths. Finally, by preliminarily applying the proposed methods to field VSP data, we demonstrate the advantages of the proposed method in the resolution and stability of attenuation estimates and these observations correlate with those of numerical tests.

Description: Capacitively Coupled Resistivity (CCR) is conventionally used to emulate DC resistivity measurements and may provide important information about the ice content of material in periglacial areas. The application of CCR theoretically enables the determination of both electrical parameters, that is, the resistivity and the electrical permittivity, by analysing magnitude and phase shift spectra. The electrical permittivity may dominate the impedance, especially in periglacial areas or regions of hydrogeological interest. However, previous theoretical work suggested that the phase shift may strongly depend on electrode height above ground, implying that electrode height must be known with great accuracy to determine electrical permittivity. Here, we demonstrate with laboratory test measurements, theoretical modelling and by analysing the Jacobian matrix of the inversion that the sensitivity towards electrode height is drastically reduced if the electrical permittivity is frequency dependent in a way that is typical for ice. For the first time, we used a novel broad-band CCR device ‘Chameleon’ for a field test located in one of the ridge galleries beneath the crest of Mount Zugspitze. A permanently ice covered bottom of a tunnel was examined. For the inversion of the measured spectra, the frequency dependence of the electrical parameters was parametrized in three different ways: A Debye Model for pure ices, a Cole–Cole Model for pure ices and a dual Cole–Cole Model including interfacial water additionally. The frequency-dependent resistivity and permittivity spectra obtained from the inversion, including low- and high-frequency limits, agree reasonably well with laboratory and field measurements reported in the literature.

Description: We evaluate the applicability and the effectiveness of GPR attribute analysis for high-resolution glacier imaging and characterization, testing this approach on 4-D GPR multifrequency data collected in a small glacier in the Eastern Alps, by repeating the acquisition along the same profiles in four different periods of the year 2013. The main objectives are to image and characterize the glacier's inner structure and to quantitatively monitor the seasonal thawing of near-surface frozen materials (snow/firn). A multiattribute approach is used to characterize the subsurface through different attribute categories, including instantaneous and textural attributes considering not only amplitude-, phase- and frequency-related attributes, but also other more complex and integrated parameters. We combine information from more than one attribute into a single image with composite displays, using overlays or mixed displays. The results demonstrate that the developed GPR attribute analysis can provide significant improvements in the discrimination of GPR signals, and obtain enhanced and more constrained data interpretations.

Description: Migration velocity analysis aims at automatically updating the large-scale components of the velocity model, called macromodel. Extended Common Image Gathers are panels used to evaluate focusing after imaging and are constructed as a function of a spatial shift introduced in the imaging condition. We investigate how transmitted waves can also be used in migration velocity analysis: instead of back-propagating the residuals associated with reflected waves, we propose to back-propagate the full wavefield. The image function, equivalent to the migrated section for reflected data, does not exhibit localized events in space along horizons but is still sensitive to the choice of the background velocity model and can thus be coupled to the same objective function defined in the image domain. In order to enhance the benefits of direct waves, we consider a cross-well configuration. Direct waves provide a large illumination between two vertical wells. Associated Common Image Gathers present different characteristics than the ones associated with reflected waves in surface acquisition. In particular, energy is spread over up to the maximum penetration depth. We invert cross-well seismic data along two lines. In the first case, the input data contain the full wavefield dominated by transmitted waves. It demonstrates the possibility to handle transmitted waves to determine the velocity model. It appears that the misfit in the data domain is largely reduced after inversion. In the second case, we use the same algorithm, but with reflected observed data only, as in a classical approach. Most of velocity updates are localized around the reflectivity, leading to an incorrect final model. This demonstrates the benefit of transmitted waves for migration velocity analysis in a cross-well configuration.

Description: An integrated approach to interpret Self-Potential (SP) anomalies based on spectral analysis and tomographic methods is presented. The Maximum Entropy Method (MEM) is used for providing accurate estimates of the depth of the anomaly source. The 2-D tomographic inversion technique, based on the underground charge occurrence probability (COP) function, is, then, used to fully characterize the anomalous body, as the MEM is not helpful in delineating the shape of the anomaly source. The proposed integrated approach is applied for the inversion of synthetic SP data generated by geometrically simple anomalous bodies, such as cylinders and inclined sheets. This numerical study has allowed the determination of mathematical relationships between zero lines of the COP distributions, the polarization angles and the positions along the profile of the causative sources, which have been of great help for interpreting the related SP anomalies. Finally, the analysis of field examples shows the high potential applicability of the proposed integrated approach for SP data inversion.

Description: Volcanic and tectonic events are the main processes involved in the generation of the oceanic crust and responsible for the seismicity associated with seafloor spreading. To monitor this activity, usually not or poorly detected by land-based seismological stations, we deployed from February 2012 to February 2013 a network of autonomous hydrophones to compare the behaviour of the ultraslow-spreading Southwest Indian ridge (SWIR) with that of the intermediate-spreading Southeast Indian ridge (SEIR). The rate of seismicity is similar for both ridges, suggesting that there is no systematic relationship between seismicity and spreading rates. The along-axis distribution of the seismic events, however, does differ, reflecting the rate dependence of accretion modes. Earthquakes are sparse and regularly spaced and scattered along the SWIR, reflecting prevailing tectonic processes. By contrast, along the SEIR, events are irregularly distributed and focus at ridge-segment ends and transforms faults, reflecting the ridge segmentation; only two swarms occurred at a segment centre and are probably caused by a magmatic event. This seismicity distribution thus looks controlled by segment-scale crustal heterogeneities along the SEIR and by regional-scale contrasting accretion processes along the SWIR, probably driven by different lithospheric and asthenospheric dynamics on either side of the Melville fracture zone. The comparison of hydroacoustic and teleseismic catalogues shows that, along these spreading ridges, the background seismicity observed in 1 yr by a hydroacoustic network is representative of the seismicity observed over two decades by land-based networks.

Description: Geothermal well GRT-1 (Rittershoffen, Alsace) was drilled in 2012. Its open-hole section (extending down to a depth of 2.6 km) penetrated fractured sandstones and granite. In 2013, the well was subjected to Thermal, Chemical and Hydraulic (TCH) stimulation, which improved the injectivity index fivefold. The goal of the study was to assess the impact of the stimulation by comparing pre- and post-stimulation well-logging (acoustic and temperature [ T ] logs) and mud-logging data. This comparison revealed modifications of almost all the natural fractures. However, not all of these fractures are associated with permeability enhancement, and the post-stimulation T logs are important for characterizing this enhancement. Chemical alteration due to mechanical erosion at the tops and bottoms of the fractures was observed in the sandstones. These zones display indications of very small new permeability after the TCH stimulation. Because a major fault zone caved extensively where it crosses the borehole, it was not imaged in the acoustic logs. However, this originally permeable zone was enhanced as demonstrated by the T logs. Based on the natural injectivity of this fault zone, hydraulic erosion and thermal microcracking of its internal quartz veins are associated with this permeability enhancement. Although local changes in the borehole wall observed in the acoustic images cannot be directly linked to the improved injectivity index, the comparison of the acoustic image logs allows for identification of fracture zones impacted by the TCH stimulation.

Description: We investigate the spatial coherence of underwater ambient noise using a yearlong time-series measured off Ascension Island. Qualitative agreement with observed cross-correlations is achieved using a simple range-dependent model, constrained by earlier, active tomographic studies in the area. In particular, the model correctly predicts the existence of two weakly dispersive normal modes in the microseism frequency range, with the group speed of one of the normal modes being smaller than the sound speed in water. The agreement justifies our interpretation of the peaks of the measured cross-correlation function of ambient noise as modal arrivals, with dispersion that is sensitive to crustal velocity structure. Our observations are consistent with Scholte to Moho head wave coupled propagation, with double mode conversion occurring due to the bathymetric variations between receivers. We thus demonstrate the feasibility of interrogating crustal properties using noise interferometry of moored hydrophone data at ranges in excess of 120 km.

Description: Recent tectonic reconstructions of the South Atlantic have partitioned the ocean basin into several segments based upon one or more proposed intraplate South American deformation zones. In several of these reconstructions, opening of the southern segment(s) by seafloor spreading prior to Aptian-Albian time is accompanied by contemporaneous strike-slip motion along an intraplate boundary extending southeastward from the Andean Cochabamba—Santa Cruz bend to the Rio Grande Fracture Zone (RGFZ). We have examined new magnetic data over the Pelotas, Santos and Campos Basins, offshore Argentina and Brazil, acquired by ION-GXT in tandem with long-offset, long record seismic reflection data, and identified seafloor spreading anomalies M4, M3, M2 and M0 (~131, ~129, ~128 and ~125 Ma). Integrating these results with our earlier work, we have been able to correlate magnetochrons M4, M3, M2 and M0 north and south of the RGFZ on the South American margin, and north and south of the Walvis Ridge on the African side. Our results are therefore inconsistent with diachronous opening models that involve substantial continental strike-slip motion north of RGFZ during M4 to M0 time. Although the ocean basin may have opened from south to north, our results indicate that seafloor spreading began north of the RGFZ earlier than previously proposed.

Description: Determining thin layer thickness is very important for reservoir characterization and CO 2 quantification. Given its high time–frequency resolution and robustness, the complex spectral decomposition method was applied on time-lapse 3-D seismic data from the Ketzin pilot site for CO 2 storage to evaluate the frequency-dependent characteristics of thin layers at the injection level. Higher temporal resolution and more stratigraphic details are seen in the all-frequency and monochromatic reflectivity amplitude sections obtained by complex spectral decomposition compared to the stacked sections. The mapped geologic discontinuities within the reservoir are consistent with the preferred orientation of CO 2 propagation. Tuning frequency mapping shows the thicknesses of the reservoir sandstone and gaseous CO 2 is consistent with the measured thickness of the sandstone unit from well logging. An attempt to discriminate between pressure effects and CO 2 saturation using the extracted tuning frequency indicates that CO 2 saturation is the main contributor to the amplitude anomaly at the Ketzin site. On the basis of determined thickness of gaseous CO 2 in the reservoir, quantitative analysis of the amount of CO 2 was performed and shows a discrepancy between the injected and calculated CO 2 mass. This may be explained by several uncertainties, like structural reservoir heterogeneity, a limited understanding of the complex subsurface conditions, error of determined tuning frequency, the presence of ambient noise and ongoing CO 2 dissolution.

Description: Comprehensive analytical solutions to 3-D axisymmetric problems for static response of multilayered thermoelastic media subjected to surface loads and containing sources are presented in this study. The solution strategy employs Laplace and Hankel transforms to the field variables. The problem is formulated in cylindrical coordinate system and in this coordinate system vector surface harmonics and generalized propagator matrix are introduced to find the solution for the problem for the behaviour of thermoelastic multilayered media subject to surface loads and containing heat sources. A high-order adaptive Gaussian quadrature method with continued fraction expansions is employed to approximate the integral solutions expressed in terms of semi-infinite Hankel-type integrals. It is the first time to apply the proposed solution method to investigate the behaviour of repository for heat-emitting high-level nuclear waste (HLW) in a geological formation where the HLW can be regarded as a decaying with time point heat source.

Description: We present a series of processes for understanding and analysing controlled-source electromagnetic (CSEM) responses for a conductive and permeable earth. To realize the CSEM response, a new 3-D CSEM forward modelling algorithm based on an edge finite element method for both electrically conductive and magnetically permeable heterogeneities is developed. The algorithm shows highly accurate results in validation tests against a semi-analytic solution for stratified earth and an integral form of the scattered field. We describe the vector behaviour of an anomalous magnetic field originating from a conductive and permeable anomaly when the loop sources are deployed over a conductive half-space. The CSEM response of the conductive and permeable anomaly is classified into three effects originating from: conductivity perturbations, permeability perturbations and the coupling of these two effects. The separated individual results and the corresponding integral equation form of the anomalous field help to better understand the physical behaviour. We confirm the characteristic features of the CSEM response from the conductive and permeable anomaly, for example, (1) the general dominance of the induction effect in the out-of-phase response accompanied by a non-negligible magnetization effect from the magnetic anomaly in a conductive half-space and (2) the dominance of near frequency-independent magnetization effects in the in-phase response at relatively low frequencies and change in ruling part of the in-phase response into the induction effect as the frequency increases. We also demonstrate the effect of coupling mode and show that its maximum contribution is limited to a few per cent level of other two modes, induction and magnetization mode, even when the heterogeneity of our model is strong. In our synthetic survey, using examples of land-based profiling surveys of low induction number and intermediate regime, we find that the effect of magnetization can be used as an indicator of the existence of magnetic material. One important point to note from this study is the importance of accurate cancelling-out or estimation of free-space responses, which can mask the magnetic responses to imaging the varying magnetic property.

Description: Extraction of fluids from subsurface reservoirs induces changes in pore pressure, leading not only to geomechanical changes, but also perturbations in seismic velocities and hence observable seismic attributes. Time-lapse seismic analysis can be used to estimate changes in subsurface hydromechanical properties and thus act as a monitoring tool for geological reservoirs. The ability to observe and quantify changes in fluid, stress and strain using seismic techniques has important implications for monitoring risk not only for petroleum applications but also for geological storage of CO 2 and nuclear waste scenarios. In this paper, we integrate hydromechanical simulation results with rock physics models and full-waveform seismic modelling to assess time-lapse seismic attribute resolution for dynamic reservoir characterization and hydromechanical model calibration. The time-lapse seismic simulations use a dynamic elastic reservoir model based on a North Sea deep reservoir undergoing large pressure changes. The time-lapse seismic traveltime shifts and time strains calculated from the modelled and processed synthetic data sets (i.e. pre-stack and post-stack data) are in a reasonable agreement with the true earth models, indicating the feasibility of using 1-D strain rock physics transform and time-lapse seismic processing methodology. Estimated vertical traveltime shifts for the overburden and the majority of the reservoir are within ±1 ms of the true earth model values, indicating that the time-lapse technique is sufficiently accurate for predicting overburden velocity changes and hence geomechanical effects. Characterization of deeper structure below the overburden becomes less accurate, where more advanced time-lapse seismic processing and migration is needed to handle the complex geometry and strong lateral induced velocity changes. Nevertheless, both migrated full-offset pre-stack and near-offset post-stack data image the general features of both the overburden and reservoir units. More importantly, the results from this study indicate that integrated seismic and hydromechanical modelling can help constrain time-lapse uncertainty and hence reduce risk due to fluid extraction and injection.

Description: A new plate motions model for the northwest Africa–North America Plate pair during the Oligocene and early Miocene is presented. The model is accompanied by a high-resolution isochron map for the central Atlantic region, resulting from a re-examination of 423 ship tracks from the NGDC data base for the area between the 15°20' FZ and the Azores triple junction. A new digital model of fracture zones for this region and a set of 309 magnetic profiles crossing the Oligocene to recent oceanic crust within the study area allowed to determine accurate finite reconstruction poles for the North America–northwest Africa conjugate plate pair between the early Miocene (Chron 6) and the early Oligocene (Chron 13). For times older than Chron 7 (~25 Ma), the finite reconstruction poles were calculated using a reliable data set coming exclusively from the region south of the Canary Islands FZ (~32°N), which allowed to test the rigidity of the northwest African oceanic lithosphere during the Oligocene–early Miocene phase of Atlas orogeny. A comparison of theoretical magnetic isochrons with observed magnetic lineations systematically shows that anomalously high spreading rates occurred in the area north of the Canary Islands FZ before Chron 7, thereby suggesting that the formation of the Atlas mountain, rather than being a localized intracontinental process, was logically linked to the central Atlantic spreading history. Thus, an independent Moroccan Plate could have existed during the Oligocene–early Miocene time interval, which included both the oceanic lithosphere north of the Canary Islands FZ and the northern Maghrebian areas of Morocco, Algeria and Tunisia. In this eventuality, the Atlas mountain belt should be reinterpreted as a giant flower structure associated with dextral transpression.

Description: Spectral induced polarization measurements are affected by temperature variations due to a variety of temperature-dependent parameters that control the complex electrical conductivity. Most important is the influence of the ion mobility, which increases with increasing temperature. It is responsible for the increase of the conductivity of the fluid in the pores with temperature and influences the electrical double layer on the mineral surface. This work is based on laboratory measurements of 13 sandstone samples from different sources with different geological and petrophysical characteristics. We measured the complex impedance in a frequency range from 0.01 to 100 Hz and a temperature range from 0 to 40 °C. The main observation is a decrease of the characteristic time (defined by the inverse of the frequency, at which the phase shift is maximum) with increasing temperature. The strength of this decrease differs from one sample to another. The temperature dependence of the phase shift magnitude cannot easily be generalized, as it depends on the particular sample. The experimental findings suggest that neglecting the influence of temperature on complex conductivity may lead to significant errors when estimating hydraulic conductivity from relaxation time. We also simulate the temperature dependence with a theoretical model of membrane polarization and review some of the model properties, with an emphasis on the temperature dependence of the parameters. The model reproduces several features characterizing the measured data, including the temperature dependence of the characteristic times. Computed tomography and microscope images of the pore structure of three samples also allow us to associate differences in the geometrical parameters used in the modelling with pore scale parameters of the actual samples.

Description: Knowledge about the pressure dependency of elastic and electrical properties is important for a variety of geophysical applications. We present a technique to invert for the stiff and compliant porosity from velocity measurements made as a function of differential pressure on saturated sandstones. A dual porosity concept is used for dry rock compressibility and a squirt model is employed for the pressure and frequency dependent elastic properties of the rocks when saturated. The total porosity obtained from inversion shows satisfactory agreement with experimental results. The electrical cementation factor was determined using the inverted porosity in combination with measured electrical conductivity. It was found that cementation factor increased exponentially with increasing differential pressure during isostatic loading. Elastic compressibility, electrical cementation factor and electrical conductivity of the saturated rocks correlate linearly with compliant porosity, and electrical cementation factor and electrical conductivity exhibit linear correlations with elastic compressibility of the saturated rocks under loading. The results show that the dual porosity concept is sufficient to explain the pressure dependency of elastic, electrical and joint elastic-electrical properties of saturated porous sandstones.

Description: A new plate motions model for the northwest Africa–North America Plate pair during the Oligocene and early Miocene is presented. The model is accompanied by a high-resolution isochron map for the central Atlantic region, resulting from a re-examination of 423 ship tracks from the NGDC data base for the area between the 15°20' FZ and the Azores triple junction. A new digital model of fracture zones for this region and a set of 309 magnetic profiles crossing the Oligocene to recent oceanic crust within the study area allowed to determine accurate finite reconstruction poles for the North America–northwest Africa conjugate plate pair between the early Miocene (Chron 6) and the early Oligocene (Chron 13). For times older than Chron 7 (~25 Ma), the finite reconstruction poles were calculated using a reliable data set coming exclusively from the region south of the Canary Islands FZ (~32°N), which allowed to test the rigidity of the northwest African oceanic lithosphere during the Oligocene–early Miocene phase of Atlas orogeny. A comparison of theoretical magnetic isochrons with observed magnetic lineations systematically shows that anomalously high spreading rates occurred in the area north of the Canary Islands FZ before Chron 7, thereby suggesting that the formation of the Atlas mountain, rather than being a localized intracontinental process, was logically linked to the central Atlantic spreading history. Thus, an independent Moroccan Plate could have existed during the Oligocene–early Miocene time interval, which included both the oceanic lithosphere north of the Canary Islands FZ and the northern Maghrebian areas of Morocco, Algeria and Tunisia. In this eventuality, the Atlas mountain belt should be reinterpreted as a giant flower structure associated with dextral transpression.

Description: Maximizing vertical resolution is a key objective in seismic data processing. Early deconvolution and spectral balancing algorithms assumed that the seismic source wavelet was temporally invariant, or stationary. In practice, seismic scattering and attenuation give rise to non-stationary seismic source wavelets. To address this issue, most conventional time-varying deconvolution wavelet shaping and spectral modelling techniques using the stationary polynomial fitting assume the wavelet to be locally stationary within a small number of overlapping analysis windows while the fitting coefficients are invariant with all the frequencies. In this paper, we show an improvement obtained by modelling smoothly varying spectra of the seismic wavelet using non-stationary polynomial fitting in the time–frequency domain. We first decompose each seismic trace using a generalized S-transform that provides a good time–frequency distribution for the estimation of the time-varying wavelet spectra. We then model the slowly varying source wavelet spectrum at each time sample by a smooth low-order polynomial. Finally, we spectrally balance the modelled wavelet to flatten the seismic response, thereby increasing vertical resolution. We calibrate the algorithm on a simple synthetic and then apply it to a 3-D land survey acquired in western China, showing the value on both vertical slices through seismic amplitude and attribute time slices. Our new algorithm significantly improves the vertical resolution of the seismic signal, while not increasing the noise.

Description: Knowledge about the pressure dependency of elastic and electrical properties is important for a variety of geophysical applications. We present a technique to invert for the stiff and compliant porosity from velocity measurements made as a function of differential pressure on saturated sandstones. A dual porosity concept is used for dry rock compressibility and a squirt model is employed for the pressure and frequency dependent elastic properties of the rocks when saturated. The total porosity obtained from inversion shows satisfactory agreement with experimental results. The electrical cementation factor was determined using the inverted porosity in combination with measured electrical conductivity. It was found that cementation factor increased exponentially with increasing differential pressure during isostatic loading. Elastic compressibility, electrical cementation factor and electrical conductivity of the saturated rocks correlate linearly with compliant porosity, and electrical cementation factor and electrical conductivity exhibit linear correlations with elastic compressibility of the saturated rocks under loading. The results show that the dual porosity concept is sufficient to explain the pressure dependency of elastic, electrical and joint elastic-electrical properties of saturated porous sandstones.

Description: In order to perform gas exploration and determine the distribution pattern of gas in the Yanchang Oil Field in the eastern part of the North Shaanxi Slope, Ordos Basin, China, gravity and magnetic survey data were systemically collated, processed and interpreted in combination with the drilling data and recent seismic data. The genesis of gravity and magnetic anomalies and the relationship between the characteristics of the gravity and magnetic fields and known gas distribution were explored in order to predict the favourable exploration targets for gas. Gravity anomalies resulted both from the lateral variation in density of the basement rock and lateral lithologic transformation in the sedimentary cover. The regional magnetic anomalies were mainly caused by the basement metamorphic rocks and the residual magnetic anomalies may reflect the amount and general location of the volcanic materials in the overlying strata. The residual gravity and magnetic anomalies generated by high-density sandstone and high content of volcanics in the gas reservoir of the upper Paleozoic distorted and deformed the anomaly curves when they were stacked onto the primary background anomaly. The gas wells were generally found to be located in the anomaly gradient zones, or the distorted part of contour lines, and the flanks of high and low anomalies, or the transitional zones between anomaly highs and lows. The characteristics of gravity and magnetic fields provide significant information that can be used for guidance when exploring the distribution of gas. Based on these characteristics, five favourable areas for gas exploration were identified; these are quasi-equally spaced like a strip extending from the southeast to the northwest.

Description: Spectral induced polarization measurements are affected by temperature variations due to a variety of temperature-dependent parameters that control the complex electrical conductivity. Most important is the influence of the ion mobility, which increases with increasing temperature. It is responsible for the increase of the conductivity of the fluid in the pores with temperature and influences the electrical double layer on the mineral surface. This work is based on laboratory measurements of 13 sandstone samples from different sources with different geological and petrophysical characteristics. We measured the complex impedance in a frequency range from 0.01 to 100 Hz and a temperature range from 0 to 40 °C. The main observation is a decrease of the characteristic time (defined by the inverse of the frequency, at which the phase shift is maximum) with increasing temperature. The strength of this decrease differs from one sample to another. The temperature dependence of the phase shift magnitude cannot easily be generalized, as it depends on the particular sample. The experimental findings suggest that neglecting the influence of temperature on complex conductivity may lead to significant errors when estimating hydraulic conductivity from relaxation time. We also simulate the temperature dependence with a theoretical model of membrane polarization and review some of the model properties, with an emphasis on the temperature dependence of the parameters. The model reproduces several features characterizing the measured data, including the temperature dependence of the characteristic times. Computed tomography and microscope images of the pore structure of three samples also allow us to associate differences in the geometrical parameters used in the modelling with pore scale parameters of the actual samples.

Description: The streaming potential phenomenon is produced by the flow of an electrolyte in a porous medium and is used for geophysical prospecting. It is quantified through an electrokinetic (EK) coefficient. The dependence of the EK coefficient on the conductivity of the electrolyte is described by the Helmholtz–Smoluchowski (HS) equation. This equation provides successful forecasts of the EK coefficient in the standard range of concentration. However, experimental measurements show deviations to this equation at extreme low and extreme high salinities. The aim of this study is to model the EK coefficient using Lattice Boltzmann simulations in a 2-D capillary channel, with a view to understanding these deviations. The effect of the constitutive parameters of the HS equation such as the permittivity and the viscosity is discussed. The validity of the HS equation using strong potentials is assessed. Finally, a model of bulk fluid conductivity is derived. This model allows to take into account the change of local ionic distribution in the vicinity of the mineral. It appears to have a significant impact on the derivation of potentials at low salinities and in the presence of polyvalent counterions.

Description: Extraction of fluids from subsurface reservoirs induces changes in pore pressure, leading not only to geomechanical changes, but also perturbations in seismic velocities and hence observable seismic attributes. Time-lapse seismic analysis can be used to estimate changes in subsurface hydromechanical properties and thus act as a monitoring tool for geological reservoirs. The ability to observe and quantify changes in fluid, stress and strain using seismic techniques has important implications for monitoring risk not only for petroleum applications but also for geological storage of CO 2 and nuclear waste scenarios. In this paper, we integrate hydromechanical simulation results with rock physics models and full-waveform seismic modelling to assess time-lapse seismic attribute resolution for dynamic reservoir characterization and hydromechanical model calibration. The time-lapse seismic simulations use a dynamic elastic reservoir model based on a North Sea deep reservoir undergoing large pressure changes. The time-lapse seismic traveltime shifts and time strains calculated from the modelled and processed synthetic data sets (i.e. pre-stack and post-stack data) are in a reasonable agreement with the true earth models, indicating the feasibility of using 1-D strain rock physics transform and time-lapse seismic processing methodology. Estimated vertical traveltime shifts for the overburden and the majority of the reservoir are within ±1 ms of the true earth model values, indicating that the time-lapse technique is sufficiently accurate for predicting overburden velocity changes and hence geomechanical effects. Characterization of deeper structure below the overburden becomes less accurate, where more advanced time-lapse seismic processing and migration is needed to handle the complex geometry and strong lateral induced velocity changes. Nevertheless, both migrated full-offset pre-stack and near-offset post-stack data image the general features of both the overburden and reservoir units. More importantly, the results from this study indicate that integrated seismic and hydromechanical modelling can help constrain time-lapse uncertainty and hence reduce risk due to fluid extraction and injection.

Description: We study the accuracy and numerical stability of three eigenvector sets for modelling the coupled poroelastic and electromagnetic layered-Earth response. We use a known eigenvector set, its flux-normalized version and a newly derived flux-normalized set. The new set is chosen such that the system is properly uncoupled when the coupling between the poroelastic and electromagnetic fields vanishes. We carry out two different numerical stability tests: the first test focuses on the internal system, eigenvector and eigenvalue consistency; the second test investigates the stability and preciseness of the flux-normalized systems by looking at identity relations. We find that the known set shows the largest deviation for both tests, whereas the new set performs best. In two additional numerical modelling experiments, these numerical inaccuracies are shown to generate numerical noise levels comparable to small signals, such as signals coming from the important interface conversion responses, especially when the coupling coefficient is small. When coupling vanishes completely, the known set does not produce proper results. The new set produces numerically stable and accurate results in all situations. We therefore strongly recommend to use this newly derived set for future layered-Earth seismo-electromagnetic modelling experiments.

Description: The Quaternary Añavieja-Dévanos tufa system is located in the northern sector of the Iberian Chain. It has been previously tackled by means sedimentological studies focused on the available outcrops and some boreholes. They have permitted the proposal of a sedimentary scenario that fits with a pool-barrage fluvial tufa model. However a better knowledge of the characteristics and internal distribution of the usually non-outcropping pool deposits as well as of its relationship with barrage deposits has not been evaluated in detail yet. Palaeoenvironmental studies on tufas are usually biased because tufas are commonly delicate facies exposed to intense erosion during water level fall stages; for this reason outcrops are usually scarce and very often coincide with the most cemented barrage deposits. In order to analyse the internal characteristics of the tufa deposits under study, but also the lateral correlation among different facies, ground penetrating radar (GPR) has been employed both for the evaluation of its applicability in such kind of environments and to improve, if possible, the sedimentary model using geophysical data in sectors without outcrops. A GPR survey including different antennas ranging from 50 to 500 MHz along different sectors and its comparison with natural outcrops has been carried out. GPR results have permitted to deduce clear differences between pool and barrage deposits and to recognise its internal structure and geometrical relationships. The survey also permitted an approach to different scales of heterogeneities in the radarfacies evaluation by using distinct antennas and therefore, reaching different resolutions and penetrations. The resulting integration from different antennas allows three different attenuant and eight reflective radarfacies to be defined permitting a better approach to the real extension of the pool areas. These results have permitted to decipher the horizontal and vertical facies changes and the identification of a scarcer development of pool deposits than expected in the studied system.

Description: Two novel transient controlled source electromagnetic methods called circular electrical dipole (CED) and differential electrical dipole (DED) are theoretically analysed for applications in shallow marine environments. 1-D and 3-D time-domain modelling studies are used to investigate the detectability and applicability of the methods when investigating resistive layers/targets representing hydrocarbon-saturated formations. The results are compared to the conventional time-domain horizontal electrical dipole (HED) and vertical electrical dipole (VED) sources. The applied theoretical modelling studies demonstrate that CED and DED have higher signal detectability towards resistive targets compared to TD-CSEM, but demonstrate significantly poorer signal amplitudes. Future CED/DED applications will have to solve this issue prior to measuring. Furthermore, the two novel methods have very similar detectability characteristics towards 3-D resistive targets embedded in marine sediments as VED while being less susceptible towards non-verticality. Due to the complex transmitter design of CED/DED the systems are prone to geometrical errors. Modelling studies show that even small transmitter inaccuracies have strong effects on the signal characteristics of CED making an actual marine application difficult at the present time. In contrast, the DED signal is less affected by geometrical errors in comparison to CED and may therefore be more adequate for marine applications.

Description: An L-configured, three-component short period seismic array was deployed on the Ross Ice Shelf, Antarctica during November 2014. Polarization analysis of ambient noise data from these stations shows linearly polarized waves for frequency bands between 0.2 and 2 Hz. A spectral peak at about 1.6 Hz is interpreted as the resonance frequency of the water column and is used to estimate the water layer thickness below the ice shelf. The frequency band from 4 to 18 Hz is dominated by Rayleigh and Love waves propagating from the north that, based on daily temporal variations, we conclude were generated by field camp activity. Frequency–slowness plots were calculated using beamforming. Resulting Love and Rayleigh wave dispersion curves were inverted for the shear wave velocity profile within the firn and ice to ~150 m depth. The derived density profile allows estimation of the pore close-off depth and the firn–air content thickness. Separate inversions of Rayleigh and Love wave dispersion curves give different shear wave velocity profiles within the firn. We attribute this difference to an effective anisotropy due to fine layering. The layered structure of firn, ice, water and the seafloor results in a characteristic dispersion curve below 7 Hz. Forward modelling the observed Rayleigh wave dispersion curves using representative firn, ice, water and sediment structures indicates that Rayleigh waves are observed when wavelengths are long enough to span the distance from the ice shelf surface to the seafloor. The forward modelling shows that analysis of seismic data from an ice shelf provides the possibility of resolving ice shelf thickness, water column thickness and the physical properties of the ice shelf and underlying seafloor using passive-source seismic data.

Description: We present an algorithm to recover the Bayesian posterior model probability density function of subsurface elastic parameters, as required by the full pressure field recorded at an ocean bottom cable due to an impulsive seismic source. Both the data noise and source wavelet are estimated by our algorithm, resulting in robust estimates of subsurface velocity and density. In contrast to purely gradient based approaches, our method avoids model regularization entirely and produces an ensemble of models that can be visualized and queried to provide meaningful information about the sensitivity of the data to the subsurface, and the level of resolution of model parameters. Our algorithm is trans-dimensional and performs model selection, sampling over a wide range of model parametrizations. We follow a frequency domain approach and derive the corresponding likelihood in the frequency domain. We present first a synthetic example of a reservoir at 2 km depth with minimal acoustic impedance contrast, which is difficult to study with conventional seismic amplitude versus offset changes. Finally, we apply our methodology to survey data collected over the Alba field in the North Sea, an area which is known to show very little lateral heterogeneity but nevertheless presents challenges for conventional post migration seismic amplitude versus offset analysis.

Description: A novel application of seismic interferometry (SI) and Marchenko imaging using both surface and borehole data is presented. A series of redatuming schemes is proposed to combine both data sets for robust deep local imaging in the presence of velocity uncertainties. The redatuming schemes create a virtual acquisition geometry where both sources and receivers lie at the horizontal borehole level, thus only a local velocity model near the borehole is needed for imaging, and erroneous velocities in the shallow area have no effect on imaging around the borehole level. By joining the advantages of SI and Marchenko imaging, a macrovelocity model is no longer required and the proposed schemes use only single-component data. Furthermore, the schemes result in a set of virtual data that have fewer spurious events and internal multiples than previous virtual source redatuming methods. Two numerical examples are shown to illustrate the workflow and to demonstrate the benefits of the method. One is a synthetic model and the other is a realistic model of a field in the North Sea. In both tests, improved local images near the boreholes are obtained using the redatumed data without accurate velocities, because the redatumed data are close to the target.

Description: The effect of the selection of different nuclear magnetic resonance (NMR) relaxation times for permeability estimation is investigated for a set of fully brine-saturated rocks acquired from Cretaceous carbonate reservoirs in the North Sea and Middle East. Estimators that are obtained from the relaxation times based on the Pythagorean means are compared with estimators that are obtained from the relaxation times based on the concept of a cumulative saturation cut-off. Select portions of the longitudinal ( T 1 ) and transverse ( T 2 ) relaxation-time distributions are systematically evaluated by applying various cut-offs, analogous to the Winland-Pittman approach for mercury injection capillary pressure (MICP) curves. Finally, different approaches to matching the NMR and MICP distributions using different mean-based scaling factors are validated based on the performance of the related size-scaled estimators. The good results that were obtained demonstrate possible alternatives to the commonly adopted logarithmic mean estimator and reinforce the importance of NMR-MICP integration to improving carbonate permeability estimates.

Description: Static time-delay corrections are frequency independent and ignore velocity variations away from the assumed vertical ray path through the subsurface. There is therefore a clear potential for improvement if the finite frequency nature of wave propagation can be properly accounted for. Such a method is presented here based on the Born approximation, the assumption of surface consistency and the misfit of instantaneous phase. The concept of instantaneous phase lends itself very well for sweep-like signals, hence these are the focus of this study. Analytical sensitivity kernels are derived that accurately predict frequency-dependent phase shifts due to P -wave anomalies in the near surface. They are quick to compute and robust near the source and receivers. An additional correction is presented that re-introduces the nonlinear relation between model perturbation and phase delay, which becomes relevant for stronger velocity anomalies. The phase shift as function of frequency is a slowly varying signal, its computation therefore does not require fine sampling even for broad-band sweeps. The kernels reveal interesting features of the sensitivity of seismic arrivals to the near surface: small anomalies can have a relative large impact resulting from the medium field term that is dominant near the source and receivers. Furthermore, even simple velocity anomalies can produce a distinct frequency-dependent phase behaviour. Unlike statics, the predicted phase corrections are smooth in space. Verification with spectral element simulations shows an excellent match for the predicted phase shifts over the entire seismic frequency band. Applying the phase shift to the reference sweep corrects for wavelet distortion, making the technique akin to surface consistent deconvolution, even though no division in the spectral domain is involved. As long as multiple scattering is mild, surface consistent finite frequency phase corrections outperform traditional statics for moderately large velocity contrasts, also in the case of a free surface with topography. The kernels could also be used to update the velocity model, if one has knowledge of (wave propagation related only) surface consistent phase operators.

Description: The seislet transform has been demonstrated to have a better compression performance for seismic data compared with other well-known sparsity promoting transforms, thus it can be used to remove random noise by simply applying a thresholding operator in the seislet domain. Since the seislet transform compresses the seismic data along the local structures, the seislet thresholding can be viewed as a simple structural filtering approach. Because of the dependence on a precise local slope estimation, the seislet transform usually suffers from low compression ratio and high reconstruction error for seismic profiles that have dip conflicts. In order to remove the limitation of seislet thresholding in dealing with conflicting-dip data, I propose a dip-separated filtering strategy. In this method, I first use an adaptive empirical mode decomposition based dip filter to separate the seismic data into several dip bands (5 or 6). Next, I apply seislet thresholding to each separated dip component to remove random noise. Then I combine all the denoised components to form the final denoised data. Compared with other dip filters, the empirical mode decomposition based dip filter is data-adaptive. One only needs to specify the number of dip components to be separated. Both complicated synthetic and field data examples show superior performance of my proposed approach than the traditional alternatives. The dip-separated structural filtering is not limited to seislet thresholding, and can also be extended to all those methods that require slope information.

Description: Noise interferometry is the process by which approximations to acoustic Green's functions, which describe sound propagation between two locations, are estimated by cross-correlating time series of ambient noise measured at those locations. Noise-interferometry-based approximations to Green's functions can be used as the basis for a variety of inversion algorithms, thereby providing a purely passive alternative to active-source ocean acoustic remote sensing. In this paper we give an overview of results from noise interferometry experiments conducted in the Florida Straits at 100 m depth in December 2012, and at 600 m depth in September/October 2013. Under good conditions for noise interferometry, estimates of cross-correlation functions are shown to allow one to perform advanced phase-coherent signal processing techniques to perform waveform inversions, estimate currents by exploiting non-reciprocity, perform time-reversal/back-propagation calculations and investigate modal dispersion using time-warping techniques. Conditions which are favourable for noise interferometry are identified and discussed.

Description: In this paper, we consider a new spectral finite volume method (FVM) for the elastic wave equations. Our new FVM is based on a piecewise constant approximation on a fine mesh and a high-order polynomial reconstruction on a coarser mesh. Our new method is constructed based on two existing techniques, the high-order FVM and the spectral FVM. In fact, we will construct a new method to take advantage of both methods. More precisely, our method has two distinctive features. The first one is that the local polynomial reconstructions are performed on the coarse triangles and the reconstruction matrices for all the coarse triangles are the same. This fact enhances the parallelization of our algorithm. We will present a parallel implementation of our method and show excellent efficiency results. The second one is that, by using a suitable number of finer triangles with a coarse triangle, we obtain an overdetermined reconstruction system, which can enhance the robustness of the reconstruction process. To derive our scheme, standard finite volume technique is applied to each fine triangle, and the high-order reconstructed polynomials, computed on coarse triangles, are used to compute numerical fluxes. We will present numerical results to show the performance of our method. Our method is presented for 2-D problems, but the same methodology can be applied to 3-D.

Description: At the heart of the full waveform inversion (FWI) implementation is wavefield extrapolation, and specifically its accuracy and cost. To obtain accurate, dispersion free wavefields, the extrapolation for modelling is often expensive. Combining an efficient extrapolation with a novel gradient preconditioning can render an FWI implementation that efficiently converges to an accurate model. We, specifically, recast the extrapolation part of the inversion in terms of its spectral components for both data and gradient calculation. This admits dispersion free wavefields even at large extrapolation time steps, which improves the efficiency of the inversion. An alternative spectral representation of the depth axis in terms of sine functions allows us to impose a free surface boundary condition, which reflects our medium boundaries more accurately. Using a newly derived perfectly matched layer formulation for this spectral implementation, we can define a finite model with absorbing boundaries. In order to reduce the nonlinearity in FWI, we propose a multiscale conditioning of the objective function through combining the different directional components of the gradient to optimally update the velocity. Through solving a simple optimization problem, it specifically admits the smoothest approximate update while guaranteeing its ascending direction. An application to the Marmousi model demonstrates the capability of the proposed approach and justifies our assertions with respect to cost and convergence.

Description: Spectral induced polarization (SIP) data are commonly analysed using phenomenological models. Among these models the Cole–Cole (CC) model is the most popular choice to describe the strength and frequency dependence of distinct polarization peaks in the data. More flexibility regarding the shape of the spectrum is provided by decomposition schemes. Here the spectral response is decomposed into individual responses of a chosen elementary relaxation model, mathematically acting as kernel in the involved integral, based on a broad range of relaxation times. A frequently used kernel function is the Debye model, but also the CC model with some other a priorly specified frequency dispersion (e.g. Warburg model) has been proposed as kernel in the decomposition. The different decomposition approaches in use, also including conductivity and resistivity formulations, pose the question to which degree the integral spectral parameters typically derived from the obtained relaxation time distribution are biased by the approach itself. Based on synthetic SIP data sampled from an ideal CC response, we here investigate how the two most important integral output parameters deviate from the corresponding CC input parameters. We find that the total chargeability may be underestimated by up to 80 per cent and the mean relaxation time may be off by up to three orders of magnitude relative to the original values, depending on the frequency dispersion of the analysed spectrum and the proximity of its peak to the frequency range limits considered in the decomposition. We conclude that a quantitative comparison of SIP parameters across different studies, or the adoption of parameter relationships from other studies, for example when transferring laboratory results to the field, is only possible on the basis of a consistent spectral analysis procedure. This is particularly important when comparing effective CC parameters with spectral parameters derived from decomposition results.

Description: Viability for the development of an engineered geothermal system (EGS) in the oilsands region near Fort McMurray, Alberta, is investigated by studying the structure of the Precambrian basement rocks with magnetotellurics (MT). MT data were collected at 94 broad-band stations on two east–west profiles. Apparent resistivity and phase data showed little variation along each profile. The short period MT data detected a 1-D resistivity structure that could be identified as the shallow sedimentary basin underlain by crystalline basement rocks to a depth of 4–5 km. At lower frequencies a strong directional dependence, large phase splits, and regions of out-of-quadrant (OOQ) phase were detected. 2-D isotropic inversions of these data failed to produce a realistic resistivity model. A detailed dimensionality analysis found links between large phase tensor skews (~15°), azimuths, OOQ phases and tensor decomposition strike angles at periods greater than 1 s. Low magnitude induction vectors, as well as uniformity of phase splits and phase tensor character between the northern and southern profiles imply that a 3-D analysis is not necessary or appropriate. Therefore, 2-D anisotropic forward modelling was used to generate a resistivity model to interpret the MT data. The preferred model was based on geological observations of outcropping anisotropic mylonitic basement rocks of the Charles Lake shear zone, 150 km to the north, linked to the study area by aeromagnetic and core sample data. This model fits all four impedance tensor elements with an rms misfit of 2.82 on the southern profile, and 3.3 on the northern. The conductive phase causing the anisotropy is interpreted to be interconnected graphite films within the metamorphic basement rocks. Characterizing the anisotropy is important for understanding how artificial fractures, necessary for EGS development, would form. Features of MT data commonly interpreted to be 3-D (e.g. out of OOQ phase and large phase tensor skew) are shown to be interpretable with this 2-D anisotropic model.

Description: Deep geological repositories, isolated from the geosphere by an engineered bentonite barrier, are currently considered the safest solution for high-level radioactive waste (HLRW) disposal. As the physical conditions and properties of the bentonite barrier are anticipated to change with time, seismic tomography was suggested as a viable technique to monitor the physical state and integrity of the barrier and to timely detect any unforeseen failure. To do so, the seismic monitoring system needs to be optimized, and this can be achieved by conducting numerical simulations of wave propagation in the repository geometry. Previous studies treated bentonite as an elastic medium, whereas recent experimental investigations indicate its pronounced viscoelastic behaviour. The aims of this contribution are (i) to numerically estimate the effective attenuation of bentonite as a function of temperature T and water content W c , so that synthetic data can accurately reproduce experimental traces and (ii) assess the feasibility and limitation of the HLRW repository monitoring by simulating the propagation of sonic waves in a realistic repository geometry. A finite difference method was utilized to simulate the wave propagation in experimental and repository setups. First, the input of the viscoelastic model was varied to achieve a match between experimental and numerical traces. The routine was repeated for several values of W c and T , so that quality factors Q p ( W c , T ) and Q s ( W c , T ) were obtained. Then, the full-scale monitoring procedure was simulated for six scenarios, representing the evolution of bentonite's physical state. The estimated Q p and Q s exhibited a minimum at W c = 20 per cent and higher sensitivity to W c , rather than T , suggesting that pronounced inelasticity of the clay has to be taken into account in geophysical modelling and analysis. The repository-model traces confirm that active seismic monitoring is, in principle, capable of depicting physical changes in the bentonite barrier. However, the locations of sources and receivers relative to the tunnel need to be optimized to achieve best sensitivity to the changes. It was also observed that first arrivals do not necessarily represent the most informative part of the signal; thus the time windows of the analysis need to be selected, accounting for the geometry of the acquisition system. This contribution should be considered as a next step towards the development of a numerical approach to aid in the design and optimization of a non-intrusive monitoring system in HLRW repositories.

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.

Description: We develop a procedure to invert time domain induced polarization (IP) data for inductive sources. Our approach is based upon the inversion methodology in conventional electrical IP (EIP), which uses a sensitivity function that is independent of time. However, significant modifications are required for inductive source IP (ISIP) because electric fields in the ground do not achieve a steady state. The time-history for these fields needs to be evaluated and then used to define approximate IP currents. The resultant data, either a magnetic field or its derivative, are evaluated through the Biot-Savart law. This forms the desired linear relationship between data and pseudo-chargeability. Our inversion procedure has three steps: (1) Obtain a 3-D background conductivity model. We advocate, where possible, that this be obtained by inverting early-time data that do not suffer significantly from IP effects. (2) Decouple IP responses embedded in the observations by forward modelling the TEM data due to a background conductivity and subtracting these from the observations. (3) Use the linearized sensitivity function to invert data at each time channel and recover pseudo-chargeability. Post-interpretation of the recovered pseudo-chargeabilities at multiple times allows recovery of intrinsic Cole-Cole parameters such as time constant and chargeability. The procedure is applicable to all inductive source survey geometries but we focus upon airborne time domain EM (ATEM) data with a coincident-loop configuration because of the distinctive negative IP signal that is observed over a chargeable body. Several assumptions are adopted to generate our linearized modelling but we systematically test the capability and accuracy of the linearization for ISIP responses arising from different conductivity structures. On test examples we show: (1) our decoupling procedure enhances the ability to extract information about existence and location of chargeable targets directly from the data maps; (2) the horizontal location of a target body can be well recovered through inversion; (3) the overall geometry of a target body might be recovered but for ATEM data a depth weighting is required in the inversion; (4) we can recover estimates of intrinsic and that may be useful for distinguishing between two chargeable targets.

Description: The extraction of spectral information in the inversion process of time-domain (TD) induced polarization (IP) data is changing the use of the TDIP method. Data interpretation is evolving from a qualitative description of the subsurface, able only to discriminate the presence of contrasts in chargeability parameters, towards a quantitative analysis of the investigated media, which allows for detailed soil- and rock-type characterization. Two major limitations restrict the extraction of the spectral information of TDIP data in the field: (i) the difficulty of acquiring reliable early-time measurements in the millisecond range and (ii) the self-potential background drift in the measured potentials distorting the shape of the late-time IP responses, in the second range. Recent developments in TDIP acquisition equipment have given access to full-waveform recordings of measured potentials and transmitted current, opening for a breakthrough in data processing. For measuring at early times, we developed a new method for removing the significant noise from power lines contained in the data through a model-based approach, localizing the fundamental frequency of the power-line signal in the full-waveform IP recordings. By this, we cancel both the fundamental signal and its harmonics. Furthermore, an efficient processing scheme for identifying and removing spikes in TDIP data was developed. The noise cancellation and the de-spiking allow the use of earlier and narrower gates, down to a few milliseconds after the current turn-off. In addition, tapered windows are used in the final gating of IP data, allowing the use of wider and overlapping gates for higher noise suppression with minimal distortion of the signal. For measuring at late times, we have developed an algorithm for removal of the self-potential drift. Usually constant or linear drift-removal algorithms are used, but these algorithms often fail in removing the background potentials present when the electrodes used for potential readings are previously used for current injection, also for simple contact resistance measurements. We developed a drift-removal scheme that models the polarization effect and efficiently allows for preserving the shape of the IP responses at late times. Uncertainty estimates are essential in the inversion of IP data. Therefore, in the final step of the data processing, we estimate the data standard deviation based on the data variability within the IP gates and the misfit of the background drift removal Overall, the removal of harmonic noise, spikes, self-potential drift, tapered windowing and the uncertainty estimation allows for doubling the usable range of TDIP data to almost four decades in time (corresponding to four decades in frequency), which will significantly advance the applicability of the IP method.

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: 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: 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: Joint inversion and time-lapse inversion techniques of geophysical data are often implemented in an attempt to improve imaging of complex subsurface structures and dynamic processes by minimizing negative effects of random and uncorrelated spatial and temporal noise in the data. We focus on the structural cross-gradient (SCG) approach (enforcing recovered models to exhibit similar spatial structures) in combination with time-lapse inversion constraints applied to surface-based electrical resistivity and seismic traveltime refraction data. The combination of both techniques is justified by the underlying petrophysical models. We investigate the benefits and trade-offs of SCG and time-lapse constraints. Using a synthetic case study, we show that a combined joint time-lapse inversion approach provides an overall improvement in final recovered models. Additionally, we introduce a new approach to reweighting SCG constraints based on an iteratively updated normalized ratio of model sensitivity distributions at each time-step. We refer to the new technique as the Automatic Joint Constraints (AJC) approach. The relevance of the new joint time-lapse inversion process is demonstrated on the synthetic example. Then, these approaches are applied to real time-lapse monitoring field data collected during a quarter-scale earthen embankment induced-piping failure test. The use of time-lapse joint inversion is justified by the fact that a change of porosity drives concomitant changes in seismic velocities (through its effect on the bulk and shear moduli) and resistivities (through its influence upon the formation factor). Combined with the definition of attributes (i.e. specific characteristics) of the evolving target associated with piping, our approach allows localizing the position of the preferential flow path associated with internal erosion. This is not the case using other approaches.

Description: Geophysical inversion should ideally produce geologically realistic subsurface models that explain the available data. Multiple-point statistics is a geostatistical approach to construct subsurface models that are consistent with site-specific data, but also display the same type of patterns as those found in a training image. The training image can be seen as a conceptual model of the subsurface and is used as a non-parametric model of spatial variability. Inversion based on multiple-point statistics is challenging due to high nonlinearity and time-consuming geostatistical resimulation steps that are needed to create new model proposals. We propose an entirely new model proposal mechanism for geophysical inversion that is inspired by texture synthesis in computer vision. Instead of resimulating pixels based on higher-order patterns in the training image, we identify a suitable patch of the training image that replace a corresponding patch in the current model without breaking the patterns found in the training image, that is, remaining consistent with the given prior. We consider three cross-hole ground-penetrating radar examples in which the new model proposal mechanism is employed within an extended Metropolis Markov chain Monte Carlo (MCMC) inversion. The model proposal step is about 40 times faster than state-of-the-art multiple-point statistics resimulation techniques, the number of necessary MCMC steps is lower and the quality of the final model realizations is of similar quality. The model proposal mechanism is presently limited to 2-D fields, but the method is general and can be applied to a wide range of subsurface settings and geophysical data types.

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: The oceanic lithosphere in the Bay of Bengal (BOB) formed 80–120 Ma following the breakup of eastern Gondwanaland. Since its formation, it has been affected by the emplacement of two long N-S trending linear aseismic ridges (85°E and Ninetyeast) and by the loading of ca . 20-km of sediments of the Bengal Fan. Here, we present the results of a combined spatial and spectral domain analysis of residual geoid, bathymetry and gravity data constrained by seismic reflection and refraction data. Self-consistent geoid and gravity modelling defined by temperature-dependent mantle densities along a N–S transect in the BOB region revealed that the depth to the lithosphere–asthenosphere boundary (LAB) deepens steeply from 77 km in the south to 127 km in north, with the greater thickness being anomalously thick compared to the lithosphere of similar-age beneath the Pacific Ocean. The Geoid-Topography Ratio (GTR) analysis of the 85°E and Ninetyeast ridges indicate that they are compensated at shallow depths. Effective elastic thickness ( T e ) estimates obtained through admittance/ coherence analysis as well as the flexural modelling along these ridges led to the conclusions: (i) 85°E Ridge was emplaced in off-ridge environment ( T e  = 10–15 km); (ii) the higher T e values of ~25 km over the Afanasy Nikitin Seamount (ANS) reflect the secondary emplacement of the seamount peaks in off-ridge environment, (iii) that the emplacement of the Ninetyeast Ridge north of 2°N occurred in an off-ridge environment as indicated by higher T e values (25–30 km). Furthermore, the admittance analysis of geoid and bathymetry revealed that the admittance signatures at wavelengths 〉800 km are compensated by processes related to upper mantle convection.

Description: Varied applications of seismic interferometry have arisen in the last decade; however, the potential of the method to improve reflection seismic processing in hardrock environments has not explicitly been regarded. Therefore, in this paper we investigate the potential of retrieving the first arrivals originally hindered by high noise level in the exploitation of controlled-source data acquired over the iron-oxide apatite-rich deposit at Grängesberg (Sweden) and its mining-induced structures. The supervirtual first arrivals retrieved using interferometry methodologies allowed first breaks to be picked more extensively than in the original data. Revised static corrections significantly improved the linearity of the first arrivals and continuity of reflections in the source gathers. Especially, reflections considerably enhanced in the source gathers stacked constructively in the final section. Comparison with geological data, supported by traveltime forward modelling, indicates that these reflections represent the unmined part of the deposit. Other reflections at shallower depth are interpreted as anthropogenic faults possibly located at lithological contacts. Tomographic inversion was also run using the augmented traveltime data. The greater resolution and penetration of this new tomographic image allowed better bridging with the results of the reflection seismic section. Velocity anomalies depict the presence of mining-induced structures and a potential ‘Brewery fault’, which is believed to put at risks an urban area. Even though the potential of first-arrival retrieval seems to be case-dependent, this study illustrates that interferometry may substantially improve the accuracy of static corrections and subsequent stack for hardrock imaging, as well as in the penetration and resolution of traveltime tomograms.

Description: The Arctic continental margin contains large amounts of methane in the form of methane hydrates. The west Svalbard continental slope is an area where active methane seeps have been reported near the landward limit of the hydrate stability zone. The presence of bottom simulating reflectors (BSRs) on seismic reflection data in water depths greater than 600 m suggests the presence of free gas beneath gas hydrates in the area. Resistivity obtained from marine controlled source electromagnetic (CSEM) data provides a useful complement to seismic methods for detecting shallow hydrate and gas as they are more resistive than surrounding water saturated sediments. We acquired two CSEM lines in the west Svalbard continental slope, extending from the edge of the continental shelf (250 m water depth) to water depths of around 800 m. High resistivities (5–12 m) observed above the BSR support the presence of gas hydrate in water depths greater than 600 m. High resistivities (3–4 m) at 390–600 m water depth also suggest possible hydrate occurrence within the gas hydrate stability zone (GHSZ) of the continental slope. In addition, high resistivities (4–8 m) landward of the GHSZ are coincident with high-amplitude reflectors and low velocities reported in seismic data that indicate the likely presence of free gas. Pore space saturation estimates using a connectivity equation suggest 20–50 per cent hydrate within the lower slope sediments and less than 12 per cent within the upper slope sediments. A free gas zone beneath the GHSZ (10–20 per cent gas saturation) is connected to the high free gas saturated (10–45 per cent) area at the edge of the continental shelf, where most of the seeps are observed. This evidence supports the presence of lateral free gas migration beneath the GHSZ towards the continental shelf.

Description: Full waveform inversion (FWI) is a promising technique for inverting a high-resolution subsurface velocity model. The success of FWI highly depends on a fairly well initial velocity model. We propose a method for building a good initial velocity model that can be put into the FWI framework for inverting a nearly perfect velocity structure. We use a well log interpolated velocity model as a high-fidelity initial model for the subsequent FWI. The interpolation problem is solved via a least-squares method with a geological structural regularization. In order to obtain the geological structure of subsurface reflectors, an initial reverse time migration (RTM) with a fairly realistic initial velocity model is conducted, and the local slope of subsurface structure is roughly calculated from the RTM image. The well log interpolated initial velocity model can be very close to the true velocity while containing a small velocity anomaly or oversmoothing caused by the imperfect velocity interpolation. The anomaly and oversmoothing effect can be compensated during the subsequent FWI iterations. We use a relatively simple-layered model and the more complicated Marmousi velocity model to demonstrate the applicability of the proposed approach. We start from a very smooth velocity model and obtain a nearly perfect FWI result which is much better than the traditional FWI result without the velocity interpolation. The migrated images from the RTM method using different velocity models are also compared to further confirm the effectiveness of the proposed framework. Regarding the field deployment, we suggest that future drilling of exploration wells can be seismic-oriented, which can help fully utilize the information of well logs for building initial subsurface velocity model and will facilitate a wide application of the proposed methodology.

Description: Series of experiments were performed to study the spectral induced polarization (SIP) response of sands fully saturated with water, into which gaseous CO 2 or N 2 was injected, in the frequency range 0.5 Hz–20 kHz. Three main observations were made. (1) SIP parameters were always most affected by gas injection when the frequency of the injected signal was in the intermediate range (1 〈 f  〈 20 kHz). This point emphasizes the interest of broadening the frequency range of SIP surveys toward intermediate frequencies. It also implies that more work is needed in order to understand and quantify the parasitic effects that occur at this frequency range (EM coupling and electrode polarization). (2) Through all the experiments performed, we were able to distinguish the parameters variations caused by a reduction of the water saturation level (invasion of a resistive gas phase in the porous space) from those resulting from dissolution processes (increase of the electrical conductivity of the saturating water). (3) The quadrature conductivity '', which is mainly linked to the capacitive properties and inductive response of the media, always shows relative variations stronger than the resistivity of the media, hence demonstrating the interest of the additional information provided by the SIP method as against the classical resistivity method.

Description: The intermediate frequency range 10 2 –10 5  Hz forms the transition range between the spectral induced polarization frequency domain and the dielectric spectroscopy frequency domain. Available experimental data showed that the spectral induced polarization response of sands fully saturated with water was particularly sensitive to variations of the saturating water electrical conductivity value in the intermediate frequency range. An empirical and a mechanistic model have been developed and confronted to this experimental data. This confrontation showed that the Maxwell Wagner polarization alone is not sufficient to explain the observed signal in the intermediate frequency range. The SIP response of the media was modelled by assigning relatively high dielectric permittivity values to the sand particle or high effective permittivity values to the media. Such high values are commonly observed in the dielectric spectroscopy literature when entering the intermediate frequency range. The physical origin of these high dielectric permittivity values is discussed (grain shape, electromagnetic coupling), and a preliminary study is presented which suggests that the high impedance values of the non-polarizable electrodes might play a significant role in the observed behaviour.