Estimations of formation velocity, permeability, and shear‐wave anisotropy using acoustic logs

Geophysics ◽  
1996 ◽  
Vol 61 (2) ◽  
pp. 437-443 ◽  
Author(s):  
Ningya Cheng ◽  
Chuen Hon Cheng

Field data sets collected by an array monopole acoustic logging tool and a shear wave logging tool are processed and interpreted. The P‐ and S‐wave velocities of the formation are determined by threshold detection with cross‐correlation correction from the full waveform and the shear‐wave log, respectively. The array monopole acoustic logging data are also processed using the extended Prony’s method to estimate the borehole Stoneley wave phase velocity and attenuation as a function of frequency. The well formation between depths of 2950 and 3150 ft (899 and 960 m) can be described as an isotropic elastic medium. The inverted [Formula: see text] from the Stoneley wave phase velocity is in excellent agreement with the shear‐wave log results in this section. The well formation between the depths of 3715 and 3780 ft (1132 and 1152 m) can be described as a porous medium with shear‐wave velocity anisotropy about 10% to 20% and with the symmetry axis perpendicular to the borehole axis. The disagreement between the shear‐wave velocity from the Stoneley wave inversion and the direct shear‐wave log velocity in this section is beyond the errors in the measurements. Estimated permeabilities from low‐frequency Stoneley wave velocity and attenuation data are in good agreement with the core measurements. Also it is proven that the formation permeability is not the cause of the discrepancy. From the estimated “shear/pseudo‐Rayleigh” phase velocities in the array monopole log and the 3-D finite‐difference synthetics in the anisotropic formation, the discrepancy can best be explained as shear‐wave anisotropy.

Geophysics ◽  
2021 ◽  
pp. 1-57
Author(s):  
Bowen Li ◽  
Alexey Stovas

Characterizing the kinematics of seismic waves in elastic vertical transversely isotropic (VTI) media involves four independent parameters. To reduce the complexity, the acoustic approximation for P-waves reduces the number of required parameters to three by setting the vertical S-wave velocity to zero. However, since only the SV-wave phase velocities parallel or perpendicular to the symmetry axis are indirectly set to zero, the acoustic approximation leads to coupled P-wave components and SV-wave artifacts. The new acoustic approximation suggests setting the vertical S-wave velocity as a phase angle-dependent variable so that the SV-wave phase velocity is zero at all phase angles. We find that manipulating this parameter is a valid way for P-wave approximation, but doing so inevitably leads to zero- or non-zero-valued spurious SV-wave components. Thus, we have developed a novel approach to efficiently approximate and thoroughly separate the two wave modes in VTI media. First, the exact P- and SV-wave phase velocity expressions are rewritten by introducing an auxiliary function. After confirming the insensitivity of this function, we construct a new expression for it and obtain simplified P- and SV-wave phase velocity expressions, which are three- and four-parameter, respectively. This approximation process leads to the same reasonable error for both wave modes. Accuracy analysis indicates that for the P-wave, the overall accuracy performance of our approach is comparable to that of some existing three-parameter approximations. We then derive the corresponding P- and SV-wave equations in tilted transversely isotropic (TTI) media and provide two available solutions, the hybrid finite-difference/pseudo-spectral scheme and the low-rank approach. Numerical examples illustrate the separability and high accuracy of the proposed P- and SV-wave simulation methods in TTI media.


2014 ◽  
Vol 88 ◽  
pp. 178-191 ◽  
Author(s):  
Hao-peng Chen ◽  
Liang-bao Zhu ◽  
Qing-dong Wang ◽  
Pan Zhang ◽  
Ying-hang Yang

2021 ◽  
Author(s):  
Prabhakar Kumar ◽  
Dibakar Ghosal

<p>The continent-continent collision between the Indian and Asian Plate formed a series of major faults from north to south along the Himalayan belt. Among these Himalayan Frontal Thrust (HFT) is the southernmost and youngest one and is tectonically very active. Any information on the shear wave velocity distribution across the fault is therefore very important. In this study, we have used the Wide Angle Multichannel Analysis of Surface Wave (WAMASW) to estimate the subsurface shear wave velocity profiles across HFT at Pawalgarh in Uttarakhand, India, using widely used stochastic global search Particle Swarm Optimization (PSO) and Grey wolf Optimization (GWO) algorithms. To gain confidence on the accuracy of the inversion results, we first generated an elastic synthetic seismic shot gather with ground rolls by using the forward modelling scheme of SOFI2D for a two-layer velocity depth model overlying a half-space. The generated gather was then processed in MATLAB to generate the experimental dispersion curve using the Phase shift method. We then extracted the fundamental mode for the gather and inverted it using the standard PSO and GWO algorithms and estimated 1D shear wave velocity profile. After getting acceptable results for the synthetic dataset, we then applied the PSO algorithm to generate the 1D S-wave velocity (Vs) profile across the Himalayan Frontal Thrust (HFT). In the study area, the Rayleigh wave phase velocity for the first shot varies from 444 to 743 m/s. We then obtained the 1D shear wave velocity profiles and a jump in Vs is observed across the HFT indicating variation in the sediment stiffness across the fault.</p><p><strong>Keywords: </strong>WAMASW, dispersion, Meta- Heuristic, PSO, GWO, 1D Shear wave velocity</p><p> </p>


Sign in / Sign up

Export Citation Format

Share Document