预测页岩储层体积改造后生产效果的方法
【技术领域】
[0001 ]本申请属于油气田开发领域，具体地说，涉及一种预测页岩储层体积改造后生产 效果的方法。
【背景技术】
[0002] 随着全球能源工业的快速发展和能源需求的与日倶增，越来越多的低渗透油气藏 投入开发利用。特低渗透气藏具有储层渗透率极低、存在平面非均质性和启动压力梯度、应 力敏感效应等特点。压裂水平井技术是特低渗透气藏增储上产和经济有效开发的关键技术 之一，但常规压裂形成单一裂缝的改造模式不能有效地沟通储层，难以实现高效商业化开 发。针对特低渗透致密气藏开展体积压裂能够形成复杂缝网，从而显著提高储层的整体渗 透率，实现大幅度增产。在特低渗透致密气藏压裂水平井生产过程中，随着地层压力不断下 降，不同尺度气孔、溶蚀孔洞、孔隙喉道等容易发生收缩变形，引起储层渗透率的应力敏感 效应。在对特低渗透气藏压裂水平井体积压裂的产量进行预测时，必须综合考虑渗透率非 均质性、启动压力梯度、储层压力敏感和体积压裂形成缝网效应等因素的综合影响。
[0003] 目前针对水平井体积压裂的产量预测模型主要有双线性流、三线性流、五线性流 等渗流模型，上述模型均没有同时考虑特低渗透气藏非均质性、特低渗透气藏流体流动时 存在启动压力梯度和压力下降时产生储层应力敏感效应的综合影响。
【发明内容】
[0004] 有鉴于此，本申请所要解决的技术问题是现有技术没有同时考虑特低渗透气藏非 均质性、特低渗透气藏流体流动时存在启动压力梯度和压力下降时产生储层应力敏感效应 的综合影响。
[0005] 为了解决上述技术问题，本申请公开了一种预测页岩储层体积改造后生产效果的 方法，主要包括以下步骤：
[0006 ] 1)收集储层、流体、水平井筒的基本参数；
[0007] 2)将非均质气藏沿所述水平井筒的长度方向划分成至少两个渗流带，所述渗流带 含有复杂体积裂缝的人工裂缝；
[0008] 3)建立每一个所述渗流带的渗流单元模拟模型；
[0009] 4)建立至少两个所述渗流单元的渗流模拟模型；
[0010] 5)耦合各个所述渗流单元的渗流的模拟模型，预测非均质页岩气藏压裂水平井的 稳态条件下生产效果；
[0011] 6)将时间离散化，重复步骤1~5,预测非均质页岩气藏非稳态条件下生产效果。
[0012] 进一步的，所述步骤1)中，所述储层的基本参数包括:原地应力方向、储层厚度、储 层温度、孔隙度、渗透率、基质渗透率变形因子;所述流体的基本参数包括:气体粘度、气体 临界压力、气体偏差因子、气体临界温度;所述水平井筒的基本参数包括:水平井方位、水平 井井筒长度。
[0013] 进一步的，所述步骤2)具体为:根据储层在水平井筒长度方向上渗透率的差异，将 渗透率相同的储层划分成同一个含有人工裂缝渗透带的渗流带。
[0014] 进一步的，所述渗流单元的基本流动过程包括:基质线性流、缝网高渗透带线性 流、当量井径向流和裂缝内的线性渗流；以第1个条带基质中的第1个渗流单元为例计算渗 流单元中各流动过程的阻力。
[0015] 进一步的，所述步骤3)具体包括：
[0016] a、根据关系式： 计算I区线性流的阻 ， 力，式中：R叫为第i渗透带中第j个渗流单元I区基质线性流阻力与缝网线性流阻力之和， MPa2/( (m3/d) (mPa · s)); T为气藏温度，单位为K; Tsc为标准大气温度，单位为K; psc为标准大 气压，单位为MPa; I^为第i渗透带中第j条裂缝与上一条裂缝渗流距离，m; Xf为裂缝半长，单 位为m; ki为第i渗透带基质渗透率，单位为mD ;h为储层厚度，单位为m;
[00Π ] b、根据关系式
，计算II区径向流阻力;式中： Rnij为第i渗透带中第j个渗流单元当量井径向流阻力;单位为MPa2/((m3/d)(mPa · S)) ;kfij 为裂缝渗透率，单位为Dc; Wfij为裂缝宽度，单位为cm;
[0018] c、根据所述的I区线性流的阻力和II区径向流阻力，建立所述渗流带的渗流单元 的模拟模型。
[0019]进一步的，所述步骤&中，所述关系式：
[0022]根据气藏启动压力梯度下的基质线性流动阶段、高渗透缝网线性流阶段流量与拟 压力，得出公式：
(1)
[0024]根据拟压力定义，气体流过距离d时消耗的启动压力为Glj · d，将Glj · d定义为启 动压力PB;因此可做以下转化：
[0025] CN 105507893 A 说明书 3/17 页 (2) t t
[0026] 式中：也ij、itdij分别为I区基质线性流外拟压力和内拟压力，Ifeuij-Ifedij为流体流过 长度为dijl区基质的启动拟压差，单位为MPa2/(mPa · s) ;Gij为I区基质的启动压力梯度，单 位为MPa/m;y为天然气黏度，单位为mPa · s;z为天然气偏差因子，无因次;dij为I区基质线性 流长度，单位为m; dnij为I区缝网线性流长度，单位为m; qij为流量，单位为m3/d; Xf为裂缝半 长，单位为m; ki为第i渗透带基质渗透率，单位为mD; kni为第i渗透带缝网渗透率，单位为mD; h为储层厚度，单位为m; T为气藏温度，单位为K; Ts。为标准大气温度，单位为K; ps。为标准大 气压，单位为MPa。
[0027]进一步的，所述步骤b中：
[0028]根据同一个渗流单元内，I区和II区边界位置处的压力相等，I区线性流的内压即 为II区径向流的外压，则径向流流量与压差关系式为：
(5)
[0030] 式中:Ikfij为II区径向流内拟压力，单位为MPa2/(mPa · s);rwij为当量井半径;单位 为m;
[0031] 所述当量直井的半径5通过当量井径模型求取，具体为：
[0032] 通过建立裂缝产量和直井产量相等的关系式从而得到当量井径与裂缝半长之间 关系式，下面分别建立求取当量直径和裂缝的产量：
[0033] 直井产量:求取定压供给边界条带气藏中一口直井六:的稳态产量，根据镜像反映 原理，将其转化为一排生产井和一排注气井的产量稳态进行求解；
[0034] 根据等倌渗流阻力法，得一口井产量：
(6)
[0036] 有限导流裂缝产量：
[0037]
(J)
[0038] 式中：Ψ为裂缝方向上x点处拟压力，单位为MPa2/(mPa · s);kfij为裂缝渗透率，单 位为Dc; Wf ij为裂缝宽度，单位为cm;
[0039] 解出裂缝内平均压力近似为裂缝内压力，流体从基质流到裂缝内满足线性流动规 律，线性流外拟压力为IK111，求出基质线性流流暈，再由质量守恒得裂缝产量：
(8)
[0041]对比当量直井的产量式(6)与有限导流裂缝的产量式(8)得当量井半径公式：
(9)
[0043] 式中，rw为当量井半径，单位为m;
[0044]将式(8)代入式(5)得到渗流单元II区当量井附近等效径向流阻力为：
C10) 〇
[0046] 进一步的，所述步骤4)具体为：
[0047] 在求得各区流动阻力之后，根据不同渗流区域间满足压力连续、流量相等的基本 原理，利用等值渗流阻力法即求得压裂水平井产能公式。
[0048]进一步的，所述步骤5)具体为：
[0049] 流体由两边向中间流动，先分析两边渗透带裂缝产能方程再分析中间裂缝产能方 程，相邻裂缝间用等值渗流阻力法连接各阻力得到
个方程组成的线性方程组；
[0050] 计算第1渗透带中裂缝产量，公式为：
(π)
[0052]计算第3渗透带中裂缝产量，公式为：
(12)
[0054]由于第2渗透带最外两个渗流单元分别与1、3渗透带相连，则其I区基质线性流阻 力为：
(13) (14)
[0057] 假设中间分流裂缝为第2渗透带中第m条裂缝所等效的井排，该井排为左右两流量 汇合处，所以可列出N2+l个方程：
[0058] (15)
)
[0059] 式中：q2m为左边流向分流井排的流量，单位为m3/d; q2U+i)为右边流向分流井排的 流量，分流裂缝的实际流量为两者之和，单位为m3/d;
[0060] 同理，当存在任意个非均质渗透带时，方程也可由此得出，此时第N个非均质渗透 带，同时I <N〈Nm，的渗流方程为：
(16)
[0062]根据拟压力定义，拟启动拟压差可采用压力平方差表示：
(17) (18)
[0065]同理可得地层中任意一点压力p的拟压力函数表达式：
(19)
[0067]基于以上推导过程，得到一个(N+l) X (N+1)的封闭线性方程组，N为裂缝条数，采 用高斯消元法求解得到各裂缝产量，最终叠加得到稳态条件下压裂水平井产能；
[0068]压裂水平井总的产量为：
(20)
[0070] 式中，Q为压裂水平井产量，单位为m3/d。
[0071] 进一步的，所述步骤6)具体包括:求解特低渗透非均质气藏在非稳态条件下