Weak-anisotropy approximation of P-wave geometrical spreading in horizontally layered anisotropic media of arbitrary symmetry: TTI specification

Geophysics ◽  
2021 ◽  
pp. 1-61
Author(s):  
Veronique Farra ◽  
Ivan Psencik
Keyword(s):  
Surface Profile ◽  
Seismic Wave Propagation ◽  
Multilayered Structure ◽  
P Wave ◽  
Seismic Exploration ◽  
Finite Limit ◽  
Geometrical Spreading ◽  
Weak Anisotropy ◽  
Reference Velocity ◽  
Relative Errors

Understanding the role of geometrical spreading and estimating its effects on seismic wave propagation play an important role in several techniques used in seismic exploration. The spreading can be estimated through dynamic ray tracing or determined from reflection traveltime derivatives. In the latter case, derivatives of non-hyperbolic moveout approximations are often used. We offer an alternative approach based on the weak-anisotropy approximation. The resulting formula is applicable to P-waves reflected from the bottom of a stack of horizontal layers, in which each layer can be of arbitrary anisotropy. At an arbitrary surface point, the formula depends, in each layer, on the thickness of the layer, on the P-wave reference velocity used for the construction of reference rays, and on nine P-wave weak-anisotropy (WA) parameters specifying the layer anisotropy. Along an arbitrary surface profile, the number of WA parameters reduces to five parameters related to the profile. WA parameters represent an alternative to the elastic moduli, and as such can be used for the description of any anisotropy. The relative error of the approximate formula for a multilayered structure consisting of layers of anisotropy between 8% and 20% is, at most, 10%. For models including layers of anisotropy stronger than 20%, the relative errors may reach, locally, even 30%. For any offset, relative errors remain under a finite limit, which varies with anisotropy strength.

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.

Analytic study of the geometrical spreading of P-waves in a layered transversely isotropic medium with a vertical symmetry axis

Geophysics ◽  
2000 ◽  
Vol 65 (4) ◽  
pp. 1305-1315 ◽  
Author(s):  
Hongbo Zhou ◽  
George A. McMechan
Keyword(s):  
Isotropic Medium ◽  
Symmetry Axis ◽  
Transverse Isotropy ◽  
Analytical Formula ◽  
Transversely Isotropic ◽  
P Wave ◽  
Geometrical Spreading ◽  
Transversely Isotropic Medium ◽  
Weak Anisotropy ◽  
Special Cases

An analytical formula for geometrical spreading is derived for a horizontally layered transversely isotropic medium with a vertical symmetry axis (VTI). With this expression, geometrical spreading can be determined using only the anisotropy parameters in the first layer, the traveltime derivatives, and the source‐receiver offset. Explicit, numerically feasible expressions for geometrical spreading are obtained for special cases of transverse isotropy (weak anisotropy and elliptic anisotropy). Geometrical spreading can be calculated for transversly isotropic (TI) media by using picked traveltimes of primary nonhyperbolic P-wave reflections without having to know the actual parameters in the deeper subsurface; no ray tracing is needed. Synthetic examples verify the algorithm and show that it is numerically feasible for calculation of geometrical spreading. For media with a few (4–5) layers, relative errors in the computed geometrical spreading remain less than 0.5% for offset/depth ratios less than 1.0. Errors that change with offset are attributed to inaccuracy in the expression used for nonhyberbolic moveout. Geometrical spreading is most sensitive to errors in NMO velocity, followed by errors in zero‐offset reflection time, followed by errors in anisotropy of the surface layer. New relations between group and phase velocities and between group and phase angles are shown in appendices.

P-wave reflection-moveout approximation for horizontally layered media of arbitrary moderate anisotropy

Geophysics ◽  
2020 ◽  
Vol 85 (2) ◽  
pp. C61-C70 ◽  
Author(s):  
Véronique Farra ◽  
Ivan Pšenčík
Keyword(s):  
Coordinate System ◽  
Transverse Isotropy ◽  
Local Coordinate System ◽  
Surface Profile ◽  
Taylor Series Expansion ◽  
Layered Media ◽  
P Wave ◽  
Orthorhombic Symmetry ◽  
Weak Anisotropy ◽  
Cartesian Coordinate

We present an approximate nonhyperbolic P-wave moveout formula applicable to horizontally layered media of moderate or weak anisotropy of arbitrary symmetry and orientation. Anisotropy symmetry and its orientation may differ from layer to layer. Instead of commonly used Taylor-series expansion of the square of the reflection traveltime in terms of the square of the offset, we use the weak-anisotropy approximation, in which the square of the reflection traveltime is expanded in terms of weak-anisotropy (WA) parameters. The resulting formula is simple, and it provides a transparent relation between the traveltimes and WA parameters. Along an arbitrarily chosen single surface profile, it depends, in each layer, on the thickness of the layer, on the reference P-wave velocity used for the construction of reference rays, and on three WA parameters specified in the Cartesian coordinate system related to the profile. In each layer, these three “profile” WA parameters depend on “local” WA parameters specifying anisotropy of a given layer in a local coordinate system and on directional cosines specifying the orientation of the local coordinate system with respect to the profile one. The number of local P-wave WA parameters may vary from three for transverse isotropy or six for orthorhombic symmetry to nine for triclinic symmetry. Our tests of the accuracy indicate that the maximum relative traveltime errors do not exceed 0.5% or 2.5% for weak or moderate P-wave anisotropy, respectively.

Geometrical spreading of P-waves in horizontally layered, azimuthally anisotropic media

Geophysics ◽  
2005 ◽  
Vol 70 (5) ◽  
pp. D43-D53 ◽  
Author(s):  
Xiaoxia Xu ◽  
Ilya Tsvankin ◽  
Andrés Pech
Keyword(s):  
Transversely Isotropic ◽  
Anisotropic Media ◽  
Azimuthal Velocity ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Geometrical Spreading ◽  
Weak Anisotropy ◽  
P Waves ◽  
Anisotropic Models ◽  
Wide Range

For processing and inverting reflection data, it is convenient to represent geometrical spreading through the reflection traveltime measured at the earth's surface. Such expressions are particularly important for azimuthally anisotropic models in which variations of geometrical spreading with both offset and azimuth can significantly distort the results of wide-azimuth amplitude-variation-with-offset (AVO) analysis. Here, we present an equation for relative geometrical spreading in laterally homogeneous, arbitrarily anisotropic media as a simple function of the spatial derivatives of reflection traveltimes. By employing the Tsvankin-Thomsen nonhyperbolic moveout equation, the spreading is represented through the moveout coefficients, which can be estimated from surface seismic data. This formulation is then applied to P-wave reflections in an orthorhombic layer to evaluate the distortions of the geometrical spreading caused by both polar and azimuthal anisotropy. The relative geometrical spreading of P-waves in homogeneous orthorhombic media is controlled by five parameters that are also responsible for time processing. The weak-anisotropy approximation, verified by numerical tests, shows that azimuthal velocity variations contribute significantly to geometrical spreading, and the existing equations for transversely isotropic media with a vertical symmetry axis (VTI) cannot be applied even in the vertical symmetry planes. The shape of the azimuthally varying spreading factor is close to an ellipse for offsets smaller than the reflector depth but becomes more complicated for larger offset-to-depth ratios. The overall magnitude of the azimuthal variation of the geometrical spreading for the moderately anisotropic model used in the tests exceeds 25% for a wide range of offsets. While the methodology developed here is helpful in modeling and analyzing anisotropic geometrical spreading, its main practical application is in correcting the wide-azimuth AVO signature for the influence of the anisotropic overburden.

scholarly journals S-wave experiments for the exploration of a deep geothermal carbonate reservoir in the German Molasse Basin

2021 ◽  
Vol 9 (1) ◽  
Author(s):  
Britta Wawerzinek ◽  
Hermann Buness ◽  
Hartwig von Hartmann ◽  
David C. Tanner
Keyword(s):  
Upper Jurassic ◽  
Carbonate Reservoir ◽  
P Wave ◽  
Seismic Survey ◽  
Seismic Exploration ◽  
Molasse Basin ◽  
S Wave ◽  
Induced Earthquakes ◽  
Conversion Point ◽  
Facies Classification

AbstractThere are many successful geothermal projects that exploit the Upper Jurassic aquifer at 2–3 km depth in the German Molasse Basin. However, up to now, only P-wave seismic exploration has been carried out. In an experiment in the Greater Munich area, we recorded S-waves that were generated by the conventional P-wave seismic survey, using 3C receivers. From this, we built a 3D volume of P- to S-converted (PS) waves using the asymptotic conversion point approach. By combining the P-volume and the resulting PS-seismic volume, we were able to derive the spatial distribution of the vp/vs ratio of both the Molasse overburden and the Upper Jurassic reservoir. We found that the vp/vs ratios for the Molasse units range from 2.0 to 2.3 with a median of 2.15, which is much higher than previously assumed. This raises the depth of hypocenters of induced earthquakes in surrounding geothermal wells. The vp/vs ratios found in the Upper Jurassic vary laterally between 1.5 and 2.2. Since no boreholes are available for verification, we test our results against an independently derived facies classification of the conventional 3D seismic volume and found it correlates well. Furthermore, we see that low vp/vs ratios correlate with high vp and vs velocities. We interpret the latter as dolomitized rocks, which are connected with enhanced permeability in the reservoir. We conclude that 3C registration of conventional P-wave surveys is worthwhile.

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.

Anisotropic geometrical-spreading correction for wide-azimuth P-wave reflections

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. D161-D170 ◽  
Author(s):  
Xiaoxia Xu ◽  
Ilya Tsvankin
Keyword(s):  
Transversely Isotropic ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Fractured Reservoir ◽  
Geometrical Spreading ◽  
Spreading Factor ◽  
Velocity Anisotropy ◽  
Reflection Data ◽  
Fitting Procedure ◽  
Avo Attributes

Compensation for geometrical spreading along a raypath is one of the key steps in AVO (amplitude-variation-with-offset) analysis, in particular, for wide-azimuth surveys. Here, we propose an efficient methodology to correct long-spread, wide-azimuth reflection data for geometrical spreading in stratified azimuthally anisotropic media. The P-wave geometrical-spreading factor is expressed through the reflection traveltime described by a nonhyperbolic moveout equation that has the same form as in VTI (transversely isotropic with a vertical symmetry axis) media. The adapted VTI equation is parameterized by the normal-moveout (NMO) ellipse and the azimuthally varying anellipticity parameter [Formula: see text]. To estimate the moveout parameters, we apply a 3D nonhyperbolic semblance algorithm of Vasconcelos and Tsvankin that operates simultaneously with traces at all offsets andazimuths. The estimated moveout parameters are used as the input in our geometrical-spreading computation. Numerical tests for models composed of orthorhombic layers with strong, depth-varying velocity anisotropy confirm the high accuracy of our travetime-fitting procedure and, therefore, of the geometrical-spreading correction. Because our algorithm is based entirely on the kinematics of reflection arrivals, it can be incorporated readily into the processing flow of azimuthal AVO analysis. In combination with the nonhyperbolic moveout inversion, we apply our method to wide-azimuth P-wave data collected at the Weyburn field in Canada. The geometrical-spreading factor for the reflection from the top of the fractured reservoir is clearly influenced by azimuthal anisotropy in the overburden, which should cause distortions in the azimuthal AVO attributes. This case study confirms that the azimuthal variation of the geometrical-spreading factor often is comparable to or exceeds that of the reflection coefficient.

scholarly journals A simple P-wave polarization analysis: its application to earthquake location

1994 ◽  
Vol 37 (5) ◽  
Author(s):  
B. Alessandrini ◽  
M. Cattaneo ◽  
M. Demartin ◽  
M. Gasperini ◽  
V. Lanza
Keyword(s):  
Ground Motion ◽  
Seismic Wave ◽  
Heterogeneous Media ◽  
Seismic Wave Propagation ◽  
Least Square ◽  
P Wave ◽  
Motion Direction ◽  
Arrival Times ◽  
Hypocentral Location ◽  
Ray Path

We present a method for hypocentral location which takes into account all three components of ground motion and not only the vertical one, as it is usually done by standard least-square techniques applied to arrival times. Assuming that P-wave particle motion direction corresponds to the propagation direction of the seismic wave, we carried out a simple statistical analysis of ground motion amplitudes, carefully using three-component records. We obtained the azimuth and the emersion angle of the seismic ray, which, added to Pg and Sg arrival times, allowed us to find reliable hypocentral coordinates of some local events by means of a ray-tracing technique. We compared our locations to those obtained using a least-square technique: our polarization method's dependence on the accuracy of the model used (on the contrary, the least-square technique proved to be quite stable with respect to changes in the model's velocity parameters) led us to conclude that polarization data provide coherent information on the true ray-path and can be successfully used for both location procedures and seismic wave propagation studizs in strongly heterogeneous media.

scholarly journals Vertical to Horizontal Spectral Ratio (VHSR) Response of Seismic Wave Propagation in a Homogeneous Elastic – Poroelastic Medium Using The Spectral Finite Element Method

2021 ◽  
Vol 11 (1) ◽  
pp. 95
Author(s):  
Sudarmaji Saroji ◽  
Budi Eka Nurcahya ◽  
Nivan Ramadhan Sugiantoro
Keyword(s):  
Finite Element Method ◽  
Elastic Medium ◽  
Seismic Waves ◽  
Low Frequency ◽  
Seismic Wave Propagation ◽  
P Wave ◽  
Poroelastic Medium ◽  
Compressional Waves ◽  
Spectral Finite Element Method ◽  
Spectral Finite Element

<p>Numerical modeling of 2D seismic wave propagation using spectral finite element method to estimate the response of seismic waves passing through the poroelastic medium from a hydrocarbon reservoir has been carried out. A hybrid simple model of the elastic - poroelastic - elastic with a mesoscopic scale element size of about 50cm was created. Seismic waves which was in the form of the ricker function are generated on the first elastic medium, propagated into the poroelastic medium and then transmitted to the second elastic medium. Pororoelastic medium is bearing hydrocarbon fluid in the form of gas, oil or water. Vertical and horizontal component of velocity seismograms are recorded on all mediums. Seismograms which are recorded in the poroelastic and second elastic medium show the existence of slow P compressional waves following fast P compressional waves that do not appear on the seismogram of the first elastic medium. The slow P wave is generated when the fast P wave enters the interface of the elastic - poroelastic boundary, propagated in the poroelastic medium and is transmited to the second elastic medium. The curves of Vertical to horizontal spectrum ratio (VHSR) which are observed from seismograms recorded in the poroelastic and the second elastic medium show that the peak of VHSR values at low frequency correlated with the fluid of poroelastic reservoir. The highest VHSR value at the low frequency which is recorded on the seismogram is above the 2.5 Hz frequency for reservoirs containing gas and oil in the second elastic medium, while for the medium containing water is the highest VHSR value is below the 2.5 Hz frequency.</p>

scholarly journals Seismic AVOA Inversion for Weak Anisotropy Parameters and Fracture Density in a Monoclinic Medium

2020 ◽  
Vol 10 (15) ◽  
pp. 5136 ◽  
Author(s):  
Zijian Ge ◽  
Shulin Pan ◽  
Jingye Li
Keyword(s):  
Fractured Rock ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
P Wave ◽  
Inversion Method ◽  
Fracture Density ◽  
Weak Anisotropy ◽  
Amplitude Versus Offset ◽  
Porosity And Permeability ◽  
Transverse Isotropic

In shale gas development, fracture density is an important lithologic parameter to properly characterize reservoir reconstruction, establish a fracturing scheme, and calculate porosity and permeability. The traditional methods usually assume that the fracture reservoir is one set of aligned vertical fractures, embedded in an isotropic background, and estimate some alternative parameters associated with fracture density. Thus, the low accuracy caused by this simplified model, and the intrinsic errors caused by the indirect substitution, affect the estimation of fracture density. In this paper, the fractured rock of monoclinic symmetry assumes two non-orthogonal vertical fracture sets, embedded in a transversely isotropic background. Firstly, assuming that the fracture radius, width, and orientation are known, a new form of P-wave reflection coefficient, in terms of weak anisotropy (WA) parameters and fracture density, was obtained by substituting the stiffness coefficients of vertical transverse isotropic (VTI) background, normal, and tangential fracture compliances. Then, a linear amplitude versus offset and azimuth (AVOA) inversion method, of WA parameters and fracture density, was constructed by using Bayesian theory. Tests on synthetic data showed that WA parameters, and fracture density, are stably estimated in the case of seismic data containing a moderate noise, which can provide a reliable tool in fracture prediction.

Sign in / Sign up

Export Citation Format

Share Document