📚 课程目录第3章递减曲线分析基础

3. 递减曲线分析基础

递减曲线分析（Decline Curve Analysis, DCA）是页岩气井产量预测与EUR估算中最经典、最常用的方法之一。它的核心思想非常朴素：把历史产量数据按时间排列，寻找一条能够描述产量随时间下降规律的数学曲线，再用这条曲线外推到经济极限产量，从而估算最终可采储量。Arps在1945年提出的三条递减方程——指数递减、双曲递减和调和递减——构成了这一方法体系的数学基础。理解它们的数学形式、参数物理意义以及适用边界，是后续所有产量预测工作的前提。

三种递减方程的数学形式

Arps递减方程统一以产量随时间的变化率来描述递减行为。设 q 为时刻 t 的产量，q_i 为初始产量，D_i 为初始递减率，b 为递减指数，则三种形式可以写成如下统一表达式：

指数递减 (b = 0):   q = q_i · exp(-D_i · t)
双曲递减 (0 < b < 1): q = q_i / (1 + b · D_i · t)^(1/b)
调和递减 (b = 1):   q = q_i / (1 + D_i · t)

可以看到，调和递减其实是双曲递减在 b = 1 时的特例，而指数递减是 b = 0 时的极限情形。因此三者并非彼此独立，而是同一族曲线在不同递减指数下的表现。实际工作中，工程师通常先拟合双曲递减，再根据拟合出的 b 值判断数据更接近哪一种递减类型。

与产量方程配套的还有累计产量方程。以双曲递减为例，从 t = 0 积分到 t，可得累计产量：

N_p = [q_i^b / (D_i · (1 - b))] · (q_i^(1-b) - q^(1-b))

当 b = 0 时该式退化为指数递减的累计产量 N_p = (q_i - q) / D_i；当 b = 1 时则退化为调和递减的 N_p = (q_i / D_i) · ln(q_i / q)。累计产量方程在EUR估算中直接用于计算从当前产量递减到经济极限产量期间的总产出。

参数的物理意义

初始产量 q_i 是递减开始时刻的产量，通常取关井压裂后稳定生产的第一段产量，或拟合曲线在 t = 0 处的截距。它并不等于测试产量，而是经过拟合后代表递减起点的等效产量。

初始递减率 D_i 表示单位时间内产量下降的比例，量纲为时间的倒数（常用1/月或1/年）。例如 D_i = 0.08/月 意味着递减初期每月产量下降约8%。需要强调的是，D_i 是递减开始时的瞬时递减率，随着时间推移，双曲递减的瞬时递减率会逐渐减小。

递减指数 b 是三者中最关键也最难确定的参数。它控制曲线的弯曲程度：b = 0 时产量按指数快速衰减，曲线在半对数坐标上呈直线；b 越接近1，递减越缓慢，曲线越"平缓"，晚期产量维持能力越强。页岩气井由于基质致密、裂缝网络复杂，实际生产数据往往呈现 b 值偏高（常在0.7～1.2之间）的特征，这也是页岩气井"初期递减快、后期递减慢"的数学体现。

在双曲递减中，瞬时递减率 D(t) = D_i / (1 + b·D_i·t) 随时间单调下降；而指数递减的瞬时递减率恒等于 D_i。这一差别是判断数据属于哪类递减的重要依据。

适用条件与局限性

Arps方程本质上是经验公式，其成立依赖若干隐含假设：产量递减由地层能量衰竭主导，生产制度相对稳定，没有大规模补能或重复压裂等人工干预。对于页岩气井，这些条件在投产初期到中期通常基本满足，因此Arps方法在页岩气DCA中应用广泛。

但局限性同样明显。第一，b 值的外推风险：当 b ≥ 1 时，双曲递减的累计产量在时间趋于无穷时发散（数学上不收敛），这意味着直接外推会高估EUR。工程上通常对 b 设上限（如不超过1.0），或在晚期切换到指数递减。第二，Arps无法区分物理机理，它只描述"产量怎么变"，不解释"为什么这样变"，因此当出现井间干扰、压裂液返排、井筒积液或生产制度调整时，拟合曲线可能失真。第三，对短生产历史（不足12～18个月）的井，b 与 D_i 存在强相关性，拟合结果不唯一，外推EUR不确定性大。

实际拟合中的操作要点

在工程实践中，拟合Arps方程通常采用以下流程。首先对产量数据做时间对齐，剔除关井、修井、返排等异常段，保留稳定生产段。然后以月为时间单位，将产量与时间绘制在半对数坐标上：若数据点近似直线，优先尝试指数递减；若明显上凸（衰减变慢），则用双曲递减拟合。拟合时以累计产量或产量残差最小为目标，同时约束 b 在合理区间内。

举一个教学示例（数据为合理假设，非实测）：某页岩气井投产初期产量 q_i = 80×10⁴ m³/d，拟合得 D_i = 0.10/月，b = 0.85。若经济极限产量取 q_lim = 5×10⁴ m³/d，代入双曲递减产量方程反解时间：

q_lim / q_i = (1 + b·D_i·t)^(-1/b)
0.0625 = (1 + 0.085·t)^(-1.176)
解得 t ≈ 34 个月

再代入累计产量方程，可估算出从投产到经济极限的累计产出。若把 b 改为1.0重新计算，得到的EUR会明显偏大，这正是 b 值敏感性带来的不确定性。因此在实际EUR评估中，通常给出 b 的低、中、高三种情形，形成EUR区间而非单一数值。

不要用不足12个月的生产数据直接外推EUR。短历史下 b 与 D_i 高度相关，拟合曲线可以穿过数据点但外推结果差异巨大。建议至少积累18～24个月稳定生产数据后再做正式DCA，并配合其他方法交叉验证。

总结来说，Arps三方程提供了描述产量递减的统一数学框架：指数递减对应恒定递减率，调和递减对应递减率随产量线性下降，双曲递减介于两者之间并由 b 控制弯曲程度。掌握 q_i、D_i、b 的物理含义与相互约束，理解 b ≥ 1 时的外推风险和生产干预导致的拟合失真，才能在页岩气井EUR估算中正确使用这一经典工具，而不是把它当作"套公式出结果"的黑箱。