本发明属于石油地震勘探数据分析处理,涉及动校拉伸畸变处理,具体地说是压制动校拉伸畸变的方法。
背景技术:
1、随着地震勘探技术的进一步发展,大偏移距地震数据采集已成为新常态,此类数据在后续处理时会出现严重的动校拉伸畸变问题。消除动校拉伸畸变对改善地震数据的分辨率,提高地震数据的保真度和浅层成像质量具有十分重要的意义。目前,通常采用动校道集切除、高阶动校正及动校正处理等方法消除动校拉伸畸变造成的影响,这些方法在一定程度上消除动校拉伸畸变的同时也有着各自的缺陷。
2、通过动校道集切除消除动校拉伸畸变,易导致浅层有效反射信息的丢失,不利于资料浅中层成像。再者,基于高阶动校正的动校拉伸畸变消除方法,采用了更多次的二项展开式,在获得较为精确的动校速度的同时,虽然能克服部分动校拉伸引起的成像畸变,但由于动校拉伸的固有特性,该方法也仅能减弱大偏移距地震数据采集所带来的动校时差高次项误差,并不能真正从本质上压制拉伸畸变。
3、常规的动校正处理会造成大偏移距数据波形畸变,降低频带,从而降低对薄互层的分辨力以及后续属性分析的精度,给高精度地震勘探带来不利影响。在动校正处理过程中,会发生子波拉伸,使子波频率降低,如何通过子波整形处理压制动校拉伸畸变也成为待解决的问题。
技术实现思路
1、本发明的目的,是要提供一种压制动校拉伸畸变的方法,通过子波提取、反褶积处理及整形算子应用,来抑制远偏移距数据的动校拉伸畸变,达到压制动校拉伸带来的成像畸变,从而达到提高地震资料成像品质和可靠性,提高中浅层资料利用率的目的。
2、为实现上述目的,本发明所采用的技术方案如下:
3、一种压制动校拉伸畸变的方法,它是对地震资料进行叠前数据预处理及共中心点域数据重排后,利用井震结合进行近远偏移距数据处理,提取得目标子波s1(t)和待整形子波s2(t);
4、对目标子波s1(t)和待整形子波s2(t)进行反褶积处理,求得滤波算子f(t);
5、将滤波算子f(t)通过褶积处理应用于叠前道集,压制动校拉伸畸变。
6、作为一种限定,所述叠前预处理数据为:获取地震资料,经叠前去噪、振幅补偿和反褶积处理得预处理数据;
7、所述共中心点域数据重排为:将预处理数据进行共中心点域数据重排,得共中心点道集数据。
8、作为进一步限定,所述近远偏移距数据处理中,近偏移距数据处理的方法为:
9、从炮检点偏移距数据中,筛选出共中心点道集数据中的近偏移距数据,经动校正和叠加处理,得近道叠加数据,再与钻测井合成记录或vsp走廊叠加结果进行波形吻合度对比:
10、若波形吻合度对比结果为吻合度高,则输出近偏移距叠加数据,从近偏移距叠加数据中提取出近道自相关子波,作为目标子波s1(t);
11、若波形吻合度对比结果为吻合度低,则重新进行数据筛选。
12、其中,若波形吻合度对比结果中,波形、频率及相位一致性高或地震数据与钻测井数据两者波形相似系数达到0.85及以上,则认为吻合度高。可根据本领域技术人员常识,依实际情况进行调整。
13、作为再进一步限定,所述近远偏移距数据处理中,远偏移距数据处理的方法为:
14、从炮检点偏移距数据中,筛选出共中心点道集数据中已动校拉伸畸变的远偏移距数据,经动校正和叠加处理,输出远偏移距叠加数据,从远偏移距叠加数据中提取出远道自相关子波,作为待整形子波s2(t)。
15、所述待整形子波s2(t)为已动校拉伸畸变的远偏移距叠加数据子波。
16、其中,已动校拉伸畸变的远偏移距数据是指的数据,其中δf为同一反射同相轴远偏移距共中心点道集主频与近偏移距共中心点道集主频的变化量,f为反射同相轴近偏移距共中心点道集主频。可根据本领域技术人员常识,依实际情况进行调整。
17、所述反褶积处理是根据最小二乘准则和超定方程求解滤波算子f(t)。
18、其中,反褶积处理过程如下:
19、(1)确定关于滤波算子f(t)的线性方程组
20、假设滤波算子f(t)作用后输出的f(t)*s2(t)与目标子波s1(t)的残差为e(t),则有
21、e(t)=f(t)*s2(t)-s1(t)。
22、在实际计算中,滤波算子随着时间t的无限增大而逐渐趋于0,滤波算子具有有限能量,为可实现的滤波算子。
23、根据最小二乘准则,令误差函数求解滤波算子f(t),利用极限的必要条件得到关于f(t)的线性方程组,其中n,m为滤波算子的长度:
24、
25、(2)通过解超定方程求得滤波算子f(t)
26、首先,解超定方程rtra=rtx,其中
27、
28、目标子波矩阵
29、滤波算子矩阵
30、可得到,
31、
32、根据自相关的定义:x(t)为在一个时间序列的不同时间上相似性的度量。
33、用下式表示
34、
35、则,
36、以此类推,得出:
37、
38、同时x(t)和y(t)的互相关定义为
39、
40、则同理从而求得整形算子f(t)(其中n=m为滤波算子长度)
41、由于采用了上述的技术方案,本发明与现有技术相比,所取得的技术进步在于:
42、本发明的方法通过子波提取、反褶积处理及整形算子应用,排除最小偏移距或零偏移距共中心点道集子波提取的随机性,并通过与钻测井资料波形吻合度对比进行质控,实现井控保真处理,通过抑制远偏移距地震数据的动校拉伸畸变,达到压制动校拉伸带来的成像畸变,从而提高地震资料成像品质和可靠性,同时提高中浅层地震资料利用率。
43、本发明中目标子波的求取是通过叠前去噪、振幅补偿、反褶积等一系列地震资料预处理工作后所得,排除了地震子波在传播过程中受干扰波以及传播路径中能量、频率衰减导致的畸变,尽量恢复地震激发子波原始状态。
44、本发明可应用于石油地震勘探数据分析处理技术领域,用于压制动校拉伸畸变,提高中浅层资料成像品质,降低钻探失误风险,为油气田高效勘探开发提供可靠数据。
45、下面结合附图及具体实施例对本发明作更进一步详细说明。
1.一种压制动校拉伸畸变的方法,其特征在于,它是对地震资料进行叠前数据预处理及共中心点域数据重排后,利用井震结合进行近远偏移距数据处理,提取得目标子波s1(t)和待整形子波s2(t);
2.根据权利要求1所述的压制动校拉伸畸变的方法,其特征在于,所述叠前预处理数据为:获取地震资料,经叠前去噪、振幅补偿和反褶积处理得预处理数据;
3.根据权利要求2所述的压制动校拉伸畸变的方法,其特征在于,所述近远偏移距数据处理中,近偏移距数据处理的方法为:
4.根据权利要求2所述的压制动校拉伸畸变的方法,其特征在于,所述近远偏移距数据处理中,远偏移距数据处理的方法为:
5.根据权利要求1或2所述的压制动校拉伸畸变的方法,其特征在于,所述反褶积处理是根据最小二乘准则和超定方程求解滤波算子f(t)。