ALBERT

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.

Beschreibung: 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.