ALBERT

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