The full acoustic wave train in a laboratory model of a borehole

Geophysics ◽  
1982 ◽  
Vol 47 (11) ◽  
pp. 1512-1520 ◽  
Author(s):  
S. T. Chen
Keyword(s):  
Acoustic Wave ◽  
Wave Velocity ◽  
Theoretical Calculations ◽  
Wave Train ◽  
High Amplitude ◽  
P Wave ◽  
Laboratory Model ◽  
Stoneley Wave ◽  
S Wave ◽  
Full Wave

We studied the characteristics of acoustic wave propagation in a fluid-filled borehole using as a laboratory model a concrete cylinder 2 ft high and 2 ft in diameter with a 1/4-inch diameter borehole along its axis. The model represents sonic logging in the field reduced by a factor of 40. We recorded the full wave train consisting of a refracted compressional P wave, a refracted shear S wave, and guided waves including a number of normal modes and a Stoneley wave. Exploiting the dispersive properties of a modal wave and the source-receiver frequency characteristics, we were able to isolate the S–wave, which contains much valuable information about the formation rock, but which has not been widely used since it is difficult to extract from the full wave train. The observed Stoneley wave had a very high amplitude at low frequency and showed little dispersion. Stoneley-wave velocity is closely related to S–wave velocity and formation density, and can be measured very accurately because the Stoneley wave generally has high amplitude and low attenuation. It can therefore be used indirectly to obtain the S–wave velocity, even when the S–wave cannot be measured directly. In general, the observed characteristics of each component wave agreed with our theoretical calculations but their relative amplitudes did not. We believe these discrepancies were caused, in part, by the fact that rock attenuation and the latitudinal angular dependence of the source radiation were not taken into account in the theoretical calculations.

Shear‐wave logging with dipole sources

Geophysics ◽  
1988 ◽  
Vol 53 (5) ◽  
pp. 659-667 ◽  
Author(s):  
S. T. Chen
Keyword(s):  
Shear Wave ◽  
Wave Velocity ◽  
Wave Train ◽  
P Wave ◽  
Dipole Source ◽  
Stoneley Wave ◽  
S Wave ◽  
S Waves ◽  
On Line ◽  
First Arrival

Laboratory measurements have verified a novel technique for direct shear‐wave logging in hard and soft formations with a dipole source, as recently suggested in theoretical studies. Conventional monopole logging tools are not capable of measuring shear waves directly. In particular, no S waves are recorded in a soft formation with a conventional monopole sonic tool because there are no critically refracted S rays when the S-wave velocity of the rock is less than the acoustic velocity of the borehole fluid. The present studies were conducted in the laboratory with scale models representative of sonic logging conditions in the field. We have used a concrete model to represent hard formations and a plastic model to simulate a soft formation. The dipole source, operating at frequencies lower than those conventionally used in logging, substantially suppressed the P wave and excited a wave train whose first arrival traveled at the S-wave velocity. As a result, one can use a dipole source to log S-wave velocity directly on‐line by picking the first arrival of the full wave train, in a process similar to that used in conventional P-wave logging. Laboratory experiments with a conventional monopole source in a soft formation did not produce S waves. However, the S-wave velocity was accurately estimated by using Biot’s theory, which required measuring the Stoneley‐wave velocity and knowing other borehole parameters.

A theoretical study of acoustic S-wave and P-wave velocity logging with conventional and dipole sources in soft formations

Geophysics ◽  
1988 ◽  
Vol 53 (10) ◽  
pp. 1334-1342 ◽  
Author(s):  
Graham A. Winbow
Keyword(s):  
Wave Velocity ◽  
Water Waves ◽  
Fluid Velocity ◽  
Soft Materials ◽  
P Wave ◽  
Stoneley Wave ◽  
Drilling Mud ◽  
S Wave ◽  
Low Frequencies ◽  
P Wave Velocity

This paper is focused on the special features of the wavetrains recorded by conventional and dipole sonic logging tools in soft formations defined to be those whose shar velocity is less than the sound velocity of drilling mud. Such formations are commonn in the Gulf Coast, the Canadian Arctic, the Bass Strait of Australia, and many other region. A conventional logging tool operating at normal frequencies [Formula: see text] records P waves, water waves, and Stoneley waves in soft formations. A dipole tool records modal waves and water waves at frequencies of order 15 kHz, but produces almost pure S-wave first arrivals at low frequencies [Formula: see text] since at 1 kHz, a mode which we refer to as a “dipole Stoneley wave” is efficiently excited. For very soft materials such as clays, where the formation P-wave velocity can be less than the fluid velocity, the formation P velocity can be logged by operating a conventional sonic tool at low frequencies [Formula: see text] so as to excite a leaky mode traveling at very close to the formation P-wave velocity. Water waves are not important for high‐velocity formations where they arrive at the trailing edge of the modal part of the wavetrain. However, in soft formalions they form a prominent part of the wavetrain at normal logging frequencies [Formula: see text] and disappear at low frequencies [Formula: see text]. Water waves are carried by leaky modes.

Local and global fluid effects on sonic wave modes

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. D369-D381 ◽  
Author(s):  
Elliot J. H. Dahl ◽  
Kyle T. Spikes
Keyword(s):  
Rayleigh Wave ◽  
Wave Velocity ◽  
Fluid Pressure ◽  
Wave Mode ◽  
P Wave ◽  
Stoneley Wave ◽  
S Wave ◽  
Local Pressure ◽  
Wave Velocities ◽  
Wave Modes

Most subsurface formations of value to exploration contain a heterogeneous fluid-filled pore space, where local fluid-pressure effects can significantly change the velocities of passing seismic waves. To better understand the effect of these local pressure gradients on borehole wave propagation, we combined Chapman’s squirt-flow model with Biot’s poroelastic theory. We applied the unified theory to a slow and fast formation with permeable borehole walls containing different quantities of compliant pores. These results are compared with those for a formation with no soft pores. The discrete wavenumber summation method with a monopole point source generates the wavefields consisting of the P-, S-, leaky-P, Stoneley, and pseudo-Rayleigh waves. The resulting synthetic wave modes are processed using a weighted spectral semblance (WSS) algorithm. We found that the resulting WSS dispersion curves closely matched the analytical expressions for the formation compressional velocity and solutions to the period equation for dispersion for the P-wave, Stoneley-wave, and pseudo-Rayleigh wave phase velocities in the slow and fast formations. The WSS applied to the S-wave part of the waveforms, however, did not correlate as well with its respective analytical expression for formation S-wave velocity, most likely due to interference of the pseudo-Rayleigh wave. To separate changes in formation P- and S-wave velocities versus fluid-flow effects on the Stoneley-wave mode, we computed the slow-P wave dispersion for the same formations. We found that fluid-saturated soft pores significantly affected the P- and S-wave effective formation velocities, whereas the slow-P wave velocity was rather insensitive to the compliant pores. Thus, the large phase-velocity effect on the Stoneley wave mode was mainly due to changes in effective formation P- and S-wave velocities and not to additional fluid mobility.

Detection of gas in sandstone reservoirs using AVO analysis: A 3-D seismic case history using the Geostack technique

Geophysics ◽  
1994 ◽  
Vol 59 (9) ◽  
pp. 1362-1376 ◽  
Author(s):  
Jan L. Fatti ◽  
George C. Smith ◽  
Peter J. Vail ◽  
Peter J. Strauss ◽  
Philip R. Levitt
Keyword(s):  
Wave Velocity ◽  
High Amplitude ◽  
Igneous Rocks ◽  
P Wave ◽  
S Wave ◽  
Gas Field ◽  
Data Set ◽  
Rock Types ◽  
Fluid Factor ◽  
Clastic Sedimentary

The Geostack technique is a method of analyzing seismic amplitude variation with offset (AVO) information. One of the outputs of the analysis is a set of direct hydrocarbon indicator traces called “fluid factor” traces. The fluid factor trace is designed to be low amplitude for all reflectors in a clastic sedimentary sequence except for rocks that lie off the (mudrock line.) The mudrock line is the line on a crossplot of P‐wave velocity against S‐wave velocity on which water‐saturated sandstones, shales, and siltstones lie. Some of the rock types that lie off the mudrock line are gas‐saturated sandstones, carbonates, and igneous rocks. In the absence of carbonates and igneous rocks, high amplitude reflections on fluid factor traces would be expected to represent gas‐saturated sandstones. Of course, this relationship does not apply exactly in nature, and the extent to which the mudrock line model applies varies from area to area. However, it is a useful model in many basins of the world, including the one studied here. Geostack processing has been done on a 3-D seismic data set over the Mossel Bay gas field on the southern continental shelf of South Africa. We found that anomalously high amplitude fluid factor reflections occurred at the top and base of the gas‐reservoir sandstone. Maps were made of the amplitude of these fluid factor reflections, and it was found that the high amplitude values were restricted mainly to the gas field area as determined by drilling. The highest amplitudes were found to be located roughly in the areas of best reservoir quality (i.e., highest porosity) in areas where the reservoir is relatively thick.

Sonic logging in deviated boreholes penetrating an anisotropic formation: Laboratory study

Geophysics ◽  
2007 ◽  
Vol 72 (4) ◽  
pp. E125-E134 ◽  
Author(s):  
Zhenya Zhu ◽  
Shihong Chi ◽  
M. Nafi Toksöz
Keyword(s):  
Transversely Isotropic ◽  
Body Waves ◽  
P Wave ◽  
Laboratory Model ◽  
Stoneley Wave ◽  
S Wave ◽  
Wave Velocities ◽  
Anisotropic Formation ◽  
Sonic Logging ◽  
Drill Holes

Development of deepwater fields requires drilling deviated or horizontal wells. Many formations are highly anisotropic, that is, the P- and S-wave velocities vary with propagation direction. Sonic logs acquired in these wells need to be corrected for anisotropy effects before the logs can be used in formation evaluation and seismic applications. In this study, we use a laboratory model made of an orthorhombic Phenolite block to study acoustic logging in deviated wells. We first measure the qP-, qSV-, and SH-wave group velocities by using body waves at angles of 0°, 15°, 30°, 45°, 60°, 75°, and 90° relative to the slowest P-wave principal axis of the Phenolite block. We then drill holes at the same angles in the block. We record monopole and dipole sonic waveforms in these holes and extract the qP-, qSV-, SH-, and Stoneley-wave velocities by using the slowness-time semblance method. The velocities measured through the use of monopole logging and dipole logging vary with borehole deviations. We find that an equivalent transversely isotropic (TI) model can fit the measured qP-, qSV-, and Stoneley-wave velocities very well. The S-wave velocities at low to medium borehole deviations can be used to differentiate an orthorhombic material from a TI one.

A reality check on full-wave inversion applied to land seismic data for near-surface modeling

2022 ◽  
Vol 41 (1) ◽  
pp. 40-46
Author(s):  
Öz Yilmaz ◽  
Kai Gao ◽  
Milos Delic ◽  
Jianghai Xia ◽  
Lianjie Huang ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
Surface Modeling ◽  
P Wave ◽  
Borehole Data ◽  
Near Surface ◽  
Traveltime Tomography ◽  
Wave Inversion ◽  
Full Wave ◽  
Elastic Inversion

We evaluate the performance of traveltime tomography and full-wave inversion (FWI) for near-surface modeling using the data from a shallow seismic field experiment. Eight boreholes up to 20-m depth have been drilled along the seismic line traverse to verify the accuracy of the P-wave velocity-depth model estimated by seismic inversion. The velocity-depth model of the soil column estimated by traveltime tomography is in good agreement with the borehole data. We used the traveltime tomography model as an initial model and performed FWI. Full-wave acoustic and elastic inversions, however, have failed to converge to a velocity-depth model that desirably should be a high-resolution version of the model estimated by traveltime tomography. Moreover, there are significant discrepancies between the estimated models and the borehole data. It is understandable why full-wave acoustic inversion would fail — land seismic data inherently are elastic wavefields. The question is: Why does full-wave elastic inversion also fail? The strategy to prevent full-wave elastic inversion of vertical-component geophone data trapped in a local minimum that results in a physically implausible near-surface model may be cascaded inversion. Specifically, we perform traveltime tomography to estimate a P-wave velocity-depth model for the near-surface and Rayleigh-wave inversion to estimate an S-wave velocity-depth model for the near-surface, then use the resulting pairs of models as the initial models for the subsequent full-wave elastic inversion. Nonetheless, as demonstrated by the field data example here, the elastic-wave inversion yields a near-surface solution that still is not in agreement with the borehole data. Here, we investigate the limitations of FWI applied to land seismic data for near-surface modeling.

Numerical and Theoretical Investigation on the Origin of the Multiple Peaks of the H/V Ratio Curve

2021 ◽  
Author(s):  
Wanbo Xiao ◽  
Siqi Lu ◽  
Yanbin Wang
Keyword(s):  
Rayleigh Wave ◽  
Wave Velocity ◽  
Subsurface Layer ◽  
P Wave ◽  
Theoretical Formula ◽  
S Wave ◽  
Multiple Peaks ◽  
Ratio Curve ◽  
Multiple Layer ◽  
Wave Resonance

<p>Despite the popularity of the horizontal to vertical spectral ratio (HVSR) method in site effect studies, the origin of the H/V peaks has been controversial since this method was proposed. Many previous studies mainly focused on the explanation of the first or single peak of the H/V ratio, trying to distinguish between the two hypotheses — the S-wave resonance and ellipticity of Rayleigh wave. However, it is common both in numerical simulations and practical experiments that the H/V ratio exhibits multiple peaks, which is essential to explore the origin of the H/V peaks.</p><p>The cause for the multiple H/V peaks has not been clearly figured out, and once was simply explained as the result of multi subsurface layers. Therefore, we adopted numerical method to simulate the ambient noise in various layered half-space models and calculated the H/V ratio curves for further comparisons. The peak frequencies of the H/V curves accord well with the theoretical frequencies of S-wave resonance in two-layer models, whose frequencies only depend on the S wave velocity and the thickness of the subsurface layer. The same is true for models with varying model parameters. Besides, the theoretical formula of the S-wave resonance in multiple-layer models is proposed and then supported by numerical investigations as in the cases of two-layer models. We also extended the S-wave resonance to P-wave resonance and found that its theoretical frequencies fit well with the V/H peaks, which could be an evidence to support the S-wave resonance theory from a new perspective. By contrast, there are obvious differences between the higher orders of the H/V ratio peaks and the higher orders of Rayleigh wave ellipticity curves both in two-layer and multiple-layer models. The Rayleigh wave ellipticity curves are found to be sensitive to the Poisson’s ratio and the thickness of the subsurface layer, so the variation of the P wave velocity can affect the peak frequencies of the Rayleigh wave ellipticity curves while the H/V peaks show slight change. The Rayleigh wave ellipticity theory is thus proved to be inappropriate for the explanation of the multiple H/V peaks, while the possible effects of the Rayleigh wave on the fundamental H/V peak still cannot be excluded.</p><p>Based on the analyses above, we proposed a new evidence to support the claim that the peak frequencies of the H/V ratio curve, except the fundamental peaks, are caused by S-wave resonance. The relationship between the P-wave resonance and the V/H peaks may also find further application.</p>

Using Seismic Data to Predict Shale Pore Pressure and Overpressure Sweet Spots: A Case Study from the Lower Silurian Longmaxi Formation in Sichuan Basin, China

2021 ◽  
Author(s):  
Sheng Chen ◽  
Qingcai Zeng ◽  
Xiujiao Wang ◽  
Qing Yang ◽  
Chunmeng Dai ◽  
...  
Keyword(s):  
Prediction Model ◽  
Wave Velocity ◽  
Shale Gas ◽  
Seismic Data ◽  
P Wave ◽  
Formation Pressure ◽  
S Wave ◽  
P Wave Velocity ◽  
New Formation ◽  
Sweet Spots

Abstract Practices of marine shale gas exploration and development in south China have proved that formation overpressure is the main controlling factor of shale gas enrichment and an indicator of good preservation condition. Accurate prediction of formation pressure before drilling is necessary for drilling safety and important for sweet spots predicting and horizontal wells deploying. However, the existing prediction methods of formation pore pressures all have defects, the prediction accuracy unsatisfactory for shale gas development. By means of rock mechanics analysis and related formulas, we derived a formula for calculating formation pore pressures. Through regional rock physical analysis, we determined and optimized the relevant parameters in the formula, and established a new formation pressure prediction model considering P-wave velocity, S-wave velocity and density. Based on regional exploration wells and 3D seismic data, we carried out pre-stack seismic inversion to obtain high-precision P-wave velocity, S-wave velocity and density data volumes. We utilized the new formation pressure prediction model to predict the pressure and the spatial distribution of overpressure sweet spots. Then, we applied the measured pressure data of three new wells to verify the predicted formation pressure by seismic data. The result shows that the new method has a higher accuracy. This method is qualified for safe drilling and prediction of overpressure sweet spots for shale gas development, so it is worthy of promotion.

Constraints on the composition of the crust and uppermost mantle in northwestern Canada: Vp/Vs variations along Lithoprobe's SNorCLE transect

2005 ◽  
Vol 42 (6) ◽  
pp. 1205-1222 ◽  
Author(s):  
Gabriela Fernández-Viejo ◽  
Ron M Clowes ◽  
J Kim Welford
Keyword(s):  
Wave Velocity ◽  
Upper Mantle ◽  
Laboratory Data ◽  
High Ratio ◽  
P Wave ◽  
S Wave ◽  
Uppermost Mantle ◽  
Crust And Upper Mantle ◽  
Slave Craton ◽  
Crustal Composition

Shear-wave seismic data recorded along four profiles during the SNoRE 97 (1997 Slave – Northern Cordillera Refraction Experiment) refraction – wide-angle reflection experiment in northwestern Canada are analyzed to provide S-wave velocity (Vs) models. These are combined with previous P-wave velocity (Vp) models to produce cross sections of the ratio Vp/Vs for the crust and upper mantle. The Vp/Vs values are related to rock types through comparisons with published laboratory data. The Slave craton has low Vp/Vs values of 1.68–1.72, indicating a predominantly silicic crustal composition. Higher values (1.78) for the Great Bear and eastern Hottah domains of the Wopmay orogen imply a more mafic than average crustal composition. In the western Hottah and Fort Simpson arc, values of Vp/Vs drop to ∼1.69. These low values continue westward for 700 km into the Foreland and Omineca belts of the Cordillera, providing support for the interpretation from coincident seismic reflection studies that much of the crust from east of the Cordilleran deformation front to the Stikinia terrane of the Intermontane Belt consists of quartzose metasedimentary rocks. Stikinia shows values of 1.78–1.73, consistent with its derivation as a volcanic arc terrane. Upper mantle velocity and ratio values beneath the Slave craton indicate an ultramafic peridotitic composition. In the Wopmay orogen, the presence of low Vp/Vs ratios beneath the Hottah – Fort Simpson transition indicates the presence of pyroxenite in the upper mantle. Across the northern Cordillera, low Vp values and a moderate-to-high ratio in the uppermost mantle are consistent with the region's high heat flow and the possible presence of partial melt.

Nonlinear two‐dimensional elastic inversion of multioffset seismic data

Geophysics ◽  
1987 ◽  
Vol 52 (9) ◽  
pp. 1211-1228 ◽  
Author(s):  
Peter Mora
Keyword(s):  
Wave Equation ◽  
Wave Velocity ◽  
A Priori ◽  
Gaussian Model ◽  
Wave Impedance ◽  
P Wave ◽  
S Wave ◽  
Elastic Wave Equation ◽  
Gradient Direction ◽  
P Wave Velocity

The treatment of multioffset seismic data as an acoustic wave field is becoming increasingly disturbing to many geophysicists who see a multitude of wave phenomena, such as amplitude‐offset variations and shearwave events, which can only be explained by using the more correct elastic wave equation. Not only are such phenomena ignored by acoustic theory, but they are also treated as undesirable noise when they should be used to provide extra information, such as S‐wave velocity, about the subsurface. The problems of using the conventional acoustic wave equation approach can be eliminated via an elastic approach. In this paper, equations have been derived to perform an inversion for P‐wave velocity, S‐wave velocity, and density as well as the P‐wave impedance, S‐wave impedance, and density. These are better resolved than the Lamé parameters. The inversion is based on nonlinear least squares and proceeds by iteratively updating the earth parameters until a good fit is achieved between the observed data and the modeled data corresponding to these earth parameters. The iterations are based on the preconditioned conjugate gradient algorithm. The fundamental requirement of such a least‐squares algorithm is the gradient direction which tells how to update the model parameters. The gradient direction can be derived directly from the wave equation and it may be computed by several wave propagations. Although in principle any scheme could be chosen to perform the wave propagations, the elastic finite‐ difference method is used because it directly simulates the elastic wave equation and can handle complex, and thus realistic, distributions of elastic parameters. This method of inversion is costly since it is similar to an iterative prestack shot‐profile migration. However, it has greater power than any migration since it solves for the P‐wave velocity, S‐wave velocity, and density and can handle very general situations including transmission problems. Three main weaknesses of this technique are that it requires fairly accurate a priori knowledge of the low‐ wavenumber velocity model, it assumes Gaussian model statistics, and it is very computer‐intensive. All these problems seem surmountable. The low‐wavenumber information can be obtained either by a prior tomographic step, by the conventional normal‐moveout method, by a priori knowledge and empirical relationships, or by adding an additional inversion step for low wavenumbers to each iteration. The Gaussian statistics can be altered by preconditioning the gradient direction, perhaps to make the solution blocky in appearance like well logs, or by using large model variances in the inversion to reduce the effect of the Gaussian model constraints. Moreover, with some improvements to the algorithm and more parallel computers, it is hoped the technique will soon become routinely feasible.

Sign in / Sign up

Export Citation Format

Share Document