一种油气藏孔隙度分布图生成方法、装置、设备、介质及程序与流程

本发明涉及油气田勘探，尤其涉及一种油气藏孔隙度分布图生成方法、装置、电子设备、介质及程序。

背景技术：

1、断控缝洞型油气藏的储层发育受构造运动、岩溶过程的影响很大，储集体空间分布特征复杂，储集体内的孔隙度非均质极强。孔隙度是油气藏地质储量估算的关键参数，地质储量是开发方案设计的基础。因此，研究一种适用于断控缝洞型油气藏的孔隙度计算方法对此类油气藏的有效开发具有重要的意义。

2、研究早期，大型洞穴的孔隙度根据产量情况赋予经验值；孔洞的孔隙度参考测井解释的结果；裂缝的孔隙度基于统计的裂缝开度，利用测井解释或离散裂缝模型计算。这种处理方法受研究人员主观性的影响较强，孔隙度模拟结果的误差较大；另外，研究人员建立了孔隙度-纵波阻抗的关系图版，拟合了孔隙度-纵波阻抗关系曲线，并基于拟合的关系曲线通过纵波阻抗计算得到孔隙度结果。由于孔隙度与纵波阻抗关系在二维坐标系中呈现出较宽的条带状，单一的拟合曲线难以准确地表征这样的关系，忽略了储层孔隙结构存在的非均质性特征，导致孔隙度的分析结果准确性较低，因此，如何提高油气藏孔隙度分布图生成的准确性成了亟待解决的问题。

技术实现思路

1、针对上述问题，本发明的实施例提供了一种油气藏孔隙度分布图生成方法、装置、电子设备、介质及程序。

2、第一方面，本发明实施例提供了一种油气藏孔隙度分布图生成方法，包括：

3、根据预设的断控缝洞储层类型对目标油气藏进行分类表征，得到所述目标油气藏的储集体轮廓模型；

4、获取所述目标油气藏的孔隙度数据，利用所述储集体轮廓模型和所述孔隙度数据生成所述目标油气藏的波阻抗约束数据；

5、根据绘制概率图法对所述孔隙度数据进行曲线拟合，得到所述目标油气藏的拟合参数，利用所述拟合参数和所述波阻抗约束数据构建所述目标油气藏的孔隙度分布模型；

6、对所述孔隙度分布模型进行参数优化，得到优化后的孔隙度分布模型，利用所述优化后的孔隙度分布模型生成所述目标油气藏的孔隙度分布图。

7、根据本发明的实施例，所述根据预设的断控缝洞储层类型对目标油气藏进行分类表征，得到所述目标油气藏的储集体轮廓模型，包括：

8、获取目标油气藏的目标地球物理属性，确定所述目标地球物理属性的属性阈值；

9、根据所述属性阈值和预设的断控缝洞储层类型对目标油气藏进行表征，得到所述目标油气藏的分类表征；

10、根据预设的断控缝洞储层类型对所述分类表征进行特征融合，得到所述目标油气藏的储集体轮廓模型。

11、根据本发明的实施例，所述利用所述储集体轮廓模型和所述孔隙度数据生成所述目标油气藏的波阻抗约束数据，包括：

12、根据所述储集体轮廓模型对所述目标油气藏进行纵波阻抗反演，得到所述目标油气藏的阻抗分布；

13、对所述阻抗分布进行正态得分变换，得到所述阻抗分布的正态阻抗分布；

14、利用所述正态阻抗分布、所述孔隙度数据和预设的同位协同克里金算法生成所述目标油气藏的波阻抗约束数据，其中，所述预设的同位协同克里金算法为：

15、

16、其中，是位置u处的待估值点的孔隙度估值误差，α是参与估值点的标识，n是参与估值点的总数，λα是第α个参与估值点的孔隙度误差权重系数，y(uα)是参与估值点孔隙度误差，λ′是待估值点的波阻抗误差权重系数，y2(u)是待估值点的波阻抗误差，u是位置标识，uα是第α个参与估值点的位置标识。

17、根据本发明的实施例，所述根据绘制概率图法对所述孔隙度数据进行曲线拟合，得到所述目标油气藏的拟合参数，包括：

18、根据所述预设的断控缝洞储层类型对所述孔隙度数据进行数据分类，得到所述孔隙度数据的分类孔隙度；

19、对所述分类孔隙度进行升序排列，得到所述分类孔隙度中孔隙度数据的数据标记；

20、利用所述数据标记和预设的分位数算法生成所述分类孔隙度的分位数，其中，所述预设的分位数算法为：

21、

22、其中，qi是所述分类孔隙度的分位数，i是所述孔隙度数据的数据标记，n是所述分类孔隙度中孔隙度数据的数据总数；

23、根据所述分位数生成所述分类孔隙度的曲线参数，确定所述曲线参数为所述目标油气藏的拟合参数。

24、根据本发明的实施例，所述根据所述分位数生成所述分类孔隙度的曲线参数，包括：

25、生成所述分位数的累积分布函数反函数，利用预设的线性表征函数和所述累积分布函数反函数确定所述分类孔隙度的曲线参数，其中，所述预设的线性表征函数为：

26、x＝μ+σg-1(q)

27、其中，x是分类孔隙度的孔隙度数据集合，μ是线性表征函数的截距，σ是线性表征函数的斜率，g-1(q)是所述分位数集合的累积分布函数反函数，q是所述分类孔隙度的分位数集合。

28、根据本发明的实施例，所述对所述孔隙度分布模型进行参数优化，得到优化后的孔隙度分布模型，包括：

29、生成所述孔隙度分布模型的正态拟合误差和权重系数误差；

30、根据所述正态拟合误差和所述权重系数误差对所述孔隙度分布模型进行参数优化，得到优化后的孔隙度分布模型。

31、第二方面，本发明实施例提供了一种油气藏孔隙度分布图生成装置，其特征在于，包括：

32、分类表征模块，用于根据预设的断控缝洞储层类型对目标油气藏进行分类表征，得到所述目标油气藏的储集体轮廓模型；

33、波阻抗约束模块，用于获取所述目标油气藏的孔隙度数据，利用所述储集体轮廓模型和所述孔隙度数据生成所述目标油气藏的波阻抗约束数据；

34、拟合参数模块，用于根据绘制概率图法对所述孔隙度数据进行曲线拟合，得到所述目标油气藏的拟合参数，利用所述拟合参数和所述波阻抗约束数据构建所述目标油气藏的孔隙度分布模型；

35、参数优化模块，用于对所述孔隙度分布模型进行参数优化，得到优化后的孔隙度分布模型，利用所述优化后的孔隙度分布模型生成所述目标油气藏的孔隙度分布图。

36、第三方面，本发明实施例提供了一种电子设备，其包括：

37、处理器；

38、用于存储所述处理器可执行指令的存储器；

39、其中，所述处理器被配置为执行所述指令，以实现如前第一方面所述的一种油气藏孔隙度分布图生成方法。

40、第四方面，本发明实施例提供了一种介质，其上存储有计算机程序，该程序被被处理器执行时，实现如前第一方面所述的一种油气藏孔隙度分布图生成方法。

41、与现有技术相比，本发明的上述技术方案具有如下有益效果：

42、本发明的实施例利用同位协同克里金算法生成所述目标油气藏的波阻抗约束数据，并利用绘制概率图法对所述孔隙度数据进行曲线拟合，得到所述目标油气藏的拟合参数，再采用高斯随机模拟算法将波阻抗约束数据和拟合参数模拟孔隙度三维空间分布，充分挖掘了测井、地震等数据中包含的孔隙度分布及空间趋势信息，可构建较为准确的孔隙度模型，为断控缝洞型油气藏的储量估算及储量不确定性分析提供依据，因此本发明提出的油气藏孔隙度分布图生成方法、装置、设备、介质及程序，可以解决生成油气藏孔隙度分布图的准确性较低的问题。