西南石油大学学报(自然科学版)  2020, Vol. 42 Issue (1): 147-154
压差时量的提出及平面单向流变量的等价性    [PDF全文]
刘海龙    
中国石油大庆油田有限责任公司勘探开发研究院, 黑龙江 大庆 163712
摘要: 针对平面单向流主动变化量与被动变化量之间函数关系有待明确的问题,提出了压差时量的概念,开展了平面单向流变量之间函数关系研究,在研究中利用压差时量与累积注入量之间的函数关系,结合Buckley Leverett方程和Welge方程,推导了压差时量与出口端含水饱和度、含水率、累积采油量、采出程度之间的函数关系,获得了以下认识:平面单向流任意两个变量之间存在由参数方程所确定的函数关系,压差时量是参数方程的参数;平面单向流中的变量具有等价性。根据以上认识,以层间矛盾研究为例,利用压差时量概念和变量等价性研究了层间矛盾产生的机理。压差时量的概念和变量等价性的认识,明确了主动变量与被动变量之间的函数关系,建立了驱替过程中各变量之间的等价关系,该认识对相关理论问题和应用问题的研究起到一定促进作用。
关键词: 压差时量    平面单向流    变量    等价性    层间矛盾    机理    函数关系    采收率    
The Theory of Pressure-difference-time and the Equivalence of Variables in Planar One-dimension Flow
LIU Hailong    
Exploration and Development Research Institute, Daqing Oilfield Company Limited, PetroChina, Daqing, Heilongjiang 163712, China
Abstract: In view of the fact that the functional relationship between the active and passive variations of planar one-dimension flow needs to be clarified, the concept of pressure difference time is proposed. The author studies the functional relationship between the active and passive variations of planar one-dimension flow. In this study, the functional relationship between pressure difference time and cumulative water injection given in the author's literature is used in Buckley-leverett equation and Welge equation to derive the functional relationship between pressure difference time and water saturation, water ratio, cumulative oil production and the degree of reserve recovery. The following understandings have been obtained:there is a functional relationship between any two variables of planar one-dimension flow, which is determined by the parametric equation. Pressure Difference Time is the parameter of the parametric equation. Variables of planar one-dimension flow are equivalent. Based on the above understanding, taking the study of interlayer interference as an example, the mechanism of interlayer interference is studied by using the concept of pressure difference time and the equivalence of variables. Using the concept of pressure difference time and the knowledge of equivalence of variables, the author clarifies the functional relationship between active variable and passive variable, and the equivalence relationship between variables in displacement process. This understanding has a certain role in promoting the studies of theoretical issues and application issues.
Keywords: pressure-difference-time    planar one-dimension flow    variable    equivalence    interlayer interference    mechanism    functional relationship    recovery ratio    
引言

在水驱油驱替过程中,通过施加压力,经过一段时间作用,产生渗流,多孔介质含水饱和度上升、含油饱和度下降,同时,采出端采出一定量的油和水。目前,平面单向流理论比较成熟和完整,是相关室内实验、现场实践、项目研究的基础[1-5]。平面单向流所涉及的物理量都是油田开发中重要的研究对象。

在平面单向流的变量中,注采压差和时间是可控制的量,即使研究者控制流速,实验设备也是通过压力调整实现流速控制。注采压差和时间是主动变化的量,其他变量随着注采压差和时间的变化而变化,有必要从主动变化、被动变化的角度,根据平面单向流基础理论,研究主动变量与被动变量之间的函数关系。

笔者在文献[6-7]中给出了压力、时间与累积注入量、注入倍数之间的函数关系,压力在时间上的积分值参与计算,本文提出以此积分值作为一个新的量,命名为压差时量。

在平面单向流主动变量与被动变量之间关系的研究中,目前渗流理论研究较多的是被动变量之间的关系[8-14]。例如,Buckley-Leverett方程、Welge方程涉及累积注入量、平均含水饱和度、出口端含水饱和度、某位置含水饱和度等,这些变量都是被动变量,或者在变量关系研究中涉及部分主动变量[15],例如在设定流速前提下,利用Buckley-Leverett方程研究时间与出口端含水饱和度之间的关系,笔者在文献[6-7]中给出了全部主动变量注采压差、时间与被动变量累积注入量之间的关系。本文提出的压差时量体现了注采压差和时间的综合作用。

主动变量在给定驱替条件下如何引起被动变量变化,是驱替过程中的因果顺序问题,在这一问题的研究上未见与压差时量相类似的内容。本文在研究主动变量与被动变量之间关系的过程中,总结出平面单向流中的变量具有等价性,各变量都能够反映驱替程度,不同的是物理意义和数值范围。

上述问题的研究未见资料报道,这些问题的提出来自变量关系研究上的需求,而非实际问题解决上的需求。

虽然研究内容不是来自于实际问题需求,但是主动变量与被动变量之间关系的研究能够为实际问题的解决提供新思路。压差时量体现了注采压差和时间的综合作用,产量递减、含水上升、井距配合、层间矛盾、开发效果差异等问题的机理研究,都直接或间接涉及同时性、时序性和变量的等价性。目前,上述实际问题的许多研究是依据拟合模型、数值试算、经验公式[16-18],例如,产量递减模型、井距组合数值模拟方案试算、层间矛盾数值模拟方案试算等。从同时性、时序性和变量等价性出发,从渗流原理上对上述实际问题进行机理研究,相关资料不多见,本文内容可为这些实际问题的机理研究提供不同的思路,在这些实际问题的研究上起到促进作用。

本文的研究对象是与压差时量有函数关系的变量,由压差时量无法确定的变量如瞬时流速等不作为研究对象。

1 压差时量的提出

在平面单向流变压力驱替过程中,注采两端压差曲线与时间轴在坐标系中围出曲边梯形,该曲边梯形的面积称为压差时量。压差时量的计算公式为

$ \begin{align} {A_{\rm{r}}} = \int_{{t_0}}^{{t_{\rm{n}}}} {\left[ {{p_{\rm{w}}}\left( t \right) - {p_{\rm{o}}}\left( t \right)} \right]} {\rm{d}}t \end{align} $ (1)

式中:${A_{\rm{r}}}$—压差时量,Pa$\cdot$s;

$t_0$—起始时间,s;

$t_{\rm{n}}$—结束时间,s;

${p_{\rm{w}}}\left(t \right)$—随时间变化的注入端压力,Pa;

${p_{\rm{o}}}\left(t \right)$—随时间变化的采出端压力,Pa;

$t$—时间,s。

在平面单向流中,压力、时间是由研究者控制的量。当研究者控制流速时,实验设备仍是通过压力控制实现流速控制。

在驱替过程中,压力和时间是主动变化的量,其他变量随着压力和时间的变化而变化,由式(1),压差时量包括压力、时间两个方面,压差时量体现压力和时间的综合作用。根据文献[6-7]有公式

采出端见水前

$ \begin{align} {Q_{{\rm{iq}}}} = \dfrac{{ - {\mu _{\rm{o}}}L\phi + {{\left\{ {{{\left( {{\mu _{\rm{o}}}L\phi } \right)}^2} + 2{A_{{\rm{r1}}}}K\phi \left[ {{\mu _{\rm{w}}}\alpha - {\mu _{\rm{o}}}f{'}\left( {{S_{{\rm{wf}}}}} \right)} \right]} \right\}}^{0.5}}}}{{{\mu _{\rm{w}}}\alpha - {\mu _{\rm{o}}}f{'}\left( {{S_{{\rm{wf}}}}} \right)}}{\rm{}} \end{align} $ (2)

式中:${Q_{{\rm{iq}}}}$—见水前累积注入量,m$^3$

$\mu _{\rm{o}}$—油的黏度,Pa$\cdot$s;

$L$—注采井距,m;

$\phi$—孔隙度,无因次;

${A_{{\rm{r1}}}}$—以水驱油开始驱替为起始时间至见水前某一时间的压差时量,Pa$\cdot$s;

$K$—渗透率,D;

$\mu _{\rm{w}}$—水的黏度,Pa$\cdot$s;

$\alpha$—常数,无因次,详见文献[6-7];

$f{'}\left({{S_{{\rm{wf}}}}} \right)$—含水率函数在前缘含水饱和度处的导数,无因次;

$S_{{\rm{wf}}}$—前缘含水饱和度,无因次。

采出端见水后

$ \begin{align} {A_{{\rm{r2}}}} = \dfrac{{{\mu _{\rm{w}}}{L^2}\phi }}{K}\varphi \left( {\dfrac{1}{{f{'}\left( {{S_{{\rm{we}}}}} \right)}}} \right){\rm{}} \end{align} $ (3)

式中:${A_{{\rm{r2}}}}$—以见水为起始时间至见水后某一时间的压差时量,Pa$\cdot$s;

$\varphi$—与相对渗透率曲线和油水黏度有关的单调增加函数,详见文献[6-7];

$S_{{\rm{we}}}$—出口端含水饱和度,无因次;

$f{'}\left({{S_{{\rm{we}}}}} \right)$—含水率函数在出口端含水饱和度处的导数,无因次;

${\dfrac{1}{{f'\left({{S_{{\rm{we}}}}} \right)}}}$—累积注水倍数,可计算累积注入量。

式(2)和式(3)给出了压差时量与累积注入量之间的函数关系。

2 平面单向流变量的等价性

油层含水饱和度分布如图 1所示,油层截面积取1~m$^2$

图1 含水饱和度分布图 Fig. 1 The distribution of water saturation

变量等价性推导过程需要用到以下两个公式

Buckley-Leverett方程

$ \begin{align} x = \dfrac{{f{'}\left( {{S_{\rm{w}}}} \right)}}{\phi }{Q_{\rm{i}}}{\rm{}} \end{align} $ (4)

式中:$x$—油层位置,m;

$S_{\rm{w}}$$x$位置处的含水饱和度,无因次;

$f{'}\left({{S_{\rm{w}}}} \right)$—含水率函数的导数,无因次;

$Q_{\rm{i}}$—累积注入量,m$^3$

Welge方程

$ \begin{align} \overline {{S_{\rm{w}}}} - {S_{{\rm{we}}}} = \dfrac{{1 - f\left( {{S_{{\rm{we}}}}} \right)}}{{f'\left( {{S_{{\rm{we}}}}} \right)}}{\rm{}} \end{align} $ (5)

式中:

$\overline {{S_{\rm{w}}}}$—平均含水饱和度,无因次;

$f\left({{S_{{\rm{we}}}}} \right)$—出口端含水率,无因次。

$\overline {{S_{\rm{w}}}}$$S_{{\rm{we}}}$求导

$ \begin{align} {\overline {{S_{\rm{w}}}} ^{'}} = \dfrac{{f{''}\left( {{S_{{\rm{we}}}}} \right)\left[ {f\left( {{S_{{\rm{we}}}}} \right) - 1} \right]}}{{[f{'}{{\left( {{S_{{\rm{we}}}}} \right)]}^2}}} \end{align} $

式中:$f{''}\left({{S_{{\rm{we}}}}} \right)$—含水率函数在出口端含水饱和度处的二阶导数,无因次。

见水后,含水率与出口端含水饱和度关系曲线一般是凸弧,因此,$f{''}\left({{S_{{\rm{we}}}}} \right) < 0$,可知,${\overline {{S_{\rm{w}}}} ^{'}} > 0$,所以$\overline {{S_{\rm{w}}}} $是出口端含水饱和度$S_{\rm{we}}$的单调增加函数。

2.1 累积注入量、累积采液量与压差时量的函数关系

式(4)中$x$取出口端位置$L$$S_{\rm{w}}$取前缘含水饱和度$S_{\rm{wf}}$,见水时累积注入量为

$ \begin{align} {Q_{{\rm{ij}}}} = \dfrac{{L\phi }}{{f{'}\left( {{S_{{\rm{wf}}}}} \right)}}{\rm{}} \end{align} $ (6)

式中:${Q_{{\rm{ij}}}}$—见水时累积注入量,m$^3$

由式(6)、式(2),见水时压差时量为

$ \begin{align} {A_{{\rm{rj}}}} = \dfrac{{{L^2}\phi \left[ {{\mu _{\rm{w}}}\alpha + {\mu _{\rm{o}}}f{'}\left( {{S_{{\rm{wf}}}}} \right)} \right]}}{{2Kf{'}{{\left( {{S_{{\rm{wf}}}}} \right)}^2}}}{\rm{}} \end{align} $ (7)

式中:${A_{{\rm{rj}}}}$—见水时压差时量,Pa$\cdot$s。

整理式(3),有

$ \begin{align} \dfrac{1}{{f{'}\left( {{S_{{\rm{we}}}}} \right)}} = {\varphi ^{ - 1}} \left({\dfrac{{{A_{r2}}K}}{{{\mu _{\rm{w}}}{L^2}\phi }}}\right) {\rm{}} \end{align} $ (8)

式(4)中,$x$取出口端位置$L$$S_{\rm{w}}$取出口端含水饱和度$S_{\rm{we}}$,有

$ \begin{align} L = \dfrac{{f{'}\left( {{S_{{\rm{we}}}}} \right)}}{\phi }{Q_{{\rm{ih}}}}{\rm{}} \end{align} $ (9)

式中:$Q_{{\rm{ih}}}$—见水后累积注入量,m$^3$

由式(8)、式(9),见水后累积注入量为

$ \begin{align} {Q_{{\rm{ih}}}} = L\phi {\varphi ^{ - 1}}\left(\dfrac{{{A_{{\rm{r2}}}}K}}{{{\mu _{\rm{w}}}{L^2}\phi }}\right) \end{align} $ (10)

$\varphi$函数是单调增加函数[6-7]$\varphi$函数导数大于0,由反函数求导方法可知,${\varphi ^{ - 1}}$函数导数大于0,${\varphi ^{ - 1}}$函数也是单调增加函数。

$A_{\rm{r}}$是以水驱油开始驱替为起始时间至任意时间的压差时量,则$A_{\rm{r}}$$A_{\rm{r1}}$$A_{\rm{r2}}$之间存在以下关系

见水前,$A_{\rm{r1}}$=$A_{\rm{r}}$,见水后,$A_{\rm{r2}}$=$A_{\rm{r}}-A_{{\rm{r}}j}$

代入式(2)、式(10),由式(2)、式(7)、式(10),累积注入量为

$ \begin{align} F_{Q_{\rm{i}}}\left( {{A_{\rm{r}}}} \right) = \left\{ {\begin{array}{l} {\dfrac{{ - {\mu _{\rm{o}}}L\phi + {{\left\{ {{{\left( {{\mu _{\rm{o}}}L\phi } \right)}^2} + 2{A_{\rm{r}}}K\phi \left[ {{\mu _{\rm{w}}}\alpha - {\mu _{\rm{o}}}f'\left( {{S_{{\rm{wf}}}}} \right)} \right]} \right\}}^{0.5}}}}{{{\mu _{\rm{w}}}\alpha - {\mu _{\rm{o}}}f'\left( {{S_{{\rm{wf}}}}} \right)}}, ~~~~{A_{\rm{r}}} \leqslant {A_{{\rm{rj}}}}}\\[11pt] {L\phi {\varphi ^{ - 1}}\left( {\dfrac{{\left( {{A_{\rm{r}}} - {A_{{\rm{rj}}}}} \right)K}}{{{\mu _{\rm{w}}}{L^2}\phi }}} \right), ~~~~{A_{\rm{r}}} > {A_{{\rm{rj}}}}} \end{array}} \right. \end{align} $ (11)

式中:

$F_{Q_{\rm{i}}}\left({{A_{\rm{r}}}} \right)$—累积注入量函数。

${{\varphi ^{ - 1}}}$函数是单调增加函数,分析式(11)可知,$F_{Q_{\rm{i}}}\left({{A_{\rm{r}}}} \right)$是单调增加函数。

在不考虑岩石、流体压缩性情况下,累积注入量与累积采液量相等

$ \begin{align} F_{Q_{\rm{p}}}\left( {{A_{\rm{r}}}} \right) = \left\{ {\begin{array}{l} {\dfrac{{ - {\mu _{\rm{o}}}L\phi + {{\left\{ {{{\left( {{\mu _{\rm{o}}}L\phi } \right)}^2} + 2{A_{\rm{r}}}K\phi \left[ {{\mu _{\rm{w}}}\alpha - {\mu _{\rm{o}}}f{'}\left( {{S_{{\rm{wf}}}}} \right)} \right]} \right\}}^{0.5}}}}{{{\mu _{\rm{w}}}\alpha - {\mu _{\rm{o}}}f{'}\left( {{S_{{\rm{wf}}}}} \right)}}, ~~~~{A_{\rm{r}}} \leqslant {A_{{\rm{rj}}}}}\\[11pt] {L\phi {\varphi ^{ - 1}}\left(\dfrac{{\left( {{A_{\rm{r}}} - {A_{{\rm{rj}}}}} \right)K}}{{{\mu _{\rm{w}}}{L^2}\phi }}\right), ~~~~{A_{\rm{r}}} > {A_{{\rm{rj}}}}} \end{array}} \right. \end{align} $ (12)

式中:

$F_{Q_{\rm{p}}}\left({{A_{\rm{r}}}} \right)$—累积采液量函数,是单调增加函数。

2.2 出口端含水饱和度与压差时量的函数关系

由式(4)、式(11)得含水饱和度分布

$ \begin{align} x = \dfrac{{f{'}\left( {{S_{\rm{w}}}} \right)}}{\phi }F_{Q_{\rm{i}}}\left( {{A_{\rm{r}}}} \right) \end{align} $ (13)

式(13)中,$x$取出口端位置$L$,出口端含水饱和度为

$ \begin{align} F_{S_{{\rm{we}}}}\left( {{A_{\rm{r}}}} \right) = {f^{{'} - 1}}\left( {\dfrac{{L\phi }}{{F_{Q_{\rm{i}}}\left( {{A_{\rm{r}}}} \right)}}} \right) \end{align} $

式中:$F_{S_{{\rm{we}}}}\left({{A_{\rm{r}}}} \right)$—出口端含水饱和度函数。

见水后,含水率与出口端含水饱和度关系曲线一般是凸弧,$f{''}\left({{S_{{\rm{we}}}}} \right) < 0$,所以,$f'\left({{S_{\rm{w}}}} \right)$是单调减小函数,由反函数求导方法可知,${f^{{'} - 1}}$的导数小于0,${f^{{'} - 1}}$是单调减小函数,由复合函数求导法则可知,$F_{S_{{\rm{we}}}}\left({{A_{\rm{r}}}} \right)$导数大于0,$F_{S_{{\rm{we}}}}\left({{A_{\rm{r}}}} \right)$是单调增加函数。

考虑见水前情况

$ \begin{align} F_{S_{{\rm{we}}}}\left( {{A_{\rm{r}}}} \right) = \left\{ {\begin{array}{l} {{S_{{\rm{wc}}}}, ~~~~{A_{\rm{r}}} < {A_{{\rm{rj}}}}}\\ {{f^{{'} - 1}}\left( {\dfrac{{L\phi }}{{F_{Q_{\rm{i}}}\left( {{A_{\rm{r}}}} \right)}}} \right), ~~~~{A_{\rm{r}}} \geqslant {A_{{\rm{rj}}}}} \end{array}} \right. \end{align} $ (14)

式中:$S_{{\rm{wc}}}$—束缚水饱和度,无因次。

见水后,$F_{S_{{\rm{we}}}}\left({{A_{\rm{r}}}} \right)$是单调增加函数。

2.3 含水率与压差时量的函数关系

含水率

$ \begin{align} f\left( {{S_{{\rm{we}}}}} \right) = \dfrac{1}{{1 + \dfrac{{{\mu _{\rm{w}}}}}{{{\mu _{\rm{o}}}}}\dfrac{{{K_{{\rm{ro}}}}\left( {{S_{{\rm{we}}}}} \right)}}{{{K_{{\rm{rw}}}}\left( {{S_{{\rm{we}}}}} \right)}}}} \end{align} $ (15)

式中:

${{K_{{\rm{ro}}}}\left({{S_{{\rm{we}}}}} \right)}$—油相相对渗透率,无因次;

${{K_{{\rm{rw}}}}\left({{S_{{\rm{we}}}}} \right)}$—水相相对渗透率,无因次。

式(14)代入式(15)得

$ \begin{align} F_f\left( {{A_{\rm{r}}}} \right) = \dfrac{1}{{1 + \dfrac{{{\mu _{\rm{w}}}}}{{{\mu _{\rm{o}}}}}\dfrac{{{K_{{\rm{ro}}}} {F_{S_{{\rm{we}}}}\left( {{A_{\rm{r}}}} \right)} }}{{{K_{{\rm{rw}}}} {F_{S_{{\rm{we}}}}\left( {{A_{\rm{r}}}} \right)} }}}} \end{align} $

式中:$F_f\left({{A_{\rm{r}}}} \right)$—含水率函数。

${f\left({{S_{{\rm{we}}}}} \right)}$${F_{S_{{\rm{we}}}}\left({{A_{\rm{r}}}} \right)}$是单调增加函数,由复合函数求导法则可知,$F_f\left({{A_{\rm{r}}}} \right)$导数大于0,$F_f\left({{A_{\rm{r}}}} \right)$函数是单调增加函数。

考虑见水前情况

$ \begin{align} F_f\left( {{A_{\rm{r}}}} \right) = \left\{ {\begin{array}{l} {0, ~~~~{A_{\rm{r}}} < {A_{{\rm{rj}}}}}\\ {\dfrac{1}{{1 + \dfrac{{{\mu _{\rm{w}}}}}{{{\mu _{\rm{o}}}}}\dfrac{{{K_{{\rm{ro}}}} {F_{S_{{\rm{we}}}}\left( {{A_{\rm{r}}}} \right)}}}{{{K_{{\rm{rw}}}} {F_{S_{{\rm{we}}}}\left( {{A_{\rm{r}}}} \right)} }}}}, ~~~~{A_{\rm{r}}} \geqslant {A_{{\rm{rj}}}}} \end{array}} \right. \end{align} $ (16)

见水后$F_f\left({{A_{\rm{r}}}} \right)$函数是单调增加函数。

2.4 累积采油量、采出程度与压差时量的函数关系

式(14)代入式(5),见水后平均含水饱和度

$ \begin{align} \overline {{S_{\rm{w}}}\left({{A_{\rm{r}}}} \right)} = \dfrac{{1 - f\left({F_{S_{{\rm{we}}}}\left({{A_{\rm{r}}}} \right)} \right)}}{{f{'}\left({F_{S_{{\rm{we}}}}\left({{A_{\rm{r}}}} \right)} \right)}} + F_{S_{{\rm{we}}}}\left({{A_{\rm{r}}}} \right){\rm{}} \end{align} $ (17)

式中:$\overline {{S_{\rm{w}}}\left({{A_{\rm{r}}}} \right)}$—见水后平均含水饱和度函数。

由前述式(5)单调性分析,平均含水饱和度$\overline {{S_{\rm{w}}}} $是含水饱和度${S_{{\rm{we}}}}$的单调增加函数,由前述式(14)单调性分析,${S_{{\rm{we}}}}$${{A_{\rm{r}}}}$的单调增加函数,由复合函数求导法则可知,式(17)$\overline {{S_{\rm{w}}}\left({{A_{\rm{r}}}} \right)} $导数大于0,$\overline {{S_{\rm{w}}}\left({{A_{\rm{r}}}} \right)}$是单调增加函数。

累积采油量

$ \begin{align} F_{N_{\rm{p}}}\left( {{A_{\rm{r}}}} \right) = L\phi \left[ {\overline {{S_{\rm{w}}}\left( {{A_{\rm{r}}}} \right)} - {S_{{\rm{wc}}}}} \right] \end{align} $

式中:

$F_{N_{\rm{p}}}\left({{A_{\rm{r}}}} \right)$—累积采油量函数。

考虑见水前情况

$ \begin{align} F_{N_{\rm{p}}}\left( {{A_{\rm{r}}}} \right) = \left\{ {\begin{array}{l} {F_{Q_{\rm{i}}}\left( {{A_{\rm{r}}}} \right), ~~~~{A_{\rm{r}}} \leqslant {A_{{\rm{rj}}}}}\\[11pt] {L\phi \left[ {\overline {{S_{\rm{w}}}\left( {{A_{\rm{r}}}} \right)} - {S_{{\rm{wc}}}}} \right], ~~~~{A_{\rm{r}}} > {A_{{\rm{rj}}}}} \end{array}} \right. \end{align} $ (18)

${\overline {{S_{\rm{w}}}\left({{A_{\rm{r}}}} \right)} }$是单调增加函数,所以$F_{N_{\rm{p}}}\left({{A_{\rm{r}}}} \right)$是单调增加函数。

由采出程度定义

$ \begin{align} F_R\left( {{A_{\rm{r}}}} \right) = \dfrac{{F_{N_{\rm{p}}}\left( {{A_{\rm{r}}}} \right)}}{{L\phi \left( {1 - {S_{{\rm{wc}}}}} \right)}} \end{align} $ (19)

式中:

$F_R\left({{A_{\rm{r}}}} \right)$—采出程度函数。

${F_{N_{\rm{p}}}\left({{A_{\rm{r}}}} \right)}$是单调增加函数,所以$F_R\left({{A_{\rm{r}}}} \right)$是单调增加函数。

2.5 变量之间的等价关系

由式(11)、式(12)、式(14)、式(16)、式(18)和式(19),可得到以下关系

$ \begin{align} \left\{ {\begin{array}{l} {{A_{\rm{r}}} = {A_{\rm{r}}}}\\ {{Q_{\rm{i}}} = F_{Q_{\rm{i}}}\left( {{A_{\rm{r}}}} \right)}\\ {{Q_{\rm{p}}} = F_{Q_{\rm{p}}}\left( {{A_{\rm{r}}}} \right)}\\ {{S_{{\rm{we}}}} = F_{S_{{\rm{we}}}}\left( {{A_{\rm{r}}}} \right)}\\ {\begin{array}{l} {{f_{}} = F_f\left( {{A_{\rm{r}}}} \right)}\\ {{N_{\rm{p}}} = F_{N_{\rm{p}}}\left( {{A_{\rm{r}}}} \right)}\\ {R = F_R\left( {{A_{\rm{r}}}} \right)} \end{array}} \end{array}} \right. \end{align} $ (20)

式中:

$Q_{\rm{p}}$—累积产液量,m$^3$

$N_{\rm{p}}$—累积采油量,m$^3$

$f$—含水率,无因次;

$R$—采出程度,无因次。

由式(20)可知,平面单向流7项变量中的任意两项变量之间,都存在由参数方程所确定的函数关系,压差时量是参数方程的参数。

由式(11)、式(12)、式(14)、式(16)、式(18)和式(19)单调性,根据参数方程求导方法,见水前,${A_{\rm{r}}}$${Q_{\rm{i}}}$${Q_{\rm{p}}}$${N_{\rm{p}}}$$R$,5项变量中的任意两项变量之间由参数方程所确定的函数是单调函数,见水前${S_{\rm{we}}}$${f_{\rm{w}}}$是不变的常量;见水后,${A_{\rm{r}}}$${Q_{\rm{i}}}$${Q_{\rm{p}}}$${S_{\rm{we}}}$${f_{\rm{w}}}$${N_{\rm{p}}}$$R$,7项变量中的任意两项变量之间由参数方程所确定的函数是单调函数。

见水前,选取${A_{\rm{r}}}$${Q_{\rm{i}}}$${Q_{\rm{p}}}$${N_{\rm{p}}}$$R$中的任意两项变量$c$$q$

如果$c$是条件,$q$是结论,由式(20)及上述单调性分析可知,$q$$c$的单调函数,所以由$c$可以推出$q$,不存在多解。

由式(20)及上述单调性分析可知,$c$$q$的单调函数,所以由$q$可以推出$c$,不存在多解。

所以,$c$$q$的充分必要条件。

由等价的定义,$c$$q$等价。

因此,平面单向流见水前${A_{\rm{r}}}$${Q_{\rm{i}}}$${Q_{\rm{p}}}$${N_{\rm{p}}}$$R$之间互相等价,这5项变量具有等价性。

同理,平面单向流见水后${A_{\rm{r}}}$${Q_{\rm{i}}}$${Q_{\rm{p}}}$${S_{\rm{we}}}$${f_{\rm{w}}}$${N_{\rm{p}}}$$R$之间互相等价,这7项变量具有等价性。

由以上推导过程可知,$c$$q$是与压差时量有单调函数关系的变量,压差时量充当参数方程中的参数,因此,在平面单向流中,见水前,与压差时量有单调函数关系的变量在单调区间内与上述5项变量等价,见水后,与压差时量有单调函数关系的变量在单调区间内与上述7项变量等价,与压差时量无单调函数关系的变量没有此结论,例如瞬时流速。

等价性应用举例:见水后,平面单向流含水率和采出程度等价,含水率98%等价于采出程度某一数值,含水率98%不变,对应采出程度某一数值不变,在不考虑岩石、流体压缩性情况下,平面单向流采收率不变。

3 层间矛盾机理研究中的应用

非均质多层油田注水开发后,由于储层物性差异,各层含水率、采出程度、水淹状况等方面产生差异,通常把这种差异概括为层间矛盾。

目前,层间矛盾问题有很多研究[19-23],本文应用压差时量概念和变量等价性对层间矛盾机理进行研究。

图 2给出了大庆油田某区块某层系内小层的动用情况,由图可知,各小层采出程度、含水、平均含水饱和度存在较大差异。

图2 小层动用情况对比图 Fig. 2 Contrast chart of producing situation

由于是一套层系,所以各小层的开发历程相同,注采压差和开发时间差别不大,可以认为各小层的压差时量相同。各小层渗透率、孔隙度、相对渗透率等性质存在不同,所以,各小层的式(20)参数方程中常量和函数曲线形态不同,因此,各小层参数方程的参数(压差时量)虽然相同,但是某一变量在不同小层的函数值却不相同,层间矛盾出现。

根据以上分析,产生层间矛盾的机理是:多个油层同时开发,在压差时量差异不大的情况下,压差时量与关注变量之间的函数关系在不同油层有所不同,导致关注变量在不同油层函数值不同,层间矛盾出现。

分析图 2,第14层与第16层比较,第16层采出程度较低,但是含水率却较高,在不进行油层改造、层内调整的情况下,如果增加第16层采出程度,减小采出程度层间矛盾,根据变量的等价性,采出程度增加等价于含水率增加,使第16层含水率更高,导致含水率层间矛盾增加。第12层、第13层含水率与含水饱和度也存在类似情况。

由以上分析可知,在不进行油层改造、层内调整的情况下,小层之间采出程度、含水率、含水饱和度3项变量不能同时一致,某一变量一致,其他两项变量就会出现不一致,这是驱替过程中变量的等价性和储层物性差异造成的。

通过油层改造、层内调整能够实现多个变量层间矛盾同时减小的效果,例如,第16层可通过实施压裂、补孔、封堵等措施既提高采出程度又控制含水,原理是这些措施改变了储层物性和驱替条件,变量之间函数关系中的常量发生改变,使变量之间不再具有等价性,从而达到多个变量同时均衡的调整效果,实现稳油控水的目的。

4 结论

(1)在平面单向流变压力驱替过程中,注采两端压差曲线与时间轴在坐标系中围出曲边梯形,该曲边梯形的面积称为压差时量。

(2)压差时量、累积注入量、累积采液量、出口端含水饱和度、含水率、累积采油量、采出程度7项变量中的任意两项变量之间,都存在由参数方程所确定的函数关系,压差时量是参数方程的参数。

(3)见水前,压差时量、累积注入量、累积采液量、累积采油量、采出程度5项变量具有等价性;见水后,压差时量、累积注入量、累积采液量、出口端含水饱和度、含水率、累积采油量、采出程度7项变量具有等价性。

(4)不考虑岩石、流体压缩性情况下,平面单向流采收率不变。

(5)在平面单向流中,见水前,与压差时量有单调函数关系的变量在单调区间内与上述5项变量等价;见水后,与压差时量有单调函数关系的变量在单调区间内与上述7项变量等价。

(6)产生层间矛盾的机理:多个油层同时开发,在压差时量差异不大的情况下,压差时量与关注变量之间的函数关系在不同油层有所不同,导致关注变量在不同油层函数值不同,层间矛盾出现。

参考文献
[1]
喻高明. 非均质油藏注水开发指标计算方法的改进[J]. 石油勘探与开发, 1997, 24(2): 80-83.
YU Gaoming. A correction to the budgetary estimating method used in the development design for waterflooding reservoir[J]. Petroleum Exploration and Development, 1997, 24(2): 80-83. doi: 10.3321/j.issn:1000-0747.1997.02.020
[2]
宋兆杰, 李治平, 赖枫鹏, 等. 高含水期油田水驱特征曲线关系式的理论推导[J]. 石油勘探与开发, 2013, 40(2): 201-208.
SONG Zhaojie, LI Zhiping, LAI Fengpeng, et al. Derivation of water flooding characteristic curve for high water-cut oilfields[J]. Petroleum Exploration and Development, 2013, 40(2): 201-208. doi: 10.11698/PED.2013.02.09
[3]
吕爱民, 王立伟, 龙涛, 等. 基于Buckley-Leverett理论的典型缝洞油藏一维驱替开发指标计算方法[J]. 西安石油大学学报(自然科学版), 2018, 33(1): 61-67, 74.
LÜ Aimin, WANG Liwei, LONG Tao, et al. Calculation method for one dimensional displacement development indexes of typical fracture-vuggy reservoir models based on buckley-leverett theory[J]. Journal of Xi'an Shiyou University (Natural Science Edition), 2018, 33(1): 61-67, 74. doi: 10.3969/j.issn.1673-064X.2018.01.010
[4]
宋考平, 吴玉树, 计秉玉. 水驱油藏剩余油饱和度分布预测的φ函数法[J]. 石油学报, 2006, 27(3): 91-95.
SONG Kaoping, WU Yushu, JI Bingyu. A (φ)-function method for estimating distribution of residual oil saturation in water drive reservoir[J]. Acta Petrolei Sinica, 2006, 27(3): 91-95. doi: 10.3321/j.issn:0253-2697.2006.03.020
[5]
胡伟, 杨胜来, 翟羽佳, 等. 油-水相对渗透率曲线优化校正新方法[J]. 石油学报, 2015, 36(7): 871-875.
HU Wei, YANG Shenglai, ZHAI Yujia, et al. A new optimization and correction method of oil-water phase relative permeability curve[J]. Acta Petrolei Sinica, 2015, 36(7): 871-875. doi: 10.7623/syxb201507011
[6]
刘海龙. 一维水驱油恒压驱替渗流过程推导[J]. 大庆石油学院学报, 2012, 36(3): 90-95.
LIU Hailong. The derivation of one-dimensional displacement of oil by water in constant pressure[J]. Journal of Daqing Petroleum Institute, 2012, 36(3): 90-95. doi: 10.3969/j.issn.2095-4107.2012.03.016
[7]
刘海龙, 由春梅, 姜雪岩, 等. 一维水驱油变压力驱替过程推导[J]. 石油地质与工程, 2014, 28(5): 80-82, 86.
LIU Hailong, YOU Chunmei, JIANG Xueyan, et al. The derivation of one-dimensional water/oil displacement in variation pressure[J]. Petroleum Geology and Engineering, 2014, 28(5): 80-82, 86. doi: 10.3969/j.issn.1673-8217.2014.05.025
[8]
金蓉蓉. 新型含水率与采出程度关系理论曲线的推导[J]. 大庆石油地质与开发, 2015, 34(3): 72-75.
JIN Rongrong. Derivation of the new theoretical relationship curve between watercut and recovery factor[J]. Petroleum Geology and Oilfield Development in Daqing, 2015, 34(3): 72-75. doi: 10.3969/J.ISSN.1000-3754.2015.03.013
[9]
杜晓康, 李治平, 田丰, 等. 水驱油田高含水期体积波及系数预测方法[J]. 大庆石油地质与开发, 2018, 37(1): 83-87.
DU Xiaokang, LI Zhiping, TIAN Feng, et al. Predicting method of the volumetric swept coefficient for high-watercut waterflooded oilfields[J]. Petroleum Geology and Oilfield Development in Daqing, 2018, 37(1): 83-87. doi: 10.19597/J.ISSN.1000-3754.201703011
[10]
邓森, 王怒涛, 孟令强, 等. 高含水期两种新型水驱特征曲线的建立与应用[J]. 大庆石油地质与开发, 2017, 36(4): 58-63.
DENG Sen, WANG Nutao, MENG Lingqiang, et al. Establishment and application of the new two-type water-flooding characteristic curves at high watercut stage[J]. Petroleum Geology and Oilfield Development in Daqing, 2017, 36(4): 58-63. doi: 10.19597/J.ISSN.1000-3754.201612012
[11]
张金庆. 水驱曲线的进一步理论探讨及童氏图版的改进[J]. 中国海上油气, 2019, 31(1): 86-93.
ZHANG Jinqing. A further theoretical discussion on water flooding curve and improvement of Tong's chart[J]. China Offshore Oil and Gas, 2019, 31(1): 86-93. doi: 10.11935/j.issn.1673-1506.2019.01.010
[12]
李跃刚, 肖峰, 徐文, 等. 基于气水相对渗透率曲线的产水气井开采效果评价——以苏里格气田致密砂岩气藏为例[J]. 天然气工业, 2015, 35(12): 27-34.
LI Yuegang, XIAO Feng, XU Wen, et al. Performance evaluation on water-producing gas wells based on gas-water relative permeability curves:A case study of tight sandstone gas reservoirs in the Sulige Gas Field, Ordos Basin[J]. Natural Gas Industry, 2015, 35(12): 27-34. doi: 10.3787/j.issn.1000-0976.2015.12.004
[13]
苏海波. 反映动态启动压力梯度的低渗透油藏渗流模型[J]. 西南石油大学学报(自然科学版), 2015, 37(6): 105-111.
SU Haibo. A flow model reflecting dynamic kick-off pressure gradient for low permeability reservoirs[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2015, 37(6): 105-111. doi: 10.11885/j.issn.1674-5086.2015.05.06.05
[14]
薛颖, 石立华, 席天德. 童氏水驱曲线的改进及应用[J]. 西南石油大学学报(自然科学版), 2015, 37(1): 141-145.
XUE Ying, SHI Lihua, XI Tiande. Tong's water drive curve improvement and its application[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2015, 37(1): 141-145. doi: 10.11885/j.issn.1674-5086.2012.12.20.02
[15]
高文君, 徐冰涛, 黄瑜, 等. 水驱油田含水率预测方法研究及拓展[J]. 石油与天然气地质, 2017, 38(5): 993-999.
GAO Wenjun, XU Bingtao, HUANG Yu, et al. Research on and development of prediction method of water cut in water flooding oilfield[J]. Oil & Gas Geology, 2017, 38(5): 993-999. doi: 10.11743/ogg20170518
[16]
陈元千, 唐玮. 广义递减模型的建立及应用[J]. 石油学报, 2016, 37(11): 1410-1413.
CHEN Yuanqian, TANG Wei. Establishment and application of generalized decline model[J]. Acta Petrolei Sinica, 2016, 37(11): 1410-1413. doi: 10.7623/syxb201611009
[17]
吴晓慧. 大庆长垣油田特高含水期水驱精细挖潜措施后产量变化规律[J]. 大庆石油地质与开发, 2018, 37(5): 71-75.
WU Xiaohui. Changed laws of the production after waterflooding finely-tapped-potential stimulations for Daqing placanticline oilfields at the stage of extra-high watercut[J]. Petroleum Geology and Oilfield Development in Daqing, 2018, 37(5): 71-75. doi: 10.19597/J.ISSN.10003754.201804038
[18]
苑志旺, 杨宝泉, 杨莉, 等. 深水浊积砂岩油田含水上升机理及优化注水技术——以西非尼日尔三角洲盆地AKPO油田为例[J]. 石油勘探与开发, 2018, 45(2): 287-296.
YUAN Zhiwang, YANG Baoquan, YANG Li, et al. Water-cut rising mechanism and optimized water injection technology for deepwater turbidite sandstone oilfield:A case study of AKPO Oilfield in Niger Delta Basin, West Africa[J]. Petroleum Exploration and Development, 2018, 45(2): 287-296. doi: 10.11698/PED.2018.02.11
[19]
张继成, 何晓茹, 周文胜, 等. 大段合采油井层间干扰主控因素研究[J]. 西南石油大学学报(自然科学版), 2015, 37(4): 101-106.
ZHANG Jicheng, HE Xiaoru, ZHOU Wensheng, et al. Main controlling factors of interlayer interference in big intervals commingled production oil wells[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2015, 37(4): 101-106. doi: 10.11885/j.issn.1674-5086.2015.03.12.02
[20]
周文胜, 李倩茹, 耿站立, 等. 合采井层间干扰现象数学模拟研究[J]. 西南石油大学学报(自然科学版), 2017, 39(6): 109-116.
ZHOU Wensheng, LI Qianru, GENG Zhanli, et al. Mathematical simulation study on interlayer interference in commingled production[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2017, 39(6): 109-116. doi: 10.11885/j.issn.1674-5086.2016.08.20.01
[21]
谢一婷, 陈朝晖, 柴小颖, 等. 利用采油速度干扰系数评价油藏层间干扰[J]. 西南石油大学学报(自然科学版), 2016, 38(6): 119-124.
XIE Yiting, CHEN Zhaohui, CHAI Xiaoying, et al. Interlayer interference evaluation by the recovery rate interference factor for the waterflooding oil reservoir in the medium-late period[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2016, 38(6): 119-124. doi: 10.11885/j.issn.16745086.2015.11.11.25
[22]
冯其红, 王波, 王相, 等. 高含水油藏细分注水层段组合优选方法研究[J]. 西南石油大学学报(自然科学版), 2016, 38(2): 103-108.
FENG Qihong, WANG Bo, WANG Xiang, et al. Study on layer combination optimization method of subdivision water injection in high water-cut reservoirs[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2016, 38(2): 103-108. doi: 10.11885/j.issn.1674-5086.2013.12.17.04
[23]
黄世军, 康博韬, 程林松, 等. 海上普通稠油油藏多层合采层间干扰定量表征与定向井产能预测[J]. 石油勘探与开发, 2015, 42(4): 488-495.
HUANG Shijun, KANG Botao, CHENG Linsong, et al. Quantitative characterization of interlayer interference and productivity prediction of directional wells in the multilayer commingled production of ordinary offshore heavy oil reservoirs[J]. Petroleum Exploration and Development, 2015, 42(4): 488-495. doi: 10.11698/PED.2015.04.10