基于测井数据确定毛管力曲线及毛管力曲线分布的方法.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN115034112A(43)申请公布日2022.09.09(21)申请号202210672799.3G06F119/14(2020.01)(22)申请日2022.06.14G06F119/22(2020.01)(71)申请人西南石油大学地址610500四川省成都市新都区新都大道8号(72)发明人彭瑀黄义涛李勇明李政澜伍翊嘉(74)专利代理机构成都其知创新专利代理事务所(普通合伙)51326专利代理师房立普(51)Int.Cl.G06F30/23(2020.01)E21B49/00(2006.01)E21B43/00(2006.01)G06F111/10(2020.01)权利要求书3页说明书7页附图8页(54)发明名称基于测井数据确定毛管力曲线及毛管力曲线分布的方法(57)摘要本发明公开了一种基于测井数据确定毛管力曲线及毛管力曲线分布的方法，确定毛管力曲线的方法包括以下步骤：S1：进行已有岩心进行毛管力测试，获得毛管力曲线；S2：确定曲线拐点，建立拐点处的进汞饱和度和相应毛管力的比值与岩心渗透率的关系式一，门槛压力、拐点处毛管力以及拐点处进汞饱和度的关系式二，拐点处毛管力与岩心渗透率的关系式三，以及毛管力与进汞体积的关系式四；S3：数据拟合得到关系式一和三中的特征参数；S4：推导得到孔吼结构几何因子、门槛压力与岩心渗透率之间的关系式五；S5：根据测井数据获得目标井的岩心渗透率，结合关系式五，确定其对应的毛管力曲线。本发明能够为油藏数值模拟提供更加准确的毛管力曲线数据。CN115034112ACN115034112A权利要求书1/3页1.一种基于测井数据确定毛管力曲线的方法，其特征在于，包括以下步骤：S1：对目标井所处目标区块现有的岩心进行毛管力测试，获得所述岩心对应取心井的毛管力曲线；S2：确定所述毛管力曲线的拐点，建立所述拐点处的进汞饱和度和相应毛管力的比值与岩心渗透率的关系式一，门槛压力、拐点处毛管力以及拐点处进汞饱和度之间的关系式二，拐点处毛管力与岩心渗透率之间的关系式三，以及毛管力与进汞体积的关系式四；S3：根据所述毛管力曲线，结合所述拐点，数据拟合得到所述关系式一中的特征参数一和所述关系式三中的特征参数二；S4：根据所述关系式一、关系式二、关系式三、关系式四以及步骤S3获得的特征参数，推导得到孔吼结构几何因子、门槛压力与岩心渗透率之间的关系式五；S5：根据测井数据获得所述目标井的岩心渗透率，根据目标井的岩心渗透率，结合所述关系式五，即可确定目标井岩心渗透率对应的毛管力曲线。2.根据权利要求1所述的基于测井数据确定毛管力曲线的方法，其特征在于，步骤S2中，所述关系式一为：bK＝a(SHgA/pcA)(1)式中：K为井位网格点处的岩样渗透率；a和b均为常数，即所述特征参数一；SHgA为拐点处的进汞饱和度；pcA为拐点处的毛管压力；所述关系式二为：pd＝pcASHgA(2)式中：pd为门槛压力；所述关系式三为：dpcA＝c/K(3)式中：c和d均为常数，即所述特征参数二；所述关系式四为：lg(Vb/Vb∞)lg(pc/pd)＝‑Fg/2.303(4)式中：Vb为进汞体积；Vb∞为完全连通的孔隙体积；pc为毛细管压力；Fg为孔吼结构几何因子。3.一种基于测井数据确定毛管力曲线分布的方法，其特征在于，首先利用权利要求1或2所述的基于测井数据确定毛管力曲线的方法确定毛管力曲线，然后利用网格插值算法确定井位网格点之间每个网格的毛管力曲线，从而确定毛管力分布。4.根据权利要求3所述的基于测井数据确定毛管力曲线分布的方法，其特征在于，利用网格插值算法确定井位网格点之间每个网格的毛管力曲线，从而确定毛管力分布，具体包括以下步骤：S1'：根据目标区块的井位分布建立包含目标区块的m×m矩阵网格面，假设网格分布区域为Xmin≤X≤Xmax，Ymin≤Y≤Ymax，则网格点H(i，j)的坐标为(xi，yj)；S2'：逐步比较找出距离网格点H(i，j)不同距离的井位网格点(xk，yk)，根据距离的远近分别给予不同的权重系数；S3'：将目标区块边界外的网格点的毛管力曲线取值为0；对目标区块边界内的网格点，以网格点H(i，j)为圆心，以R为搜索半径，搜索得到位于R范围内的所有井位网格点，并将其2CN115034112A权利要求书2/3页按照距离网格点H(i，j)的距离大小进行排列；S4'：根据搜索得到的所有井位网格点的结果，确定各网格点的毛管力曲线：若与网格点H(i，j)距离最近的井位网格点，其与网格点H(i，j)的网格距离d1≤R/10，则网格点H(i，j)的毛管力曲线直接等于该井位网格点毛管压力曲线；若搜索半径R范围内的井位网格数n≥3，则网格点H(i，j)的毛管力采用距离加权平均公式