多子波克希霍夫地震数据偏移方法
【技术领域】
[0001] 本发明涉及地震资源（石油，煤炭，矿产等）勘探中的地震资料处理和解释技术领 域，属于一种多子波克希霍夫地震数据偏移方法。
【背景技术】
[0002] 地震数据偏移是地震数据处理过程中的一个重要环节。目前地震数据都是以某个 地球平面振动的振幅按等时间隔采样来表述的，一个地震道代表了平面上一个点在一个时 间段的振动，实际表述为一个浮点数数组。因此，目前常规的地震资料处理都是基于，或建 立在振动样点之上的。在基于样点的偏移计算中，不可避免地用到对振动场的样点插值，从 而引起频散噪音和地震波动或子波的拉伸变形，造成进一步地震资料解释的困难，也有可 能造成错误的解释。
[0003] 近年来在国际上新出现了基于小波变换的偏移计算方法，例如YuZhou等作者于 2004年提出了基于小波变换的叠前多尺度（multiscale)克希霍夫偏移的方法。这种方法 首先应用小波变换将地震道分解成不同尺度（scale)的小波，然后应用常规的克希霍夫偏 移算法移动小波的系数，然后根据小波系数重构共反射点道集，叠加得到偏移结果。虽然这 种方法可降低频散噪音，但也存在以下问题：
【主权项】
1. 一种多子波克希霍夫地震数据偏移方法，其特征在于：它包括下列步骤： (1) 、生成多地震子波序列，即不同主频或最大振幅频率的地震子波序列； (2) 、利用上述地震子波序列对每个地震道进行分解，得到I1范数优化的稀疏地震子波 系数序列，提高对地震道的解析程度； (3) 、选定所要偏移的子波的主频率，或所要偏移的子波的主频率范围； (4) 、对波克希霍夫偏移算法中，将给定成像时间计算输入样点时间的计算方法改成给 定输入子波时间计算该子波在成像中的时间，然后将输入地震子波准确地移动到成像中的 位置，然后重构叠加得到精确的偏移结果； (5) 、对移动后的共反射点地震道进行重构叠加得到精确的偏移结果，得到所给定子波 的偏移成像；理论上，对给定子波序列中的任何一个子波，都可以得到一个成像输出，该输 出代表地震数据中与该子波对应的信息的成像；同样，得到不同给定主频区间的子波的成 像；将部分主频区间的成像叠加，就等于该主频区间的子波的成像；所有主频不重叠的成 像叠加，就是所有输入信息的成像。
2. 根据权利要求1所述的多子波克希霍夫地震数据偏移方法，其特征在于： 在步骤（1)中生成多地震子波序列中给定一个子波，通过三次样条函数插值和二进制 搜寻计算得到给定子波最大振幅频率或主频的子波的计算方法。
3. 根据权利要求2所述的多子波克希霍夫地震数据偏移方法，其特征在于： 在步骤（2)中利用一个子波序列，求满足I1范数优化的稀疏地震子波系数序列的方 法。
4. 根据权利要求3所述的多子波克希霍夫地震数据偏移方法，其特征在于： 在步骤（4)中由输入子波，移动子波到共反射点道集中的地震道，从而达到偏移的结 果。
5. 根据权利要求4所述的多子波克希霍夫地震数据偏移方法，其特征在于： 在第四步中由输入子波时间t，计算移动到共反射点道集中的地震道的时间t(i的方法。
6. 根据权利要求5所述的多子波克希霍夫地震数据偏移方法，其特征在于： 在步骤（1)至步骤（5)中得到对应于给定子波主频或主频范围的偏移结果的方法。
【专利摘要】一种多子波克希霍夫地震数据偏移方法，它包括下列步骤：(1)、生成多地震子波序列；(2)、利用多地震子波序列对地震道进行分解，得到l1范数优化的稀疏地震子波系数序列；(3)、选定所要偏移的子波的主频率；(4)、对波克希霍夫偏移算法中，将给定成像时间计算输入样点时间的计算方法改成给定输入子波时间计算该子波在成像中的时间；(5)、对移动后的共反射点地震道进行重构叠加，得到所给定子波的偏移成像；本发明的有益效果是：它提高地震数据偏移精度，以得到多地震子波不同波段的，更加精确的地下岩层和物理特征的分布和变化的偏移结果，为下步对资源(油气，煤炭，矿产等)的解释，提供更可靠的资料。
【IPC分类】G01V1-28
【公开号】CN104570100
【申请号】CN201510056165
【发明人】安平, 林禄春
【申请人】杰奥世博（北京）技术有限公司
【公开日】2015年4月29日
【申请日】2015年2月3日