(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211495714.5 (22)申请日 2022.11.28 (71)申请人 中国石油大 学 （华东） 地址 257061 山东省东营市东营区北一路 739号 (72)发明人 王森　王昆　冯其红　李顺明　 杜宜静　吴政　王相　张纪远　 张先敏　 (74)专利代理 机构 北京高沃 律师事务所 1 1569 专利代理师 万慧华 (51)Int.Cl. G06F 30/20(2020.01) E21B 49/00(2006.01) G06F 111/10(2020.01) (54)发明名称 一种水驱油藏储层参数时变规律预测方法 及系统 (57)摘要 本发明涉及一种水驱油藏储层参数时变规 律预测方法及系统， 属于油气田开发技术领域。 首先， 根据实测孔径分布曲线和实测T2谱， 建立 T2谱与孔径关系模型； 获取水驱物理模拟实验过 程中不同注水量下岩样的T2谱； 然后将获得不同 注水量下的T2谱， 转化为不同注水量下的孔径分 布， 通过孔径分布建立不同注水量下的孔隙网络 模型， 开展油水两相渗流模拟， 最终得到不同PV 倍数下岩样所对应的油水两相相对渗透率曲线； 所述油水两相相对渗透率曲线即为所述水驱油 藏储层参数， 从而有效解决了 现有技术中难以同 时测量给定岩样渗透率变化规律和相对渗透率 曲线变化规律的问题， 从而准确预测储层参数时 变规律。 权利要求书3页 说明书15页 附图7页 CN 115526067 A 2022.12.27 CN 115526067 A 1.一种水驱油藏储层参数时变规 律预测方法， 其特 征在于， 所述方法包括： 获取水驱油藏储层岩 样的实测孔径分布曲线和实测T2谱； 根据所述实测孔径分布曲线和所述实测T2谱， 建立T2谱与孔径关系实测模型； 获取水驱物理模拟实验过程中不同注水量下的所述岩 样的T2谱； 利用所述T2谱与孔径关系实测模型， 根据不同注水量下的所述岩样的T2谱， 得到不同注 水量下的孔径分布曲线； 根据不同注水量下的孔径分布曲线， 构建所述孔径分布曲线随岩样注水参数变化的数 学模型； 所述岩样注水参数包括： 注水量、 岩样初始渗透率、 平均注水速度、 初始孔隙度、 泥 质类型、 初始泥质含量和粒度中值； 根据所述数 学模型构建不同注水量下的孔隙网络模型； 利用不同注水量下的孔隙网络模型进行油水两相渗流模拟， 计算不同注水量下岩样所 对应的油水两相 相对渗透率曲线； 所述油水两相 相对渗透率曲线即为水驱油藏储层参数。 2.根据权利要求1所述的水驱油藏储层参数时变规律预测方法， 其特征在于， 所述获取 水驱油藏储层岩 样的实测孔径分布曲线和实测T2谱， 具体包括： 利用高压 压汞方法获取岩 样的实测孔径分布曲线； 利用核磁共 振方法获取岩 样的实测T2谱。 3.根据权利要求1所述的水驱油藏储层参数时变规律预测方法， 其特征在于， 所述根据 所述实测孔径分布曲线和所述实测T2谱， 建立T2谱与孔径关系实测模型， 具体包括： 设置转换系数C； 根据所述实测T2谱和所述转换系数C， 利用T2谱与孔径关系模型， 计算孔径分布曲线； 所 述T2谱与孔径关系模型为： r= CT2， 其中， C为 转换系数； r为 孔隙半径； T2为横向弛豫时间； 判断所述孔径分布曲线和所述实测孔径分布曲线的第一个峰值是否重合， 若不重合， 则调节所述转换系数C的值， 并返回步骤 “根据所述实测T2谱和所述转换系数C， 利用T2谱与 孔径关系模型， 计算孔径分布曲线 ”， 直到所述孔径分布曲线和所述 实测孔径分布曲线的第 一个峰值重合， 得到T2谱与孔径关系实测模型。 4.根据权利要求1所述的水驱油藏储层参数时变规律预测方法， 其特征在于， 所述根据 不同注水量下 的孔径分布曲线， 构建所述孔径分布曲线随岩样注水参数变化的数学模型， 具体包括： 将所述不同注水量下的孔径分布曲线作为若干个正态分布函数的叠加； 正态分布函数 的个数与孔径分布曲线上的峰值个数保持一 致； 计算每个所述正态分布函数的数 学期望和标准差； 构建所述正态分布函数的数学期望、 标准差与岩样注水参数之间的函数关系式， 得到 所述孔径分布曲线随岩 样注水参数变化的数 学模型。 5.根据权利要求4所述的水驱油藏储层参数时变规律预测方法， 其特征在于， 所述计算 每个所述正态分布函数的数 学期望和标准差， 具体包括： 利用EM算法计算每 个所述正态分布函数的数 学期望和标准差 。 6.根据权利要求1所述的水驱油藏储层参数时变规律预测方法， 其特征在于， 在所述获 取水驱物理模拟实验过程中不同注水量下的所述岩 样的T2谱之后， 还 包括： 获取水驱物理模拟实验过程中不同注水量下的所述岩 样的实际渗透率。权　利　要　求　书 1/3 页 2 CN 115526067 A 27.根据权利要求6所述的水驱油藏储层参数时变规律预测方法， 其特征在于， 在所述根 据所述数 学模型构建不同注水量下的孔隙网络模型之后， 还 包括： 利用所述实际渗透率校正所述 孔隙网络模型。 8.根据权利要求7所述的水驱油藏储层参数时变规律预测方法， 其特征在于， 所述利用 所述实际渗透率校正所述 孔隙网络模型， 具体包括： 利用所述 孔隙网络模型进行 单相水的孔隙网络模拟； 获取模拟得到的孔隙网络 两端的压力差和流 量； 根据所述压力差和流 量利用达西公式计算模型渗透率； 比较所述模型渗透率与所述实 际渗透率是否相同； 若不同， 则调节所述孔隙网络模型 的物理参数， 并返回 “利用所述孔隙网络模型进行单相水 的孔隙网络模拟 ”步骤， 直至所述 模型渗透率与所述实际渗透率相同。 9.根据权利要求1所述的水驱油藏储层参数时变规律预测方法， 其特征在于， 所述利用 不同注水量下的孔隙网络模型进 行油水两相渗流模拟， 计算不同注水量下岩样所对应的油 水两相相对渗透率曲线， 具体包括： 进行油水两相渗流模拟， 采用油驱完全饱和水的孔隙网络模型， 直到孔隙网络模型不 出水为止， 建立束缚水状态下的油相有效渗透率 ； ； 式中， 表示油相的粘度； 表示油相通过孔隙网络模型的流量； L表示孔隙网络模型的长度； A表示孔隙网络模型的截面积； ∆P表示孔隙网络模型入口与出口间的压力差； 进行水驱油模拟 实验， 获取油水两相在孔隙网络模型中的流量和孔隙网络模型两端的 压力差， 计算油水两相共存时各相的有效渗透率， 计算公式为： ； ； 其中， 、 分别表示油相和水相的有效渗透率， 表示水相的粘 度； 表示水相通过孔隙网络模型的流 量； 当含水率达到99.95%时或注水量达到30倍PV后， 计算残余油状态下的水相有效渗透 率， 计算方程为： ； 式中： 表示残余油状态下的水相有效渗 透率。 10.一种水驱油藏储层参数时变规 律预测系统， 其特 征在于， 所述系统包括： 初始数据获取模块， 用于获取 水驱油藏储层岩 样的实测孔径分布曲线和实测T2谱； 模型建立模块， 用于根据所述实测孔径分布曲线和所述实测T2谱， 建立T2谱与孔径关系 实测模型； T2谱获取模块， 用于获取 水驱物理模拟实验过程中不同注水量下的所述岩 样的T2谱； 孔径分布曲线建立模块， 用于利用所述T2谱与孔径关系实测模型， 根据不同注水量下的权　利　要　求　书 2/3 页 3 CN 115526067 A 3