《径向流条件下油水相对渗透率曲线的计算方法.pdf》由会员分享，可在线阅读，更多相关《径向流条件下油水相对渗透率曲线的计算方法.pdf（11页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102915406 A (43)申请公布日 2013.02.06 C N 1 0 2 9 1 5 4 0 6 A *CN102915406A* (21)申请号 201210418648.1 (22)申请日 2012.10.26 G06F 19/00(2006.01) (71)申请人中国石油大学（华东） 地址 266580 山东省青岛市经济技术开发区 长江西路66号 (72)发明人侯健 李振泉 王代刚 邴绍献 罗福全 杜庆军 于波 苏映宏 (74)专利代理机构青岛高晓专利事务所 37104 代理人杨大兴 (54) 发明名称 径向流条件下油水相对渗透率曲线的计算方 法 (。
2、57) 摘要 本发明公开了一种径向流条件下油水相 对渗透率曲线的计算方法，特征是在径向流岩 心驱替实验的基础上，选取不同时刻的驱替压 差和累积产量作为动态数据建立最小二乘目 标函数，利用三次均匀B样条模型表征油水相 对渗透率曲线，结合径向流数值模拟器，通过 Levenberg-Marquardt算法不断调整相渗表征模 型的控制参数向量，使动态数据观测值与预测值 的误差平方和达到最小，最终反演计算油水相对 渗透率曲线。本发明实用性强，在处理径向流岩心 驱替实验数据的基础上可对油水相对渗透率曲线 进行准确计算，为开展油、水两相流体在多孔介质 中的渗流机理研究提供了有效的工具。 (51)Int.Cl。
3、. 权利要求书2页 说明书4页 附图4页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 4 页 附图 4 页 1/2页 2 1.一种径向流条件下油水相对渗透率曲线的计算方法，其特征在于：在径向流岩心驱 替实验的基础上，选取不同时刻的驱替压差和累积产量作为动态数据建立最小二乘目标 函数，利用三次均匀B样条模型表征油水相对渗透率曲线，结合径向流数值模拟器，通过 Levenberg-Marquardt算法不断调整相渗表征模型的控制参数向量，使动态数据观测值与 预测值的误差平方和达到最小，最终反演计算油水相对渗透率曲线，具体实现过程如下： 步骤1：通过径向流岩。
4、心驱替实验采集不同时刻的驱替压差、累积产油量以及累积产 水量数据； 步骤2：基于模型预测值应与实际观测值相吻合的理论，选取不同时刻驱替压差和累 积产量作为动态数据建立最小二乘目标函数，表达式为： 式中，O(m)为目标函数；m为m1阶模型参数向量；T为表征向量或矩阵转置的符号； d obs 为n1阶动态数据观测值向量；g(m)为n1阶动态数据预测值向量；C D 为nn阶权 重协方差矩阵； 步骤3：利用三次均匀B样条模型表征油水相对渗透率曲线； 步骤4：结合径向流数值模拟器，利用Levenberg-Marquardt算法不断调整拟控制参数 向量x和y，每次迭代结束，通过求解上述线性方程组计算控制参。
5、数向量C w 和C o ，得到满足 单调上凸性要求的三次均匀B样条形式的油水相对渗透率曲线，不断循环迭代，使动态数 据观测值与预测值的误差平方和达到最小，最终反演计算油水相对渗透率曲线。 2.根据权利要求1所述的径向流条件下油水相对渗透率曲线的计算方法，其特征在 于，所述步骤3进一步包括： 步骤3a.定义无因次含水饱和度，构建三次均匀B样条形式的油、水相对渗透率， 无因次含水饱和度为 相应的三次均匀B样条形式为 l=o,w 式中，S wD 无因次含水饱和度；S wc 、S or 分别为束缚水饱和度、残余油饱和度；k rl 为l相 的相对渗透率；n为控制节点个数；分别为油、水相相对渗透率的控制节。
6、点；B j， p (S wD )为四阶(3次)的B样条基函数； 步骤3b.特殊控制节点处理：室内处理相渗实验数据时，通常将束缚水饱和度下的 油相渗透率作为基准渗透率，并假定束缚水饱和度和残余油饱和度已知，因此有 及成立；为保证三次均匀B样条形式的油、水相对渗透率曲线通过端点和 需引入映射点和并建立以下关系： 且 步骤3c.保证油水相对渗透率曲线的单调上凸性：基于控制节点的单调上凸性等价 于对应B样条曲线单调上凸性的理论，引入对数变换将控制参数向量转化为拟控制参数向 量，以保证油水相对渗透率曲线的单调上凸性： 权 利 要 求 书CN 102915406 A 2/2页 3 对于水相相对渗透率曲线 。
7、1in 对于油相相对渗透率曲线 1in-1。 权 利 要 求 书CN 102915406 A 1/4页 4 径向流条件下油水相对渗透率曲线的计算方法 技术领域 ： 0001 本发明涉及石油开发领域，特别是对径向流条件下油水相对渗透率曲线的计算方 法。 背景技术 ： 0002 油水相对渗透率曲线是油藏开发中的重要数据，它能够反映油、水相流体在多孔 介质中的渗流规律。目前，油水相对渗透率曲线主要通过单向流岩心驱替实验获取，计算方 法以JBN等解析方法为主，而解析方法由于假设条件较为理想，往往导致相对渗透率曲线 的计算精度较低。基于自动历史拟合技术的数值反演方法近十年来发展较快。与解析方法 相比，数。
8、值反演方法应用于室内时，既可综合利用见水前、后动态数据，又能够考虑毛细管 压力、非均质性等因素对相渗反演结果的影响，计算得到的相对渗透率曲线完整且精度较 高。 0003 油田开发实践表明，实际油藏在近井地带以径向渗流为主，单向流岩心驱替实验 中流体为线性渗流，两者渗流方式不同，导致将目前相渗曲线计算理论向矿场推广应用时 存在偏差，影响应用效果。因此，开展径向流岩心驱替实验，研究径向流条件下油水相对渗 透率曲线的计算方法具有重要意义。但由于径向流岩心驱替实验操作难度大、耗费时间长， 流体径向渗流的理论研究不足，使得现有的径向流油水相对渗透率曲线计算方法较少。 发明内容 ： 0004 本发明要解决。
9、的技术问题是提供一种径向流条件下油水相对渗透率曲线的计算 方法，用于准确计算径向流油水相对渗透率曲线。 0005 本发明解决其技术问题所采取的技术方案是在径向流岩心驱替实验的基础上，选 取不同时刻的驱替压差和累积产量作为动态数据建立最小二乘目标函数，利用三次均匀B 样条模型表征油水相对渗透率曲线，结合径向流数值模拟器，通过Levenberg-Marquardt 算法不断调整相渗表征模型的控制参数向量，使动态数据观测值与预测值的误差平方和达 到最小，最终反演计算油水相对渗透率曲线，具体按如下步骤实现： 0006 步骤1：通过径向流岩心驱替实验采集不同时刻的驱替压差、累积产油量以及累 积产水量数据。
10、； 0007 步骤2：基于模型预测值应与实际观测值相吻合的理论，选取不同时刻压力和累 积产量作为动态数据建立最小二乘目标函数，表达式为： 0008 0009 式中，O(m)为目标函数；m为m1阶模型参数向量；T为表征向量或矩阵转置的符 号；d obs 为n1阶动态数据观测值向量；g(m)为n1阶动态数据预测值向量；C D 为nn阶 权重协方差矩阵； 0010 步骤3：利用三次均匀B样条模型表征油水相对渗透率曲线。该步骤还进一步包 说 明 书CN 102915406 A 2/4页 5 括： 0011 步骤3a.定义无因次含水饱和度，构建三次均匀B样条形式的油、水相对渗透率： 0012 无因次含水。
11、饱和度为 0013 0014 相应的三次均匀B样条形式为 0015 l=o,w 0016 式中，S w D为无因次含水饱和度；S wc 、S or 分别为束缚水饱和度、残余油饱和度；k rl 为 l相的相对渗透率；n为控制节点个数；分别为油、水相相对渗透率的控制节点； B j，p (S w D)为四阶(3次)的B样条基函数； 0017 步骤3b.特殊控制节点处理：室内处理相渗实验数据时，通常将束缚水饱和度下 的油相渗透率作为基准渗透率，并假定束缚水饱和度和残余油饱和度已知，因此有 及成立。为保证三次均匀B样条形式的油、水相对渗透率曲线通过端点和 需引入映射点和并建立以下关系： 0018 且 0。
12、019 步骤3c.保证油水相对渗透率曲线的单调上凸性：基于控制节点的单调上凸性等 价于对应B样条曲线单调上凸性的理论，引入对数变换将控制参数向量转换为拟控制参数 向量，以保证油水相对渗透率曲线的单调上凸性： 0020 对于水相相对渗透率曲线 0021 1in 0022 对于油相相对渗透率曲线 0023 1in-1。 0024 步骤4：结合径向流数值模拟器，利用Levenberg-Marquardt算法不断调整拟控制 参数向量x和y，每次迭代结束，通过求解线性方程组计算控制参数向量C w 和C o ，并生成满 足单调上凸性要求的三次均匀B样条形式的油水相对渗透率曲线，不断循环迭代，使动态 数据观。
13、测值与预测值的误差平方和达到最小，最终反演计算油水相对渗透率曲线。 0025 本发明与现有技术相比具体的有益效果是：在径向流岩心驱替实验的基础上，引 入三次均匀B样条模型作为油水相对渗透率曲线表征模型，通过Levenberg-Marquardt算 法对驱替压差和累积产量数据进行自动历史拟合，实现了径向流油水相对渗透率曲线的准 确计算，克服了传统方法的缺陷。 0026 下面结合附图和实施例来进一步描述本发明。 附图说明 ： 说 明 书CN 102915406 A 3/4页 6 0027 图1为径向流岩心驱替实验装置原理图。 0028 图2为圆盘状岩样中的径向渗流特征。 0029 图3为三次均匀B。
14、样条模型的局部拟合特性示意图。 0030 图4为基于三次均匀B样条模型生成油水相对渗透率曲线的效果展示图。 0031 图5为驱替压差的拟合效果图。 0032 图6为累积产油量的拟合效果图。 0033 图7为累积产水量的拟合效果图。 0034 图8为估算得到的径向流油水相对渗透率曲线。 0035 图9为径向流条件下油水相对渗透率曲线计算方法的流程框图。 具体实施方式 ： 0036 以下通过对某油田天然水湿砂岩油藏取心圆盘状岩样开展径向流岩心驱替实验 并结合附图来详细说明本发明的内容和实现原理。 0037 1.开展径向流岩心驱替实验，采集相渗实验数据 0038 按图1所示设计径向流岩心驱替实验装置。
15、，对某油田天然水湿砂岩油藏取心圆盘 状岩样开展径向流岩心驱替实验，其径向渗流特征如图2所示，模拟圆盘状岩样的基本参 数如表1所示。实验过程中，选用的注采方式为中间注、四周采，实验控制条件为定液注入、 定压产出。径向流岩心驱替实验的具体步骤为： 0039 将岩样清洗、烘干，测定岩样孔隙度和气测渗透率，并利用氮气法测定孔隙体 积； 0040 将烘干后的岩样抽真空饱和模拟地层水； 0041 将圆盘状岩样放置在图1所示的驱替系统上，首先用油驱水实验模拟油藏的形 成(驱替过程)并测算束缚水饱和度，然后选择合适的驱替速度或驱替压差进行水驱油实 验(吸吮过程)，准确记录不同时刻的累积产油量、累积产水量和岩样。
16、两端的驱替压差，当 含水率达到99.95%时或注水30倍孔隙体积后，测算残余油饱和度，结束实验； 0042 表1模拟圆盘状岩样基本参数 0043 0044 2.建立最小二乘目标函数 0045 综合采用不同时刻的驱替压差和累积产油量、累积产水量作为动态数据，建立最 小二乘目标函数，表征动态数据观测值与模型预测值的误差平方和。 0046 3.利用三次均匀B样条模型表征油水相对渗透率曲线 说 明 书CN 102915406 A 4/4页 7 0047 由于三次均匀B样条模型具有如图3所示的局部拟合特性，使得它更具有普遍意 义，自由度高，故选取三次均匀B样条模型表征油水相对渗透率曲线，具体包括： 00。
17、48 定义无因次含水饱和度，构建三次均匀B样条形式的油、水相对渗透率； 0049 特殊控制节点处理：将束缚水饱和度下的油相渗透率作为基准渗透率，并假定 束缚水饱和度和残余油饱和度已知，并引入两个控制节点建立映射关系使得三次均匀B样 条形式的油水相对渗透率曲线通过控制参数向量的端点； 0050 引入对数变换将控制参数向量转换为拟控制参数向量，保证油水相对渗透率曲 线的单调上凸性； 0051 最终通过三次均匀B样条模型描述油水相对渗透率曲线，效果如图4所示。 0052 4.径向流条件下油水相对渗透率曲线的计算 0053 结合径向流数值模拟器，利用Levenberg-Marquardt算法不断调整三。
18、次均匀B样 条模型的拟控制参数向量，每次迭代结束，通过求解线性方程组计算相应的控制参数向量， 并生成满足单调上凸性要求的油水相对渗透率曲线，不断循环迭代，直至动态数据观测值 与模型预测值的误差平方和达到最小，进而反演计算径向流油水相对渗透率曲线。驱替压 差、累积产油量和累积产水量的拟合效果如图5图7所示，估算得到的径向流油水相对渗 透率曲线如图8所示。 0054 图9为本发明径向流条件下油水相对渗透率曲线计算方法的流程框图。本发明基 于自动历史拟合技术，建立了径向流油水相对渗透率曲线计算方法，可有效地用于油、水相 流体在多孔介质中的渗流机理研究。 说 明 书CN 102915406 A 1/4页 8 图1 图2 图3 说 明 书 附 图CN 102915406 A 2/4页 9 图4 图5 说 明 书 附 图CN 102915406 A 3/4页 10 图6 图7 说 明 书 附 图CN 102915406 A 10 4/4页 11 图8 图9 说 明 书 附 图CN 102915406 A 11 。