彭仁艳,徐振旺,刘文峰,余 锋,董旭光
1.中国石化石油工程地球物理有限公司科技研发中心,江苏 南京 210000
2.中国石油辽河油田分公司勘探开发研究院,辽宁 盘锦 124010
3.中国石化石油工程地球物理有限公司华东分公司,江苏 南京 210000
4.中国石化石油工程地球物理有限公司国际业务发展中心,江苏 南京 210000
地震同相轴自动识别与拾取在速度分析、资料解释等方面具有重要的应用价值,但对于大规模三维地震数据,人工拾取工作量巨大,效率低下,无法满足工业生产的需要[1-8]。在地震剖面上自动识别与拾取同相轴方面,前人已进行了大量研究,基本思想为基于数据的空间分布特征,利用局部相似性寻找线性相关的同相轴[9-14]。根据计算局部相似性算法的不同对拾取方法进行分类,目前,工业界中使用最广泛的是基于倾角信息的拾取算法,即根据地震数据的空间局部相似性,估计整个数据空间的倾角场,然后,根据参考道或种子点,沿着估计的倾角方向依次追踪相邻道,最后,将拾取的点进行模式分类,构造出相应的反射层位。其中,描述局部相关性,计算倾角的方法是整个自动拾取的关键步骤,其计算方法也有很多种[15],如利用瞬时波数方向估计反射层的倾角[16]、利用平面波分解或预测计算局部空间的倾角[17]、利用结构张量的方法估计倾角[18]等。平面波分解假设同相轴是局部线性的,计算精度不高;结构张量对数据比较敏感,易受噪音的干扰。另……