弹性波方程的数值模拟

来源 :同济大学理学院 同济大学 | 被引量 : 0次 | 上传用户:qq77880066
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
地震波传播数值模拟是研究地震波传播规律的重要手段,在石油和天然气勘探研究中受到越来越多的关注。它可以为地震数据采集、数据处理和解释提供理论指导。 地震波数值模拟中的一个重要问题是如何控制工作量的问题。有限元和有限元差分相结合的方法虽然可以大量的缩短计算时间,但是长时间的计算会出现非物理现象。本文提出有限元和有限元差分形式计算相结合的方法不仅保留了上述方法的优点,而且长时间计算没有产生非物理现象。 在对PML边界条件的方程进行有限元离散时,我们采用两种方法处理PML方程的边界积分项。一种方法是省略边界积分项。另一种方法是保留边界的积分项,本文将对这两种方法分别进行了讨论。 如何处理自由边界的地表条件是地震波数值模拟中的一个重要问题。在实际中,研究单位根据给定计算区域的实际介质传播数值划定地表形状,这样得到的是锯齿状的地表形状,会引发数据记录的信噪比问题。本文使用样条插值逼近地表形状,这样大大减少了人为的波场散射现象。 此外,本文还讨论了计算区域的角点问题。 最后,结合新疆地区实际模型的数值例子,我们使用本文采用的地表处理和控制工作量的方法进行模拟,并用波场快照和时间序列图的形式将数值结果表示出来。结果表明本文讨论的方法是行之有效的。
其他文献
本文讨论反演下列抛物型偏微分方程定解问题{ut=a2△u-q(x)u(x,t)∈Q×(0,T)u(x,t)=0(x,t)∈Ω×[0,T]u(x,0)=u0(x)x∈Ω未知系数q(x)的数值优化方法,其中初始条件u0(x)是已知的。所
在工程项目建筑的过程中因为业主的需求以及设计方面的原因等等经常会需要对于工程项目进行一定的变更,在工程变更的过程中存在着较多的风险问题,一旦控制不好就会影响整个工
许多产品在制作之前要进行几何曲线与曲面的设计,例如汽车车壳设计、船体设计与放样、飞机机身、机舱及机翼的设计,甚至刀片、刀架、鞋模等的外型设计以及服装设计,都对曲线与曲
基于一般自回归条件异方差原理的金融变量回报估计方法,即GARCH类模型已成为描述金融市场波动性的有效工具。而GARCH模型参数估计的主要方法BHHH算法和广义矩方法在实际运算
Voronoi图是一个关于空间划分的基础数据结构。作为计算几何的一个重要分支,它的应用越来越为人们所关注。100年来,它被应用在与几何信息相关的各个领域。随着计算机技术的普及