频率域逆时偏移的优化算法

来源 :哈尔滨工程大学 | 被引量 : 0次 | 上传用户:xuhuangyun1118
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
石油和天燃气是国家的战略物质,存储量的减少迫使人们去寻找更难发现的油气藏,对勘探技术的要求也越来越高。地震勘探、地震数据处理是寻找石油和天燃气的主要手段,叠前深度偏移是地震数据成像技术中的重要环节,主要研究方法分为两大类:一类是克希霍夫积分法,另一类是波动方程方法即基于波动理论的波场递归延拓方法。克希霍夫方法在处理不规则数据、偏移速度分析、计算效率和局部成像孔径控制方面具有优越性和更强的灵活性,缺点是有高频近似,难以处理复杂波现象。波动方程方法在处理复杂波现象如多旅行时,焦散和振幅保持方面比克希霍夫方法更自然更简单,在复杂构造下成像精度更高。基于波动方程的逆时偏移成像由于不受地下构造倾角的限制,精度更高,振幅更精确,近年来逐渐成为主流的偏移技术。逆时偏移的缺点是效率低,对有限差分法存在频散,所以该方法的研究重点在于提高计算效率和去噪。  本论文选择在频率-空间域实现逆时偏移,这样可以避免因时间差分带来的频散,理论上讲比时间域偏移具有更好的成像精度。在空间域分别采用五点和九点差分格式来逼近二阶导数,从而得到一个系数稀疏的频率域波动方程组。通过和共轭梯度法的比较发现LU分解法尽管更加消耗内存,但计算速度更快。所以本文采用了LU分解法来求解此方程组得到延拓波场,最后根据最大相关成像准则进行成像。由数值算例可以看出九点差分格式的频散明显好于五点差分,为了进一步提高计算效率,论文使用了合成震源记录技术,数值算例验证了该方法正确、有效。
其他文献
有理函数Julia集的拓扑是复解析动力系统研究的重要问题之一,多项式Julia集的连通性由于Branner-Hubbard猜想的证明[47]已得到较为完整的刻画.对于有理函数动力系统,二次有理
填充与覆盖问题是图论中非常重要而又基本的问题,在物理学、计算机网络及组合优化等领域都有十分重要的意义。填充和覆盖是一对具有对偶性质的概念。作为填充与覆盖问题中一个