ALBERT

Description: ABSTRACT Up–down wavefield decomposition is effectuated by a scaled addition or subtraction of the pressure and vertical particle velocity, generally on horizontal or vertical surfaces, and works well for data given on such surfaces. The method, however, is not applicable to decomposing a wavefield when it is given at one instance in time, i.e. on snapshots. Such situations occur when a wavefield is modelled with methods like finite‐difference techniques, for the purpose of, for example, reverse time migration, where the entire wavefield is determined per time instance. We present an alternative decomposition method that is exact when working on snapshots of an acoustic wavefield in a homogeneous medium, but can easily be approximated to heterogeneous media, and allows the wavefield to be decomposed in arbitrary directions. Such a directional snapshot wavefield decomposition is achieved by recasting the acoustic system in terms of the time derivative of the pressure and the vertical particle velocity, as opposed to the vertical derivative in up–down decomposition for data given on a horizontal surface. As in up–down decomposition of data given at a horizontal surface, the system can be eigenvalue decomposed and the inverse of the eigenvector matrix decomposes the wavefield snapshot into fields of opposite directions, including up–down decomposition. As the vertical particle velocity can be rotated at will, this allows for decomposition of the wavefield into any spatial direction; even spatially varying directions are possible. We show the power and effectiveness of the method by synthetic examples and models of increasing complexity.

Description: ABSTRACT This article addresses the question whether time-lapse seismic reflection techniques can be used to follow and quantify the effects of solution salt mining. Specifically, the production of magnesium salts as mined in the north of the Netherlands is considered. The use of seismic time-lapse techniques to follow such a production has not previously been investigated. For hydrocarbon production and CO 2 storage, time-lapse seismics are used to look at reservoir changes mainly caused by pressure and saturation changes in large reservoirs, while for solution mining salt is produced from caverns with a limited lateral extent, with much smaller production volumes and a fluid (brine) replacing a solid (magnesium salt). In our approach we start from the present situation of the mine and then study three different production scenarios, representing salt production both in vertical and lateral directions of the mine. The present situation and future scenarios have been transformed into subsurface models that were input to an elastic finite-difference scheme to create synthetic seismic data. These data have been analysed and processed up to migrated seismic images, such that time-lapse analyses of intermediate and final results could be done. From the analyses, it is found that both vertical and lateral production is visible well above the detection threshold in difference data, both at pre-imaging and post-imaging stages. In quantitative terms, an additional production of the mine of 6 m causes time-shifts in the order of 2 ms (pre-imaging) and 4 ms (post-imaging) and amplitude changes of above 20% in the imaged sections. A laterally oriented production causes even larger amplitude changes at the edge of the cavern due to replacement of solid magnesium salt with brine introducing a large seismic contrast. Overall, our pre-imaging and post-imaging time-lapse analysis indicates that the effects of solution salt mining can be observed and quantified on seismic data. The effects seem large enough to be observable in real seismic data containing noise.

Description: ABSTRACT When a seismic source is placed in the water at a height less than a wavelength from the water–solid interface, a prominent S-wave arrival can be observed. It travels kinematically as if it was excited at the projection point of the source on the interface. This non-geometric S-wave has been investigated before, mainly for a free-surface configuration. However, as was shown in a field experiment, the non-geometric S-wave can also be excited at a fluid–solid configuration if the S-wave speed in the solid is less than the sound speed in the water. The amplitude of this wave exponentially decreases when the source is moved away from the interface revealing its evanescent character in the fluid. In the solid, this particular converted mode is propagating as an ordinary S-wave and can be transmitted and reflected as such. There is a specific region of horizontal slownesses where this non-geometric wave exists, depending on the ratio of the S-wave velocity and the sound speed of water. Only for ratios smaller than 1, this wave appears. Lower ratios result in a wider region of appearance. Due to this property, this particular P-S converted mode can be identified and filtered from other events in the Radon domain.

Description: ABSTRACT Wavefield decomposition forms an important ingredient of various geophysical methods. An example of wavefield decomposition is the decomposition into upgoing and downgoing wavefields and simultaneous decomposition into different wave/field types. The multi-component field decomposition scheme makes use of the recordings of different field quantities (such as particle velocity and pressure). In practice, different recordings can be obscured by different sensor characteristics, requiring calibration with an unknown calibration factor. Not all field quantities required for multi-component field decomposition might be available, or they can suffer from different noise levels. The multi-depth-level decomposition approach makes use of field quantities recorded at multiple depth levels, e.g., two horizontal boreholes closely separated from each other, a combination of a single receiver array combined with free-surface boundary conditions, or acquisition geometries with a high-density of vertical boreholes. We theoretically describe the multi-depth-level decomposition approach in a unified form, showing that it can be applied to different kinds of fields in dissipative, inhomogeneous, anisotropic media, e.g., acoustic, electromagnetic, elastodynamic, poroelastic, and seismoelectric fields. We express the one-way fields at one depth level in terms of the observed fields at multiple depth levels, using extrapolation operators that are dependent on the medium parameters between the two depth levels. Lateral invariance at the depth level of decomposition allows us to carry out the multi-depth-level decomposition in the horizontal wavenumber–frequency domain. We illustrate the multi-depth-level decomposition scheme using two synthetic elastodynamic examples. The first example uses particle velocity recordings at two depth levels, whereas the second example combines recordings at one depth level with the Dirichlet free-surface boundary condition of zero traction. Comparison with multi-component decomposed fields shows a perfect match in both amplitude and phase for both cases. The multi-depth-level decomposition scheme is fully customizable to the desired acquisition geometry. The decomposition problem is in principle an inverse problem. Notches may occur at certain frequencies, causing the multi-depth-level composition matrix to become uninvertible, requiring additional notch filters. We can add multi-depth-level free-surface boundary conditions as extra equations to the multi-component composition matrix, thereby overdetermining this inverse problem. The combined multi-component–multi-depth-level decomposition on a land data set clearly shows improvements in the decomposition results, compared with the performance of the multi-component decomposition scheme.

Description: Seismic waves converted from compressional to shear mode in the shallow subsurface can be useful not only for obtaining shear-wave velocity information but also for improved processing of deeper reflection data. These waves generated at deep seas have been used successfully in hydrocarbon exploration; however, acquisition of good-quality converted-wave data in shallow marine environments remains challenging. We have looked into this problem through field experiments and synthetic modeling. A high-resolution seismic survey was conducted in a shallow-water canal using different types of seismic sources; data were recorded with a four-component water-bottom cable. Observed events in the field data were validated through modeling studies. Compressional waves converted to shear waves at the water bot-tom and at shallow reflectors were identified. The shear waves showed distinct linear polarization in the horizontal plane and low velocities in the marine sediments. Modeling results indicated the presence of a nongeometric shear-wave arrival excited only when the dominant wavelength exceeded the height of the source with respect to the water/sediment interface, as observed in air-gun data. This type of shear wave has a traveltime that corresponds to the raypath originating not at the source but at the interface directly below the source. Thus, these shear waves, excited by the source/water-bottom coupled system, kinematically behave as if they were generated by an S-wave source placed at the water bottom. In a shallow-water environment, the condition appears to be favorable for exciting such shear waves with nongeometric arrivals. These waves can provide useful information of shear-wave velocity in the sediments.