2. 油藏物质平衡方程推导与物理意义

物质平衡法的全部推导都建立在一个朴素而严格的物理事实上：在油藏开采过程中，储层孔隙体积的变化量，必然等于其中流体（油、气、水）体积变化量之和。换句话说，地层是一个封闭的"容器"，采出的流体由三部分来源补充——原油本身的膨胀、溶解气脱出的膨胀、以及束缚水与岩石的膨胀，再加上可能侵入的外来水。只要把这几项在地下条件下的体积变化逐一写清楚，并令其与采出流体的地下体积相等，就得到了通用物质平衡方程。

一、符号约定与参照条件

推导前必须先固定参照系，否则正负号会混乱。本节统一采用如下约定：以原始油藏压力 p_i 下的状态为基准，所有体积均折算到地下条件；凡是从油藏中"采出"的量取正号，凡是"侵入"或"膨胀补充"的量在方程中移到等号另一侧时取负号。常用符号为：N 为原始地质储量（地面原油体积，m³）；N_p 为累产油量（地面体积）；R_p 为累计生产气油比；R_si 为原始溶解气油比；R_s 为当前压力下溶解气油比；B_o、B_oi 为当前与原始原油体积系数；B_g 为天然气体积系数；W_e 为累计水侵量；W_p 为累产水量；c_w、c_f 为地层水与岩石压缩系数；S_wi 为束缚水饱和度；Δp = p_i − p 为压降。

二、从体积平衡到通用方程

设油藏原始孔隙体积为 V_p，其中原油的地下体积为 N·B_oi。当压力由 p_i 降到 p 时，油藏内剩余流体占据的地下体积应等于原始孔隙体积加上岩石与束缚水的膨胀量，再减去水侵占据的部分。把"采出"与"地下剩余"两侧对齐，可写出通用形式：

N·B<sub>oi</sub> + N·(R<sub>si</sub> − R<sub>s</sub>)·B<sub>g</sub> + W<sub>e</sub>
   = N·B<sub>o</sub> + (R<sub>p</sub> − R<sub>s</sub>)·N<sub>p</sub>·B<sub>g</sub> + W<sub>p</sub>·B<sub>w</sub>
     + N·B<sub>oi</sub>·(c<sub>w</sub>·S<sub>wi</sub> + c<sub>f</sub>)·Δp / (1 − S<sub>wi</sub>)

整理后即得教材中常见的通用物质平衡方程：

N·(B<sub>o</sub> − B<sub>oi</sub>) + N·(R<sub>si</sub> − R<sub>s</sub>)·B<sub>g</sub>
  + N·B<sub>oi</sub>·(c<sub>w</sub>·S<sub>wi</sub> + c<sub>f</sub>)·Δp/(1 − S<sub>wi</sub>)
  = N<sub>p</sub>·[B<sub>o</sub> + (R<sub>p</sub> − R<sub>s</sub>)·B<sub>g</sub>] + W<sub>p</sub>·B<sub>w</sub> − W<sub>e</sub>

等号左侧三项依次是：原油膨胀项、溶解气脱出膨胀项、岩石与束缚水膨胀项；等号右侧是采出物的地下体积（油、伴生气、水）减去水侵补充量。这就是"地下体积平衡"的完整表达。

三、各体积项的物理含义

原油膨胀项 N(B_o − B_oi)：压力下降使原油体积系数增大，地下原油"胀大"，这部分胀出的体积正是驱动原油流向井底的弹性能量之一。它只与压降和原油压缩性有关，与生产时间无关。

溶解气脱出项 N(R_si − R_s)B_g：压力低于饱和压力后，原本溶在油中的天然气析出，在地层中形成自由气。由于气体体积系数 B_g 远大于 1（地下体积远大于地面体积），这一项往往是弹性能量中贡献最大的部分，也是溶解气驱油藏的主要驱动来源。

岩石与束缚水膨胀项：孔隙压力降低使岩石骨架被上覆压力压缩、孔隙体积缩小，同时束缚水膨胀。两者合起来用综合压缩系数 c_t = (c_wS_wi + c_f)/(1 − S_wi) 表示。该项数值通常只占总量百分之几，但在高含水、低饱和压力或封闭小油藏中不可忽略。

水侵项 W_e：边底水在压差作用下侵入油藏，相当于从外部补充了地下体积，因此移到方程右侧时取负号。它是水驱油藏物质平衡计算中最难确定的量，通常需要结合水侵模型（如 Fetkovich、Van Everdingen-Hurst）拟合。

采出项：N_p[B_o + (R_p − R_s)B_g] 表示累产油及其携带的伴生气折算到地下的体积；W_pB_w 为产出水的地下体积。注意 R_p 是累计生产气油比，它包含了溶解气与自由气两部分，这正是"自由气"在方程中的体现方式。

四、一个教学示例

假设某未饱和油藏：N = 2.0×10⁶ m³，B_oi = 1.20，B_o = 1.24，R_si = 80 m³/m³，当前 R_s = 80（未饱和，无气脱出），c_w = 4.5×10⁻⁴ MPa⁻¹，c_f = 6.0×10⁻⁴ MPa⁻¹，S_wi = 0.25，Δp = 8 MPa，W_e = 0，W_p = 0。则原油膨胀项 = 2.0×10⁶×(1.24−1.20) = 8.0×10⁴ m³；岩石与水膨胀项 = 2.0×10⁶×1.20×(4.5×10⁻⁴×0.25 + 6.0×10⁻⁴)×8/(1−0.25) ≈ 1.75×10⁴ m³。两项合计约 9.75×10⁴ m³，即在该压降下理论可采出约 9.75×10⁴ m³ 地下体积的油，对应地面累产油 N_p ≈ 9.75×10⁴/1.24 ≈ 7.86×10⁴ m³。这一结果说明：未饱和油藏的弹性能量主要来自原油膨胀，岩石与水膨胀只占约 18%。

五、使用方程时的判断要点

• 压力高于饱和压力时，R_s = R_si，溶解气项为零，方程退化为弹性驱动形式。
• 存在气顶时，需在方程中单独加入气顶膨胀项 m·N·B_oi·(B_g/B_gi − 1)，本节通用式未含气顶，属简化形式。
• 水侵量 W_e 未知时，方程含 N 与 W_e 两个未知数，需用多压力点数据联立或引入水侵模型求解。
• 所有体积系数必须取自同一压力、同一温度条件下的 PVT 数据，避免不同来源数据混用。

掌握这一方程的推导逻辑后，后续无论是弹性驱、溶解气驱还是水驱的储量计算，都只是在本通用式上做取舍与变形。理解每一项"从哪来、到哪去"，比记住公式本身更重要。