共查询到17条相似文献，搜索用时 125 毫秒
1.
2.
本文主要研究了黏弹性HTI和EDA介质中地震波的波动参数(包括相速度、慢度、偏振向量和群速度),并基于摄法,推导了P、SV、SH波波动参数的弱各向异性近似公式.文章提出了慢度向量的三种定义形式,分析对比了各种定义方法在求解christo-ffel方程时的具体方法,指出特殊分量法为各向异性黏弹性介质提供了一种研究均匀和非均匀波的更简单、使用更普遍的方法.基于特殊分量法,通过求解christoffel方程,推导出黏弹性HTI介质中均匀、非均匀波的精确相速度、慢度和群速度计算公式,并通过模型计算研究了SH波的相速度特征及其随相角和不均匀参数D的变化规律,结果表明参数D对地震波的相速度大小有一定影响,但对其方位特性无影响,在EDA介质中相速度随方位角变化的规律仍然可指示介质的对称轴方向和裂隙的走向.基于摄动法,以弹性EDA介质为背景介质,通过模型计算对均匀SH波的近似公式的正确性和精度进行验证,结果证明其最大相对误差为1.15%. 相似文献
3.
对Christoffel公式进行Bond变换得到EDA介质的Christoffel方程,并由其非零解推导出EDA介质中视横波(qSV)、横波(SH)、视纵波(qP)的相速度、群速度、偏振向量(质点的振动方向)的三维计算公式.通过模型计算分析了具有水平对称轴的各向异性(HTI)介质和EDA介质中介质对称轴的极角和方位角对相速度、群速度及偏振向量的影响,对其随极角、方位角的变化特征进行了分析,并采用Matlab进行了数值计算,对其特征采用三维显示.通过取极角或方位角为零简化得到HTI介质和具有垂直对称轴的各向异性(VTI)介质中地震波的相速度、群速度,对EDA介质中的三维计算结果进行退化验证. 通过数值计算进一步验证了地震波相速度与EDA介质对称轴的相互关系. 结果表明,通过广角地震勘探可探明地下介质的裂隙走向及密度,从而确定灾害体产状. 相似文献
4.
相速度和偏振方向是研究地震波传播规律和描述介质特性的重要参数,在理论研究和实际应用中有重要作用.本文假定倾斜横向各向同性(TTI)介质对称轴位于观测坐标系XOZ面内,在此观测坐标系下直接推导了TTI介质弹性波相速度和偏振方向的解析表达式,再进一步利用Thomsen弱各向异性理论,推导了弱各向异性近似条件下弹性波相速度以及qP波和qSV波偏振方向表达式.理论分析和数值试例表明,在相速度方面,随着各向异性介质参数改变,qP波和qSH波速度变化较为平缓,qSV波速度变化较为剧烈.弹性波相速度近似式误差均较小,能较好地近似精确相速度.在偏振方向方面,SH波偏振方向只是传播方向和对称轴倾角的函数,而与各向异性参数无关,SH波偏振方向既垂直于传播方向,又垂直于TTI介质对称轴方向.除特定方向外,qP波和qSV波的偏振方向与传播方向均成一定角度,并且随TTI介质对称轴倾角的改变而改变;在精确和近似情况下,qP波和qSV波的偏振方向始终垂直;在精度允许范围内,偏振方向的弱各向异性近似式与理论解析式吻合较好. 相似文献
5.
由于构造运动等作用,TI介质对称轴往往沿空间任意方向分布,具有任意空间取向对称轴的TI（ATI）介质更符合实际地质情况.VTI介质与ATI介质的相速度在形式上具有一致性,VTI介质中地震波的相角对应ATI介质对称轴与地震波传播方向的夹角.本文基于Tsvankin的VTI介质精确相速度公式,利用TI介质对称轴和地震波传播方向上单位向量的数量积和向量积来计算ATI介质的精确相速度.根据弱各向异性假设,导出qP波和qSV波的近似相速度,分析了近似公式的误差,讨论总结了ATI介质qP波和qSV波的相速度特征.本文中的单位向量采用观测坐标系表示,通过相角关系,可以较为方便地由ATI介质近似相速度导出频散关系,然后借助傅里叶逆变换推导出时间-波数域qP波和qSV波解耦的波动方程.数值算例表明本文的波动方程是qP波和qSV波解耦的,波场计算结果稳定,未出现明显的数值频散,验证了本文方法的有效性. 相似文献
6.
相速度和群速度是研究地震波传播规律和描述介质特性的重要参数,是弹性波传播理论中的核心内容,在理论研究和实际应用中有重要作用.本文根据VTI介质的刚度矩阵,利用Bond变换建立了TTI介质刚度矩阵.再利用TTI介质刚度矩阵,结合弹性动力学的本构方程、牛顿运动微分方程和几何方程,得到了三维TTI介质弹性波波动方程和Christoffel方程.通过本征值方法求解Christoffel方程,推导了三维TTI介质弹性波相速度的解析表达式.利用Berryman和Crampin推导各向异性介质群速度公式,根据三维TTI介质的相速度解析式推导了三维TTI介质群速度解析表达式.数值试例表明,随着各向异性介质参数改变,TI介质弹性波相速度变化较为平缓,群速度变化较为剧烈,qP波和SH波速度变化较为平缓,qSV波速度变化较为剧烈. 相似文献
7.
8.
实际地层中地震波传播普遍存在速度和衰减各向异性现象,研究黏弹各向异性介质中高频地震波传播理论有助于揭示地震波的传播特征.本文针对黏弹性VTI介质,从Christoffel矩阵的解析特征值出发推导出qP、qSV和qSH波的复相速度和复射线速度的解析表达式,并应用实射线追踪方法确定出均匀复射线速度矢量,由此计算出实射线速度和实射线衰减以及实射线品质因子.基于非均匀复相速度和均匀复射线速度的解析表达式,推导了实射线慢度和实射线衰减关于黏弹性模量(包括弹性模量和Q值)的敏感度核函数,该敏感度核函数反映各个黏弹性模量对地震波复走时的影响程度.不同岩石样本的数值计算结果显示,实走时对弹性模量更为敏感,而射线衰减(虚走时)对弹性模量和Q值的敏感程度相当.本研究可为黏弹性VTI介质中地震射线追踪和复走时层析成像提供理论基础. 相似文献
9.
煤层中存在的裂隙会导致介质表现为各向异性,本文以HTI型煤层为例,结合各向异性介质弹性矩阵和各向异性裂隙理论,推导出不同充填物的垂直裂隙中各向异性参数表达式,将其应用于地震波响应分析;通过改进的交错网格差分法和各向异性Christoffel方程波场分解法,得到地震波合成记录和分解后的P波和SV波记录;将Thomsen群速度与相速度公式,经过坐标轴旋转变换,得到HTI型煤层中不同各向异性参数的地震波速度响应表达式;建立不同类型煤层地质模型,分析了裂隙密度、裂隙充填物以及煤层厚度等参数变化时的地震波响应特征.研究结果为分析垂向裂隙各向异性薄煤层地震波传播规律提供工具,为选用相应地震数据进行地震波各向异性参数反演提供依据. 相似文献
10.
Kelvin-Voigt均匀黏弹性介质中传播的地震波 总被引：3，自引：0，他引：3
研究了在小扰动的情况下, 地震波在Kelvin-Voigt均匀黏弹性介质中传播时的瞬态响应以及对应的三阶偏微分方程的解析解. 获得了在脉冲震源条件下Kelvin-Voigt均匀黏弹性三阶波动方程的平面波解. 运用脉冲叠加原理, 获得了在任意震源条件下Kelvin-Voigt均匀黏弹性三阶波动方程的平面波解. 讨论了地震波在Kelvin-Voigt均匀黏弹性介质中的传播速度和衰减, 推导出了地震波的衰减系数和传播速度与介质的密度、弹性模量和黏滞系数之间的精确关系. 这些结果能够用于地震勘探中黏弹性岩性参数的反演. 相似文献
11.
12.
13.
Anisotropy in subsurface geological models is primarily caused by two factors: sedimentation in shale/sand layers and fractures. The sedimentation factor is mainly modelled by vertical transverse isotropy (VTI), whereas the fractures are modelled by a horizontal transversely isotropic medium (HTI). In this paper we study hyperbolic and non‐hyperbolic normal reflection moveout for a package of HTI/VTI layers, considering arbitrary azimuthal orientation of the symmetry axis at each HTI layer. We consider a local 1D medium, whose properties change vertically, with flat interfaces between the layers. In this case, the horizontal slowness is preserved; thus, the azimuth of the phase velocity is the same for all layers of the package. In general, however, the azimuth of the ray velocity differs from the azimuth of the phase velocity. The ray azimuth depends on the layer properties and may be different for each layer. In this case, the use of the Dix equation requires projection of the moveout velocity of each layer on the phase plane. We derive an accurate equation for hyperbolic and high‐order terms of the normal moveout, relating the traveltime to the surface offset, or alternatively, to the subsurface reflection angle. We relate the azimuth of the surface offset to its magnitude (or to the reflection angle), considering short and long offsets. We compare the derived approximations with analytical ray tracing. 相似文献
14.
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. 相似文献
15.
Average elastic properties of a fluid‐saturated fractured rock are discussed in association with the extremely slow and dispersive Krauklis wave propagation within individual fractures. The presence of the Krauklis wave increases P‐wave velocity dispersion and attenuation with decreasing frequency. Different laws (exponential, power, fractal, and gamma laws) of distribution of the fracture length within the rock show more velocity dispersion and attenuation of the P‐wave for greater fracture density, particularly at low seismic frequencies. The results exhibit a remarkable difference in the P‐wave reflection coefficient for frequency and angular dependency from the fractured layer in comparison with the homogeneous layer. The biggest variation in behaviour of the reflection coefficient versus incident angle is observed at low seismic frequencies. The proposed approach and results of calculations allow an interpretation of abnormal velocity dispersion, high attenuation, and special behaviour of reflection coefficients versus frequency and angle of incidence as the indicators of fractures. 相似文献
16.
对于倾斜叠层和非垂直裂隙岩层,用具有任意空间取向对称轴的TI (ATI)模型来描述更符合实际观测. 本文基于坐标变换的方法,研究任意强弱ATI介质中体波速度的角散和方位变化特征. 研究结果表明,ATI介质中体波速度随传播方向变化的速度图案相对TI对称轴确定,此确定的速度图案与TI的Thomsen参数相关;速度特征只依赖于传播矢量与对称轴的夹角. 因此,随着TI对称轴取向的空间变化和测线方位的变化,体波速度图案呈现多样性变化,并具有一定的对称性、渐变和重复性. 研究结果有助于进一步的理论研究和各向异性资料处理解释. 相似文献
17.
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. 相似文献