CN105160146A - 水驱特征关系图版生成方法及装置 - Google Patents

Abstract

本发明提供了一种水驱特征关系图版生成方法及装置，所述的水驱特征关系图版生成方法包括：获取油粘度、水粘度、油藏产出液含水率、残余油饱和度、束缚水含水饱和度及原始状态下含油饱和度；根据所述油粘度、水粘度及残余油饱和度建立油藏产出液含水率的关系式；根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度与平均含水饱和度的表达式；对所述油藏产出液含水率进行微分运算，得到微分方程；求解得到水驱油藏含水率与采出程度关系的通解公式；根据所述通解公式中计算油藏产出液含水率与采出程度的关系值，并根据所述关系值绘制水驱特征关系图版。

Description

水驱特征关系图版生成方法及装置

技术领域

本发明涉及油田开发的油藏开采工程方法，特别是关于一种水驱特征关系图版生成方法及装置。

背景技术

水驱曲线是水驱油藏的特征曲线，利用水驱曲线不但可以预测水驱油藏的开发动态，而且还可预测水驱油藏的可采储量及采收率，在国内外已广泛应用于水驱油田的水驱开发效果评价。现有技术中，计算含水率与采出程度关系的方法有很多种，举例说明如下：

1)WelgeHJ.Asimplifiedmethodforcomputingoilrecoverybygasorwaterdrive[J].Trans.，AIME，1952，195:91-98。

该方法基于达西定律和物质守恒原理，推导了一维线性驱替***中平均含水饱和度与出口端含水饱和度关系式：

${\stackrel{\‾}{S}}_w=S_w+\frac{1-f_w}{f_w^{\′}}---\left(7\right)$

式中：为油藏中平均含水饱和度；S_w为油藏中油井处的含水饱和度；f_w为油藏产出液含水率；f_w'为含水率关于含水饱和度的导数。

2)童宪章.油井产状和油藏动态分析[M].石油工业出版社，1981，P37-41。

该方法基于乙型水驱曲线及国内外25个油藏数据，建立含水率与采出程度关系式：

$\lg \frac{f}{1-f}=7.5(R-R_m)+1.69---\left(8\right)$

关系式(2)称为童氏水驱特征关系式，式中：f为油藏含水率；R为油藏采出程度；R_m为油藏最终采收率。利用上式可绘制油田不同最终采收率下含水率与采出程度关系曲线图版，如图1所示。

3)陈元千.水驱曲线关系式的推导[J].石油学报，1985，6(2):69-78。

该方法基于油水两相的驱替理论和实验研究成果，近似地推导了含水率与采出程度关系式为：

$\lg \frac{f}{1-f}=A+BR---\left(9\right)$

其中： $A=log\frac{{\mathrm{\μ}}_o{B}_o{\mathrm{\γ}}_w}{{\mathrm{n\μ}}_w{B}_w{\mathrm{\γ}}_o}+\frac{m}{4.606}(3S_{wi}+S_{or}-1),$ B＝3mS_oi/4.606

式中：f为油藏含水率；R为油藏采出程度；μ_o,μ_w分别为油、水粘度，mPa·s；B_o,B_w分别为油、水体积系数；γ_o,γ_w分别为油、水的比重，kg/m³；m,n为与岩石及流体有关的常数；S_wi为束缚水含水饱和度；S_oi为原始状态下含油饱和度；S_or为残余油饱和度。

该方法指出当油田含水率达到一定值之后，水油比与采出程度在半对数轴***中成一直线关系。

现有的含水率与采出程度关系曲线图版，主要存在以下两个问题，一是水驱特征曲线上没有体现无水采油期原油的采出程度或是含水率要达到一定值后才适用，这与水驱曲线是水驱油藏固有的特征曲线相矛盾；二是将童氏水驱特征关系式(2)中的系数7.5或1.96人为地修改为其他常数，以满足目标研究区域实际水驱规律，缺乏理论性和通用性；三是水驱关系式的理论推导过程中考虑了较多的经验值或是近似值，对于实际油田拟合度较差，不能合理的分析油田的水驱开发特征和效果。

发明内容

本发明提供一种水驱特征关系图版生成方法及装置，以将油田实际生产数据绘制到理论图版上进行对比进而分析水驱开发效果。

为了实现上述目的，本发明实施例提供一种水驱特征关系图版生成方法，所述的水驱特征关系图版生成方法包括：

获取油粘度μ_o、水粘度μ_w、油藏产出液含水率f_w、残余油饱和度S_or、束缚水含水饱和度S_wi及原始状态下含油饱和度S_oi；

根据所述油粘度μ_o、水粘度μ_w及残余油饱和度S_or建立油藏产出液含水率f_w的关系式：

$f_w=1-\frac{50{\mathrm{\μ}}_w}{{\mathrm{\μ}}_o}{(1-S_{or}-S_{we})}^3---\left(1\right);$

根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_we与平均含水饱和度的表达式：

$S_{we}={\stackrel{\‾}{S}}_w-\frac{1-f_w}{f_w^{\′}\left(S_{we}\right)}---\left(2\right);$

其中， ${\stackrel{\‾}{S}}_w=S_{wi}+R(1-S_{wi})---\left(3\right);$

对所述油藏产出液含水率f_w进行微分运算，得到微分方程：

$f_w^{\′}\left(R\right)=\frac{\∂f_w}{\∂S_w}\frac{\∂S_{we}}{\∂{\stackrel{\‾}{S}}_{we}}\frac{\∂{\stackrel{\‾}{S}}_w}{\∂R}---\left(4\right);$

根据式(1)至(4)求解得到水驱油藏含水率与采出程度关系的通解公式：

$f_w=1-e^C{\[S_{oi}(1-S_{wi}-{\mathrm{RS}}_{oi}-S_{or})\]}^{\frac{3}{S_{oi}^2}}---\left(5\right);$

根据所述通解公式中计算油藏产出液含水率与采出程度的关系值，并根据所述关系值绘制水驱特征关系图版；

其中，S_we为出口端含水饱和度，常数C的表达式为：

$C=-\frac{3}{S_{oi}^2}ln\[S_{oi}(1-S_{wi}-R_m{S}_{oi}-S_{or})\]-3.912---\left(6\right).$

在一实施例中，根据所述油粘度μ_o、水粘度μ_w、出口端含水饱和度S_we及残余油饱和度S_or建立油藏产出液含水率f_w的关系式，包括：当地下油水粘度介于在1～10之间时，根据所述油粘度μ_o、水粘度μ_w、出口端含水饱和度S_we及残余油饱和度S_or建立油藏产出液含水率f_w的关系式。

在一实施例中，根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_w与平均含水饱和度的表达式，包括：在非活塞式水驱油条件下，根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_w与平均含水饱和度的表达式。

在一实施例中，根据质量守恒原理计算平均含水饱和度

为了实现上述目的，本发明实施例提供一种水驱特征关系图版生成装置，所述的水驱特征关系图版生成装置包括：

参数获取单元，用于获取油粘度μ_o、水粘度μ_w、油藏产出液含水率f_w、残余油饱和度S_or、束缚水含水饱和度S_wi、油藏采出程度R及原始状态下含油饱和度S_oi；

含水率关系式建立单元，用于根据所述油粘度μ_o、水粘度μ_w及残余油饱和度S_or建立油藏产出液含水率f_w的关系式：

$f_w=1-\frac{50{\mathrm{\μ}}_w}{{\mathrm{\μ}}_o}{(1-S_{or}-S_{we})}^3---\left(1\right);$

饱和度计算单元，用于根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_we与平均含水饱和度的表达式：

$S_{we}={\stackrel{\‾}{S}}_w-\frac{1-f_w}{f_w^{\′}\left(S_{we}\right)}---\left(2\right);$

其中， ${\stackrel{\‾}{S}}_w=S_{wi}+R(1-S_{wi})---\left(3\right);$

微分运算单元，用于对所述油藏产出液含水率f_w进行微分运算，得到微分方程：

$f_w^{\′}\left(R\right)=\frac{\∂f_w}{\∂S_w}\frac{\∂S_{we}}{\∂{\stackrel{\‾}{S}}_{we}}\frac{\∂{\stackrel{\‾}{S}}_w}{\∂R}---\left(4\right);$

通解公式求解单元，用于根据式(1)至(4)求解得到水驱油藏含水率与采出程度关系的通解公式：

$f_w=1-e^C{\[S_{oi}(1-S_{wi}-{\mathrm{RS}}_{oi}-S_{or})\]}^{\frac{3}{S_{oi}^2}}---\left(5\right);$

图版绘制单元，用于根据所述通解公式中计算油藏产出液含水率与采出程度的关系值，并根据所述关系值绘制水驱特征关系图版；

其中，S_we为出口端含水饱和度，常数C的表达式为：

$C=-\frac{3}{S_{oi}^2}ln\[S_{oi}(1-S_{wi}-R_m{S}_{oi}-S_{or})\]-3.912---\left(6\right).$

在一实施例中，所述含水率关系式建立单元具体用于：当地下油水粘度介于在1～10之间时，根据所述油粘度μ_o、水粘度μ_w、出口端含水饱和度S_we及残余油饱和度S_or建立油藏产出液含水率f_w的关系式。

在一实施例中，所述饱和度计算单元具体用于：在非活塞式水驱油条件下，根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_w与平均含水饱和度的表达式。

在一实施例中，饱和度计算单元还用于根据质量守恒原理计算平均含水饱和度

本发明利用微分方程理论推导含水率及采出程度的解析式，减少了经验公式及近似取值带来的误差。根据实际油田相渗曲线数据，确定不同的初边界条件，从而利用微分方程通解确定油田的水驱特征图版。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为传统的童氏水驱特征曲线图版及实际油田数据对比示意图；

图2为本发明实施例的水驱特征关系图版生成方法流程图；

图3为本发明实施例的水驱特征曲线图版与实际油田数据对比示意图；

图4为本发明实施例的水驱特征关系图版生成装置的结构示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

水驱曲线是水驱油藏的特征曲线，利用水驱曲线不但可以预测水驱油藏的开发动态，而且还可预测水驱油藏的可采储量及采收率，在国内外已广泛应用于水驱油田的水驱开发效果评价。注水或者天然边底水开发油田过程中，水体不断将原油往井底推进，油藏中含油饱和度随着含水饱和度的增加而减少，当含水饱和度达到一定值之后，油井处含水饱和度大于束缚水饱和度，油井见水，含水率开始增大，并且含水率与含水饱和度成一定关系式。由质量守恒原理及油田相渗曲线数据，可以得到含水率与原油采出程度的理论关系及曲线图版，将油田实际生产数据绘制到理论图版进行匹配，从而分析油田的开发效果并指导未来的开发调整，因此只有建立合理的水驱特征曲线关系式及图版，才能更好的分析指导油田高效开发。

因此，本发明实施例提供一种水驱特征关系图版生成方法，如图2所示，该水驱特征关系图版生成方法包括：

S201：获取油粘度μ_o、水粘度μ_w、油藏产出液含水率f_w、残余油饱和度S_or、束缚水含水饱和度S_wi及原始状态下含油饱和度S_oi；

S202：根据所述油粘度μ_o、水粘度μ_w及残余油饱和度S_or建立油藏产出液含水率f_w的关系式：

$f_w=1-\frac{50{\mathrm{\μ}}_w}{{\mathrm{\μ}}_o}{(1-S_{or}-S_{we})}^3---\left(1\right);$

S203：根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_we与平均含水饱和度的表达式：

$S_{we}={\stackrel{\‾}{S}}_w-\frac{1-f_w}{f_w^{\′}\left(S_{we}\right)}---\left(2\right);$

其中， ${\stackrel{\‾}{S}}_w=S_{wi}+R(1-S_{wi})---\left(3\right);$

S204：对所述油藏产出液含水率f_w进行微分运算，得到微分方程：

$f_w^{\′}\left(R\right)=\frac{\∂f_w}{\∂S_w}\frac{\∂S_{we}}{\∂{\stackrel{\‾}{S}}_{we}}\frac{\∂{\stackrel{\‾}{S}}_w}{\∂R}---\left(4\right);$

S205：根据式(1)至(4)求解得到水驱油藏含水率与采出程度关系的通解公式：

$f_w=1-e^C{\[S_{oi}(1-S_{wi}-{\mathrm{RS}}_{oi}-S_{or})\]}^{\frac{3}{S_{oi}^2}}---\left(5\right);$

S206：根据所述通解公式中计算油藏产出液含水率与采出程度的关系值，并根据所述关系值绘制水驱特征关系图版；

其中，S_we为出口端含水饱和度，常数C的表达式为：

$C=-\frac{3}{S_{oi}^2}ln\[S_{oi}(1-S_{wi}-R_m{S}_{oi}-S_{or})\]-3.912---\left(6\right).$

上述方程中，给定原始含油饱和度、束缚水饱和度、残余油饱和度以及不同的最终原油采收率代入到通解公式(5)中，从而计算得到不同最终采收率下一系列含水率与采出程度关系值，绘制关系曲线图，即得新型水驱特征曲线图版，如图3所示。

步骤S202中，根据艾富罗斯的实验理论研究表明，建立油藏产出液含水率f_w的关系式(1)需要满足的条件为：地下油水粘度比在1～10之间。

步骤S202中，根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_w与平均含水饱和度的表达式，一般需要在非活塞式水驱油条件下进行。

在一实施例中，计算平均含水饱和度可以根据质量守恒原理(即物质平衡原理)。

下面结合具体的实例说明上述水驱特征关系图版生成方法：

根据艾富罗斯的实验理论研究表明，当地下油水粘度比在1～10之间时，油水两相流动的出口端含水率为：

$f_w=1-\frac{50{\mathrm{\μ}}_w}{{\mathrm{\μ}}_o}{(1-S_{or}-S_{we})}^3---\left(1\right)$

在非活塞式水驱油条件下，由Buckley-Leverett水驱替理论和Welge方程求得出口端含水饱和度和平均含水饱和度表达式为S_we出口端含水饱和度

$S_w={\stackrel{\‾}{S}}_w-\frac{1-f_w}{f_w^{\′}\left(S_w\right)}$ 现有的(2)

由物质平衡原理知油层中平均含水饱和度可表示为

${\stackrel{\‾}{S}}_w=S_{wi}+R(1-S_{wi})---\left(3\right)$

同时由微分方程得

$f_w^{\′}\left(R\right)=\frac{\∂f_w}{\∂S_w}\frac{\∂S_w}{\∂{\stackrel{\‾}{S}}_w}\frac{\∂{\stackrel{\‾}{S}}_w}{\∂R}---\left(4\right)$

联立以上(4)～(7)求解得到水驱油藏含水率与采出程度关系的通解公式为

$f_w=1-e^C{\[S_{oi}(1-S_{wi}-{\mathrm{RS}}_{oi}-S_{or})\]}^{\frac{3}{S_{oi}^2}}---\left(5\right)$

其中常数C与极限含水率及最终采收率有关，当极限含水取0.98时，C的表达式为

$C=-\frac{3}{S_{oi}^2}ln\[S_{oi}(1-S_{wi}-R_m{S}_{oi}-S_{or})\]-3.912---\left(6\right)$

上述方程(6)中，给定原始含油饱和度、束缚水饱和度、残余油饱和度以及不同的最终原油采收率代入到本专利得到的通解中，从而计算得到不同最终采收率下一系列含水率与采出程度关系值，绘制关系曲线图，即得新型水驱特征曲线图版，如图1所示。

本发明基于前任实验研究成果及Welge平均含水饱和度方程，利用微分方程理论，推导出含水率与采出程度的通解公式，结合实际油田相渗及生产数据，绘制符合不同油田的理论水驱特征曲线。本次创立的理论图版(图1)与实际油田数据集童氏曲线(如图2)进行对比，可以看出童氏曲线并不能解释油田后期开发数据点，不能合理的预测油田最终采出程度(预测采收率大于60％)，而本发明绘制的图版(图1)不仅理论性更强，而且能较好的解释油田水驱开发特征及预测最终水驱采收率(预测采收率为45％与油藏标定采收43％接近)。

从两图中可以看出，随着油田的开采，含水率与采出程度的关系向横轴(采出程度)偏向，油田水驱效率都是越来越好。然而由Buckley-Leverett水驱理论知，当油藏均质性足够好时，无水采油期采出程度将较高，而从童氏理论图版可以看出，理论曲线中并没有出现无水采油期的值，同时童氏曲线并不能解释说明油田后期开发数据点的变化；而本发明的图版(图1)结合油田的实际数据可以绘制考虑无水采油期的实际油田水驱特征曲线，能够较好地分析油田开发效果及开发指标预测。

本发明实施例提供一种水驱特征关系图版生成装置，如图4所示，该水驱特征关系图版生成装置包括：参数获取单元401，含水率关系式建立单元402，饱和度计算单元403，微分运算单元404，通解公式求解单元405及图版绘制单元406。

参数获取单元401用于获取油粘度μ_o、水粘度μ_w、油藏产出液含水率f_w、残余油饱和度S_or、束缚水含水饱和度S_wi、油藏采出程度R及原始状态下含油饱和度S_oi；

含水率关系式建立单元402用于根据所述油粘度μ_o、水粘度μ_w及残余油饱和度S_or建立油藏产出液含水率f_w的关系式：

$f_w=1-\frac{50{\mathrm{\μ}}_w}{{\mathrm{\μ}}_o}{(1-S_{or}-S_{we})}^3---\left(1\right);$

饱和度计算单元403用于根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_we与平均含水饱和度的表达式：

$S_{we}={\stackrel{\‾}{S}}_w-\frac{1-f_w}{f_w^{\′}\left(S_{we}\right)}---\left(2\right);$

其中， ${\stackrel{\‾}{S}}_w=S_{wi}+R(1-S_{wi})---\left(3\right);$

微分运算单元404用于对所述油藏产出液含水率f_w进行微分运算，得到微分方程：

$f_w^{\′}\left(R\right)=\frac{\∂f_w}{\∂S_w}\frac{\∂S_{we}}{\∂{\stackrel{\‾}{S}}_{we}}\frac{\∂{\stackrel{\‾}{S}}_w}{\∂R}---\left(4\right);$

通解公式求解单元405用于根据式(1)至(4)求解得到水驱油藏含水率与采出程度关系的通解公式：

$f_w=1-e^C{\[S_{oi}(1-S_{wi}-{\mathrm{RS}}_{oi}-S_{or})\]}^{\frac{3}{S_{oi}^2}}---\left(5\right);$

图版绘制单元406用于根据所述通解公式中计算油藏产出液含水率与采出程度的关系值，并根据所述关系值绘制水驱特征关系图版；

其中，S_we为出口端含水饱和度，常数C的表达式为：

$C=-\frac{3}{S_{oi}^2}ln\[S_{oi}(1-S_{wi}-R_m{S}_{oi}-S_{or})\]-3.912---\left(6\right).$

在一实施例中，含水率关系式建立单元402具体用于：当地下油水粘度介于在1～10之间时，根据所述油粘度μ_o、水粘度μ_w、出口端含水饱和度S_we及残余油饱和度S_or建立油藏产出液含水率f_w的关系式。

在一实施例中，所述饱和度计算单元403具体用于：在非活塞式水驱油条件下，根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_w与平均含水饱和度的表达式。

在一实施例中，饱和度计算单元403还用于根据质量守恒原理计算平均含水饱和度

本发明针对目前推导及应用含水率与采出程度关系式存在的问题，基于前人实验成果，即油藏出口端含水率与含水饱和度满足一定的关系式；同时结合质量守恒原理以及Welge平均含水饱和度方程，建立含水率关于采出程度的微分方程，利用微分方程求解公式得到采出程度与含水率的通解关系式，从而根据实际相渗数据绘制相应的水驱特征曲线图版。将油田实际生产数据绘制到理论图版上进行对比进而分析水驱开发效果。

本发明利用微分方程理论推导含水率及采出程度的解析式，减少了经验公式及近似取值带来的误差。根据实际油田相渗曲线数据，确定不同的初边界条件，从而利用微分方程通解确定油田的水驱特征图版。

本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种水驱特征关系图版生成方法，其特征在于，所述的水驱特征关系图版生成方法包括：

获取油粘度μ_o、水粘度μ_w、油藏产出液含水率f_w、残余油饱和度S_or、束缚水含水饱和度S_wi及原始状态下含油饱和度S_oi；

根据所述油粘度μ_o、水粘度μ_w及残余油饱和度S_or建立油藏产出液含水率f_w的关系式：

$f_w=1-\frac{50{\mathrm{\μ}}_w}{{\mathrm{\μ}}_o}{(1-S_{or}-S_{we})}^3---\left(1\right);$

根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_we与平均含水饱和度的表达式：

$S_{we}={\stackrel{\‾}{S}}_w-\frac{1-f_w}{f_w^{\′}\left(S_{we}\right)}---\left(2\right);$

其中， ${\stackrel{\‾}{S}}_w=S_{wi}+R(1-S_{wi})---\left(3\right);$

对所述油藏产出液含水率f_w进行微分运算，得到微分方程：

$f_w^{\′}\left(R\right)=\frac{\∂f_w}{\∂S_w}\frac{\∂S_{we}}{\∂{\stackrel{\‾}{S}}_{we}}\frac{\∂{\stackrel{\‾}{S}}_w}{\∂R}---\left(4\right);$

根据式(1)至(4)求解得到水驱油藏含水率与采出程度关系的通解公式：

$f_w=1-e^C{\[S_{oi}(1-S_{wi}-{\mathrm{RS}}_{oi}-S_{or})\]}^{\frac{3}{S_{oi}^2}}---\left(5\right);$

根据所述通解公式中计算油藏产出液含水率与采出程度的关系值，并根据所述关系值绘制水驱特征关系图版；

其中，S_we为出口端含水饱和度，常数C的表达式为：

$C=-\frac{3}{S_{oi}^2}ln\[S_{oi}(1-S_{wi}-R_m{S}_{oi}-S_{or})\]-3.912---\left(6\right).$

2.根据权利要求1所述的水驱特征关系图版生成方法，其特征在于，根据所述油粘度μ_o、水粘度μ_w、出口端含水饱和度S_we及残余油饱和度S_or建立油藏产出液含水率f_w的关系式，包括：当地下油水粘度介于在1～10之间时，根据所述油粘度μ_o、水粘度μ_w、出口端含水饱和度S_we及残余油饱和度S_or建立油藏产出液含水率f_w的关系式。

3.根据权利要求1所述的水驱特征关系图版生成方法，其特征在于，根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_w与平均含水饱和度的表达式，包括：在非活塞式水驱油条件下，根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_w与平均含水饱和度的表达式。

4.根据权利要求1所述的水驱特征关系图版生成方法，其特征在于，根据质量守恒原理计算平均含水饱和度

5.一种水驱特征关系图版生成装置，其特征在于，所述的水驱特征关系图版生成装置包括：

参数获取单元，用于获取油粘度μ_o、水粘度μ_w、油藏产出液含水率f_w、残余油饱和度S_or、束缚水含水饱和度S_wi、油藏采出程度R及原始状态下含油饱和度S_oi；

含水率关系式建立单元，用于根据所述油粘度μ_o、水粘度μ_w及残余油饱和度S_or建立油藏产出液含水率f_w的关系式：

$f_w=1-\frac{50{\mathrm{\μ}}_w}{{\mathrm{\μ}}_o}{(1-S_{or}-S_{we})}^3---\left(1\right);$

饱和度计算单元，用于根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_we与平均含水饱和度的表达式：

$S_{we}={\stackrel{\‾}{S}}_w-\frac{1-f_w}{f_w^{\′}\left(S_{we}\right)}---\left(2\right);$

其中， ${\stackrel{\‾}{S}}_w=S_{wi}+R(1-S_{wi})---\left(3\right);$

微分运算单元，用于对所述油藏产出液含水率f_w进行微分运算，得到微分方程：

$f_w^{\′}\left(R\right)=\frac{\∂f_w}{\∂S_w}\frac{\∂S_{we}}{\∂{\stackrel{\‾}{S}}_{we}}\frac{\∂{\stackrel{\‾}{S}}_w}{\∂R}---\left(4\right);$

通解公式求解单元，用于根据式(1)至(4)求解得到水驱油藏含水率与采出程度关系的通解公式：

$f_w=1-e^C{\[S_{oi}(1-S_{wi}-{\mathrm{RS}}_{oi}-S_{or})\]}^{\frac{3}{S_{oi}^2}}---\left(5\right);$

图版绘制单元，用于根据所述通解公式中计算油藏产出液含水率与采出程度的关系值，并根据所述关系值绘制水驱特征关系图版；

其中，S_we为出口端含水饱和度，常数C的表达式为：

$C=-\frac{3}{S_{oi}^2}\ln \[S_{oi}(1-S_{wi}-R_m{S}_{oi}-S_{or})\]-3.912---\left(6\right).$

6.根据权利要求5所述的水驱特征关系图版生成装置，其特征在于，所述含水率关系式建立单元具体用于：当地下油水粘度介于在1～10之间时，根据所述油粘度μ_o、水粘度μ_w、出口端含水饱和度S_we及残余油饱和度S_or建立油藏产出液含水率f_w的关系式。

7.根据权利要求5所述的水驱特征关系图版生成装置，其特征在于，所述饱和度计算单元具体用于：在非活塞式水驱油条件下，根据Buckley-Leverett水驱替理论及Welge方程求得出口端含水饱和度S_w与平均含水饱和度的表达式。

8.根据权利要求5所述的水驱特征关系图版生成装置，其特征在于，饱和度计算单元还用于根据质量守恒原理计算平均含水饱和度