Borehole seismic‐source radiation pattern in transversely isotropic media

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 29-42 ◽  
Author(s):  
Wenjie Dong ◽  
M. Nafi Toksöz

We extend previous discussions on crosswell tomography in anisotropic formations by deriving the radiation patterns of three typical downhole seismic sources (impulsive air gun or dynamite, wall‐clamped vertical vibrators, and cylindrical bender) inside a fluid‐filled borehole embedded in a transversely isotropic (TI) formation. The method of steepest descents, in conjuncture with the low‐frequency and far‐field assumptions, is applied to the exact displacement integrals of these sources to obtain their radiation patterns asymptotically. In spite of complications caused by quasi‐P‐ and quasi‐SV‐wave coupling and wavefront triplication in homogeneous TI media, the final results can still be expressed in slowness components determined by a ray direction, which is desired when source radiation effects are to be accounted for by ray‐based tomography techniques. Tests with the radiation patterns show that while the effect of anisotropy on P‐waves is moderate, its effect on the S‐wave pattern is significant even for slightly anisotropic formations. One can predict the S‐wave pattern from the sign of the Thomsen’s measure δ*.

Geophysics ◽  
1994 ◽  
Vol 59 (11) ◽  
pp. 1780-1785 ◽  
Author(s):  
Richard L. Gibson ◽  
Chengbin Peng

An accurate characterization of borehole seismic sources is necessary to model and interpret waveforms observed in crosshole and reverse vertical seismic profiling (VSP) surveys, since the radiation pattern of a source will directly influence the amplitudes of elastic wave arrivals at receiver locations. Any attempt to study these data or perform inversions of amplitude data without incorporating the borehole effects will have serious limitations. Most previous studies of borehole seismic source radiation patterns have applied low‐frequency approximations to develop expressions for the radiation patterns of volume injection or stress sources (Heelan, 1953; White, 1960; White and Senghush, 1963; Lee and Balch, 1982; Lee, 1986; Kurkjian, 1986; Meredith, 1990; Winbow, 1991; Ben‐Menahem and Kostek, 1991). For example, Lee and Balch (1982) used this approach, along with a steepest descent solution, to derive closed‐form analytic expressions for the asymptotic far‐field radiation from sources located in uncased boreholes. Meredith (1990) applied the same methodology to study the radiation patterns of a variety of types of sources, though he also computed full waveform synthetic seismograms using the discrete wavenumber method. Likewise, Greenfield (1978) used full waveform numerical solutions to compute seismograms for force sources applied to the wall of a cylindrical cavity.


2022 ◽  
Vol 41 (1) ◽  
pp. 9-18
Author(s):  
Andrew Brenders ◽  
Joe Dellinger ◽  
Imtiaz Ahmed ◽  
Esteban Díaz ◽  
Mariana Gherasim ◽  
...  

The promise of fully automatic full-waveform inversion (FWI) — a (seismic) data-driven velocity model building process — has proven elusive in complex geologic settings, with impactful examples using field data unavailable until recently. In 2015, success with FWI at the Atlantis Field in the U.S. Gulf of Mexico demonstrated that semiautomatic velocity model building is possible, but it also raised the question of what more might be possible if seismic data tailor-made for FWI were available (e.g., with increased source-receiver offsets and bespoke low-frequency seismic sources). Motivated by the initial value case for FWI in settings such as the Gulf of Mexico, beginning in 2007 and continuing into 2021 BP designed, built, and field tested Wolfspar, an ultralow-frequency seismic source designed to produce seismic data tailor-made for FWI. A 3D field trial of Wolfspar was conducted over the Mad Dog Field in the Gulf of Mexico in 2017–2018. Low-frequency source (LFS) data were shot on a sparse grid (280 m inline, 2 to 4 km crossline) and recorded into ocean-bottom nodes simultaneously with air gun sources shooting on a conventional dense grid (50 m inline, 50 m crossline). Using the LFS data with FWI to improve the velocity model for imaging produced only incremental uplift in the subsalt image of the reservoir, albeit with image improvements at depths greater than 25,000 ft (approximately 7620 m). To better understand this, reprocessing and further analyses were conducted. We found that (1) the LFS achieved its design signal-to-noise ratio (S/N) goals over its frequency range; (2) the wave-extrapolation and imaging operators built into FWI and migration are very effective at suppressing low-frequency noise, so that densely sampled air gun data with a low S/N can still produce useable model updates with low frequencies; and (3) data density becomes less important at wider offsets. These results may have significant implications for future acquisition designs with low-frequency seismic sources going forward.


Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1471-1480 ◽  
Author(s):  
Thomas M. Daley ◽  
Dale Cox

A recently developed borehole seismic source, the orbital vibrator, was successfully deployed in a crosswell survey in a fractured basalt aquifer. This seismic source uses a rotating eccentric mass to generate seismic energy. Source sweeps with clockwise and counter‐clockwise rotations are recorded at each source location. Because this source generates circularly polarized waves, unique processing algorithms are used to decompose the recordings into two equivalent linearly oscillating, orthogonally oriented seismic sources. The orbital vibrator therefore generates P‐ and S‐waves simultaneously for all azimuths. A coordinate rotation based on P‐wave particle motion is used to align the source components from various depths. In a field experiment, both P‐ and S‐wave arrivals were recorded using fluid‐coupled hydrophone sensors. The processed field data show clear separation of P‐ and S‐wave arrivals for in‐line and crossline source components, respectively. A tensor convolutional description of the decomposition process allows for extension to multicomponent sensors.


Geophysics ◽  
1995 ◽  
Vol 60 (3) ◽  
pp. 748-757 ◽  
Author(s):  
Wenjie Dong ◽  
M. Nafi Toksöz

The source and receiver boreholes in crosshole seismology are usually considered unimportant except for their effects on body wave radiation and reception patterns. We present counter examples by analyzing a real crosswell data set from Buckhorn, Illinois, using computer simulations. The algorithm used is a combination of the boundary element method (for the source borehole) and the borehole coupling theory (for the receiver borehole) in transversely isotropic media. We find that most of the strong events in the data are inexplicable unless both boreholes are included in the modeling. The importance of the boreholes stems from the local geology which consists of highly contrasted sedimentary rocks. At a high‐contrast interface, wave conversion is no longer a negligible secondary effect. In fact, converted waves can be stronger than the primaries.


Geophysics ◽  
1995 ◽  
Vol 60 (3) ◽  
pp. 735-747 ◽  
Author(s):  
Wenjie Dong ◽  
Michel Bouchon ◽  
M. Nafi Toksöz

In modeling waves radiated from a borehole seismic source in layered isotropic or anisotropic media, the commonly used numerical methods (e.g., finite difference and finite element) encounter difficulties because of the large scale difference between the borehole diameter and the formation extent. To get around this problem, we apply the indirect boundary element method to establish a general algorithm for modeling source radiation from open and cased boreholes in layered transversely isotropic (TI) media. The essence of the algorithm is to use discrete secondary sources (unknowns) on both sides of the borehole wall (formation/cement interface) to represent the influence of the interface on wave scattering, so that wave propagation inside and outside the borehole can be carried out by Green’s functions. The discrete distribution of the secondary sources is determined by matching boundary conditions on the borehole wall. Comparison with the discrete wavenumber method validates the implementation. Applications to fluid‐filled open and cased boreholes in three‐layer media demonstrate the creation of guided modes in low velocity layers. Presence of anisotropy complicates the guided modes as a result of dispersion and P‐ and SV‐waves coupling in homogeneous TI media. Presence of casing and cement enhances the visibility of the guided modes.


2019 ◽  
Vol 38 (9) ◽  
pp. 680-690 ◽  
Author(s):  
Benoît Teyssandier ◽  
John J. Sallas

Ten years ago, CGG launched a project to develop a new concept of marine vibrator (MV) technology. We present our work, concluding with the successful acquisition of a seismic image using an ocean-bottom-node 2D survey. The expectation for MV technology is that it could reduce ocean exposure to seismic source sound, enable new acquisition solutions, and improve seismic data quality. After consideration of our objectives in terms of imaging, productivity, acoustic efficiency, and operational risk, we developed two spectrally complementary prototypes to cover the seismic bandwidth. In practice, an array composed of several MV units is needed for images of comparable quality to those produced from air-gun data sets. Because coupling to the water is invariant, MV signals tend to be repeatable. Since far-field pressure is directly proportional to piston volumetric acceleration, the far-field radiation can be well controlled through accurate piston motion control. These features allow us to shape signals to match precisely a desired spectrum while observing equipment constraints. Over the last few years, an intensive validation process was conducted at our dedicated test facility. The MV units were exposed to 2000 hours of in-sea testing with only minor technical issues.


Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 2082-2091 ◽  
Author(s):  
Bjørn Ursin ◽  
Ketil Hokstad

Compensation for geometrical spreading is important in prestack Kirchhoff migration and in amplitude versus offset/amplitude versus angle (AVO/AVA) analysis of seismic data. We present equations for the relative geometrical spreading of reflected and transmitted P‐ and S‐wave in horizontally layered transversely isotropic media with vertical symmetry axis (VTI). We show that relatively simple expressions are obtained when the geometrical spreading is expressed in terms of group velocities. In weakly anisotropic media, we obtain simple expressions also in terms of phase velocities. Also, we derive analytical equations for geometrical spreading based on the nonhyperbolic traveltime formula of Tsvankin and Thomsen, such that the geometrical spreading can be expressed in terms of the parameters used in time processing of seismic data. Comparison with numerical ray tracing demonstrates that the weak anisotropy approximation to geometrical spreading is accurate for P‐waves. It is less accurate for SV‐waves, but has qualitatively the correct form. For P waves, the nonhyperbolic equation for geometrical spreading compares favorably with ray‐tracing results for offset‐depth ratios less than five. For SV‐waves, the analytical approximation is accurate only at small offsets, and breaks down at offset‐depth ratios less than unity. The numerical results are in agreement with the range of validity for the nonhyperbolic traveltime equations.


Sign in / Sign up

Export Citation Format

Share Document