E-mail Alert
-
致密油藏具有渗透率低、渗流阻力大、存在启动压力梯度等特点，传统的开发方式难以实现高效开发[1-2]。水平井、长缝压裂技术的发展，促进了致密油藏的进一步开发。长缝压裂技术通过油、水井的大型压裂制造长缝，使裂缝在井排方向上延伸，进而形成一种类似于水平井开发的定向、定量造长缝的压裂完井技术。通过长缝压裂，可以形成较大面积的渗流通道，大幅提高致密油藏产能和开发效果[3-4]。
朱维耀等人[5]的研究结果表明致密油藏压裂裂缝中流体的流动存在高速非达西渗流。苏玉亮等人[6-7]的室内实验及现场研究表明，当裂缝内流体的渗流速度较大时，会呈现出高速非达西渗流，达西定律将无法准确地描述其渗流规律。Holditch(1976)等人[8]最早将裂缝中流体的高速非达西渗流的影响考虑到现场压裂设计中。Mcdaniel(1989)和Vincent(2000)等人[9-10]的研究结果表明，裂缝内的高速非达西渗流将严重的影响油藏最终的采收率及经济效益。Gil(2003)及Hernandez(2004)等人[11-12]在前人研究的基础上分析得出，高速非达西系数的阈值为0.1，即高速非达西渗流引起的压力降小于或等于总压降的10%，则高速非达西渗流的影响可以忽略不计且高速非达西系数随着采油速度和裂缝半长的增加而增加。目前国内外致密油藏长缝压裂渗流规律和产能模型研究中，没有综合考虑致密油藏的储层特征及裂缝中高速非达西渗流的特点，具有较大的局限性。因此，在致密油藏存在启动压力梯度的基础上考虑了裂缝中高速非达西渗流的影响，建立了长缝压裂直井基质-裂缝复合流动模型，并计算分析了压力动态和最佳裂缝半长变化规律。
-
致密油藏液体的流动规律不符合常规的达西定律，在渗流时除粘滞力外，还要克服吸附层的阻力。当外加压力梯度大于启动压力梯度时，液体才开始流动。故可以应用运动方程描述致密油藏基质中流体的流动为
$$v = - \frac{k}{\mu }\left( {\frac{{\partial p}}{{\partial x}} - G} \right)$$ (1) 式中，v为流体流动速度，m/s; k为基质渗透率，10−3 μm2；μ为原油黏度，mPa · s；
$\dfrac{{\partial p}}{{\partial x}}$ 为压力梯度，MPa/m；G为启动压力梯度，MPa/m。长缝压裂技术是致密油藏油井增加产量、提高经济效益的有效措施。致密油藏经过水力压裂后，裂缝渗透率大大提高。此时，流体首先从储层基质流向裂缝，再沿着裂缝流向井底。当裂缝的面积较小时，裂缝内的液体流速较大。对于高速流动的液体，其流动规律遵循的是Forchheimer二项式方程，如式(2)所示[13]。分析二项式可以看出，当渗流速度很小时，平方项可以忽略不计，下式就转化为达西渗流公式。由此可以看出式(2)中的第一项表示黏滞阻力引起的压力损失，第二项则表示由惯性力而引起的压力损失。当渗流速度较小时，第一项占优势；当渗流速度较大时，第二项占优势。
$$ - \frac{{{\rm{d}}p}}{{{\rm{d}}x}} = \frac{\mu }{k}\bar v{\rm{ + }}\rho \beta \bar v\left| v \right|$$ (2) 式中，
$\rho $ 为流体密度，kg/m3；$\beta $ 为惯性因子，m−1；$\bar v$ 为平均流体流动速度，m/s。 -
致密油藏基质-裂缝耦合流动模型的基本假设：无限、均质地层中心1口垂直裂缝井；有限导流能力双翼对称裂缝，裂缝半长xf，缝宽wf，裂缝垂向完全穿透储层，裂缝渗透率kf为常数；流体首先从基质流到裂缝，进而从裂缝流到井筒，模型分为基质系统和裂缝系统，为了简化计算将裂缝均分n段；基质中流体流动考虑启动压力梯度，裂缝中流体流动考虑非达西渗流系数的影响。
-
$$ {p_{\rm{D}}} = \frac{{kh\left[ {{p_i} - p\left( {x,y,t} \right)} \right]}}{{Q{B_{\rm{o}}}\mu }} $$ (3) $$ {p_{{\rm{fD}}}} = \frac{{{k_{\rm{f}}}h\left[ {{p_i} - {p_{\rm{f}}}\left( {x,t} \right)} \right]}}{{Q{B_{\rm{o}}}\mu }} $$ (4) $$ F = \frac{{\beta \rho {k_{\rm{f}}}}}{\mu }\frac{{Q{B_{\rm{o}}}}}{{{w_{\rm{f}}}h}} $$ (5) $${q_{{\rm{fD}}}} = \frac{{{q_{{\rm{f}}\left( {{\rm{x}},{\rm{t}}} \right)}}{x_{\rm{f}}}}}{{Q{B_{\rm{o}}}}}$$ (6) $${t_{\rm{D}}} = \frac{{kt}}{{\phi \mu {c_{\rm{t}}}x_{\rm{f}}^2}}$$ (7) $${x_{\rm{D}}} = \frac{x}{{{x_{\rm{f}}}}}$$ (8) $$ {c_{{\rm{fD}}}} = \frac{{{k_{\rm{f}}}{w_{\rm{f}}}}}{{k{x_{\rm{f}}}}} $$ (9) $$ {c_\text{η}} = \frac{{{k_{\rm{f}}}}}{{{\phi _{\rm{f}}}{c_{{\rm{tf}}}}}}\frac{{\phi {c_{\rm{t}}}}}{k} $$ (10) $${G_{\rm{D}}} = {c_{\rm{f}}}{x_{\rm{f}}}G$$ (11) 式中，pD为无因次地层压力，下标D表示无因次属性；k为基质渗透率，10−3 μm；h为地层厚度，m；pi为第i个网格的地层压力，MPa；p为地层压力，MPa; Q为井底流量，cm3/s；Bo为原油体积系数，m3/m3；pfD为无因次裂缝压力；pf为裂缝压力，MPa；
$F$ 为无因次裂缝非达西渗流系数；wf为裂缝宽度，m；qfD为无因次裂缝流量；qf为裂缝流量，cm3/s；xf为裂缝半长，m；tD为无因次时间；t为真实时间，s；$\phi $ 为地层孔隙度；ct为综合压缩系数，MPa−1；xD为无因次长度；cfD为无因次裂缝导流系数；cη为无因次裂缝扩散系数；GD为无因次启动压力梯度。 -
(1)基质系统。由式(1)得考虑启动压力梯度的基质系统流体流动的无因次渗流方程为
$$\left\{ {\begin{array}{l} {\dfrac{{{\partial ^2}{p_{\rm{D}}}}}{{\partial x_{\rm{D}}^2}} + \dfrac{{{\partial ^2}{p_{\rm{D}}}}}{{\partial y_{\rm{D}}^2}} - {G_{\rm{D}}}\dfrac{{\partial {p_{\rm{D}}}}}{{\partial {x_{\rm{D}}}}} - {G_{\rm{D}}}\dfrac{{\partial {p_{\rm{D}}}}}{{\partial {y_{\rm{D}}}}} = \dfrac{{\partial {p_{\rm{D}}}}}{{\partial {t_{\rm{D}}}}}} \\ {{{\left. {\dfrac{{\partial {p_{\rm{D}}}}}{{\partial {y_{\rm{D}}}}}} \right|}_{{y_{\rm{D}}} = 0}} = {q_{{\rm{fD}}}}\left( {{x_{{\rm{D}},}}{t_{\rm{D}}}} \right) - \dfrac{{kh}}{{\mu {c_{\rm{f}}}Q{B_{\rm{o}}}}}{G_{\rm{D}}}} \\ p\left( {{x_{\rm{D}}} \to \infty ,{y_{\rm{D}}} \to \infty ,{t_{\rm{D}}}} \right) = 0 \\ p\left( {{x_{\rm{D}}},{y_{\rm{D}}},t = 0} \right) = 0 \end{array}} \right.$$ (12) (2)裂缝系统。由式(2)得考虑裂缝中非达西渗流系数的裂缝系统流动的无因次渗流方程为
$$\left\{ \begin{array}{l} {\dfrac{{\partial {p_{\rm{D}}}}}{{\partial {x_{\rm{D}}}}}\dfrac{\partial }{{\partial {x_{\rm{D}}}}}\left( {\delta \dfrac{{\partial {p_{{\rm{fD}}}}}}{{\partial {x_{\rm{D}}}}}} \right) - \dfrac{{{q_{{\rm{fD}}}}\left( {{x_{\rm{D}}},{t_{\rm{D}}}} \right)}}{{{c_{{\rm{fD}}}}}} = \dfrac{1}{{{c_\text{η}}}}\dfrac{{\partial {p_{{\rm{fD}}}}}}{{\partial {t_{\rm{D}}}}}} \\ {\left. {\delta \dfrac{{\partial {p_{{\rm{fD}}}}}}{{\partial {x_{\rm{D}}}}}} \right|_{{x_{\rm{D}}} = 0}} = \dfrac{1}{{2{c_{{\rm{fD}}}}}} \\ {\left. {\dfrac{{\partial {p_{\rm{D}}}}}{{\partial {x_{\rm{D}}}}}} \right|_{{x_{\rm{D}}} = 1}} = 0 \\ {p\left( {{x_{\rm{D}}},{t_{\rm{D}}} = 0} \right) = 0} \\ \delta = \dfrac{1}{{1 + F \times {q_{n{\rm{D}}}}}} \\ \end{array} \right.$$ (13) 式中，cf为裂缝导流系数；qnD为无因次节点流量。
进一步得到裂缝系统中任意一个网格i的无因次渗流方程为
$$\left\{ \begin{array}{l} \dfrac{{{\partial ^2}{p_{{\rm{fD}}}}}}{{\partial x_{\rm{D}}^2}} - \dfrac{{{q_{{\rm{fD}}i}}\left( {{x_{\rm{D}}},{t_{\rm{D}}}} \right)}}{{{\delta _i}{c_{{\rm{fD}}}}}} = \dfrac{1}{{{c_\text{η}}{\delta _i}}}\dfrac{{\partial {p_{{\rm{fD}}}}}}{{\partial {t_{\rm{D}}}}}\\ {\left. {\dfrac{{\partial {p_{{\rm{fD}}}}}}{{\partial {x_{\rm{D}}}}}} \right|_{{x_{\rm{D}}}_{i - 1}}} = \dfrac{1}{{{c_{{\rm{fD}}}}}}{\left. {\left( {\dfrac{{{q_{n{\rm{D}}}}}}{\delta }} \right)} \right|_{{x_{{\rm{D}}i - 1}}}}\\ {\left. {\dfrac{{\partial {p_{{\rm{fD}}}}}}{{\partial {x_{\rm{D}}}}}} \right|_{{x_{\rm{D}}}_i}} = \dfrac{1}{{{c_{{\rm{fD}}}}}}{\left. {\left( {\dfrac{{{q_{n{\rm{D}}}}}}{\delta }} \right)} \right|_{{x_{{{\rm{D}}_i}}}}}\\ {p_{{\rm{fD}}}}\left( {{x_{\rm{D}}},{t_{\rm{D}}} = 0} \right) = 0 \end{array} \right.$$ (14) -
对上述基质及裂缝的渗流方程进行Laplace变换并应用点源函数求解[14]，可得拉式空间内基质系统、裂缝系统的解为
$${\bar p_{\rm{D}}}\left( {{x_{{\rm{D}}j}},u} \right) = \mathop \sum \limits_1^n \left[ {\bar q_{{\rm{fD}}i}^{\rm{*}}\left( u \right) - \frac{{kh}}{{\mu {c_{\rm{f}}}Q{B_{\rm{o}}}}}{G_{\rm{D}}}} \right]{a_{ij}}$$ (15) $${\bar p_{{\rm{fD}}i}}\left( {{x_{\rm{D}}},{u_i}} \right) = {b_i}\left( {{x_{\rm{D}}}} \right)\bar q_{n{\rm{D}}i - 1}^* + {c_i}\left( {{x_{\rm{D}}}} \right)\bar q_{n{\rm{D}}i}^{\rm{*}} + {d_i}\bar q_{{\rm{fD}}i}^{\rm{*}}$$ (16) 式(15)、(16)中，
$$\begin{split} {a_{ij}} =& \frac{1}{u} {\int\limits_{{x_{{\rm{D}}i{\rm{ - }}1}}}^{{x_{{\rm{D}}i}}} {{K_{\rm{o}}}} \left( {\left| {{x_{{\rm{D}}j}} - {x_{\rm{D}}}^\prime } \right|\sqrt u } \right){\rm{d}}{x_{\rm{D}}}^\prime + } \frac{1}{u} {\int\limits_{{x_{{\rm{D}}i{\rm{ - }}1}}}^{{x_{{\rm{D}}i}}} {{K_{\rm{o}}}} \left( {\left| {{x_{{\rm{D}}j}} + {x_{\rm{D}}}^\prime } \right|\sqrt u } \right){\rm{d}}{x_{\rm{D}}}^\prime } \end{split}$$ (17) $$\begin{split}{b_i}\left( {{x_{\rm{D}}}} \right) = \frac{1}{{{c_{{\rm{fD}}}}\sqrt {{u_i} \cdot {c_\text{η} }^{ - 1}} }}\left\{ {\frac{{2\cos {\rm{h}}\left[ {\left( {{x_{\rm{D}}} - {x_{{\rm{D}}i{\rm{ - 1}}}}} \right)\sqrt {{u_i} \cdot {c_\text{η} }^{ - 1}} } \right]}}{{{\rm{exp}}\left[ {2\left( {{x_{{\rm{D}}i}} - {x_{{\rm{D}}i - 1}}} \right)\sqrt {{u_i} \cdot {c_\text{η} }^{ - 1}} } \right] - 1}} + } {{\rm{exp}}\left[ {\left( {{x_{{\rm{D}}i - 1}} - {x_{\rm{D}}}} \right)\sqrt {{u_i} \cdot {c_\text{η} }^{ - 1}} } \right]} \right\} \end{split}$$ (18) $$\begin{split} {c_i}\left( {{x_{\rm{D}}}} \right) = \frac{{ - 1}}{{{c_{{\rm{fD}}}}\sqrt {{u_i} \cdot {c_\text{η} }^{ - 1}} }}\left\{ {\frac{{2\cos {\rm{h}}\left[ {\left( {{x_{{\rm{D}}i}} - {x_{\rm{D}}}} \right)\sqrt {{u_i} \cdot {c_\text{η} }^{ - 1}} } \right]}}{{{\rm{exp}}\left[ {2\left( {{x_{{\rm{D}}i}} - {x_{{\rm{D}}i - 1}}} \right)\sqrt {{u_i} \cdot {c_\text{η} }^{ - 1}} } \right] - 1}} + } {{\rm{exp}}\left( {{x_{\rm{D}}} - {x_{{\rm{D}}i}}} \right)\sqrt {{u_i} \cdot {c_\text{η} }^{ - 1}} } \right\} \end{split}$$ (19) $${d_i} = - \frac{{{c_\text{η}}}}{{{c_{{\rm{fD}}}}{u_i}}}$$ (20) 式中，Ko为零阶贝塞尔函数；u为拉普拉斯变量；下标i表示节点序号或者网格序号。
由式(15)~式(16)将裂缝系统相邻两个网格的解沿网格之间的界面进行耦合求解，将基质系统和裂缝系统沿裂缝界面进行耦合求解。由此得到复合流动模型的耦合解为
$$\left[ {\begin{array}{*{20}{c}} {\dfrac{{{c_1}}}{{{\tau _1}}} - \dfrac{{{b_1}}}{{{\tau _2}}}}& \cdots &0&0&{\dfrac{{{d_1}}}{{{\tau _{\rm{1}}}}}}& \cdots &0 \\ \vdots & \ddots & \vdots & \vdots & \vdots & \ddots & \vdots \\ 0& \cdots &{\dfrac{{{c_{n - 1}}}}{{{\tau _{n - 1}}}} - \dfrac{{{b_{n - 1}}}}{{{\tau _n}}}}&{ - \dfrac{{{c_{n - 1}}}}{{{\tau _n}}}}&0& \cdots &{ - \dfrac{{{d_{n - 1}}}}{{{\tau _n}}}} \\ {{c_1}}& \cdots &0&0&{{d_1} - \dfrac{{{a_{11}}}}{{{\delta _1}}}}& \cdots &{ - \dfrac{{{a_{n1}}}}{{{\delta _n}}}} \\ \vdots & \ddots & \vdots & \vdots & \vdots & \ddots & \vdots \\ 0& \cdots &{{b_n}}&{{c_n}}&{ - \dfrac{{{a_{{\rm{1}}n}}}}{{{\delta _1}}}}& \cdots &{\dfrac{{{a_{nn}}}}{{{\delta _n}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} \begin{gathered} \bar q_{n{\rm{D1}}}^* \\ \bar q_{n{\rm{D2}}}^* \\ \end{gathered} \\ {\bar q_{n{\rm{D}}3}^*} \\ \begin{gathered} \vdots \\ \bar q_{n{\rm{D}}n{\rm{ - }}1}^* \\ \end{gathered} \\ \begin{gathered} \bar q_{n{\rm{D}}n}^* \\ \begin{array}{*{20}{c}} {\bar q_{{\rm{fD}}1}^*} \\ {\bar q_{{\rm{fD}}2}^*} \\ \vdots \\ {\bar q_{{\rm{fD}}n}^*} \end{array} \\ \end{gathered} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} { - \dfrac{{{b_1}}}{{{\tau _1}}}\bar q_{n{\rm{D0}}}^{\rm{*}}} \\ 0 \\ 0 \\ \vdots \\ 0 \\ { - {b_1}\bar q_{n{\rm{D0}}}^{\rm{*}} - {G_{\rm{D}}}\left( {{a_{11}} + \cdots + {a_{n{\rm{1}}}}} \right)} \\ { - {G_{\rm{D}}}\left( {{a_{12}} + \cdots + {a_{{\rm{2}}n}}} \right)} \\ \vdots \\ { - {G_{\rm{D}}}\left( {{a_{{\rm{1}}n}} + \cdots + {a_{nn}}} \right)} \end{array}} \right]$$ (21) 式中，
$${\delta _i} = \delta ({x_{\rm{D}}},{t_{\rm{D}}})$$ (22) $$ {x_{{\rm{D}}i - 1}} \leqslant {x_{\rm{D}}} \leqslant {x_{{\rm{D}}i}} $$ (23) $${\tau _{{\rm{D}}i}} = \int\limits_0^{{t_{\rm{D}}}} {\delta _i}\left( u \right){\rm{d}}u$$ (24) 通过编程求解以上线性方程组，得出裂缝内的流量分布及每1节点的流量。代入裂缝的第1个网格方程求解得到拉普拉斯域的井底压力为
$$ {\bar p_{\rm{D}}} = {b_1}\bar q_{n{\rm{D0}}}^{\rm{*}} + {c_1}\bar q_{n{\rm{D}}1}^{\rm{*}} + {d_1}\bar q_{{\rm{fD1}}}^{\rm{*}} $$ (25) Van Everdingen等人[15]给出了拉式空间内定产条件下井底流动压力和定压条件下井底产量之间的关系为
$$ \overline q{{_{\rm{D}}}} = \frac{1}{{{u^2}\overline \varphi {{_{\rm{D}}}} }} $$ (26) 式中，
$\overline q{{_{\rm{D}}}}$ 为定压生产时拉式空间无因次井底产量；$\overline \varphi{{ _{\rm{D}}}}$ 为定产量生产时拉式空间无因次井底压力。对式(25)、(26)进行stehfest数值反演[16]，可以得到实时域空间的无因次井底压力及无因次井底产量。
-
为验证解析模型的可靠性，利用COMSOL有限元模拟器，建立一个致密油藏长缝压裂数值模型，设定模型边界为10 000 m，通过偏微分方程系数进行非达西渗流系数定义，数值模型中其他致密油储层的相关参数与解析模型所用参数保持一致，具体参数见表1。数值模拟结果显示，当生产时间为6 000 d时，渗流所波及的半径远小于模型边界，这时可以认为模型是无限大的。结果对比如图1所示，研究建立的解析模型与数值模型拟合度较好。同时相对于数值模型，解析模型具有计算时间短、能够直接应用于现场等优点。
表 1 致密油藏储层参数
Table 1. Reservoir parameters of tight oil reservoir
地层参数 数值 地层参数 数值 地层孔隙度$\phi $ 0.1 综合压缩系数
${c_{\rm{t}}}$/MPa−10.003 裂缝孔隙度${\phi _{\rm{f}}}$ 0.001 原油体积系数
${B_{\rm{o}}}$/(m3 · m−3)1.1 地层渗透率
$k$/10−3 μm20.1 裂缝长度${x_{\rm{f}}}$/m 200 裂缝渗透率
${k_{\rm{f}}}$/10−3 μm21 000 裂缝宽度${w_{\rm{f}}}$/m 0.005 地层厚度h/m 15 井底产油量Q/(cm3∙s−1) 200 原油密度$\rho $/(g · cm−3) 0.85 原油黏度$\mu $/(mPa∙s) 2.5 油藏初始压力${p_i}$/MPa 30 井底压力$p$/MPa 15 
图 1 解析模型与数值模型
Figure 1. Analytical model and numerical model
-
根据储层参数(表1)，计算定产条件不同无因次启动压力梯度下无因次井底压力及压力导数随无因次时间变化曲线如图2所示。无因次井底压力及压力导数随无因次启动压力梯度增大而增大，但启动压力梯度对生产动态的后期影响较大，这是因为生产后期，压降传播的范围增大，随启动压力梯度增加，地层中消耗的能量增大，因而生产压差较大。
图 2 不同启动压力下压力及压力导数曲线
Figure 2. Pressure and pressure derivative at different threshold pressures
-
根据储层参数(表1)，计算定产条件不同无因次非达西渗流系数F下无因次井底压力及压力导数随无因次时间变化曲线如图3所示。由图可以看出，无因次井底压力及压力导数随着非达西渗流系数的增大而增大，且无因次非达西渗流系数对生产动态的前期影响较大。这是因为生产前期，裂缝中流体的渗流速度较大，由式(2)可知，Forchheimer渗流方程的第二项不可忽略，相当于增加了流体流动的渗流阻力，故此时非达西渗流对压力及压力导数的影响很大；而生产后期，生产压差扩大，裂缝中流体的渗流速度较小，Forchheimer渗流方程第二项可以忽略，因此非达西渗流对压力曲线的影响较小。
图 3 不同非达西渗流系数下压力及压力导数曲线
Figure 3. Pressure and pressure derivative at different non-Darcy flow coefficient
-
定压生产条件下，计算生产时间为10年时不同非达西渗流系数下累积产油量随裂缝半长的变化曲线如图4所示。由图可以看出，累积产油量随非达西渗流系数F的增大而减小；累积产油量随裂缝半长的增大而增大，但增大幅度越来越小，即存在最佳裂缝半长。
图 4 不同非达西渗流系数下累积产量随裂缝 半长的变化曲线
Figure 4. Variation of cumulative production with half fracture length at different non-Darcy flow coefficient
由图4得出致密油藏长缝压裂的最佳裂缝半长随裂缝中非达西渗流系数的变化曲线如图5所示。由图可以看出，致密油藏水力压裂直井的最佳裂缝半长随裂缝中非达西渗流系数的增大而减小。这是因为高速非达西渗流在裂缝中产生的附加压降，增加了流体的渗流阻力，从而需要的裂缝半长越短。
图 5 最佳裂缝半长随非达西渗流系数的变化曲线
Figure 5. Variation of optimal half fracture length with non-Darcy flow coefficient
-
(1)考虑致密油藏中启动压力梯度及裂缝中非达西渗流的影响，建立了致密油藏基质-裂缝复合流动半解析模型，通过拉普拉斯变换，应用点源函数和Stehest数值反演计算了真实空间域内的井底压力及井底产量。
(2)通过分析启动压力梯度、非达西渗流系数对生产动态的影响发现，无因次压力及压力导数随着启动压力梯度、非达西渗流系数的增大而增大，且启动压力梯度对生产后期的压力曲线影响较大，而非达西渗流系数对生产前期的压力曲线影响较大。
(3)对不同非达西渗流系数下致密油藏长缝压裂的裂缝半长进行优化分析，结果表明裂缝中非达西渗流越严重，累积井底产量越低，最佳裂缝半长越短。研究成果为致密油藏长缝压裂方案设计、试井分析提供理论指导。
Matrix-fracture coupling flow model of long-crack fracturing vertical well in tight oil reservoirs
-
摘要: 致密油藏长缝压裂压力动态和产量变化规律目前尚未明确，开发方案设计缺少理论模型指导。针对长缝压裂基质-裂缝复合流动问题，建立了考虑储层启动压力梯度、裂缝内高速非达西渗流影响的长缝压裂直井基质-裂缝复合流动模型，并运用Laplace变换、点源函数、Stehest数值反演等方法进行求解，明确了致密油藏长缝压裂井压力响应特征和影响因素。研究结果表明，无因次压力及压力导数随启动压力梯度的增加而增大，但启动压力梯度对无因次压力及压力导数的前期影响较小，后期影响较大；非达西渗流系数对无因次压力及压力导数曲线的前期影响较大，随着非达西渗流系数的增大，无因次压力及压力导数增大；裂缝中非达西渗流系数越大，压裂井的产量越低，且压裂井的最佳裂缝半长越小。研究成果为致密油藏长缝压裂方案设计、试井分析等工作提供理论指导。Abstract: Long-crack fracturing pressure behavior and production change laws of tight oil reservoirs have not been clarified clearly, and the development scheme is designed without the guidance of theoretical model. To deal with the matrix-fracture combined flow of long-crack fracturing, this paper established a matrix-fracture combined flow model of long-crack fracturing vertical well which considers the influences of threshold reservoir pressure gradient and intra-fracture high speed non-Darcy flow. Then, it was solved by means of Laplace transformation, point source function and Stehest numerical inversion. Finally, the pressure respond behaviors of long-crack fracturing well in tight oil reservoirs and their influential factors were determined. It is indicated that dimensionless pressure and pressure derivative increase with the increase of threshold pressure gradient, but the influence of threshold pressure gradient on dimensionless pressure and pressure derivative is less in the early stage and greater in the late stage. Non-Darcy flow coefficient has greater influence on dimensionless pressure and pressure derivative in the early stage. And with the increase of non-Darcy flow coefficient, dimensionless pressure and pressure derivative increase. The higher the non-Darcy flow coefficient of fracture is, the lower the production rate of fracturing well is and the shorter the optimal half fracture length of fracturing well is. The research results provide the theoretical guidance for long-crack fracturing scheme design and well test analysis of tight oil reservoir.
-
Key words:
- tight oil reservoir /
- long-crack fracturing /
- non-Darcy flow /
- coupling flow model
-
图 2 不同启动压力下压力及压力导数曲线
Figure 2. Pressure and pressure derivative at different threshold pressures
图 3 不同非达西渗流系数下压力及压力导数曲线
Figure 3. Pressure and pressure derivative at different non-Darcy flow coefficient
图 4 不同非达西渗流系数下累积产量随裂缝 半长的变化曲线
Figure 4. Variation of cumulative production with half fracture length at different non-Darcy flow coefficient
图 5 最佳裂缝半长随非达西渗流系数的变化曲线
Figure 5. Variation of optimal half fracture length with non-Darcy flow coefficient
表 1 致密油藏储层参数
Table 1. Reservoir parameters of tight oil reservoir
地层参数 数值 地层参数 数值 地层孔隙度 $\phi $ 0.1 综合压缩系数 ${c_{\rm{t}}}$ /MPa−10.003 裂缝孔隙度 ${\phi _{\rm{f}}}$ 0.001 原油体积系数 ${B_{\rm{o}}}$ /(m3 · m−3)1.1 地层渗透率 $k$ /10−3 μm20.1 裂缝长度 ${x_{\rm{f}}}$ /m200 裂缝渗透率 ${k_{\rm{f}}}$ /10−3 μm21 000 裂缝宽度 ${w_{\rm{f}}}$ /m0.005 地层厚度h/m 15 井底产油量Q/(cm3∙s−1) 200 原油密度 $\rho $ /(g · cm−3)0.85 原油黏度 $\mu $ /(mPa∙s)2.5 油藏初始压力 ${p_i}$ /MPa30 井底压力 $p$ /MPa15 
下载: 导出CSV
-
[1] 姜瑞忠, 韩光伟, 汪洋, 等. 低渗气藏特殊渗流机理及稳定产能预测方法[J]. 石油钻采工艺, 2015, 37(4):67-71. JIANG Ruizhong, HAN Guangwei, WANG Yang, et al. Research on special seepage mechanism and prediction method for stable productivity of low permeability gas pools[J]. Oil Drilling & Production Technology, 2015, 37(4): 67-71. [2] 尹虎, 王新海, 刘洪, 等. 考虑启动压力梯度的页岩气藏数值模拟[J]. 天然气与石油, 2012, 30(4):43-45. doi: 10.3969/j.issn.1006-5539.2012.04.012 YIN Hu, WANG Xinhai, LIU Hong, et al. Numerical simulation of shale gas reservoir considering startup pressure gradient[J]. Natural Gas & Oil, 2012, 30(4): 43-45. doi: 10.3969/j.issn.1006-5539.2012.04.012 [3] 赵益忠, 程远方, 刘钰川, 等. 启动压力梯度对低渗透油藏微观渗流及开发动态的影响[J]. 油气地质与采收率, 2013, 20(1):67-69. doi: 10.3969/j.issn.1009-9603.2013.01.017 ZHAO Yizhong, CHENG Yuanfang, LIU Yuchuan, et al. Study on influence of start-up pressure gradient to micro-seepage in low permeability reservoirs and development trends[J]. Petroleum Geology and Recovery Efficiency, 2013, 20(1): 67-69. doi: 10.3969/j.issn.1009-9603.2013.01.017 [4] 任宗孝, 吴晓东, 金岗, 等. 致密油藏体积压裂水平井半解析渗流模型[J]. 石油钻采工艺, 2016, 38(5):633-639. REN Zongxiao, WU Xiaodong, JIN Gang, et al. Semi-analytical seepage model for horizontal wells with volume fracturing in tight oil reservoirs[J]. Oil Drilling & Production Technology, 2016, 38(5): 633-639. [5] 朱维耀, 王明, 刘合, 等. 特低渗油层非达西渗流压裂水平井非定常渗流产能预测新方法[J]. 特种油气藏, 2010, 17(5):63-66. doi: 10.3969/j.issn.1006-6535.2010.05.017 ZHU Weiyao, WANG Ming, LIU He, et al. A new productivity model of non-Darcy flow for hydraulic fracturing horizontal wells in extra-low permeability reservoirs with unsteady porous flow[J]. Special Oil & Gas Reservoirs, 2010, 17(5): 63-66. doi: 10.3969/j.issn.1006-6535.2010.05.017 [6] 苏玉亮, 盛广龙, 王文东, 等. 页岩气藏多重介质耦合流动模型[J]. 天然气工业, 2016, 36(2):52-59. doi: 10.3787/j.issn.1000-0976.2016.02.007 SU Yuliang, SHENG Guanglong, WANG Wendong, et al. A multi-media coupling flow model for shale gas reservoirs[J]. Natural Gas Industry, 2016, 36(2): 52-59. doi: 10.3787/j.issn.1000-0976.2016.02.007 [7] YAO Y, LI G, QIN P. Seepage features of high-velocity non-Darcy flow in highly productive reservoirs[J]. Journal of Natural Gas Science and Engineering, 2015, 27: 1732-1738. doi: 10.1016/j.jngse.2015.10.039 [8] HOLDITCH S A, MORSE R A. The effects of non-Darcy flow on the behavior of hydraulically fractured gas wells[J]. Journal of Petroleum Technology, 1976, 28(10): 1169-1179. doi: 10.2118/5586-PA [9] MCDANIEL B W. Design of fracturing treatments for gas wells considering gel filter cake residue and non-Darcy flow effects[R]. SPE 19095, 1989. [10] VINCENT M C, PEARSON C M, KULLMAN J. Non-Darcy and multiphase flow in propped fractures: case studies illustrate the dramatic effect on well productivity[R]. SPE/CERI Gas Technology Symposium, 2000: 71-84. [11] GIL J A, OZKAN E, RAGHAVAN R. Fractured-well-test design and analysis in the presence of non-Darcy flow[J]. SPE Reservoir Evaluation & Engineering, 2003, 6(3): 185-196. [12] LOPEZ-HERNANDEZ H D, VALKO P P, PHAM T T. Optimum fracture treatment design minimizes the impact of non-Darcy flow effects[R]. SPE 90195, 2004: 1725-1740. [13] GEERTSMA J. Estimating the coefficient of inertial resistance in fluid flow through porous media[J]. Society of Petroleum Engineers Journal, 1974, 14(14): 445-450. [14] VAN KRUYSDIJK C P J W. Semianalytical modeling of pressure transients in fractured reservoirs[R]. SPE 18169, 1988. [15] VAN EVERDINGEN A F, HURST W. The application of the Laplace transformation to flow problems in reservoirs[J]. Journal of Petroleum Technology, 1949, 1(12): 305-324. doi: 10.2118/949305-G [16] STEHFEST H. Numerical inversion of Laplace transforms[J]. Communications of the ACM, 1970, 13(1): 47-49. doi: 10.1145/361953.361969 -
点击查看大图
图(5) / 表 (1)
计量
- 文章访问数: 3
- HTML全文浏览量: 117
- PDF下载量: 56
- 被引次数: 0
