E-mail Alert
下载:
-
缝洞型碳酸盐岩油气藏储量丰富，天然裂缝普遍存在于储层中，在压裂增产过程中天然裂缝对水力裂缝延伸方向和扩展形态具有重要影响，能否利用天然裂缝沟通溶洞储集体是储层改造成功与否的关键[1-4]。
许多专家学者针对储层天然裂缝对油气田开发的影响问题进行了理论和实验研究。在油气藏开采方面，樊建明等[5]采用高斯函数对地应力、天然裂缝、人工裂缝各自优势方向进行耦合，将利用自然能量开发井网作为研究对象建立数值模型，研究不同优势方向下油藏的开采方法。姚军等[6]研究了油藏中大裂缝分布、长度、方位对水驱油的影响。苏玉亮等[7]通过对天然裂缝导流能力、长度、密度、方位如何影响油藏开采展开分析，得出天然裂缝的密度、方位对采收率有较大影响。李博[8]利用嵌入式离散裂缝模型对缝洞型油藏进行模拟研究，模拟结果表明天然裂缝数量越多，储层连通性越好，充分连通水力裂缝和天然裂缝对增产具有重要意义。孙若凡[9]在水平井压裂产能模型中考虑了天然裂缝，发现天然裂缝所占比例相对于它的渗透率对产能影响更大。朱圣举等[10]结合现场工程实际，对影响天然裂缝开启压力的因素进行了研究，结果表明对注水压力的控制可以优化开发效果。刘熙远[11]建立了三维地质模型，在储层注水开发过程中考虑天然裂缝，认为沿着裂缝方向注采井连通性更强。吴忠宝等[12]基于等值渗流阻力法原理和现场施工资料，建立微裂缝发育的地质模型和模拟方程，模拟结果和现场结果吻合较好。王宁[13]对缝洞型油藏中天然微裂缝发育程度进行量化分析，发现裂缝发育程度对油藏压后产能有较大影响。在储层改造方面，曾青冬等[14]建立井筒和裂缝流体流动模型，讨论天然裂缝逼近角与水力裂缝扩展的关系。马耕等[15]进行了存在天然裂缝的水力压裂实验，实验结果表明天然裂缝形态对水力裂缝扩展有不同影响。宋晨鹏等[16]建立二维模型对水力裂缝和天然裂缝相交扩展进行模拟，发现小尺寸的天然裂缝对水力裂缝影响不大。林鹤等[17]利用微震检测和施工数据研究发现活动性天然裂缝相对于非活动性天然裂缝对压裂改造效果的影响更明显。赵金洲等[18]研究了水力裂缝沟通天然裂缝后不同因素对裂缝网络的影响，发现天然裂缝发育、地应力差小的环境有利于形成缝网。赵海洋等[19]、荣元帅等[20]对塔河缝洞型碳酸盐岩高效酸压改造进行了大量的现场论证和分析。
上述研究主要集中在天然裂缝连通对储层产能的影响，对于缝洞型碳酸盐岩油藏中天然裂缝如何影响压裂裂缝扩展及沟通溶洞的机理尚不明确，没有提出明确的技术思想。本文基于缝洞型碳酸盐岩储层特点，提出了循缝找洞改造储层思想，同时通过建立天然裂缝和溶洞储集体同时存在的压裂模型，研究单条天然裂缝及天然裂缝网络对压裂裂缝扩展轨迹的影响，对缝洞型油藏中压裂裂缝如何沟通溶洞进行模拟研究，并对研究成果进行推广应用。
-
考虑一个包含裂缝面
${\Gamma _{\rm{f}}}$ 的三维体$\Omega $ ，拉力F施加在外边界${\Gamma _{\rm{h}}}$ 上，边界的法向量为$\nu $ ，加载在准静态情况下，控制方程的有限元强形式为$$ \left\{ {\begin{array}{l}\nabla \sigma +\rho b=0{\text{，在}}\Omega {\text{上}}\\ \sigma \nu =t{\text{，在}}{\Gamma }_{\rm{h}}{\text{上}}\end{array}} \right.$$ (1) 式中，
$\sigma $ 为应力，Pa；$\;\rho $ 为材料密度，kg/m3；$b$ 为三维体的体力，N；$t$ 为边界单位法向载荷，Pa。在水力压裂过程中，在裂缝表面
${\Gamma _{\rm{f}}}$ 上，拉力满足在裂缝面两侧连续条件$${\rm{ - }}\sigma {f^ + }{\rm{ = }}\sigma {f^{\rm{ - }}} = p$$ (2) 式中，
$f$ 为裂缝面局部的单位法向量，无量纲；$p$ 为作用在裂缝面表面的水压，Pa。将总位移
$u$ 分解为连续部分${u_{\rm{c}}}$ 和非连续部分${u_{\rm{d}}}$ ，可以得到总位移为$$u(x) = {u_{\rm{c}}}(x){\rm{ + }}{u_{\rm{d}}}(x)$$ (3) 将裂缝面从正面到负面的位移差定义为
$$w = u({x^ + }) - u({x^ - })$$ (4) 代入公式(2)、(3)、(4)计算，得到裂缝表面能为
$\displaystyle\int_{{\Gamma _{\rm{f}}}} {pw{\rm{d}}\Gamma } $ 。根据变形体虚功原理，得到控制方程的有限元弱形式最终表示为
$$\int_\Omega {\delta \varepsilon \sigma {\rm{d}}\Omega } = \int_\Omega {\delta ub{\rm{d}}\Omega } {\rm{ + }}\int_{{\Gamma _{\rm{h}}}} {\delta uF{\rm{d}}\Gamma } + \int_{{\Gamma _{\rm{f}}}} {pw{\rm{d}}\Gamma } $$ (5) 式中，
$\delta $ 为开度，m；$\varepsilon $ 为应变；$u$ 为位移，m。 -
TOUGH2流体模拟采用基于积分有限差分的空间离散化技术(IFD)和Newton-Raphson方法，对流体各分量的质量平衡方程和各时刻的能量平衡方程进行线性化。在TOUGH2中，每一流体组分的总质量平衡由每一阶段的相关贡献进行累积。对于任意组分得到
$$\frac{\partial }{{\partial t}}{M^{\rm{K}}} - {Q^{\rm{K}}} = - \nabla \cdot \left( {q_{\rm{l}}^{\rm{K}} + q_{\rm{g}}^{\rm{K}}} \right)$$ (6) 式中，
${M^{\rm{K}}}$ 为组分K的单位体积质量，kg/m3；${Q^{\rm{K}}}$ 为组分K的能量通量，kg/(m3 · s)；$q_\Psi ^{\rm{K}}$ 为$\Psi$ 相中K组分的流量(液体$\Psi$ 用l表示，气体$\Psi$ 用g表示；空气K用a表示，水K用w表示)，kg/(m3 · s)。${M^{\rm{K}}}$ 的计算公式为$${M^{\rm{K}}} = \phi {S_{\rm{l}}}{\rho _{\rm{l}}}X_{\rm{l}}^{\rm{K}} + \phi {S_{\rm{g}}}{\rho _{\rm{g}}}X_{\rm{g}}^{\rm{K}}$$ (7) 式中，
$\phi $ 为孔隙度；$S$ 为饱和度；$\;{\rho}_{\Psi}$ 为$\Psi$ 相密度，kg/m3；$X_\Psi ^{\rm{K}}$ 为在$\Psi$ 相中K组分的质量分数。在
$\Psi$ 相中，每个组分的质量流量为达西流量和扩散流量之和$$q_\Psi ^{\rm{K}} = - {\rho _\Psi }X_\Psi ^{\rm{K}}\frac{{{k_{\rm{a}}}{k_{{\rm{r}}\Psi }}}}{{{\mu _\Psi }}}\left( {\nabla {p_\Psi } - {\rho _\Psi }a\nabla {\textit{z}}} \right) + i_\Psi ^{\rm{K}}$$ (8) 式中，
${k_{\rm{a}}}$ 为绝对渗透率，m2；${k_{{\rm{r}}\Psi }}$ 为$\Psi $ 相相对渗透率，%；${\mu _\Psi }$ 为$\Psi $ 相黏度，Pa · s；${p_\Psi }$ 为$\Psi $ 相压力，Pa；$a$ 为重力加速度，m/s2；$i_\Psi ^{\rm{K}}$ 为$\Psi$ 相中K组分的扩散通量，kg/(m3 · s)。用Fick定律计算扩散通量
$$i_\Psi ^{\rm{K}} = - \phi {S_\Psi }\tau {\rho _\Psi }{D_\Psi }\nabla X_\Psi ^{\rm{K}}$$ (9) 式中，
${S_\Psi }$ 为$\Psi $ 相饱和度；${D_\Psi }$ 为介质中依赖于温度、气压、介质曲折率和气体饱和度的有效分子扩散系数，在模拟过程中假设温度是恒定的；$\tau $ 为介质曲折率。利用上述方程及时间的一阶有限差分离散化，在时间
${e^{k + 1}} = {e^k} + \Delta e$ ，平衡方程重新定义为$$R_{}^{k + 1} = M_{}^{k + 1} - M_{}^k - \frac{{\Delta t}}{{{V_{}}}}\left\{ {\sum\limits_{{m}} {{A_{{m}}}q_{{m}}^{k + 1} + {V_{}}Q_{}^{k + 1}} } \right\}$$ (10) 式中，R为残差，kg/m3；M为单位体积质量，kg/m3；
$V$ 为单元体积，m3；A为单元面面积，m2；$q$ 为通量，kg/(m2 · s)；Q为源项，kg/(m3 · s)；k为迭代步数；m为离散单元号。对于非线性，采用牛顿-拉弗逊迭代法生成线性方程，在每个时间步中迭代，直到满足某个收敛条件。
-
AiFrac模拟器基于有限无网格法(FEMM)，有限无网格法主要用来模拟非连续模型[21-22]。如图1所示的横截面存在3种不同的单元类型：红色的裂缝单元、黄色的过渡单元和绿色的FEMM单元。裂缝穿过的单元称为裂缝单元，与裂缝单元相邻的单元称为过渡单元，剩下的单元为FEMM单元。存在两种结点：黑色实心的FE结点和黑色空心的PU结点，裂缝单元包含的所有结点为PU结点，FEMM单元包含的所有结点为FE结点，过渡单元同时包含PU结点和FE结点。对于任意点
$e = $ $ \left\{ {x,y,{\textit{z}}} \right\}$ ，在单元$\Omega $ 上，位移场$u(e)$ 可以由节点形函数和节点位移插值近似为
图 1 裂缝单元、过渡单元和FEMM单元的定义
Figure 1. Definition of fracture unit, transition unit and FEMN unit
$$u(e) = \sum\limits_{i = 1}^4 {{\omega _i}(e){u_i}(e)} $$ (11) 式中，
$u(e)$ 为FEMM中的全局近似，${\omega _i}(e)$ 为非负权重函数且权重函数总和等于1，${u_i}(e)$ 为结点$i$ 的局部近似。与FE结点相关的局部近似设为1，用最小二乘法构造在节点上的局部近似。对于断裂单元，应正确反映断裂表面的不连续位移场。根据Shepard公式建立分形单元的权函数，对于四面体单元
$\Omega $ 的节点检索集为$ {\zeta }_{\Omega }=\left\{1,2,3,4\right\}$ ，按照可视化准则，假设裂缝面是不可见的，已知单元内计算点$e$ ，定义节点对应的一组函数为$\varphi \left( e \right) = $ $ \left\{ {{\varphi _1}\left( e \right),{\varphi _2}\left( e \right),{\varphi _3}\left( e \right),{\varphi _4}\left( e \right)} \right\}$ ，对于计算点$e$ 对应的四面体单元每个节点的有限元形函数${\varphi _i}(e)$ 可以由体积坐标得到$$\left\{ \begin{gathered} {\varphi _1}(e) = \frac{{v_{ol}(E(e){E_2}{E_3}{E_4})}}{{v_{ol}({E_1}{E_2}{E_3}{E_4})}} \\ {\varphi _2}(e) = \frac{{v_{ol}(E(e){E_3}{E_4}{E_1})}}{{v_{ol}({E_1}{E_2}{E_3}{E_4})}} \\ {\varphi _3}(e) = \frac{{v_{ol}(E(e){E_4}{E_2}{E_1})}}{{v_{ol}({E_1}{E_2}{E_3}{E_4})}} \\ {\varphi _4}(e) = \frac{{v_{ol}(E(e){E_1}{E_2}{E_3})}}{{v_{ol}({E_1}{E_2}{E_3}{E_4})}} \\ \end{gathered} \right.$$ (12) 式中，
$v_{ol}({E_1}{E_2}{E_3}{E_4})$ 为四面体单元体积，$v_{ol}(E(e){E_i}{E_j}{E_k})$ 为任意点$E(e)$ 和其他3个顶点$\left\{ {{E_i},{E_j},{E_k}} \right\}$ 构成的四面体单元体积。当计算点$e$ 与单元节点的连线与裂缝面相交时，将其${\varphi _i}(e)$ 设为零函数。定义
${\zeta _\Omega }$ 为FEMM中与单元域$\Omega $ 相关的一个节点检索集，定义可见区$\zeta _\Omega ^{{\rm{vis}}}(e) \subset {\zeta _\Omega }$ 为计算点对应的所有非零加权函数的集合。对于包含裂缝面的任意单元$\Omega $ 中的计算点$e$ ，${\varphi _i}(e)$ 可以表示为$${\varphi '_i}(e) = \left\{ \begin{gathered} {\varphi _i}(e),e \in \zeta _\Omega ^{{\rm{vis}}} \\ 0,e \notin \zeta _\Omega ^{{\rm{vis}}} \\ \end{gathered} \right.$$ (13) 利用下式得出包含裂缝面与结点相关的有限元法单元权函数为
$${\omega _i}(e) = \frac{{{{\varphi '_i}}(e)}}{{{{\varphi '_1}}(e) + {{\varphi '_2}}(e) + {{\varphi '_3}}(e) + {{\varphi '_4}}(e)}}$$ (14) 过渡单元的体积权函数与标准有限元公式相同，传统的有限元单元形函数用标准有限元形函数计算。
-
TOUGH2-AiFrac耦合求解算法，将TOUGH2流体求解器[23-24]与AiFrac固体求解器[24]相结合，耦合过程见图2。流体流动和流体压力计算采用TOUGH2求解器，岩石变形和裂缝扩展计算采用AiFrac固体求解器。
图 2 TOUGH2-FEMM仿真原理
Figure 2. TOUGH2-FEMM simulation principle
TOUGH2和AiFrac建模是在不同的网格上进行的，但是几何结构和节点编号是相同的。如图3所示，TOUGH2分析每个单元的中心点，AiFrac需要在每个单元的8个节点上进行三维分析。
图 3 TOUGH和FEMM计算节点示意图
Figure 3. Sketch of TOUGH and FEMM calculate nodes
建模数据将TOUGH2单元的中心点插入到AiFrac单元的顶点，生成TOUGH2和AiFrac网格的过程见图4。首先，用TOUGH2建立了一套用于模拟裂缝和围岩基质中流体流动的六面体单元，然后，在之前六面体单元的基础上，使用AiFrac生成了一组四面体单元，用于模拟岩石变形和压裂，。
图 4 TOUGH网格构造出AiFrac网格
Figure 4. AiFrac network constructed on the basis of TOUGH network
-
研究缝洞型储层中单条天然裂缝对缝洞沟通的影响，建立如图5所示的模型。在35 m×35 m×1.5 m的假三维平板岩体模型中存在长10 m的水力裂缝，x方向为水平最大地应力方向，假设模拟岩体各向同性、线弹性，模型参数设置如下：岩石弹性模量30 GPa，泊松比0.2，孔隙度0.05，岩石渗透系数
$1 \times {10^{ - 18}}$ ，裂缝孔隙度0.25，裂缝渗透系数$3 \times {10^{ - 10}}$ 。在水力裂缝表面施加20 MPa的均匀压力。在假三维平板岩体模型中，将天然溶洞几何形状简化为圆柱体，圆柱体上下面圆心坐标为(25,17,0.75)，圆柱体厚度为1.5 m，溶洞内充满流体，压力为20 MPa。
图 5 单条天然裂缝-溶洞的沟通模拟模型
Figure 5. Simulation model of single natural fracture-vug connection
考虑到岩体网格尺寸引起的裂缝扩展敏感性，当岩体网格尺寸越小时，裂缝扩展模拟结果越接近于解析解，但同时会造成计算量大幅度增加。通过多组网格尺寸的模拟对比，选取h=0.5 m作为模拟网格尺寸。岩体四面体网格总数117 600，将圆柱体上下面圆心坐标(25,17,0.75)，半径4 m内框选的所有网格视为溶洞网格。
-
建立无天然裂缝的缝洞模型(图6)，基于TOUGH2-AiFrac耦合求解器计算水力裂缝扩展及缝洞沟通情况。模拟设置6个计算步，结果表明无天然裂缝时水力裂缝在压力作用下水平扩展，未与天然溶洞连通。
图 6 无天然裂缝的缝洞模型
Figure 6. Fracture-vug model without natural fracture
-
基于无天然裂缝的缝洞模型，在压裂裂缝扩展路径上增加一条与水平最大地应力方向夹角50°的天然裂缝，但裂缝走向偏离溶洞(图7)，研究天然裂缝对缝洞沟通的影响。假设天然裂缝未与水力裂缝连通时，缝内压力忽略不计。模拟设置3个计算步，结果表明压裂裂缝与天然裂缝连通，但未与溶洞连通(图8)。
图 7 天然裂缝偏离溶洞的缝洞模型
Figure 7. Fracture-vug model with natural fractures deviating from vugs
图 8 天然裂缝偏离溶洞时缝洞沟通结果
Figure 8. Fracture-vug connection result of natural fractures deviating from vugs
-
对比于上述天然裂缝偏离溶洞的情况，在压裂裂缝扩展路径上增加一条与水平最大地应力方向夹角43°的天然裂缝，裂缝走向指向溶洞(图9)。模拟设置3个计算步，结果表明压裂裂缝与天然裂缝连通后，与附近的溶洞储集体连通(图10)。
图 9 天然裂缝指向溶洞的缝洞模型
Figure 9. Fracture-vug model of natural fractures deviating towards vugs
图 10 天然裂缝指向溶洞时缝洞沟通结果
Figure 10. Fracture-vug connection result of natural fractures deviating towards vugs
-
天然裂缝作为缝洞型油藏中沟通溶洞储集体的重要通道，在储层中广泛存在。研究天然裂缝网络影响水力裂缝扩展及缝洞沟通，对于提高缝洞型油藏油气产量具有实际意义。储层中天然裂缝的特征表述包括长度L、面密度α、走向和倾角，考虑到研究区储层的特性，设置一组长度5±1 m，面密度2%，走向平行于z轴，倾角45±5°的天然裂缝网络。如图11所示在50 m×50 m×1.5 m的假三维平板岩体模型中存在长7 m的水力裂缝，x方向为水平最大地应力方向，在水力裂缝表面施加20 MPa的均匀压力，将天然溶洞几何形状简化为圆柱体，圆柱体圆心坐标(40,13,0.75)，厚度为1.5 m，溶洞内部充满流体，压力20 MPa，模型边界位移设置为0，模型参数取值与第2节相同。模型选取h=0.5 m作为岩体模拟网格尺寸，与第2节模型模拟网格相同。
图 11 天然裂缝网络-溶洞的缝洞沟通模型
Figure 11. Fracture-vug connection model of natural fracture network-vug
基于无天然裂缝网络时缝洞沟通模型，加入指向溶洞的天然裂缝网络(图12)。假设每条天然裂缝是独立存在的，将多条已经贯通的天然裂缝视为一条裂缝，不考虑天然裂缝中的流体压力。
图 12 天然裂缝网络指向溶洞时缝洞沟通结果
Figure 12. Fracture-vug connection result of natural fracture network deviating towards vugs
存在天然裂缝网络时，水力裂缝在压力作用下先沿着水平向扩展，之后与天然裂缝沟通，并将流体注入天然裂缝内，天然裂缝在压力作用下扩展，连通更多的天然裂缝，形成相互连通的裂缝网络，最终连通远端的溶洞储集体。
-
建立如图13所示的三维缝洞型油藏模型，模拟岩体尺寸50 m×50 m×50 m，位于岩体中心的圆形水力裂缝长5 m，x方向为水平最大地应力方向。设置两组长度5±2 m，走向平行于z轴，倾角分别为45°和135°的天然裂缝网络，裂缝面密度2%。在水力裂缝表面施加20 MPa的均匀压力，天然溶洞圆心坐标分别为(10,30,25)、(38,14,25)，溶洞内部充满流体，压力20 MPa，模型边界位移设置为0，模型参数设置与第2节相同。
图 13 三维天然裂缝网络-溶洞的缝洞沟通模型
Figure 13. Fracture-vug connection model of 3D natural fracture network-vug
计算结果显示，水力裂缝在压力作用下向四周扩展，连通上、下两条倾斜的天然裂缝，并注入流体，上、下两条天然裂缝在压力作用下扩展，连通更多的天然裂缝，形成相互连通的裂缝导流网络，最终沟通远端的溶洞储集体(图14)。
图 14 三维天然裂缝网络-溶洞的缝洞沟通结果
Figure 14. Fracture-vug connection result of 3D natural fracture network-vug
上述模拟研究表明，对于不在水平最大地应力方向上的溶洞，人工主裂缝难以对其进行有效沟通。为增加溶洞储集的沟通几率和数量，需充分利用储层中的天然裂缝，在不同方向形成相互连通的裂缝通道，实现对不同方向的多个溶洞进行沟通，即通过控制施工参数，使注入流体遵循原生天然裂缝的展布形态进行流动、连通，利用缝洞型碳酸盐岩储层缝-洞伴生发育的特点，沟通井周不同方向的溶洞储集体，即循缝找洞。
-
以循缝找洞思想为指导，制定储层改造技术方案，现场应用85井次，有效率100%，与常规酸压邻井相比，2020年增加原油产量16.12万t，增油效果显著。以W1井为例，该井区域油气富集，但生产无产量，拟通过储层改造实现增产，改造层段为奥陶系6 413.98~6 543.00 m，层位为
${{\rm O}_{2}yj} + {{\rm O}_{1 - 2}y}$ ，岩性为石灰岩，缝洞储层发育，且分布在井筒附近。施工过程中，通过在压裂初期注入压裂液形成一定长度的人工主裂缝，之后注入酸液充分激活人工主裂缝侧面及未端的天然裂缝，形成相互连通的裂缝通道，以增加沟通溶洞的几率。注入普通胶凝酸后在注入排量不变的条件下，压裂压力呈现较大幅度下降，显示人工裂缝与溶洞储集体存在沟通。W1井在压裂前基本无产量，压裂改造后产量呈现大幅度上升，平均日产量37.7 t，且持续稳产。
-
基于TOUGH2-AiFrac耦合求解算法建立了缝洞型油藏压裂裂缝扩展模型，研究了天然裂缝对缝洞沟通的影响机制,研究表明：缝洞型碳酸盐岩循缝找洞的储层改造思想是可行的；天然裂缝网络对压裂裂缝扩展轨迹具有重要影响，可充分利用天然裂缝的导向作用沟通溶洞；基于缝洞型碳酸盐岩循缝找洞储层改造思想，可实现不同方位溶洞储集体的沟通；现场应用证实缝洞型碳酸盐岩循缝找洞储层改造可大幅度提升井周油气动用程度，提高采收率。
Fracturing technology of searching for vugs along fractures in fractured-vuggy carbonate reservoirs
-
摘要: 深层及超深层缝洞型碳酸盐岩溶洞是油气的主要储集空间，压裂裂缝受地应力控制主要沿水平最大地应力方向扩展，并沟通该方向上的溶洞，其他方向的溶洞无法与井眼连通，导致井周储量动用程度低。根据碳酸盐岩储层伴生缝洞发育特点，提出循缝找洞的储层改造思想，即压裂过程中通过控制泵注压力注入流体，流体遵循原生天然裂缝的展布形态进行流动，形成相互连通的人工裂缝通道，沟通不同方向的多个溶洞，减少地应力对压裂裂缝形态的影响。通过构建含天然裂缝网络和溶洞的模型，采用TOUGH2-AiFrac耦合求解算法，研究天然裂缝数量、走向对压裂裂缝沟通溶洞的影响规律。结果显示天然裂缝走向对压裂裂缝扩展轨迹有较大影响，天然裂缝走向偏向溶洞有利于压裂裂缝沟通溶洞；不同数量、不同走向的天然裂缝网络能够有效连通不同方向的多个溶洞，说明能通过循缝找洞实现井周储量的高效动用。以循缝找洞思想为指导形成的技术方案现场应用85井次，2020年增加原油产量16.12万t，为深层及超深层缝洞型碳酸盐岩油气藏的改造增产提供了方向。Abstract: Deep and ultradeep fractured-vuggy carbonate vugs are the main oil and gas reservoir spaces. Under the control of in-situ stress, hydraulic fractures propagate mainly along the direction of the maximum horizontal stress to connect the vugs in this direction, while the vugs in the other directions cannot be connected with the boreholes, so the circumferential reserve producing degree is low. The reservoir stimulation idea of searching for vugs along fractures was put forward according to the development characteristics of the associated fractures and vugs of carbonate reservoirs. That is to inject the fluid by controlling the pumping pressure in the process of fracturing, so that the fluid can flow along the distribution forms of primary natural fractures to form interconnected hydraulic fracture pathways to connect multiple vugs in different directions and reduce the influences of in-situ stress on the morphology of hydraulic fractures. Then, the influential laws of the quantity and strike of natural fractures on the connection of vugs by hydraulic fractures were studied by using the TOUGH2-AiFrac coupling algorithm and constructing the model containing natural fracture network and vug. It is indicated that the strike of natural fractures has greater influence on the propagation tracks of hydraulic fractures, and when the strike of natural fractures is deviated towards vugs, it is favorable for the hydraulic fractures to connect vugs. In addition, natural fracture networks of different quantities and strikes can effectively connect multiple vugs in different directions, indicating that efficient production of circumferential reserves can be realized by searching for vugs along fractures. The technical scheme formulated according to the idea of searching for vugs along fractures has been applied on site 85 well times and the incremental oil in 2020 is 16.12×104 t. The research results provide the direction for the stimulation of deep and ultradeep fractured-vuggy carbonate oil and gas reservoirs.
-
Key words:
- ultradeep /
- carbonate rock /
- fracture and vug /
- fracturing /
- searching for vugs along fractures
-
图 1 裂缝单元、过渡单元和FEMM单元的定义
Figure 1. Definition of fracture unit, transition unit and FEMN unit
图 4 TOUGH网格构造出AiFrac网格
Figure 4. AiFrac network constructed on the basis of TOUGH network
图 5 单条天然裂缝-溶洞的沟通模拟模型
Figure 5. Simulation model of single natural fracture-vug connection
图 7 天然裂缝偏离溶洞的缝洞模型
Figure 7. Fracture-vug model with natural fractures deviating from vugs
图 8 天然裂缝偏离溶洞时缝洞沟通结果
Figure 8. Fracture-vug connection result of natural fractures deviating from vugs
图 9 天然裂缝指向溶洞的缝洞模型
Figure 9. Fracture-vug model of natural fractures deviating towards vugs
图 10 天然裂缝指向溶洞时缝洞沟通结果
Figure 10. Fracture-vug connection result of natural fractures deviating towards vugs
图 11 天然裂缝网络-溶洞的缝洞沟通模型
Figure 11. Fracture-vug connection model of natural fracture network-vug
图 12 天然裂缝网络指向溶洞时缝洞沟通结果
Figure 12. Fracture-vug connection result of natural fracture network deviating towards vugs
图 13 三维天然裂缝网络-溶洞的缝洞沟通模型
Figure 13. Fracture-vug connection model of 3D natural fracture network-vug
-
[1] 苏建政, 张素青. 碳酸盐岩非均质性对水力裂缝起裂影响的有限元法研究[J]. 石油钻采工艺, 2008, 30(4):63-66. doi: 10.3969/j.issn.1000-7393.2008.04.016 SU Jianzheng, ZHANG Suqing. Studying the effect of carbonatite heterogeneity on hydraulic fracture initiation by using finite element method[J]. Oil Drilling & Production Technology, 2008, 30(4): 63-66. doi: 10.3969/j.issn.1000-7393.2008.04.016 [2] 邢静, 李勇明. 碳酸盐岩油藏酸压模型研究及应用[J]. 石油钻采工艺, 2009, 31(3):62-66. doi: 10.3969/j.issn.1000-7393.2009.03.016 XING Jing, LI Yongming. Model study of acid fracturing model for carbonate reservoirs[J]. Oil Drilling & Production Technology, 2009, 31(3): 62-66. doi: 10.3969/j.issn.1000-7393.2009.03.016 [3] 罗源, 杨乾龙, 刘平礼, 等. 裂缝性碳酸盐岩酸化模拟新方法[J]. 石油钻采工艺, 2015, 37(5):73-77, 94. doi: 10.13639/j.odpt.2015.05.018 LUO Yuan, YANG Qianlong, LIU Pingli, et al. A new simulation method for fractured carbonate acidizing[J]. Oil Drilling & Production Technology, 2015, 37(5): 73-77, 94. doi: 10.13639/j.odpt.2015.05.018 [4] LIU Z, TANG X, TAO S, et al. Mechanism of connecting natural caves and wells through hydraulic fracturing in fracture-cavity reservoirs[J]. Rock Mechanics and Rock Engineering, 2020, 53: 5511-5530. doi: 10.1007/s00603-020-02225-w [5] 樊建明, 陈小东, 雷征东, 等. 鄂尔多斯盆地致密油藏天然裂缝与人工裂缝特征及开发意义[J]. 中国石油大学学报(自然科学版), 2019, 43(3):98-106. doi: 10.3969/j.issn.1673-5005.2019.03.011 FAN Jianming, CHEN Xiaodong, LEI Zhengdong, et al. Characteristics of natural and hydraulic fractures in tight oil reservoir in Ordos Basin and its implication to field development[J]. Journal of China University of Petroleum(Edition of Natural Science), 2019, 43(3): 98-106. doi: 10.3969/j.issn.1673-5005.2019.03.011 [6] 姚军, 王子胜, 张允, 等. 天然裂缝性油藏的离散裂缝网络数值模拟方法[J]. 石油学报, 2010, 31(2):284-288. doi: 10.7623/syxb201002018 YAO Jun, WANG ZI Sheng, ZHANG Yun, et al. Numerical simulation method of discrete fracture network for naturally fractured reservoirs[J]. Acta Petrolei Sinica, 2010, 31(2): 284-288. doi: 10.7623/syxb201002018 [7] 苏玉亮, 王霞, 李涛, 等. 天然裂缝对油藏开采特征的影响[J]. 钻采工艺, 2006, 29(6):41-43. doi: 10.3969/j.issn.1006-768X.2006.06.014 SU Yuliang, WANG Xia, LI Tao, et al. Effect of natural fracture on the low permeability reservoir production characteristics[J]. Drilling & Production Technology, 2006, 29(6): 41-43. doi: 10.3969/j.issn.1006-768X.2006.06.014 [8] 李博. 致密油藏分段压裂水平井嵌入式离散裂缝模型数值模拟研究[D]. 成都: 西南石油大学, 2018. LI Bo. Numerical study of multiple-stage fractured horizontal well in tight oil reservoirs based on embedded discrete fracture model[D]. Chengdu: Southwest Petroleum University, 2018. [9] 孙若凡. 考虑天然裂缝影响的致密油藏压裂水平井产能模型研究[D]. 北京: 中国地质大学(北京), 2020. SUN Ruofan. Study on productivity model of fracturing horizontal well in tight reservoirs considering the effect of natural fractures[D]. Beijing: China University of Geosciences (BeiJing), 2020. [10] 朱圣举, 赵向原, 张皎生, 等. 低渗透砂岩油藏天然裂缝开启压力及影响因素[J]. 西北大学学报(自然科学版), 2016, 46(4):573-578. doi: 10.16152/j.cnki.xdxbzr.2016-04-020 ZHU Shengju, ZHAO Xiangyuan, ZHANG Jiaosheng, et al. Fracture opening pressure and its influence factors in low-permeability sandstone reservoirs[J]. Journal of Northwest University(Natural Science Edition), 2016, 46(4): 573-578. doi: 10.16152/j.cnki.xdxbzr.2016-04-020 [11] 刘熙远. 考虑天然裂缝的特低渗透油藏时变数值模拟研究及应用[D]. 青岛: 中国石油大学(华东), 2018. LIU Xiyuan. Research and application of ultra-low permeability reservoirs time-varying numerical simulation considering natural fractures[D]. Qingdao: China University of Petroleum (East China), 2018. [12] 吴忠宝, 胡文瑞, 宋新民, 等. 天然微裂缝发育的低渗透油藏数值模拟[J]. 石油学报, 2009, 30(5):727-730, 734. doi: 10.3321/j.issn:0253-2697.2009.05.016 WU Zhongbao, HU Wenrui, SONG Xinmin, et al. Numerical simulation on low-permeability reservoir with abundant natural micro-fractures[J]. Acta Petrolei Sinica, 2009, 30(5): 727-730, 734. doi: 10.3321/j.issn:0253-2697.2009.05.016 [13] 王宁. Erdos盆地低渗透油藏天然裂缝发育程度量化分析及应用[D]. 西安: 西安石油大学, 2017. WANG Ning. Quantitative analysis and application of natural fracture development degree in low permeability reservoirs in Erdos Basin[D]. Xi’an: Xi’an Shiyou University, 2017. [14] 曾青冬, 姚军, 孙致学. 页岩气藏压裂缝网扩展数值模拟[J]. 力学学报, 2015, 47(6):994-999. doi: 10.6052/0459-1879-15-014 ZENG Qingdong, YAO Jun, SUN Zhixue. Numerical modeling of fracture network propagation in shale reservoirs[J]. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(6): 994-999. doi: 10.6052/0459-1879-15-014 [15] 马耕, 张帆, 刘晓, 等. 裂缝性储层中水力裂缝扩展规律的试验研究[J]. 采矿与安全工程学报, 2017, 34(5):993-999. doi: 10.13545/j.cnki.jmse.2017.05.025 MA Geng, ZHANG Fan, LIU Xiao, et al. Experimental study on hydraulic fracture propagation in fractured reservoir[J]. Journal of Mining and Safety Engineering, 2017, 34(5): 993-999. doi: 10.13545/j.cnki.jmse.2017.05.025 [16] 宋晨鹏, 卢义玉, 夏彬伟, 等. 天然裂缝对煤层水力压裂裂缝扩展的影响[J]. 东北大学学报(自然科学版), 2014, 35(5):756-760. doi: 10.3969/j.issn.1005-3026.2014.05.034 SONG Chenpeng, LU Yiyu, XIA Binwei, et al. Effects of natural fractures on hydraulic fractures propagation of coal seams[J]. Journal of Northeastern University (Natural Science), 2014, 35(5): 756-760. doi: 10.3969/j.issn.1005-3026.2014.05.034 [17] 林鹤, 李德旗, 周博宇, 等. 天然裂缝对压裂改造效果的影响[J]. 石油地球物理勘探, 2018, 53(增刊2):156-161, 167. doi: 10.13810/j.cnki.issn.1000-7210.2018.S2.024 LIN He, LI Deqi, ZHOU Boyu, et al. Influences of natural cracks on fracturing[J]. Oil Geophysical Prospecting, 2018, 53(S2): 156-161, 167. doi: 10.13810/j.cnki.issn.1000-7210.2018.S2.024 [18] 赵金洲, 李勇明, 王松, 等. 天然裂缝影响下的复杂压裂裂缝网络模拟[J]. 天然气工业, 2014, 34(1):68-73. doi: 10.3787/j.issn.1000-0976.2014.01.010 ZHAO Jinzhou, LI Yongming, WANG Song, et al. Simulation of a complex fracture network influenced by natural fractures[J]. Natural Gas Industry, 2014, 34(1): 68-73. doi: 10.3787/j.issn.1000-0976.2014.01.010 [19] 赵海洋, 林涛, 陈朝刚. 塔河油田碳酸盐岩油藏重复酸压可行性研究[J]. 石油钻采工艺, 2009, 31(1):113-116, 120. doi: 10.3969/j.issn.1000-7393.2009.01.029 ZHAO Haiyang, LIN Tao, CHEN Chaogang. Study on acid-refracturing feasibility in carbonate reservoirs of Tahe oilfield[J]. Oil Drilling & Production Technology, 2009, 31(1): 113-116, 120. doi: 10.3969/j.issn.1000-7393.2009.01.029 [20] 荣元帅, 何新明. 塔河油田碳酸盐岩缝洞型油藏大型复合酸压选井优化论证[J]. 石油钻采工艺, 2011, 33(4):84-87. doi: 10.3969/j.issn.1000-7393.2011.04.022 RONG Yuanshuai, HE Xinming. Well selection optimization verification of large scale compound acid fracturing in carbonate fracture-cavity reservoir in Tahe oilfield[J]. Oil Drilling & Production Technology, 2011, 33(4): 84-87. doi: 10.3969/j.issn.1000-7393.2011.04.022 [21] LIU Q, SUN L, TANG X, et al. Simulate intersecting 3D hydraulic cracks using a hybrid “FE-Meshfree” method[J]. Engineering Analysis With Boundary Elements, 2018, 91: 24-43. doi: 10.1016/j.enganabound.2018.03.005 [22] TANG X H, RUTQVIST J, HU M S, et al. Modeling three-dimensional fluid-driven propagation of multiple fractures using TOUGH-FEMM[J]. Rock Mechanics and Rock Engineering, 2019, 52(2): 611-627. doi: 10.1007/s00603-018-1715-7 [23] 薛志波. 井筒液位监控系统中数控声源部分的研究[D]. 成都: 西南石油大学, 2007. XUE Zhibo. The research of the fluid position supervisory system acoustic source[D]. Chengdu: Southwest Petroleum University, 2007. [24] XU J, TANG X, WANG Z, et al. Investigating the softening of weak interlayers during landslides using nanoindentation experiments and simulations[J]. Engineering Geology, 2020, 277: 105801. doi: 10.1016/j.enggeo.2020.105801 -
点击查看大图
图(14)
计量
- 文章访问数: 11
- HTML全文浏览量: 103
- PDF下载量: 36
- 被引次数: 0
