[发明专利]径向流条件下油水相对渗透率曲线的计算方法

摘要

本发明公开了一种径向流条件下油水相对渗透率曲线的计算方法，特征是在径向流岩心驱替实验的基础上，选取不同时刻的驱替压差和累积产量作为动态数据建立最小二乘目标函数，利用三次均匀B样条模型表征油水相对渗透率曲线，结合径向流数值模拟器，通过Levenberg-Marquardt算法不断调整相渗表征模型的控制参数向量，使动态数据观测值与预测值的误差平方和达到最小，最终反演计算油水相对渗透率曲线。本发明实用性强，在处理径向流岩心驱替实验数据的基础上可对油水相对渗透率曲线进行准确计算，为开展油、水两相流体在多孔介质中的渗流机理研究提供了有效的工具。

主权项

一种径向流条件下油水相对渗透率曲线的计算方法，其特征在于：在径向流岩心驱替实验的基础上，选取不同时刻的驱替压差和累积产量作为动态数据建立最小二乘目标函数，利用三次均匀B样条模型表征油水相对渗透率曲线，结合径向流数值模拟器，通过Levenberg‑Marquardt算法不断调整相渗表征模型的控制参数向量，使动态数据观测值与预测值的误差平方和达到最小，最终反演计算油水相对渗透率曲线，具体实现过程如下：步骤1：通过径向流岩心驱替实验采集不同时刻的驱替压差、累积产油量以及累积产水量数据；步骤2：基于模型预测值应与实际观测值相吻合的理论，选取不同时刻驱替压差和累积产量作为动态数据建立最小二乘目标函数，表达式为： O ( m ) = 1 2 ( g ( m ) - d obs ) T C D - 1 ( g ( m ) - d obs ) 式中，O(m)为目标函数；m为m×1阶模型参数向量；T为表征向量或矩阵转置的符号；dobs为n×1阶动态数据观测值向量；g(m)为n×1阶动态数据预测值向量；CD为n×n阶权重协方差矩阵；步骤3：利用三次均匀B样条模型表征油水相对渗透率曲线；步骤4：结合径向流数值模拟器，利用Levenberg‑Marquardt算法不断调整拟控制参数向量x和y，每次迭代结束，通过求解上述线性方程组计算控制参数向量Cw和Co，得到满足单调上凸性要求的三次均匀B样条形式的油水相对渗透率曲线，不断循环迭代，使动态数据观测值与预测值的误差平方和达到最小，最终反演计算油水相对渗透率曲线。