>【主权项】
1. 一种获取=相相对渗透率曲线的方法，其特征在于，包括： 51, 基于目标区域的气-水交替非混相驱替实验，获取不同时刻的驱替压差、含水率W 及气油比，所述不同时刻的驱替压差、含水率W及气油比构成了动态数据观测值向量； 52, 基于所构建的表征模型所表征的先验=相相对渗透率曲线W及所构建的组分数值 模拟模型，获取当前动态数据预测值向量； 53, 利用所述动态数据观测值向量和所述当前动态数据预测值向量，对所构建的目标 函数进行迭代计算，直到所述目标函数的迭代计算满足预设的迭代收敛条件； 54, 根据满足所述迭代收敛条件的目标函数，获取所述目标区域的油-气-水=相相对 渗透率曲线。2. 根据权利要求1所述的方法，其特征在于，所述组分数值模拟模型基于所确定的流体 物性参数分布特征来构建，所述流体物性参数特征包括饱和压力、原油密度、粘度和体积系 数的变化规律。3. 根据权利要求1所述的方法，其特征在于，所述步骤S2包括： S21，获取先验控制参数向量集合，所述先验控制参数向量集合包括多个相对渗透率先 验模型所对应的先验控制参数向量； 522, 基于所构建的组分数值模拟模型，利用所述先验控制参数向量集合进行组分数值 模拟，W得到当前时间步的所述多个相对渗透率先验模型所对应的先验动态数据预测值向 量； 523, 计算每个所述先验动态数据预测值向量与所述动态数据观测值向量之间的差值 或误差平方和，W得到符合预设要求的校正控制参数向量集合； 524, 计算所述校正控制参数向量集合所对应的先验动态数据预测值向量的平均值，W 得到当前动态数据预测值向量， 相应的，所述步骤S3包括： S31，利用所述动态数据观测值向量和所述当前动态数据预测值向量，计算所述目标函 数的当前数值； S32，判断当前迭代计算是否满足迭代收敛条件； 533, 在判断出当前迭代计算不满足迭代收敛条件时，将所述校正控制参数向量集合作 为先验控制参数向量集合，重复步骤S22-S24,得到当前动态数据预测值向量； 534, 在得到所述校正控制参数向量集合所对应的当前动态数据预测值向量后，重复步 骤S31-32,直到得到满足所述迭代收敛条件的目标函数。4. 根据权利要求1或3所述的方法，其特征在于，所构建的目标函数的表达式如下：其中，0(m)为目标函数;m为所述表征模型中的控制参数向量;T为表征向量或矩阵转置 的符号;cUs为动态数据观测值向量;g(m)为动态数据预测值向量;Cd为动态数据测量误差 的协方差矩阵;mprinr为所述控制参数向量m的先验模型信息的平均值;Cm为先验模型的协方 差矩阵。5. 根据权利要求4所述的方法，其特征在于，所述根据满足所述迭代收敛条件的目标函 数，获取所述目标区域的油-气-水=相相对渗透率曲线包括： 根据满足所述迭代收敛条件的目标函数，确定出所述控制参数向量中各控制参数的当 前数值； 根据所述控制参数向量中各控制参数的当前数值，计算油-水两相系统的油相相对渗 透率、水相相对渗透率W及油-气两相系统的气相相对渗透率W及油相相对渗透率； 利用所计算的油-水两相系统的油相相对渗透率、水相相对渗透率W及油-气两相系统 的气相相对渗透率W及油相相对渗透率，计算=相渗流条件下的油相相对渗透率。6. 根据权利要求5所述的方法，其特征在于，通过W下公式来计算=相渗流条件下的油 相相对渗透率： Kro =化iw+Krw K Krog+Krg ) - ( Krw+Krg ) 式中，私。为；相渗流条件下的油相相对渗透率;Kr?为油-水两相系统的油相相对渗透 率;Krw为油-水两相系统的水相相对渗透率;KrDg为油-气两相系统的油相相对渗透率;Krg为 油-气两相系统的气相相对渗透率。 7 . -种获取=相相对渗透率曲线的装置，其特征在于，包括： 第一获取单元，用于基于目标区域的气-水交替非混相驱替实验，获取不同时刻的驱替 压差、含水率W及气油比，所述不同时刻的驱替压差、含水率W及气油比构成了动态数据观 测值向量； 第二获取单元，用于基于所构建的表征模型所表征的先验=相相对渗透率曲线W及所 构建的组分数值模拟模型，获取当前动态数据预测值向量； 迭代计算单元，用于利用所述动态数据观测值向量和所述当前动态数据预测值向量， 对所构建的目标函数进行迭代计算，直到所述目标函数的迭代计算满足预设的迭代收敛条 件； 第=获取单元，用于根据满足所述迭代收敛条件的目标函数，获取所述目标区域的油-气-水=相相对渗透率曲线。8.根据权利要求7所述的装置，其特征在于，所述第二获取单元包括： 获取子单元，用于获取先验控制参数向量集合，所述先验控制参数向量集合包括多个 相对渗透率先验模型所对应的先验控制参数向量； 组分数值模拟子单元，用于基于所构建的组分数值模拟模型，利用所述先验控制参数 向量集合进行组分数值模拟，W得到当前时间步的所述多个相对渗透率先验模型所对应的 先验动态数据预测值向量； 第一计算子单元，用于计算每个所述先验动态数据预测值向量与所述动态数据观测值 向量之间的差值或误差平方和，W得到校正控制参数向量集合； 第二计算子单元，用于计算所述校正控制参数向量集合所对应的动态数据预测值向量 的平均值，W得到当前动态数据预测值向量， 相应的，所述迭代计算单元包括： 第=计算子单元，用于利用所述动态数据观测值向量和所述当前动态数据预测值向 量，计算所述目标函数的当前数值； 判断子单元，用于判断当前迭代计算是否满足迭代收敛条件； 第一控制子单元，用于在判断出当前迭代计算不满足迭代收敛条件时，将所述校正控 制参数向量集合作为先验控制参数向量集合，并控制所述组分数值模拟子单元、所述第一 计算子单元W及所述第二计算子单元进行相应的操作； 第二控制子单元，用于控制所述第=计算子单元W及所述判断子单元的操作，直到得 到满足迭代收敛条件的目标函数。9. 根据权利要求7或8所述的装置，其特征在于，所述迭代计算单元中所构建的目标函 数的表达式如下：其中，0(m)为目标函数;m为所述表征模型中的控制参数向量;T为表征向量或矩阵转置 的符号;cUs为动态数据观测值向量;g(m)为动态数据预测值向量;Cd为动态数据测量误差 的协方差矩阵;mprinr为所述控制参数向量m的先验模型信息的平均值;Cm为先验模型的协方 差矩阵。10. 根据权利要求9所述的装置，其特征在于，所述第=获取单元包括： 确定子单元，用于根据所述符合迭代收敛条件的目标函数，确定出所述控制参数向量 中各控制参数的当前数值. 第四计算子单元，用于根据所述控制参数向量中各控制参数的当前数值，计算油-水两 相系统的油相相对渗透率、水相相对渗透率W及油-气两相系统的气相相对渗透率W及油 相相对渗透率； 第五计算子单元，用于利用所计算的油-水两相系统的油相相对渗透率、水相相对渗透 率W及油-气两相系统的气相相对渗透率W及油相相对渗透率，计算=相渗流条件下的油 相相对渗透率。
【专利摘要】本申请公开了一种获取三相相对渗透率曲线的方法及装置。该方法包括：基于目标区域的气-水交替非混相驱替实验，获取不同时刻的驱替压差、含水率以及气油比，所述不同时刻的驱替压差、含水率以及气油比构成了动态数据观测值向量；基于所构建的表征模型所表征的先验三相相对渗透率曲线以及所构建的组分数值模拟模型，获取当前动态数据预测值向量；利用所述动态数据观测值向量和所述动态数据预测值向量，对所构建的目标函数进行迭代计算，直到所述目标函数的迭代计算满足预设的迭代收敛条件；根据满足所述迭代收敛条件的目标函数，获取目标区域的油-气-水三相相对渗透率曲线。通过本申请所公开的技术方案，可以获取适用于气-水交替非混相驱油技术的三相相对渗透率曲线。
【IPC分类】G06F19/00
【公开号】CN105653845
【申请号】
【发明人】王代刚, 李勇, 胡永乐, 李保柱, 侯健, 刘晓彤, 王伟, 王 琦, 谭柱
【申请人】中国石油天然气股份有限公司
【公开日】2016年6月8日
【申请日】2015年12月24日