论文部分内容阅读
根据地震资料的CMP道集中同相轴经部分动校正后为抛物线的假设,本文给出了加权抛物Radon变换叠前地震数据重建方法以及结合波动方程预测的抛物Radon变换表面多次波衰减方法。
地震资料的重建和数据规则化是地震资料处理的重要环节,本文建议了一种在数据域加权的Radon变换数据规则化方法:加权抛物Radon变换法(WPRT)。WPRT法通过在迭代过程中引入变化着的权系数,拓展和改进了传统抛物Radon变换方法,它可同时完成地震数据的不规则采样的规则化和空道及近偏移距道重建,且有更高的计算效率,由于加权后所形成的矩阵仍具有Toeplitz结构,因此仍可使用Levinson递推法来进行求解。本文地震数据重建的算法是后面多次波衰减方法对比研究的基础,因为基于波动理论预测的表面多次波衰减方法需要用到近偏移距的地震数据。
针对现行抛物Radon变换地震数据重建方法,本文建议了引入道均衡的提高计算效率的措施,该方法在保证计算精度的条件下可以提高地震道的插值、外推和反假频重采样的计算效率。
表面多次波是一种能量很强的相干干扰,在海洋地震勘探中表面多次波尤为明显。现在的多次波衰减方法大致可分为两类:基于波动方程预测的衰减方法(预测法)和利用地震波时距曲线特征的滤波衰减法(滤波法)。滤波法使用的一个共同假设是一次波和多次波在动校正速度上存在较大的差异。当假设条件不满足(如地层速度逆转)时滤波法不能得到理想的多次波衰减效果。预测法依据地震波传播和多次波形成的机理,在全波场的条件下能有效预测出所有阶的多次波,在衰减时通过能量最小准则用一维优化算法来求取逆源子波,进而得到表面多次波衰减的结果剖面,该方法在应用中有很多限制,如近偏移距不能缺失、震源和检波器的特性要一致、道间距和炮间距要相等。本文结合滤波法和预测法,提出了结合波动方程粗略预测和Radon域切除的表面多次波衰减方法(预测切除法),该方法利用表面多次波形成机理来预测多次波的旅行时,基于预测的多次波的CMP道集和初始地震剖面的CMP道集的对应的抛物Radon变换结果设计自适应滤波函数,可准确得到抛物Radon变换域的多次波剖面,经反变换后从原始地震记录中减去就得到衰减多次波的地震剖面。
本文将预测切除的抛物Radon变换法用于含速度逆转的模拟地震数据的表面多次波衰减的试算,结果表明Radon域的自适应滤波衰减了所有的表面多次波,而有效波的幅值得到保持。本文方法不用优化算法来消除震源子波效应,同时滤波函数是逐一CMP道集自适应来设计的,所以预测切除的抛物Radon变换法消除了预测法的优化求取子波所引入的较大的计算量和误差,也不象Radon变换法那样需要人为干预来设计切除函数。本文算法有如下优点:无震源和检波器方向一致特性的限制、计算效率较高、能在不损害有效波的前提下有效衰减所有阶表面多次波,另外,本文方法无需地下介质的任何信息并且对介质的速度没有任何假设。
本文给出了复杂SMAART地震数据的表面多次波衰减的处理流程和处理结果,通过近偏移距和中偏移距剖面以及炮记录和CMP道集的结果显示本文方法可以衰减全部偏移距的表面多次波,对比滤波衰减掉的多次波,本文方法对于复杂的数据仍能保持住有效波的振幅。另外本文给出了近偏移距缺失时表面多次波预测和衰减试算结果,该结果说明本文对表面多次波衰减预处理中的近偏移距外推没有严格的要求,所以本文方法的实用性很强。