储层孔隙度预测方法及装置.pdf

《储层孔隙度预测方法及装置.pdf》由会员分享，可在线阅读，更多相关《储层孔隙度预测方法及装置.pdf（20页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 202010452820.X (22)申请日 2020.05.26 (71)申请人 中国石油天然气股份有限公司 地址 100007 北京市东城区东直门北大街9 号 (72)发明人 王磊雷明徐中华杜炳毅 陈彬滔石兰亭方乐华史忠生 薛罗马轮史江龙 (74)专利代理机构 北京三友知识产权代理有限 公司 11127 代理人 赵平周永君 (51)Int.Cl. G01V 1/30(2006.01) G01V 1/28(2006.01) (54)发明名称 储层孔隙度预测方法及装置 (57)摘。

2、要 本申请实施例提供一种储层孔隙度预测方 法及装置， 方法包括： 根据目标工区的测井曲线 数据， 确定纵波阻抗与孔隙度的对应关系和横波 阻抗与孔隙度的对应关系； 根据所述纵波阻抗与 孔隙度的对应关系和所述横波阻抗与孔隙度的 对应关系， 确定泊松阻抗属性与孔隙度的对应关 系； 根据所述泊松阻抗属性与孔隙度的对应关系 和所述目标工区的三维地震数据体， 确定所述目 标工区的孔隙度分布； 本申请能够通过增加横波 约束降低孔隙流体对于孔隙度预测的不利影响， 提高储层孔隙度定量预测的精度。 权利要求书2页 说明书11页 附图6页 CN 111751879 A 2020.10.09 CN 111751879。

3、 A 1.一种储层孔隙度预测方法， 其特征在于， 所述方法包括： 根据目标工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关系和横波阻抗与孔隙 度的对应关系； 根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对应关系， 确定泊 松阻抗属性与孔隙度的对应关系； 根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工区的三维地震数据体， 确定 所述目标工区的孔隙度分布。 2.根据权利要求1所述的储层孔隙度预测方法， 其特征在于， 所述根据目标工区的测井 曲线数据， 确定纵波阻抗与孔隙度的对应关系和横波阻抗与孔隙度的对应关系， 包括： 根据测井曲线数据中的纵波速度、 横波速度以及地层密度， 确。

4、定纵波阻抗和横波阻抗； 根据预设岩石物理分析模型和预设拟合系数分别对所述纵波阻抗与孔隙度和所述横 波阻抗与孔隙度进行拟合分析， 确定所述纵波阻抗与所述孔隙度的对应关系和所述横波阻 抗与所述孔隙度的对应关系。 3.根据权利要求1所述的储层孔隙度预测方法， 其特征在于， 所述根据所述纵波阻抗与 孔隙度的对应关系和所述横波阻抗与孔隙度的对应关系， 确定泊松阻抗属性与孔隙度的对 应关系， 包括： 根据纵波阻抗与孔隙度的对应关系、 横波阻抗与孔隙度的对应关系以及与所述目标工 区对应的预设泊松阻抗调节因子， 确定所述目标工区中泊松阻抗属性与孔隙度的对应关 系。 4.根据权利要求1所述的储层孔隙度预测方法，。

5、 其特征在于， 所述根据所述泊松阻抗属 性与孔隙度的对应关系和所述目标工区的三维地震数据体， 确定所述目标工区的孔隙度分 布， 包括： 根据泊松阻抗属性与孔隙度的对应关系和所述目标工区的泊松阻抗属性地震数据体， 确定所述目标工区的孔隙度分布。 5.一种储层孔隙度预测装置， 其特征在于， 包括： 阻抗确定模块， 用于根据目标工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关 系和横波阻抗与孔隙度的对应关系； 预测模型构建模块， 用于根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔 隙度的对应关系， 确定泊松阻抗属性与孔隙度的对应关系； 孔隙度预测模块， 用于根据所述泊松阻抗属性与孔隙度的对应。

6、关系和所述目标工区的 三维地震数据体， 确定所述目标工区的孔隙度分布。 6.根据权利要求5所述的储层孔隙度预测装置， 其特征在于， 所述阻抗确定模块包括： 横纵阻抗确定单元， 用于根据测井曲线数据中的纵波速度、 横波速度以及地层密度， 确 定纵波阻抗和横波阻抗； 横纵阻抗孔隙度关系确定单元， 用于根据预设岩石物理分析模型和预设拟合系数分别 对所述纵波阻抗与孔隙度和所述横波阻抗与孔隙度进行拟合分析， 确定所述纵波阻抗与所 述孔隙度的对应关系和所述横波阻抗与所述孔隙度的对应关系。 7.根据权利要求5所述的储层孔隙度预测装置， 其特征在于， 所述预测模型构建模块包 括： 权利要求书 1/2 页 2 。

7、CN 111751879 A 2 泊松阻抗孔隙度关系确定单元， 用于根据纵波阻抗与孔隙度的对应关系、 横波阻抗与 孔隙度的对应关系以及与所述目标工区对应的预设泊松阻抗调节因子， 确定所述目标工区 中泊松阻抗属性与孔隙度的对应关系。 8.根据权利要求5所述的储层孔隙度预测装置， 其特征在于， 所述孔隙度预测模块包 括： 孔隙度分布确定单元， 用于根据泊松阻抗属性与孔隙度的对应关系和所述目标工区的 泊松阻抗属性地震数据体， 确定所述目标工区的孔隙度分布。 9.一种电子设备， 包括存储器、 处理器及存储在存储器上并可在处理器上运行的计算 机程序， 其特征在于， 所述处理器执行所述程序时实现权利要求1。

8、至4任一项所述的储层孔 隙度预测方法的步骤。 10.一种计算机可读存储介质， 其上存储有计算机程序， 其特征在于， 该计算机程序被 处理器执行时实现权利要求1至4任一项所述的储层孔隙度预测方法的步骤。 权利要求书 2/2 页 3 CN 111751879 A 3 储层孔隙度预测方法及装置 技术领域 0001 本申请涉及地质勘探领域， 具体涉及一种储层孔隙度预测方法及装置。 背景技术 0002 在石油地球物理勘探过程中， 储层岩石的物理特性决定了储层圈闭的储量规模， 直接影响到石油勘探的经济效益。 储层岩石的物理特性通常包括孔隙度、 含烃饱和度以及 渗透率等信息， 其中孔隙度表征了地下岩石的有效。

9、储集空间的大小， 含烃饱和度表征了单 位储集空间内石油和天然气的含量， 渗透率表征了储层岩石孔隙中流体的流动性。 通常情 况下， 对于常规储层岩石， 孔隙度信息决定了储层的可采规模以及开发潜力， 是影响储层资 源量计算的重要参数。 在地震勘探实践中， 通过构建孔隙度与弹性参数之间的关系式， 可以 间接计算储层孔隙度分布情况， 现阶段常规的方法是利用叠后地震反演的纵波阻抗数据体 计算目的层孔隙度体。 0003 随着勘探程度的深入， 基于弹性信息预测储层孔隙度已经获得许多有益应用， 然 而由于横波信息的缺失以及测量精度等原因使得其在储层孔隙度预测应用方面受到限制。 近年来， 随着岩石物理模型的不断。

10、完善， 使得基于双相介质的横波预测精度越来越高， 越来 越能够有效反映储层岩石的孔隙特征以及流体特征， 导致基于弹性信息的孔隙度真实值定 量预测成为重要发展趋势以及研究热点。 0004 发明人发现， 常规的基于纵波阻抗的孔隙度定量预测由于受到储层流体的影响， 导致孔隙度与纵波阻抗之间的关系具有多解性， 使得孔隙度预测结果准确性有所降低， 不 利于储层精细描述。 发明内容 0005 针对现有技术中的问题， 本申请提供一种储层孔隙度预测方法及装置， 能够通过 增加横波约束降低孔隙流体对于孔隙度预测的不利影响， 提高储层孔隙度定量预测的精 度。 0006 为了解决上述问题中的至少一个， 本申请提供以。

11、下技术方案： 0007 第一方面， 本申请提供一种储层孔隙度预测方法， 包括： 0008 根据目标工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关系和横波阻抗与 孔隙度的对应关系； 0009 根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对应关系， 确 定泊松阻抗属性与孔隙度的对应关系； 0010 根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工区的三维地震数据体， 确定所述目标工区的孔隙度分布。 0011 进一步地， 所述根据目标工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关 系和横波阻抗与孔隙度的对应关系， 包括： 0012 根据测井曲线数据中的纵波速度、 横波速度以及。

12、地层密度， 确定纵波阻抗和横波 说明书 1/11 页 4 CN 111751879 A 4 阻抗； 0013 根据预设岩石物理分析模型和预设拟合系数分别对所述纵波阻抗与孔隙度和所 述横波阻抗与孔隙度进行拟合分析， 确定所述纵波阻抗与所述孔隙度的对应关系和所述横 波阻抗与所述孔隙度的对应关系。 0014 进一步地， 所述根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度 的对应关系， 确定泊松阻抗属性与孔隙度的对应关系， 包括： 0015 根据纵波阻抗与孔隙度的对应关系、 横波阻抗与孔隙度的对应关系以及与所述目 标工区对应的预设泊松阻抗调节因子， 确定所述目标工区中泊松阻抗属性与孔隙度的对。

13、应 关系。 0016 进一步地， 所述根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工区的三 维地震数据体， 确定所述目标工区的孔隙度分布， 包括： 0017 根据泊松阻抗属性与孔隙度的对应关系和所述目标工区的泊松阻抗属性地震数 据体， 确定所述目标工区的孔隙度分布。 0018 第二方面， 本申请提供一种储层孔隙度预测装置， 包括： 0019 阻抗确定模块， 用于根据目标工区的测井曲线数据， 确定纵波阻抗与孔隙度的对 应关系和横波阻抗与孔隙度的对应关系； 0020 预测模型构建模块， 用于根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗 与孔隙度的对应关系， 确定泊松阻抗属性与孔隙度的对应关。

14、系； 0021 孔隙度预测模块， 用于根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工 区的三维地震数据体， 确定所述目标工区的孔隙度分布。 0022 进一步地， 所述阻抗确定模块包括： 0023 横纵阻抗确定单元， 用于根据测井曲线数据中的纵波速度、 横波速度以及地层密 度， 确定纵波阻抗和横波阻抗； 0024 横纵阻抗孔隙度关系确定单元， 用于根据预设岩石物理分析模型和预设拟合系数 分别对所述纵波阻抗与孔隙度和所述横波阻抗与孔隙度进行拟合分析， 确定所述纵波阻抗 与所述孔隙度的对应关系和所述横波阻抗与所述孔隙度的对应关系。 0025 进一步地， 所述预测模型构建模块包括： 0026 泊松。

15、阻抗孔隙度关系确定单元， 用于根据纵波阻抗与孔隙度的对应关系、 横波阻 抗与孔隙度的对应关系以及与所述目标工区对应的预设泊松阻抗调节因子， 确定所述目标 工区中泊松阻抗属性与孔隙度的对应关系。 0027 进一步地， 所述孔隙度预测模块包括： 0028 孔隙度分布确定单元， 用于根据泊松阻抗属性与孔隙度的对应关系和所述目标工 区的泊松阻抗属性地震数据体， 确定所述目标工区的孔隙度分布。 0029 第三方面， 本申请提供一种电子设备， 包括存储器、 处理器及存储在存储器上并可 在处理器上运行的计算机程序， 所述处理器执行所述程序时实现所述的储层孔隙度预测方 法的步骤。 0030 第四方面， 本申请。

16、提供一种计算机可读存储介质， 其上存储有计算机程序， 该计算 机程序被处理器执行时实现所述的储层孔隙度预测方法的步骤。 0031 由上述技术方案可知， 本申请提供一种储层孔隙度预测方法及装置， 通过根据目 说明书 2/11 页 5 CN 111751879 A 5 标工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关系和横波阻抗与孔隙度的对应关 系， 再根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对应关系， 确定 泊松阻抗属性与孔隙度的对应关系， 进而根据所述泊松阻抗属性与孔隙度的对应关系和所 述目标工区的三维地震数据体， 确定所述目标工区的孔隙度分布， 本申请通过拟合分析分 别。

17、得到纵波阻抗与孔隙度以及横波阻抗与孔隙度之间的线性关系， 然后根据泊松阻抗属性 的计算公式， 综合利用纵波阻抗与横波阻抗信息构建泊松阻抗属性与孔隙度之间的定量关 系， 通过增加横波约束降低孔隙流体对于孔隙度预测的不利影响， 提高储层孔隙度定量预 测的精度。 附图说明 0032 为了更清楚地说明本申请实施例或现有技术中的技术方案， 下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍， 显而易见地， 下面描述中的附图是本申请 的一些实施例， 对于本领域普通技术人员来讲， 在不付出创造性劳动的前提下， 还可以根据 这些附图获得其他的附图。 0033 图1为本申请实施例中的储层孔隙度预测方法。

18、的流程示意图之一； 0034 图2为本申请实施例中的储层孔隙度预测方法的流程示意图之二； 0035 图3为本申请实施例中的储层孔隙度预测装置的结构图之一； 0036 图4为本申请实施例中的储层孔隙度预测装置的结构图之二； 0037 图5为本申请实施例中的储层孔隙度预测装置的结构图之三； 0038 图6为本申请实施例中的储层孔隙度预测装置的结构图之三； 0039 图7为本申请一实施例中纵波速度、 横波速度、 地层密度、 孔隙度曲线以及计算得 到的泊松阻抗曲线的示意图； 0040 图8为本申请一实施例中纵波阻抗与孔隙度交汇分析以及拟合结果的示意图； 0041 图9为本申请一实施例中横波阻抗与孔隙度。

19、交汇分析以及拟合结果的示意图； 0042 图10为本申请一实施例中构建的泊松阻抗与孔隙度线性关系式以及交汇分析结 果的示意图； 0043 图11为本申请一实施例中基于泊松阻抗属性预测得到的孔隙度分布结果的示意 图； 0044 图12为本申请实施例中的电子设备的结构示意图。 具体实施方式 0045 为使本申请实施例的目的、 技术方案和优点更加清楚， 下面将结合本申请实施例 中的附图， 对本申请实施例中的技术方案进行清楚、 完整的描述， 显然， 所描述的实施例是 本申请一部分实施例， 而不是全部的实施例。 基于本申请中的实施例， 本领域普通技术人员 在没有作出创造性劳动前提下所获得的所有其他实施例。

20、， 都属于本申请保护的范围。 0046 考虑到常规的基于纵波阻抗的孔隙度定量预测由于受到储层流体的影响， 导致孔 隙度与纵波阻抗之间的关系具有多解性， 使得孔隙度预测结果准确性有所降低， 不利于储 层精细描述的问题， 本申请提供一种储层孔隙度预测方法及装置， 通过根据目标工区的测 井曲线数据， 确定纵波阻抗与孔隙度的对应关系和横波阻抗与孔隙度的对应关系， 再根据 说明书 3/11 页 6 CN 111751879 A 6 所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对应关系， 确定泊松阻抗属 性与孔隙度的对应关系， 进而根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工区 的三维地震。

21、数据体， 确定所述目标工区的孔隙度分布， 本申请通过拟合分析分别得到纵波 阻抗与孔隙度以及横波阻抗与孔隙度之间的线性关系， 然后根据泊松阻抗属性的计算公 式， 综合利用纵波阻抗与横波阻抗信息构建泊松阻抗属性与孔隙度之间的定量关系， 通过 增加横波约束降低孔隙流体。 0047 为了能够通过增加横波约束降低孔隙流体对于孔隙度预测的不利影响， 提高储层 孔隙度定量预测的精度， 本申请提供一种储层孔隙度预测方法的实施例， 参见图1， 所述储 层孔隙度预测方法具体包含有如下内容： 0048 步骤S101： 根据目标工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关系和 横波阻抗与孔隙度的对应关系。 00。

22、49 可选的， 所述测井曲线数据包括但不限于： 纵波速度、 横波速度、 地层密度以及孔 隙度曲线。 0050 可选的， 本申请接收到目标工区的测井曲线数据后， 基于输入的测井曲线数据分 别计算纵波阻抗曲线与横波阻抗属性。 0051 可选的， 交汇分析纵波阻抗与孔隙度以及横波阻抗与孔隙度之间的关系， 拟合得 到纵波阻抗与孔隙度以及横波阻抗与孔隙度对应的线性关系式。 0052 步骤S102： 根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对 应关系， 确定泊松阻抗属性与孔隙度的对应关系。 0053 可选的， 基于上述拟合分析结果可以构建泊松阻抗属性与孔隙度之间的关系式， 以此确定泊松阻。

23、抗属性与孔隙度的对应关系。 0054 步骤S103： 根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工区的三维地 震数据体， 确定所述目标工区的孔隙度分布。 0055 可选的， 可以通过输入纵波阻抗以及横波阻抗地震数据体并基于计算得到的泊松 阻抗属性体定量计算储层孔隙度分布。 0056 从上述描述可知， 本申请实施例提供的储层孔隙度预测方法， 能够通过根据目标 工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关系和横波阻抗与孔隙度的对应关 系， 再根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对应关系， 确定 泊松阻抗属性与孔隙度的对应关系， 进而根据所述泊松阻抗属性与孔隙度的对。

24、应关系和所 述目标工区的三维地震数据体， 确定所述目标工区的孔隙度分布， 本申请通过拟合分析分 别得到纵波阻抗与孔隙度以及横波阻抗与孔隙度之间的线性关系， 然后根据泊松阻抗属性 的计算公式， 综合利用纵波阻抗与横波阻抗信息构建泊松阻抗属性与孔隙度之间的定量关 系， 通过增加横波约束降低孔隙流体。 0057 为了能够准确确定横波和纵波的阻抗以及其与孔隙度的对应关系， 在本申请的储 层孔隙度预测方法的一实施例中， 参见图2， 还可以具体包含如下内容： 0058 步骤S201： 根据测井曲线数据中的纵波速度、 横波速度以及地层密度， 确定纵波阻 抗和横波阻抗。 0059 步骤S202： 根据预设岩石。

25、物理分析模型和预设拟合系数分别对所述纵波阻抗与孔 隙度和所述横波阻抗与孔隙度进行拟合分析， 确定所述纵波阻抗与所述孔隙度的对应关系 说明书 4/11 页 7 CN 111751879 A 7 和所述横波阻抗与所述孔隙度的对应关系。 0060 可选的， 基于输入的测井曲线数据分别计算纵波阻抗曲线与横波阻抗属性。 0061 具体的， 纵波阻抗计算公式为： 0062 AIVpDen； 0063 横波阻抗计算公式为： 0064 SIVsDen； 0065 其中， AI为纵波阻抗， Vp为纵波速度， Den为密度， SI为横波阻抗， Vs为横波速度。 0066 可选的， 分别对已知井纵波阻抗与孔隙度以及。

26、横波阻抗与孔隙度进行拟合分析， 建立对应的线性关系式， 其中纵波阻抗与孔隙度线性关系式为： AIk1POR+d1； 横波阻抗 与孔隙度线性关系式为： SIk2POR+d2， 其中， AI为纵波阻抗， POR为孔隙度， SI为横波阻 抗， k1、 d1、 k2、 d2分别为拟合系数。 0067 为了能够根据横波阻抗和纵波阻抗与孔隙度的对应关系确定泊松阻抗属性与孔 隙度的对应关系， 在本申请的储层孔隙度预测方法的一实施例中， 还可以具体包含如下内 容： 0068 根据纵波阻抗与孔隙度的对应关系、 横波阻抗与孔隙度的对应关系以及与所述目 标工区对应的预设泊松阻抗调节因子， 确定所述目标工区中泊松阻抗。

27、属性与孔隙度的对应 关系。 0069 可选的， 根据泊松阻抗计算公式PIAI-cSI， 将拟合得到的纵波阻抗、 横波阻抗 与孔隙度的线性关系式代入进去， 得到泊松阻抗属性与孔隙度之间的关系式， 即： 0070 PI(k1POR+d1)-c(k2POR+d2)(k1-ck2)POR+(d1-cd2)； 0071 其中， c为泊松阻抗调节因子， 可以根据已知井岩石物理分析获取。 0072 为了能够根据泊松阻抗属性与孔隙度的对应关系对目标工区的孔隙度分布进行 精准预测， 在本申请的储层孔隙度预测方法的一实施例中， 还可以具体包含如下内容： 0073 根据泊松阻抗属性与孔隙度的对应关系和所述目标工区的。

28、泊松阻抗属性地震数 据体， 确定所述目标工区的孔隙度分布。 0074 可选的， 输入的纵波阻抗以及横波阻抗地震数据体基于常规的叠前地震反演得 到， 泊松阻抗属性计算公式为： PIAI-cSI。 基于泊松阻抗属性与孔隙度之间的关系得到 孔隙度定量计算公式： 0075 0076 利用上述公式， 结合泊松阻抗属性以及拟合系数可以定量计算储层孔隙度分布， 相对于传统的单纯利用纵波阻抗计算孔隙度的方法， 该方法通过引入横波阻抗约束降低了 孔隙流体对于孔隙度预测结果的不利影响， 提高了预测精度。 0077 为了能够通过增加横波约束降低孔隙流体对于孔隙度预测的不利影响， 提高储层 孔隙度定量预测的精度， 本。

29、申请提供一种用于实现所述储层孔隙度预测方法的全部或部分 内容的储层孔隙度预测装置的实施例， 参见图3， 所述储层孔隙度预测装置具体包含有如下 内容： 0078 阻抗确定模块10， 用于根据目标工区的测井曲线数据， 确定纵波阻抗与孔隙度的 对应关系和横波阻抗与孔隙度的对应关系。 说明书 5/11 页 8 CN 111751879 A 8 0079 预测模型构建模块20， 用于根据所述纵波阻抗与孔隙度的对应关系和所述横波阻 抗与孔隙度的对应关系， 确定泊松阻抗属性与孔隙度的对应关系。 0080 孔隙度预测模块30， 用于根据所述泊松阻抗属性与孔隙度的对应关系和所述目标 工区的三维地震数据体， 确定。

30、所述目标工区的孔隙度分布。 0081 从上述描述可知， 本申请实施例提供的储层孔隙度预测装置， 能够通过根据目标 工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关系和横波阻抗与孔隙度的对应关 系， 再根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对应关系， 确定 泊松阻抗属性与孔隙度的对应关系， 进而根据所述泊松阻抗属性与孔隙度的对应关系和所 述目标工区的三维地震数据体， 确定所述目标工区的孔隙度分布， 本申请通过拟合分析分 别得到纵波阻抗与孔隙度以及横波阻抗与孔隙度之间的线性关系， 然后根据泊松阻抗属性 的计算公式， 综合利用纵波阻抗与横波阻抗信息构建泊松阻抗属性与孔隙度之间的。

31、定量关 系， 通过增加横波约束降低孔隙流体。 0082 为了能够准确确定横波和纵波的阻抗以及其与孔隙度的对应关系， 在本申请的储 层孔隙度预测装置的一实施例中， 参见图4， 所述阻抗确定模块10包括： 0083 横纵阻抗确定单元11， 用于根据测井曲线数据中的纵波速度、 横波速度以及地层 密度， 确定纵波阻抗和横波阻抗。 0084 横纵阻抗孔隙度关系确定单元12， 用于根据预设岩石物理分析模型和预设拟合系 数分别对所述纵波阻抗与孔隙度和所述横波阻抗与孔隙度进行拟合分析， 确定所述纵波阻 抗与所述孔隙度的对应关系和所述横波阻抗与所述孔隙度的对应关系。 0085 为了能够根据横波阻抗和纵波阻抗与孔。

32、隙度的对应关系确定泊松阻抗属性与孔 隙度的对应关系， 在本申请的储层孔隙度预测装置的一实施例中， 参见图5， 所述预测模型 构建模块20包括： 0086 泊松阻抗孔隙度关系确定单元21， 用于根据纵波阻抗与孔隙度的对应关系、 横波 阻抗与孔隙度的对应关系以及与所述目标工区对应的预设泊松阻抗调节因子， 确定所述目 标工区中泊松阻抗属性与孔隙度的对应关系。 0087 为了能够根据泊松阻抗属性与孔隙度的对应关系对目标工区的孔隙度分布进行 精准预测， 在本申请的储层孔隙度预测装置的一实施例中， 参见图6， 所述孔隙度预测模块 30包括： 0088 孔隙度分布确定单元31， 用于根据泊松阻抗属性与孔隙度。

33、的对应关系和所述目标 工区的泊松阻抗属性地震数据体， 确定所述目标工区的孔隙度分布。 0089 为了更进一步说明本方案， 本申请还提供一种应用上述储层孔隙度预测装置实现 储层孔隙度预测方法的具体应用实例， 具体包含有如下内容： 0090 步骤1： 输入目的层测井曲线数据。 0091 具体的， 获取研究工区内已知井的测井曲线数据， 包括： 纵波速度、 横波速度、 密度 曲线以及孔隙度曲线， 测井曲线如图7所示， 已知井目的层段深度范围为4615-4632米。 0092 步骤2： 基于输入的测井曲线数据分别计算纵波阻抗曲线与横波阻抗属性。 0093 具体的， 纵波阻抗计算公式为： 0094 AIV。

34、pDen； 0095 横波阻抗计算公式为： 说明书 6/11 页 9 CN 111751879 A 9 0096 SIVsDen； 0097 其中， AI为纵波阻抗， Vp为纵波速度， Den为密度， SI为横波阻抗， Vs为横波速度。 0098 步骤3： 交汇分析纵波阻抗与孔隙度以及横波阻抗与孔隙度之间的关系， 拟合得到 纵波阻抗与孔隙度以及横波阻抗与孔隙度对应的线性关系式。 0099 具体的， 根据岩石物理分析方法， 交汇分析纵波阻抗与孔隙度以及横波阻抗与孔 隙度之间的关系， 并拟合得到纵波阻抗与孔隙度以及横波阻抗与孔隙度对应的线性关系。 0100 如图8所示， 为纵波阻抗与孔隙度交汇分析。

35、结果， 图中横坐标为孔隙度POR， 纵坐标 为纵波阻抗AI， 拟合得到的线性关系式为AI-20250POR-19720， 相关系数R20.902， 可 以看出纵波阻抗与孔隙度具有很好的线性关系， 随着孔隙度增大， 纵波阻抗属性线性减小。 0101 如图9所示， 为横波阻抗与孔隙度交汇分析结果， 图中横坐标为孔隙度POR， 纵坐标 为横波阻抗SI， 拟合得到的线性关系为SI-8276POR+10440， 相关系数R20.807， 可以 看出相对于纵波阻抗， 横波阻抗与孔隙度之间的线性关系相对较弱， 但趋势稳定。 0102 步骤4： 基于拟合分析结果构建泊松阻抗属性与孔隙度之间的关系式。 0103。

36、 根据泊松阻抗计算公式PIAI-cSI， 将拟合得到的纵波阻抗、 横波阻抗与孔隙 度的线性关系式代入进去， 得到泊松阻抗属性与孔隙度之间的关系式： 0104 PI(-20250POR-19720)-c(-8276POR+10440)； 0105 根据岩石物理分析结果以及工区背景资料分析， 确定目的层段泊松阻抗调节因子 c0.7。 将c值代入到泊松阻抗属性PI公式， 则建立泊松阻抗属性与孔隙度之间的关系式 为： 0106 PI-14450POR+12410。 0107 如图10展示了泊松阻抗属性与孔隙度交汇分析结果， 图中横坐标为孔隙度POR， 纵 坐标为泊松阻抗PI， 泊松阻抗与孔隙度线性关系。

37、式对应的拟合线如图中实线所示， 相关系 数R20.946， 可以看出相对于纵波阻抗以及横波阻抗， 泊松阻抗属性与孔隙度具有更好的 线性相关性， 尤其当孔隙度大于0.1时， 基于泊松阻抗计算的孔隙度值与实测孔隙度误差很 小。 0108 步骤5： 输入纵波阻抗以及横波阻抗地震数据体并基于计算得到的泊松阻抗属性 体定量计算储层孔隙度分布。 0109 输入的纵波阻抗以及横波阻抗地震数据体基于常规的叠前地震反演得到， 根据泊 松阻抗计算公式得到泊松阻抗属性地震数据体： PIAI-0.7SI。 然后步骤4得到的泊松阻 抗属性与孔隙度之间的关系得到孔隙度定量计算公式： 0110 POR-(PI-12410)。

38、/14450。 0111 图11展示了研究区基于泊松阻抗属性预测的孔隙度分布图， 图11中井W-1在目的 层段为高孔气层， 根据钻井以及岩心测试结果， 在目的层段W-1井孔隙度为0.18左右。 本实 施例中， 根据泊松阻抗属性预测的孔隙度在W-1井点处的POR0.19左右， 与岩心测试结果 相对误差为5， 表明该方法具有很高的预测精度。 0112 由上述描述可知， 本申请至少还可以实现如下的技术效果： 0113 通过拟合分析分别得到纵波阻抗与孔隙度以及横波阻抗与孔隙度之间的线性关 系， 然后根据泊松阻抗属性的计算公式， 综合利用纵波阻抗与横波阻抗信息构建泊松阻抗 属性与孔隙度之间的定量关系， 。

39、通过增加横波约束降低孔隙流体对于孔隙度预测的不利影 说明书 7/11 页 10 CN 111751879 A 10 响， 提高储层孔隙度定量预测的精度。 0114 从硬件层面来说， 为了能够通过增加横波约束降低孔隙流体对于孔隙度预测的不 利影响， 提高储层孔隙度定量预测的精度， 本申请提供一种用于实现所述储层孔隙度预测 方法中的全部或部分内容的电子设备的实施例， 所述电子设备具体包含有如下内容： 0115 处理器(processor)、 存储器(memory)、 通信接口(Communications Interface)和 总线； 其中， 所述处理器、 存储器、 通信接口通过所述总线完成相互。

40、间的通信； 所述通信接口 用于实现储层孔隙度预测装置与核心业务系统、 用户终端以及相关数据库等相关设备之间 的信息传输； 该逻辑控制器可以是台式计算机、 平板电脑及移动终端等， 本实施例不限于 此。 在本实施例中， 该逻辑控制器可以参照实施例中的储层孔隙度预测方法的实施例， 以及 储层孔隙度预测装置的实施例进行实施， 其内容被合并于此， 重复之处不再赘述。 0116 可以理解的是， 所述用户终端可以包括智能手机、 平板电子设备、 网络机顶盒、 便 携式计算机、 台式电脑、 个人数字助理(PDA)、 车载设备、 智能穿戴设备等。 其中， 所述智能穿 戴设备可以包括智能眼镜、 智能手表、 智能手环。

41、等。 0117 在实际应用中， 储层孔隙度预测方法的部分可以在如上述内容所述的电子设备侧 执行， 也可以所有的操作都在所述客户端设备中完成。 具体可以根据所述客户端设备的处 理能力， 以及用户使用场景的限制等进行选择。 本申请对此不作限定。 若所有的操作都在所 述客户端设备中完成， 所述客户端设备还可以包括处理器。 0118 上述的客户端设备可以具有通信模块(即通信单元)， 可以与远程的服务器进行通 信连接， 实现与所述服务器的数据传输。 所述服务器可以包括任务调度中心一侧的服务器， 其他的实施场景中也可以包括中间平台的服务器， 例如与任务调度中心服务器有通信链接 的第三方服务器平台的服务器。。

42、 所述的服务器可以包括单台计算机设备， 也可以包括多个 服务器组成的服务器集群， 或者分布式装置的服务器结构。 0119 图12为本申请实施例的电子设备9600的系统构成的示意框图。 如图12所示， 该电 子设备9600可以包括中央处理器9100和存储器9140； 存储器9140耦合到中央处理器9100。 值得注意的是， 该图12是示例性的； 还可以使用其他类型的结构， 来补充或代替该结构， 以 实现电信功能或其他功能。 0120 一实施例中， 储层孔隙度预测方法功能可以被集成到中央处理器9100中。 其中， 中 央处理器9100可以被配置为进行如下控制： 0121 步骤S101： 根据目标工。

43、区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关系和 横波阻抗与孔隙度的对应关系。 0122 步骤S102： 根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对 应关系， 确定泊松阻抗属性与孔隙度的对应关系。 0123 步骤S103： 根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工区的三维地 震数据体， 确定所述目标工区的孔隙度分布。 0124 从上述描述可知， 本申请实施例提供的电子设备， 通过根据目标工区的测井曲线 数据， 确定纵波阻抗与孔隙度的对应关系和横波阻抗与孔隙度的对应关系， 再根据所述纵 波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对应关系， 确定泊松阻抗属性与。

44、孔 隙度的对应关系， 进而根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工区的三维 地震数据体， 确定所述目标工区的孔隙度分布， 本申请通过拟合分析分别得到纵波阻抗与 说明书 8/11 页 11 CN 111751879 A 11 孔隙度以及横波阻抗与孔隙度之间的线性关系， 然后根据泊松阻抗属性的计算公式， 综合 利用纵波阻抗与横波阻抗信息构建泊松阻抗属性与孔隙度之间的定量关系， 通过增加横波 约束降低孔隙流体。 0125 在另一个实施方式中， 储层孔隙度预测装置可以与中央处理器9100分开配置， 例 如可以将储层孔隙度预测装置配置为与中央处理器9100连接的芯片， 通过中央处理器的控 制来。

45、实现储层孔隙度预测方法功能。 0126 如图12所示， 该电子设备9600还可以包括： 通信模块9110、 输入单元9120、 音频处 理器9130、 显示器9160、 电源9170。 值得注意的是， 电子设备9600也并不是必须要包括图12 中所示的所有部件； 此外， 电子设备9600还可以包括图12中没有示出的部件， 可以参考现有 技术。 0127 如图12所示， 中央处理器9100有时也称为控制器或操作控件， 可以包括微处理器 或其他处理器装置和/或逻辑装置， 该中央处理器9100接收输入并控制电子设备9600的各 个部件的操作。 0128 其中， 存储器9140， 例如可以是缓存器、 。

46、闪存、 硬驱、 可移动介质、 易失性存储器、 非 易失性存储器或其它合适装置中的一种或更多种。 可储存上述与失败有关的信息， 此外还 可存储执行有关信息的程序。 并且中央处理器9100可执行该存储器9140存储的该程序， 以 实现信息存储或处理等。 0129 输入单元9120向中央处理器9100提供输入。 该输入单元9120例如为按键或触摸输 入装置。 电源9170用于向电子设备9600提供电力。 显示器9160用于进行图像和文字等显示 对象的显示。 该显示器例如可为LCD显示器， 但并不限于此。 0130 该存储器9140可以是固态存储器， 例如， 只读存储器(ROM)、 随机存取存储器 (。

47、RAM)、 SIM卡等。 还可以是这样的存储器， 其即使在断电时也保存信息， 可被选择性地擦除 且设有更多数据， 该存储器的示例有时被称为EPROM等。 存储器9140还可以是某种其它类型 的装置。 存储器9140包括缓冲存储器9141(有时被称为缓冲器)。 存储器9140可以包括应用/ 功能存储部9142， 该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央 处理器9100执行电子设备9600的操作的流程。 0131 存储器9140还可以包括数据存储部9143， 该数据存储部9143用于存储数据， 例如 联系人、 数字数据、 图片、 声音和/或任何其他由电子设备使用的数据。 。

48、存储器9140的驱动程 序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如 消息传送应用、 通讯录应用等)的各种驱动程序。 0132 通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。 通信模 块(发送机/接收机)9110耦合到中央处理器9100， 以提供输入信号和接收输出信号， 这可以 和常规移动通信终端的情况相同。 0133 基于不同的通信技术， 在同一电子设备中， 可以设置有多个通信模块9110， 如蜂窝 网络模块、 蓝牙模块和/或无线局域网模块等。 通信模块(发送机/接收机)9110还经由音频 处理器9130耦合到扬声器9131。

49、和麦克风9132， 以经由扬声器9131提供音频输出， 并接收来 自麦克风9132的音频输入， 从而实现通常的电信功能。 音频处理器9130可以包括任何合适 的缓冲器、 解码器、 放大器等。 另外， 音频处理器9130还耦合到中央处理器9100， 从而使得可 以通过麦克风9132能够在本机上录音， 且使得可以通过扬声器9131来播放本机上存储的声 说明书 9/11 页 12 CN 111751879 A 12 音。 0134 本申请的实施例还提供能够实现上述实施例中的执行主体为服务器或客户端的 储层孔隙度预测方法中全部步骤的一种计算机可读存储介质， 所述计算机可读存储介质上 存储有计算机程序，。

50、 该计算机程序被处理器执行时实现上述实施例中的执行主体为服务器 或客户端的储层孔隙度预测方法的全部步骤， 例如， 所述处理器执行所述计算机程序时实 现下述步骤： 0135 步骤S101： 根据目标工区的测井曲线数据， 确定纵波阻抗与孔隙度的对应关系和 横波阻抗与孔隙度的对应关系。 0136 步骤S102： 根据所述纵波阻抗与孔隙度的对应关系和所述横波阻抗与孔隙度的对 应关系， 确定泊松阻抗属性与孔隙度的对应关系。 0137 步骤S103： 根据所述泊松阻抗属性与孔隙度的对应关系和所述目标工区的三维地 震数据体， 确定所述目标工区的孔隙度分布。 0138 从上述描述可知， 本申请实施例提供的计算。