区域分解频率域波动方程波场模拟及反演

来源 :中国科学院大学 | 被引量 : 0次 | 上传用户:maomao1983520
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文具体研究了,基于有限差分方法,如何应用区域分解方法来求解频率域波动方程。所求解的波动方程包括声波方程和弹性波方程。   数值求解频率域声波方程,即Helmholtz方程时,当方程中的参数波数k很大时,通过各种离散得到的离散矩阵的条件数很大,而对于这种问题这时用直接法或迭代法求解一般所得到结果都不理想。我们知道对于这类坏条件数问题,通常需构造预条件子去求解。构造预条件子的方式通常有两类,一类是基于矩阵的构造方法,如通过ILU方式构造,另外一类是基于偏微分算子的构造方法。   本文第三章简单回顾了第二类基于偏微分算子的预条件子的构造方法,主要是复共轭移位拉普拉斯预条件(CSL)方法,并将这种预条件方法用到了Helmholtz方程的区域分解求解,包括非重叠型和重叠型区域分解算法,对均匀模型、层状模型和Marmousi模型做了数值计算。   在对Helmholtz方程区域分解求解的基础上,第四章研究了频率域弹性波方程的区域分解求解,包括非重叠和重叠区域分解方法。详细阐述了如何基于有限差分方法去构造界面条件,以及如何转化为子区域问题,并进行有限差分离散,最后通过预条件GMRES法和Bi-CGSTAB迭代法来求解,也进行较多了的数值计算。   最后,在第五章,初步探讨了频率域波场反演问题。该问题通过在求解区域中或边界上的某部分的观测值,来反演出最好的模型使得计算值和观测值之间的误差尽量小,这是一个非线性最小二乘问题。我们结合最优化计算方法中的Gauss-Newton法、Levenberg-Marquardt方法、拟牛顿法等来求解这个最优化问题,并作了数值计算。  
其他文献
多项式系统求解问题是数学中一个既古老又经典的问题,并在科学与工程计算中有着广泛的应用,例如机器人技术,编码理论,最优化理论,数学生物学,计算机视觉,博弈论,统计学,机器学习,控制理
学位
本文研究了几类二阶椭圆偏微分方程(组)的Dirichlet边值问题.其中涉及两类微分算子,即Monge-Ampère算子和Laplace算子.前者是典型且重要的完全非线性二阶微分算子,与这类算子相
本文首先对利率期限结构的经典理论和模型进行了介绍,然后对利率期限结构实证研究相关文献做了综述,特别是利率期限结构与宏观经济变量之间相互作用的理论依据与实证进展。在实
有限维代数可分为两大类:有限表示型和无限表示型,而无限表示型代数又可分为驯顺表示型和野表示型这互不相交的两类.目前,驯顺型代数的研究是代数表示论研究的重点.对驯顺型代
文章简略地介绍了多目标规划问题的发展过程、基本理论、求解方法以及基于最小二乘原理的娄据拟合法,详细分析并比较了目前关于多目标规划求解的四种方法及其优缺点,特别对今
本文主要研究的是关于Stein流形上的李群作用的几个问题.设K是一个紧李群,KC是它的复化,证明了一个n维的具有全纯的群KC作用的连通Stein流形,如果具有K-不动点且它上的所有K不变
本文介绍了Kyle(1985)内部交易模型,理性预期价格内部交易模型及其相关扩展.在理性预期价格模型的框架下,构建了资产服从对数正态分布的风险中性单内部交易者模型,并将高频交易
学位
利用已构建的“云粳17号/云粳20号”重组自交系(RILs)群体为试验材料,在低纬高原中国云南和高纬平原韩国水原两地间开展了粳稻主要农艺性状的异地穿梭鉴定,并分析了粳稻RILs