3. 容积法计算公式解析

容积法(Volumetric Method)是油气资源评价中计算圈闭原始地质储量最基础、应用最广的方法。它的核心思想非常朴素:先算出圈闭中能够储存油气的岩石总体积,再乘以岩石的孔隙度得到孔隙体积,然后乘以含油饱和度得到含油孔隙体积,最后把地下条件下的油气体积折算到地面标准条件。整个过程可以浓缩为一个公式,但公式中每一个参数的物理含义、取值依据和单位约定,才是实际评价工作中真正决定结果可靠性的部分。

原始地质储量公式的完整形式

对于油藏,容积法原始地质储量(地面标准条件下的原油体积)的常用表达式为:

N = 100 · A · h · φ · So / Boi

N   —— 原始地质储量(地质储量),10^4 m³(万吨需再乘地面原油密度)
A   —— 含油面积,km²
h   —— 平均有效厚度,m
φ   —— 平均有效孔隙度,小数
So  —— 平均原始含油饱和度,小数
Boi —— 原始原油体积系数,m³/m³(地下体积/地面体积,无量纲比值)
100 —— 单位换算系数(km²·m → 10^4 m³)

对于气藏,对应的表达式为:

G = 0.01 · A · h · φ · Sg / Bgi

G   —— 原始天然气地质储量,10^8 m³
Sg  —— 平均原始含气饱和度,小数
Bgi —— 原始天然气体积系数,m³/m³
0.01 —— 单位换算系数(km²·m → 10^8 m³)

两个公式结构完全一致,差别只在储量单位约定和换算系数。理解这一点,比死记两个系数更重要:换算系数的本质是把“面积用 km²、厚度用 m”算出的体积(km²·m = 10⁶ m³)转换成行业习惯的 10⁴ m³ 或 10⁸ m³ 计量单位。

容积法储量计算:参数链与量纲传递 含油面积 A km²(平面投影) 有效厚度 h m(垂直厚度) 有效孔隙度 φ 小数(0~1) 含油饱和度 So 小数(0~1) 体积系数 Boi m³/m³(>1) 岩石总体积 A · h → km²·m 孔隙体积 A · h · φ 地下含油体积 A · h · φ · So 地面原油体积 ÷ Boi(收缩折算) 单位换算 ×100 10⁶ m³ → 10⁴ m³ 原始地质储量 N N = 100 · A · h · φ · So / Boi 单位:10⁴ m³ 量纲:km²·m·(小数)³ ÷ (m³/m³) × 系数 量纲校核要点 ① A 用 km²,h 用 m,不可混用 m² ② φ、So 必须为小数,30% 写作 0.30 ③ Boi 为地下/地面体积比,恒大于 1 ④ 系数 100 与 0.01 对应不同储量单位

五个参数的物理含义与取值逻辑

含油面积 A 指圈闭中油气聚集的平面投影面积,单位 km²。它由构造图上的油水界面(或气水界面)与圈闭闭合等高线共同圈定,不是圈闭的闭合面积。实际取值时通常采用圈闭高点与油水界面所限定的含油边界在构造图上的投影面积,对于断层遮挡圈闭还要扣除断层带影响。面积是储量公式中对结果影响最敏感的参数之一,边界每外扩一点,储量按面积比例线性放大。

有效厚度 h 指储层中具有储集能力、能被油气饱和并参与渗流的岩层厚度,单位 m。它不等于地层厚度,需要依据岩性、物性(孔隙度、渗透率)和含油性下限逐井逐层剔除泥质夹层、致密段和水层。实际工作中常用测井解释成果配合岩心标定,先确定有效厚度下限标准,再在单井上统计有效厚度,最后用算术平均或面积权衡平均得到评价单元的平均有效厚度。

有效孔隙度 φ 指岩样中连通孔隙体积与岩石总体积之比,用小数表示。它只计入相互连通、能容纳并允许流体流动的孔隙,不包括孤立孔洞。取值来源包括岩心分析、测井解释和地震反演,通常以岩心分析为基准标定测井解释,再按厚度加权平均。

含油饱和度 So 指孔隙体积中被原油占据的比例,用小数表示。它与束缚水饱和度互补:So = 1 − Swi。原始含油饱和度应取油藏原始状态下的值,不能使用开发后期的剩余油饱和度。常用测井解释、岩心分析或毛管压力资料确定。

体积系数 Boi 是原油在地下油藏条件下的体积与地面脱气后体积之比,单位为 m³/m³,是一个大于 1 的无量纲比值。原油从地下到地面会因压力降低、溶解气逸出而收缩,因此地下体积大于地面体积,公式中必须除以 Boi 才能得到地面储量。Boi 一般由高压物性(PVT)分析或经验相关式确定,轻质油、挥发油数值较高,重油接近 1.0~1.1。

单位换算与量纲分析

量纲分析是检查公式使用是否正确的最有效手段。以油藏公式为例,把各参数单位代入:

A[km²] · h[m] · φ[小数] · So[小数] ÷ Boi[m³/m³]
= km²·m = 10⁶ m³(地下孔隙中油的体积,折算到地面后仍为 10⁶ m³ 量级)
× 100 → 10⁸ m³? 不对,应为 10⁶ m³ = 100 × 10⁴ m³

关键在最后一步:1 km²·m = 10⁶ m³,而 10⁶ m³ = 100 × 10⁴ m³,所以乘以 100 后结果单位就是 10⁴ m³。同理,气藏公式中 1 km²·m = 10⁶ m³ = 0.01 × 10⁸ m³,故系数取 0.01,结果单位为 10⁸ m³。这就是两个系数来源的全部逻辑——它们不是经验数字,而是纯粹的 unit conversion。

做量纲校核时,建议先把所有参数统一到“km²、m、小数、m³/m³”这套单位,再套系数。若面积误用 m²、厚度误用 cm,或孔隙度误填百分数(如填 25 而非 0.25),结果会偏差几个数量级,且往往在与其他方法对比时才被发现。

一个可复算的教学示例

假设某砂岩圈闭评价单元参数为(示例假设值,非实测):含油面积 A = 3.2 km²,平均有效厚度 h = 12.5 m,平均有效孔隙度 φ = 0.22,平均原始含油饱和度 So = 0.65,原始原油体积系数 Boi = 1.15。代入公式:

N = 100 × 3.2 × 12.5 × 0.22 × 0.65 ÷ 1.15
  = 100 × 3.2 × 12.5 × 0.143 ÷ 1.15
  = 100 × 40 × 0.143 ÷ 1.15
  = 572 ÷ 1.15
  ≈ 497.4(10⁴ m³)

即该圈闭原始地质储量约 497×10⁴ m³。若需要以质量(万吨)表示,再乘以地面原油密度(如 0.85 t/m³),得约 423 万吨。这个示例的价值不在于数字本身,而在于演示计算顺序:先算 A·h 得岩石体积,再乘 φ 得孔隙体积,再乘 So 得地下含油体积,最后除以 Boi 并乘系数。每一步都可以单独检查数量级是否合理,例如 A·h = 40 km²·m 对应 4×10⁷ m³ 岩石体积,乘以 0.22×0.65 后地下含油体积约 5.7×10⁶ m³,除以 1.15 后地面体积约 5.0×10⁶ m³,与 497×10⁴ m³ 一致。

参数敏感性:哪个参数最值得花力气

从公式结构看,N 与 A、h、φ、So 均成正比,与 Boi 成反比。因此各参数的相对误差会以相同比例传递到储量结果:面积误差 10%,储量就误差约 10%;孔隙度误差 10%,储量同样误差约 10%。但在实际评价中,各参数的获取难度和不确定性差异很大。面积受圈闭落实程度和油水界面精度控制,往往是最大的不确定性来源;有效厚度受井控程度和下限标准影响,井少时误差可达 20%~30%;孔隙度和饱和度在取心井多的区块相对可控;体积系数由 PVT 样品代表性决定。因此,提高储量精度最有效的途径通常是加强构造解释与油水界面论证、增加评价井控制程度,而不是反复微调孔隙度小数位。

常见错误有三类:一是把圈闭闭合面积直接当作含油面积;二是有效厚度中混入泥质夹层或水层;三是体积系数方向用反(乘以 Boi 而非除以 Boi),导致储量偏大。第三类错误在手工计算中尤其常见,建议每次计算后用“地下体积应大于地面体积”这一常识做一次方向校验。