论文部分内容阅读
当下许多油气藏都进入了精细勘探与开发阶段,这就要求地球物理工作者对地下介质物性参数的变化做出更加精确的反馈与成像,而地球物理勘探本质实际上是一个地球物理反问题求解过程,在求解过程中,所采用的反演方法对反演结果的好坏影响很大。地震波阻抗反演技术是其中最具代表性的反问题之一,根据反演过程中使用的地震资料类型的不同可以分为叠前波阻抗反演和叠后波阻抗反演,现阶段,叠后波阻抗反演方法种类繁多,技术也较为成熟,研究不同反演方法的优缺点、适用性、影响反演效果的各种因素和反演中的控制参数的选取,不仅能为实际应用提供指导,而且能避免应用的盲目性。最优化理论在求解地球物理反问题中的应用历史可以追溯到二十世纪八十年代,它通常是通过设置目标函数,在全局或局部范围搜索以寻求目标函数极小点或极大点所对应的变量值,从而达到求解的目的。由于地球物理学中的待求解方程组往往是非线性和高度病态的,且具有严重的多解性,在求解这类方程组时,通常有两种方案,一是通过把地球物理非线性问题转化为线性问题,用线性化方法求解非线性问题,另一种思路是采用完全非线性方法求解。与线性反演理论相比,非线性反演无论在理论上还是在实现方法上都比线性反演困难得多,而线性反演在方法原理上、可操作性等方面都更为成熟。阻尼最小二乘法,因其阻尼系数的引入,在保证了矩阵正定性的同时,矩阵的病态性得以改善,使得该算法在寻优策略上优于高斯牛顿法和最速下降法,但是如何给定阻尼系数一直以来是影响该算法稳定性和收敛性的重点问题。本文以稀疏约束脉冲反演方法为基础,研究了信赖域方法下的阻尼系数给定方式,并建立最小二乘意义下的波阻抗反演目标函数,降低了传统最小二乘法中雅克比矩阵的计算代价,并且解决了雅克比矩阵的非正定问题和奇异问题,优化了阻尼系数的迭代方式,使得算法具有更好的稳定性和全局收敛性。分别对不同子波、不同噪音、不同初始波阻抗的理论模型进行了测试,结合二维Marmousi模型,进行反演均取得较好反演结果。应用实际资料进行反演,反演结果较好,证明了本文提出方法的可行性和有效性。