[发明专利]补偿吸收衰减的叠前时间偏移方法

摘要

补偿吸收衰减的叠前时间偏移方法，应用于地震勘探中反射地震资料处理。该方法在偏移过程中补偿地球介质粘性、薄层散射导致的高频地震波幅值衰减，恢复被衰减的高频成份，使得中、深层构造的成像分辨率达到与浅层接近的程度，且可保证稳定性并避免噪音放大。该方法不对Q值的空间分布做层状或均匀的假设，通过引入描述地震波幅值吸收衰减的等效Q值，将补偿吸收衰减与叠前时间偏移有效地结合到一起。该方法能基于地面采集的反射地震资料确定等效Q值，规避了现行各类吸收衰减补偿方法在确定Q值上的困难。该方法能更好地指示中、深层构造的断裂和地层沉积样式，对油气、矿产资源勘探有重要应用价值。

主权项

1.一种补偿吸收衰减的叠前时间偏移方法，其特征在于：包括以下步骤：步骤A、用单条拖缆或测线记录人工震源激发的反射地震信号，记录到磁带上；步骤B、从磁带上读取地震信号，对叠前地震资料做常规的噪音衰减处理；针对部分共中心点，抽取共中心点道集，对抽取的道集作常规的动校正速度拾取，对所得到速度的倒数做横向平均，将得到的横向均匀的速度场作为初始偏移速度；步骤C、将噪音衰减处理后的叠前地震资料按偏移距排序、分组，将不同组地震资料放置到集群计算机的不同计算节点上；在每个计算节点，依据初始偏移速度，用常规的叠前时间方法进行偏移计算，收集各计算节点的偏移结果，形成共反射点道集；步骤D、对共反射点道集，利用初始偏移速度做反动校，再做动校正得到更新后的速度，对各个CDP点更新后的速度的倒数做空间平滑，所得到的平滑速度场即是最终的偏移速度场；步骤E、依据偏移速度场，对噪音衰减处理后的叠前地震资料做常规叠前时间偏移计算；收集各计算节点的偏移结果，形成共反射点道集，对共反射点道集做剩余动校和拉伸切除，得到偏移叠加剖面；在经增益处理后的偏移叠加剖面上选定清晰、连续的同相轴；根据给定的采样间距，确定等间距分布的拾取Q值的CDP点，依据选定的同相轴，在这些CDP点上确定从浅到深的幅值对比点和浅层的参考点，并在各点处选定应用短时傅立叶变换的时窗；从浅到深的幅值对比点处的同相轴需展现由细变粗的变化趋势；步骤F、对噪音衰减处理后的叠前地震资料应用时间域快速傅立叶变换，求半导数、施加频率窗衰减带、计算地震道能量、剔除奇异道；步骤G、确定等效Q值的取值范围，按1/Q等间距选取系列Q值，用每个Q值，对步骤F处理后的频率域叠前地震资料，做稳定的压噪反Q滤波计算，得到新的时间域地震资料；步骤H、依据偏移速度场，对每个Q值对应的反Q滤波后的地震资料做常规叠前时间偏移计算；收集各计算节点的偏移结果，形成共反射点道集，对共反射点道集做剩余动校，切除剩余动校后的拉伸部分；将需拾取Q值的CDP点的相邻共反射点道集叠加，形成各CDP点的超道集；定义小、中、大三个射线参数p₀、p₁、p₂，对全部需拾取Q值的CDP点的幅值对比点和参考点计算偏移距其中i＝0，1，2，T为该点的单程旅行时，单位秒，v_rms为该点的偏移速度，单位米/秒，射线参数p_i的单位是秒/米，偏移距d_i的单位是米；在各点选定的时窗内，按偏移距范围[0，d₀]、[d₀，d₁]和[d₁，d₂]对超道集分别进行叠加，在每点处得到三个短时时间序列；步骤I、重复步骤G和步骤H，收集各幅值对比点和参考点的对应小、中、大射线参数的三个短时时间序列，按1/Q值大小排列，在各点形成三个随1/Q值变化的小、中、大射线参数的短时Q道集；步骤J、根据小、中、大射线参数的短时Q道集拾取等效Q值，得到成像区域的等效Q值场；步骤K、在每个计算节点，对步骤F处理后的地震资料的全部地震道循环，应用稳定的变频带频率域算法计算各成像点的频率域波场实部；步骤L、由各成像点频率域波场的实部，用变频带、保能量成像条件求得各成像点的偏移幅值；步骤M、将偏移幅值按偏移距大小累加到存放偏移结果的数组中对应的偏移距上；收集各计算节点的偏移结果，形成全部CDP点的共反射点道集；对该道集做剩余动校正和拉伸切除，将不同偏移距的偏移结果叠加，形成偏移叠加剖面；步骤N、通过显示软件将偏移叠加剖面数值转换为地下反射构造的剖面图像，剖面图像将指示中、深层构造的微小断裂体系和精细的地层沉积样式，用于确定地下生、储油构造和识别油气储层。