(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210672799.3 (22)申请日 2022.06.14 (71)申请人 西南石油大 学 地址 610500 四川省成 都市新都区新都大 道8号 (72)发明人 彭瑀　黄义涛　李勇明　李政澜　 伍翊嘉　 (74)专利代理 机构 成都其知创新专利代理事务 所(普通合伙) 51326 专利代理师 房立普 (51)Int.Cl. G06F 30/23(2020.01) E21B 49/00(2006.01) E21B 43/00(2006.01) G06F 111/10(2020.01)G06F 119/14(2020.01) G06F 119/22(2020.01) (54)发明名称 基于测井数据确定毛管力 曲线及毛管力 曲 线分布的方法 (57)摘要 本发明公开了一种基于测井数据确定毛管 力曲线及 毛管力曲线分布的方法， 确定毛管力曲 线的方法包括以下步骤： S1： 进行已有岩心进行 毛管力测试， 获得毛管力曲线； S2： 确定曲线拐 点， 建立拐点处的进汞饱和度和相应毛管力的比 值与岩心渗透率的关系式一， 门槛压力、 拐点处 毛管力以及拐点处进汞饱和度的关系式二， 拐点 处毛管力与岩心渗透率的关系式三， 以及毛管力 与进汞体积 的关系式四； S3： 数据拟合得到关系 式一和三中的特征参数； S4： 推导得到孔吼结构 几何因子、 门槛压力与 岩心渗透率之间的关系式 五； S5： 根据测井数据获得目标井的岩心渗透率， 结合关系 式五， 确定其对应的毛管力曲线。 本发 明能够为油藏数值模拟提供更加准确的毛管力 曲线数据。 权利要求书3页 说明书7页 附图8页 CN 115034112 A 2022.09.09 CN 115034112 A 1.一种基于测井数据确定毛管力曲线的方法， 其特 征在于， 包括以下步骤： S1： 对目标井所处目标区块现有的岩心进行毛管力测试， 获得所述岩心对应取心井的 毛管力曲线； S2： 确定所述毛管力曲线的拐点， 建立所述拐点处的进汞饱和度和相应毛管力的比值 与岩心渗透率的关系式一， 门槛压力、 拐点处毛管力以及拐点处进汞饱和度之间的关系式 二， 拐点处毛管力与岩心渗透率之间的关系式三， 以及毛管力与进汞体积的关系式四； S3： 根据所述毛管力曲线， 结合所述拐点， 数据拟合得到所述关系式一中的特征参数一 和所述关系式三中的特 征参数二； S4： 根据所述关系式一、 关系式二、 关系式三、 关系式四以及步骤S3获得的特征参数， 推 导得到孔吼结构几何 因子、 门槛压力与岩心渗透率之间的关系式五； S5： 根据测井数据获得所述目标井的岩心渗透率， 根据目标井的岩心渗透率， 结合所述 关系式五， 即可确定目标井岩心渗透率对应的毛管力曲线。 2.根据权利要求1所述的基于测井数据确定毛管力曲线的方法， 其特征在于， 步骤S2 中， 所述关系式一 为： K＝a(SHgA/pcA)b                             (1) 式中： K为井位网格点处的岩样渗透率； a和b均为常数， 即所述特征参数一； SHgA为拐点 处的进汞饱和度； pcA为拐点处的毛管压力； 所述关系式二 为： pd＝pcASHgA                               (2) 式中： pd为门槛压力； 所述关系式三 为： pcA＝c/Kd                                (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范围内的所有井位网格点， 并将其权　利　要　求　书 1/3 页 2 CN 115034112 A 2按照距离网格点H(i， j)的距离大小 进行排列； S4'： 根据搜索得到的所有井位网格点的结果， 确定各网格点的毛管力曲线： 若与网格点H(i， j)距 离最近的井位网格点， 其与网格点H(i， j)的网格距 离d1≤R/10， 则 网格点H(i， j)的毛管力曲线直接等于该井位网格点毛管压力曲线； 若搜索半径R范围内的井位网格数n≥3， 则网格点H(i， j)的毛管力采用距离加权平均 公式进行求 解获得； 若搜索半径R范围内的井位网格数n＜3， 则网格点H(i， j)的毛管力采用逼近函数式插 值求得； S5'： 以地质模型纵向分层为基础， 重复步骤S1' ‑S4'， 依次求得所述目标区块纵向第1 层、 第2层、 ……、 第e层的毛管力曲线分布。 5.根据权利要求4所述的基于测井数据确定毛管力曲线分布的方法， 其特征在于， 步骤 S1'中， 坐标xi和yj分别通过 下式进行计算： xi＝Xmin+(Xmax‑Xmin)×(i‑1)/(m‑1) i＝1,2,……,m              (5) yj＝Ymin+(Ymax‑Ymin)×(j‑1)/(m‑1) j＝1,2,……,m              (6) 式中： xi、 yi分别为网格点H(i,j)的横坐标、 纵坐标； Xmin、 Xmax分别为矩阵网格面上横坐 标的最小值、 最大值； Ymin、 Ymax分别为矩阵网格面上纵坐标的最小值、 最大值； i、 j分别 为网 格点H(i,j)所处网格的横向网格数、 纵向网格数； m为矩阵网格面的横向或纵向的网格总 数。 6.根据权利要求4所述的基于测井数据确定毛管力曲线分布的方法， 其特征在于， 步骤 S2'中， 给予权重系数时： 距离 网格点H(i， j)越近的井位网格点(xk， yk)权重越大， 距离网格 点H(i， j)越远的井位网格点(xk， yk)权重越小。 7.根据权利要求6所述的基于测井数据确定毛管力曲线分布的方法， 其特征在于， 步骤 S2'中， 给予权 重系数时， 采用的权函数为： 其中： 式中： w为距 离权重； α 为经验系数； dk为网格点H(i， j)与井位网格点(xk， yk)的距离； xk为 井位网格点(xk， yk)的横坐标； yk为井位网格点(xk， yk)的纵坐标。 8.根据权利要求4所述的基于测井数据确定毛管力曲线分布的方法， 其特征在于， 步骤 S4'中， 所述距离加权平均公式为： 式中： z(i， j)为网格点H(i， j)的毛管力值； n为搜索 半径R范围内的井位网格点总数； pc (i'， j')为与网格点H(i， j)距离最近的井位网格点的毛管力值； wk为第k个井位网格点的距 离权重； k为井位网格点累加数。权　利　要　求　书 2/3 页 3 CN 115034112 A 3