1. 引言

天然气作为一种清洁低碳的化石能源，是化石能源向新能源转换的重要过渡能源，在全球能源绿色低碳转型中发挥重要作用 [1] 。仅2022年世界天然气消费量仍高达3.94万亿立方米。页岩气作为一种非常规天然气资源，是一种清洁、高效的能源资源，资源储量丰富。其大规模开发利用有助于降低我国对于石油、天然气等传统能源的进口依赖，提高能源自主保障能力、保障国家能源安全。

目前我国技术可采页岩气储量可达8.82 × 10⁸ m³ [2] ，在涪陵、长宁、威远等地区，以实现较为完善的工业化开发 [3] 。然而，页岩储层通常具有低孔、超低渗等特征，导致页岩气在储层内流动十分困难，目前水平井分段压裂技术已成为页岩气开发的一种关键手段，须采用“体积压裂 + 水平井”技术才能对页岩气进行商业开发 [4] 。国内外专业学者对页岩气生产规律、递减特征、产能主控因素等开展了大量工作研究，发现页岩气生产压力具有初期递减快、后期递减慢的特征 [5] ；产能影响因素有优质储层发育厚度、钻遇长度、储层改造程度、Langmuir体积、应力敏感系数等，导致如何准确预测页岩气井最终可采储量(EUR)目前仍是业界难点 [6] 。目前，常用的页岩气井产量计算方法有解析法、数值模拟法、经验法、物质平衡法。其中，数值模拟法是一种基于对连续性方程进行数值化求解，从而模拟地下流体渗流规律并进行油气藏动态预测的方法，模型推导难度大，涉及的物理参数多，计算难度高，而且无法表征产量与压力间的函数关系 [7] 。解析法是一种基于对连续性方程进行线性化处理，采用数学手段进行解析求解，从而实现产量预测的方法 [8] [9] 。模型解析求解难度通常较大，且通常须对非线性问题进行线性化处理。物质平衡法则是基于质量守恒定律建立的，即原始地质储采量等于采出量与地下剩余量之和，通常只研究采出气量与地层压力之间的关系 [10] 。经验法是一种基于矿场实际生产数据变化特征，采用经验模型进行回归的方法，但因为参数较多且以试算得到，没有累计气量的直接计算公式，且误差较大 [11] 。

基于上述分析，现有的页岩气产量计算方法均有一定不足。因此，本文将基于页岩气渗流阻力理论，结合物质平衡方法，建立了页岩气压裂水平井产量预测模型。该模型物理意义明确，推导简便，且易于计算，能准确地表征页岩气井生产动态特征。

2. 页岩气物质平衡方程的建立

由于页岩气藏孔隙度和渗透率极低，且具有自生自储特征，因此，可将页岩气藏视为封闭气藏。根据物质守恒原理，页岩气藏的物质平衡方程(气体体积守恒)的基本形式为：

(1)

式中： $G_f^i$ 为页岩储层原始游离气量，m³； $G_a^i$ 为页岩储层原始吸附气量，m³；G_f为页岩储层当前游离气量，m³；G_a为页岩储层当前吸附气量，m³；G_p为页岩储层累计产出气量，m³。

通过考虑页岩气体积、孔隙体积和束缚水的体积，可得原始条件下页岩气藏容积为：

$V_{ci}^{}=V_p^{}-V_{wc}^{}=A_g^{}h\varphi \left(1-S_{wc}\right)$ (2)

页岩气藏原始吸附气体积为：

$G_a^i=\frac{{\rho }_s^{}{V}_b^{}}{1000}{V}_a^i$ (3)

式中：A_g为气藏面积，m²；h为气藏厚度，m； $\varphi$ 为孔隙度，无因次；S_wc为束缚水饱和度，无因次；V_p为气藏孔隙体积，m³；V_wc为气藏束缚水体积，m³。

当气藏产出一定气量G_p之后，地层压力从p_i下降到某一压力p时，页岩气藏吸附气体积：

$G_a^{}=\frac{{\rho }_s^{}{V}_b^{}}{1000}{V}_a^{}$ (4)

其中，页岩气吸附量用Langmuir模型进行表示：

$V_a^{}=\frac{V_L^{}p}{p_L^{}+p}$ (5)

气藏孔隙体积减小量：

$\Delta V_p^{}=V_p^{}{c}_p^{}\Delta p$ (6)

束缚水体积膨胀量：

$\Delta V_{wc}^{}=V_p^{}{c}_w^{}\Delta p$ (7)

其中

$\Delta p=p_i^{}-p$ (8)

因此，页岩气藏当前体积为：

$V_c^{}=V_{ci}^{}\left(1-\frac{c_p^{}+c_w^{}{S}_{wc}^{}}{1-S_{wc}^{}}\Delta p\right)$ (9)

式中：ρ_s为岩石密度，kg/m³；V_b为气藏体积，m³；V_a为页岩储层当前吸附气量，m³；c_p为岩石压缩系数，1/MPa；c_w为地层水压缩系数，1/MPa。

封闭气藏中的原始气量为：

$G_f^i=\frac{V_f^i}{B_{gi}^{}}=\frac{V_{ci}^{}}{B_{gi}^{}}$ (10)

其中，体积系数可表示为：

$B_{gi}^{}=\frac{p_{sc}^{}}{Z_{sc}^{}{T}_{sc}^{}}\frac{Z_i^{}{T}_i^{}}{p_i^{}}$ (11)

当前地面游离气量体积为：

$G_f^{}=\frac{V_c^{}}{B_g^{}}=\frac{V_{ci}^{}}{B_g^{}}\left(1-\frac{c_p^{}+c_w^{}{S}_{wc}^{}}{1-S_{wc}^{}}\Delta p\right)$ (12)

其中，体积系数可表示为：

$B_{gi}^{}=\frac{p_{sc}^{}}{Z_{sc}^{}{T}_{sc}^{}}\frac{Z_i^{}{T}_i^{}}{p}$ (13)

将式(3)、式(4)、式(10)和式(12)代入式(1)，并整理得：

$\frac{p}{Z}\left(1-c_c^{}\Delta p\right)=\frac{p_i^{}}{Z_i^{}}\left(1-\frac{G_p^f}{G_f^i}\right)$ (14)

其中，

$G_p^f=G_p^{}-\frac{{\rho }_s^{}{V}_b^{}}{1000}\left(V_a^i-V_a^{}\right)$ (15)

$c_c^{}=\frac{c_p^{}+c_w^{}{S}_{wc}^{}}{1-S_{wc}^{}}$ (16)

式中：B_gi为原始条件下气体体积压缩系数，1/MPa； $V_f^i$ 为原始地下游离气量；Z_sc为标准状态压因子，无因次；T_sc为标准状态温度，K；P_sc为标准状态压力，MPa；Z_i为气藏当前地层压缩因子，无因次；T_i为气藏原始地层温度，K；P_i为气藏原始地层压力，MPa。

3. 页岩气压裂水平井产量预测模型研究

根据渗流阻力法，并利用面积等值原则，建立压裂水平井产能公式为：

$q_{sc}^{}=\frac{\pi Nh{Z}_{sc}^{}{T}_{sc}^{}\left[\psi \left(P\right)-\psi \left(P_{wf}\right)\right]}{P_{sc}^{}T\left[\frac{1}{k_m^{}}\ln \sqrt{\frac{x_e^{}}{2x_f^{}}+\frac{1}{k_f^{}}\ln \left(\frac{1}{r_{wef}^{}}\sqrt{\frac{2L{x}_f^{}}{\pi N}}\right)}\right]}$ (17)

其中，

$\psi \left(P\right)=2{\displaystyle {\int }_{P_i^{}}^P\frac{P}{\mu Z}\text{d}P}$ (18)

式中：L为压裂水平井有效长度，m；N为压裂段数，段；x_e为储层宽度，m；k_f为裂缝渗透率，mD；k_m为储层渗透率，mD；x_f为裂缝半长，m；q_sc为地面条件下日产气量，10⁴ m³/d；ψ为地层压力对应的拟压力，MPa²/(mPa·s)；ψ (P_wf)为井底流压对应的拟压力，MPa²/(mPa·s)。

4. 页岩气压裂水平井产量预测模型校验与分析

4.1. 页岩气压裂水平井产量预测模型的校验

选取四川盆地某页岩气井真实生产数据，首先通过该井实际日产气量，利用本文所建的产量预测模型反算出井底流压，然后根据井底流压的拟合情况完成所建模型的可靠性分析。拟合结果如图1所示，新建的产量预测模型能够较好地描述出真实气井的井底流压变化特征。因此，本文所建的产量预测模型是可靠的。

4.2. 页岩气压裂水平井产能影响因素分析

1) 吸附作用

页岩气主要以游离态和吸附态两种形式赋存于页岩储层中，不同地质条件下，页岩气的赋存状态有较大差异 [12] 。吸附气是页岩气重要的赋存形式之一。页岩气吸附规律常采用Langmuir模型进行描述，因此，本文研究了Langmuir体积和Langmuir压力对于页岩气井产量影响。由图2可知，页岩气井稳产期和累产气量均随着Langmuir体积和Langmuir压力的增加而增加，不考虑吸附作用会显著低估气井产量，并且Langmuir体积对页岩气井产量的影响要高于Langmuir压力。(由图2和图3可知)

2) 裂缝段数

多级压裂水平井技术已是页岩气藏成功开发的核心技术 [13] ，但由于地层条件和压裂工艺的复杂性导致压裂裂缝呈现出不同形态，对页岩气产量造成影响 [14] 。由图4可知，当裂缝段数由10段分别增至13段和15段时，稳产期由202天分别增至942天和1895天，气井累计产量由1.27 × 10⁸ m³分别增至1.89 × 10⁸ m³和2.48 × 10⁸ m³，增大倍数依次为0.49倍和0.96倍。因此，压裂段数的增加可促进页岩气井产量的增加，但随着气井持续生产，相较于压裂段数对气井早期产量的影响，压裂段数对气井中后期产量的影响将逐渐变小。

3) 裂缝半长

由图5可知，当裂缝半长由30 m分别提高至60 m、90 m和120 m时，稳产期由407天分别增至942天、1648天和2179天，增大倍数分别为1.31倍、3.05倍和4.35倍，气井累计产量由1.25 × 10⁸ m³分别增至1.89 × 10⁸ m³、2.41 × 10⁸ m³和2.95 × 10⁸ m³，增大倍数依次为0.51倍、0.93倍和1.35倍。因此，在不考虑经济约束条件下，增大裂缝半长将有利于提高气井产能。

5. 结论

本文基于物质平衡方法建立了页岩气压裂水平井产量预测模型，并且通过真实页岩气井生产数据对该模型进行了校验。在此基础上，完成了吸附作用、压裂半长和压裂段数对页岩气井产量的影响规律研究。计算分析表明：

1) Langmuir体积对页岩气井产量的影响要高于Langmuir压力，若不考虑页岩气吸附作用则会显著低估气井产量。因此，考虑页岩气吸附作用并准确获取吸附参数是提高页岩气井产量预测结果可靠性的关键因素之一。

2) 页岩气井产量与裂缝半长和压裂段数都呈正相关，但裂缝半长对于页岩气井整个生产阶段都有较大的影响，而压裂段数对气井产量的影响会随生产时间的增加而逐渐减小，因此在页岩气储层改造时应充分发挥二者的协同效应。

基金资助

重庆科技学院大学生科技创新训练计划项目(项目编号：2023144)；重庆市自然科学基金面上项目(项目编号：CSTB2023NSCQ-MSX0264)；重庆市教委科学技术研究项目青年项目(项目编号：KJQN202201517)；重庆科技学院人才引进科研启动资助项目(项目编号：ckrc2022025)。

参考文献