Determination of formation shear attenuation from dipole sonic log data

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. D73-D79 ◽  
Author(s):  
Qiaomu Qi ◽  
Arthur C. H. Cheng ◽  
Yunyue Elita Li
Keyword(s):  
Reservoir Characterization ◽  
Wave Attenuation ◽  
Synthetic Data ◽  
Low Frequency ◽  
P Wave ◽  
Flexural Wave ◽  
S Wave ◽  
Additional Energy ◽  
Log Data ◽  
Sonic Log

ABSTRACT Formation S-wave attenuation, when combined with compressional attenuation, serves as a potential hydrocarbon indicator for seismic reservoir characterization. Sonic flexural wave measurements provide a direct means for obtaining the in situ S-wave attenuation at log scale. The key characteristic of the flexural wave is that it propagates at the formation shear slowness and experiences shear attenuation at low frequency. However, in a fast formation, the dipole log consists of refracted P- and S-waves in addition to the flexural wave. The refracted P-wave arrives early and can be removed from the dipole waveforms through time windowing. However, the refracted S-wave, which is often embedded in the flexural wave packet, is difficult to separate from the dipole waveforms. The additional energy loss associated with the refracted S-wave results in the estimated dipole attenuation being higher than the shear attenuation at low frequency. To address this issue, we have developed a new method for accurately determining the formation shear attenuation from the dipole sonic log data. The method uses a multifrequency inversion of the frequency-dependent flexural wave attenuation based on energy partitioning. We first developed our method using synthetic data. Application to field data results in a shear attenuation log that is consistent with lithologic interpretation of other available logs.

Extraction of reflected events from sonic-log waveforms using the Karhunen-Loève transform

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. D265-D277 ◽  
Author(s):  
Junxiao Li ◽  
Kristopher A. Innanen ◽  
Guo Tao
Keyword(s):  
Principal Components ◽  
Laboratory Data ◽  
Synthetic Data ◽  
P Wave ◽  
Stoneley Wave ◽  
S Wave ◽  
Data Set ◽  
Sonic Log ◽  
Karhunen Loeve Transform ◽  
Wave Modes

Sonic-reflection logging, a recently developed borehole geophysical scheme, is in principle capable of providing a clear view of outside the well bore. In this type of acoustic well logging, a key technical obstacle is that the reflected wave signal is almost entirely obscured by the directly arriving P-, S-, and Stoneley wave modes. Effective extraction of these reflection signals from the full acoustic waveforms is therefore a critical data-processing step. We have examined the use of the Karhunen-Loève (KL) transform, combined with a band-limiting filter, as a technique for the extraction of reflections of interest from a mixture with directly arriving wave modes of much higher amplitude. Under the assumption that large energy (squared-amplitude) differences exist between each wave component, the direct Stoneley wave, S-wave, and the P-wave are eliminated sequentially by subtracting the most significant principal components, after which the remaining signal is seen to be dominated by reflected events. Thereafter, the extracted reflections can be used in migration to provide interpretable images of the structures outside the borehole. Synthetic data are used to develop and justify our procedure for subtraction of appropriate KL principal components. Laboratory data are used to demonstrate in detail the suppression of unwanted modes. For comparison, the multiscale slowness-time-coherence method is applied to extract reflections from the same data set. The procedure is exemplified on a field data case with attention paid in particular to the consequences to imaging of near-borehole structures.

Estimation of shear-wave interval attenuation from mode-converted data

Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. D11-D19 ◽  
Author(s):  
Bharath Shekar ◽  
Ilya Tsvankin
Keyword(s):  
Shear Wave ◽  
Attenuation Coefficient ◽  
Reservoir Characterization ◽  
Wave Attenuation ◽  
Pure Shear ◽  
Synthetic Data ◽  
Spectral Ratio ◽  
Ratio Method ◽  
S Wave ◽  
Target Layer

Interval attenuation measurements provide valuable information for reservoir characterization and lithology discrimination. We extend the attenuation layer-stripping method of Behura and Tsvankin to mode-converted (PS) waves with the goal of estimating the S-wave interval attenuation coefficient. By identifying PP and PS events with shared ray segments and applying the [Formula: see text] method, we first perform kinematic construction of pure shear (SS) events in the target layer and overburden. Then, the modified spectral-ratio method is used to compute the effective shear-wave attenuation coefficient for the target reflection. Finally, application of the dynamic version of velocity-independent layer stripping to the constructed SS reflections yields the interval S-wave attenuation coefficient in the target layer. The attenuation coefficient estimated for a range of source-receiver offsets can be inverted for the interval attenuation parameters. The method is tested on multicomponent synthetic data generated with the anisotropic reflectivity method for layered VTI (transversely isotropic with a vertical symmetry axis) and orthorhombic media.

P- and S-wave attenuation logs from monopole sonic data

Geophysics ◽  
2000 ◽  
Vol 65 (3) ◽  
pp. 755-765 ◽  
Author(s):  
Xinhua Sun ◽  
Xiaoming Tang ◽  
C. H. (Arthur) Cheng ◽  
L. Neil Frazer
Keyword(s):  
Wave Attenuation ◽  
Fluid Velocity ◽  
P Wave ◽  
Reservoir Rock ◽  
Rock Properties ◽  
S Wave ◽  
Intrinsic Attenuation ◽  
Modified Method ◽  
Receiver Array ◽  
Attenuation Profiles

In this paper, a modification of an existing method for estimating relative P-wave attenuation is proposed. By generating synthetic waveforms without attenuation, the variation of geometrical spreading related to changes in formation properties with depth can be accounted for. With the modified method, reliable P- and S-wave attenuation logs can be extracted from monopole array acoustic waveform log data. Synthetic tests show that the P- and S-wave attenuation values estimated from synthetic waveforms agree well with their respective model values. In‐situ P- and S-wave attenuation profiles provide valuable information about reservoir rock properties. Field data processing results show that this method gives robust estimates of intrinsic attenuation. The attenuation profiles calculated independently from each waveform of an eight‐receiver array are consistent with one another. In fast formations where S-wave velocity exceeds the borehole fluid velocity, both P-wave attenuation ([Formula: see text]) and S-wave attenuation ([Formula: see text]) profiles can be obtained. P- and S-wave attenuation profiles and their comparisons are presented for three reservoirs. Their correlations with formation lithology, permeability, and fractures are also presented.

Anisotropic correction of sonic logs in wells with large relative dip

Geophysics ◽  
2008 ◽  
Vol 73 (1) ◽  
pp. E1-E5 ◽  
Author(s):  
Lev Vernik
Keyword(s):  
Reservoir Characterization ◽  
Low Frequency ◽  
P Wave ◽  
Amplitude Variation With Offset ◽  
Seismic Amplitude ◽  
Avo Inversion ◽  
Pore Pressure Prediction ◽  
Siliciclastic Rocks ◽  
Dip Angle ◽  
Sonic Logs

Seismic reservoir characterization and pore-pressure prediction projects rely heavily on the accuracy and consistency of sonic logs. Sonic data acquisition in wells with large relative dip is known to suffer from anisotropic effects related to microanisotropy of shales and thin-bed laminations of sand, silt, and shale. Nonetheless, if anisotropy parameters can be related to shale content [Formula: see text] in siliciclastic rocks, then I show that it is straightforward to compute the anisotropy correction to both compressional and shear logs using [Formula: see text] and the formation relative dip angle. The resulting rotated P-wave sonic logs can be used to enhance time-depth ties, velocity to effective stress transforms, and low-frequency models necessary for prestack seismic amplitude variation with offset (AVO) inversion.

Direct shear-wave seismic survey in Sanhu area, Qaidam Basin, west China

2022 ◽  
Vol 41 (1) ◽  
pp. 47-53
Author(s):  
Zhiwen Deng ◽  
Rui Zhang ◽  
Liang Gou ◽  
Shaohua Zhang ◽  
Yuanyuan Yue ◽  
...  
Keyword(s):  
Shear Wave ◽  
Reservoir Characterization ◽  
Qaidam Basin ◽  
P Wave ◽  
Data Sets ◽  
Direct Shear ◽  
Subsurface Structure ◽  
S Wave ◽  
West China ◽  
Wave Data

The formation containing shallow gas clouds poses a major challenge for conventional P-wave seismic surveys in the Sanhu area, Qaidam Basin, west China, as it dramatically attenuates seismic P-waves, resulting in high uncertainty in the subsurface structure and complexity in reservoir characterization. To address this issue, we proposed a workflow of direct shear-wave seismic (S-S) surveys. This is because the shear wave is not significantly affected by the pore fluid. Our workflow includes acquisition, processing, and interpretation in calibration with conventional P-wave seismic data to obtain improved subsurface structure images and reservoir characterization. To procure a good S-wave seismic image, several key techniques were applied: (1) a newly developed S-wave vibrator, one of the most powerful such vibrators in the world, was used to send a strong S-wave into the subsurface; (2) the acquired 9C S-S data sets initially were rotated into SH-SH and SV-SV components and subsequently were rotated into fast and slow S-wave components; and (3) a surface-wave inversion technique was applied to obtain the near-surface shear-wave velocity, used for static correction. As expected, the S-wave data were not affected by the gas clouds. This allowed us to map the subsurface structures with stronger confidence than with the P-wave data. Such S-wave data materialize into similar frequency spectra as P-wave data with a better signal-to-noise ratio. Seismic attributes were also applied to the S-wave data sets. This resulted in clearly visible geologic features that were invisible in the P-wave data.

Retrieval of source-extent parameters and the interpretation of corner frequency

1983 ◽  
Vol 73 (6A) ◽  
pp. 1499-1511
Author(s):  
Paul Silver
Keyword(s):  
Body Wave ◽  
Amplitude Spectrum ◽  
Low Frequency ◽  
P Wave ◽  
Corner Frequency ◽  
S Wave ◽  
Dislocation Source ◽  
Statistical Measures ◽  
Dislocation Models ◽  
Planar Dislocation

Abstract A method is proposed for retrieving source-extent parameters from far-field body-wave data. At low frequency, the normalized P- or S-wave displacement amplitude spectrum can be approximated by |Ω^(r^,ω)| = 1 − τ2(r^)ω2/2 where r^ specifies a point on the focal sphere. For planar dislocation sources, τ2(r^) is linearly related to statistical measures of source dimension, source duration, and directivity. τ2(r^) can be measured as the curvature of |Ω^(r^,ω)| at ω = 0 or the variance of the pulse Ω^(r^,t). The quantity ωc=2τ−1(r^) is contrasted with the traditional corner frequency ω0, defined as the frequency at the intersection of the low- and high-frequency trends of |Ω^(r^,ω)|. For dislocation models without directivity, ωc(P) ≧ ωc(S) for any r^. A mean corner frequency defined by averaging τ2(r^) over the focal sphere, ω¯c=2<τ2(r^)>−1/2, satisfies ωc(P) > ωc(S) for any dislocation source. This behavior is not shared by ω0. It is shown that ω0 is most sensitive to critical times in the rupture history of the source, whereas ωc is determined by the basic parameters of source extent. Evidence is presented that ωc is the corner frequency measured on actual seismograms. Thus, the commonly observed corner frequency shift (P-wave corner greater than the S-wave corner), now viewed as a shift in ωc is simply a result of spatial finiteness and is expected to be a property of any dislocation source. As a result, the shift cannot be used as a criterion for rejecting particular dislocation models.

Local and global fluid-flow effects on dipole flexural waves

Geophysics ◽  
2020 ◽  
Vol 85 (1) ◽  
pp. D1-D11
Author(s):  
Elliot J. H. Dahl ◽  
Kyle T. Spikes
Keyword(s):  
Fluid Flow ◽  
Low Frequency ◽  
Summation Method ◽  
Sensitivity Analyses ◽  
Flexural Wave ◽  
Biot Theory ◽  
S Wave ◽  
Local Flow ◽  
Flow Effects ◽  
Wave Induced

Wave-induced fluid flow (WIFF) can significantly alter the effective formation velocities and cause increasing waveform dispersion and attenuation. We have used modified frame moduli from the theory of Chapman together with the classic Biot theory to improve the understanding of local- and global-flow effects on dipole flexural wave modes in boreholes. We investigate slow and fast formations with and without compliant pores, which induce local flow. The discrete wavenumber summation method generates the waveforms, which are then processed with the weighted spectral semblance method to compare with the solution of the period equation. We find compliant pores to decrease the resulting effective formation P- and S-wave velocities, that in turn decrease the low-frequency velocity limit of the flexural wave. Furthermore, depending on the frequency at which the local-flow dispersion occurs, different S-wave velocity predictions from the flexural wave become possible. This issue is investigated through changing the local-flow critical frequency. Sensitivity analyses of the flexural-wave phase velocity to small changes in WIFF parameters indicate the modeling to be mostly sensitive to compliant pores in slow and fast formations.

scholarly journals A method of geo-acoustic parameter inversion in shallow sea by the Bayesian theory and the acoustic pressure field

2019 ◽  
Vol 283 ◽  
pp. 06003
Author(s):  
Guangxue Zheng ◽  
Hanhao Zhu ◽  
Jun Zhu
Keyword(s):  
Sound Pressure ◽  
Pressure Field ◽  
Wave Attenuation ◽  
Acoustic Parameter ◽  
Bayesian Theory ◽  
P Wave ◽  
Shallow Sea ◽  
S Wave ◽  
Wave Velocities ◽  
Parameter Inversion

A method of geo-acoustic parameter inversion based on the Bayesian theory is proposed for the acquisition of acoustic parameters in shallow sea with the elastic seabed. Firstly, the theoretical prediction value of the sound pressure field is calculated by the fast field method (FFM). According to the Bayesian theory, we establish the misfit function between the measured sound pressure field and the theoretical pressure field. It is under the assumption of Gaussian data errors which are in line with the likelihood function. Finally, the posterior probability density (PPD) of parameters is given as the result of inversion. Our research is conducted in the light of Metropolis sample rules. Apart from numerical simulations, a scaled model experiment has been taken in the laboratory tank. The results of numerical simulations and tank experiments show that sound pressure field calculated by the result of inversion is consistent with the measured sound pressure field. Besides, s-wave velocities, p-wave velocities and seafloor density have fewer uncertainties and are more sensitive to complex sound pressure than s-wave attenuation and p-wave attenuation. The received signals calculated by inversion results are keeping with received signals in the experiment which verify the effectiveness of this method.

Near-surface scattering phenomena and implications for tunnel detection

2015 ◽  
Vol 3 (1) ◽  
pp. SF43-SF54 ◽  
Author(s):  
Shelby L. Peterie ◽  
Richard D. Miller
Keyword(s):  
Synthetic Data ◽  
P Wave ◽  
S Wave ◽  
Seismic Line ◽  
Data Set ◽  
Near Surface ◽  
Tunnel Detection ◽  
Impedance Contrast ◽  
Diffraction Imaging ◽  
Imaging Algorithms

Tunnel locations are accurately interpreted from diffraction sections of focused mode converted P- to S-wave diffractions from a perpendicular tunnel and P-wave diffractions from a nonperpendicular (oblique) tunnel. Near-surface tunnels are ideal candidates for diffraction imaging due to their small size relative to the seismic wavelength and large acoustic impedance contrast at the tunnel interface. Diffraction imaging algorithms generally assume that the velocities of the primary wave and the diffracted wave are approximately equal, and that the diffraction apex is recorded directly above the scatterpoint. Scattering phenomena from shallow tunnels with kinematic properties that violate these assumptions were observed in one field data set and one synthetic data set. We developed the traveltime equations for mode-converted and oblique diffractions and demonstrated a diffraction imaging algorithm designed for the roll-along style of acquisition. Potential processing and interpretation pitfalls specific to these diffraction types were identified. Based on our observations, recommendations were made to recognize and image mode-converted and oblique diffractions and accurately interpret tunnel depth, horizontal location, and azimuth with respect to the seismic line.

Rock-physics modeling of ultrasonic P- and S-wave attenuation in partially frozen brine and unconsolidated sand systems and comparison with laboratory measurements

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. MR153-MR171 ◽  
Author(s):  
Linsen Zhan ◽  
Jun Matsushima
Keyword(s):  
Ultrasonic Wave ◽  
Wave Attenuation ◽  
Freezing Point ◽  
Rock Physics ◽  
Ultrasonic Waves ◽  
P Wave ◽  
Dominant Factor ◽  
S Wave ◽  
High Attenuation ◽  
Unconsolidated Sand

The nonintuitive observation of the simultaneous high velocity and high attenuation of ultrasonic waves near the freezing point of brine was previously measured in partially frozen systems. However, previous studies could not fully elucidate the attenuation variation of ultrasonic wave propagation in a partially frozen system. We have investigated the potential attenuation mechanisms responsible for previously obtained laboratory results by modeling ultrasonic wave transmission in two different partially frozen systems: partially frozen brine (two phases composed of ice and unfrozen brine) and unconsolidated sand (three phases composed of ice, unfrozen brine, and sand). We adopted two different rock-physics models: an effective medium model for partially frozen brine and a three-phase extension of the Biot model for partially frozen unconsolidated sand. For partially frozen brine, our rock-physics study indicated that squirt flow caused by unfrozen brine inclusions in porous ice could be responsible for high P-wave attenuation around the freezing point. Decreasing P-wave attenuation below the freezing point can be explained by the gradual decrease of squirt flow due to the gradual depletion of unfrozen brine. For partially frozen unconsolidated sand, our rock-physics study implied that squirt flow between ice grains is a dominant factor for P-wave attenuation around the freezing point. With decreasing temperature lower than the freezing point, the friction between ice and sand grains becomes more dominant for P-wave attenuation because the decreasing amount of unfrozen brine reduces squirt flow between ice grains, whereas the generation of ice increases the friction. The increasing friction between ice and sand grains caused by ice formation is possibly responsible for increasing the S-wave attenuation at decreasing temperatures. Then, further generation of ice with further cooling reduces the elastic contrast between ice and sand grains, hindering their relative motion; thus, reducing the P- and S-wave attenuation.

Sign in / Sign up

Export Citation Format

Share Document