共查询到20条相似文献，搜索用时 15 毫秒
1.
Although it is widely recognized that anisotropy can have a significant influence on the focusing and positioning of migrated reflection events, conventional depth imaging methods still operate with isotropic velocity fields. Here, we present an application of a 2D migration velocity analysis (MVA) algorithm, designed for factorized v(x, z) VTI (transversely isotropic with a vertical symmetry axis) media, to an offshore data set from West Africa. By approximating the subsurface with factorized VTI blocks, it is possible to decouple the spatial variations in the vertical velocity from the anisotropic parameters with minimal a priori information. Since our method accounts for lateral velocity variation, it produces more accurate estimates of the anisotropic parameters than those previously obtained with time‐domain techniques. The values of the anellipticity parameter η found for the massive shales exceed 0.2, which confirms that ignoring anisotropy in the study area can lead to substantial imaging distortions, such as mis‐stacking and mispositioning of dipping events. While some of these distortions can be removed by using anisotropic time processing, further marked improvement in image quality is achieved by prestack depth migration with the estimated factorized VTI model. In particular, many fault planes, including antithetic faults in the shallow part of the section, are better focused by the anisotropic depth‐migration algorithm and appear more continuous. Anisotropic depth migration facilitates structural interpretation by eliminating false dips at the bottom of the section and improving the images of a number of gently dipping features. One of the main difficulties in anisotropic MVA is the need to use a priori information for constraining the vertical velocity. In this case study, we successfully reconstructed the time–depth curve from reflection data by assuming that the vertical velocity is a continuous function of depth and estimating the vertical and lateral velocity gradients in each factorized block. If the subsurface contains strong boundaries with jumps in velocity, knowledge of the vertical velocity at a single point in a layer is sufficient for our algorithm to determine all relevant layer parameters. 相似文献
2.
由所建立的三维qP波相速度表示式出发,导出并解析求解各向异性介质中的频散方程,得到三维各向异性介质中的相移算子,进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质,发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化;文中给出的算例虽然是针对二维VTI介质的,但所提出的算法同样适用于三维TI介质. 相似文献
3.
In real strata anisotropy and viscosity extensively exists. They degraded waveforms in amplitude, resulting in which reducing of image resolution. To obtain high-precision imaging of deep reservoirs, we extended the separated viscous and anisotropic reverse time migration (RTM) to a stable viscoacoustic anisotropic RTM for vertical transverse isotropic (VTI) media, based on single generalized standard and linear solid (GSLS) media theory.. We used a pseudo-spectral method to develop the numerical simulation. By introducing a regularization operator to eliminate the high-frequency instability problem, we built a stable inverse propagator and achieved viscoacoustic VTI media RTM. High-resolution imaging results were obtained after correcting for the effects of anisotropy and viscosity. Synthetic tests verify the validity and accuracy of algorithm. 相似文献
4.
Modeling and RTM for arbitrary-order pure acoustic wave equation in VTI media using normalized pseudo-analytical method 
下载免费PDF全文
下载免费PDF全文 The conventional pseudo-acoustic wave equations(PWEs) in vertical transversely isotropic(VTI)media may generate SV-wave artifacts and propagation instabilities when anisotropy parameters cannot satisfy the pseudo-acoustic assumption. One solution to these issues is to use pure acoustic anisotropic wave equations, which can produce stable and pure P-wave responses without any SVwave pollutions. The commonly used pure acoustic wave equations(PAWEs) in VTI media are mainly derived from the decoupled P-SV dispersion relation based on first-order Taylor-series expansion(TE), thus they will suffer from accuracy loss in strongly anisotropic media. In this paper, we adopt arbitrary-order TE to expand the square root term in Alkhalifah's accurate acoustic VTI dispersion relation and solve the corresponding PAWE using the normalized pseudoanalytical method(NPAM) based on optimized pseudodifferential operator. Our analysis of phase velocity errors indicates that the accuracy of our new expression is perfectly acceptable for majority anisotropy parameters. The effectiveness of our proposed scheme also can be demonstrated by several numerical examples and reverse-time migration(RTM) result. 相似文献
5.
Václav Bucha 《Studia Geophysica et Geodaetica》2012,56(2):533-552
We use Kirchhoff prestack depth migration to calculate migrated sections in 3-D simple anisotropic homogeneous velocity models in order to demonstrate the impact of anisotropy on migrated images. The recorded wave field is generated in models composed of two homogeneous layers separated by one either non-inclined or inclined curved interface. The anisotropy in the upper layer is triclinic. We apply Kirchhoff prestack depth migration to velocity models with different types of anisotropy: a triclinic anisotropic medium, an isotropic medium, transversely isotropic media with a horizontal (HTI) and vertical (VTI) symmetry axis. We observe asymmetry in migration caused by triclinic anisotropy and we show the errors of the migrated interface caused by inaccurate velocity models used for migration. The study is limited to P-waves. 相似文献
6.
Imaging diffracted waves can provide useful information about complex subsurface geology and fracture networks. Separation of diffractions from typically more intensive reflected events can be done based on specularity, which measures deviation from Snell’s law. Here, we analyze two formulations of specularity and their applicability to diffraction processing in the presence of anisotropy. We show that the most common definition of specularity, originally introduced for pure modes in isotropic media, remains valid for both pure and converted waves in arbitrarily anisotropic models. The other formulation operates directly with the difference between the slowness projections onto the reflector for the incident and reflected waves. Testing on a VTI (transversely isotropic with a vertical symmetry axis) diffraction ramp model demonstrates that both formulations produce satisfactory results for anisotropic media with appropriate tapering of the specularity gathers. Then separation and imaging of diffractions is performed for the structurally complex VTI Marmousi model. We also analyze the sensitivity of diffractions in the specularity gathers to errors in the symmetry-direction velocity and anellipticity parameter η. 相似文献
7.
Analysis of Thomsen parameters for finely layered VTI media 总被引：2，自引：0，他引：2
James G. Berryman Vladimir Y. Grechka & Patricia A. Berge 《Geophysical Prospecting》1999,47(6):959-978
Since the work of Postma and Backus, much has been learned about elastic constants in vertical transversely isotropic (VTI) media when the anisotropy is due to fine layering of isotropic elastic materials. Nevertheless, there has continued to be some uncertainty about the possible range of Thomsen's anisotropy parameters ε and δ for such media. We use both Monte Carlo studies and detailed analysis of Backus' equations for both two- and three-component layered media to establish the results presented. We show that ε lies in the range ?3/8 ε ½[〈v2p〉〈v?2p〉?1], for finely layered media; smaller positive and all negative values of ε occur for media with large fluctuations in the Lamé parameter λ in the component layers. We show that δ can also be either positive or negative, and that for constant density media, sign (δ) = sign (〈v?2p〉 ? 〈v?2s〉〈v2s/v2p〉). Monte Carlo simulations show that among all theoretically possible random media, positive and negative δ are equally likely in finely layered media. (Of course, the δs associated with real earth materials may span some smaller subset of those that are theoretically possible, but answering this important question is beyond our present scope.) Layered media having large fluctuations in λ are those most likely to have positive δ. This is somewhat surprising since ε is often negative or a small positive number for such media, and we have the general constraint that ε ? δ > 0 for layered VTI media. Since Gassmann's results for fluid-saturated porous media show that the mechanical effects of fluids influence only the Lamé parameter λ, not the shear modulus μ, these results suggest that small positive δ occurring together with small positive ε (but somewhat larger than δ) may be indicative of changing fluid content in a layered earth. 相似文献
8.
A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity VP0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces. 相似文献
9.
Pure-mode wave propagation is important for applications ranging from imaging to avoiding parameter tradeoff in waveform inversion. Although seismic anisotropy is an elastic phenomenon, pseudo-acoustic approximations are routinely used to avoid the high computational cost and difficulty in decoupling wave modes to obtain interpretable seismic images. However, such approximations may result in inaccuracies in characterizing anisotropic wave propagation. We propose new pure-mode equations for P- and S-waves resulting in an artefact-free solution in transversely isotropic medium with a vertical symmetry axis. Our approximations are more accurate than other known approximations as they are not based on weak anisotropy assumptions. Therefore, the S-wave approximation can reproduce the group velocity triplications in strongly anisotropic media. The proposed approximations can be used for accurate modelling and imaging of pure P- and S-waves in transversely isotropic media. 相似文献
10.
A new wave equation is derived for modelling viscoacoustic wave propagation in transversely isotropic media under acoustic transverse isotropy approximation. The formulas expressed by fractional Laplacian operators can well model the constant-Q (i.e. frequency-independent quality factor) attenuation, anisotropic attenuation, decoupled amplitude loss and velocity dispersion behaviours. The proposed viscoacoustic anisotropic equation can keep consistent velocity and attenuation anisotropy effects with that of qP-wave in the constant-Q viscoelastic anisotropic theory. For numerical simulations, the staggered-grid pseudo-spectral method is implemented to solve the velocity–stress formulation of wave equation in the time domain. The constant fractional-order Laplacian approximation method is used to cope with spatial variable-order fractional Laplacians for efficient modelling in heterogeneous velocity and Q media. Simulation results for a homogeneous model show the decoupling of velocity dispersion and amplitude loss effects of the constant-Q equation, and illustrate the influence of anisotropic attenuation on seismic wavefields. The modelling example of a layered model illustrates the accuracy of the constant fractional-order Laplacian approximation method. Finally, the Hess vertical transversely isotropic model is used to validate the applicability of the formulation and algorithm for heterogeneous media. 相似文献
11.
从含Thomsen各向异性参数的qP波相速度表示式出发,建立并求解三维VTI介质中的频散方程,得到三维VTI介质中的相移算子,进而将以相移算子为基础的最优分裂Fourier方法推广到三维VTI介质,发展了一个三维VTI介质的深度偏移方法.文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中.文中提出的偏移算法对相移法引入了高阶校正项来补偿介质横向变化的影响,使该方法可应用于横向非均匀VTI介质的陡角度成像,文中给出的偏移脉冲响应很好地证明了这一点. 相似文献
12.
On shear-wave triplication in transversely isotropic media 总被引：1，自引：0，他引：1
The exact solution to the problem of qSV triplication in homogenous transversely isotropic media has been long known, but the result is algebraically complex and is seldom applied in practice. We present an appropriate approximation (not assuming weak qSV-anisotropy) that simplifies the conditions for the onset of off-axis triplication as anisotropy is increased, identifying the anisotropy parameter σ as the controlling parameter. It follows that commonly reported surface-seismic P-wave move-out measurements imply that many formations in the earth's sedimentary crust support off-axis qSV triplications. For typical Vp/Vs velocity ratios and a horizontally stratified earth, however, off-axis qSV triplications appear to only occur for shear-wave incidence angles too far from the vertical to be sampled by surface-seismic converted-wave survey geometries. 相似文献
13.
14.
In recent years, wave‐equation imaged data are often presented in common‐image angle‐domain gathers as a decomposition in the scattering angle at the reflector, which provide a natural access to analysing migration velocities and amplitudes. In the case of anisotropic media, the importance of angle gathers is enhanced by the need to properly estimate multiple anisotropic parameters for a proper representation of the medium. We extract angle gathers for each downward‐continuation step from converting offset‐frequency planes into angle‐frequency planes simultaneously with applying the imaging condition in a transversely isotropic with a vertical symmetry axis (VTI) medium. The analytic equations, though cumbersome, are exact within the framework of the acoustic approximation. They are also easily programmable and show that angle gather mapping in the case of anisotropic media differs from its isotropic counterpart, with the difference depending mainly on the strength of anisotropy. Synthetic examples demonstrate the importance of including anisotropy in the angle gather generation as mapping of the energy is negatively altered otherwise. In the case of a titled axis of symmetry (TTI), the same VTI formulation is applicable but requires a rotation of the wavenumbers. 相似文献
15.
We have developed new basic theories for calculating the conversion point and the travel time of the P-SV converted wave (C-wave)
in anisotropic, inhomogeneous media. This enables the use of conventional procedures such as semblance analysis, Dix-type
model building and Kirchhoff summation, to implement anisotropic processing, and makes anisotropic processing affordable.
Here we present these new developments in two parts: basic theory and application to velocity analysis and parameter estimation.
This part deals with the basic theory, including both conversion-point calculation and moveout analysis.
Existing equations for calculating the PS-wave (C-wave) conversion point in layered media with vertical transverse isotropy
(VTI) are strictly limited to offsets about half the reflector depth (an offset-depth ratio, xlz, of 0.5), and those for calculating the C-wave traveltimes are limited to offsets equal to the reflector depth (x/z=l.0). In contrast, the new equations for calculating the conversion-point extend into offsets about three-times the reflector
depth (x/z=3.0), those for calculating the C-wave traveltimes extend into offsets twice the reflector depth (x/z=2.0). With the improved accuracy, the equations can help in C-wave data processing and parameter estimation in anisotropic,
inhomogeneous media.
This work is funded by the Edinburgh Anisotropy Project (EAP) of the British Geological Survey.
First author:
Xiangyang Li, Mr. Li is currently a professorial research seismologist (Grade 6) and technical director of the Edinburgh Anisotropy Project
in the British Geological Survey. He also holds a honorary professorship in multicomponent seismology at the School of Geosciences,
University of Edinburgh. He received his BSc(1982) in Geophysics from Changchun Geological Institute, China, an MSc (1984)
in applied geophysics from East China Petroleum Institute (now known as the China University of Petroleum), and a PhD (1992)
in seismology from the University of Edinburgh. During 1984–1987, he worked as a lecturer with the East China Petroleum Institute.
Since 1991, he has been employed by the British Geological Survey. His research interests include seismic anisotropy and multicomponent
seismology. 相似文献
16.
假设横向各向同性(TI)介质的对称轴是垂直的(VTI)或者水平的(HTI)能给实际资料处理带来便利,然而实际TI介质的对称轴往往是倾斜的(TTI),忽略对称轴倾角可能给各向异性参数提取和成像带来偏差,因此需要研究是否能、以及什么条件下能忽略TTI介质对称轴倾角.本文通过理论研究和数值分析研究了与TTI介质弹性性质最接近的VTI介质(OAVTI)的弹性常数和各向异性参数与原TTI介质的弹性常数和各向异性参数之间的联系与差别.结果表明:OAVTI介质各向异性参数与原TTI介质各向异性参数之间的差别可统一表示成F(α0/β0,ε,δ,γ)ξ2的形式,其中F(α0/β0,ε,δ,γ)是无量纲各向异性参数(ε, δ, γ)的线性函数,ξ是对称轴倾角;ξ的大小对各参数的误差起主导作用,一般不建议忽略20°~25°以上的对称轴倾角;当ξ较小时,即使是对强各向异性的TTI介质作VTI近似,引起的P波各向异性参数误差也很小,因此在纵波资料处理中忽略TTI介质对称轴倾角通常是可行的;即使在小ξ条件下,倾斜对称轴对SV波也有显著影响,因此在转换波资料处理中,不建议忽略TTI介质的对称轴倾角.本文的研究为分析忽略TTI介质对称轴倾角的可行性提供了理论依据和简便的判据. 相似文献
17.
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. 相似文献
18.
Tariq Alkhalifah 《Geophysical Prospecting》2002,50(4):373-382
A linearized eikonal equation is developed for transversely isotropic (TI) media with a vertical symmetry axis (VTI). It is linear with respect to perturbations in the horizontal velocity or the anisotropy parameter η. An iterative linearization of the eikonal equation is used as the basis for an algorithm of finite-difference traveltime computations. A practical implementation of this iterative technique is to start with a background model that consists of an elliptically anisotropic, inhomogeneous medium, since traveltimes for this type of medium can be calculated efficiently using eikonal solvers, such as the fast marching method. This constrains the perturbation to changes in the anisotropy parameter η (the parameter most responsible for imaging improvements in anisotropic media). The iterative implementation includes repetitive calculation of η from traveltimes, which is then used to evaluate the perturbation needed for the next round of traveltime calculations using the linearized eikonal equation. Unlike isotropic media, interpolation is needed to estimate η in areas where the traveltime field is independent of η, such as areas where the wave propagates vertically.
Typically, two to three iterations can give sufficient accuracy in traveltimes for imaging applications. The cost of each iteration is slightly less than the cost of a typical eikonal solver. However, this method will ultimately provide traveltime solutions for VTI media. The main limitation of the method is that some smoothness of the medium is required for the iterative implementation to work, especially since we evaluate derivatives of the traveltime field as part of the iterative approach. If a single perturbation is sufficient for the traveltime calculation, which may be the case for weak anisotropy, no smoothness of the medium is necessary. Numerical tests demonstrate the robustness and efficiency of this approach. 相似文献
Typically, two to three iterations can give sufficient accuracy in traveltimes for imaging applications. The cost of each iteration is slightly less than the cost of a typical eikonal solver. However, this method will ultimately provide traveltime solutions for VTI media. The main limitation of the method is that some smoothness of the medium is required for the iterative implementation to work, especially since we evaluate derivatives of the traveltime field as part of the iterative approach. If a single perturbation is sufficient for the traveltime calculation, which may be the case for weak anisotropy, no smoothness of the medium is necessary. Numerical tests demonstrate the robustness and efficiency of this approach. 相似文献
19.
Marijan Dravinski 《地震工程与结构动力学》2003,32(5):653-670
Scattering of plane harmonic waves by a three‐dimensional basin of arbitrary shape embedded within elastic half‐space is investigated by using an indirect boundary integral equation approach. The materials of the basin and the half‐space are assumed to be the most general anisotropic, homogeneous, linearly elastic solids without any material symmetry (i.e. triclinic). The unknown scattered waves are expressed in terms of three‐dimensional triclinic time harmonic full‐space Green's functions. The results have been tested by comparing the surface response of semi spherical isotropic and transversely isotropic basins for which the numerical solutions are available. Surface displacements are presented for a semicircular basin subjected to a vertical incident plane harmonic pseudo‐P‐, SV‐, or SH‐wave. These results are compared with the motion obtained for the corresponding equivalent isotropic models. The results show that presence of the basin may cause significant amplification of ground motion when compared to the free‐field displacements. The peak amplitude of the predominant component of surface motion is smaller for the anisotropic basin than for the corresponding isotropic one. Anisotropic response may be asymmetric even for symmetric geometry and incidence. Anisotropic surface displacement generally includes all three components of motion which may not be the case for the isotropic results. Furthermore, anisotropic response strongly depends upon the nature of the incident wave, degree of material anisotropy and the azimuthal orientation of the observation station. These results clearly demonstrate the importance of anisotropy in amplification of surface ground motion. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献