ALBERT

Description: Passive interferometry technology is based on the relation between the reflection and the transmission responses of the subsurface. The transmission response can be received at surface in the presence of the ambient noise source in the subsurface with the cross-correlation (CC) or multidimensional deconvolution methods. We investigate the feasibility of electromagnetic (EM) wave passive interferometry with CC method. We design a 2-D finite-difference time domain (FDTD) algorithm to simulate the long-duration ground penetrating radar (GPR) measurements with random distribution of passive EM sources. The noise sources have random duration time, waveform and spatial distribution. We test the FDTD GPR passive interferometry code with above source characteristics and apply the method to light non-aqueous phase liquid (LNAPL) monitoring. Based on the model simulation data, by using common midpoint velocity analysis and normal move out correction to process the interferometry retrieve record, we can accurately obtain the dynamic changing characteristics of the target's permittivity. The LNAPL dynamic leakage model can be imaged as well. The synthetic results demonstrate that the GPR passive interferometry is feasible in subsurface LNAPL monitoring. Our work provides a foundation for a passive interferometry field application using GPR.

Description: 3-D electrical resistivity surveys and inversion models are required to accurately resolve structures in areas with very complex geology where 2-D models might suffer from artefacts. Many 3-D surveys use a grid where the number of electrodes along one direction ( x ) is much greater than in the perpendicular direction ( y ). Frequently, due to limitations in the number of independent electrodes in the multi-electrode system, the surveys use a roll-along system with a small number of parallel survey lines aligned along the x -direction. The ‘Compare R’ array optimization method previously used for 2-D surveys is adapted for such 3-D surveys. Offset versions of the inline arrays used in 2-D surveys are included in the number of possible arrays (the comprehensive data set) to improve the sensitivity to structures in between the lines. The array geometric factor and its relative error are used to filter out potentially unstable arrays in the construction of the comprehensive data set. Comparisons of the conventional (consisting of dipole-dipole and Wenner–Schlumberger arrays) and optimized arrays are made using a synthetic model and experimental measurements in a tank. The tests show that structures located between the lines are better resolved with the optimized arrays. The optimized arrays also have significantly better depth resolution compared to the conventional arrays.

Description: We apply a reversible-jump Markov chain Monte Carlo method to sample the Bayesian posterior model probability density function of 2-D seafloor resistivity as constrained by marine controlled source electromagnetic data. This density function of earth models conveys information on which parts of the model space are illuminated by the data. Whereas conventional gradient-based inversion approaches require subjective regularization choices to stabilize this highly non-linear and non-unique inverse problem and provide only a single solution with no model uncertainty information, the method we use entirely avoids model regularization. The result of our approach is 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. We represent models in 2-D using a Voronoi cell parametrization. To make the 2-D problem practical, we use a source–receiver common midpoint approximation with 1-D forward modelling. Our algorithm is transdimensional and self-parametrizing where the number of resistivity cells within a 2-D depth section is variable, as are their positions and geometries. Two synthetic studies demonstrate the algorithm's use in the appraisal of a thin, segmented, resistive reservoir which makes for a challenging exploration target. As a demonstration example, we apply our method to survey data collected over the Scarborough gas field on the Northwest Australian shelf.

Description: Compressional and shear wave seismic measurements were performed in an old railway tunnel and in galleries excavated in a 250-m-thick Toarcian claystone formation in the Tournemire experimental station (France). Three component (3C) geophones and three orthogonal orientations of the vibroseismic force source were used. Additionally, vertical seismic profiling (VSP) measurements were recorded with a 3C borehole geophone, a hydrophone and a microphone in a 159 m deep borehole (ID180) in the tunnel. The seismic data show that Toarcian claystone has strong transverse isotropy (TI) with a vertical symmetry axis. The qP , SH and qSV wave propagation velocities in horizontal directions—the plane of isotropy of the TI medium—are measured as 3550, 1850 and 1290 m s –1 , respectively. The zero-offset VSP reveals that only one shear wave propagates in the vertical (depth) direction and the P - and S -wave velocities are 3100 and 1375 m s –1 , respectively. Four elastic moduli of the TI medium are determined from the seismic velocities and from the bulk density of 2.53 g cm –3 : c 11 = 31.9 GPa, c 33 = 24.3 GPa, c 44 = 4.5 GPa and c 66 = 8.7 GPa. A walkaway VSP with the borehole geophone at 50 m depth in borehole ID180 and shot points in the galleries leads to oblique seismic ray paths which allow us to determine the fifth elastic modulus of the TI medium to c 13 = 16 GPa. The tube wave recorded by a hydrophone in the water filled lower part of the borehole propagates with 1350 m s –1 , which confirms the estimate of the elastic constant c 66 . The analysis of body wave and surface wave data from a seismic experiment in Galerie Est shows reflections from several fracture zones in the gallery floor. The thickness of the excavation damaged zone (EDZ) in the floor of Galerie Est is estimated to 0.7 m.

Description: The Tonga-Kermadec forearc is deforming in response to on-going subduction of the Pacific Plate beneath the Indo-Australian Plate. Previous research has focussed on the structural development of the forearc where large bathymetric features such as the Hikurangi Plateau and Louisville Ridge seamount chain are being subducted. Consequently, knowledge of the ‘background’ forearc in regions of normal plate convergence is limited. We report on an ~250-km-long multichannel seismic reflection profile that was shot perpendicular to the Tonga-Kermadec trench at ~28°S to determine the lateral and temporal variations in the structure, stratigraphy and deformation of the Kermadec forearc resulting solely from Pacific Plate subduction. Interpretation of the seismic profile, in conjunction with regional swath bathymetry data, shows that the Pacific Plate exhibits horst and graben structures that accommodate bending-induced extensional stresses, generated as the trenchward dip of the crust increases. Trench infill is also much thicker than expected at 1 km which, we propose, results from increased sediment flux into and along the trench. Pervasive normal faulting of the mid-trench slope most likely accommodates the majority of the observed forearc extension in response to basal subduction erosion, and a structural high is located between the mid- and upper-trench slopes. We interpret this high as representing a dense and most likely structurally robust region of crust lying beneath this region. Sediment of the upper-trench slope documents depositional hiatuses and on-going uplift of the arc. Strong along-arc currents appear to erode the Kermadec volcanic arc and distribute this sediment to the surrounding basins, while currents over the forearc redistribute deposits as sediment waves. Minor uplift of the transitional Kermadec forearc, observed just to the north of the profile, appears to relate to an underlying structural trend as well as subduction of the Louisville Ridge seamount chain 250 km to the north. Relative uplift of the Kermadec arc is observed from changes in the tilt of upper-trench slope deposits and extensional faulting of the basement immediately surrounding the Louisville Ridge.

Description: New methods are required to combine the information contained in the passive electrical and seismic signals to detect, localize and monitor hydromechanical disturbances in porous media. We propose a field experiment showing how passive seismic and electrical data can be combined together to detect a preferential flow path associated with internal erosion in a Earth dam. Continuous passive seismic and electrical (self-potential) monitoring data were recorded during a 7-d full-scale levee (earthen embankment) failure test, conducted in Booneschans, Netherlands in 2012. Spatially coherent acoustic emissions events and the development of a self-potential anomaly, associated with induced concentrated seepage and internal erosion phenomena, were identified and imaged near the downstream toe of the embankment, in an area that subsequently developed a series of concentrated water flows and sand boils, and where liquefaction of the embankment toe eventually developed. We present a new 4-D grid-search algorithm for acoustic emissions localization in both time and space, and the application of the localization results to add spatially varying constraints to time-lapse 3-D modelling of self-potential data in the terms of source current localization. Seismic signal localization results are utilized to build a set of time-invariant yet spatially varying model weights used for the inversion of the self-potential data. Results from the combination of these two passive techniques show results that are more consistent in terms of focused ground water flow with respect to visual observation on the embankment. This approach to geophysical monitoring of earthen embankments provides an improved approach for early detection and imaging of the development of embankment defects associated with concentrated seepage and internal erosion phenomena. The same approach can be used to detect various types of hydromechanical disturbances at larger scales.

Description: We propose a two-phase damage theory in a viscoelastic medium to study the pressure and porosity diffusion in fractured near-surface porous rocks. The key ingredient in the viscoelastic theory is that the pressure difference between solid and fluid is divided into three parts, which contribute to reversible elastic potential energy, irreversible viscous entropy production and surface energy stored during deformation. The resulting continuum description of weakening and failure (distributed void generation and microcracking) in a linear Kelvin body accounts for surface energy being created by both viscous and elastic deformational work. The model shows that while non-linear permeability models leads to an enhanced diffusivity, damage makes the matrix more compressible if we assume the geometry/size of cracks remain unchanged. The net effect is that the porosity diffusivity is reduced causing fluid infiltration to accumulate closer to the injection source, leading to a slower fluid diffusion during hydraulic fracturing with a fixed porosity boundary condition. However if a constant overpressure boundary condition is applied, a weakened matrix with damage leads to greater pressure diffusivity than for porosity.

Description: We define an algorithm in the time domain for computing the unwrapped instantaneous phase and its derivative, the instantaneous frequency, using only derivatives and integrals. It does not require user-defined parameters, like most algorithms proposed so far. We validate and compare its performance with respect to open-source and commercial software by synthetic and real data examples.

Description: Simulating electromagnetic fields in the quasi-static regime by solving Maxwell's equations is a central task in many geophysical applications. In most cases, geophysical targets of interest exhibit complex topography and bathymetry as well as layers and faults. Capturing these effects with a sufficient level of detail is a huge challenge for numerical simulations. Standard techniques require a very fine discretization that can result in an impracticably large linear system to be solved. A remedy is to use locally refined and adaptive meshes, however, the potential coarsening is limited in the presence of highly heterogeneous and anisotropic conductivities. In this paper, we discuss the application of multiscale finite volume (MSFV) methods to Maxwell's equations in frequency domain. Given a partition of the fine mesh into a coarse mesh the idea is to obtain coarse-to-fine interpolation by solving local versions of Maxwell's equations on each coarsened grid cell. By construction, the interpolation accounts for fine scale conductivity changes, yields a natural homogenization, and reduces the fine mesh problem dramatically in size. To improve the accuracy for singular sources, we use an irregular coarsening strategy. We show that using MSFV methods we can simulate electromagnetic fields with reasonable accuracy in a fraction of the time as compared to state-of-the-art solvers for the fine mesh problem, especially when considering parallel platforms.

Description: Wave-equation tomography (WT) and full waveform inversion (FWI) are combined through a hybrid misfit function to estimate high-resolution subsurface structures starting from a poorly constrained initial velocity model. Both methods share the same wavefield forward modelling and inversion schemes, while they differ only on the ways to calculate misfit functions and hence the ways to sample in the model space. Aiming at minimizing the cross-correlation phase delay between synthetic and real data, WT can be used to retrieve the long- and middle-wavelength model components, which are essential to FWI. Compared to ray-based traveltime tomography that is based on asymptotic high-frequency approximation, WT provides a better resolution by exploring the band-limited feature of seismic wavefield. On the other hand, FWI is capable of resolving the short-wavelength model component, complementing the WT. In this study, we apply WT to surface first-arrival refraction data, and apply FWI to both refraction and reflection data. We assign adaptive weights to the two different misfit measurements and build a progressive inversion strategy. To illustrate the advantage of our strategy over conventional ‘ray tomography + FWI’ approach, we show in a synthetic lens test that WT can provide extra subsurface information that is critical for a successful FWI application. To further show the efficiency, we test our strategy on the 2-D Marmousi model where satisfactory inversion results are achieved without much manual intervention. Finally, we apply the inversion strategy to a deep-water seismic data set acquired offshore Sumatra with a 12-km-long streamer. In order to alleviate several practical problems posed by the deep-water setting, we apply downward continuation (DC) to generate a virtual ocean bottom experiment data set prior to inversion. The new geometry after DC boosts up the shallow refractions, as well as avoiding cumbersome modelling through the thick water column, thus reducing the computation cost by 85 per cent. The inversion result from the new data set shows high-resolution shallow sediment structures and the migration images prove the superiority of the inverted model over a conventional tomography model.

Description: The activation of Late Quaternary faults in the Central Apennines (Italy) could generate earthquakes with magnitude of about 6.5, and the Monte Vettore fault system probably belongs to the same category of seismogenetic faults. Such structure has been defined ‘silent’, because of its geological and geomorphological evidences of past activation, but the absence of historical records in the seismic catalogues to be associated with its activation. The ‘Piano di Castelluccio’ intramountain basin, resulting from the Quaternary activity of normal faults, is characterized by a secondary fault strand highlighted by a NW–SE fault scarp: it has been already studied through palaeoseismological trenches, which highlighted evidences of Quaternary shallow faulting due to strong earthquakes, and through a 2-D ground penetrating radar (GPR) survey, showing the first geophysical signature of faulting for this site. Within the same place, a 3-D GPR volume over a 20 20 m area has been collected. The collection of radar echoes in three dimensions allows to map both the vertical and lateral continuity of shallow geometries of the fault zone (Fz), imaging features with high resolution, ranging from few metres to centimetres and therefore imaging also local variations at the microscale. Several geophysical markers of faulting, already highlighted on this site, have been taken as reference to plan the 3-D survey. In this paper, we provide the first 3-D subsurface imaging of an active shallow fault belonging to the Umbria-Marche Apennine highlighting the subsurface fault geometry and the stratigraphic sequence up to a depth of about 5 m. From our data, geophysical faulting signatures are clearly visible in three dimensions: diffraction hyperbolas, truncations of layers, local attenuated zones and varying dip of the layers have been detected within the Fz. The interpretation of the 3-D data set provided qualitative and quantitative geological information in addition to the fault location, like its geometry, boundaries and an estimation of the fault throw.

Description: To estimate the seismic hazard, the geometry (dip, length and orientation) and the dynamics (type of displacements and amplitude) of the faults in the area of interest need to be understood. In this paper, in addition to geomorphologic observations, we present the results of two ground penetrating radar (GPR) campaigns conducted in 2010 and 2011 along the Emeelt fault in the vicinity of Ulaanbaatar, capital of Mongolia, located in an intracontinental region with low deformation rate that induces long recurrence time between large earthquakes. As the geomorphology induced by the fault activity has been highly smoothed by erosion processes since the last event, the fault location and geometry is difficult to determine precisely. However, by using GPR first, a non-destructive and fast investigation, the fault and the sedimentary deposits near the surface can be characterized and the results can be used for the choice of trench location. GPR was performed with a 50 MHz antenna over 2-D lines and with a 500 MHz antenna for pseudo-3-D surveys. The 500 MHz GPR profiles show a good consistency with the trench observations, dug next to the pseudo-3-D surveys. The 3-D 500 MHz GPR imaging of a palaeochannel crossed by the fault allowed us to estimate its lateral displacement to be about 2 m. This is consistent with a right lateral strike-slip displacement induced by an earthquake around magnitude 7 or several around magnitude 6. The 2-D 50 MHz profiles, recorded perpendicular to the fault, show a strong reflection dipping to the NE, which corresponds to the fault plane. Those profiles provided complementary information on the fault such as its location at shallow depth, its dip angle (from 23° to 35°) and define its lateral extension.

Description: A new method to obtain the statistics of a geostatistical model is introduced. The method elicits the statistical information from a geological expert directly, by iteratively updating a population of vectors of statistics, based on the expert's subjective opinion of the corresponding geological simulations. Thus, it does not require the expert to have knowledge of the mathematical and statistical details of the model. The process uses a genetic algorithm to generate new vectors. We demonstrate the methodology for a particular geostatistical model used to model rock pore-space, which simulates the spatial distribution of matrix and pores over a 2-D grid, using multipoint statistics specified by conditional probabilities. Experts were asked to use the algorithm to estimate the statistics of a given target pore-space image with known statistics; thus, their numerical rates of convergence could be calculated. Convergence was measured for all experts, showing that the algorithm can be used to find appropriate probabilities given the expert's subjective input. However, considerable and apparently irreducible residual misfit was found between the true statistics and the estimates of statistics obtained by the experts, with the root-mean-square error on the conditional probabilities typically 〉0.1. This is interpreted as the limit of the experts’ abilities to distinguish between realizations of different spatial statistics using the algorithm. More accurate discrimination is therefore likely to require complementary elicitation techniques or sources of information independent of expert opinion.

Description: We show analytically that a well-known transfer function previously derived for the scalar acoustic problem that converts measurements from a 3-D (real-world) setting to a 2-D equivalent is directly applicable to the vector electromagnetic borehole ground penetrating radar problem. We also show that the transfer function's precision is improved for the low-loss case through the use of complex velocity. The transfer function has a strong effect on amplitude, and is therefore a critical preprocessing step for 2-D full-wave inversion when finding conductivity is of concern. We demonstrate the effectiveness of the transfer function through various numerical experiments and a synthetic frequency-domain full-wave inversion. We also compare the effectiveness of this curved-ray transfer function to a quasi-straight-ray transfer function. The inversion demonstrates the positive effect the transfer functions have on recovering conductivity and also that they are effective even when there are sharp velocity contrasts.

Description: We report and analyse the tsunami recorded in the northwestern Indian Ocean at the Makran region following the M w 7.7 Pakistan inland strike-slip earthquake on 2013 September 24. We analyse eleven tide gauge records as well as one DART record of this tsunami and perform numerical modelling of the tsunami that would be triggered by a range of possible sources. The tsunami registered a maximum wave height of 109 cm at the Qurayat tide gauge station (Oman). The dominant period of the tsunami was around 12 min, although wavelet analysis showed that parts of the tsunami energy were partitioned into a slightly wider period range of 7 and 16 min. Tsunami backward ray tracing showed that the tsunami source was possibly located offshore Jiwani (Pakistan) and that the tsunami was most likely triggered by the main shock. The aftershocks are distributed in the inland region and the coseismic vertical and horizontal displacements are also limited inland implying that the tsunami was generated by secondary sources triggered by the earthquake. Different possible tsunami sources including a mud volcano at the location of the newly generated island, and a mud volcano or diapir at offshore deep water were examined through numerical modelling and all failed to reproduce the observed waveforms. Numerical modelling showed that a submarine slump with a source dimension of about 10–15 km and a thickness of about 100 m located at 61.49°E and 24.62°N, that is, about 60–70 km off the Jiwani coast (Pakistan), seems capable of reasonably reproducing the wave amplitudes and periods of the observed tsunami waveforms. This event was the second instrumentally recorded tsunami in the region, after the Makran tsunami of 1945 November, and provides evidence for a hazard from landslide/slump-generated waves following seismic activity in the area.

Description: Joint inversion of different geophysical data sets is becoming a more popular and powerful tool, and it has been performed on data sensitive both to the same physical parameter and to different physical parameters. Joint inversion is undertaken to reduce acceptable model space and to increase sensitivity to model parameters that one method alone is unable to resolve adequately. We examine and implement a novel hybrid joint inversion approach. In our inversion scheme a model—the reference model—is fixed, and the information shared with the subsurface structure obtained from another method will be maximized; in our case conductivity structures from magnetotelluric (MT) inversion. During inversion, the joint probability distribution of the MT and the specified reference model is estimated and its entropy minimized in order to guide the inversion result towards a solution that is statistically compatible with the reference model. The powerful feature of this technique is that no explicit relationships between estimated model parameters and reference model ones are presumed: if a link exists in data then it is highlighted in the estimation of the joint probability distribution, if no link is required, then none is enforced. Tests performed verify the robustness of this method and the advantages of it in a 1-D anisotropic scenario are demonstrated. A case study was performed with data from Central Germany, effectively fitting an MT data set from a single station within as minimal an amount of anisotropy as required.

Description: Several slope failures are observed near the deformation front on the frontal ridges of the northern Cascadia accretionary margin off Vancouver Island. The cause for these events is not clear, although several lines of evidence indicate a possible connection between the occurrence of gas hydrate and submarine landslide features. The presence of gas hydrate is indicated by a prominent bottom-simulating reflector (BSR), at a depth of ~265–275 m beneath the seafloor (mbsf), as interpreted from vertical-incidence and wide-angle seismic data beneath the ridge crests of the frontal ridges. For one slide, informally called Slipstream Slide, the velocity structure inferred from tomography analyses shows anomalous high velocities values of about 2.0 km s –1 at shallow depths of 100 mbsf. The estimated depth of the glide plane (100 ± 10 m) closely matches the depth of these shallow high velocities. In contrast, at a frontal ridge slide just to the northwest (informally called Orca Slide), the glide plane occurs at the same depth as the current BSR. Our new results indicate that the glide plane of the Slipstream slope failure is associated with the contrast between sediments strengthened by gas hydrate and overlying sediments where little or no hydrate is present. In contrast, the glide plane of Orca Slide is between sediment strengthened by hydrate underlain by sediments beneath the gas hydrate stability zone, possibly containing free gas. Additionally, a set of margin perpendicular normal faults are imaged from seafloor down to BSR depth at both frontal ridges. As inferred from the multibeam bathymetry, the estimated volume of the material lost during the slope failure at Slipstream Slide is about 0.33 km 3 , and ~0.24 km 3 of this volume is present as debris material on the ocean basin floor. The 20 per cent difference is likely due to more widely distributed fine sediments not easily detectable as bathymetric anomalies. These volume estimates on the Cascadia margin are approaching the mass failure volume for other slides that have generated large tsunamis—for example 1–3 km 3 for a 1998 Papua New Guinea slide.

Description: The method presented here aims to assess the tsunami threat very rapidly after the occurrence of a large earthquake, using as input the parameters of the seismic source, and an approach based on Green's summation. We show that the main weakness of the approach (the need to consider only linear shallow water propagation) is largely compensated by the advantages in terms of computing performance and independence with respect to pre-computed scenarios. To test the approach and to illustrate its implementation in a real environment, we focus on the Sea of Oman, a tsunamigenic area characterized by Makran subduction zone which detailed structure is partially unknown and where secondary tsunami sources must also be taken into account, both for hazard studies and warning purposes. The potential source area is partitioned into a grid of unity water sources. A shallow water (SW) numerical model is used to pre-compute the corresponding empirical Green's functions on several points of interest located on the coasts of Iran, Pakistan and Oman. The comparison between Green's summation and the direct SW computation using the full resolution of the bathymetric grid shows that the accuracy is good enough for practical applications.

Description: Clay minerals as products of hydrothermal alteration significantly influence the hydraulic and mechanical properties of crystalline rock. Therefore, the localization and characterization of alteration zones by downhole measurements is a great challenge for the development of geothermal reservoirs. The magnetite bearing granite of the geothermal site in Soultz-sous-Forêts (France) experienced hydrothermal alteration during several tectonic events and clay mineral formation is especially observed in alteration halos around fracture zones. During the formation of clay minerals, magnetite was oxidized into hematite, which significantly reduces the magnetic susceptibility of the granite from ferrimagnetic to mostly paramagnetic values. The aim of this study was to find out if there exists a correlation between synthetic clay content logs (SCCLs) and measurements of magnetic susceptibility on cuttings in the granite in order to characterize their alteration mineralogy. Such a correlation has been proven for core samples of the EPS1 reference well. SCCLs were created from gamma ray and fracture density logs using a neural network. These logs can localize altered fracture zones in the GPK1-4 wells, where no core material is available. Mass susceptibility from 261 cutting samples of the wells GPK1–GPK4 was compared with the neural network derived synthetic logs. We applied a combination of temperature dependent magnetic susceptibility measurements with optical and electron microscopy, and energy dispersive X-ray spectroscopy to discriminate different stages of alteration. We found, that also in the granite cuttings an increasing alteration grade is characterized by an advancing oxidation of magnetite into hematite and a reduction of magnetic susceptibility. A challenge to face for the interpretation of magnetic susceptibility data from cuttings material is that extreme alteration grades can also display increased susceptibilities due to the formation of secondary magnetite. Low magnetic susceptibility can also be attributed to primary low magnetite content, if the granite facies changes. In order to interpret magnetic susceptibility from cuttings, contaminations with iron from wear debris of the drilling tools must be eliminated. Provided that the magnetic mineralogy of the granite is known in detail, this method in combination with petrographic investigations is suited to indicate and characterize hydrothermal alteration and the appearance of clay.

Description: Elastic and mechanical weakening from water saturation are widely known to occur in sedimentary rocks, and particularly in carbonate rocks. To improve our understanding of the physics underlying this phenomenon, ultrasonic ( f  ~ 0.5 MHz) elastic properties are measured on a large suite of clean limestones and sandstones at very low saturations from relative humidity (RH) variations at ambient conditions. Measurements clearly highlight an elastic weakening (i.e. decrease in elastic wave velocity) from moisture adsorption. P - and S -wave velocities are similarly affected by adsorption, but in a different way for limestones and sandstone samples. While the elastic properties of limestone samples show almost no RH dependence, a large weakening is observed for samples of Fontainebleau sandstone that increases with the samples’ porosity. The main elastic weakening effect is likely to result from adsorption of fluid at grain contacts. It thus affects particularly granular rocks such as sandstones while well-cemented limestones are not affected. The granular model from Murphy et al. , accounting for surface energy effects, proves to be appropriate. Applying this model, it is shown that (i) P - and S -wave velocities have the same dependence on surface energy, which is consistent with the measurements and (ii) surface energy values obtained from the ultrasonic data using this model correlate with RH, and are consistent with the expected value for quartz crystals at vapour pressure. Yet, porosity, which relates to degree of cementation in the particular case of Fontainebleau sandstone, appears to be an additional parameter. A modified model is thus derived using the cementation model from Digby, accounting for a bonding radius at grain contact. It proves to apply well to the measured data. The fundamental difference between limestones’ and sandstones’ dependence to RH appears to be related to a microstructural difference. Saturation variations from RH increase depend on specific surface area, which is particularly low in Fontainebleau sandstones and large in microporous limestones. However elastic weakening from RH is more important in sandstones owing to their granular microstructure.

Description: A novel experimental method is introduced to estimate the Thomsen's elastic anisotropy parameters and of a transversely isotropic shale under variable stress and saturation conditions. The method consists in recording P -wave velocities along numerous paths on a cylindrical specimen using miniature ultrasonic transducers. Such an overdetermined set of measurements is specifically designed to reduce the uncertainty associated with the determination of Thomsen's parameter compared to the classical method for which a single off-axis measurement is used (usually at 45° to the specimen's axis). This method is applied to a specimen of Opalinus Clay recovered from the Mont-Terri Underground Research Laboratory in Switzerland. The specimen is first saturated with brine at low effective pressure and then subjected to an effective pressure cycle up to 40 MPa, followed by a triaxial loading up to failure. During saturation and deformation, the evolution of P -wave velocities along a maximum of 240 ray paths is monitored and Thomsen's parameters α , and are computed by fitting Thomsen's weak anisotropy model to the data. The values of and obtained at the highest confining pressures reached during the experiment are comparable with those predicted from X-ray diffraction texture analysis and modelling for Opalinus Clay reported in the literature. These models neglect the effect of soft-porosity on elastic properties, but become relevant when soft porosity is closed at high effective pressure.

Description: We present a novel accurate and efficient goal-oriented adaptive finite-element method solution for complex multi-electrodes resistivity system with arbitrary smooth surface topographies. A simple Green's function of a half-space model is adopted to eliminate the singularity. A unified boundary value problem for the regular potential is formulated for a multi-electrodes system so that it shares a common system matrix. In addition, a goal-oriented error estimation technique is developed to generate an optimal common grid so that highly accurate solutions are obtained with minimum computation cost. Synthetic models are used to verify our algorithm and excellent agreements are obtained by comparing with other methods.

Description: The Geological Survey of Sweden has been collecting airborne tensor very low frequency data (VLF) over several decades, covering large parts of the country. The data has been an invaluable source of information for identifying conductive structures that can among other things be related to water-filled fault zones, wet sediments that fill valleys or ore mineralizations. Because the method only uses two differently polarized plane waves of very similar frequency, vertical resolution is low and interpretation is in most cases limited to maps that are directly derived from the data. Occasionally, 2-D inversion is carried out along selected profiles. In this paper, we present for the first time a 3-D inversion for tensor VLF data in order to further increase the usefulness of the data set. The inversion is performed using a non-linear conjugate gradient scheme (Polak-Ribière) with an inexact line-search. The gradient is obtained by an algebraic adjoint method that requires one additional forward calculation involving the adjoint system matrix. The forward modelling is based on integral equations with an analytic formulation of the half-space Green's tensor. It avoids typically required Hankel transforms and is particularly amenable to singularity removal prior to the numerical integration over the volume elements. The system is solved iteratively, thus avoiding construction and storage of the dense system matrix. By using fast 3-D Fourier transforms on nested grids, subsequently farther away interactions are represented with less detail and therefore with less computational effort, enabling us to bridge the gap between the relatively short wavelengths of the fields (tens of metres) and the large model dimensions (several square kilometres). We find that the approximation of the fields can be off by several per cent, yet the transfer functions in the air are practically unaffected. We verify our code using synthetic calculations from well-established 2-D methods, and trade modelling accuracy off against computational effort in order to keep the inversion feasible in both respects. Our compromise is to limit the permissible resistivity to not fall below 100 m to maintain computational domains as large as 10 10 km 2 and computation times on the order of a few hours on standard PCs. We investigate the effect of possible local violations of these limits. Even though the conductivity magnitude can then not be recovered correctly, we do not observe any structural artefacts related to this in our tests. We invert a data set from northern Sweden, where we find an excellent agreement of known geological features, such as contacts or fault zones, with elongated conductive structures, while high resistivity is encountered in probably less disturbed geology, often related to topographic highs, which have survived predominantly glacial erosion processes. As expected from synthetic studies, the resolution is laterally high, but vertically limited down to the top of conductive structures.

Description: Induced polarization is a geophysical method looking to image and interpret low-frequency polarization mechanisms occurring in porous media. Below 10 kHz, the quadrature conductivity of metal-free sandy and clayey materials exhibits a distribution of relaxation times, which can be related to the pore size distribution of these porous materials. When the polarization spectra are fitted with a Cole–Cole model, we first observe that the main relaxation time is controlled by the main pore size of the material and that the Cole–Cole exponent c is never much above 0.5, a value corresponding to a Warburg function. The complex conductivity is then obtained through a convolution product between the pore size distribution and such Warburg function. We also provide a way to recover the pore size distribution by performing a deconvolution of measured spectra using the Warburg function. A new dataset of mercury porosimetry and induced polarization data of six siliciclastic materials supports the hypothesis that the Cole–Cole relaxation time is strongly controlled by the pore size, and especially the characteristic pore size corresponding to the peak of the pore size distribution from mercury porosimetry. The distribution of the pore throat sizes of these materials seems fairly well recovered using the Warburg decomposition of the spectral induced polarization spectra but additional data will be needed to confirm this finding.

Description: We propose a new, simple and efficient method to image electrical resistivity between a set of wells. Our procedure consists of two steps: first, we map the interfaces between various subsurface formations using seismoelectric conversions; second, we derive the formation resistivity using image-guided cross-well electric tomography. In the first step, we focus seismic energy at a set of points located on a regular grid between wells, which enables us to map the geological formations in terms of heterogeneities in electrical, hydraulic and/or seismic properties. The density of the scanning points (i.e. the seismoelectric image resolution) is related to the wavelength of the seismic impulse used to scan the formations. Each time the seismic energy is focused at a point, the resulting electrical potential burst (equivalent to the one generated by a volumetric seismic source) is recorded remotely at a set of electrodes positioned in wells (the reference electrode can be located on the ground surface or far enough to be considered at infinity). We construct a high-resolution ‘seismoelectric’ image by assigning the electrical potential simulated at these fixed electrodes to the location of the seismic focus. In a follow-up step, the structure of this image is used in image-guided inversion to improve electrical resistivity tomography between the two wells. The structural information from the seismoelectric image is used to impose constraints on the model covariance matrix used in the inversion of the electrical resistivity data. This approach offers new perspectives in recovering fine structure of resistivity (high definition resistivity tomography) between the wells, which cannot be resolved through conventional cross-well resistivity or from seismic tomography alone.

Description: In subduction zones, shallow subsurface structures are the manifestation of the plate interactions at depth. However, significant water depths, rough bathymetry and presence of heavily deformed accretionary wedge materials hamper imaging of the near-surface features to a great extent using conventional imaging techniques. In this study, we show results using an integrated processing technique to a multichannel seismic data set acquired in 2006 from the northwestern offshore Sumatra. We start with first downward continuing the 12-km-long surface streamer data to the seafloor, followed by a high-resolution traveltime tomography of refracted phases to determine a detailed velocity–depth model of subsurface, which in turns, is used for pre-stack depth migration in order to delineate the shallow subsurface structures beneath the trench, subduction front and outer accretionary wedge. Our velocity–depth model and the depth migrated image depict variation of sediment properties across the front and structures of uppermost sedimentary sequence with an unprecedented high resolution providing the precise location of the frontal and conjugate thrusts, highly folded sedimentary sequences, which in turns describe their relationship with the top of the subducting plate and factors that control rupture propagation to the trench. Furthermore, we estimate the porosity distribution across the front, where we find a 12 and 18 per cent decrease in porosity beneath the deformation front and the inner accretionary plateau at 500 m below the seafloor, respectively, which we interpret to be due to the compaction. A significant decrease in porosity at the plate interface below 5–6-km thick sediments near the deformation front would increase the coupling, leading to the rupture propagation up to the trench, uplifting 4.5 km water and producing large tsunami.

Description: We present a scale- and parameter-adaptive method to pre-condition the gradient of the parameters to be inverted in time-domain 2-D elastic full-waveform inversion (FWI). The proposed technique, which relies on a change of variables of the model parameters, allows to balance the value of the gradient of the Lamé parameters and density throughout the model in each step of the multiscale inversion. The main difference compared to existing gradient pre-conditioners is that the variables are automatically selected based on a least-squares minimization criteria of the gradient weight, which corresponds to the product of the gradient by a power of the parameter to be inverted. Based on numerical tests made with (1) a modified version of the Marmousi-2 model, and (2) a high-velocity and density local anomaly model, we illustrate that the value of the power helps to balance the gradient throughout the model. In addition, we show that a particular value exists for each parameter that optimizes the inversion results in terms of accuracy and efficiency. For the two models, the optimal power is ~2.0–2.5 and ~1.5 for the first and second Lamé parameters, respectively; and between 3 and 6, depending on the inverted frequency, for density. These power values provide the fastest and most accurate inversion results for the three parameters in the framework of multiscale and multishooting FWI using three different optimization schemes.