基于压缩傅里叶插值算法的地震数据规则化技术

基于压缩傅里叶插值算法的地震数据规则化技术

随着地震勘探的深入，地震数据的保幅要求也越来越高。在进行叠前属性分析时,如果叠前道集数据振幅不均匀、一致性差,则会严重影响属性分析结果的准确性。可以说,如今地震资料处理的目标不仅仅是提高最终成像的质量,还应该获得振幅相对保真、一致性较好的叠前道集。然而,在实际地震资料采集过程中,由于地表条件复杂,变观现象普遍存在,炮检点分布不规则导致覆盖次数不均匀,影响了地震数据振幅的一致性关于插值技术的研究,国内外学者已经做了很多工作针对上述问题,本文采用计算效率较高的压缩傅里叶插值算法,插值时考虑了地震数据的主测线(Inline)、联络测线(Crossline)、时间、偏移距和方位角等五个维度的信息,同时针对压缩傅里叶算法不抗假频的问题,增加了局部的倾角校正处理以提高压缩傅里叶插值的抗假频能力,实现了对不规则地震数据的五维插值规则化处理,利用模型数据和实际地震资料对方法进行了测试分析。1空间带宽波数数据分析根据传统的Shannon采样定理,对于任意最高频率为f式中:基函数sinc(x)=sin(x)/x。式中:s(x式中:xCHEN等以上插值研究的内容都基于固定带宽的信号重建方法,而实际地震数据在不同插值位置的邻近地震道分布各不相同,用相同带宽插值显然不够精确。为此,MOORE式中:y代表插值期望输出位置。对任意指定的空间位置插值时,首先根据数据空间采样间隔确定最大奈奎斯特波数,以此奈奎斯特波数为基础,按一定间隔不断缩小空间波数,并根据公式(4)分别计算这些波数对应的最小均方根误差,优选出满足设定的最小均方根误差门槛值(例如10%)的最大波数,然后根据公式(5)计算该波数对应的插值权系数。这一过程也可以视为对空间带宽波数k的扫描,通过波数扫描确定最优的插值权系数,具体的波数扫描流程如图1所示。目前大部分插值技术利用横向、纵向和时间3个维度的三维数据插值,忽略了数据在不同偏移距和不同方位角之间的联系,插值之后的数据不能很好地保真、保幅以及保持偏移距和方位角信息,影响了插值结果的振幅一致性。本文采用了五维插值技术,将插值数据分选为共中心点道集,同时考虑不同面元数据本身包含的主测线、联络线、时间、偏移距和方位角等五维信息。计算插值权系数的过程与二维、三维数据插值类似,每一维都相对独立地进行波数扫描,确定各自最优插值波数,最终插值公式如下:式中:s为期望插值位置的输出道;i代表不同的插值维度;s为避免空间假频影响,常规的压缩傅里叶算法在插值前需要对数据进行动校正处理。但是在实际应用中,如果数据含有多次波或动校正速度出现偏差,则部分同相轴无法完全校平,插值时仍然会受到空间假频的干扰。因此,本文利用五维压缩傅里叶算法插值时,增加了倾角校正处理,倾角的求取思想与生成速度分析的速度谱方法类似,通过设定倾角范围,对期望插值位置的每个时间样点进行倾角扫描,计算出不同倾角的相干值,优选出相干值最大的倾角。增加倾角校正的插值算法保证了插值时同相轴的局部水平,插值公式如下。式中:Δt2数据分析结果为了检验增加倾角校正的五维压缩傅里叶插值算法抗假频的能力,建立一个含有陡倾角的三维速度模型,如图2a所示,层速度分别为1500,3000,6000m/s。对该模型进行三维地震波有限差分正演模拟,上边界为自由边界条件。任意选取一炮地震数据,如图2b所示,对该炮点位置进行五维插值。插值的主要参数:波数扫描范围为最大奈奎斯特波数(由插值空间采样间隔确定)的100%~0,扫描间隔5%;最小均方误差10%;倾角扫描范围12ms,扫描间隔2ms(数据采样间隔2ms)。对比分析了常规压缩傅里叶算法与含倾角校正压缩傅里叶算法对含假频数据的插值效果。图2b为实际模拟的炮集记录,包括第一、第二个反射界面的一次反射波,以及来自第一个界面的一阶和二阶全程多次波,此外还包括来自倾斜界面的反射波。图2c为采用常规只经过动校正处理的五维压缩傅里叶算法得到的插值数据,可以看到,第一个反射界面反射波近偏移距部分和第二个反射界面的反射波插值效果较好,但是第一个反射界面的远偏移距部分反射波、多次波以及陡倾角界面的反射波插值效果不理想。这是因为利用速度模型的一次波速度对原始数据进行动校正时,无法校平多次波和倾斜界面的反射波,同时浅层反射波的远偏移距还存在动校拉伸现象,插值时受到空间假频的影响。图2d为增加了倾角校正的五维压缩傅里叶算法插值结果。由于进行了倾角校正,保证了插值时同相轴局部水平,因此取得了较好的插值效果。对比图2b和图2d可以看出,增加倾角校正的五维压缩傅里叶插值算法获得的插值数据与理论数据差别不大,波场干净,反射波、多次波以及倾斜界面反射波与原始炮集非常接近,插值信息逼真,空间假频影响得到了较好压制。3推进各维插值空间分布的规则化某工区地表条件复杂,观测系统变观位置较多,地震数据采集的炮点、检波点分布如图3a所示,红色代表检波点位置,蓝色代表炮点位置。受地表因素影响,炮、检波点分布很不规则,设计满覆盖次数为42次,但为了避让障碍物而采取了变观措施,实际采集的地震数据覆盖次数从10多次到70多次不等,如图3b所示。最小均方根误差门槛值是压缩傅里叶插值的重要参数,要根据实际数据的信噪比特点进行选取,原始数据信噪比高时,最小均方根的门槛值可以相对较大,原始数据信噪比低时,最小均方根的门槛值要相对较小,本次实验选取的门槛值为10%。波数扫描间隔对计算效率有很大的影响,扫描间隔过小时计算时间较长,扫描间隔过大时影响插值精度,本次实验对计算效率和插值效果进行了权衡,选取的波数扫描间隔为10%。其它插值参数与模拟数据测试参数一致。图4a为五维压缩傅里叶插值方法规则化后数据的炮、检点位置分布图,可见规则化后数据炮、检点分布更加规则。图4b为规则化后数据的覆盖次数属性图,可见数据规则化后覆盖次数的均匀性更好,部分缺失炮点得到弥补。图5为插值前后的CMP道集数据,红色实线代表炮检距分布范围,对比可以发现,五维插值规则化处理后炮检距分布更加均匀,原始数据近偏缺失得到弥补。为了更加准确地分析五维压缩傅里叶插值对炮检距的规则化作用,对该CMP道集数据作了炮检距分布属性分析,如同6所示。可以看出,原始数据存在大量变观现象,炮检距分布稀疏、不均匀,近炮检距数据缺失,而经过五维压缩傅里叶插值后的炮点位置分布规则,炮检距分布更加均匀。分别对五维压缩傅里叶插值前后数据进行克希霍夫叠前时间偏移处理,图7对比了插值前后叠前时间偏移处理得到的CRP道集。图7a是五维压缩傅里叶插值前CRP道集,由于原始数据炮点、检波点空间分布不规则,导致CRP道集振幅出现跳动,一致性较差,浅中层同相轴杂乱、不连续,不利于后续的叠前属性分析及AVO处理(图中红色椭圆线圈内)。图7b是五维压缩傅里叶插值规则化处理后的CRP道集,浅中层同相轴连续性增强,道与道之间的能量更加均衡、自然,振幅相对变化关系趋于合理。这表明经过五维压缩傅里叶插值处理后的数据炮点、检波点空间分布规则,CRP道集的振幅一致性得到改善,道集质量得到提高。图8对比了五维压缩傅里叶插值前后的叠前时间偏移剖面。图8a是五维压缩傅里叶插值前偏移叠加剖面,由于炮检点分布不规则,存在变观、空洞等现象,浅层同相轴受到严重影响,表现为部分同相轴连续性较差(图中红色椭圆线圈内)。图8b是五维压缩傅里叶插值后的叠前时间偏移叠加剖面,对比可以发现,利用五维插值技术对炮检点位置进行规则化处理后,由于部分缺失炮点得到弥补,浅层同相轴的连续性得到明显改善。同时,由于五维插值技术改善了CRP道集质量(图7),五维插值后的偏移叠加剖面信噪比得到提高,一些复杂的断层构造更为清晰(图中黄绿色方框内)。为了进一步考察五维压缩傅里叶插值规则化处理对振幅一致性的影响,对五维压缩傅里叶插值前后数据体进行了振幅属性切片对比分析(图9)。图9a是插值前300ms处的振幅切片,由于炮检点分布不规则的影响,振幅切片存在空洞。图9b是五维压缩傅里叶插值后300ms处的振幅切片,经过五维插值处理后,浅层空洞得到了一定弥补,对地震属性边界的刻画更加清晰。图9c和图9d是五维压缩傅里叶插值前后1300ms处的振幅切片,对比可见,插值前由炮检点分布不规则引起的采集脚印在插值后得到了很好的消除(图中红色线圈内)。对主要目的层T4维插值规则化空间采样点分布本文采用压缩傅里叶插值算法,考虑了地震数据的主测线、联络测线、时间、偏移距和方位角等五个维度的信息,同时增加了局部倾角校正处