波动方程时间域全波形反演方法研究

来源 :中国科学院大学 | 被引量 : 0次 | 上传用户:womjun
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
地震波形反演是一种利用地震反射数据得到地下介质特征的方法。全波形反演(FWI)是一种用全波长模拟数据与实际资料数据进行数据拟合的过程,数学上,它相当于求解一个非线性的最小二乘问题。FWI能得出地下介质物理参数的高精度高分辨率的成像。   本文将着重在三个方面的研究:1、FWI方法的本质是非线性优化问题,文中将讨论如何有效求解非线性极小化问题,比较各种方法的优劣性。2、FWI方法的非线性性会使目标函数收敛到一个错误的局部极小值,文中将研究一些策略来避免陷入局部极小值的情况。3、如果模型参数的数量庞大,那么优化方法的收敛速度将会很慢,它直接影响到算法的效率,文中将研究加速迭代收敛的方法。   本文以2D声波方程为数学模型,使用有限差分法进行波场模拟,再进行时间域的全波形反演。本文将共轭梯度法、拟牛顿法和高斯牛顿法应用于声波方程反演,并比较了各种方法的收敛性、计算量和存储量需求。由于传统的线搜索方法只利用了目标函数一阶导数信息,为了充分利用高斯牛顿法中已求得的二阶导数信息,本文将Levenberg-Marquardt法和信赖域法应用于波形反演,这两种方法都是利用二阶导数信息调整步长,构造搜索方向。文中分析并阐明了这两种步长搜索策略在应用中的优势和劣处。   对于大规模复杂模型的计算,本文以Marmousi模型为例,构造了单尺度上求解FWI的算法。为了处理局部极小值问题,在单尺度算法基础上提出了频率多尺度算法。从数值算例可以看出,频率多尺度有效利用低频信息为高频反演构造出一个可靠的初始模型,有效地解决了初值离真解较远时,反演易陷入局部极小值的问题。为了节约计算时间,本文将频域多尺度全波形反演算法区域并行化,将区域划分为多块交给不同的进程计算。此外,还结合多重网格的思想进行了频率多尺度全波形反演,充分发挥了低频反演的优势,大大地减少了计算量。   小规模反演结果的分析表明了各优化方法用于非线件问题的效率和效果,可以为其它领域求解非线性优化问题作为参考,也给以后推广到大规模算法奠定了基础;而大规模的并行反演计算结果表明了文中所构造算法的有效件,为进一步应用提供了基础。  
其他文献
本文用统计理论与方法对描述检出能力的特征量(包括响应变量的临界值、净状态变量的临界值及净状态变量的最小可检出值)进行了系统的研究.  本论文的研究内容分四个部分.第一
在英语中,副词是指在句子中表示行为和状态特征的词。长期以来,副词都是小学英语教学的一个难点。在传统的英语教学中,对副词的教学往往以模仿和机械性记忆的方式展开,这样的
本论文主要研究了一类退化的反应扩散方程的行波解的存在性稳定性.  本文主要分为三章.第一章主要介绍本文研究的模型以及反应扩散方程行波解的研究背景和国内外研究现状,并
随机时滞微分方程被广泛的用于对物理、生物、航空、材料科学等许多领域中的不确定现象的建模。在理论上求解此类方程是一件十分困难的事情。因此,数值研究此类方程变成一个十
全文共分为两章,第一章主要讨论了具负反馈的双阈值的二元神经网络模型(Ⅰ)解的收敛性与极限环的存在唯一性问题.在所设定的初始函数空间内,对于阈值σ和σ的不同取值范围给
非负矩阵分解是实现大规模数据处理与数据分析的一种非常有效的数据挖掘方法.与其它传统的矩阵分解方法(如PCA,SVD等)不同的是,非负矩阵分解算法是在矩阵中所有元素均为非负的条
约束矩阵不等式及其最小二乘问题是数值代数领域中的一个重要研究课题,它在图像恢复、控制论以及组合优化等领域中都有重要应用.此外,在很多实际应用中,不仅要求约束矩阵具有对
请下载后查看,本文暂不支持在线获取查看简介。 Please download to view, this article does not support online access to view profile.
耦合神经网络系统的同步是当今研究的主要内容之一,已经渗透到通信系统、电子系统、生物系统等各个领域,因此,许多学者致力于这个领域的研究并取得了一些重要的理论成果。特别是
随着我国教育改革的不断深入和发展,我国中职院校的教学也在改革的过程中不断取得新的进展和突破.我们在进行教育教学的过程中,经常提到要让学生德智体美劳全面发展,由此可见