一种基于三元动态参数的水驱油相渗时变规律表征方法技术

本发明专利技术公开了一种基于三元动态参数的水驱油相渗时变规律表征方法，包括如下步骤：划分稳定生产阶段；阶段递减分析；确定相渗曲线的端点值与形态；回归不同瞬时水驱强度下残余油饱和度与累积水驱强度的关系；回归残余油下水相相对渗透率与瞬时水驱强度的关系；回归油相指数m、水相指数n分别与累积水驱强度、瞬时水驱强度的关系；联立得到一种基于有效润湿相饱和度、累积水驱强度、瞬时水驱强度的三元动态参数水驱油相渗时变表征通式。经过实际生产效果检验证明，利用该方法可以显著提高中高渗水驱砂岩油藏开发历史拟合精度与剩余油饱和度分布预测精度，进而提高生产动态指标变化趋势预测准确度，具有较强的通用性。

【技术实现步骤摘要】
一种基于三元动态参数的水驱油相渗时变规律表征方法
本专利技术涉及一种油气勘探开发技术中，中高渗水驱砂岩油藏水驱油相渗曲线时变规律表征技术。
技术介绍
水驱条件下，储层的孔隙结构和润湿性发生了较大的变化，因此水驱油效率与水驱油相渗曲线并不是一成不变的。然而目前的表征方法多局限于单独考虑累积水驱强度或瞬时水驱强度对相渗曲线的影响，难以对同时考虑累积水驱强度与瞬时水驱强度对水驱油相渗曲线动态变化规律的影响进行综合表征，无法实时预测开发工作制度调整后的相渗曲线变化，导致油田开发中高含水期剩余油分布预测精度不高。为了实现油田稳油控水、有效调整挖潜与提高原油采收率，需要对水驱油藏相渗曲线动态变化规律有更加深入的认识。然而传统表征方法具备以下固有的缺陷：(1)室内试验所用岩心不一定能够代表地下真实的渗流特征变化规律；(2)只能建立相渗曲线端点值与累积水驱强度或瞬时水驱强度的单因素关系，考虑因素不够全面，并且无法考虑累积水驱强度与瞬时水驱强度对相渗曲线形态变化的影响；(3)类似沉积单元相渗曲线回归分析得到的相渗曲线变化规律地质局限性大，难以推广使用。以上缺陷造成剩余油饱和度分布预测精度不高，无法准确指导开发调整。
技术实现思路
为了克服上述之不足，本专利技术的目的在于提供一种准确评价中高渗水驱砂岩油藏相渗曲线动态变化规律的基于三元动态参数的水驱油相渗时变规律表征方法,该表征方法同时考虑了累积水驱强度与瞬时水驱强度对水驱油相渗曲线端点值以及形态动态变化规律影响。本专利技术所采用的技术方案是：一种基于三元动态参数的水驱油相渗时变规律表征方法，包括如下步骤：S1，基于中高渗水驱砂岩油藏生产数据分析，结合张金庆水驱曲线Wp/Np＝a+b*Wp/Np2，划分稳定生产阶段，即选取Wp/Np2与Wp/Np呈近似直线关系的各阶段动态数据；其中，Wp——累积产水量，104m3；Np——累积产油量，104m3；a——回归系数；b——可动油储量，104m3；S2，结合阶段递减分析，得到各生产阶段的可动油储量，以此计算对应的驱油效率和残余油饱和度；S3，在各稳定生产阶段分别拟合采出程度-含水率关系曲线，得到残余油下水相相对渗透率、油相指数和水相指数，从而确定相渗曲线的端点值与形态；S4，使用水相面通量表征累积水驱强度，使用月产液量表征瞬时水驱强度，将瞬时水驱强度接近的相邻稳定生产阶段组合为单元组，对各单元组分别回归残余油饱和度与水相面通量的关系，得到残余油饱和度与水相面通量、月产液量的通式，不同瞬时水驱强度下残余油饱和度与累积水驱强度明显呈现不同的变化趋势，通式如式(5)所示：Sor＝f1(CS,TS)(5)其中，Sor——残余油饱和度，CS——水相面通量，m；TS——月产液量，m3；S5，回归残余油下水相相对渗透率与瞬时水驱强度的关系，得到各单元组残余油下水相相对渗透率Krw(Sor)与瞬时水驱强度的通式；S6，对各单元组分别回归油相指数m与累积水驱强度、瞬时水驱强度的关系，得到油相指数m是与累积水驱强度、瞬时水驱强度相关的函数；S7，对各单元组分别回归水相指数n与累积水驱强度、瞬时水驱强度的关系，得到水相指数n与累积水驱强度、瞬时水驱强度相关的函数；S8，将步骤S4、S5、S6、S7中的通式联立得到一种基于有效润湿相饱和度、累积水驱强度、瞬时水驱强度的三元动态参数水驱油相渗时变表征通式，从而能够计算得到实时变化的水驱油相渗曲线。进一步地，所述步骤S2中,各生产阶段的可动油储量为以当前生产趋势外推废弃产量达到0时的可动油储量。进一步地，所述步骤S3中,具体步骤如下:对参数a赋一初始值，计算不同平均含水饱和度Sw下经过a、b值校正的水驱油相渗渗透率比值Kro/Krw：其中，Kro——油相相对渗透率；Krw——水相相对渗透率；μo——油相粘度，mPa·s；μw——水相粘度，mPa·s；Bo——油相体积系数；Bw——水相体积系数；N——地质储量，104m3；Sw——平均含水饱和度；Swc——束缚水饱和度；结合水驱油相渗曲线Willhite幂函数经验通式：其中，Sor——残余油饱和度；Kro(Swc)——束缚水饱和度下油相相对渗透率；Krw(Sor)——残余油饱和度下水相相对渗透率；m——油相指数；n——水相指数；式(2)除以式(3)，两边取对数，得式(4)：y＝mx1-nx2+t(4)其中通常Kro(Swc)数值取1，由式(4)经二元线性回归，同时通过调整参数a值，在各稳定生产阶段分别拟合采出程度-含水率关系曲线，得到残余油下水相相对渗透率Krw(Sor)、油相指数m和水相指数n，从而确定相渗曲线的端点值与形态。进一步地，所述步骤S5中，做出各单元组残余油下水相相对渗透率Krw(Sor)与瞬时水驱强度的通式如式(7)所示：Krw(Sor)＝f2(TS)(7)TS——月产液量，使用月产液量表征瞬时水驱强度。进一步地，所述步骤S6中，在瞬时水驱强度接近的各单元组内，分别回归油相指数m与累积水驱强度的关系式，通式如式(9)所示：m＝f3(CS)(9)特别的，当油相指数m与累积水驱强度呈线性关系式时，有：m＝f3(CS)＝J·(CS)+K(10)其中，斜率J、截距K为回归系数；分别回归斜率J、截距K与瞬时水驱强度的关系式，通式如式(11)、式(12)所示：J＝f4(TS)(11)K＝f5(TS)(12)特别的，当斜率J、截距K与瞬时水驱强度呈线性关系式时，有：J＝f4(TS)＝V·(TS)+W(13)K＝f5(TS)＝X·(TS)+Y(14)油相指数m是与累积水驱强度、瞬时水驱强度相关的函数：m＝f6(CS,TS)(15)其中，V、W、X、Y为回归系数。进一步地，所述步骤S7中，在瞬时水驱强度接近的各单元组内，分别回归水相指数n与累积水驱强度的关系式，通式如式(16)所示：n＝f7(CS)(16)特别的，当水相指数n与累积水驱强度呈线性关系式时，有：n＝f7(CS)＝-P·(CS)+Q(17)其中，斜率P、截距Q为回归系数；分别回归斜率P、截距Q与瞬时水驱强度的关系式(图11、图12)，通式如式(18)、式(19)所示：P＝f8(TS)(18)Q＝f9(TS)(19)；特别的，当斜率P、截距Q与瞬时水驱强度呈线性关系式时，有：P＝f8(TS)＝A·(TS)+B(20)Q＝f9(TS)＝R·(TS)+U(21)水相指数n是与累积水驱强度、瞬时水驱强度相关的函数：n＝f10(CS,TS)(22)其中，A、B、R、U为回归系数。进一步地，所述步骤S8中，联立通式(5)、(7)、(15)、(22)，得到水驱油相渗曲线端点值与形态动态变化表征通式，如式(23)、式(24)所示:对于水湿储层，有效润湿相饱和度为：对应的基于三元动态参数(有效润湿相饱和度、累积水驱强度、瞬时水驱强度)的水湿储层水驱油相渗曲线时变规律表征通式为：其中Krnwp、Krwp分别为非湿相、润湿相相对渗透率；对于油湿储层，有效润湿相饱和度为：对应的基于三元动态参数(有效润湿相饱和度、累积水驱强度、瞬时水驱强度)的油湿储层水驱油相渗曲线时变规律表征通式为：本专利技术的有益效果：经过实际生产效果检验证明，利用该方法可以显著提高中高渗水驱砂岩油藏开发历史拟合精度与剩余油饱和度分布预测精度，进而提高生产动态本文档来自技高网...

【技术保护点】
1.一种基于三元动态参数的水驱油相渗时变规律表征方法，其特征在于,包括如下步骤：S1，基于中高渗水驱砂岩油藏生产数据分析，结合张金庆水驱曲线Wp/Np＝a+b*Wp/Np2，划分稳定生产阶段，即选取Wp/Np2与Wp/Np呈近似直线关系的各阶段动态数据；其中，Wp——累积产水量，104m3；Np——累积产油量，104m3；a——回归系数；b——可动油储量，104m3；S2，结合阶段递减分析，得到各生产阶段的可动油储量，以此计算对应的驱油效率和残余油饱和度；S3，在各稳定生产阶段分别拟合采出程度‑含水率关系曲线，得到残余油下水相相对渗透率、油相指数和水相指数，从而确定相渗曲线的端点值与形态；S4，使用水相面通量表征累积水驱强度，使用月产液量表征瞬时水驱强度，将瞬时水驱强度接近的相邻稳定生产阶段组合为单元组，对各单元组分别回归残余油饱和度与水相面通量的关系，得到残余油饱和度与水相面通量、月产液量的通式，不同瞬时水驱强度下残余油饱和度与累积水驱强度明显呈现不同的变化趋势，通式如式(5)所示：Sor＝f1(CS,TS)       (5)，其中，Sor——残余油饱和度，CS——水相面通量，m；TS——月产液量，m3；S5，回归残余油下水相相对渗透率与瞬时水驱强度的关系，得到各单元组残余油下水相相对渗透率Krw(Sor)与瞬时水驱强度的通式；S6，对各单元组分别回归油相指数m与累积水驱强度、瞬时水驱强度的关系，得到油相指数m是与累积水驱强度、瞬时水驱强度相关的函数；S7，对各单元组分别回归水相指数n与累积水驱强度、瞬时水驱强度的关系，得到水相指数n与累积水驱强度、瞬时水驱强度相关的函数；S8，将步骤S4、S5、S6、S7中的通式联立得到一种基于有效润湿相饱和度、累积水驱强度、瞬时水驱强度的三元动态参数水驱油相渗时变表征通式，从而能够计算得到实时变化的水驱油相渗曲线。...

【技术特征摘要】
1.一种基于三元动态参数的水驱油相渗时变规律表征方法，其特征在于,包括如下步骤：S1，基于中高渗水驱砂岩油藏生产数据分析，结合张金庆水驱曲线Wp/Np＝a+b*Wp/Np2，划分稳定生产阶段，即选取Wp/Np2与Wp/Np呈近似直线关系的各阶段动态数据；其中，Wp——累积产水量，104m3；Np——累积产油量，104m3；a——回归系数；b——可动油储量，104m3；S2，结合阶段递减分析，得到各生产阶段的可动油储量，以此计算对应的驱油效率和残余油饱和度；S3，在各稳定生产阶段分别拟合采出程度-含水率关系曲线，得到残余油下水相相对渗透率、油相指数和水相指数，从而确定相渗曲线的端点值与形态；S4，使用水相面通量表征累积水驱强度，使用月产液量表征瞬时水驱强度，将瞬时水驱强度接近的相邻稳定生产阶段组合为单元组，对各单元组分别回归残余油饱和度与水相面通量的关系，得到残余油饱和度与水相面通量、月产液量的通式，不同瞬时水驱强度下残余油饱和度与累积水驱强度明显呈现不同的变化趋势，通式如式(5)所示：Sor＝f1(CS,TS)(5)，其中，Sor——残余油饱和度，CS——水相面通量，m；TS——月产液量，m3；S5，回归残余油下水相相对渗透率与瞬时水驱强度的关系，得到各单元组残余油下水相相对渗透率Krw(Sor)与瞬时水驱强度的通式；S6，对各单元组分别回归油相指数m与累积水驱强度、瞬时水驱强度的关系，得到油相指数m是与累积水驱强度、瞬时水驱强度相关的函数；S7，对各单元组分别回归水相指数n与累积水驱强度、瞬时水驱强度的关系，得到水相指数n与累积水驱强度、瞬时水驱强度相关的函数；S8，将步骤S4、S5、S6、S7中的通式联立得到一种基于有效润湿相饱和度、累积水驱强度、瞬时水驱强度的三元动态参数水驱油相渗时变表征通式，从而能够计算得到实时变化的水驱油相渗曲线。2.根据权利要求1所述的基于三元动态参数的水驱油相渗时变规律表征方法，其特征在于：所述步骤S2中,各生产阶段的可动油储量为以当前生产趋势外推废弃产量达到0时的可动油储量。3.根据权利要求1所述的基于三元动态参数的水驱油相渗时变规律表征方法，其特征在于,所述步骤S3中,具体步骤如下:对参数a赋一初始值，计算不同平均含水饱和度Sw下经过a、b值校正的水驱油相渗渗透率比值Kro/Krw：，其中，Kro——油相相对渗透率；Krw——水相相对渗透率；μo——油相粘度，mPa·s；μw——水相粘度，mPa·s；Bo——油相体积系数；Bw——水相体积系数；N——地质储量，104m3；Sw——平均含水饱和度；Swc——束缚水饱和度；结合水驱油相渗曲线Willhite幂函数经验通式：，其中，Sor——残余油饱和度；Kro(Swc)——束缚水饱和度下油相相对渗透率；Krw(Sor)——残余油饱和度下水相相对渗透率；m——油相指数；n——水相指数；式(2)除以式(3)，两边取对数，得式(4)：y＝mx1-nx2+t(4)，其中通常Kro(Swc)数值取1，由式(4...

【专利技术属性】
技术研发人员：雷霄，王雯娟，王世朝，张乔良，黄冬梅，薛国庆，洪楚侨，
申请(专利权)人：中国海洋石油集团有限公司，中海石油中国有限公司湛江分公司，
类型：发明
国别省市：北京,11