张 钋,汪道柳,徐 昇,周洪波,刘金俊

(Statoil Gulf Service,Houston,77042)

对于复杂地质构造区域的叠前深度偏移成像来说,除了选取适当的偏移算法之外,高精度的速度模型是获得良好成像质量的另一关键因素,而如何准确估算地层的各向异性参数则一直是速度建模中的难点及研究重点。传统的基于层析成像反演的地震速度建模方法的基本原理是通过数学优化过程调整初始速度模型使得实际观测数据与模拟数据之间的差异(即优化过程的目标函数)最小化,进而求出某种数学意义下与实际数据吻合的最佳速度模型。对于不同的层析成像算法及不同的已知数据,目标函数的选择往往也不尽相同,例如,如果已知数据是拾取的地震波旅行时,则目标函数可定义为拾取的旅行时与正演模拟的旅行时之差;如果数据是叠前深度偏移后输出的共成像点道集,则目标函数可定义为相应成像道集中同相轴的拉平程度。众所周知,与大多数反问题类似,地震层析成像反演通常是欠定和病态的[1-2],而且该问题在反演各向异性介质模型参数的情况下尤为严重。另外,由于地震资料采集的技术、环境以及成本等诸多因素的限制,实际地震数据往往只是空间有限采集范围内的地表数据,这也造成层析成像反演不同模型参数之间,以及模型参数和地层深度之间存在一定的不确定性,即存在多种模型参数组合可以使得观测数据与模拟数据在一定误差范围内相吻合[3-5]。……