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

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.

Geophysics ◽  
1988 ◽  
Vol 53 (9) ◽  
pp. 1207-1218 ◽  
Author(s):  
L. J. Baker ◽  
G. A. Winbow

Wave trains produced by conventional and multipole sonic logging tools may be expected to depend upon whether the formation has been altered by the drilling process. Such alteration may include invasion by drilling fluids and/or drilling damage. We show, by theoretical modeling, that a dipole or quadrupole tool operating at conventional logging frequencies (in the range 10–20 kHz) detects P waves from the virgin formation with a much higher signal‐to‐noise ratio than does a monopole tool. This permits the multipole tool to measure formation P‐wave velocities two to three times farther away from the borehole than a conventional monopole tool. This larger radius of investigation typically extends beyond the altered zone for most situations, even for sources and receivers spaced several meters apart. Our conclusion is valid only if the velocity of the altered zone is less than the velocity of the virgin formation. If the velocity of the altered zone is larger than that of the virgin formation, no appreciable benefit is obtained by using a multipole tool. Similar results are demonstrated in measuring S‐wave velocities.


Geophysics ◽  
1982 ◽  
Vol 47 (11) ◽  
pp. 1512-1520 ◽  
Author(s):  
S. T. Chen

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.


Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. D369-D381 ◽  
Author(s):  
Elliot J. H. Dahl ◽  
Kyle T. Spikes

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.


Geophysics ◽  
1980 ◽  
Vol 45 (1) ◽  
pp. 3-17 ◽  
Author(s):  
Franklyn K. Levin

P‐wave, SV‐wave, and SH‐wave velocities are computed for transversely isotropic solids formed from two isotropic solids. The combinations are shale‐sandstone and shale‐limestone solids of an earlier paper (Levin, 1979), but one velocity of the nonshale component is allowed to vary over the range of Poisson’s ratios σ = 0 to σ = 0.45, i.e., from a rigid solid to a near‐liquid. When the S‐wave velocity of either the sandstone or limestone is varied, the ratio of horizontal P‐wave velocity to vertical P‐wave velocity goes through a maximum as σ increases and subsequently falls to values less than unity as σ approaches 0.5. The P‐wave velocity that would be found with a short surface spread also goes through a maximum and, at σ = 0.5, is less than the P‐wave velocity of either isotropic component. SV‐wave velocities found for data from a short spread are unreasonably large; SH‐wave velocities decrease monotonically as σ increases, but the ratio of horizontal SH‐wave velocity to vertical SH‐wave velocity goes through a minimum of unity.


1996 ◽  
Vol 86 (6) ◽  
pp. 1704-1713 ◽  
Author(s):  
R. D. Catchings ◽  
W. H. K. Lee

Abstract The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (<70 m) P-and S-wave velocity structure. Our seismic velocity models for the Tarzana site indicate that the local velocity structure may have contributed significantly to the observed shaking. P-wave velocities range from 0.9 to 1.65 km/sec, and S-wave velocities range from 0.20 and 0.6 km/sec for the upper 70 m. We also found evidence for a local S-wave low-velocity zone (LVZ) beneath the top of the hill. The LVZ underlies a CDMG strong-motion recording site at depths between 25 and 60 m below ground surface (BGS). Our velocity model is consistent with the near-surface (<30 m) P- and S-wave velocities and Poisson's ratios measured in a nearby (<30 m) borehole. High Poisson's ratios (0.477 to 0.494) and S-wave attenuation within the LVZ suggest that the LVZ may be composed of highly saturated shales of the Modelo Formation. Because the lateral dimensions of the LVZ approximately correspond to the areas of strongest shaking, we suggest that the highly saturated zone may have contributed to localized strong shaking. Rock sites are generally considered to be ideal locations for site response in urban areas; however, localized, highly saturated rock sites may be a hazard in urban areas that requires further investigation.


Geophysics ◽  
2021 ◽  
pp. 1-78
Author(s):  
Da Shuai ◽  
Alexey Stovas ◽  
Jianxin Wei ◽  
Bangrang Di ◽  
Yang Zhao

The linear slip theory is gradually being used to characterize seismic anisotropy. If the transversely isotropic medium embeds vertical fractures (VFTI medium), the effective medium becomes orthorhombic. The vertical fractures, in reality, may exist in any azimuth angle which leads the effective medium to be monoclinic. We apply the linear slip theory to create a monoclinic medium by only introducing three more physical meaning parameters: the fracture preferred azimuth angle, the fracture azimuth angle, and the angular standard deviation. First, we summarize the effective compliance of a rock as the sum of the background matrix compliance and the fracture excess compliance. Then, we apply the Bond transformation to rotate the fractures to be azimuth dependent, introduce a Gaussian function to describe the fractures' azimuth distribution assuming that the fractures are statistically distributed around the preferred azimuth angle, and average each fracture excess compliance over azimuth. The numerical examples investigate the influence of the fracture azimuth distribution domain and angular standard deviation on the effective stiffness coefficients, elastic wave velocities, and anisotropy parameters. Our results show that the fracture cluster parameters have a significant influence on the elastic wave velocities. The fracture azimuth distribution domain and angular standard deviation have a bigger influence on the orthorhombic anisotropy parameters in the ( x2, x3) plane than that in the ( x1, x3) plane. The fracture azimuth distribution domain and angular standard deviation have little influence on the monoclinic anisotropy parameters responsible for the P-wave NMO ellipse and have a significant influence on the monoclinic anisotropy parameters responsible for the S1- and S2-wave NMO ellipse. The effective monoclinic can be degenerated into the VFTI medium.


2019 ◽  
Vol 283 ◽  
pp. 06003
Author(s):  
Guangxue Zheng ◽  
Hanhao Zhu ◽  
Jun Zhu

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.


2002 ◽  
Vol 42 (1) ◽  
pp. 627
Author(s):  
R.G. Williams ◽  
G. Roberts ◽  
K. Hawkins

Seismic energy that has been mode converted from pwave to s-wave in the sub-surface may be recorded by multi-component surveys to obtain information about the elastic properties of the earth. Since the energy converted to s-wave is missing from the p-wave an alternative to recording OBC multi-component data is to examine p-wave data for the missing energy. Since pwave velocities are generally faster than s-wave velocities, then for a given reflection point the converted s-wave signal reaches the surface at a shorter offset than the equivalent p-wave information. Thus, it is necessary to record longer offsets for p-wave data than for multicomponent data in order to measure the same information.A non-linear, wide-angle (including post critical) AVO inversion has been developed that allows relative changes in p-wave velocities, s-wave velocities and density to be extracted from long offset p-wave data. To extract amplitudes at long offsets for this inversion it is necessary to image the data correctly, including correcting for higher order moveout and possibly anisotropy if it is present.The higher order moveout may itself be inverted to yield additional information about the anisotropy of the sub-surface.


2015 ◽  
Vol 76 (2) ◽  
Author(s):  
Badee Alshameri ◽  
Ismail Bakar ◽  
Aziman Madun ◽  
Edy Tonnizam Mohamad

One of the main geophysical tools (seismic tools) in the laboratory is the bender element. This tool can be used to measure some dynamic soil properties (e.g. shear and Young’s modulus). However, even if it relatively simple to use the bender element, inconsistent testing procedures can cause poor quality in the bender element data. One of the bender element procedure that always neglected is the alignment (different positions of bender element receiver to the transmitter in the vertical axis). The alignment effect was evaluated via changing the horizontal distance between transmitter and receiver starting from 0 to 110 mm for two sizes of the sample's thickness (i.e. 63.17 mm and 91.51 mm). Five methods were applied to calculate the travel times. Those methods were as the following: visually, first-peak, maximum-peak, CCexcel and CCGDS. In general, the experiments indicated uncertain results for both of the P-wave (primary wave) and S-wave (secondary wave) velocities at zone of Dr:D above 0.5:1 (where Dr is the horizontal distance of the receiver from the vertical axis and D is the thickness of the sample). On the other hand, both the visual and first-peak methods show the wave velocities results are higher than obtained from other methods. However, the ratio between the amplitude of transmitter signals to receiver amplitude signal was taken to calculate the damping-slope of the P-wave and S-wave. Thus the results from damping slope show steeply slope when the ratio of  Dr:D is above 0.5:1 compare with gentle slope below ratio 0.5:1 at the sample with thickness equal to 91.51 mm, while there is no variation at a slope in sample with thickness equal to 63.17 mm.


Sign in / Sign up

Export Citation Format

Share Document