三维波动方程的高精度正演方法研究

来源 :哈尔滨工程大学 | 被引量 : 0次 | 上传用户:victor530505911
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
地下介质中波的传播情况的研究,对于地质勘探具有十分重要的应用价值,为了更好的模拟地下波场,需要不断地发展和完善解决问题的方法.而波动方程正演作为正问题研究的重要分支,在众多领域有着重要的应用.例如,我们可以用弹性动力学中的各种波动方程来描述地震波在地下介质的传播,并且通过波动方程的正演模拟来获取地震资料.对于波动方程正问题而言,它面临着构造的显格式对步长有一定的限制、数值解不适定等困难,并且实际应用中会有相当大的计算量和存储量.因此,波动方程正演的研究具有重要的理论意义及广泛的应用价值.对于二维波动方程的研究为更高维度波动方程的研究奠定了理论基础.  首先,本文介绍正演模拟的主要方法,包括有限差分法、有限元法以及伪谱法.其中,对波动方程的初边值问题以及有限差分法三种常用差分格式进行简单介绍,包括古典显格式、古典隐格式以及VonNeumann格式;对于有限元法,介绍其找近似解的主要六步过程;对于伪谱法,介绍中心差分和时间导数,以及傅里叶变换和空间导数.  其次,以二维波动方程作为数学模型.然后应用显式格式与隐式格式的有限差分法对模型进行离散,构造十一点隐式格式和十五点隐式格式,并且采用不同的方式来处理边界条件,然后验证算法的收敛性.接下来,使用MATLAB对构造的显式差分计算格式进行数值模拟,分析模拟效果,说明第二种和第三种边界条件处理要优于一般的边界条件处理方法.  最后,在二维波动方程的基础上,研究三维波动方程的数学模型.其中对三维波动方程进行离散,构造十五点隐式格式和二十一点隐式格式,并采用不同的方式来处理边界条件,然后验证算法的收敛性.接下来,使用MATLAB对构造的显式差分计算格式进行数值模拟,分析模拟效果,说明在三种模型下,虚拟辅助边界层法处理要优于其它两种方法.
其他文献
本文研究带有工件组的单机继列批在线排序问题,批容量分无限和有限,目标函数为最小化最大完工时间。工件的到达时刻是任意的,每个工件都有各自的安装时间。机器在每一时刻最
切换时滞系统是一类重要的混杂系统。这类系统的连续动态、离散动态和时滞相互作用,使得系统的动态特性极为复杂。切换时滞系统有着广泛的实际背景和理论研究意义,目前对切换时
在本文,我们主要考虑一般的线性算子s在满足‖Sz1‖P或‖Sz1‖p,α有界的条件下,给出s有界或紧的刻画,文章主要包含以下几个部分:  第一章,主要介绍了文章的背景,及在整个复平面
随着中国经济的快速发展、国民生产总值的提高带动了人民的生活质量,人均寿命不断延长,1950-1970年我国平均生育率高达6.11(2000年生育率降低至1.22),那个年代出生的人群正逐渐步入老年,这使得我国人口老龄化、老年抚养比等问题越发严峻;近些年由于城镇经济的迅猛发展、粮食价格的降低等因素,使得农村劳动力大多选择进城务工而放弃耕地,这导致耕地利用率严重降低甚至出现弃种的现象,虽然2009年建