共查询到20条相似文献，搜索用时 156 毫秒
1.
2.
3.
随着多分量采集技术的发展,弹性波逆时偏移技术在三维各向异性介质复杂地质构造成像中得到了广泛的应用.然而耦合的P波场和S波场,会在传播过程中产生串扰噪声,降低弹性波逆时偏移的成像精度.为了解决这一问题,本研究针对具有倾斜各向异性对称轴的三维横向各向同性(Transverse Isotropy, TI)介质,提出了一种矢量弹性波场快速解耦方法,可以有效提高偏移剖面的成像质量.该方法首先通过坐标转换,将观测系统坐标系的垂直轴旋转到TI介质的对称轴方向,在新坐标系下,根据具有垂直对称轴的三维横向各向同性(Vertical Transverse Isotropy, VTI)介质中的分解算子,推导出三维TI介质解耦算子表达式.接着引入一种在空间域快速计算分解波场的方法,来实现空间域矢量P波场和S波场分离,极大地提高了计算效率.最后,通过点积成像条件,将提出的P/S波分解方法引入到三维TI介质弹性波逆时偏移中,得到高精度的PP和PS成像.与以往的波场分解方法相比,本文方法具有数值稳定和计算效率高的特点.数值算例表明,应用上述三维TI分解算子得到的偏移剖面有效压制了噪声,提高了成像质量. 相似文献
4.
本文基于弱各向异性（WA）介质的正反演公式和qP波的坐标变换,推导了利用qP波反演任意倾斜对称轴的横向各向同性（TTI）介质的各向异性参数和对称轴方向的公式.理论和数值实验表明,利用2个相互正交的变井源距垂直地震剖面（walkaway VSP）可以完全确定钻井中TTI介质qP波的3个WA参数和对称轴的2个方向参数.我们完成了几个由不同数量剖面组成的walkaway VSP模拟实验,使用TTI模型和一般各向异性模型对模拟数据进行了反演,证明了反演公式的正确性和可靠性.使用这些公式,对来自Java Sea的由3条剖面组成的walkaway VSP观测数据进行了各向异性反演,获得了钻井中接收点处介质的WA参数. 相似文献
5.
相速度和偏振方向是研究地震波传播规律和描述介质特性的重要参数,在理论研究和实际应用中有重要作用.本文假定倾斜横向各向同性(TTI)介质对称轴位于观测坐标系XOZ面内,在此观测坐标系下直接推导了TTI介质弹性波相速度和偏振方向的解析表达式,再进一步利用Thomsen弱各向异性理论,推导了弱各向异性近似条件下弹性波相速度以及qP波和qSV波偏振方向表达式.理论分析和数值试例表明,在相速度方面,随着各向异性介质参数改变,qP波和qSH波速度变化较为平缓,qSV波速度变化较为剧烈.弹性波相速度近似式误差均较小,能较好地近似精确相速度.在偏振方向方面,SH波偏振方向只是传播方向和对称轴倾角的函数,而与各向异性参数无关,SH波偏振方向既垂直于传播方向,又垂直于TTI介质对称轴方向.除特定方向外,qP波和qSV波的偏振方向与传播方向均成一定角度,并且随TTI介质对称轴倾角的改变而改变;在精确和近似情况下,qP波和qSV波的偏振方向始终垂直;在精度允许范围内,偏振方向的弱各向异性近似式与理论解析式吻合较好. 相似文献
6.
《地震地质》2020,(3)
完全匹配层吸收边界条件通常可以很好地吸收模型边界的地震反射波,但对于横向各向同性介质的模拟效果欠佳,且尚在发展阶段。为此,文中推导了横向各向同性介质中弹性动力学波动方程,给出了施加完全匹配层(PML)吸收边界条件下时间域二阶、空间域十阶精度的高阶交错网格的有限差分形式,并分别建立了均匀的垂直向对称轴的横向各向同性介质(VTI介质)和倾斜向对称轴的横向各向同性介质(TTI介质)模型。计算结果表明,对于对称轴为任意角度的横向各向同性介质,当PML边界层厚度达到一定的数值时,可以很好地抑制人工边界所产生的地震波反射效应,且PML的吸收效果不会被入射角与入射波频率影响。 相似文献
7.
用偏微分方程组特征理论研究层状并以垂直轴为对称轴的横向各向同性（TIV）介质的参数反演问题，首先从弹性波运动方程与TIV介质的应力－应变关系导出了平面波耦合方程组的特征型．根据奇性分析与特征积分法给出了连续情形下的特征线边界条件，连续清形下的波场延拓方程即平面波耦合方程组的特征型与特征线边界条件组成了层状TIV平滑介质、弱间断介质参数反演问题的基本方程组．并导出了间断情形下的波场延拓方程与特征线边界条件，这些方程组可用于层状介质（间断情况）的参数反演．基于这些基本方程组，探讨了利用地面多分量地震资料反演层状TIV介质多个弹性参数的问题． 相似文献
8.
地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征．各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用．为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程．本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程．文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势． 相似文献
9.
二维横各向同性弹性随机介质中的波场特征 总被引：9，自引：4，他引：5
本文通过交错网格有限差分正演．模拟了平面地震波在二维横各向同性弹性随机介质模型中的传播及其自激自收时间记录．为研究横各向同性弹性随机介质模型中的波场特征,我们在五个不同的时间区段上,分别计算剖面的三个统计特征(横向中心频率、纵向中心频率、波场能量相对值)．这样,对应每一个横各向同性弹性随机介质模型．均可计算得到15个不同的波场特征量．我们通过在二维横各向同性弹性随机介质中的正演模拟．研究当自相关长度以及介质的各向异性系数变化时,对应的上述波场特征量的变化特点．证实了在随机介质模型中．各向异性系数的变化会引起波场记录上的某些统计特征的变化,归纳得出了若干结论． 相似文献
10.
我们发展了一种模拟复杂地表下含裂缝介质地震波场的方法,这对于解释山地地区的地震资料具有重要意义。基于Coates-Schoenberg方法,把裂缝引入到有限差分法(FD)中,从而使包含裂缝的单元里的弹性介质就具有了局部的各向异性。为了模拟起伏的地表地形,我们借助于贴体网格,将笛卡尔坐标系的具有水平对称轴的横向各向同性介质(HTI)的弹性波方程和自由边界条件变换到曲线坐标系中,采用一种稳定的、显式的二阶精度的有限差分方法离散(曲线坐标系)HTI介质中的弹性波方程。数值实例充分地展现了在不规则地球表面的影响下裂缝介质中地震波传播的复杂性。合成地震记录和波场快照表明裂缝端点产生的散射波在地表处会受不规则地表地形的作用,再次被散射;同理,地表地形产生的散射波,经过裂缝端点时也会被再次散射,尤其是瑞利面波产生的散射波,因其能量很强,严重污染了地震记录,使得识别地下裂缝等产生的有效信息变得异常困难。这对山地地震勘探中资料的解释具有重要意义。 相似文献
11.
由于构造运动等作用,TI介质对称轴往往沿空间任意方向分布,具有任意空间取向对称轴的TI（ATI）介质更符合实际地质情况.VTI介质与ATI介质的相速度在形式上具有一致性,VTI介质中地震波的相角对应ATI介质对称轴与地震波传播方向的夹角.本文基于Tsvankin的VTI介质精确相速度公式,利用TI介质对称轴和地震波传播方向上单位向量的数量积和向量积来计算ATI介质的精确相速度.根据弱各向异性假设,导出qP波和qSV波的近似相速度,分析了近似公式的误差,讨论总结了ATI介质qP波和qSV波的相速度特征.本文中的单位向量采用观测坐标系表示,通过相角关系,可以较为方便地由ATI介质近似相速度导出频散关系,然后借助傅里叶逆变换推导出时间-波数域qP波和qSV波解耦的波动方程.数值算例表明本文的波动方程是qP波和qSV波解耦的,波场计算结果稳定,未出现明显的数值频散,验证了本文方法的有效性. 相似文献
12.
Luděk Klimeš 《Studia Geophysica et Geodaetica》2016,60(3):391-402
For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we estimate to what extent the medium is transversely isotropic (uniaxial) and determine the direction of its reference symmetry axis expressed in terms of the unit reference symmetry vector. If the medium is exactly transversely isotropic (exactly uniaxial), we obtain the direction of its symmetry axis. We can also calculate the first–order and second–order spatial derivatives of the reference symmetry vector which may be useful in tracing the reference rays for the coupling ray theory. The proposed method is tested using various transversely isotropic (uniaxial) and approximately transversely isotropic (approximately uniaxial) media. 相似文献
13.
假设横向各向同性(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介质对称轴倾角的可行性提供了理论依据和简便的判据. 相似文献
14.
Characterizing the expressions of seismic waves in elastic anisotropic media depends on multiparameters. To reduce the complexity, decomposing the P-mode wave from elastic seismic data is an effective way to describe the considerably accurate kinematics with fewer parameters. The acoustic approximation for transversely isotropic media is widely used to obtain P-mode wave by setting the axial S-wave phase velocity to zero. However, the separated pure P-wave of this approach is coupled with undesired S-wave in anisotropic media called S-wave artefacts. To eliminate the S-wave artefacts in acoustic waves for anisotropic media, we set the vertical S-wave phase velocity as a function related to propagation directions. Then, we derive a pure P-wave equation in transversely isotropic media with a horizontal symmetry axis by introducing the expression of vertical S-wave phase velocity. The differential form of new expression for pure P-wave is reduced to second-order by inserting the expression of S-wave phase velocity as an auxiliary operator. The results of numerical simulation examples by finite difference illustrate the stability and accuracy of the derived pure P-wave equation. 相似文献
15.
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. 相似文献
16.
Luděk Klimeš 《Studia Geophysica et Geodaetica》2018,62(2):255-260
For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we can estimate the extent to which the medium is transversely isotropic, and determine the direction of its reference symmetry axis. In this paper, we rotate the given stiffness tensor about this reference symmetry axis, and determine the reference transversely isotropic (uniaxial) stiffness tensor as the average of the rotated stiffness tensor over all angles of rotation. The obtained reference transversely isotropic (uniaxial) stiffness tensor represents an analytically differentiable approximation of the given generally anisotropic stiffness tensor. The proposed analytic method is compared with a previous numerical method in two numerical examples. 相似文献
17.
Two-dimensional seismic processing is successful in media with little structural and velocity variation in the direction perpendicular to the plane defined by the acquisition direction and the vertical axis. If the subsurface is anisotropic, an additional limitation is that this plane is a plane of symmetry. Kinematic ray propagation can be considered as a two-dimensional process in this type of medium. However, two-dimensional processing in a true-amplitude sense requires out-of-plane amplitude corrections in addition to compensation for in-plane amplitude variation. We provide formulae for the out-of-plane geometrical spreading for P- and S-waves in transversely isotropic and orthorhombic media. These are extensions of well-known isotropic formulae.
For isotropic and transversely isotropic media, the ray propagation is independent of the azimuthal angle. The azimuthal direction is defined with respect to a possibly tilted axis of symmetry. The out-of-plane spreading correction can then be calculated by integrating quantities which describe in-plane kinematics along in-plane rays. If, in addition, the medium varies only along the vertical direction and has a vertical axis of symmetry, no ray tracing need be carried out. All quantities affecting the out-of-plane geometrical spreading can be derived from traveltime information available at the observation surface.
Orthorhombic media possess no rotational symmetry and the out-of-plane geometrical spreading includes parameters which, even in principle, are not invertible from in-plane experiments. The exact and approximate formulae derived for P- and S-waves are nevertheless useful for modelling purposes. 相似文献
For isotropic and transversely isotropic media, the ray propagation is independent of the azimuthal angle. The azimuthal direction is defined with respect to a possibly tilted axis of symmetry. The out-of-plane spreading correction can then be calculated by integrating quantities which describe in-plane kinematics along in-plane rays. If, in addition, the medium varies only along the vertical direction and has a vertical axis of symmetry, no ray tracing need be carried out. All quantities affecting the out-of-plane geometrical spreading can be derived from traveltime information available at the observation surface.
Orthorhombic media possess no rotational symmetry and the out-of-plane geometrical spreading includes parameters which, even in principle, are not invertible from in-plane experiments. The exact and approximate formulae derived for P- and S-waves are nevertheless useful for modelling purposes. 相似文献
18.
The offset‐midpoint travel‐time pyramid for P‐wave in 2D transversely isotropic media with a tilted symmetry axis 
下载免费PDF全文
下载免费PDF全文 For pre‐stack phase‐shift migration in homogeneous isotropic media, the offset‐midpoint travel time is represented by the double‐square‐root equation. The travel time as a function of offset and midpoint resembles the shape of Cheops’ pyramid. This is also valid for transversely isotropic media with a vertical symmetry axis. In this study, we extend the offset‐midpoint travel‐time pyramid to the case of 2D transversely isotropic media with a tilted symmetry axis. The P‐wave analytical travel‐time pyramid is derived under the assumption of weak anelliptical property of the tilted transverse isotropy media. The travel‐time equation for the dip‐constrained transversely isotropic model is obtained from the depth‐domain travel‐time pyramid. The potential applications of the derived offset‐midpoint travel‐time equation include pre‐stack Kirchhoff migration, anisotropic parameter estimation, and travel‐time calculation in transversely isotropic media with a tilted symmetry axis. 相似文献
19.
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. 相似文献
20.