scholarly journals Phased array compaction cell for measurement of the transversely isotropic elastic properties of compacting sediments

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. WA113-WA123 ◽  
Author(s):  
Kurt T. Nihei ◽  
Seiji Nakagawa ◽  
Frederic Reverdy ◽  
Larry R. Myer ◽  
Luca Duranti ◽  
...  
Keyword(s):  
Elastic Properties ◽  
Phased Array ◽  
Transversely Isotropic ◽  
P Wave ◽  
S Wave ◽  
P Waves ◽  
Array Measurements ◽  
Experimental Apparatus ◽  
Marine Sediment Sample ◽  
Plane P Waves

Sediments undergoing compaction typically exhibit transversely isotropic (TI) elastic properties. We present a new experimental apparatus, the phased array compaction cell, for measuring the TI elastic properties of clay-rich sediments during compaction. This apparatus uses matched sets of P- and S-wave ultrasonic transducers located along the sides of the sample and an ultrasonic P-wave phased array source, together with a miniature P-wave receiver on the top and bottom ends of the sample. The phased array measurements are used to form plane P-waves that provide estimates of the phase velocities over a range of angles. From these measurements, the five TI elastic constants can be recovered as the sediment is compacted, without the need for sample unloading, recoring, or reorienting. This paper provides descriptions of the apparatus, the data processing, and an application demonstrating recovery of the evolving TI properties of a compacting marine sediment sample.

SCATTERING OF COMPRESSION WAVES BY SPHERICAL OBSTACLES

Geophysics ◽  
1959 ◽  
Vol 24 (1) ◽  
pp. 30-39 ◽  
Author(s):  
Leon Knopoff
Keyword(s):  
Wave Beam ◽  
Wave Length ◽  
Incident Beam ◽  
P Wave ◽  
Rigid Sphere ◽  
S Wave ◽  
P Waves ◽  
Compression Waves ◽  
Plane P Waves ◽  
Special Case

The scattering of plane P waves by a spherical obstacle is formulated. A calculation is given for the special case of scattering by a perfectly rigid sphere in which the medium outside has a Poisson’s ratio of [Formula: see text]. The range of sizes of obstacles used in the calculation includes radii very small compared with wave length and radii comparable to the wave length. For incident P waves, scattered P and S are computed with shifts in time phase occurring in both with respect to the incident beam. For small obstacles, the scattered S wave is generally broadside to the scattered P‐wave beam.

Seismic vertical transversely isotropic parameter inversion from P- and S-wave cross-borehole measurements in an aquifer environment

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. D101-D116
Author(s):  
Julius K. von Ketelhodt ◽  
Musa S. D. Manzi ◽  
Raymond J. Durrheim ◽  
Thomas Fechner
Keyword(s):  
Transversely Isotropic ◽  
P Wave ◽  
Perturbation Equation ◽  
S Wave ◽  
P Waves ◽  
Data Set ◽  
Near Surface ◽  
S Waves ◽  
Wave Measurements ◽  
Wave Splitting

Joint P- and S-wave measurements for tomographic cross-borehole analysis can offer more reliable interpretational insight concerning lithologic and geotechnical parameter variations compared with P-wave measurements on their own. However, anisotropy can have a large influence on S-wave measurements, with the S-wave splitting into two modes. We have developed an inversion for parameters of transversely isotropic with a vertical symmetry axis (VTI) media. Our inversion is based on the traveltime perturbation equation, using cross-gradient constraints to ensure structural similarity for the resulting VTI parameters. We first determine the inversion on a synthetic data set consisting of P-waves and vertically and horizontally polarized S-waves. Subsequently, we evaluate inversion results for a data set comprising jointly measured P-waves and vertically and horizontally polarized S-waves that were acquired in a near-surface ([Formula: see text]) aquifer environment (the Safira research site, Germany). The inverted models indicate that the anisotropy parameters [Formula: see text] and [Formula: see text] are close to zero, with no P-wave anisotropy present. A high [Formula: see text] ratio of up to nine causes considerable SV-wave anisotropy despite the low magnitudes for [Formula: see text] and [Formula: see text]. The SH-wave anisotropy parameter [Formula: see text] is estimated to be between 0.05 and 0.15 in the clay and lignite seams. The S-wave splitting is confirmed by polarization analysis prior to the inversion. The results suggest that S-wave anisotropy may be more severe than P-wave anisotropy in near-surface environments and should be taken into account when interpreting cross-borehole S-wave data.

The joint nonhyperbolic moveout inversion of PP and PS data in VTI media

Geophysics ◽  
2002 ◽  
Vol 67 (6) ◽  
pp. 1929-1932 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Vertical Velocity ◽  
Single Layer ◽  
Transversely Isotropic ◽  
Accurate Estimate ◽  
P Wave ◽  
S Wave ◽  
P Waves ◽  
Layer Cake ◽  
Two Parameters ◽  
Vti Media

Nonhyperbolic moveout of P‐waves in horizontally layered transversely isotropic media with a vertical symmetry axis (VTI) can be used to estimate the anellipticity coefficient η in addition to the NMO velocity Vnmo,P. Those two parameters are sufficient for time processing of P‐wave data (despite a certain instability in the inversion for η), but they do not constrain the vertical velocity VP0 and the depth scale of the model. It has been suggested in the literature that this ambiguity in the depth‐domain velocity analysis for layer‐cake VTI media can be resolved by combining long‐spread reflection traveltimes of P‐waves and mode‐converted PSV‐waves. Here, we show that reflection traveltimes of horizontal PSV events help to determine the ratio of the P‐ and S‐wave vertical velocities and the NMO velocity of SV‐waves, and they give a more accurate estimate of η. However, nonhyperbolic moveout of PSV‐waves turns out to be mostly controlled by wide‐angle P‐wave traveltimes and does not provide independent information for the inversion. As a result, even for a single‐layer model and uncommonly large offsets, traveltimes of P‐ and PSV‐waves cannot be inverted for the vertical velocity and anisotropic parameters ε and δ. To reconstruct the horizontally layered VTI model from surface data, it is necessary to combine long‐spread traveltimes of pure P and SV reflections.

A finite-difference approach for solving pure quasi-P-wave equations in transversely isotropic and orthorhombic media

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. C161-C172 ◽  
Author(s):  
Xueyan Li ◽  
Hejun Zhu
Keyword(s):  
Finite Difference ◽  
Wave Equations ◽  
Differential Operators ◽  
Transversely Isotropic ◽  
Coupled Systems ◽  
P Wave ◽  
Finite Difference Schemes ◽  
S Wave ◽  
P Waves ◽  
Pseudo Differential Operators

Starting from the dispersion relation and setting S-wave velocity along symmetry axes to zero, pseudoacoustic-wave equations have been proposed to describe the kinematics of acoustic wavefields in transversely isotropic (TI) and orthorhombic media. To date, the numerical stability of the pseudoacoustic-wave equations has been improved by developing coupled systems of wave equations; however, most simulations still suffer from S-wave artifacts that are the fundamental solutions of the fourth- and sixth-order partial differential equations. Pure quasi-P-wave equations accurately describe the traveltimes of P-waves in TI and orthorhombic media and are free of S-wave artifacts. However, it is difficult to directly solve the pure quasi-P-wave equations using conventional finite-difference schemes due to the presence of pseudo-differential operators. We approximated these pseudo-differential operators by algebraic expressions, whose coefficients can be determined by minimizing differences between the true and approximated values of the pseudo-differential operators in the wavenumber domain. The derived new coupled systems involve modified acoustic-wave equations and a Poisson’s equation that can be solved by conventional finite-difference stencils and fast Poisson’s solver. Several 2D and 3D numerical examples demonstrate that the simulations based on the new systems are free of S-wave artifacts and have correct kinematics of quasi-P-waves in TI and orthorhombic media.

scholarly journals New acoustic approximation for transversely isotropic media with a vertical symmetry axis

Geophysics ◽  
2019 ◽  
Vol 85 (1) ◽  
pp. C1-C12 ◽  
Author(s):  
Shibo Xu ◽  
Alexey Stovas ◽  
Tariq Alkhalifah ◽  
Hitoshi Mikada
Keyword(s):  
Wave Propagation ◽  
Seismic Data ◽  
Eikonal Equation ◽  
Symmetry Axis ◽  
Transversely Isotropic ◽  
P Wave ◽  
Seismic Data Processing ◽  
S Wave ◽  
P Waves ◽  
Acoustic Approximation

Seismic data processing in the elastic anisotropic model is complicated due to multiparameter dependency. Approximations to the P-wave kinematics are necessary for practical purposes. The acoustic approximation for P-waves in a transversely isotropic medium with a vertical symmetry axis (VTI) simplifies the description of wave propagation in elastic media, and as a result, it is widely adopted in seismic data processing and analysis. However, finite-difference implementations of that approximation are plagued with S-wave artifacts. Specifically, the resulting wavefield also includes artificial diamond-shaped S-waves resulting in a redundant signal for many applications that require pure P-wave data. To derive a totally S-wave-free acoustic approximation, we have developed a new acoustic approximation for pure P-waves that is totally free of S-wave artifacts in the homogeneous VTI model. To keep the S-wave velocity equal to zero, we formulate the vertical S-wave velocity to be a function of the model parameters, rather than setting it to zero. Then, the corresponding P-wave phase and group velocities for the new acoustic approximation are derived. For this new acoustic approximation, the kinematics is described by a new eikonal equation for pure P-wave propagation, which defines the new vertical slowness for the P-waves. The corresponding perturbation-based approximation for our new eikonal equation is used to compare the new equation with the original acoustic eikonal. The accuracy of our new P-wave acoustic approximation is tested on numerical examples for homogeneous and multilayered VTI models. We find that the accuracy of our new acoustic approximation is as good as the original one for the phase velocity, group velocity, and the kinematic parameters such as vertical slowness, traveltime, and relative geometric spreading. Therefore, the S-wave-free acoustic approximation could be further applied in seismic processing that requires pure P-wave data.

scholarly journals An efficient Helmholtz solver for acoustic transversely isotropic media

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. C75-C83 ◽  
Author(s):  
Zedong Wu ◽  
Tariq Alkhalifah
Keyword(s):  
Wave Equation ◽  
Transversely Isotropic ◽  
Anisotropic Media ◽  
P Wave ◽  
Seismic Exploration ◽  
S Wave ◽  
P Waves ◽  
Transversely Isotropic Media ◽  
Isotropic Media ◽  
Separable Approximation

The acoustic approximation, even for anisotropic media, is widely used in current industry imaging and inversion algorithms mainly because P-waves constitute most of the energy recorded in seismic exploration. The resulting acoustic formulas tend to be simpler, resulting in more efficient implementations, and they depend on fewer medium parameters. However, conventional solutions of the acoustic-wave equation with higher-order derivatives suffer from S-wave artifacts. Thus, we separate the quasi-P-wave propagation in anisotropic media into the elliptic anisotropic operator (free of the artifacts) and the nonelliptic anisotropic components, which form a pseudodifferential operator. We then develop a separable approximation of the dispersion relation of nonelliptic-anisotropic components, specifically for transversely isotropic media. Finally, we iteratively solve the simpler lower-order elliptical wave equation for a modified source function that includes the nonelliptical terms represented in the Fourier domain. A frequency-domain Helmholtz formulation of the approach renders the iterative implementation efficient because the cost is dominated by the lower-upper decomposition of the impedance matrix for the simpler elliptical anisotropic model. In addition, the resulting wavefield is free of S-wave artifacts and has a balanced amplitude. Numerical examples indicate that the method is reasonably accurate and efficient.

Geometrical spreading in a layered transversely isotropic medium with vertical symmetry axis

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 2082-2091 ◽  
Author(s):  
Bjørn Ursin ◽  
Ketil Hokstad
Keyword(s):  
Ray Tracing ◽  
Seismic Data ◽  
Symmetry Axis ◽  
Transversely Isotropic ◽  
Analytical Approximation ◽  
Geometrical Spreading ◽  
S Wave ◽  
Weak Anisotropy ◽  
P Waves ◽  
Amplitude Versus

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.

Reflection moveout approximations for P-waves in a moderately anisotropic homogeneous tilted transverse isotropy layer

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. C175-C185 ◽  
Author(s):  
Ivan Pšenčík ◽  
Véronique Farra
Keyword(s):  
Transverse Isotropy ◽  
Transversely Isotropic ◽  
P Wave ◽  
Weak Anisotropy ◽  
P Waves ◽  
3D Space ◽  
Axis Of Symmetry ◽  
Relative Errors ◽  
Symmetry Tests ◽  
Ti Media

We have developed approximate nonhyperbolic P-wave moveout formulas applicable to weakly or moderately anisotropic media of arbitrary anisotropy symmetry and orientation. Instead of the commonly used Taylor expansion of the square of the reflection traveltime in terms of the square of the offset, we expand the square of the reflection traveltime in terms of weak-anisotropy (WA) parameters. No acoustic approximation is used. We specify the formulas designed for anisotropy of arbitrary symmetry for the transversely isotropic (TI) media with the axis of symmetry oriented arbitrarily in the 3D space. Resulting formulas depend on three P-wave WA parameters specifying the TI symmetry and two angles specifying the orientation of the axis of symmetry. Tests of the accuracy of the more accurate of the approximate formulas indicate that maximum relative errors do not exceed 0.3% or 2.5% for weak or moderate P-wave anisotropy, respectively.

Numerical simulation of azimuthal acoustic logging in a borehole penetrating a rock formation boundary

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D283-D291 ◽  
Author(s):  
Peng Liu ◽  
Wenxiao Qiao ◽  
Xiaohua Che ◽  
Xiaodong Ju ◽  
Junqiang Lu ◽  
...  
Keyword(s):  
Domain Model ◽  
P Wave ◽  
S Wave ◽  
Angle Variation ◽  
P Waves ◽  
Acoustic Logging ◽  
S Waves ◽  
Rock Formation ◽  
Logging Tool ◽  
Difference Time

We have developed a new 3D acoustic logging tool (3DAC). To examine the azimuthal resolution of 3DAC, we have evaluated a 3D finite-difference time-domain model to simulate a case in which the borehole penetrated a rock formation boundary when the tool worked at the azimuthal-transmitting-azimuthal-receiving mode. The results indicated that there were two types of P-waves with different slowness in waveforms: the P-wave of the harder rock (P1) and the P-wave of the softer rock (P2). The P1-wave can be observed in each azimuthal receiver, but the P2-wave appears only in the azimuthal receivers toward the softer rock. When these two types of rock are both fast formations, two types of S-waves also exist, and they have better azimuthal sensitivity compared with P-waves. The S-wave of the harder rock (S1) appears only in receivers toward the harder rock, and the S-wave of the softer rock (S2) appears only in receivers toward the softer rock. A model was simulated in which the boundary between shale and sand penetrated the borehole but not the borehole axis. The P-wave of shale and the S-wave of sand are azimuthally sensitive to the azimuth angle variation of two formations. In addition, waveforms obtained from 3DAC working at the monopole-transmitting-azimuthal-receiving mode indicate that the corresponding P-waves and S-waves are azimuthally sensitive, too. Finally, we have developed a field example of 3DAC to support our simulation results: The azimuthal variation of the P-wave slowness was observed and can thus be used to reflect the azimuthal heterogeneity of formations.

On the Green function of the Lord–Shulman thermoelasticity equations

2019 ◽  
Vol 220 (1) ◽  
pp. 393-403 ◽  
Author(s):  
Zhi-Wei Wang ◽  
Li-Yun Fu ◽  
Jia Wei ◽  
Wanting Hou ◽  
Jing Ba ◽  
...  
Keyword(s):  
Green’S Function ◽  
Green's Function ◽  
Parabolic Type ◽  
Second Order ◽  
P Wave ◽  
Thermal Mode ◽  
S Wave ◽  
Order Tensor ◽  
P Waves ◽  
Thermal Conductivities

SUMMARY Thermoelasticity extends the classical elastic theory by coupling the fields of particle displacement and temperature. The classical theory of thermoelasticity, based on a parabolic-type heat-conduction equation, is characteristic of an unphysical behaviour of thermoelastic waves with discontinuities and infinite velocities as a function of frequency. A better physical system of equations incorporates a relaxation term into the heat equation; the equations predict three propagation modes, namely, a fast P wave (E wave), a slow thermal P wave (T wave), and a shear wave (S wave). We formulate a second-order tensor Green's function based on the Fourier transform of the thermodynamic equations. It is the displacement–temperature solution to a point (elastic or heat) source. The snapshots, obtained with the derived second-order tensor Green's function, show that the elastic and thermal P modes are dispersive and lossy, which is confirmed by a plane-wave analysis. These modes have similar characteristics of the fast and slow P waves of poroelasticity. Particularly, the thermal mode is diffusive at low thermal conductivities and becomes wave-like for high thermal conductivities.

Sign in / Sign up

Export Citation Format

Share Document