一种基于约束稀疏脉冲反演的砂泥岩预测方法及装置与流程

本发明涉及油气田勘探，特别涉及一种基于约束稀疏脉冲反演的砂泥岩预测方法及装置。

背景技术：

1、油气开发从最初的粗犷开发逐步转化为精细开发模式，水平井多分支水平井开发、注水注气开发方式应运而生，随之而来的是对地下地质情况的精细油藏描述要求，弄清地下储层的分布细节是精细油藏描述与精准开发的基础，储层预在油藏地质研究中的作用愈发重要。

2、目前油藏储层的刻画主要通以下3种技术方案来解决：精细地层对比，通过邻井精细地层对比，确定砂岩储层的展布方向和分布规模，但是很难精细刻画某一层砂体的展布细节、延伸长短和联通关系。约束稀疏脉冲反演，敏感性分析可行的前提下，开展精细井震标定，提取地震综合子波，融合低频模型，进行地震数据的反褶积运算，得到地层绝对波阻抗数据，进而确定砂泥岩分布，此方法对于薄储层分辨能力有限。地质统计学反演，该反演方法分辨率取决了网格的大小，但该方法具有一定的随机性，不确定性较大，预测准确性有待进一步提高。

技术实现思路

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、根据所述砂泥岩储层的纵波阻抗均方根振幅值，确定砂岩的纵波阻抗均方根振幅值和泥岩的纵波阻抗均方根振幅值；

26、根据砂岩的纵波阻抗均方根振幅值和泥岩的纵波阻抗均方根振幅值，确定砂岩与泥岩的分布概率。

27、在一个实施例中，所述根据所述地震数据综合子波和砂岩与泥岩的分布概率，进行地质统计学反演，包括：

28、根据所述所有单井在各个层位的反射系数和地震数据、地震数据综合子波和砂泥岩分布概率，进行地质统计学反演，得到统计反演结果；

29、根据新增单井在各个层位的反射系数和地震数据，对所述统计反演结果进行验证，若所述统计反演结果与所述新增单井在各个层位的反射系数和地震数据的偏差率大于预设的最大偏差率，则调整所述稀疏反演参数，重新进行稀疏反演和统计反演，并验证，直至反演结果的数据与实际井数据一致的比例大于预设的阈值。

30、第二方面，本发明实施例提供了一种约束基于稀疏脉冲反演的砂泥岩预测装置，包括：

31、采集模块，用于采集砂泥岩储层的地震数据；

32、敏感模块，用于根据所述地震数据，确定砂泥岩储层的纵波阻抗对砂泥岩储层区分的敏感效应；

33、子波模块，用于进行井震精细标定，确定地震数据综合子波；

34、模型模块，用于建立时间域构造格架模型和纵波阻抗模型；

35、稀疏模块，用于根据所述地震数据和地震数据综合子波，确定约束稀疏脉冲反演的约束稀疏脉冲反演参数，进行约束稀疏脉冲反演，得到砂岩与泥岩的分布概率；

36、统计模块，用于根据所述地震数据综合子波和所述砂岩与泥岩的分布概率，进行地质统计学反演。

37、第三方面，本发明实施例提供了一种计算设备，包括：存储器、处理器及存储于存储器上并可在处理器上运行的计算机程序，所述处理器执行的所属程序时实现所述的一种基于约束稀疏脉冲反演的砂泥岩预测方法。

38、第四方面，本发明实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述的一种基于约束稀疏脉冲反演的砂泥岩预测方法。

39、本发明实施例提供的上述技术方案的有益效果至少包括：

40、本发明实施例提供了一种基于约束稀疏脉冲反演的砂泥岩预测方法，包括：采集砂泥岩储层的地震数据；根据地震数据，确定砂泥岩储层的纵波阻抗对砂泥岩储层区分的敏感效应；进行井震精细标定，确定地震数据综合子波；根据地震数据，建立时间域构造格架模型和纵波阻抗模型；根据地震数据和地震数据综合子波，确定约束稀疏脉冲反演的约束稀疏脉冲反演参数，进行约束稀疏脉冲反演，得到砂岩与泥岩的分布概率；根据地震数据综合子波和砂岩与泥岩的分布概率，进行地质统计学反演。本发明实施例在先进行敏感性分析，确定反演可行性的前提下，开展精细井震标定、提取地震数据综合子波，同时利用构造解释的层位、断层数据和地震数据建立纵波阻抗模型，开展约束稀疏脉冲反演，得到纵波阻抗均方根振幅值。然后基于纵波阻抗均方根振幅值，结合单井敏感性分析和已钻井认识，进行砂泥岩平面分布概率刻度。利用刻度好的砂泥岩平面分布概率分布作为约束，在地质统计学反演参数分析和确定的基础上，开展地质统计学反演，并利用多次实现取平均值压制随机性。通过该技术创新，实现了油藏储层平面规模和纵向分布的精细刻画，为油藏的精细描述、水平井部署纵向选层及注水开发的奠定了有力的基础。本发明实施例有效指导了研究区有利目标层位、目标区域优选、井位部署，进而实现了研究区的高效勘探开发，应用前景广阔。

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

42、下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。