共查询到20条相似文献，搜索用时 33 毫秒
1.
—Approximate PP plane wave displacement coefficients of reflection and transmission for weak contrast interfaces separating weakly but arbitrarily anisotropic elastic media are presented. The PP reflection coefficient for such an interface has been derived recently by Vavry?uk and P?en?ík (1997). The PP transmission coefficient presented in this paper was derived by the same approach. The coefficients are given as a sum of the coefficient for the weak contrast interface separating two nearby isotropic media and a term depending linearly on contrasts of the so-called weak anisotropy (WA) parameters (parameters specifying deviation of properties of the medium from isotropy), across the interface. While the reflection coefficient depends only on 8 of the complete set of the WA parameters describing P-wave phase velocity in weakly anisotropic media, the transmission coefficient depends on their complete set. The PP reflection coefficient depends on "shear-wave splitting parameter" γ. Tests of accuracy of the approximate formulae are presented on several models. 相似文献
2.
Weak‐anisotropy approximation for P‐wave reflection coefficient at the boundary between two tilted transversely isotropic media 
下载免费PDF全文
下载免费PDF全文 Existing and commonly used in industry nowadays, closed‐form approximations for a P‐wave reflection coefficient in transversely isotropic media are restricted to cases of a vertical and a horizontal transverse isotropy. However, field observations confirm the widespread presence of rock beds and fracture sets tilted with respect to a reflection boundary. These situations can be described by means of the transverse isotropy with an arbitrary orientation of the symmetry axis, known as tilted transversely isotropic media. In order to study the influence of the anisotropy parameters and the orientation of the symmetry axis on P‐wave reflection amplitudes, a linearised 3D P‐wave reflection coefficient at a planar weak‐contrast interface separating two weakly anisotropic tilted tranversely isotropic half‐spaces is derived. The approximation is a function of the incidence phase angle, the anisotropy parameters, and symmetry axes tilt and azimuth angles in both media above and below the interface. The expression takes the form of the well‐known amplitude‐versus‐offset “Shuey‐type” equation and confirms that the influence of the tilt and the azimuth of the symmetry axis on the P‐wave reflection coefficient even for a weakly anisotropic medium is strong and cannot be neglected. There are no assumptions made on the symmetry‐axis orientation angles in both half‐spaces above and below the interface. The proposed approximation can be used for inversion for the model parameters, including the orientation of the symmetry axes. Obtained amplitude‐versus‐offset attributes converge to well‐known approximations for vertical and horizontal transverse isotropic media derived by Rüger in corresponding limits. Comparison with numerical solution demonstrates good accuracy. 相似文献
3.
— A P-wave tomographic method for 3-D complex media (3-D distribution of elastic parameters and curved interfaces) with orthorhombic symmetry is presented in this paper. The technique uses an iterative linear approach to the nonlinear travel-time inversion problem. The hypothesis of orthorhombic anisotropy and 3-D inhomogeneity increases the set of parameters describing the model dramatically compared to the isotropic case. Assuming a Factorized Anisotropic Inhomogeneous (FAI) medium and weak anisotropy, we solve the forward problem by a perturbation approach. We use a finite element approach in which the FAI medium is divided into a set of elements with polynomial elastic parameter distributions. Inside each element, analytical expressions for rays and travel times, valid to first-order, are given for P waves in orthorhombic inhomogeneous media. More complex media can be modeled by introducing interfaces separating FAI media with different elastic properties. Simple formulae are given for the Fréchet derivatives of the travel time with respect to the elastic parameters and the interface parameters. In the weak anisotropy hypothesis the P-wave travel times are sensitive only to a subset of the orthorhombic parameters: the six P-wave elastic parameters and the three Euler angles defining the orientation of the mirror planes of symmetry. The P-wave travel times are inverted by minimizing in terms of least-squares the misfit between the observed and calculated travel times. The solution is approached using a Singular Value Decomposition (SVD). The stability of the inversion is ensured by making use of suitable a priori information and/or by applying regularization. The technique is applied to two synthetic data sets, simulating simple Vertical Seismic Profile (VSP) experiments. The examples demonstrate the necessity of good 3-D ray coverage when considering complex anisotropic symmetry. 相似文献
4.
The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995. 相似文献
5.
6.
In an acoustic transversely isotropic medium, there are two waves that propagate. One is the P-wave and another one is the S-wave (also known as S-wave artefact). This paper is devoted to analyse the S-wave in two-dimensional acoustic transversely isotropic media with a tilted symmetry axis. We derive the S-wave slowness surface and traveltime function in a homogeneous acoustic transversely isotropic medium with a tilted symmetry axis. The S-wave traveltime approximations in acoustic transversely isotropic media with a tilted symmetry axis can be mapped from the counterparts for acoustic transversely isotropic media with a vertical symmetry axis. We consider a layered two-dimensional acoustic transversely isotropic medium with a tilted symmetry axis to analyse the S-wave moveout. We also illustrate the behaviour of the moveout for reflected S-wave and converted waves. 相似文献
7.
Tilted transversely isotropic formations cause serious imaging distortions in active tectonic areas (e.g., fold‐and‐thrust belts) and in subsalt exploration. Here, we introduce a methodology for P‐wave prestack depth imaging in tilted transversely isotropic media that properly accounts for the tilt of the symmetry axis as well as for spatial velocity variations. For purposes of migration velocity analysis, the model is divided into blocks with constant values of the anisotropy parameters ε and δ and linearly varying symmetry‐direction velocity VP0 controlled by the vertical (kz) and lateral (kx) gradients. Since determination of tilt from P‐wave data is generally unstable, the symmetry axis is kept orthogonal to the reflectors in all trial velocity models. It is also assumed that the velocity VP0 is either known at the top of each block or remains continuous in the vertical direction. The velocity analysis algorithm estimates the velocity gradients kz and kx and the anisotropy parameters ε and δ in the layer‐stripping mode using a generalized version of the method introduced by Sarkar and Tsvankin for factorized transverse isotropy with a vertical symmetry axis. Synthetic tests for several models typical in exploration (a syncline, uptilted shale layers near a salt dome and a bending shale layer) confirm that if the symmetry‐axis direction is fixed and VP0 is known, the parameters kz, kx, ε and δ can be resolved from reflection data. It should be emphasized that estimation of ε in tilted transversely isotropic media requires using nonhyperbolic moveout for long offsets reaching at least twice the reflector depth. We also demonstrate that application of processing algorithms designed for a vertical symmetry axis to data from tilted transversely isotropic media may lead to significant misfocusing of reflectors and errors in parameter estimation, even when the tilt is moderate (30°). The ability of our velocity analysis algorithm to separate the anisotropy parameters from the velocity gradients can be also used in lithology discrimination and geologic interpretation of seismic data in complex areas. 相似文献
8.
James Gaiser 《Geophysical Prospecting》2019,67(9):2414-2431
Utilizing shear-wave (S-wave) data acquired with compressional waves (P-waves) is becoming more common as joint imaging and inversion techniques improve. Interest in S-waves radiated from vertical sources and buried explosives exploits conversion to P-waves as primary reflections (SP-waves) for reducing acquisition costs and for application to legacy data. However, recent investigations overstate the extent of SP-wave illumination and show isotropic processing results with narrow bandwidth frequency and wavenumber data. I demonstrate that illumination with SP-waves is limited in general to near vertical polar angles up to around 30° or 35° for VP/VS of 2 or 3, respectively. At greater angles, S-waves are typically in the P-wave evanescent range and cannot excite SP-wave reflections. Contrary to recent claims, these sources for P-wave do not radiate SH-waves polarized in horizontal planes in all azimuths. I show these properties for isotropic media with radiation expressions for amplitude derived in vector slowness coordinates. Also, I extend these expressions to transversely isotropic media with a vertical symmetry axis to show agreement with synthetic seismic data that only quasi SV-waves are radiated and become more narrowly focused towards 45°. Furthermore, in orthorhombic media, synthetic data show that fast S1- and slow S2-waves polarized parallel and perpendicular to fractures may appear as SV- and SH-waves. For the partially saturated fracture model studied here, S1-wave radiation has broader azimuthal illumination than slow S2-waves, which are more narrowly focused in azimuth. These produce SP-wave splitting signatures on vertical component reflection data that are nearly identical to PS-wave signatures on radial horizontal component data. Separating these fast and slow SP-waves is an additional processing challenge. 相似文献
9.
Propagation in the plane of mirror symmetry of a monoclinic medium, with displacement normal to the plane, is the most general circumstance in anisotropic media for which pure shear-wave propagation can occur at all angles. Because the pure shear mode is uncoupled from the other two modes, its slowness surface in the plane is an ellipse. When the mirror symmetry plane is vertical the pure shear waves in this plane are SH waves and the elliptical SH sheet of the slowness surface is, in general, tilted with respect to the vertical axis. Consider a half-space of such a monoclinic medium, called medium M, overlain by a half-space of isotropic medium I with plane SH waves incident on medium M propagating in the vertical symmetry plane of M. Contrary to the appearance of a lack of symmetry about the vertical axis due to the tilt of the SH-wave slowness ellipse, the reflection and transmission coefficients are symmetrical functions of the angle of incidence, and further, there exists an isotropic medium E with uniquely determined density and shear speed which gives exactly the same reflection and transmission coefficients underlying medium J as does monoclinic medium M. This means that the underlying monoclinic medium M can be replaced by isotropic medium E without changing the reflection and transmission coefficients for all values of the angle of incidence. Thus no set of SH seismic experiments performed in the isotropic medium in the symmetry plane of the underlying half-space can reveal anything about the monoclinic anisotropy of that underlying half-space. Moreover, even when the underlying monoclinic half-space is stratified, there exists a stratified isotropic half-space that gives the identical reflection coefficient as the stratified monoclinic half-space for all angles of incidence and all frequencies. 相似文献
10.
We obtain the reflection and transmission coefficients at a fracture in transversely isotropic media, whose symmetry axis is perpendicular to the fracture surface. We consider dissimilar upper and lower media. The fracture is modeled as an interface with boundary discontinuities in the displacement and the particle velocity. The stress components are proportional to the displacement and velocity discontinuities through the specific stiffnesses and specific viscosities, respectively. The stiffness introduces frequency-dependence and phase changes in the fracture response and the viscosity is related to the energy loss. We also calculate the energy balance at the fracture and the dissipated energy. The theory is illustrated by computing the reflection coefficient of a fracture present in the Antarctic ice cap. In this case, the reflection coefficient decreases with increasing incidence angle and then approaches 1 at grazing angle. 相似文献
11.
The arrival-time curve of a reflection from a horizontal interface, beneath a homogeneous isotropic layer, is a hyperbola in the x - t-domain. If the subsurface is one-dimensionally inhomogeneous (horizontally layered), or if some or all of the layers are transversely isotropic with vertical axis of symmetry, the statement is no longer strictly true, though the arrival-time curves are still hyperbola-like. In the case of transverse isotropy, however, classical interpretation of these curves fails. Interval velocities calculated from t2 - x2-curves do not always approximate vertical velocities and therefore cannot be used to calculate depths of reflectors. To study the relationship between velocities calculated from t2 - x2-curves and the true velocities of a transversely isotropic layer, we approximate t2 - x2-curves over a vertically inhomogeneous transversely isotropic medium by a three-term Taylor series and calculate expressions for these terms as a function of the elastic parameters. It is shown that both inhomogeneity and transverse isotropy affect slope and curvature of t2 - x2-curves. For P-waves the effect of transverse isotropy is that the t2 - x2-curves are convex upwards; for SV-waves the curves are convex downwards. For SH-waves transverse isotropy has no effect on curvature. 相似文献
12.
Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures. 相似文献
13.
James G. Berryman 《Geophysical Prospecting》2009,57(2):193-208
Certain crack-influence parameters of Sayers and Kachanov are shown to be directly related to Thomsen's weak-anisotropy seismic parameters for fractured reservoirs when the crack/fracture density is small enough. These results are then applied to the problem of seismic wave propagation in polar reservoirs, i.e., those anisotropic reservoirs having two axes that are equivalent but distinct from the third axis), especially for horizontal transversely isotropic seismic wave symmetry due to the presence of aligned vertical fractures and resulting in azimuthal seismic wave symmetry at the Earth's surface. The approach presented suggests one method of inverting for fracture density from wave speed data. A significant fraction of the technical effort in the paper is devoted to showing how to predict the angular location of the true peak (or trough) of the quasi-SV-wave for polar media and especially how this peak is related to another angle that is very easy to compute. The axis of symmetry is always treated here as the x 3 -axis for either vertical transversely isotropic symmetry (due, for example, to horizontal cracks), or horizontal transversely isotropic symmetry (due to aligned vertical cracks). Then, the meaning of the stiffnesses is derived from the fracture analysis in the same way for vertical transversely isotropic and horizontal transversely isotropic media, but for horizontal transverse isotropy the wave speeds relative to the Earth's surface are shifted by 90o in the plane perpendicular to the aligned vertical fractures. Skempton's poroelastic coefficient B is used as a general means of quantifying the effects of fluids inside the fractures. Explicit Biot-Gassmann-consistent formulas for Thomsen's parameters are also obtained for either drained or undrained fractures resulting in either vertical transversely isotropic or horizontal transversely isotropic symmetry of the reservoir. 相似文献
14.
15.
Many rocks possess electrical properties with a clearly expressed anisotropy. The anisotropic character of the rocks is often overlooked in forming the Fréchet derivatives or sensitivity functions for parameter updating during the inversion of DC resistivity data. In this study we have compared the sensitivity patterns for an isotropic, homogeneous model with that for a transversely isotropic (i.e. anisotropic) model having a tilted axis of symmetry using a pole–pole array. The sensitivity functions are expressed in terms of the derivatives of the electric potential U with respect to the average conductivity σm (geometric mean of the longitudinal and transverse conductivities) and the coefficient of anisotropy λ. Results are plotted in both cross-section form and plan view for various dip and strike angles of the axis of symmetry. The derivative dU/dλ decreases more rapidly than the isotropic value dU/dσ, and shows pronounced asymmetry and weakening of magnitude with increasing dip of the plane of symmetry. The derivative dU/dσm also exhibits the asymmetric pattern (except for vertical and horizontal dip cases). The positive region between the electrodes only extends to a small depth compared to the isotropic derivative, even in the case of a vertical axis of symmetry (VTI medium). The ratio of this anisotropic derivative to the isotropic derivative, when plotted as a function of position and depth shows prominent differences in both the sign and the magnitude of the sensitivities, especially for steep dips and for strongly anisotropic rocks. The plot highlights the dangers of an isotropic assumption. Even for mildly anisotropic rocks (λ < 1.2) the possibility for error in interpretation is considerable. Combined borehole and surface measurements are needed to diagnose anisotropy. Further work is needed to design optimal electrode configurations in anisotropic situations. 相似文献
16.
Transverse isotropy with a vertical axis of symmetry is a common form of anisotropy in sedimentary basins, and it has a significant influence on the seismic amplitude variation with offset. Although exact solutions and approximations of the PP-wave reflection coefficient for the transversely isotropic media with vertical axis of symmetry have been explicitly studied, it is difficult to apply these equations to amplitude inversion, because more than three parameters need to be estimated, and such an inverse problem is highly ill-posed. In this paper, we propose a seismic amplitude inversion method for the transversely isotropic media with a vertical axis of symmetry based on a modified approximation of the reflection coefficient. This new approximation consists of only three model parameters: attribute A, the impedance (vertical phase velocity multiplied by bulk density); attribute B, shear modulus proportional to an anellipticity parameter (Thomsen's parameter ε−δ); and attribute C, the approximate horizontal P-wave phase velocity, which can be well estimated by using a Bayesian-framework-based inversion method. Using numerical tests we show that the derived approximation has similar accuracy to the existing linear approximation and much higher accuracy than isotropic approximations, especially at large angles of incidence and for strong anisotropy. The new inversion method is validated by using both synthetic data and field seismic data. We show that the inverted attributes are robust for shale-gas reservoir characterization: the shale formation can be discriminated from surrounding formations by using the crossplot of the attributes A and C, and then the gas-bearing shale can be identified through the combination of the attributes A and B. We then propose a rock-physics-based method and a stepwise-inversion-based method to estimate the P-wave anisotropy parameter (Thomsen's parameter ε). The latter is more suitable when subsurface media are strongly heterogeneous. The stepwise inversion produces a stable and accurate Thomsen's parameter ε, which is proved by using both synthetic and field data. 相似文献
17.
18.
19.
In this paper, we have considered the reflection and refraction of a plane wave at an interface between two half-spaces. The
lower half-spaces is composed of highly anisotropic triclinic crystalline material and the upper half-space is homogeneous
and isotropic. It has been assumed that due to incidence of a plane quasi-P (qP) wave, three types of waves, namely, quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH), will be generated in the lower half space whereas P and S waves will be generated in the upper half space. The phase velocities of all the quasi waves have been calculated. It has
been assumed that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation.
Some specific relations have been established between directions of motion and propagation, respectively. The expressions
for reflection coefficients of qP, qSV, qSH and refracted coefficients of P and SV waves are obtained. Results of reflection and refraction coefficients are presented. 相似文献
20.
Ilya Tsvankin 《Geophysical Prospecting》1997,45(3):479-512
The transversely isotropic (TI) model with a tilted axis of symmetry may be typical, for instance, for sediments near the flanks of salt domes. This work is devoted to an analysis of reflection moveout from horizontal and dipping reflectors in the symmetry plane of TI media that contains the symmetry axis. While for vertical and horizontal transverse isotropy zero-offset reflections exist for the full range of dips up to 90°, this is no longer the case for intermediate axis orientations. For typical homogeneous models with a symmetry axis tilted towards the reflector, wavefront distortions make it impossible to generate specular zero-offset reflected rays from steep interfaces. The ‘missing’ dipping planes can be imaged only in vertically inhomogeneous media by using turning waves. These unusual phenomena may have serious implications in salt imaging. In non-elliptical TI media, the tilt of the symmetry axis may have a drastic influence on normal-moveout (NMO) velocity from horizontal reflectors, as well as on the dependence of NMO velocity on the ray parameter p (the ‘dip-moveout (DMO) signature’). The DMO signature retains the same character as for vertical transverse isotropy only for near-vertical and near-horizontal orientation of the symmetry axis. The behaviour of NMO velocity rapidly changes if the symmetry axis is tilted away from the vertical, with a tilt of ±20° being almost sufficient to eliminate the influence of the anisotropy on the DMO signature. For larger tilt angles and typical positive values of the difference between the anisotropic parameters ε and δ, the NMO velocity increases with p more slowly than in homogeneous isotropic media; a dependence usually caused by a vertical velocity gradient. Dip-moveout processing for a wide range of tilt angles requires application of anisotropic DMO algorithms. The strong influence of the tilt angle on P-wave moveout can be used to constrain the tilt using P-wave NMO velocity in the plane that includes the symmetry axis. However, if the azimuth of the axis is unknown, the inversion for the axis orientation cannot be performed without a 3D analysis of reflection traveltimes on lines with different azimuthal directions. 相似文献