共查询到19条相似文献，搜索用时 203 毫秒
1.
2.
利用结构力学与最优控制的模拟理论,研究弹性波在层状介质中传播的数值计算方法. 将弹性波传播问题导向哈密顿(Hamilton)体系,在哈密顿体系中,推导出一种新的半解析单元,称之为动力-部分杂交元,由此导出一套哈密顿体系下的半解析数值计算方法. 本文给出了该方法在层状正交各向异性材料介质的弹性波传播问题的数值算例,分析了一定频率的弹性波在层状介质中传播时的位移、应力的模式. 计算结果展现了Hamilton体系和辛几何在弹性波传播问题研究的应用前景. 相似文献
3.
二维频率空间域25点优化系数差分格式弹性波数值模拟(英文) 总被引：2，自引：0，他引：2
频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播.但是该方法的最大瓶颈是对于计算机内存的巨大需求.我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量.同时为了减少短筹分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数.为了克服边界反射,采用了最佳匹配层吸收边界条件.数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度.使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率.该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础. 相似文献
4.
数值计算方法是考察非线性弹性波在介质中（如岩石）传播特征的重要手段.非线性弹性波的数值模拟存在陡峭间断面（点）、数值振荡以及误差的指数级增长等现象而破坏数值解的稳定性、收敛性，能否消减上述现象的不利影响成为制约数值方法有效与否的重要因素.文中同时引入了FCT算子和幅值限制器，采用中心差分格式对具有垂直对称轴的横向各向同性（VTI）介质中的二维非线性弹性波进行数值模拟，从而克服了上述困难；介绍了适用于非线性弹性波的吸收边界条件，给出了差分方程的稳定性条件.在验证了方法的有效性后成功地获取了二维VTI介质中非线性弹性波的三分量地震正演记录，表明非线性波在传播过程中会发生波形畸变等现象. 相似文献
5.
文章提出一种基于谱元法的区域分解算法(DDM-SEM),并将其用于频域弹性波问题的计算.该方法结合了谱元法的高精度以及区域分解法的高并行特性,适合于大尺度弹性波问题的高精度和高效率计算.在DDMSEM算法中,原始的大尺度问题被分解为若干个子区域,每个子区域采用谱元法独立求解,子区域之间的交界面通过频域黎曼传输条件连接.本文提出的区域分解方法在不同的子区域不仅可采用非匹配网格,而且还可以选取不同的插值阶数,实现计算效率最大化.通过对每个子区域的内部自由度和边界自由度进行分离,结合高效和天然并行性的分块LDU直接求解器能快速求解系统矩阵.数值结果验证了该方法的精度和效率,同时也表明DDMSEM算法在计算大规模、多尺度的弹性波模型中展现出数值方法的优越性.同时,本文提出的方法对于频域下多源地震波反问题的计算具有潜在的应用价值. 相似文献
6.
波动问题有限元离散后会引起数值误差, 数值频散的本质就是数值误差传播引起的非物理解. 数值频散不仅没有实际意义, 而且还会影响对真实波动现象的认识. 为厘清有限元三角网格中波动数值频散的影响因素, 本文推导了集中质量矩阵和一致质量矩阵的频散函数, 同时给出了组合质量矩阵的频散函数, 并对不同质量矩阵的数值频散进行了对比研究. 理论分析和数值计算结果表明: 有限元三角网格中波动的数值频散受网格布局、 波传播方向、 单元网格纵横比以及质量矩阵的影响; 一致质量矩阵的数值频散比集中质量矩阵更易受到波传播方向的影响; 不合理的三角网格单元会对数值相速度(数值频散)产生不良影响; 正三角网格中波动的数值频散几乎不受波传播方向的影响; 一致质量矩阵与集中质量矩阵的线性组合能够有效地压制数值频散. 相似文献
7.
相速度和群速度是研究地震波传播规律和描述介质特性的重要参数,是弹性波传播理论中的核心内容,在理论研究和实际应用中有重要作用.本文根据VTI介质的刚度矩阵,利用Bond变换建立了TTI介质刚度矩阵.再利用TTI介质刚度矩阵,结合弹性动力学的本构方程、牛顿运动微分方程和几何方程,得到了三维TTI介质弹性波波动方程和Christoffel方程.通过本征值方法求解Christoffel方程,推导了三维TTI介质弹性波相速度的解析表达式.利用Berryman和Crampin推导各向异性介质群速度公式,根据三维TTI介质的相速度解析式推导了三维TTI介质群速度解析表达式.数值试例表明,随着各向异性介质参数改变,TI介质弹性波相速度变化较为平缓,群速度变化较为剧烈,qP波和SH波速度变化较为平缓,qSV波速度变化较为剧烈. 相似文献
8.
当地下介质存在各向异性时,在观测坐标系下的弹性参数与自然坐标系下的弹性参数不一定相同．首先,根据势能密度和耗散能密度与坐标轴无关的原理,推导出了双相各向异性介质中观测坐标系下弹性参数与自然坐标系下弹性参数之间的关系;然后,从任意双相各向异性中弹性波波动方程出发,得出了该方程的伪谱法数值解法;最后,通过数值模拟,观测到了存在于双相各向异性介质中的４类波,即快纵波、慢纵波、快横波和慢横波．在双相各向异性介质中,ＳＶ波传播的波前面上仍然存在波面尖角,这些尖角在界面上要发生反射和透射．另外,数值模拟结果中可见转换慢纵波和慢纵波的转换波． 相似文献
9.
10.
孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen（LAC）广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段. 相似文献
11.
准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础. 在二维加权近似解析离散化(WNAD)算法的基础上, 本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法. 我们首先介绍三维WNAD方法的构造过程, 然后与经典的差分格式--交错网格(SG)算法进行了比较. 理论分析和数值算例表明, WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟. 同时, 本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律, 观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等. 在TTI介质分界面处, 弹性波产生更加复杂的折射、反射和波型转化, 使得波场非常复杂, 研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识. 相似文献
12.
常见弹性波动理论的建立是基于介质均匀这一基本假设,实际介质的非均匀性非常普遍.为研究连续介质中波的传播特征,本文从弹性力学中建立弹性波动方程的三个基本方程出发,考虑连续介质弹性参数的空变特征,建立非均匀介质的弹性波动方程,利用Alkhalifah声学近似思想建立位移表征的纵波波动方程,利用本征值问题求解方法建立标量波频率-波数域传播算子,从而建立描述纵波传播的标量波方程,其中波函数为纵波位移的散度,不同于均匀介质标量波方程的波函数为位移势.随后推导含PML边界波动方程差分格式并建立不同模型数值模拟进行数值试算,与均匀假设标量波方程和变密度方程对比证明本方法的准确性和稳定性. 相似文献
13.
Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations 总被引：2，自引：0，他引：2
The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second-order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite-element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media. 相似文献
14.
15.
Myung Hoon Kim Yunseok Choi Young Ho Cha Changsoo Shin 《Pure and Applied Geophysics》2009,166(12):1967-1985
In order to correctly interpret marine exploration data, which contain many elastic signals such as S waves, surface waves and converted waves, we have developed both a frequency-domain modeling algorithm for acoustic-elastic coupled media with an irregular interface, and the corresponding waveform inversion algorithm. By applying the continuity condition between acoustic (fluid) and elastic (solid) media, wave propagation can be properly simulated throughout the coupled domain. The arbitrary interface is represented by tessellating square and triangular finite elements. Although the resulting complex impedance matrix generated by finite element methods for the acoustic-elastic coupled wave equation is asymmetric, we can exploit the usual back-propagation algorithm used in the frequency domain through modern sparse matrix technology. By running numerical experiments on a synthetic model, we demonstrate that our inversion algorithm can successfully recover P- and S-wave velocity and density models from marine exploration data (pressure data only). 相似文献
16.
Terumi Touhei 《地震工程与结构动力学》1995,24(7):937-949
Generally, waves in layered media include several surface wave modes showing dispersive properties. These properties often complicate the understanding of the propagation characteristics of elastic waves. In this paper, not only a method for decomposing transient elastic waves into normal modes is presented but also investigations are applied to decomposed normal modes. This method is developed by means of DWFE and the modal superposition method. In the formulation process, transfer functions describing the relationship between the amplitude of normal modes and that of force density are defined. It was found that the propagation properties of the normal modes were well explained by the dispersion curves and transfer functions. From investigations of the normal mode propagation, properties of each normal mode can be discussed. 相似文献
17.
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
ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｅｄｉａｉｎｅａｒｔｈａｒｅｑｕｉｔｅｃｏｍｐｌｅｘ．Ｔｈｅｒｅｅｘｉｓｔｓｅｖｅｒａｌｕｎｃｏｎｔｉｎｕｏｕｓｐｌａｉｎｓ．Ｎｏｒｍａｌｙ,ｉｔｉｓｕｓｅｄｔｏａｐｐｒｏｘｉｍａｔｅｒｅａｌｍｅｄｉｕｍｗｉｔｈｌａｙ... 相似文献
18.
Modelling the wave phenomena in acoustic and elastic media with sharp variations of physical properties using the grid‐characteristic method 
下载免费PDF全文
下载免费PDF全文 Alena V. Favorskaya Michael S. Zhdanov Nikolay I. Khokhlov Igor B. Petrov 《Geophysical Prospecting》2018,66(8):1485-1502
This paper introduces a novel method of modelling acoustic and elastic wave propagation in inhomogeneous media with sharp variations of physical properties based on the recently developed grid‐characteristic method which considers different types of waves generated in inhomogeneous linear‐elastic media (e.g., longitudinal, transverse, Stoneley, Rayleigh, scattered PP‐, SS‐waves, and converted PS‐ and SP‐waves). In the framework of this method, the problem of solving acoustic or elastic wave equations is reduced to the interpolation of the solutions, determined at earlier time, thus avoiding a direct solution of the large systems of linear equations required by the FD or FE methods. We apply the grid‐characteristic method to compare wave phenomena computed using the acoustic and elastic wave equations in geological medium containing a hydrocarbon reservoir or a fracture zone. The results of this study demonstrate that the developed algorithm can be used as an effective technique for modelling wave phenomena in the models containing hydrocarbon reservoir and/or the fracture zones, which are important targets of seismic exploration. 相似文献
19.
模拟弹性波在孔隙介质中传播，需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上，首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法，并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先，本文通过均匀孔隙模型的数值解与解析解的对比，验证所提出的数值方法的正确性；然后，本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能，将完全匹配层与廖氏和阻尼吸收边界进行了对比，研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明，在孔隙介质中，完全匹配层作为吸收边界能十分有效地吸收衰减外行波，无论对体波还是面波，是一种高效边界吸收算法. 相似文献