Inhomogeneous waves in isotropic anelastic media: explicit expressions for Q

Author(s):  
Xu Liu ◽  
Stewart Greenhalgh ◽  
Bing Zhou ◽  
Huijian Li

<p>Seismic waves propagating in attenuative materials are generally inhomogenous waves which, unlike homogeneous waves, have different directions of propagation and attenuation. The degree of wave inhomogeneity can be represented by the inhomogeneity parameter D which varies from 0 to infinity (Cerveny & Psencik, 2005). The dissipation (1/Q)  factors of   inhomogeneous waves vary according to the different definitions. Based on the complex energy balance equations (Carcione, 2001) and the mixed specification of the slowness vector (Cerveny & Psencik, 2005), explicit formulas for the dissipation factors of P- and SV-waves are developed under the two different definitions, (1) 1/Q<sub>V</sub>, the ratio of the time-averaged dissipated energy density to the time-averaged strain-energy density, and (2) 1/Q<sub>T</sub>, the time-averaged dissipated energy density to the time-averaged energy density. By setting the degree of wave inhomogeneity D as zero, these dissipation factor expressions are reduced to their special case versions as homogeneous waves, i.e., 1/Q<sub>VH</sub> = -Im(M)/Re(M) and 1/Q<sub>TH</sub> = 2α<em>v</em>/ω , where, M is the wave modulus, α the attenuation coefficient, <em>v</em> the phase velocity and ω the frequency. An example viscoelastic material is chosen to represent the dissipative features of a reservoir for which P-waves are normally more dissipative than S-waves. The calculated dissipation factors of P-waves under the two definitions (i.e. 1/Q<sub>PV</sub> and 1/Q<sub>PT</sub>) decrease with increasing degree of wave inhomogeneity. For the counterpart S waves, 1/Q<sub>SV</sub> is independent of the degree of wave inhomogeneity and 1/Q<sub>ST</sub> shows the trend of increasing with increasing degree of wave inhomogeneity. These findings can be explained by the limiting dissipation factors (defined at the infinite degree of inhomogeneity) which all depend only on the shear modulus.  To ensure the correctness of our results, we repeated each step of the investigation  in a parallel way based on Buchen’s (1971) classic real value energy balance equation, including derivation of explicit formulas for 1/Q<sub>PV</sub> and 1/Q<sub>PT</sub> , with inhomogeneity angle γ  ( -π/2 < γ < π/2) representing the degree of inhomogeneity of the plane wave. We also obtain the inhomogeneity-independent formula for 1/Q<sub>SV</sub>, and  exactly the same phase velocity and dissipation factor dispersion results for the example material.</p><p><strong>Acknowledgements</strong></p><p>We are grateful to the College of Petroleum Engineering & Geosciences, King Fahd University of Petroleum and Minerals, Kingdom of Saudi Arabia for supporting this research.</p><p><strong>References</strong></p><p>Buchen, P.W., 1971. Plane waves in linear viscoelastic media, Geophysical Journal of the Royal Astronomical Society, 23, 531-542.</p><p>Carcione, J. M., 2001. Wave fields in real media:Wave propagation in anisotropic, anelastic and porous media: Pergamon Press, Inc.</p><p>Cerveny, V. & Psencik, I., 2005. Plane waves in viscoelastic anisotropic media—I. Theory, Geophysical Journal International, 161, 197–212.</p>

2019 ◽  
Vol 17 (2) ◽  
pp. 300-312 ◽  
Author(s):  
Xu Liu ◽  
Stewart Greenhalgh ◽  
Bing Zhou ◽  
Huijian Li

Abstract We derive explicit expressions for the dissipation factors of inhomogeneous P and SV-waves in isotropic viscoelastic media. The Q−1 values are given as concise and simple functions of material parameters and the wave inhomogeneity parameter using two different definitions. Unlike homogenous waves, inhomogeneous waves may have significant differences in the values of dissipation factors because of different definitions. For example, under one of the three dissipation factor definitions that Q−1 is equal to the time-averaged dissipated-energy density divided by twice the time-averaged strain-energy density, it is found and proved that the dissipation factor of SV-waves is totally independent of the inhomogeneity parameter. For materials in which P-waves are normally more dissipative than S-waves (e.g. a porous reservoir), the dissipation factors of P-waves tend to decrease with increasing degree of inhomogeneity. Based on Buchan's classic real value energy balance equation, a parallel investigation is conducted for each step similar to that based on the Carcione equations, including derivation of explicit formulas (with inhomogeneity angle representing the degree of inhomogeneity of a plane wave), and dissipation curves calculations. We also obtain an inhomogeneity independent formula of $Q_{\, SV}^{ - 1}$, and exactly the same phase velocity and attenuation dispersion results for the example material.


Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. T209-T224 ◽  
Author(s):  
Xu Liu ◽  
Stewart Greenhalgh ◽  
José M. Carcione

Seismic waves in an attenuative porous medium are generally inhomogeneous waves, which have different directions of propagation and attenuation. The dissipation factors (1/ Q) of inhomogeneous waves are strongly dependent on the degree of wave inhomogeneity and cannot be expressed correctly with the usual 1/ Q expressions valid only for homogeneous waves. We have used the differing definitions of 1/ Q for inhomogeneous waves (i.e., the ratio of the time-averaged dissipated energy density to the time-averaged strain energy density or time-averaged total energy density) and the complex form of the energy balance equations of poroviscoelastic media to derive concise and explicit expressions for the dissipation factors. They are given as simple functions of the material parameters and the wave inhomogeneity parameter for inhomogeneous SV-waves and fast and slow P-waves. The isotropic, poroviscoelastic medium under consideration is upscaled from effective Biot theory for a double-porosity solid, which is the most general theory to describe wave propagation in a reservoir. We find that, if the inhomogeneity parameter is infinite (i.e., the inhomogeneity angle is 90°) for all three Biot waves, then the dissipation factors only depend on the ratio of the imaginary to the real part of the complex shear modulus. Our explicit expressions for the dissipation factors of poroviscoelastic materials also are reduced to obtain their counterparts for viscoelastic media as a special case. The inhomogeneous waves in an example poroviscoelastic material are used to demonstrate that the 1/ Q values of the three Biot waves strongly depend on the inhomogeneity parameters and furthermore the different definitions may cause significant differences of 1/ Q values. We find that the dissipation factor of fast P-waves may decrease with the increasing degree of inhomogeneity, which contradicts previously published results.


2010 ◽  
Vol 2010 ◽  
pp. 1-11 ◽  
Author(s):  
Inder Singh ◽  
Dinesh Kumar Madan ◽  
Manish Gupta

3D solutions of the dynamical equations in the presence of external forces are derived for a homogeneous, prestressed medium. 2D plane waves solutions are obtained from general solutions and show that there exist two types of plane waves, namely, quasi-P waves and quasi-SV waves. Expressions for slowness surfaces and apparent velocities for these waves are derived analytically as well as numerically and represented graphically.


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D283-D291 ◽  
Author(s):  
Peng Liu ◽  
Wenxiao Qiao ◽  
Xiaohua Che ◽  
Xiaodong Ju ◽  
Junqiang Lu ◽  
...  

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.


Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 955-967 ◽  
Author(s):  
Hiroshi Amano

An analytical solution to separate P‐waves and S‐waves in vertical seismic profile (VSP) wavefields is derived using combinations of certain terms of the formal solution for forward VSP modeling. Some practical applications of this method to synthetic seismograms and field data are investigated and evaluated. Little wave distortion is recognized, and the weak wavefield masked by dominant wavetrains can be extracted with this method. The decomposed wavefield is expressed in the frequency‐depth (f-z) domain as a linear combination of up to the third‐order differential of traces, which is approximated by trace differences in the practical separation process. In general, five traces with single‐component data are required in this process, but the same process is implemented with only three traces in the acoustic case. Two‐trace extrapolation is applied to each edge of the data gather to enhance the accuracy of trace difference. Since the formulas are developed in the f-z domain, the influence of anelasticity can be taken into account, and the calculation is carried out fast enough with the benefit of the fast Fourier transform (FFT).


Author(s):  
N. H. Scott

Inhomogeneous small-amplitude plane waves of (complex) frequency ω are propagated through a linear dissipative material which displays hereditary viscoelasticity. The energy density, energy flux and dissipation are quadratic in the small quantities, namely, the displacement gradient, velocity and velocity gradient, each harmonic with frequency ω , and so give rise to attenuated constant terms as well as to inhomogeneous plane waves of frequency 2 ω . The quadratic terms are usually removed by time averaging but we retain them here as they are of comparable magnitude with the time-averaged quantities of frequency ω . A new relationship is derived in hereditary viscoelasticity that connects the amplitudes of the terms of the energy density, energy flux and dissipation that have frequency 2 ω . It is shown that the complex group velocity is related to the amplitudes of the terms with frequency 2 ω rather than to the attenuated constant terms as it is for homogeneous waves in conservative materials.


Sign in / Sign up

Export Citation Format

Share Document