首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
任意空间取向TI弹性张量解析表述   总被引:6,自引:4,他引:2       下载免费PDF全文
姚陈  蔡明刚 《地球物理学报》2009,52(9):2345-2348
本文理论给出任意空间取向TI(ATI)四阶弹性张量的解析表述,其以VTI弹性常数及其简单组合为系数,包括各向同性项、TI对称轴方向矢量分量的二次项和四次项,其中TI对称轴方向矢量可以在固定坐标系定义, 也可以相对三维倾斜界面甚至相对波传播方向.相比四阶张量变换法和Bond变换法,ATI弹性张量能简洁而透明地为本构关系和波动方程提供四阶张量的所有元素. ATI弹性张量能为诸多方面的理论研究提供支撑.  相似文献   

2.
介质的弹性常数为三维四阶张量的分量,共有81个,由于应力张量和应变张量的对称性及能量密度是应变的二次函数,一般各向异常性介质的独立弹性常数可减为21个,如果介质具有较高的对称性,独立弹性常数的数目会更少。 对于地壳和上地幔,具有5个独立弹性常数的横向各向同性介质是一个非常好的近似,本研究中横向各向同性介质的对称轴方向可以是任意的(即对称轴可以不平等于铅直方向),在此情况下,需要进行坐标变换,如果已知介质在某一坐标系(其坐标轴平行或垂直于介质的对称轴)中的弹性常数,我们能够容易地利用变换公式得到变换后新坐标系中的弹性常数。 本文提出了一种方案,利用伪谱法既能模拟横向各向同性介质中的平面波,也能模拟点源激发的波场。在勘探地球物理和地震学中,模拟横向各向同性介拮中传播的平面波及区域源产生的波是最重要的研究课题之一。然而在一般各向异性介质中,很难或不可能确定弹性波的相速度和偏振方向,但在横向各向同性介质中,则可以通过坐标变换来实现,这里我们所提出的方法可以用于横向各向同性介质中弹性波的模拟。  相似文献   

3.
Although synthetic seismograms can, in many cases, be generated by direct numerical methods, approximate analytical ray-theory solutions are often very useful in forward and especially inverse problems of seismic wave propagation in complex 3-D media. In this paper, the frequency-domain zero-order ray-theory Green tensor in a heterogeneous anisotropic elastic medium is derived from the zero-order ray-theory approximation using the representation theorem applied in ray-centred coordinates.  相似文献   

4.
Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement.  相似文献   

5.
A general definition of seismic wave impedance is proposed as a matrix differential operator transforming the displacement boundary conditions into traction ones. This impedance is proportional to the standard acoustic impedance at all incidence angles and allows extensions to attenuative media and to the full elastic case. In all cases, reflection amplitudes at the contact of two media are uniquely described by the ratios of their impedances. Here, the anelastic acoustic impedance is studied in detail and attenuation contrasts are shown to produce phase‐shifted reflections. Notably, the correspondence principle (i.e., the approach based on complex‐valued elastic modules in the frequency domain) leads to incorrect phase shifts of the impedance due to attenuation and consequently to wrong waveforms reflected from attenuation contrasts. Boundary conditions and the Lagrange formulation of elastodynamics suggest that elastic constants should remain real in the presence of attenuation and the various types of energy dissipation should be described by their specific mechanisms. The correspondence principle and complex‐valued elastic moduli appear to be applicable only to homogeneous media and therefore they should be used with caution when applied to heterogeneous cases.  相似文献   

6.
Equations are presented to determine the phase shift of the amplitude of the elastic Green tensor due to both simple (line) and point caustics in anisotropic media. The phase-shift rules for the Green tensor are expressed in terms of the paraxial-ray matrices calculated by dynamic ray tracing. The phase-shift rules are derived both for 2×2 paraxial-ray matrices in ray-centred coordinates and for 3×3 paraxial-ray matrices in general coordinates. The reciprocity of the phase shift of the Green tensor is demonstrated. Then a simple example is given to illustrate the positive and negative phase shifts in anisotropic media, and also to illustrate the reciprocity of the phase shift of the Green tensor.  相似文献   

7.
本文应用交错网格高阶有限差分方法模拟弹性波在三维各向同性介质中的传播。采用时间上二阶、空间上高阶近似的交错网格高阶差分公式求解三维弹性波位移-应力方程,并在计算边界处应用基于傍轴近似法得到的三维弹性波方程吸收边界条件。在此基础上进行了三维盐丘地质模型的地震波传播数值模拟试算。试算结果表明该方法模拟精度高,在很大程度上减小了数值频散,绕射波更加丰富,而且适用于介质速度具有纵向变化和横向变化的情况。  相似文献   

8.
基于临界孔隙度模型的地震波传播   总被引:1,自引:1,他引:0       下载免费PDF全文
基于岩石物里学中临界孔隙度模型,建立一种简洁的均匀弹性流体饱和孔隙介质模型,进行地震波传播研究.首先定义了构建目标模型的基本力学模型:介绍了全孔隙度区间内基本力学模型和目标孔隙介质的含义,其中基本力学模型除了完全弹性固体模型S和完全弹性流体模型F还包括临界孔隙模型C.然后通过等效力学模型推出了目标力学模型介质本构关系的组分表达形式.文中分别通过直接求取弹性参数的表达形式和运用应力应变关系两种方法得到介质模型的本构关系,进而得到该模型波动方程的组分表达形式.最后对这种介质模型进行了地震波传播的数值模拟,结合模拟结果分析孔隙对地震波传播的影响.  相似文献   

9.
李磊  郝重涛 《地球物理学报》2011,54(11):2819-2830
Thomsen提出的横向各向同性(TI)介质各向异性参数(α0、β0、ε、δ、γ)是各向异性理论研究和实际资料处理中的常用参数,Thomsen参数的取值必须符合物理学定律和实际地学情况,随意的取值可能导致无意义乃至错误的结果.本文根据热力学定律和弹性常数的物理意义,结合大量的实测数据,提出常见TI介质的Thomsen参...  相似文献   

10.
Although clay is composed of disconnected anisotropic clay platelets, many rock physics models treat the clay platelets in shale as interconnected. However, the clay matrix in shales can be modelled as anisotropic clay platelets embedded within a soft isotropic interplatelet region, allowing the influence of disconnected clay platelets on the elastic properties of the clay matrix to be analysed. In this model, properties of the interplatelet region are governed by its effective bulk and shear moduli, whereas the effective properties of the clay platelets are governed by their volume fraction, aspect ratio and elastic stiffness tensor. Together, these parameters implicitly account for variations in clay and fluid properties, as well as fluid saturation. Elastic stiffnesses of clay platelets are obtained from the literature, including both experimental measurements and first-principles calculations of the full anisotropic (monoclinic or triclinic) elastic stiffness tensors of layered silicates. These published elastic stiffness tensors are used to compile a database of equivalent transverse isotropic elastic stiffness tensors, and other physical properties, for eight common varieties of layered silicates. Clay matrix anisotropy is then investigated by examining the influence of these different elastic stiffnesses, and of varying model parameters, upon the effective transverse isotropic elastic stiffness tensor of the clay matrix. The relationship between the different clay minerals and their associated anisotropy parameters is studied, and their impact on the resulting anisotropy of the clay matrix is analysed.  相似文献   

11.
任意空间取向TI介质中体波速度特征分析   总被引:2,自引:0,他引:2       下载免费PDF全文
郝重涛  姚陈  王迅 《地震地质》2005,27(2):252-259
基于任意空间取向TI介质(简称ATI)中体波速度和偏振解析解,通过模型数值计算给出ATI介质中体波群速度和相速度的变化特征,说明TI空间取向与测线方位对速度的影响。研究表明,体波群速度图案和相速度图案相对TI对称轴固定,随TI对称轴倾角及其相对测线方位角的变化呈现出一定的对称性和重复性;可以针对ATI地区的地质情况,给出体波群速度和相速度变化图案,为进一步的理论研究提供便捷。此结果也可以直接用于VSP(垂直地震剖面)和井间地震资料的分析研究  相似文献   

12.
准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础. 在二维加权近似解析离散化(WNAD)算法的基础上, 本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法. 我们首先介绍三维WNAD方法的构造过程, 然后与经典的差分格式--交错网格(SG)算法进行了比较. 理论分析和数值算例表明, WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟. 同时, 本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律, 观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等. 在TTI介质分界面处, 弹性波产生更加复杂的折射、反射和波型转化, 使得波场非常复杂, 研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识.  相似文献   

13.
This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, temporal-and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial resuits for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.  相似文献   

14.
忽略TTI介质对称轴倾角的可行性   总被引:1,自引:1,他引:0       下载免费PDF全文
李磊  郝重涛 《地球物理学报》2012,55(6):2004-2013
假设横向各向同性(TI)介质的对称轴是垂直的(VTI)或者水平的(HTI)能给实际资料处理带来便利,然而实际TI介质的对称轴往往是倾斜的(TTI),忽略对称轴倾角可能给各向异性参数提取和成像带来偏差,因此需要研究是否能、以及什么条件下能忽略TTI介质对称轴倾角.本文通过理论研究和数值分析研究了与TTI介质弹性性质最接近的VTI介质(OAVTI)的弹性常数和各向异性参数与原TTI介质的弹性常数和各向异性参数之间的联系与差别.结果表明:OAVTI介质各向异性参数与原TTI介质各向异性参数之间的差别可统一表示成F(α00,ε,δ,γ)ξ2的形式,其中F(α00,ε,δ,γ)是无量纲各向异性参数(ε, δ, γ)的线性函数,ξ是对称轴倾角;ξ的大小对各参数的误差起主导作用,一般不建议忽略20°~25°以上的对称轴倾角;当ξ较小时,即使是对强各向异性的TTI介质作VTI近似,引起的P波各向异性参数误差也很小,因此在纵波资料处理中忽略TTI介质对称轴倾角通常是可行的;即使在小ξ条件下,倾斜对称轴对SV波也有显著影响,因此在转换波资料处理中,不建议忽略TTI介质的对称轴倾角.本文的研究为分析忽略TTI介质对称轴倾角的可行性提供了理论依据和简便的判据.  相似文献   

15.
完全匹配层吸收边界在孔隙介质弹性波模拟中的应用   总被引:20,自引:6,他引:14       下载免费PDF全文
模拟弹性波在孔隙介质中传播,需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上,首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法,并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先,本文通过均匀孔隙模型的数值解与解析解的对比,验证所提出的数值方法的正确性;然后,本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能,将完全匹配层与廖氏和阻尼吸收边界进行了对比,研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明,在孔隙介质中,完全匹配层作为吸收边界能十分有效地吸收衰减外行波,无论对体波还是面波,是一种高效边界吸收算法.  相似文献   

16.
17.
Elastic wave propagation in inhomogeneous anisotropic media   总被引:1,自引:0,他引:1  
IntroductionThemediaineartharequitecomplex.Thereexistseveraluncontinuousplains.Normaly,itisusedtoapproximaterealmediumwithlay...  相似文献   

18.
An equivalent medium model for wave simulation in fractured porous rocks   总被引:3,自引:0,他引:3  
Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane.  相似文献   

19.
孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen(LAC)广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段.  相似文献   

20.
Bos  L.  Gibson  P.  Kotchetov  M.  Slawinski  M. 《Studia Geophysica et Geodaetica》2004,48(1):265-287
The purpose of the present article is to give a precise definition and analysis from first principals of anisotropy, as the term applies to elastic media, taking care to avoid unnecessary assumptions. Two fundamental concepts, material invariance and symmetry group of a material, are defined purely in terms of the stress-strain relation. The implications of material symmetry, or in other words, of anisotropy, for the structure of the stiffness tensor are then investigated. Using the reduced notation of Voigt, these results are presented as the well-known simplifications in the form taken by the six-by-six stiffness matrix that represents the material's stiffness tensor. A new, simple proof is given for the remarkable fact that an elastic medium cannot have rotational symmetry by an angle of less than 90° without being transversely isotropic. In addition, the mutual relation that the notions of elastic symmetry and crystal symmetry have with respect to the so-called orthogonal group is sketched. Despite the historical association between anisotropic elastic materials and the study of crystals, the given presentation shows that conceptually the notion of anisotropy in elastic media is entirely independent of that of crystal symmetry.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号