3D simultaneous joint PP‐PS prestack seismic inversion at Schiehallion field,United Kingdom Continental Shelf

A workflow for simultaneous joint PP‐PS prestack inversion of data from the Schiehallion field on the United Kingdom Continental Shelf is presented and discussed. The main challenge, describing reasonable PS to PP data registration before any prestack or joint PP‐PS inversion, was overcome thanks to a two‐stage process addressing the signal envelope, then working directly on the seismic data to estimate appropriate time‐variant time‐shift volumes. We evaluated the benefits of including PS along with PP prestack seismic data in a joint inversion process to improve the estimated elastic property quality and also to enable estimation of density compared with other prestack and post‐stack inversion approaches. While the estimated acoustic impedance exhibited a similar quality independent of the inversion used (PP post‐stack, PP prestack or joint PP‐PS prestack inversion) the shear impedance estimation was noticeably improved by the joint PP‐PS prestack inversion when compared to the PP prestack inversion. Finally, the density estimated from joint PP and PS prestack data demonstrated an overall good quality, even where not well‐controlled. The main outcome of this study was that despite several data‐related limitations, inverting jointly correctly processed PP and PS data sets brought extra value for reservoir delineation as opposed to PP‐only or post‐stack inversion.     

2D and 3D amplitude-preserving elastic reverse time migration based on the vector-decomposed P- and S-wave records

The elastic reverse time migration approach based on the vector-wavefield decomposition generally uses the scalar product imaging condition to image the multicomponent seismic data. However, the resulting images contain the crosstalk artefacts and the polarity reversal problems, which are caused by the nonphysical wave modes and the angle-dependent reduction of image amplitudes, respectively. To overcome these two problems, we develop an amplitude-preserving elastic reverse time migration approach based on the vector-decomposed P- and S-wave seismic records. This approach includes two key points. The first is that we employ the vector-decomposed P- and S-wave multicomponent records to independently reconstruct the PP and PS reflection images to mitigate the crosstalk artefacts. The second is that we propose two schemes in addressing the issue of polarity reversal problem in the conventional PP image. One solution is to adopt the angle-dependent equation. Another one is to reconstruct an amplitude-preserving PP image with a separated scalar P-wave particle velocity, which has a clear physical meaning. Numerical examples using two-dimensional and three-dimensional models demonstrate that the proposed elastic reverse time migration approach can provide the images with better amplitude-preserving performance and fewer crosstalk artefacts, compared with the conventional elastic reverse time migration approach based on the scalar product imaging condition.     

A discontinuous Galerkin method for seismic wave propagation in coupled elastic and poroelastic media

Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement.     

Rock elasticity as a function of the uniaxial stress: laboratory measurements and theoretical modelling of vertical transversely isotropic and orthorhombic shales

We apply a rock-physics model that describes the relationship between the effective stress and rock elasticity. We experimentally obtain and analyse a data set containing one vertical transversely isotropic and one orthorhombic shale sample. The vertical transversely isotropic symmetry of the first sample is caused by the layered structure of the rock. The seismic orthorhombicity of the second sample could be explained after microscopic analysis of thin section, which demonstrates an imperfect disorder of inhomogeneities. Both samples were loaded uniaxially in a quasi-static regime. During the loading, we measured stress-dependent seismic velocities and sample deformations. For the analysis of the stress-dependent velocities and stiffnesses, we modelled the measured data set using a recent generalization of the porosity deformation approach. Comparison of the experimentally determined and numerically modelled data supports the applicability of the theory and helps in the interpretation of experimentally obtained data. In agreement with the theory, uniaxial stress increases the elliptic component of the seismic anisotropy and does not impact the anellipticity parameter. We demonstrate the distinct influence of the stiff and compliant porosities on the stress sensitivity of the elastic properties.     

Static and dynamic stiffness measurements with Opalinus Clay

In this work, an experimental study was carried out with the aim of reconciling static and dynamic stiffness of Opalinus Clay. The static and dynamic stiffness of core plugs from a shaly and a sandy facies of Opalinus Clay were characterized at two different stress states. The measurements included undrained quasi-static loading–unloading cycles from which the static stiffness was derived, dynamic stiffness measurement at seismic frequencies (0.5–150 Hz) and ultrasonic velocity measurements (500 kHz) probing the dynamic stiffness at ultrasonic frequencies. The experiments were carried out in a special triaxial low-frequency cell. The obtained results demonstrate that the difference between static and dynamic stiffness is due to both dispersion and non-elastic effects: Both sandy and shaly facies of Opalinus Clay exhibit large dispersion, that is, a large frequency dependence of dynamic stiffness and acoustic velocities. Especially dynamic Young's moduli exhibit very high dispersion; between seismic and ultrasonic frequencies they may change by more than a factor 2. P-wave velocities perpendicular to bedding are by more than 200 m/s higher at ultrasonic frequencies than at seismic frequencies. The static undrained stiffness of both sandy and shaly facies is strongly influenced by non-elastic effects, resulting in significant softening during both loading and unloading with increasing stress amplitude. The zero-stress extrapolated static undrained stiffness, however, reflects the purely elastic response and agrees well with the dynamic stiffness at seismic frequency.     

Measurements of the elastic and anelastic properties of sandstone flooded with supercritical CO2

The aim of this study was to investigate the effects of supercritical CO₂ (scCO₂) injection on the elastic and anelastic properties of sandstone at seismic and ultrasonic frequencies. We present the results of the low‐frequency and ultrasonic experiments conducted on water‐saturated sandstone (Donnybrook, Western Australia) flooded with scCO₂. The sandstone was cut in the direction perpendicular to a formation bedding plane and tested in a Hoek triaxial pressure cell. During the experiments with scCO_2, the low‐frequency and ultrasonic systems and the pump dispensing scCO₂ were held at a temperature of 42°C. The elastic parameters obtained for the sandstone with scCO₂ at seismic (0.1 Hz–100 Hz) and ultrasonic (～0.5 MHz) frequencies are very close to those for the dry rock. The extensional attenuation was also measured at seismic frequencies for the dry, water‐saturated, and scCO₂‐injected sandstones. The applicability of Gassmann's fluid substitution theory to obtained results was also tested during the experiments.     

Effect of micro‐inhomogeneity on the effective stress coefficients and undrained bulk modulus of a poroelastic medium: a double spherical shell model

Although most rocks are complex multi‐mineralic aggregates, quantitative interpretation workflows usually ignore this complexity and employ Gassmann equation and effective stress laws that assume a micro‐homogeneous (mono‐mineralic) rock. Even though the Gassmann theory and effective stress concepts have been generalized to micro‐inhomogeneous rocks, they are seldom if at all used in practice because they require a greater number of parameters, which are difficult to measure or infer from data. Furthermore, the magnitude of the effect of micro‐heterogeneity on fluid substitution and on effective stress coefficients is poorly understood. In particular, it is an open question whether deviations of the experimentally measurements of the effective stress coefficients for drained and undrained elastic moduli from theoretical predictions can be explained by the effect of micro‐heterogeneity. In an attempt to bridge this gap, we consider an idealized model of a micro‐inhomogeneous medium: a Hashin assemblage of double spherical shells. Each shell consists of a spherical pore surrounded by two concentric spherical layers of two different isotropic minerals. By analyzing the exact solution of this problem, we show that the results are exactly consistent with the equations of Brown and Korringa (which represent an extension of Gassmann's equation to micro‐inhomogeneous media). We also show that the effective stress coefficients for bulk volume α, for porosity n_? and for drained and undrained moduli are quite sensitive to the degree of heterogeneity (contrast between the moduli of the two mineral components). For instance, while for micro‐homogeneous rocks the theory gives n_? = 1, for strongly micro‐inhomogenous rocks, n_? may span a range of values from –∞ to ∞ (depending on the contrast between moduli of inner and outer shells). Furthermore, the effective stress coefficient for pore volume (Biot–Willis coefficient) α can be smaller than the porosity ?. Further studies are required to understand the applicability of the results to realistic rock geometries.     

Application of the nearly perfectly matched layer for seismic wave propagation in 2D homogeneous isotropic media

Numerical modelling plays an important role in helping us understand the characteristics of seismic wave propagation. The presence of spurious reflections from the boundaries of the truncated computational domain is a prominent problem in finite difference computations. The nearly perfectly matched layer has been proven to be a very effective boundary condition to absorb outgoing waves in both electromagnetic and acoustic media. In this paper, the nearly perfectly matched layer technique is applied to elastic isotropic media to further test the method's absorbing ability. The staggered‐grid finite‐difference method (fourth‐order accuracy in space and second‐order accuracy in time) is used in the numerical simulation of seismic wave propagation in 2D Cartesian coordinates. In the numerical tests, numerical comparisons between the nearly perfectly matched layer and the convolutional perfectly matched layer, which is considered the best absorbing layer boundary condition, is also provided. Three numerical experiments demonstrate that the nearly perfectly matched layer has a similar performance to the convolutional perfectly matched layer and can be a valuable alternative to other absorbing layer boundary conditions.     

Orthotropic anisotropy: on contributions of elasticity parameters to a difference in quasi-P-wave-squared velocities resulted from propagation in two orthogonal symmetry planes

We investigate the dependence of quasi P-wave phase velocity propagating in orthotropic media on particular elasticity parameters. Specifically, due to mathematical facilitation, we consider the squared-velocity difference, , resulted from propagation in two mutually perpendicular symmetry planes. In the context of the effective medium theory, may be viewed as a parameter evaluating the influence of cracks – embedded in the background medium – parallel to one or both aforementioned planes. Our investigation is both theoretical and numerical. Based on Christoffel's equations, we propose two accurate approximations of . Due to them, we interpret the aforementioned squared-velocity difference as being twice more dependent on , than on . To describe the magnitude of the dependence, we consider the proportions between the partial derivatives of . Further, it occurs that is influenced by the ratio of vertically propagating quasi P-wave to vertically propagating quasi S-wave. Anomalously high might be caused by the low P/S ratio, which in turn can be an indicator of the presence of gas in natural fractures or aligned porosity. Also, we carry out numerical sensitivity study, according to which is approximately twice more dependent on than on , twice more sensitive to than to , and equally dependent on as on . The dependence on and can be neglected, especially for small phase angles. We verify the approximations and perform the sensitivity study, using eight examples of the elasticity tensors.