一种时移地震互约束频率一致性处理方法
【技术领域】
[0001] 本发明属于油气勘探地震资料处理领域，特别是是一种时移地震频率一致性处理 方法。 现有技术
[0002] 时移地震主要目的是W地震的手段监测油气储层在油气开发过程中储层参数发 生一系列的变化。但是在实际处理过程中，由于不同时期资料采集过程中各种外部因素的 影响，会使时移地震剖面上非油藏部分带有不应该有的差异。时移地震互约束频率一致性 处理的目的是通过对同一区块不同时间采集的地震资料的处理，消除各种因素的影响，使 两期资料的频率有合理的同一性和差异性，保证不同时期地震响应的频率差异仅来自于 油气储层的变化。
[0003] 目前，时移地震互约束频率一致性处理方法不是很多且各有优缺点，目前，检索到 的时移地震互约束频率一致性处理技术有W下几种：
[0004] 1、不同步长反權积处理方法；首先对两期资料都进行不同反權积步长的测试处 理；对两期资料的测试结果进行分析、对比；通过分析、对比结果选取一组能够使两期资料 频率最为接近的反權积参数对资料进行反權积处理。此方法的缺点：参数的选取更多是靠 经验判断；处理中没有很好的体现两期资料的互约束。
[0005] 2、匹配滤波技术；把作为参考标准的一组资料称为基础资料，要校正的资料称为 监测资料，郝么匹配滤波技术可W描述如下：首先，在基础资料和监测资料中不包含地下 油藏属性变化的区域内取一个时间窗；然后，在此时间窗内让两种资料尽可能匹配，计算 出匹配滤波器；最后，用此滤波器作用于监测资料。此方法的缺点；全区用同一个校正因 子，无法满足实际资料的复杂变化情况；处理中没有很好的体现两期资料的互约束；匹配 滤波器不仅仅受频率控制，还包含相位、能量等因素，一致性处理因素不单一。
[0006] 3、滤波法；其做法就是设计一个滤波器，对两期资料进行滤波处理。此方法的缺 点：虽然可W统一两期资料的频宽，但无法统一两期资料的频率构成；参数选取缺乏科学 监控手段；易损伤储层变化引起的频率差异。
【发明内容】
[0007] 本发明的目的是针对现有技术存在的问题，提供一种能够实现空变的互约束频率 一致性处理消除非储层因素引起频率变化，从而提高时移地震互约束频率一致性处理的精 度，继而提高时移地震对储层变化的监控能力的时移地震互约束频率一致性处理方法。
[000引本发明的处理方法包括：
[0009] ①输入地震数据；输入时移地震两期数据，数据是经过时差校正、能量校正处理后 的共中必点道集或偏移后的共反射点道集。
[0010] ②数据匹配抽取；按照共中必反射点或共反射点为标识道头，对时移地震两期资 料进行分选，将两期地震资料具有相同共中必反射点或共反射点的数据放在同一个共中必 反射点或共反射点域内。
[0011] ⑨自相关计算：设计一个资料信噪比较高且不包含具有储层变化的反射信息位置 的计算时窗，在设计的自相关时窗内对所有地震道进行自相关计算处理，得到每一个地震 道的自相关结果。
[0012] ④自相关结果互均化计算；将同一个共中必反射点或共反射点域的两期地震数据 的自相关结果用η次方根的方法进行互均化叠加，使每一个共中必反射点或共反射点只有 一个互均化后的自相关结果。
[0013] ⑤求取滤波算子；利用步骤④计算出的自相关结果，求取其反算子，将送个反算子 做为滤波算子。
[0014] ⑧空变滤波；利用步骤⑤得出的滤波算子，用每一个滤波算子对相应共中必反射 点或共反射点的两期地震道进行權积处理，就可W达到两期地震资料互约束频率一致性。
[0015] ⑦输出：将滤波处理后的数据输出，并对两期时移地震数据进行分选输出，完成互 约束频率一致性处理。
[0016] 本发明的技术方案进一步包括：
[0017] (1)时移地震两期资料数据分选
[0018] 按照共中必反射点或共反射点为标识道头，对时移地震两期资料进行分选，将两 期地震资料具有相同共中必反射点或共反射点的数据放在同一个共中必反射点或共反射 点域内。
[001引 0)自相关计算
[0020] 自相关计算的结果取决于计算时窗内是否包含储层变化信息及时窗内地震数据 信噪比高低，所W在自相关计算时，计算时窗选择资料信噪比较高且不包含具有储层变化 的反射信息位置。在设计的自相关时窗内对所有地震道进行自相关计算处理，得到每一个 地震道的自相关结果。
[0021] 设e(t)为地震记录，其自相关函数为Ητ)，τ为时间延迟。自相关是一个时间 序列在不同时间点上波形相似性的量度，送是一个移动求和过程，用下式表示：
[0022]
[0023] (3)自相关结果互均化计算
[0024] 将同一个共中必反射点或共反射点域的两期地震数据的自相关结果进行互均化 叠加，使每一个共中必反射点或共反射点只有一个互均化后的自相关结果。设m为一个共 中必点或共反射点自相关结果个数，U为自相关结果，η为设定的次方根参数，均化结果为：
[00巧]
[002引 （4)求取滤波算子
[0027] 对互均化后的每一个自相关结果求取其一个反算子，将送个反算子做为滤波算 子。设存在一个滤波算子f，滤波算子与（2)式自相关结果權积可W得到一个尖脉冲函数δ，即： ^0028] 2 ···························································· (3)
[002引则f可W表示为：
[0030]
[0031] f为自相关结果S的反算子，就是所求取滤波算子。
[00础 妨空变滤波
[0033] 用每一个滤波算子对相应共中必反射点或共反射点的两期地震道进行權积处理， 就可W达到两期地震资料互约束频率一致性。
[0034] 本发明有着其他方法不具备的优势，其具体优势和特点表现在W下几个方面：
[0035] 第一、方法效果的可靠性。本方法通过对时移地震两期资料进行多道统计的方法 求取两期资料共同的滤波算子，通过滤波算子的应用实现两期资料的频率匹配，通过处理 实例证明效果明显且稳定。
[0036] 第二、实现了两期资料的互约束处理。由于本方法在滤波算子求取过程中采用两 期资料多道统计的方法，从而在频率一致性处理过程中真正的实现了互约束处理的目的。
[0037] 第H、可实现两期资