规则网格有限差分解声波方程个人总结报告.doc

规则网格有限差分解声波方程个人总结报告.doc

ID:51793202

大小:1.05 MB

页数:19页

时间:2020-03-15

规则网格有限差分解声波方程个人总结报告.doc_第1页
规则网格有限差分解声波方程个人总结报告.doc_第2页
规则网格有限差分解声波方程个人总结报告.doc_第3页
规则网格有限差分解声波方程个人总结报告.doc_第4页
规则网格有限差分解声波方程个人总结报告.doc_第5页
资源描述:

《规则网格有限差分解声波方程个人总结报告.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、地球探测科学与技术学院总结报告学校:吉林大学学院:地球探测科学与技术学院专业:勘查技术与工程(应用地球物理)科目:科学计算方法--有限差分解声波方程姓名:学号:目录一.相关理论基础31.地震波场模拟32.波动方程类型及其局限性33.数值算法类型及其优缺点4二.有限差分解声波方程基础理论知识61.需要的已知条件包括:62.弹性波方程63.声波方程的有限差分法数值模拟64.稳定性条件75.频散关系式86.有限差分参数8三.程序及结果成图8四.通过实验所发现的问题和认识12五.他人所做的有限差分解波动方程

2、程序及结果成图12参考文献及资料19有限差分解声波方程总结报告一.相关理论基础1.地震波场模拟地震波场模拟即地震正演,是指已知模型结构,通过物理或数值计算的方法模拟该地质结构下的地震波的传播,最终合成地震记录,也可以认为其是野外数据采集过程的室内再现。物理模拟花费昂贵,人们一般采用比较经济的数值模拟技术。地震波场数值模拟是在给定数学模型(如弹性波方程,声波方程等)、震源和地下几何界面、物性参数(岩层密度、速度等)情况下,研究弹性波或声波的传播规律。2.波动方程类型及其局限性(1)声波方程:二阶标量声

3、波方程:一阶压力-速度方程组:波动方程能够描述且只能描述纵波的传播规律,包括直达波、反射波、透射波、折射波等,但不能描述转换波传播规律。需要的已知条件包括:震源函数、地层速度、密度边界条件(1)弹性波方程:弹性波方程能够描述纵、横波的传播规律,包括直达波、反射波、透射波、折射波以及转换波等。需要的已知条件包括:震源函数、地层速度或根据方程的类型需要提供的地层的其它弹性参数、边界条件。(2)粘声波/弹性波方程前面讨论的是理想弹性介质,波在其中传播时,没有能量的损耗,介质中应力和应变关系严格遵循胡克定律

4、(这种理想介质称虎克固体),但波在实际介质中传播时,是有能量损耗的,这就是所谓的弹性波吸收。波在传播过程中,实际介质的不同部位之间会出现某种摩擦力,称为内摩擦力或粘滞力。这种力导致机械能向其他形式能量转换,最终转化为热能消耗掉。在地震勘探中,地震波传播的实际介质是十分复杂的。在一定条件下,即震源作用时间短,作用力微小,地球介质可以看作完全弹性模型,但随着地震勘探技术的发展,勘探精度要求提高,面临复杂地质目标时,要求地震勘探采用更加符合实际的介质模型进行研究。粘弹性介质模型更符合实际。但是到目前为止,

5、在地震资料反演处理中应用最多的还是声波方程,弹性波以及粘弹性波方程的应用还只是停留在模拟层次上。2.数值算法类型及其优缺点地震波波动方程数值模拟方法主要包括克希霍夫积分法、傅里叶变换法、有限元法和有限差分法等。(1)克希霍夫积分法①本质是波动方程积分解的一个数值计算,某种程度上相当于绕射叠加。②特点该方法计算速度较快,但由于射线追踪中存在着焦散、多重路径等问题,故其一般只能适合于较简单的模型,难以模拟复杂地层的波场信息。(1)傅里叶变换法①原理利用空间的全部信息对波场函数进行三角函数插值,能更加精确

6、地模拟地震波的传播规律,同时,利用快速傅里叶变换(FFT)进行计算,还可以提高运算效率②特点精度高,占用内存小,但缺点是计算速度较慢,对模型的适用性差,尤其是不适应于速度横向变化剧烈的模型。(2)波动方程有限元法①原理将变分法用于单元分析,得到单元矩阵,然后将单元矩阵总体求和得到总体矩阵,最后求解总体矩阵得到波动方程的数值解②优点理论上可适宜于任意地质体形态的模型,保证复杂地层形态模拟的逼真性,达到很高的计算精度③缺点占用内存和运算量均较大,不适用于大规模模拟,因此该方法在地震波勘探中尚未得到广泛地

7、应用。(3)有限差分法①基本原理该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分方法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的函数值为未知数的代数方程组。②特点该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。③主要内容v建立地球物理问题的离散有限差分模型*根据问题的特点将定解区域做网格划分*在所有网格节点上用有限差分格式对导数求近似,对函数、初始条

8、件和边界条件求近似*把原方程离散化为代数方程组v保证计算过程的可行性和计算结果的正确性:解的相容性、稳定性、收敛性v数值求解差分方程组二.有限差分解声波方程基础理论知识1.需要的已知条件包括:(1)震源函数(2)地层速度(波速)(3)边界条件2.弹性波方程:3.声波方程的有限差分法数值模拟对于二维速度-深度模型,地下介质中地震波的传播规律可以近似地用声波方程描述:(4-1)是介质在点(x,z)处的纵波速度,为描述速度位或者压力的波场。为求式(4-1)的数值解,必须将此

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。