【技术实现步骤摘要】
本公开涉及地震勘探数据解释性处理，尤其涉及一种非均质储层叠后地震精细刻画方法及装置。
技术介绍
1、在漫长的地质年代过程中，由于受到构造力、重力及其他力的作用，地层中坚硬岩石或脆性较高的地层例如碳酸盐岩通常极易产生裂缝等，造成地层的非均质性空间变化，形成非均质性地震响应。据国内外文献报道，在致密砂岩、泥页岩、碳酸盐岩、火山岩等岩石的裂缝型非均质储层中已获得高产工业油气流，主要原因在于地层中的裂缝体系可以实现孔隙型储层的沟通，从而形成油气排出途径，提供了有效的渗流通道和储集空间，增加了油气井的泄油面积，提高了钻井的油气产能，因此，寻找储层段中裂缝发育部位并确定有利裂缝发育带是油气勘探工作的重要组成部分。现有的裂缝预测方法包括基于叠前地震资料的裂缝预测方法以及基于叠后地震资料的裂缝预测方法。其中，基于叠前地震资料的裂缝地震预测技术和流程复杂，受地表条件、采集纵横比、速度精度、各向异性、资料信噪比等控制因素多，勘探成本高；而基于叠后地震资料的裂缝预测技术流程简单、快捷、稳定，可以有效克服基于叠前地震资料裂缝预测技术的上述缺陷，对于窄方位资料裂缝识别具有显著的技术优势，但因基于叠后地震资料的裂缝预测过程中，通过滤除低频，保留高频数据可以突出地层细节，处理结果保幅性差，存在逐道处理痕迹引入裂缝假象，从而不能够从全频带的地震数据体中充分挖掘反映非均质性储层地震响应问题，无法满足裂缝型均质油气藏精细刻画需求。
技术实现思路
1、本公开提出了一种非均质储层叠后地震精细刻画方法及装置，以解决以往针对
2、根据本公开的一方面，提供了一种非均质储层叠后地震精细刻画方法，包括：
3、获取工区叠前时间偏移或叠前深度偏移后输出的成像道集；
4、对成像道集进行精细切除及加权同相叠加处理，得到叠加地震数据体；
5、采用迭代滤波法，对叠加地震数据体进行迭代滤波处理，当迭代滤波次数达到预定范围次数以内，叠加地震数据体的骨架特征及信噪比满足预定要求后，完成迭代滤波，得到反射地震数据体；
6、利用反射地震数据体，在迭代滤波前的叠加地震数据体中分离出散射地震数据体；
7、基于散射地震数据体，确定非均质储层的空间展布规律以及/或剖面分布规律，完成非均质储层叠后地震精细刻画。
8、优选地，所述获取的成像道集，为保幅叠前时间偏移或叠前深度偏移后输出的成像道集；
9、其中，所述成像道集为共成像点偏移距道集或共成像点角道集。
10、优选地，所述对成像道集进行精细切除及加权同相叠加处理之前，还包括：
11、对成像道集进行数据净化处理，其中，所述数据净化处理，包括；
12、若所述成像道集存在特定频带的噪声能量，则对成像道集进行带通滤波处理；
13、若所述成像道集存在多次波能量干扰时，对成像道集中的多次波干扰进行保幅压制处理；
14、若所述成像道集的信噪比不等于预定信噪比时，对所述成像道集进行随机噪声压制处理。
15、优选地，在所述采用迭代滤波法，对叠加地震数据体进行迭代滤波处理之前，还包括：
16、若所述叠加地震数据体空间分布不符合预定要求，则对叠加地震数据体进行规则化处理，其中，所述规则化处理的方法，包括：
17、选择包含所述叠加地震数据体的最大规则化范围；
18、在所述最大规则化范围内，对叠加地震数据体之外的区域的每个面元位置的每个采样点值进行充零处理，得到规则化后的叠加地震数据体。
19、优选地，所述采用迭代滤波法，对叠加地震数据体进行迭代滤波处理的方法，包括：
20、若所述叠加地震数据体为二维地震数据，则利用二维地震工区数据迭代滤波处理公式，对所述叠加地震数据体进行迭代滤波处理；
21、若所述叠加地震数据体为三维地震数据，则利用三维地震工区数据迭代滤波处理公式，对所述叠加地震数据体进行迭代滤波处理；
22、其中，在每一次迭代滤波完成后，对得到的叠加地震数据体进行一次边界镶边处理。
23、优选地，所述二维地震工区数据迭代滤波处理公式，包括：
24、
25、式中：α为迭代滤波系数，t为当前迭代滤波次数，m、n分别代表x、z两个方向的空间位置序号，m＝1，2，……m；n＝1,2，……，n，u0为初始叠加地震数据体，ut为第t次迭代后的叠加地震数据体，ut+1为第t+1次迭代后的叠加地震数据体。
26、优选地，所述三维地震工区数据迭代滤波处理公式，包括：
27、
28、式中，α为迭代滤波系数，t为当前迭代滤波次数，m、n、k分别代表x、z、y三个方向的空间位置序号，m＝1，2，……m；n＝1,2，……，n；k＝1，2，……k，u0为初始叠加地震数据体，ut为第t次迭代后的叠加地震数据体，ut+1为第t+1次迭代后的叠加地震数据体。
29、优选地，所述边界镶边处理的方法，包括：
30、若所述叠加地震数据体为二维地震数据，则利用二维地震工区数据边界处理公式，对迭代滤波后的叠加地震数据体进行镶边处理；
31、其中，所述二维地震工区数据边界处理公式，包括：二维边处理公式以及二维角点处理公式；
32、其中，所述二维边处理公式为：
33、
34、式中：m、n分别为x、z两个方向的空间位置序号，m、n分别为x、z两个方向的最大空间位置序号，t为当前迭代滤波次数；
35、其中，所述二维角点处理公式为：
36、
37、式中：m、n分别为x方向和z方向的最大空间位置序号,t为当前迭代滤波次数；
38、以及/或；
39、若所述叠加地震数据体为三维地震数据，则利用三维地震工区数据边界处理公式，对迭代滤波后的叠加地震数据体进行镶边处理；
40、其中，所述三维地震工区数据边界处理公式，包括：三维面处理公式、三维棱边处理公式以及三维角点处理公式；
41、其中，所述三维面处理公式为：
42、
43、
44、式中：m、n、k分别为x、z、y三个方向的空间位置序号，m、n、k分别为x、z、y三个方向的最大空间位置序号，t为当前迭代滤波次数；
45、其中，所述三维棱边处理公式为：
46、
47、
48、
49、
50、式中：m、n、k分别为x、z、y三个方向的空间位置序号，m、n、k分别为x、z、y三个方向的最大空间位置序号，t为当前迭代滤波次数；
51、其中，所示三维角点处理公式为：
52、
53、
54、式中：m、n、k分别为x、z、y三个方向的最大空间位置序号，t本文档来自技高网...
【技术保护点】
1.一种非均质储层叠后地震精细刻画方法，其特征在于，包括：
2.根据权利要求1所述的非均质储层叠后地震精细刻画方法，其特征在于：
3.根据权利要求1所述的非均质储层叠后地震精细刻画方法，其特征在于，所述对成像道集进行精细切除及加权同相叠加处理之前，还包括：
4.根据权利要求1所述的非均质储层叠后地震精细刻画方法，其特征在于，在所述采用迭代滤波法，对叠加地震数据体进行迭代滤波处理之前，还包括：
5.根据权利要求1所述的非均质储层叠后地震精细刻画方法，其特征在于，所述采用迭代滤波法，对叠加地震数据体进行迭代滤波处理的方法，包括：
6.根据权利要求5所述的非均质储层叠后地震精细刻画方法，其特征在于，所述二维地震工区数据迭代滤波处理公式，包括：
7.根据权利要求5所述的非均质储层叠后地震精细刻画方法，其特征在于，所述三维地震工区数据迭代滤波处理公式，包括：
8.根据权利要求5所述的非均质储层叠后地震精细刻画方法，其特征在于，所述边界镶边处理的方法，包括：
9.根据权利要求1-8任一项所述的非均
10.一种非均质储层叠后地震精细刻画装置，其特征在于，包括：
...【技术特征摘要】
1.一种非均质储层叠后地震精细刻画方法，其特征在于，包括：
2.根据权利要求1所述的非均质储层叠后地震精细刻画方法，其特征在于：
3.根据权利要求1所述的非均质储层叠后地震精细刻画方法，其特征在于，所述对成像道集进行精细切除及加权同相叠加处理之前，还包括：
4.根据权利要求1所述的非均质储层叠后地震精细刻画方法，其特征在于，在所述采用迭代滤波法，对叠加地震数据体进行迭代滤波处理之前，还包括：
5.根据权利要求1所述的非均质储层叠后地震精细刻画方法，其特征在于，所述采用迭代滤波法，对叠加地震数据体进行迭代滤波处理的方法，包括：
6.根...
【专利技术属性】
技术研发人员：陈可洋，潘文庆，杨微，张元高，李强，陆加敏，刘洋，孙立东，戴世立，鞠林波，
申请(专利权)人：大庆油田有限责任公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。