改进的郭涛算法在剩余静校正中的应用.pdfVIP

石油天然气学报 (江汉石油学院学报)2014年2月第36卷第2期 JournalofOilandGasTechnology(J．JPI)Feb．2014Vo1．36No．2 改进的郭涛算法在剩余静校正中的应用 曹中林 ，周强，吕文彪 张华，陈爱萍 (中国石油集团川庆钻探工程有限公司地球物理勘探公司，四川成都。’。’。) [摘要]在低信噪比地震资料处理中，剩余静校正是其中的关键处理环节。传统的剩余静校正方法容易陷 入局部极值 ，难 以实现全局寻优 ，从而造成叠加剖面成像效果不佳 。为此，提 出了一种高精度 的非线性 剩余静校正方法，该方法以互相关值最大为 目标 函数 ，采用 了逐步缩小范 围的郭涛算法进行剩余静校正 量的全局寻优。通过理论模型和低信噪比资料的测试表 明，该方法优于传统的剩余静校正方法，能够快 速逼近全局最优解 ，并且能够有效改善低信噪比地震 资料的成像效果 。 [关键词]改进的郭涛算法；反演；剩余静校正 [中图分类号]P631．44 [文献标志码]A [文章编号]i000—9752 (2014)02—0062一O4 我国西部山地地 区，地形起伏多变 ，地表结构复杂 ，山体部分老地层出露地表，导致采集的地震数 据静校正问题特别严重，很大程度上影响了地震数据 的最终成像效果 。地震资料在做了野外静校正之 后，有效反射波的同相轴应该呈双曲线分布，但实际上由于静校正计算中难免存在误差，使地震记录中 还残 留一定的时差 ，有效反射波的同相轴就不是理想的双曲线 ，随后 的叠加成像效果也就不理想 。因 此 ，有必要将剩余时差在速度分析和叠加之前予以去除，这就是剩余静校正的主要任务 。 目前 ，常规的剩余静校正方法有 ：共地面点法、Satan算法和最大能量法口 ]，都不能较好地解决 低信噪比地震资料的剩余静校正问题 。因此大量 的学者引入 了非线性反演方法进行剩余静校正量的求 取 ，例如遗传算法_3]，模拟退火方法_6 以及蒙特卡罗方法[1O~11]，这些非线性反演方法在一定程度 上改善了地震资料的最终成像效果，但是都存在反演效率较慢以及容易陷入局部极值的问题 。而郭涛算 法l】纠是一种基于子空间搜索 (多父体杂交)的群体随机搜索算法 ，可 以用来求解复杂函数优化问题 。 笔者结合了郭涛算法的特性 ，提 出了逐步缩小搜索范围的改进策略，并将其应用到剩余静校正量的求取 中，实现了剩余静校正量的全局寻优 。理论模型和实际地震数据的测试表明，该方法优于传统的剩余静 校正方法 ，针对低信噪比地震资料处理具有一定的效果。 郭涛算法特点及程序 郭涛算法是一种精确的全局寻优算法_l】”]，该算法有以下 3个特点： 1)采用了演化计算中群体搜索策略，可以实现全局寻优 ，适合于复杂的非线性函数优化的问题 。 2)进行了多父体的线性重组 ，特别是子空间中随机搜索的非 凸性，保证 了随机空间中随机搜索 的 遍历性 。从N个父体向量P一 {x ，X ，…，x }中随机选择M个个体组成 向量x 一 {X ，X ，…，XM}， 构成子空间V= {X 1口·X}作为其搜索空间，其中n一 {n，“，…，nM}是M维向量，满足条件： M ∑ “一1 —0．5≤ai≤1．5 (1) z一 】 3)采用了 “劣汰策略”，每次只把群体中适应性最差 (目标函数值最大)的个体淘汰出局 ，淘汰压力最 小 ，既保证了群体的多样性，也保证 了适应性最好 (目标 函数值最小)的个体可 以 “万寿无疆”。这种 “群体 [收稿 日期]2013—06—04 [基金项 目]中国石油天然气集团公司科学研究与技术开发项目 (201113-37O6)。 [作者简介]曹中林 (I981一)，男，2004年大学毕业 ，硕士 ，工程师，现主要从事地震资料处理方面的方法研究与软件开发工作。 第 36卷第 2期 曹中林 等 ：改进的郭涛算法在剩余静校正中的应用 爬 山策略”，保证 了整个群体最后集体登上最高峰(深谷)。 郭涛算法的具体程序如下 ： Algorithm GT ： Be