共查询到20条相似文献，搜索用时 592 毫秒
1.
Elastic wave propagation in inhomogeneous anisotropic media 总被引：1，自引：0，他引：1
XIU CHENG WEI ) MIN YU DONG ) YUN YAI CHEN ) ) Geoscience Department University of Petroleum Beijing China ) Institute of Geophysics China Seismological Bureau Beijing China 《地震学报(英文版)》1998,11(6):655-667
ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｅｄｉａｉｎｅａｒｔｈａｒｅｑｕｉｔｅｃｏｍｐｌｅｘ．Ｔｈｅｒｅｅｘｉｓｔｓｅｖｅｒａｌｕｎｃｏｎｔｉｎｕｏｕｓｐｌａｉｎｓ．Ｎｏｒｍａｌｙ,ｉｔｉｓｕｓｅｄｔｏａｐｐｒｏｘｉｍａｔｅｒｅａｌｍｅｄｉｕｍｗｉｔｈｌａｙ... 相似文献
2.
The earth media friction attenuates the energy of seismic wave. Obviously, that earth medium is treated as viscoelasticity is more appropriate than as elasticity. Seismic wave in visco-elastic media contains abundant lithology information of the medium. The investiga-tions of the solutions of visco-elastic wave equations, velocities of seismic wave propagating and the at-tenuation of seismic wave in the visco-elastic media are very important for geophysical prospecting tech-nology. Under a sma… 相似文献
3.
相速度和群速度是研究地震波传播规律和描述介质特性的重要参数,是弹性波传播理论中的核心内容,在理论研究和实际应用中有重要作用.本文根据VTI介质的刚度矩阵,利用Bond变换建立了TTI介质刚度矩阵.再利用TTI介质刚度矩阵,结合弹性动力学的本构方程、牛顿运动微分方程和几何方程,得到了三维TTI介质弹性波波动方程和Christoffel方程.通过本征值方法求解Christoffel方程,推导了三维TTI介质弹性波相速度的解析表达式.利用Berryman和Crampin推导各向异性介质群速度公式,根据三维TTI介质的相速度解析式推导了三维TTI介质群速度解析表达式.数值试例表明,随着各向异性介质参数改变,TI介质弹性波相速度变化较为平缓,群速度变化较为剧烈,qP波和SH波速度变化较为平缓,qSV波速度变化较为剧烈. 相似文献
4.
常见弹性波动理论的建立是基于介质均匀这一基本假设,实际介质的非均匀性非常普遍.为研究连续介质中波的传播特征,本文从弹性力学中建立弹性波动方程的三个基本方程出发,考虑连续介质弹性参数的空变特征,建立非均匀介质的弹性波动方程,利用Alkhalifah声学近似思想建立位移表征的纵波波动方程,利用本征值问题求解方法建立标量波频率-波数域传播算子,从而建立描述纵波传播的标量波方程,其中波函数为纵波位移的散度,不同于均匀介质标量波方程的波函数为位移势.随后推导含PML边界波动方程差分格式并建立不同模型数值模拟进行数值试算,与均匀假设标量波方程和变密度方程对比证明本方法的准确性和稳定性. 相似文献
5.
时间域的波场延拓方法在本质上都可以归结为对一个空间-波数域算子的近似.本文基于一阶波数-空间混合域象征,提出一种新的方法求解解耦的二阶位移弹性波方程.该方法采用交错网格,连续使用两次一阶前向和后向拟微分算子,推导得到了解耦的二阶位移弹性波方程的波场延拓算子.由于该混合域象征在伪谱算子的基础上增加了一个依赖于速度模型的补偿项,可以补偿由于采用二阶中心差分计算时间微分项带来的误差,有效地减少模拟结果的数值频散,提高模拟精度.然而,在非均匀介质中,直接计算该二阶的波场延拓算子,每一个时间步上需要做N次快速傅里叶逆变换,其中N是总的网格点数.为了减少计算量,提出了交错网格低秩分解方法;针对常规有限差分数值频散问题,本文将交错网格低秩方法与有限差分法结合,提出了交错网格低秩有限差分法.数值结果表明,交错网格低秩方法和交错网格低秩有限差分法具有较高的精度,对于复杂介质的地震波数值模拟和偏移成像具有重要的价值. 相似文献
6.
应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法，将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质，给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度。 相似文献
7.
8.
George D. Manolis Tsviatko V. Rangelov 《Soil Dynamics and Earthquake Engineering》2006,26(10):952-959
Geological media are invariably non-homogeneous, which complicates considerably the analysis of seismically induced wave propagation phenomena. Thus, closed-form solutions in the form of Green's functions are difficult to construct, but are quite valuable in their own right and often play the role of kernels in boundary integral equation formulations that are used for the solution of complex boundary-value problems of engineering importance. In this work, we examine in some detail the types of wave-like equations that result from vector decomposition of the equations of motion for the infinitely extending non-homogeneous continuum, which would be a first step for evaluating Green's functions. Specifically, an eigenvalue analysis is first performed, followed by computations using the finite difference method for a specific example involving a soil layer with quadratically varying material parameters. The aforementioned wave-like equations, defined in terms of dilatational and rotational strains, are originally coupled. Their uncoupling involves use of algebraic transformations, which are in turn valid for certain restricted categories of non-homogeneous materials. Numerical solution of these equations clearly shows attenuation patterns and phase changes that are manifested as the incoming wave disturbance is continuously scattered by non-constant material stiffness values encountered along the propagation path. 相似文献
9.
Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method 总被引：4，自引：0，他引：4
Seismic wave propagation in transversely isotropic (TI) media is commonly described by a set of coupled partial differential equations, derived from the acoustic approximation. These equations produce pure P‐wave responses in elliptically anisotropic media but generate undesired shear‐wave components for more general TI anisotropy. Furthermore, these equations suffer from instabilities when the anisotropy parameter ε is less than δ. One solution to both problems is to use pure acoustic anisotropic wave equations, which can produce pure P‐waves without any shear‐wave contaminations in both elliptical and anelliptical TI media. In this paper, we propose a new pure acoustic transversely isotropic wave equation, which can be conveniently solved using the pseudospectral method. Like most other pure acoustic anisotropic wave equations, our equation involves complicated pseudo‐differential operators in space which are difficult to handle using the finite difference method. The advantage of our equation is that all of its model parameters are separable from the spatial differential and pseudo‐differential operators; therefore, the pseudospectral method can be directly applied. We use phase velocity analysis to show that our equation, expressed in a summation form, can be properly truncated to achieve the desired accuracy according to anisotropy strength. This flexibility allows us to save computational time by choosing the right number of summation terms for a given model. We use numerical examples to demonstrate that this new pure acoustic wave equation can produce highly accurate results, completely free from shear‐wave artefacts. This equation can be straightforwardly generalized to tilted TI media. 相似文献
10.
S. S. Kevorkyants 《Izvestiya Physics of the Solid Earth》2012,48(11-12):798-805
The system of Biot vector equations in the frequency space includes two elliptic-type vector partial differential equations with unknown displacement vectors in the solid and liquid phases. Considering the Biot equations, alongside with Pride??s equations, the key approaches to the theoretical study of the elastic waves in the two-phase fluid-saturated media, the author suggests an analytical solution for the inhomogeneous Biot equations in the frequency space, which is reduced to finding its fundamental solution (Green??s function). The solution of this problem consists of solutions for two systems of Biot equations. In the first system, only the first equation is inhomogeneous, while in the second system, only the second equation is inhomogeneous and, as it is shown, its right-hand side is exclusively a potential function. The fundamental solution of the full system of inhomogeneous Biot equations (in which both equations are inhomogeneous) is represented in the form of Green??s matrix-tensor, for the scalar elements of which the analytical relations are presented. The obtained formulas describing the elastic displacements of both the solid and liquid phases reflect three wave types, namely, compressional waves of the first and the second kind (the fast and the slow waves, respectively) and shear waves. Similar terms (those describing the same type of the elastic waves in the solid and liquid phases) in the expressions for Green??s functions are linked with each other through the coefficient that links the components of the displacement vectors of the solid and liquid phases corresponding to the given wave type. 相似文献
11.
The system of Biot vector equations in the frequency space includes two elliptic-type vector partial differential equations with unknown displacement vectors in the solid and liquid phases. Considering the Biot equations, alongside with Pride’s equations, the key approaches to the theoretical study of the elastic waves in the two-phase fluid-saturated media, the author suggests an analytical solution for the inhomogeneous Biot equations in the frequency space, which is reduced to finding its fundamental solution (Green’s function). The solution of this problem consists of solutions for two systems of Biot equations. In the first system, only the first equation is inhomogeneous, while in the second system, only the second equation is inhomogeneous and, as it is shown, its right-hand side is exclusively a potential function. The fundamental solution of the full system of inhomogeneous Biot equations (in which both equations are inhomogeneous) is represented in the form of Green’s matrix-tensor, for the scalar elements of which the analytical relations are presented. The obtained formulas describing the elastic displacements of both the solid and liquid phases reflect three wave types, namely, compressional waves of the first and the second kind (the fast and the slow waves, respectively) and shear waves. Similar terms (those describing the same type of the elastic waves in the solid and liquid phases) in the expressions for Green’s functions are linked with each other through the coefficient that links the components of the displacement vectors of the solid and liquid phases corresponding to the given wave type.
相似文献12.
用偏微分方程组特征理论研究层状并以垂直轴为对称轴的横向各向同性（TIV）介质的参数反演问题，首先从弹性波运动方程与TIV介质的应力－应变关系导出了平面波耦合方程组的特征型．根据奇性分析与特征积分法给出了连续情形下的特征线边界条件，连续清形下的波场延拓方程即平面波耦合方程组的特征型与特征线边界条件组成了层状TIV平滑介质、弱间断介质参数反演问题的基本方程组．并导出了间断情形下的波场延拓方程与特征线边界条件，这些方程组可用于层状介质（间断情况）的参数反演．基于这些基本方程组，探讨了利用地面多分量地震资料反演层状TIV介质多个弹性参数的问题． 相似文献
13.
The strong coupling of applied stress and pore fluid pressure, known as poroelasticity, is relevant to a number of applied problems arising in hydrogeology and reservoir engineering. The standard theory of poroelastic behavior in a homogeneous, isotropic, elastic porous medium saturated by a viscous, compressible fluid is due to Biot, who derived a pair of coupled partial differential equations that accurately predict the existence of two independent dilatational (compressional) wave motions, corresponding to in-phase and out-of-phase displacements of the solid and fluid phases, respectively. The Biot equations can be decoupled exactly after Fourier transformation to the frequency domain, but the resulting pair of Helmholtz equations cannot be converted to partial differential equations in the time domain and, therefore, closed-form analytical solutions of these equations in space and time variables cannot be obtained. In this paper we show that the decoupled Helmholtz equations can in fact be transformed to two independent partial differential equations in the time domain if the wave excitation frequency is very small as compared to a critical frequency equal to the kinematic viscosity of the pore fluid divided by the permeability of the porous medium. The partial differential equations found are a propagating wave equation and a dissipative wave equation, for which closed-form solutions are known under a variety of initial and boundary conditions. Numerical calculations indicate that the magnitude of the critical frequency for representative sedimentary materials containing either water or a nonaqueous phase liquid is in the kHz–MHz range, which is generally above the seismic band of frequencies. Therefore, the two partial differential equations obtained should be accurate for modeling elastic wave phenomena in fluid-saturated porous media under typical low-frequency conditions applicable to hydrogeological problems. 相似文献
14.
15.
从各向同性介质中波场数值模拟的褶积微分算子法出发,推导出了各向异性双相介质中波场传播数值计算的褶积新算法.将常见的二阶微分Biot波动方程用等效的一阶速度—应力双曲方程表示,其中未知的波场向量包括固相和流体的速度分量和应力分量,由此对方程的时间项使用交错网格差分方法计算,而对空间项则采用褶积微分算法进行求解.对各向异性双相介质在单层介质模型和双层介质模型中的波场特征进行了研究.研究的结果显示,在两层介质分界面上当地震波产生反射时能观测到两类纵波和横波,并且在衰减系数大的介质里慢纵波很难见到. 相似文献
16.
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. 相似文献
17.
18.
采用规则网格有限差分方法对二维平面弹性波动方程进行差分离散,得到相应的弹性波动方程的有限差分方程,再将弹性波动方程的差分格式与吸收边界、自由边界的离散形式结合形成弹性波动方程有限差分方程解决问题的主体,将其应用于含方形凹陷半无限非均匀介质的模型中进行数值模拟,得到此离散化模型中不同时刻不同节点的位移值。针对具体算例,运用上述方法结合科学计算软件MATLAB和结果后处理软件DIFEM ISOLINE PLOTER得到不同时刻的水平方向位移等值线图与接收器测量点处的合成位移记录,讨论非均匀介质、吸收边界、方形凹陷等对波动特性的影响。 相似文献
19.
两相饱和多孔介质的动力响应问题在地震工程领域具有重要的研究意义,由于涉及到固相和液相的动力耦合,使得该问题的求解尤为复杂。本文利用Comsol在求解多场耦合问题上的优点,针对Biot饱和多孔介质u-U耦合形式下的波动方程特征,经过一系列微分算子运算和矩阵变换得到导数形式下的波动方程,基于Comsol Multiphysics提供的广义偏微分方程模式对变形后的波动方程进行求解,并把改进后的无限元边界应用到无限域动力问题的模拟中。通过与饱和多孔介质动力响应的解析解进行对比,验证模型求解技术的可行性和正确性,并在此基础上讨论饱和土地基中空沟隔振效果与饱和土体参数孔隙率、泊松比的关系。通过研究分析,可以为饱和土地基中空沟隔振设计提供一些有价值的参考。 相似文献
20.
IntroductionMoreandmoreevidencesshowthatstratapropertiesareanisotropic.Itisveryimportanttoutilizegeophysicalwaystouncoverthepropertiesofstratainpetroleumsurveyandcoalbedmethaneexplorationaswellasearthquakeforecasting.Becauseofdetectablecharacteristicsofseismicvelocityanisotropy,itissuccessfulinpastyearsthatlayersareshiftedcorrectlytotheoriginalpositioninmigrationprocessionofseismicwaves.Anewwayofdirectinversionofstrataelasticparametersisputforthinthepaper.Itusesmultiple-wavesdatatodiscoverorp… 相似文献