自适应网格细化算法模拟地震波传播

来源 :清华大学 | 被引量 : 0次 | 上传用户:cloudzhu429
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
数值方法求解波动方程一直都是地球物理学界和计算数学界重要的研究课题。数值方法模拟地震波传播现象,可以为天然地震预测和地震波勘探提供理论依据。现有的数值方法使用非常细化的网格,可以得到很好的地震波传播模拟效果,但是采用非常细化的网格计算时,必然会增加计算时间和计算存储量。研究高效的数值方法模拟地震波传播现象,成为求解地震波动方程的关键问题。由于地下介质分布比较复杂,介质之间存在着波速变化,给地震波传播模拟造成了很大困难。高分辨率有限体积法可以得到波在不连续介质中传播的高分辨率的数值解,本文使用高分辨率有限体积法求解声波方程和弹性波方程,构造了波传播的算子,并对波传播算子进行修正限制;数值算例表明:该方法不但可以得到高分辨率的波场快照,而且有效地控制了数值振荡。通过数值解和精确解比较、频散分析,证明了高分辨率有限体积法模拟震波传播的有效性。为了提高地震波传播模拟的计算效率,本文提出了使用自适应网格细化算法结合高分辨率有限体积法求解地震声波和弹性波方程。首先,使用高分辨率有限体积法得到基网格上的数值解,接着,对数值解进行截断误差估计,根据预先设定的阈值,确定需要进行网格细化的区域,将得到的细化网格嵌套在基网格上。自适应网格细化算法自动生成层层嵌套的细化网格,时间步长与空间步长同时细化,采用树状数据结构存储网格点。本文分别使用自适应网格和非自适应网格求解地震声波和弹性波方程,并对计算时间和计算存储量进行比较,结果表明自适应网格细化算法的高效性。提升算法得到的第二代小波变换相对于小波变换运算速度更快,并且保持了多尺度分解的特性。本文利用第二代小波的快速变换和多尺度分解性质,构造了多尺度层层嵌套的自适应网格求解二维声波方程,网格点的空间导数使用有限差分法得到,并给出了声波数值传播的波场快照、计算精度的分析和计算效率的比较。数值实验结果,说明了基于第二代小波的自适应网格求解声波方程的可行性。
其他文献
太原晋阳大佛开凿于北齐年间,是迄今为止我国发现的最早的大型石刻佛像,具有极高的历史和文化价值。大佛在元末毁于战火,在1980年被重新发现时已是四肢残缺、破碎不堪。文物
新课改背景下的教育教学之中,要求各学科的老师都要积极运用人文关怀思想理念,以此来更好的培养学生的学习情感,提升学生的学科素养.作文教学作为初中语文教学中的难点内容,
中国的传统文化由于一些客观因素出现了断层,而中国家具在一味模仿西方与形式化的传统复古中均难寻发展之路。笔者从传统文化的认识误区开始,摒弃了传统文化表面化的低层次元
三塘湖盆地马朗凹陷上石炭统卡拉岗组(C2k)、哈尔加乌组(C2h)火山岩油藏勘探相继获得突破,获得了较大的油气储量规模。但火山岩油藏形成的控制因素较多,对火山岩油藏形成条件、主
数学史是研究数学起源与发展的历史,是理解数学的一种途径,其数学教育价值逐渐受到人们的重视。在小学数学教学中渗透数学史,有利于提高学生学习数学的兴趣,培养学生的创新精
【摘 要】本文分析了GW16型隔离开关检修培训过程中存在的问题,进而介绍GW16型隔离开关技能培训专用检修平台的构成、各主要部件的关键参数及技术要求,最后简单描述了该平台在实际培训中的应用。  【关键词】隔离开关;检修;检修平台  前言  GW16型隔离开关是一种单柱双臂垂直伸缩式隔离开关,其钢支架部分(不含水泥基础)的高度为2130mm-3700mm,再加上槽钢的高度,培训教师或学员在改隔离开关
从企业危机营销管理的内涵、原则入手,探讨企业危机营销管理体系建设。
本文紧密结合苏里格探区苏20区块主要目的层石盒子组油藏精细表征描述与系统地质建模研究工作,以相震约束精细储层表征应用研究为主线,以构造精细解释、地震微相识别、地震分
小组合作学习对于提高初中语文课堂效率,激发学生学习兴趣都有较好的效果.本文强调了初中语文教学中采用小组合作学习的优势,但也说明了现阶段课堂小组合作中容易出现的问题.