一种断控缝洞型储集体的刻画方法与流程

本发明属于地震勘探，更具体地，涉及一种断控缝洞型储集体的刻画方法、电子设备及介质。

背景技术：

1、目前，断控缝洞型油气藏在全世界范围内并不常见，但该类油藏储存了大量的油气资源。相对于国外，我国断控缝洞型油气藏埋深大、后期改造作用复杂，储层非均质性极强，勘探开发难度较大。“溶蚀洞”以及“缝洞带”是该类油藏的重要储集空间，对于该类储集空间的识别，目前的方法主要包括三种：第一种为利用地震反射振幅、分频数据地震属性进行溶蚀带的刻画；第二种为基于稀疏脉冲或者地质统计学反演方法，利用反演得到的参数以及敏感参数门槛值进行溶洞雕刻；第三种为利用断溶体轮廓属性对反演所用到的低频模型进行约束重构，在此基础上，基于反演得到低阻体异常进行溶洞的刻画。

2、现有方法中，第一类属性计算的方法简单方便，但受到子波的影响，无法准确地反映地下溶蚀孔洞的真实发育位置。第二、三类方法中涉及的阻抗反演虽然可以降低子波旁瓣效应，恢复溶蚀洞核部发育位置，对溶蚀洞体积进行瘦身，但是在溶蚀程度较弱的断控型油气藏中，反演的低阻溶洞异常响应容易被背景反射淹没。目前，大部分的研究均聚焦在溶蚀洞以及断溶体轮廓带的预测，鲜有涉及“缝洞带”突显识别的研究。

3、因此，期待发明一种断控缝洞型储集体的刻画方法，能够解决现有技术中缝洞带的刻画精度较低的问题。

技术实现思路

1、本发明的目的是提出一种断控缝洞型储集体的刻画方法，以解决现有技术中缝洞带的刻画精度较低的问题。

2、为了实现上述目的，本发明提供一种断控缝洞型储集体的刻画方法，包括：

3、步骤1：获取三维地震反射振幅数据d，并基于所述三维地震反射振幅数据d，得到三维相对阻抗数据；

4、步骤2：基于所述三维相对阻抗数据，获取井点处的相对阻抗伪井曲线，并基于所述相对阻抗伪井曲线和纵波阻抗测井曲线中的中-低频数据，得到纵波阻抗与相对阻抗之间的转换关系公式，进而求取绝对阻抗数据ip；

5、步骤3：基于所述绝对阻抗数据ip和所述三维地震反射振幅数据d，在构造层位三维趋势控制下，选取三维滑动窗口v，并根据选取的三维滑动窗口v，分别计算所述绝对阻抗数据ip的阻抗背景趋势ip_bg和所述三维地震反射振幅数据d的振幅背景趋势dbg；

6、步骤4：计算所述绝对阻抗数据ip和所述阻抗背景趋势ip_bg之间的差值，得到背景趋势压制后的阻抗数据ip_diff，并计算所述三维地震反射振幅数据d和所述振幅背景趋势dbg之间的差值，得到背景趋势压制后的地震振幅数据ddiff；

7、步骤5：基于所述背景趋势压制后的阻抗数据ip_diff，计算得到溶洞响应，基于所述背景趋势压制后的地震振幅数据ddiff，计算得到缝洞带响应。

8、可选地，所述步骤1包括：

9、步骤11：获取三维地震反射振幅数据d，以及测井速度数据或地震解释的速度数据；

10、步骤12：利用所述测井速度数据或所述地震解释的速度数据进行时间域和深度域的标定，得到时深标定结果；

11、步骤13：基于所述时深标定结果，估算所述目的层段对应的子波w；

12、步骤14：基于所述三维地震反射振幅数据d和所述子波w，求取相对阻抗，进而得到所述三维相对阻抗数据。

13、可选地，所述步骤2包括：

14、步骤21：基于所述三维相对阻抗数据，获取井点处的相对阻抗伪井曲线，并基于所述相对阻抗伪井曲线和纵波阻抗测井中的中-低频曲线，得到纵波阻抗与相对阻抗之间的转换关系公式，进而得到低频趋势模型ip_trend；

15、步骤22：根据所述低频趋势模型ip_trend，构建低频约束表达式trend；

16、步骤23：基于所述低频约束表达式trend，构建反演目标函数f(rp)；

17、步骤24：将所述目的层对应的子波矩阵w代入所述反演目标函数f(rp)中，得到纵波阻抗反射率rp；

18、步骤25：将所述纵波阻抗反射率rp代入以下公式中，得到绝对阻抗数据ip，

19、

20、其中，t为采样时间，ip1为低频趋势模型中起始时刻的纵波阻抗绝对值，rp(t)为采样时间为t时的纵波阻抗反射率。

21、可选地，所述步骤21包括：

22、基于所述三维相对阻抗数据，获取井点处的相对阻抗伪井曲线；

23、获取测井的原始纵波阻抗数据；

24、对所述测井的原始纵波阻抗数据依次进行低通滤波和重采样处理，得到所述井点处的实测纵波阻抗测井曲线；

25、将所述实测纵波阻抗测井曲线和所述相对阻抗伪井曲线进行交汇，得到所述纵波阻抗和所述相对阻抗之间的转换关系公式；

26、基于所述转换关系公式，将所述三维相对阻抗数据转化为所述低频趋势模型ip_trend。

27、可选地，所述低频约束表达式trend为

28、trend＝ωip||lipmrp-ip_trend||2，

29、其中，ωip为调节低频约束项权重的超参数，lip为低通滤波矩阵，m为积分矩阵，rp为纵波阻抗反射率；

30、所述反演目标函数f(rp)为

31、

32、其中，w为所述目的层对应的子波矩阵，为稀疏约束项，α为所述稀疏约束项的权重调节参数；

33、可选地，所述步骤3包括：

34、基于所述绝对阻抗数据ip和所述三维地震反射振幅数据d，沿层位构造趋势，在构造层位的发育趋势控制下，选取所述三维滑动窗口v，其中，所述三维滑动窗口v内包括i个采样点；

35、沿所述层位构造趋势，以任一个采样点为中心点，将所述三维滑动窗口v内的所述i个采样点进行重排，使得v1＜v2＜v3＜...＜vi，得到所述三维滑动窗口v内样本的中值(v1+vi)/2并将所述中值作为所述三维滑动窗口v的中心，其中，vi为窗口v内第i个采样点；

36、基于所述三维滑动窗口v内样本的中值，分别计算所述绝对阻抗数据ip的阻抗背景趋势ip_bg和所述三维地震反射振幅数据d的振幅背景趋势dbg。

37、可选地，

38、ip_diff＝ip-ip_bg，其中，ip表示所述绝对阻抗数据，ip_bg表示所述阻抗背景趋势，ip_diff表示所述背景趋势压制后的阻抗数据；

39、ddiff＝d-dbg，

40、其中，d表示所述三维地震反射振幅数据，dbg表示所述振幅背景趋势。

41、可选地，所述步骤5包括：

42、将所述背景趋势压制后的阻抗数据ip_diff代入所述溶蚀洞响应的增强计算公式中，得到溶洞响应ipe，所述溶蚀洞响应的增强计算公式为

43、

44、其中，e为自然指数，β为缩放因子且β的取值大于1；

45、将所述背景趋势压制后的地震振幅数据ddiff进行相位转换，得到归位后的缝洞带异常振幅ddiff_h；

46、基于所述归位后的缝洞带异常振幅ddiff_h，通过能量包络计算公式得到缝洞带响应env，所述能量包络的计算公式为

47、

48、其中，t为采样时间，ddiff_hf(t)为ddiff_h经过希尔伯特变换得到的实部，ddiff_hg(t)为ddiff_h经过希尔伯特变换得到的虚部。

49、一种电子设备，所述电子设备包括：

50、存储器，存储有可执行指令；

51、处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的断控缝洞型储集体的刻画方法。

52、一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的断控缝洞型储集体的刻画方法。

53、本发明的有益效果在于：

54、本发明首先获取三维地震反射振幅数据并求取三维相对阻抗数据，然后获取中-低频纵波阻抗测井曲线与相对阻抗伪井曲线之间的转换关系公式，并进行纵波阻抗反演，进而得到绝对纵波阻抗数据，再次基于反演得到的纵波绝对阻抗数据和三维地震反射振幅数据，分别求取其对应的背景趋势值，进行背景趋势干扰的压制处理，最后在背景趋势压制的基础上，进行溶蚀洞以及缝洞带响应增强优化处理，得到溶洞响应和缝洞带响应；本发明从反演出发，综合测井信息好和地震信息，在保持地下真实反射特征不受破坏的条件下，引入约束项，约束项的引入可使反演结果的空间非均质性得到有效提高，同时，本发明给出了背景趋势压制以及异常响应增强的优化方法，从而提高了溶蚀洞以及缝洞带的识别精度，再次，对实际研究区进行储集带三维刻画时，综合溶蚀洞以及缝洞带两种不同的相带共同约束，可以有效降低储集带预测的多解性，进一步提升地震数据对断控型油气藏的刻画能力。

55、本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。