1 

Hydraulic fracture characterization with dispersion measurements of seismic waves

[PDF]

2 

Experimental verification of stressinduced anisotropy

[PDF]

3 

Biangular decomposition of seismic data

[PDF]

4 

A theoretical and experimental approach to the geophoneground coupling problem based on acoustic reciprocity

[PDF]

5 

A new method to convert unleveled marine seismic data to leveled splitspread data

[PDF]

6 

The reflectivity operator for curved interfaces

[PDF]

7 

Nonlinear imaging condition and the effect of source illumination: Imaging fractures as nonwelded interfaces
Fluid flow through a fractured reservoir is often controlled by the large fractures. Seismically imaging these large fractures has the potential to illuminate their hydraulic properties.We derived a nonlinear imaging condition considering a medium containing nonwelded interfaces such as fractures for highresolution fracture imaging. This was achieved by using the general correlationtype representation theorem relating the wavefield between two different states representing different fracture compliances. We numerically tested the imaging condition to investigate the effect of the nonlinear term and that of the onesided source illumination. Assuming a dry fracture, we calculated the wavefield from a nonwelded interface. We obtained a Pwave imaging result from Pwave sources. In the case of perfect source illumination, we found that introducing the nonlinear term in the imaging condition enhances the image because the nonwelded interface is imaged as a thin layer in an otherwise homogeneous medium. Investigating the effect of onesided illumination by horizontally aligned sources revealed interesting limitations and possibilities. The imaging result for a horizontal fracture showed a volumetric distribution of nonzero amplitudes around the polarity change at the fracture that can be misinterpreted as a welded thick layer boundary. However, when the fracture was not horizontal, the imaging result was quite good and was closer to that with a perfect source illumination. These led to a new possibility of imaging subvertical fractures from surface seismic measurements, or subhorizontal fractures from vertical seismic profiling data, assuming that we successfully estimated the perturbed wavefield from the receiver responses.

[PDF]
[Abstract]

8 

A seismic vertical vibrator driven by linear synchronous motors
A linear synchronous motor (LSM) is an electric motor that can produce large controllable forces and is therefore suitable as a driving engine for a seismic vibrator. This motor consists of two independent elements, a magnet track and a coil track, allowing practically unlimited motor displacements. This makes the LSM very suitable for expanding the source frequency band to the lower frequencies in which larger strokes are needed. In contrast to hydraulic engines, the LSM performs equally well over the whole frequency range, making possible a smaller amount of signal distortion, especially at the low frequencies. To find the feasibility of an LSMdriven vibrator, we successfully designed and built a multiLSM prototype vibrator of some 1200 kg. We addressed the synchronization between the individual motor tracks and the different motors. To lower the energy consumption, a spring mechanism was implemented that delivered the force needed to lift the vibrator mass to its neutral position. The resonance belonging to this spring mechanism was successfully suppressed with the help of a position feedback control that also suppressed the temperature effects. The seismic data acquired in the field tests proved that the prototype LSM vibrator acted very well as a seismic source. It has no trouble generating pseudorandom sweeps, and even given its limited size, it generated signals within the lowfrequency regime, down to 2 Hz, rather easily.

[PDF]
[Abstract]

9 

Studying CO2 storage with ambientnoise seismic interferometry: A combined numerical feasibility study and fielddata example for Ketzin, Germany
Seismic interferometry applied to ambientnoise measurements allows the retrieval of the seismic response between pairs of receivers. We studied ambientnoise seismic interferometry (ANSI) to retrieve timelapse reflection responses from a reservoir during CO2 geologic sequestration, using the case of the experimental site of Ketzin, Germany. We applied ANSI to numerically modeled data to retrieve base and repeat reflection responses characterizing the impedances occurring at the reservoir both with and without the injection of CO2. The modeled data represented global transmission responses from bandlimited noise sources randomly triggered in space and time. We found that strong constraints on the spatial distribution of the passive sources were not required to retrieve the timelapse signal as long as sufficient sourcelocation repeatability was observed between the base and the repeat passive survey. To illustrate the potential of the technique, ANSI was applied to three days of passive field data recorded in 2012 at Ketzin. Comparison with the modeled results illustrated the potential to retrieve key reflection events using ANSI on field data from Ketzin. This study supports the idea that the geologic setting and characteristics of ambient noise at Ketzin may be opportune to monitor CO2 sequestration.

[PDF]
[Abstract]

10 

An approximate 3D computational method for realtime computation of induction logging responses
Over many years, induction logging systems have been used to create well formation logs. The major drawback for the utilization of these tools is the long simulation time for a single forward computation. We proposed an efficient computational method based on a contrasttype of integralequation formulation, in which we applied an approximation for the 3D electromagnetic field. We assumed that the dominant contribution in the integral equation is obtained by the contribution around the singularity of Green’s kernel. It is expected that the approximation yields reliable results when the (homogeneous) background conductivity around the logging tool is close to the actual conductivity at the location of the tool. We have developed a datadriven method to determine this background conductivity from the dominant part of the measured coaxial magnetic fields, which are mainly influenced by the conductivity at the tool sensors. For a synthetic model, the results were compared to the ones of a rigorous solution of the integral equation and show a good simulation response to smallscale variations in the medium. Further, the method was used to simulate the response of a realistic reservoir model. Such a model is created by a geological modeling program. We concluded that our approximate method was able to improve the approximation results in highly heterogeneous structures compared to the Born approximation and provide an effective mediumgradient around the tool. Our method, based on the wavefield approximation, also estimates the error, and hence yields a warning when the method becomes unreliable.

[PDF]
[Abstract]

11 

Retrieval of reflections from ambient noise recorded in the Mizil area, Romania
We applied seismic interferometry (SI) by crosscorrelation to ambientnoise panels recorded in the Mizil area, Romania, aiming to retrieve bodywave reflections. To achieve this goal, surface waves in the noise panels input to SI should be suppressed. We did this by selecting for input to SIonly noise panels that are not dominated by surface waves; the selection was either after visual inspection in the time domain or after automatic slowness evaluation on crosscorrelated panels. The latter used the slowness of arrivals passing through the virtualsource position at time 0 s. We discovered that the automatic slownessevaluation method allows better retrieval of reflections. From the retrieved reflection gathers, we obtained stacked images of the subsurface. Comparing the SI images to a stacked image from activesource data, we concluded that some retrieved events correspond to reflectors in the active seismic section, including known geologic markers. In a previous application of SI to ambient noise, the retrieved reflections exhibited frequency content lower than that of active data. In our results, the frequency content of the SI retrieved data was comparable to the one of the active data.

[PDF]
[Abstract]

12 

Datadriven wavefield focusing and imaging with multidimensional deconvolution: Numerical examples for reflection data with internal multiples
Standard imaging techniques rely on the single scattering assumption. This requires that the recorded data do not include internal multiples, i.e., waves that have bounced multiple times between reflectors before reaching the receivers at the acquisition surface. When multiple reflections are present in the data, standard imaging algorithms incorrectly image them as ghost reflectors. These artifacts can mislead interpreters in locating potential hydrocarbon reservoirs. Recently, we introduced a new approach for retrieving the Green’s function recorded at the acquisition surface due to a virtual source located at depth. We refer to this approach as datadriven wavefield focusing. Additionally, after applying sourcereceiver reciprocity, this approach allowed us to decompose the Green’s function at a virtual receiver at depth in its downgoing and upgoing components. These wavefields were then used to create a ghostfree image of the medium with either crosscorrelation or multidimensional deconvolution, presenting an advantage over standard prestack migration. We tested the robustness of our approach when an erroneous background velocity model is used to estimate the firstarriving waves, which are a required input for the datadriven wavefield focusing process. We tested the new method with a numerical example based on a modification of the Amoco model.

[PDF]
[Abstract]

13 

Marchenko imaging
Traditionally, the Marchenko equation forms a basis for 1D inverse scattering problems. A 3D extension of the Marchenko equation enables the retrieval of the Green’s response to a virtual source in the subsurface from reflection measurements at the earth’s surface. This constitutes an important step beyond seismic interferometry. Whereas seismic interferometry requires a receiver at the position of the virtual source, for the Marchenko scheme it suffices to have sources and receivers at the surface only. The underlying assumptions are that the medium is lossless and that an estimate of the direct arrivals of the Green’s function is available. The Green’s function retrieved with the 3D Marchenko scheme contains accurate internal multiples of the inhomogeneous subsurface. Using sourcereceiver reciprocity, the retrieved Green’s function can be interpreted as the response to sources at the surface, observed by a virtual receiver in the subsurface. By decomposing the 3D Marchenko equation, the response at the virtual receiver can be decomposed into a downgoing field and an upgoing field. By deconvolving the retrieved upgoing field with the downgoing field, a reflection response is obtained, with virtual sources and virtual receivers in the subsurface. This redatumed reflection response is free of spurious events related to internal multiples in the overburden. The redatumed reflection response forms the basis for obtaining an image of a target zone. An important feature is that spurious reflections in the target zone are suppressed, without the need to resolve first the reflection properties of the overburden.

[PDF]
[Abstract]

14 

Combining intersource seismic interferometry and sourcereceiver interferometry for deep local imaging
The virtual source method has been applied successfully to retrieve the impulse response between pairs of receivers in the subsurface. This method is further improved by an updown separation prior to the crosscorrelation to suppress the reflections from the overburden and the free surface. In a reversed situation where the sources are in the subsurface and receivers are on the surface, in principle, one can apply the same logic to retrieve the virtual response between pairs of sources by sourcereceiver reciprocity, turning the physical borehole sources into virtual receivers. However, since the updown separation is not applicable on the source side, the simple crosscorrelation of the total fields results in spurious events due to the incomplete receiver coverage around the sources. We show with a numerical example that for this configuration of borehole sources and surface receivers, one can replace such an updown separation at the source side by that of the direct and reflected waves as a first order approximation. This procedure produces the virtual receiver data that is adequate for local imaging below the source depth and is completely independent of the accuracy of the overburden velocity model. We implement this intersource type of interferometry by multidimensional deconvolution (MDD). Further, if the conventional surface survey data is available, we test the methodology from sourcereceiver interferometry (SRI) for this reverse configuration with borehole sources to retrieve the virtual receiver data with reflections coming from above, using also only the separation of the direct and reflected waves. By migrating the two sets of virtual receiver data, one can create a local image around the borehole sources in a deep area with better focusing and localization without a sophisticated velocity model.

[PDF]
[Abstract]

15 

3D acquisition geometry analysis: Incorporating information from multiples
Recent advances in survey design have led to conventional commonmidpointbased analysis being replaced by the subsurfacebased seismic acquisition analysis and design, with the emphasis on advance techniques of illumination analysis. Amongst them are waveequationbased seismic illumination analyses such as the socalled focal beam method. The method's objective is to provide quantitative insight into the combined influence of acquisition geometry, overburden structure, and migration operators on image resolution and angledependent amplitude fidelity. So far, this method has only addressed illumination by primaries. However, multiples also contain information and may illuminate the target point from other angles than primaries, resulting in higher resolution and better illumination. In this paper, the focal beam analysis concept for multiplereflected waves is introduced. The method demonstrates how the acquisitionrelated amplitude footprint can be corrected and improved using multiples.

[PDF]
[Abstract]

16 

Characterization of a heterogeneous landfill using seismic and electrical resistivity data
Understanding the processes occurring inside a landfill is important for improving the treatment of landfills. Irrigation and recirculation of leachate are widely used in landfill treatments. Increasing the efficiency of such treatments requires a detailed understanding of the flow inside the landfill. The flow depends largely on the heterogeneous distribution of density. It is, therefore, of great practical interest to determine the density distribution affecting the flow paths inside a landfill. Studies in the past have characterized landfill sites but have not led to highresolution, detailed quantitative results. We performed an Swave reflection survey, multichannel analysis of surface waves (MASW), and electrical resistivity survey to investigate the possibility of delineating the heterogeneity distribution in the body of a landfill. We found that the highresolution Swave reflection method offers the desired resolution. However, in the case of a very heterogeneous landfill and a high noise level, the processing of highresolution, shallow reflection data required special care. In comparison, MASW gave the general trend of the changes inside the landfill, whereas the electrical resistivity (ER) survey provides useful clues for interpretation of seismic reflection data. We found that it is possible to localize finescale heterogeneities in the landfill using the Swave reflection method using a highfrequency vibratory source. Using empirical relations specific to landfill sites, we then estimated the density distribution inside the landfill, along with the associated uncertainty considering different methods. The final interpretation was guided by supplementary information provided by MASW and ER tomography.

[PDF]
[Abstract]

17 

Sensitivity of the nearsurface vertical electric field land ControlledSource Electromagnetic monitoring
We investigate potential benefits of measuring the vertical electric field component in addition to the routinely measured horizontal electric field components in onshore timelapse controlledsource electromagnetics. Synthetic electromagnetic data based on a model of the Schoonebeek onshore oil field are used. We confirm that the vertical electric field component is more sensitive to small changes in the reservoir than the horizontal components, yet its amplitudes are small. Accordingly, optimal sourcereceiver geometry and precise knowledge of the verticality of the receiver dipole will be required for successful utilization of the vertical electric field.

[PDF]
[Abstract]

18 

3D Surfacewave estimation and separation: A closedloop approach
Surfacewaves are often dominant in seismic data in a shallow water and land environment. Separating them out from the seismic data is of great importance for either removing them as noise for reservoir characterization, or extracting them as signal for nearsurface characterization. However, their complex properties make the surfacewave separation significantly challenging in seismic processing. To address the challenges and adopt recent advances in seismic processing, we propose a methodology of surfacewave estimation and separation using a closedloop approach. The methodology was successfully demonstrated on real 3D seismic data with mudroll that is often dominant in a shallow water environment.

[PDF]
[Abstract]

19 

The electromagnetic response in a layered vertical transverse isotropic medium: A new look at an old problem
We determined that the electromagnetic vertical transverse isotropic response in a layered earth can be obtained by solving two equivalent scalar equations, which were for the vertical electric field and for the vertical magnetic field, involving only a scalar global reflection coefficient. Besides the complete derivation of the full electromagnetic response, we also developed the corresponding computer code called EMmod, which models the full electromagnetic fields including internal multiples in the frequencywavenumber domain and obtains the frequencyspace domain solutions through a Hankel transformation by computing the Hankel integral using a 61point GaussKronrod integration routine. The code is able to model the 3D electromagnetic field in a 1D earth for diffusive methods such as controlled source electromagnetics as well as for wave methods such as ground penetrating radar. The user has complete freedom to place the source and the receivers in any layer. The modeling is illustrated with three examples, which aim to present the different capabilities of EMmod, while assessing its correctness.

[PDF]
[Abstract]

20 

An illustration of adaptive Marchenko imaging
In Marchenko imaging, wavefields are retrieved at specified focal points in the subsurface through an iterative scheme derived from the multidimensional Marchenko equation. The method requires seismicreflection data at the earth’s surface (after freesurface multiple elimination) and an estimate of the direct wavefield from the surface to each focal point, which can be computed, for instance, in a macrovelocity model. In the first iteration, the direct wavefield is crosscorrelated with the reflection data. This operation is identical to inversewavefield extrapolation as is applied commonly in various imaging schemes, for instance, in reverse time migration (RTM). At each succeeding iteration, the result of the previous iteration is truncated in time and crosscorrelated with the reflection data again. To obtain a seismic image, a multidimensional deconvolutionbased imaging condition can be applied to the retrieved wavefields. By this approach, both primary reflections and internal multiples contribute to the construction of the image. Alternatively, a crosscorrelationbased imaging condition can be used in which only the primary reflections are imaged and the contributions of internal multiples are subtracted. The latter strategy offers more flexibility because the subtraction of redatumed internal multiples can be implemented adaptively. Through this approach, the artifacts from internal multiples can be removed effectively from a conventional RTM image.

[PDF]
[Abstract]
