一种获取三相相对渗透率曲线的方法及装置的制造方法
【技术领域】
[0001 ]本申请涉及油气藏开发技术领域，特别涉及一种获取油-气-水三相相对渗透率曲 线的方法及装置。
【背景技术】
[0002] 相对渗透率曲线是油气藏开发中的一项重要数据，它能够精细描述孔隙介质中的 多相渗流特征。目前，油-水或气-液两相相对渗透率曲线主要通过岩心驱替实验获取，计算 方法主要以Johnson-Bossler-Naumann(JBN)等解析方法以及数值反演方法为主。而对获取 油-气-水三相相对渗透率曲线的方法的研究相对较少。
[0003] 随着油藏纵向及平面非均质矛盾的加剧，中国注水开发油田主体已普遍进入高含 水、高采出程度的"双高"开发阶段，剩余油呈现"整体高度分散、局部相对富集"的分布格 局。室内实验及矿场实践表明，气-水交替非混相驱是水驱油藏后期进一步提高原油采收率 的有效技术。因此，研究适用于气-水交替非混相驱油技术的获取油-气-水三相相对渗透率 曲线的方法显得意义重大，这可以为准确描述孔隙介质中气-水交替非混相驱替的多相渗 流特征提供理论依据。
【发明内容】
[0004] 本申请实施例的目的是提供一种获取三相相对渗透率曲线的方法及装置，以实现 获取适用于气-水交替非混相驱油技术的油-气-水三相相对渗透率曲线的目的。
[0005] 为解决上述技术问题，本申请实施例提供一种获取三相相对渗透率曲线的方法及 装置是这样实现的：
[0006] 本申请实施例提供了一种获取三相相对渗透率曲线的方法，包括：
[0007] S1，基于目标区域的气-水交替非混相驱替实验，获取不同时刻的驱替压差、含水 率以及气油比，所述不同时刻的驱替压差、含水率以及气油比构成了动态数据观测值向量；
[0008] S2，基于所构建的表征模型所表征的先验三相相对渗透率曲线以及所构建的组分 数值模拟模型，获取当前动态数据预测值向量；
[0009] S3,利用所述动态数据观测值向量和所述当前动态数据预测值向量，对所构建的 目标函数进行迭代计算，直到所述目标函数的迭代计算满足预设的迭代收敛条件；
[0010] S4,根据满足所述迭代收敛条件的目标函数，获取所述目标区域的油-气-水三相 相对渗透率曲线。
[0011] 在至少一实施例中，所述组分数值模拟模型基于所确定的流体物性参数分布特征 来构建，所述流体物性参数特征包括饱和压力、原油密度、粘度和体积系数的变化规律。
[0012] 在至少一实施例中，所述步骤S2包括：
[0013] S21，获取先验控制参数向量集合，所述先验控制参数向量集合包括多个相对渗透 率先验模型所对应的先验控制参数向量；
[0014] S22,基于所构建的组分数值模拟模型，利用所述先验控制参数向量集合进行组分 数值模拟，以得到当前时间步的所述多个相对渗透率先验模型所对应的先验动态数据预测 值向量；
[0015] S23,计算每个所述先验动态数据预测值向量与所述动态数据观测值向量之间的 差值或误差平方和，以得到符合预设要求的校正控制参数向量集合；
[0016] S24,计算所述校正控制参数向量集合所对应的先验动态数据预测值向量的平均 值，以得到当前动态数据预测值向量，
[0017]相应的，所述步骤S3包括：
[0018] S31，利用所述动态数据观测值向量和所述当前动态数据预测值向量，计算所述目 标函数的当前数值；
[0019] S32,判断当前迭代计算是否满足迭代收敛条件；
[0020] S33,在判断出当前迭代计算不满足迭代收敛条件时，将所述校正控制参数向量集 合作为先验控制参数向量集合，重复步骤S22-S24,得到当前动态数据预测值向量；
[0021 ] S34，在得到所述校正控制参数向量集合所对应的当前动态数据预测值向量后，重 复步骤S31-32,直到得到满足所述迭代收敛条件的目标函数。
[0022] 在至少一实施例中，所构建的目标函数的表达式如下：
[0023]
[0024] 其中，0(m)为目标函数;m为所述表征模型中的控制参数向量;T为表征向量或矩阵 转置的符号;cUs为动态数据观测值向量;g(m)为动态数据预测值向量;C D为动态数据测量 误差的协方差矩阵;mpricir为所述控制参数向量m的先验模型信息的平均值;C m为先验模型的 协方差矩阵。
[0025] 在至少一实施例中，所述根据满足所述迭代收敛条件的目标函数，获取所述目标 区域的油-气-水三相相对渗透率曲线包括：
[0026] 根据满足所述迭代收敛条件的目标函数，确定出所述控制参数向量中各控制参数 的当前数值；
[0027] 根据所述控制参数向量中各控制参数的当前数值，计算油-水两相系统的油相相 对渗透率、水相相对渗透率以及油-气两相系统的气相相对渗透率以及油相相对渗透率；
[0028] 利用所计算的油-水两相系统的油相相对渗透率、水相相对渗透率以及油-气两相 系统的气相相对渗透率以及油相相对渗透率，计算三相渗流条件下的油相相对渗透率。 [0029]在至少一实施例中，通过以下公式来计算三相渗流条件下的油相相对渗透率：
[0030] Kro= (Krow+Krw) (Krog+Krg)_(Krw+Krg)
[0031] 式中，Kr。为三相渗流条件下的油相相对渗透率;Kr?为油-水两相系统的油相相对 渗透率;Krw为油-水两相系统的水相相对渗透率;K rcig为油-气两相系统的油相相对渗透率； Krg为油-气两相系统的气相相对渗透率。
[0032]本申请实施例还提供了一种获取三相相对渗透率曲线的装置，包括：
[0033]第一获取单元，用于基于目标区域的气-水交替非混相驱替实验，获取不同时刻的 驱替压差、含水率以及气油比，所述不同时刻的驱替压差、含水率以及气油比构成了动态数 据观测值向量；
[0034] 第二获取单元，用于基于所构建的表征模型所表征的先验三相相对渗透率曲线以 及所构建的组分数值模拟模型，获取当前动态数据预测值向量；
[0035] 迭代计算单元，用于利用所述动态数据观测值向量和所述当前动态数据预测值向 量，对所构建的目标函数进行迭代计算，直到所述目标函数的迭代计算满足预设的迭代收 敛条件；
[0036] 第三获取单元，用于根据满足所述迭代收敛条件的目标函数，获取所述目标区域 的油-气-水三相相对渗透率曲线。
[0037] 在至少一实施例中，所述第二获取单元包括：
[0038] 获取子单元，用于获取先验控制参数向量集合，所述先验控制参数向量集合包括 多个相对渗透率先验模型所对应的先验控制参数向量；
[0039] 组分数值模拟子单元，用于基于所构建的组分数值模拟模型，利用所述先验控制 参数向量集合进行组分数值模拟，以得到当前时间步的所述多个相对渗透率先验模型所对 应的先验动态数据预测值向量；
[0040] 第一计算子单元，用于计算每个所述先验动态数据预测值向量与所述动态数据观 测值向量之间的差值或误差平方和，以得到校正控制参数向量集合；
[0041] 第二计算子单元，用于计算所述校正控制参数向量集合所对应的先验动态数据预 测值向量的平均值，以得到当前动态数据预测值向量，
[0042]相应的，迭代计算单元包括：
[0043]第三计算子单元，用于利用所述动态数据观测值向量和所述当前动态数据预测值 向量，计算所述目标函数的当前数值；
[0044] 判断子单元，用于判断当前迭代计算是否满足迭代收敛条件；
[0045] 第一控制子单元，用于在判断出当前迭代计算不满足迭代收敛条件时，将所述校 正控制参数向量集合作为先验控制参数向量集合，并控制所述组分数值模拟子单元、所述 第一计算子单元以及所述第二计算子单元进行相应的操作，即