CN104727790A - 水驱油藏剩余油再富集成藏周期预测方法 - Google Patents

Abstract

本发明提供一种水驱油藏剩余油再富集成藏周期预测方法，该水驱油藏剩余油再富集成藏周期预测方法包括步骤1，求取所研究油藏的含水饱和度值；步骤2，根据该含水饱和度值选择相应的水驱油藏油水运移速度公式；以及步骤3，利用所选择的该水驱油藏油水运移速度公式计算剩余油再富集成藏需要的最长时间即剩余油再富集成藏周期。该水驱油藏剩余油再富集成藏周期预测方法为特高含水阶段低效开发油藏的战略调整、周期注采开发油藏开发、近废弃油藏的产能恢复及再度开发等方面指明了方向，其推广应用前景广阔，经济社会效益显著。

Description

水驱油藏剩余油再富集成藏周期预测方法

技术领域

本发明涉及水驱油藏开发后期的挖潜与提高采收率领域，特别是涉及到一种水驱油藏剩余油再富集成藏周期预测方法。

背景技术

剩余油一般是指经一、二次采油后，油藏进入高含水的开发后期，油藏剩余的可动资源量，是一个动态的变化量。胜利油田已进入高含水-特高含水开发期。随着开采程度加深，剩余油分布更加分散，油田稳产和调整挖潜面临越来越大的挑战。因此，准确描述和预测油层中的剩余油，特别是富集部位的分布状态是特高含水期改善油田开发效果的基础。开发实践表明，原来已经强采强注的强水淹或特强水淹而停采的油层或区域，由于地下动态条件的变化，经过一段时间后，剩余油在该油层或区域可以再次聚集，重新开采而取得显著效果。这种情况可以称之为油气重新运移聚集或再聚集或再富集，但对于剩余油再富集需要多长时间即再富集成藏周期尚缺乏简便实用的预测方法，为此我们针对水驱油藏特高含水后期剩余油分布更加分散、油田稳产和调整挖潜面临的难点，发明了一种新的水驱油藏剩余油再富集成藏周期预测方法，解决了以上技术问题。

发明内容

本发明的目的是提供一种水驱油藏剩余油再富集成藏周期预测方法，其利用不同含水饱和度条件下的油水运移速度公式计算水驱油藏剩余油再富集成藏周期之目的。

本发明的目的可通过如下技术措施来实现：水驱油藏剩余油再富集成藏周期预测方法，该水驱油藏剩余油再富集成藏周期预测方法包括：步骤1，求取所研究油藏的含水饱和度值；步骤2，根据该含水饱和度值选择相应的水驱油藏油水运移速度公式；以及步骤3，利用所选择的该水驱油藏油水运移速度公式计算剩余油再富集成藏需要的最长时间即剩余油再富集成藏周期。

本发明的目的还可通过如下技术措施来实现：

在步骤3中，根据油藏的原始含油高度、含油条带宽度情况，计算出原油上浮及水下沉所需要的总时间，即可得到该油藏剩余油再富集成藏所需的最长时间亦即剩余油再富集成藏周期。

当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(1\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(2\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}---\left(3\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1-\frac{1}{\sin \mathrm{\α}})]---\left(4\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

S_w、S_o—含水饱和度、含油饱和度，%；

g—重力加速度，m/s²；

α—地层倾角，°。

当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_w}{S_o}-{\mathrm{\ρ}}_o)g+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(5\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_w}{S_o}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(6\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}---\left(3\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1-\frac{1}{\sin \mathrm{\α}})]---\left(4\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

S_w、S_o—含水饱和度、含油饱和度，%；

g—重力加速度，m/s²；

α—地层倾角，°。

当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+2({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(7\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1+\frac{1}{\sin \mathrm{\α}})]---\left(8\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(9\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(10\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

S_w、S_o—含水饱和度、含油饱和度，%；

g—重力加速度，m/s²。

当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+2({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(7\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1+\frac{1}{\sin \mathrm{\α}})]---\left(8\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(9\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(10\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

S_w、S_o—含水饱和度、含油饱和度，%；

g—重力加速度，m/s²；

c—高含水阶段油水相对渗透率的比值与含水饱和度的函数关系式中回归系数，其中a、b也是回归系数，常数。

当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}---\left(11\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(\frac{1}{\sin \mathrm{\α}}-1)]---\left(12\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g-({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(13\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}-({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(14\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

v_⊥、v_∠—竖直、地层方向剩余油再富集成藏速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

g—重力加速度，m/s²；

S_w、S_o—含水饱和度、含油饱和度，%；

c—高含水阶段油水相对渗透率的比值与含水饱和度的函数关系式中回归系数，其中a、b也是回归系数，常数。

本发明中的水驱油藏剩余油再富集成藏周期预测方法，其利用不同含水饱和度条件下的油水运移速度公式计算水驱油藏剩余油再富集成藏周期。该预测方法解决了现有预测方法不能考虑浮力、毛管力及预测结果不理想等问题，可预测不同含水阶段、不同构造位置剩余油再富集成藏周期，为油藏剩余油潜力评价提供依据，具有预测成功率高、简便实用的特点。该方法为水驱油藏开发后期的挖潜与提高采收率，尤其是为特高含水阶段低效开发油藏的战略调整、周期注采开发油藏开发、近废弃油藏的产能恢复及再度开发等方面指明了方向，其推广应用前景广阔，经济社会效益显著。

附图说明

图1为本发明的水驱油藏剩余油再富集成藏周期预测方法的一具体实施例的流程图。

具体实施方式

为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合所附图式，作详细说明如下。

如图1所示，图1为本发明的水驱油藏剩余油再富集成藏周期预测方法的结构图。

在步骤101，根据油藏的地质及动态状况，求取该时刻油藏的含水饱和度。流程进入到步骤102。

在步骤102，根据该时刻油藏的含水饱和度大小，选择相对应的油水运移速度公式。

在一实施例中，油水运移速度公式如下所示：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(1\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(2\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}---\left(3\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1-\frac{1}{\sin \mathrm{\α}})]---\left(4\right)$

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_w}{S_o}-{\mathrm{\ρ}}_o)g+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(5\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_w}{S_o}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(6\right)$

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+2({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(7\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1+\frac{1}{\sin \mathrm{\α}})]---\left(8\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(9\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(10\right)$

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}---\left(11\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(\frac{1}{\sin \mathrm{\α}}-1)]---\left(12\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g-({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(13\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}-({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(14\right)$

式（1）～（4）为本发明的一具体实施例中当时水驱油藏油水运移过程中的渗流速度公式；

式（5）、（6）、（3）、（4）为本发明的一具体实施例中当时水驱油藏油水运移过程中的渗流速度公式；

式（7）～（10）为本发明的一具体实施例中当时水驱油藏油水运移过程中的渗流速度公式；

式（7）～（10）为本发明的一具体实施例中当时水驱油藏油水运移过程中的渗流速度公式；

式（11）～（14）为本发明的一具体实施例中当时水驱油藏油水运移过程中的渗流速度公式。

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

v_⊥、v_∠—竖直、地层方向剩余油再富集成藏速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_ro∠—竖直、地层方向油相相对渗透率，无因次；

K_rw⊥、K_rw∠—竖直、地层方向水相相对渗透率，无因次；

K_ro、K_rw—油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm³；

g—重力加速度，m/s²；

S_w、S_o—含水饱和度、含油饱和度，%；

c—高含水阶段油水相对渗透率的比值与含水饱和度的函数关系式中回归系数，其中a、b也是回归系数，常数。

流程进入到步骤103。

在步骤103，根据油藏的原始含油高度、含油条带宽度情况，计算出原油上浮及水下沉所需要的总时间，即可得到该油藏剩余油再富集成藏所需的最长时间亦即剩余油再富集成藏周期。流程结束。

在应用本发明的一具体实施例1中：

油田1为一典型条带状窄屋脊断块油藏，油藏埋深1950m，含油面积0.92km²，油层有效厚度9.3m，孔隙度27.5%，空气渗透率673×10^-3μm²，地层原油粘度9.3mPa.s，地质储量119×10⁴t，地层倾角12.5°，含油条带宽度长3.6km、宽100m～250m。试验前开油井1口，日产油0.4t/d，综合含水97.7%，动液面800m，采出程度31.5%，进入特高含水、近技术废弃阶段。

2008年4月，先后在原始含油边界之外利用邻块8口报废的油水井上返注水，最高单井日注达到500m³以上，在累计注水量达15×10⁴m³后，2008年10月扶停DXX1-22井，2009年2月扶停DXX104井，自喷生产且含水分别下降到55%和88%，其中DXX1-22井累计停产7.5年，停产前综合含水96.6%；DXX1-22井累计停产8年，停产前综合含水97.6%。之后又先后在构造高部位补孔上返DXX1-34、DXX100X44、DXX9X78、DXX100X54，扶停DXX18-1，均获自喷生产且含水大幅下降。

（1）根据相渗曲线，可知油藏停产时含水率97.1%对应的含水饱和度为62.5%。

（2）由停产时刻含水饱和度62.5%，确认选择公式（7）～（10）。

（3）考虑到辛1沙一4在重新开放前基本处于近技术废弃状态，因此，可以认为水动力压力梯度为零，即油藏处于静置状态。据公式计算，竖直方向上剩余油再富集成藏速度为5.65×10^-5cm/s，地层方向上剩余油再富集成藏速度为7.71×10^-5cm/s。于是，以油藏油水边界处起算，竖直方向上剩余油再富集成藏周期约为4.42～3.04年，平均2.13年，地层方向剩余油再富集成藏周期约为4.11～10.28年，平均7.20年。这就是说，在油藏静置时间内，油层大部分区域已基本完成剩余油再富集成藏过程，已在构造顶部形成了新的剩余油再富集油藏。这也是近技术废弃油藏辛1沙一4重新获得效益开发的根本原因。

本发明中的水驱油藏剩余油再富集成藏周期预测方法，针对水驱油藏特高含水后期剩余油分布更加分散、油田稳产和调整挖潜面临的难点，基于流体势原理和油气运移成藏理论，运用达西定律提出了水驱油藏剩余油再富集成藏周期预测方法。本发明利用不同时刻含水饱和度条件下的油水运移速度公式预测水驱油藏剩余油再富集成藏周期。该预测方法解决了现有预测方法不能考虑浮力、毛管力、重力及预测结果不理想等问题，可预测不同含水阶段、不同构造位置剩余油再富集成藏周期，为油藏剩余油潜力评价提供依据，具有预测成功率高、简便实用的特点，为特高含水阶段低效开发油藏的战略调整、周期注采油藏开发、近废弃油藏的产能恢复及再度开发等方面指明了方向。

Claims (7)

1.水驱油藏剩余油再富集成藏周期预测方法，其特征在于，该水驱油藏剩余油再富集成藏周期预测方法包括：

步骤1，求取所研究油藏的含水饱和度值；

步骤2，根据该含水饱和度值选择相应的水驱油藏油水运移速度公式；以及

步骤3，利用所选择的该水驱油藏油水运移速度公式计算剩余油再富集成藏需要的最长时间即剩余油再富集成藏周期。

2.根据权利要求1所述的水驱油藏剩余油再富集成藏周期预测方法，其特征在于，在步骤3中，根据油藏的原始含油高度、含油条带宽度情况，计算出原油上浮及水下沉所需要的总时间，即可得到该油藏剩余油再富集成藏所需的最长时间亦即剩余油再富集成藏周期。

3.根据权利要求1所述的水驱油藏剩余油再富集成藏周期预测方法，其特征在于，当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(1\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(2\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}---\left(3\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1-\frac{1}{\sin \mathrm{\α}})]---\left(4\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

S_w、S_o—含水饱和度、含油饱和度，%；

g—重力加速度，m/s²；

α—地层倾角，°。

4.根据权利要求1所述的水驱油藏剩余油再富集成藏周期预测方法，其特征在于，当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_w}{S_o}-{\mathrm{\ρ}}_o)g+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(5\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_w}{S_o}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(6\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}---\left(3\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1-\frac{1}{\sin \mathrm{\α}})]---\left(4\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

S_w、S_o—含水饱和度、含油饱和度，%；

g—重力加速度，m/s²；

α—地层倾角，°。

5.根据权利要求1所述的水驱油藏剩余油再富集成藏周期预测方法，其特征在于，当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+2({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(7\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1+\frac{1}{\sin \mathrm{\α}})]---\left(8\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(9\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(10\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

S_w、S_o—含水饱和度、含油饱和度，%；

g—重力加速度，m/s²。

6.根据权利要求1所述的水驱油藏剩余油再富集成藏周期预测方法，其特征在于，当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+2({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(7\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(1+\frac{1}{\sin \mathrm{\α}})]---\left(8\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(9\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(10\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

S_w、S_o—含水饱和度、含油饱和度，%；

g—重力加速度，m/s²；

c—高含水阶段油水相对渗透率的比值与含水饱和度的函数关系式中回归系数，其中a、b也是回归系数，常数。

7.根据权利要求1所述的水驱油藏剩余油再富集成藏周期预测方法，其特征在于，当时，该水驱油藏油水运移速度公式为：

$v_{o\⊥}=\frac{k_{\⊥}{K}_{\mathrm{ro}\⊥}}{{\mathrm{\μ}}_o}\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}---\left(11\right)$

$v_{o\∠}=\frac{k_{\∠}{K}_{\mathrm{ro}\∠}}{{\mathrm{\μ}}_o}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g(\frac{1}{\sin \mathrm{\α}}-1)]---\left(12\right)$

$v_{w\⊥}=\frac{k_{\⊥}{K}_{\mathrm{rw}\⊥}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}\sin \mathrm{\α}+(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g-({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(13\right)$

$v_{w\∠}=\frac{k_{\∠}{K}_{\mathrm{rw}\∠}}{{\mathrm{\μ}}_w}[\frac{\mathrm{\Δp}}{\mathrm{\Δl}}+\frac{(\frac{{\mathrm{\ρ}}_w{S}_o}{S_w}-{\mathrm{\ρ}}_o)g}{\sin \mathrm{\α}}-({\mathrm{\ρ}}_w-{\mathrm{\ρ}}_o)g]---\left(14\right)$

其中，v_o⊥、v_o∠—竖直、地层方向原油渗流速度，cm/s；

v_w⊥、v_w∠—竖直、地层方向原油渗流速度，cm/s；

v_⊥、v_∠—竖直、地层方向剩余油再富集成藏速度，cm/s；

k_⊥、k_∠—竖直、地层方向油层空气渗透率，10^-3μm²；

K_ro⊥、K_rw⊥—垂直方向上的油相、水相相对渗透率，无因次；

K_ro∠、K_rw∠—地层方向上的油相、水相相对渗透率，无因次；

μ_o、μ_w—地层原油、水粘度，mPa.s；

—水动力压力梯度，MPa/m；

α—地层倾角，°；

ρ_w、ρ_o—地层原油、水密度，g/cm3；

g—重力加速度，m/s²；

S_w、S_o—含水饱和度、含油饱和度，%；

c—高含水阶段油水相对渗透率的比值与含水饱和度的函数关系式中回归系数，其中a、b也是回归系数，常数。