共查询到18条相似文献，搜索用时 126 毫秒
1.
2.
3.
自然界的裂隙排列并非都是直立的，在复杂地壳应力作用下裂隙排列很可能是倾斜的.本文从各向异性弹性动力学基本方程出发，在特征矩阵方法基础上提出了计算任意倾斜裂隙介质水平界面上的反射、透射系数的简单方法，并利用该方法计算了倾斜裂隙介质P-P波、P-SV波和P-SH波的反透射系数.结果表明，P-P波、P-SV波和P-SH波的反透射系数与倾斜裂隙的倾角、测线方位以及偏移距这些物理参数之间存在一定的依赖关系.此外我们计算了C波的零入射相角时的反射系数，发现该反射系数跟裂隙倾角、裂隙的方位角有着明显的对应规律，且各向异性越强时，该反射系数越大.这些研究对于深入了解倾斜裂隙的动力学特征具有一定的意义. 相似文献
4.
以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征. 相似文献
5.
基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波. 相似文献
6.
7.
准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础. 在二维加权近似解析离散化(WNAD)算法的基础上, 本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法. 我们首先介绍三维WNAD方法的构造过程, 然后与经典的差分格式--交错网格(SG)算法进行了比较. 理论分析和数值算例表明, WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟. 同时, 本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律, 观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等. 在TTI介质分界面处, 弹性波产生更加复杂的折射、反射和波型转化, 使得波场非常复杂, 研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识. 相似文献
8.
本文基于广义线性黏弹性体，通过引入黏弹性拉梅常量，提出拉梅差异矩阵，给出了非均匀各向同性介质的黏弹性波简化方程.拉梅差异矩阵具有和物性矩阵相似的形式，与弹性、黏弹性物性矩阵之间具有特定的数量关系，体现了黏弹性与弹性物性参数之间的差异程度，可以大大简化方程.并利用该简化方程对SEG/EAGE二维盐丘模型进行了正演数值模拟，将模拟得到的单炮记录和声波、Carcione黏弹性波方程模拟得到的单炮记录进行了对比分析.结果表明该简化方程模拟得到的单炮记录不仅具有较好的效果，而且比Carcione黏弹性波方程计算效率高，与理论分析基本相符. 相似文献
9.
10.
11.
12.
—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. 相似文献
13.
地震波在各向异性介质中以一个准P波(qP)和两个准S波(qS1和qS2)的形式传播.研究三种波的相速度、群速度以及偏振方向等传播性质能够为各向异性介质中的正反演问题提供有效支撑.具有比横向各向同性(TI)介质更一般对称性的正交各向异性介质通常需要9个独立参数对其进行描述,这使得对传播特征的计算更为复杂.当两个准S波速度相近时具有耦合性,从而令慢度的计算产生奇异性.因此,奇异点(慢度面的鞍点和交叉点)附近的反射与透射(R/T)系数的求解不稳定,会导致波场振幅不准确.本文首次通过结合耦合S波射线理论和基于迭代的各向异性相速度与偏振矢量的高阶近似解,得到了适用于正交各向异性介质以qP波入射所产生的二阶R/T系数的计算方法.与基于一阶近似的结果相比,基于二阶近似的方法提高了qP波R/T系数的精度,能得到一阶耦合近似无法表达的准确的qP-qS转换波的R/T系数解,且方法适用于较强的各向异性介质. 相似文献
14.
A problem of reflection and transmission of elastic waves at a plane interface between a uniform elastic solid half-space and a porous elastic half-space containing two immiscible fluids is investigated. The theory developed by Lo, Sposito and Majer for porous media containing two immiscible fluids is employed to find out the reflection and transmission coefficients. The incident wave is assumed to propagate through the uniform elastic half-space and two cases are considered. In the first case, a beam of plane longitudinal wave is assumed to be incident and in the second case, a beam of transverse wave is assumed to be incident at the interface. By taking granite as impervious elastic medium and columbia fine sandy loam containing air-water mixture as porous medium, reflection and transmission coefficients are obtained. By neglecting the inertial coupling coefficients, these coefficients are reduced to those obtained by Tomar and Arora using the theory of Tuncay and Corapcioglu. It is found that the inertial coupling parameters significantly affect the phase speeds and the amplitude ratios of the transmitted waves. 相似文献
15.
16.
Summary An account of possible anomalous effects in reflection and refraction of elastic waves at an interface between anisotropic media is presented. These effects are due to anisotropy and they cannot occur at an interface between isotropic media. The shape of the slowness surface (its local deviations from spherical symmetry) is the decisive factor for appearance of these effects. A numerical example of such anomalous behaviour of elastic waves at a free boundary of the crystal of spinel is presented. 相似文献
17.
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. 相似文献
18.
An account of possible anomalous effects in reflection and refraction of elastic waves at an interface between anisotropic media is presented. These effects are due to anisotropy and they cannot occur at an interface between isotropic media. The shape of the slowness surface (its local deviations from spherical symmetry) is the decisive factor for appearance of these effects. A numerical example of such anomalous behaviour of elastic waves at a free boundary of the crystal of spinel is presented. 相似文献