建立三维渗透率模型的方法和装置制造方法及图纸

本申请实施例提供了一种建立三维渗透率模型的方法和装置，其中，该方法包括以下步骤：获取目标油藏的结果数据；根据目标油藏的结果数据，建立目标油藏的应力场模拟地质模型；根据目标油藏的结果数据和目标油藏成像测井资料，通过应力场模拟地质模型，计算得到目标油藏应力场分布数据；根据目标油藏应力场分布数据、目标油藏成像测井资料和目标油藏试井资料，求解得到校正渗透率；根据校正渗透率建立三维渗透率模型。由于该方案考虑了微裂缝的影响，并根据校正渗透率建立三维渗透率模型，因此解决了现有建模方法中存在的三维渗透率模型精度低、误差大的技术问题，实现了提高所建立的三维渗透率模型精度的技术效果。

【技术实现步骤摘要】

本申请涉及油藏开发
，特别涉及一种建立三维渗透率模型的方法和装置。
技术介绍
在油藏开发领域中，对目标油藏进行研究分析时，往往会先根据目标油藏的测井资料和实验结果数据，建立关于目标油藏的三维渗透率模型，再通过该三维渗透率模型对目标油藏展开具体的油藏工程研究和油藏数值模拟研究。目前，主要采用的建模方法普遍是根据测井渗透率，建立关于目标油藏的三维渗透率模型。然而，在具体实施时，目标油藏的地质储层中往往会发育有许多微裂缝，这些微裂缝会影响多孔介质中的流体的运移，进而导致在实际开发中，地质储层的实际渗透率会比测井渗透率高很多。因此，导致通过现有的建模方法，直接根据测井渗透率建立的三维渗透率模型往往存在精度低、误差大的技术问题。在对特低渗储层进行研究时，现有的建模方法所产生的上述技术问题尤为明显。针对上述问题，目前尚未提出有效的解决方案。
技术实现思路
本申请实施例提供了一种建立三维渗透率模型的方法和装置，以解决现有建模方法存在的模型精度低、误差大的技术问题。本申请实施例提供了一种建立三维渗透率模型的方法，包括：获取目标油藏的结果数据；根据所述目标油藏的结果数据，建立目标油藏的应力场模拟地质模型；根据所述目标油藏的结果数据和目标油藏的测井资料，通过所述应力场模拟地质模型，计算得到目标油藏应力场分布数据；根据所述目标油藏应力场分布数据、所述目标油藏的测井资料和目标油藏的试井资料，求解得到校正渗透率；根据所述校正渗透率，建立所述目标油藏的三维渗透率模型。在一个优选实施方式中，所述结果数据包括：岩石力学相关参数、应力场大小和应力场水平最大主应力方向。在一个优选实施方式中，所述获取目标油藏的结果数据，包括：获取目标油藏的钻井岩心样品；对所述目标油藏的钻井岩心样品进行岩石三轴应力实验，得到所述岩石力学相关参数；对所述目标油藏的钻井岩心样品进行岩石声发射实验，得到所述应力场大小；根据所述目标油藏的测井资料，得到所述应力场水平最大主应力方向。在一个优选实施方式中，根据所述目标油藏的结果数据，建立目标油藏的应力场模拟地质模型，包括：根据所述目标油藏的结果数据，通过映射网格划分，建立所述目标油藏的应力场模拟地质模型。在一个优选实施方式中，根据所述目标油藏应力场分布数据、所述目标油藏的测井资料和目标油藏的试井资料，求解得到校正渗透率，包括：根据所述目标油藏应力场分布数据，计算得到目标油藏中各个测点的裂缝破裂率分布参数值；根据所述目标油藏应力场分布数据和所述目标油藏中各个测点的裂缝破裂率分布参数值，计算得到目标油藏中各个测点的裂缝密度预测值；根据所述目标油藏中各个测点的裂缝密度预测值，将所述目标油藏划分为多个裂缝发育密度区；根据所述目标油藏的测井资料和试井资料，求解得到所述目标油藏中各个测点的的校正渗透率。在一个优选实施方式中，所述根据所述目标油藏应力场分布数据，计算得到目标油藏中各个测点的裂缝破裂率分布参数值，包括：根据所述目标油藏应力场分布数据，按照以下公式，计算得到所述目标油藏中各个测点的裂缝破裂率分布参数值：In(j)＝τn(j)/[τn](j)It(j)＝σt(j)/[σt](j)I(j)＝α(j)It(j)+β(j)In(j)其中，[τn](j)为编号为j的测点的岩石的抗剪强度，C(j)为编号为j的测点的正应力为零时岩石的抗剪强度，σn(j)为编号为j的测点的破裂面剪应力大小，为编号为j的测点的内摩擦系数，In(j)为编号为j的测点的裂缝剪破裂率，τn(j)为编号为j的测点的剪应力值，It(j)为编号为j的测点的裂缝张破裂率，σt(j)为编号为j的测点的破裂面张应力大小，[σt](j)为编号为j的测点的岩石抗张强度，α(j)为编号为j的测点的张裂缝比例，β(j)为编号为j的测点的剪裂缝比例，I(j)为编号为j的测点的裂缝破裂率分布参数值。在一个优选实施方式中，所述根据所述目标油藏应力场分布数据和所述目标油藏中各个测点的裂缝破裂率分布参数值，计算得到目标油藏中各个测点的裂缝密度预测值，包括：根据所述目标油藏应力场分布数据和所述目标油藏中各个测点的裂缝破裂率分布参数值，按照以下公式，计算得到所述目标油藏中各个测点的裂缝密度预测值：I(j)＜Ic,ξ(j)＝a1(j)I(j)2+a2(j)I(j)≥Ic,ξ(j)＝a1(j)I(j)2+a2(j)+a3(j)(A(j)+a4(j))其中，I(j)为编号为j的测点的裂缝破裂率分布参数值，Ic为I(j)的临界值，ξ(j)为编号为j的测点的裂缝密度预测值，a1(j)、a2(j)、a3(j)、a4(j)为编号为j的测点的待定常数，A(j)为编号为j的测点的能量值。在一个优选实施方式中，所述根据所述目标油藏的测井资料和试井资料，求解得到所述目标油藏中各个测点的校正渗透率，包括：根据所述目标油藏的测井资料，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中的测井渗透率的平均值；根据所述目标油藏试井资料，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中的试井渗透率；根据所述多个裂缝发育密度区中各个裂缝发育密度区中的测井渗透率的平均值和试井渗透率值，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中校正系数；根据所述目标油藏的测井资料和所述多个裂缝发育密度区中各个裂缝发育密度区中校正系数，求解得到所述目标油藏中各个测点的校正渗透率。在一个优选实施方式中，根据所述目标油藏的测井资料，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中的测井渗透率的平均值，包括：根据所述目标油藏的测井资料，按照以下公式，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中的测井渗透率的平均值：其中，为目标油藏中第F个裂缝密度发育密度区中的测井渗透率平均值，K(φ(j),MZ(j),Vsh(j))(j)为第F个裂缝密度发育密度区中第i小层中的编号为j的测点的测井渗透率，φ(j)、MZ(j)、Vsh(j)、Vsh分别为第F个裂缝密度发育密度区中第i小层中的编号为j的测点的孔隙度、粒度中值、泥质含量，hi(F)为第F个裂缝密度发育密度区中第i小层厚度，h(F)为第F个裂缝密度发育密度区中油藏实际生产层段厚度，ni为第F个裂缝密度发育密度区中第i小层内所包含的测点数，N(F)为第F个裂缝密度发育密度区中所包含小层数。在一个优选实施方式中，根据所述目标油藏的试井资料，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中的试井渗透率，包括：根据所述目标油藏的试井资料，按照以下公式，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中的试井渗透率：其中，KT(F)为目标油藏中第F个裂缝密度发育密度区中的试井渗透率，Q(F)为目标油藏中第F个裂缝密度发育密度区中的试井产量，μ(F)为第F个裂缝密度发育密度区中的地层原油粘度，re(F)为第F个裂缝密度发育密度区中的供给半径，rw(F)为第F个裂缝密度发育密度区中的油井半径，ΔP(F)为第F个裂缝密度发育密度区中的生产压差，h(F)为第F个裂缝密度发育密度区中的试井油藏厚度。在一个优选实施方式中，根据所述多个裂缝发育密度区中各个裂缝发育密度区中的测井渗透率的平均值和试井渗透率值，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中校正系数，包括：根据所本文档来自技高网...

【技术保护点】
一种建立三维渗透率模型的方法，其特征在于，包括：获取目标油藏的结果数据；根据所述目标油藏的结果数据，建立目标油藏的应力场模拟地质模型；根据所述目标油藏的结果数据和目标油藏的测井资料，通过所述应力场模拟地质模型，计算得到目标油藏应力场分布数据；根据所述目标油藏应力场分布数据、所述目标油藏的测井资料和目标油藏的试井资料，求解得到校正渗透率；根据所述校正渗透率，建立目标油藏的三维渗透率模型。

【技术特征摘要】
1.一种建立三维渗透率模型的方法，其特征在于，包括：获取目标油藏的结果数据；根据所述目标油藏的结果数据，建立目标油藏的应力场模拟地质模型；根据所述目标油藏的结果数据和目标油藏的测井资料，通过所述应力场模拟地质模型，计算得到目标油藏应力场分布数据；根据所述目标油藏应力场分布数据、所述目标油藏的测井资料和目标油藏的试井资料，求解得到校正渗透率；根据所述校正渗透率，建立目标油藏的三维渗透率模型。2.根据权利要求1所述的方法，其特征在于，所述结果数据包括：岩石力学相关参数、应力场大小和应力场水平最大主应力方向。3.根据权利要求2所述的方法，其特征在于，所述获取目标油藏的结果数据，包括：获取目标油藏的钻井岩心样品；对所述目标油藏的钻井岩心样品进行岩石三轴应力实验，得到所述岩石力学相关参数；对所述目标油藏的钻井岩心样品进行岩石声发射实验，得到所述应力场大小；根据所述目标油藏的测井资料，得到所述应力场水平最大主应力方向。4.根据权利要求2所述的方法，其特征在于，根据所述目标油藏的结果数据，建立目标油藏的应力场模拟地质模型，包括：根据所述目标油藏的结果数据，通过映射网格划分，建立所述目标油藏的应力场模拟地质模型。5.根据权利要求2所述的方法，其特征在于，根据所述目标油藏应力场分布数据、所述目标油藏的测井资料和目标油藏的试井资料，求解得到校正渗透率，包括：根据所述目标油藏应力场分布数据，计算得到目标油藏中各个测点的裂缝破裂率分布参数值；根据所述目标油藏应力场分布数据和所述目标油藏中各个测点的裂缝破裂率分布参数值，计算得到目标油藏中各个测点的裂缝密度预测值；根据所述目标油藏中各个测点的裂缝密度预测值，将所述目标油藏划分为多个裂缝发育密度区；根据所述目标油藏的测井资料和所述目标油藏的试井资料，求解得到所述目标油藏中各个测点的校正渗透率。6.根据权利要求5所述的方法，其特征在于，所述根据所述目标油藏应力场分布数据，计算得到目标油藏中各个测点的裂缝破裂率分布参数值，包括：根据所述目标油藏应力场分布数据，按照以下公式，计算得到所述目标油藏中各个测点的裂缝破裂率分布参数值：In(j)＝τn(j)/[τn](j)It(j)＝σt(j)/[σt](j)I(j)＝α(j)It(j)+β(j)In(j)其中，[τn](j)为编号为j的测点的岩石的抗剪强度，C(j)为编号为j的测点的正应力为零时岩石的抗剪强度，σn(j)为编号为j的测点的破裂面剪应力大小，为编号为j的测点的内摩擦系数，In(j)为编号为j的测点的裂缝剪破裂率，τn(j)为编号为j的测点的剪应力值，It(j)为编号为j的测点的裂缝张破裂率，σt(j)为编号为j的测点的破裂面张应力大小，[σt](j)为编号为j的测点的岩石抗张强度，α(j)为编号为j的测点的张裂缝比例，β(j)为编号为j的测点的剪裂缝比例，I(j)为编号为j的测点的裂缝破裂率分布参数值。7.根据权利要求5所述的方法，其特征在于，所述根据所述目标油藏应力场分布数据和所述目标油藏中各个测点的裂缝破裂率分布参数值，计算得到目标油藏中各个测点的裂缝密度预测值，包括：根据所述目标油藏应力场分布数据和所述目标油藏中各个测点的裂缝破裂率分布参数值，按照以下公式，计算得到所述目标油藏中各个测点的裂缝密度预测值：I(j)＜Ic,ξ(j)＝a1(j)I(j)2+a2(j)I(j)≥Ic,ξ(j)＝a1(j)I(j)2+a2(j)+a3(j)(A(j)+a4(j))其中，I(j)为编号为j的测点的裂缝破裂率分布参数值，Ic为I(j)的临界值，ξ(j)为编号为j的测点的裂缝密度预测值，a1(j)、a2(j)、a3(j)、a4(j)为编号为j的测点的待定常数，A(j)为编号为j的测点的能量值。8.根据权利要求5所述的方法，其特征在于，所述根据所述目标油藏的测井资料和所述目标油藏的试井资料，求解得到所述目标油藏中各个测点的校正渗透率，包括：根据所述目标油藏的测井资料，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中的测井渗透率的平均值；根据所述目标油藏试井资料，求解得到所述多个裂缝发育密度区中各个裂缝发育密度区中的试井渗透率；根据所述多个裂缝发育密度区中各个裂缝发育密度区中的测井渗透率的平均值和试井渗透率值，求解得到所述多个裂缝发育密度区中各...

【专利技术属性】
技术研发人员：薛永超，邸鹏伟，毛文辉，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：北京;11