裂缝-溶蚀孔洞型碳酸盐岩储层高渗透带预测方法及装置与流程

本发明涉及地球物理勘探，特别涉及一种裂缝-溶蚀孔洞型碳酸盐岩储层高渗透带预测方法及装置。

背景技术：

1、碳酸盐岩储层，尤其是裂缝-溶蚀孔洞型碳酸盐岩储层，具有非常强的非均质性，这些非均质性对储层的渗透性具有重要影响，研究高渗透带的分布对油田开发具有重要的意义。

2、渗透率在储层研究中是最重要的参数，特别是针对强非均质碳酸盐岩储层来说，其空间分布一直都是该类储层研究的重点和难点，也是预测高渗透带分布的关键。渗透率通常具有两种，一种是基质的渗透率，另一种是裂缝的渗透率，裂缝渗透率一般是使用dfn(离散裂缝网络模型)通过工程算法进行模拟，高渗透带是这二者的综合影响结果。对于裂缝性储层，基质的渗透性较低，通常裂缝渗透性具有决定性作用。对于溶蚀孔洞发育的储层，基质的渗透性有时会很高，甚至于超出裂缝对渗透性的贡献，所以，对溶蚀孔洞发育的储层而言，二者都非常重要，不能忽视任何一方。

3、溶蚀孔洞型碳酸盐岩的渗透率具有非常强的非均质性，使用常规的渗透率变异系数(c)、突变系数(tk)和级差(jk)等参数难以评价高渗透带的分布。目前还没有有效方法对渗透率空间分布进行估计，插值方法对溶蚀孔洞型碳酸盐岩基本上不适用。并且，裂缝渗透率难以估计，使用数值模拟算法得到的渗透性往往具有较大的不确定性。

4、研究表明，岩心物性分析渗透率被视为最准确的基质渗透率测量结果；生产数据可以定性地反映裂缝性质。此外，地震属性数据可与相应的渗透率结合和反演，对渗透率进行预测。然而，现有技术中缺少这种方法。

技术实现思路

1、为了解决现有技术存在问题，本发明提供了一种裂缝-溶蚀孔洞型碳酸盐岩储层高渗透带预测方法，所述方法包括：

2、根据裂缝-溶蚀孔洞型碳酸盐岩储层的振幅地震属性和岩心渗透率获取基质渗透率体；

3、根据裂缝-溶蚀孔洞型碳酸盐岩储层的分频地震属性和试井渗透率获取裂缝渗透率体；

4、根据基质渗透率体和裂缝渗透率体融合确定裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体；

5、划分裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体的高渗透区域，预测裂缝-溶蚀孔洞型碳酸盐岩储层的高渗透带。

6、进一步地，根据裂缝-溶蚀孔洞型碳酸盐岩储层的振幅地震属性和岩心渗透率获取基质渗透率体包括：

7、根据裂缝-溶蚀孔洞型碳酸盐岩岩心分析，获取岩心孔隙度和渗透率；

8、将岩心渗透率分别进行对数化和bootstrap方法样本抽样，获取服从正态分布的岩心渗透率，建立服从正态化分布的岩心渗透率和岩心孔隙度的关系；

9、将裂缝-溶蚀孔洞型碳酸盐岩储层的振幅地震属性归一化，建立归一化后的裂缝-溶蚀孔洞型碳酸盐岩储层的振幅地震属性和岩心孔隙度的关系；

10、根据服从正态分布的岩心渗透率与岩心孔隙度的关系和归一化后的裂缝-溶蚀孔洞型碳酸盐岩储层的振幅地震属性和岩心孔隙度的关系，将归一化后的裂缝-溶蚀孔洞型碳酸盐岩储层的振幅地震属性转化为基质渗透率，获取基质渗透率体。

11、进一步地，根据裂缝-溶蚀孔洞型碳酸盐岩储层的分频地震属性和试井渗透率获取裂缝渗透率体包括：

12、将裂缝-溶蚀孔洞型碳酸盐岩储层试井渗透率分别经过对数化和bootstrap方法样本抽样，获取服从正态分布的裂缝-溶蚀孔洞型碳酸盐岩储层试井渗透率；

13、对裂缝-溶蚀孔洞型碳酸盐岩储层的地震数据进行广义s变换，将时间域的地震数据转化为频率域，确定品质因子属性，建立分频地震属性；

14、将裂缝-溶蚀孔洞型碳酸盐岩储层的分频地震属性归一化，获取归一化后的分频地震属性；

15、根据服从正态分布的裂缝-溶蚀孔洞型碳酸盐岩储层试井渗透率，将归一化后的分频地震属性进行分布转换，获取裂缝渗透率体。

16、进一步地，对裂缝-溶蚀孔洞型碳酸盐岩储层的地震数据进行广义s变换，将时间域的地震数据转化为频率域，确定品质因子属性，建立分频地震属性包括：

17、将初始振幅谱a0(f)表示为，

18、|a1(f)|＝|a0(f)|e-πft/q

19、其中，t为从始发到目标的双程旅行时，a1(f)为目标反射界面的振幅谱，f为频率，q为品质因子；

20、

21、其中，即为裂缝-溶蚀孔洞型碳酸盐岩储层的分频地震属性，t0为初始时间，t1为到达目标反射界面的时间。

22、进一步地，根据服从正态分布的裂缝-溶蚀孔洞型碳酸盐岩储层试井渗透率，将归一化后的分频地震属性进行分布转换，获取裂缝渗透率体，包括，

23、根据服从正态分布的裂缝-溶蚀孔洞型碳酸盐岩储层试井渗透率获取服从正态分布的裂缝-溶蚀孔洞型碳酸盐岩储层试井渗透率数据的累积概率分布；

24、根据归一化后的分频地震属性数据获取归一化后的分频地震属性数据的累积概率分布；

25、根据服从正态分布的裂缝-溶蚀孔洞型碳酸盐岩储层试井渗透率数据的累积概率分布和归一化后的分频地震属性数据的累积概率分布进行分位数逆变换，获取裂缝渗透率体：

26、y＝fy-1[fz(z)]

27、其中，y为裂缝渗透率体，fy为归一化后的分频地震属性数据的累积概率分布，fz(z)为服从正态分布的裂缝-溶蚀孔洞型碳酸盐岩储层试井渗透率数据的累积概率分布。

28、进一步地，根据基质渗透率体和裂缝渗透率体融合确定裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体之后，还包括提高裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体的连续性：

29、计算裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体的均值和方差；

30、剔除落在95％统计检验的置信区间之外的裂缝-溶蚀孔洞型碳酸盐岩储层总渗透率体，获取结构性连续性裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体。

31、进一步地，划分裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体的高渗透区域，预测裂缝-溶蚀孔洞型碳酸盐岩储层的高渗透带包括：

32、划分出结构性连续性裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体大于500md的区域，对裂缝-溶蚀孔洞型碳酸盐岩储层的高渗透带进行预测。

33、本发明的另一方面还提供了一种裂缝-溶蚀孔洞型碳酸盐岩储层高渗透带预测装置，所述装置包括：

34、基质渗透率体获取单元，用于根据裂缝-溶蚀孔洞型碳酸盐岩储层的振幅地震属性和岩心渗透率获取基质渗透率体；

35、裂缝渗透率体获取单元，用于根据裂缝-溶蚀孔洞型碳酸盐岩储层的分频地震属性和试井渗透率获取裂缝渗透率体；

36、总渗透率体确定单元，用于根据基质渗透率体和裂缝渗透率体融合确定裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体；

37、高渗透带预测单元，用于划分裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体的高渗透区域，预测裂缝-溶蚀孔洞型碳酸盐岩储层的高渗透带。

38、进一步地，所述总渗透率体确定单元用于根据基质渗透率体和裂缝渗透率体融合确定裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体之后，还用于提高裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体的连续性：

39、计算裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体的均值和方差；

40、剔除落在95％统计检验的置信区间之外的裂缝-溶蚀孔洞型碳酸盐岩储层总渗透率体，获取结构性连续性裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体。

41、进一步地，所述高渗透带预测单元用于划分裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体的高渗透区域，预测裂缝-溶蚀孔洞型碳酸盐岩储层的高渗透带包括：

42、划分出结构性连续性裂缝-溶蚀孔洞型碳酸盐岩储层的总渗透率体大于500md的区域，预测裂缝-溶蚀孔洞型碳酸盐岩储层的高渗透带。

43、相对于现有技术，本发明具有以下的有益效果：

44、本发明提供的裂缝-溶蚀孔洞型碳酸盐岩储层高渗透带预测方法，可以有效的预测油气藏开发的高渗透带，为油气藏开发提供有意义的目标区，降低油藏模拟过程中的不确定性，对于开发布井具有重要的指导作用。

45、本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书以及附图中所指出的步骤或流程来实现和获得。