4、评价与储量计算:储量分类标准(SPE/AAPG/WPC)、容积法、物质平衡法、不确定性分析
储量评价是油气田开发价值链中承上启下的核心环节——它既是对勘探发现的量化确认,又是开发方案编制、投资决策和经济评价的基础。本章围绕储量计算的三根支柱展开:分类标准解决“怎么定义储量”、计算方法解决“怎么算准储量”、不确定性分析解决“算出来的储量有多可信”。理解这三者的逻辑关系,是开发地质工程师和油藏工程师的基本功。
4.1 储量分类标准:从“地质认识”到“商业语言”
核心原理
储量分类的本质是将地下油气资源按地质确定性和经济可行性两个维度进行分级。目前国际通行的SPE/AAPG/WPC分类体系(2007年修订版,2018年补充)将资源分为三大类:储量(Reserves)、概算资源(Contingent Resources)和远景资源(Prospective Resources)。其中储量又细分为证实储量(Proved, 1P)、概算储量(Probable, 2P)和可能储量(Possible, 3P)。
分类的核心逻辑是:商业成熟度决定是否属于储量,不确定性程度决定属于哪个子类。一个项目必须同时满足“已发现”、“可开发”、“经济可行”三个条件才能计入储量。SEC(美国证券交易委员会)的规则更严格,只承认“证实储量”,且要求基于实际生产或可靠测试数据。
分类标准对比表
| 分类维度 | 证实储量(1P) | 概算储量(2P) | 可能储量(3P) |
|---|---|---|---|
| 确定性水平 | ≥90%概率(确定性法) | 1P+2P ≥50%概率 | 1P+2P+3P ≥10%概率 |
| 数据基础 | 生产数据、DST测试、岩心/测井 | 测井解释+类比油藏 | 地质类比+地球物理间接证据 |
| 经济条件 | 当前价格与成本下可开发 | 合理预期价格下可开发 | 乐观价格情景下可开发 |
| SEC认可 | 是(唯一认可) | 否 | 否 |
案例拆解:某海上油田储量申报
某渤海湾海上油田,含油面积12.5 km²,平均有效厚度18m,孔隙度28%,含油饱和度65%,原油体积系数1.15。开发方案设计水平井+注水。
- 初始评估:基于3口探井+2口评价井,采用确定性法计算,证实储量(1P)为 5200×10⁴t,概算储量(2P)为 7800×10⁴t。
- SEC审核问题:由于注水开发尚未投产,SEC要求将1P储量下调至仅基于一次采油(溶解气驱+弹性驱)的 3800×10⁴t。
- 应对策略:加快注水先导试验,取得注水受效数据后重新申报,最终1P储量恢复至 4900×10⁴t。
这个案例说明:分类标准不是纯技术问题,而是技术+商业+监管的复合决策。开发节奏会影响储量申报口径,进而影响公司资产估值和融资能力。
常见误区
- 误区1:把“地质储量”当“可采储量”上报——地质储量是地下总量,可采储量要考虑采收率。
- 误区2:忽略经济边界——油价下跌时,高成本储量应降级为概算资源。
- 误区3:概率法与确定性法混用——同一项目不能部分参数用概率法、部分用确定性法。
4.2 容积法:静态储量的基石
核心原理
容积法通过计算油藏孔隙体积中烃类所占的体积来估算原始地质储量(OOIP),公式为:
OOIP = 100 × A × h × φ × (1-Swi) × ρo / Boi
其中:A为含油面积(km²),h为有效厚度(m),φ为有效孔隙度(小数),Swi为束缚水饱和度(小数),ρo为地面原油密度(t/m³),Boi为原始原油体积系数(m³/m³)。对于气藏,公式变为:
OGIP = 0.01 × A × h × φ × (1-Swi) × (Tsc×Pi) / (Psc×T×Zi)
容积法的核心假设是储层参数在控制体积内具有代表性。因此,参数取值方法决定了计算精度:
- 含油面积:由构造图+油水界面确定,受断层、尖灭、边底水影响。
- 有效厚度:通过测井解释的油层下限标准(孔隙度、饱和度截止值)确定。
- 孔隙度/饱和度:取岩心分析、测井解释的算术平均或厚度加权平均。
- 体积系数:由PVT实验获得,或通过相关经验公式估算。
操作流程
- 建立油藏地质模型:构造层面、断层模型、储层展布。
- 确定油水界面:结合压力梯度、测井解释、试油数据。
- 圈定含油面积:在构造图上勾绘油水界面等高线。
- 分层计算:按小层或流动单元分别计算,避免非均质性被平均化。
- 参数加权:采用面积加权或体积加权,而非简单算术平均。
- 汇总校验:与物质平衡法、递减法结果交叉验证。
案例拆解:某陆相砂岩油藏容积法计算
某断块油藏,含油面积3.8 km²,有效厚度12.5m,平均孔隙度24%,束缚水饱和度30%,原油密度0.86 t/m³,体积系数1.12。计算如下:
OOIP = 100 × 3.8 × 12.5 × 0.24 × (1-0.30) × 0.86 / 1.12
= 100 × 3.8 × 12.5 × 0.24 × 0.70 × 0.86 / 1.12
= 100 × 3.8 × 12.5 × 0.24 × 0.70 × 0.768
= 613.2 × 10⁴ t
该区块后期通过数值模拟验证,实际动态储量约580×10⁴t,误差约5.7%,在可接受范围内。误差主要来源于含油面积边界的不确定性(西北侧断层位置存在±200m摆动)。
常见误区
- 误区1:用全油藏平均参数代替分层参数——厚层与薄层物性差异大时,平均化会掩盖高渗层贡献。
- 误区2:忽略有效厚度截止值——将泥质砂岩、钙质夹层计入有效厚度会虚增储量。
- 误区3:不区分地面与地下体积——体积系数换算错误是新手常见问题。
4.3 物质平衡法:动态验证的利器
核心原理
物质平衡法基于质量守恒定律:原始烃类体积 = 剩余烃类体积 + 采出烃类体积 + 侵入水体积。对于油藏,其通式为:
N = [Np×Bo + (Gp-Ginj)×Bg - Wp×Bw + Winj×Bw] / [Boi×(Bo-Boi) + Bg×(Rsi-Rs) + m×Boi×(Bg/Bgi-1) + (1+m)×Boi×Ce×ΔP]
其中N为原始地质储量,Np为累计产油,Gp为累计产气,Wp为累计产水,Winj为累计注水,m为气顶指数,Ce为有效压缩系数。
物质平衡法的优势在于利用动态生产数据反推储量,不依赖地质参数的精确取值。但它的前提是油藏压力系统统一(即各井点压力连通),且需要准确的PVT数据和压力历史。
方法对比
| 对比维度 | 容积法 | 物质平衡法 |
|---|---|---|
| 数据来源 | 静态地质参数 | 动态生产+压力数据 |
| 适用阶段 | 勘探评价期、开发早期 | 开发中后期(有足够生产史) |
| 对非均质性 | 敏感(需精细建模) | 不敏感(体积平均效应) |
| 对压力数据 | 不依赖 | 高度依赖(需定期测压) |
| 结果含义 | 原始地质储量 | 动态连通储量(可能小于地质储量) |
案例拆解:边水油藏物质平衡诊断
某边水油藏投产5年,累计产油85×10⁴t,累计产水12×10⁴m³,地层压力从25.8 MPa降至22.3 MPa。容积法计算地质储量为1200×10⁴t。但物质平衡法计算结果显示:
- 不考虑水侵:计算储量仅680×10⁴t,明显低于容积法。
- 考虑稳态水侵:计算储量1150×10⁴t,与容积法接近。
- 诊断结论:存在活跃边水补给,压力下降速度低于封闭油藏预期。
该案例说明:物质平衡法不仅能算储量,还能诊断油藏驱动类型。若计算储量远低于容积法,往往提示存在未识别的水侵或能量补给;若远高于容积法,则可能提示含油面积比预期大或存在未动用的连通储集体。
常见误区
- 误区1:忽略压力数据质量——压力计漂移、测压时机不一致会导致计算偏差。
- 误区2:将物质平衡法用于强非均质油藏——层间渗透率差异大时,各层压力不统一,物质平衡假设失效。
- 误区3:不区分溶解气驱与气顶驱——气顶指数m取值错误会显著影响结果。
4.4 不确定性分析:从“一个数”到“一个分布”
核心原理
储量计算中的每个输入参数都有不确定性,传统确定性法只给一个“最佳估计值”,无法回答“储量可能低到什么程度”或“高到什么程度”。不确定性分析通过概率分布描述每个参数的不确定性,然后通过蒙特卡洛模拟传播这些不确定性,最终得到储量的概率分布(P10、P50、P90)。
关键参数通常采用三角分布或正态分布:
- 含油面积:三角分布(最小值/最可能/最大值),受构造解释和油水界面不确定性影响。
- 有效厚度:正态分布或对数正态分布,受测井解释截止值影响。
- 孔隙度:正态分布,由岩心分析统计得到。
- 含油饱和度:正态分布,受电阻率解释模型影响。
- 体积系数:窄分布,PVT实验误差较小。
操作流程
- 参数分布定义:基于岩心、测井、试井、类比油藏确定每个参数的P10/P50/P90值。
- 相关性设定:孔隙度与饱和度通常正相关(高孔往往高饱),需设定相关系数。
- 蒙特卡洛模拟:随机抽样5000-10000次,每次计算一个储量值。
- 结果统计:输出储量累积概率分布曲线,读取P10(低值)、P50(中值)、P90(高值)。
- 敏感性分析:计算各参数对储量方差的贡献率,识别主要风险源。
案例拆解:某致密气藏储量不确定性评估
某致密砂岩气藏,含气面积45 km²,有效厚度20-35m,孔隙度6%-10%,含气饱和度55%-70%。由于储层强非均质,各参数不确定性大:
| 参数 | P10 | P50 | P90 | 分布类型 |
|---|---|---|---|---|
| 含气面积(km²) | 38 | 45 | 52 | 三角分布 |
| 有效厚度(m) | 22 | 28 | 34 | 正态分布 |
| 孔隙度(%) | 6.5 | 8.0 | 9.5 | 正态分布 |
| 含气饱和度(%) | 58 | 63 | 68 | 正态分布 |
蒙特卡洛模拟5000次后,OGIP的P10=280×10⁸m³,P50=410×10⁸m³,P90=560×10⁸m³。敏感性分析显示:有效厚度贡献了52%的方差,孔隙度贡献28%。因此,后续评价应优先加密取心井,降低厚度解释不确定性。
常见误区
- 误区1:参数分布范围拍脑袋——必须有数据支撑或类比依据。
- 误区2:忽略参数间相关性——孔隙度与饱和度强相关时,独立抽样会夸大不确定性。
- 误区3:把P50当“最可能值”——P50是中位数,不是众数;对于偏态分布,两者差异明显。
4.5 三种方法协同与综合决策
实际工作中,三种方法不是互相替代,而是互相校验、互为补充。推荐的工作流程为:
- 勘探评价期:以容积法为主,建立初步储量范围,配合不确定性分析给出P10/P50/P90。
- 开发早期:容积法+物质平衡法交叉验证,若两者差异>20%,需排查地质模型或压力数据问题。
- 开发中后期:以物质平衡法和递减曲线法为主,容积法结果作为上限参考。
- 年度储量复核:结合新钻井、新生产数据,更新参数分布,重新计算储量。
本章内容覆盖了储量评价的完整链条:从分类标准建立“商业语言”,到容积法提供“静态骨架”,再到物质平衡法注入“动态血液”,最后用不确定性分析给出“风险边界”。四者结合,才能形成对油气田开发价值链条中这一关键节点的可靠量化认知。