2. CO2驱油藏工程基础回顾

CO2驱方案评价的第一步不是直接套用商业模拟器,而是把油藏工程的基本方程重新“校准”到CO2驱的物理条件上。常规黑油模型所依赖的若干假设——岩石不可压缩、流体为两组分、相间瞬时平衡、达西渗流线性成立——在注CO2场景下会不同程度地失效。本节按“岩石—流体—相态—渗流方程”的顺序梳理这些基础,并说明每一项在CO2驱中需要如何修正、修正到什么程度。

一、储层岩石物性:孔隙度、渗透率与压缩性

孔隙度φ与渗透率K是储量计算和产能预测的基础参数。CO2驱中需要特别关注两点。其一是有效应力敏感性:注气压力升高会降低有效应力,使裂缝和孔隙通道张开,渗透率上升;而生产井附近压力下降又可能引起压实。对于低渗、裂缝性储层,渗透率变化可达数倍,若在模型中把K视为常数,历史拟合阶段往往出现“注气压力对不上”的系统偏差。其二是岩石压缩系数Cr:在CO2驱中,溶解作用与压力变化共同引起孔隙体积变化,物质平衡方程中的综合压缩系数应写成 Ct = Cr + Sw·Cw + So·Co + Sg·Cg,其中气相压缩系数Cg在高压下近似等于1/p,量级远大于Cr,不能忽略。

此外,CO2与地层水接触后可能溶蚀碳酸盐胶结物,长期注气会引起孔隙度和渗透率的结构性变化。方案评价阶段通常先用岩心驱替实验测定不同围压、不同CO2饱和压力下的K变化曲线,再以压力相关的乘子形式写入模型,而不是简单取一个平均值。

二、流体物性:黏度、密度与界面张力

CO2驱最关键的流体物性差异是黏度比。以示例条件(温度80 ℃、压力20 MPa、原油黏度5 mPa·s)估算,CO2在该条件下的黏度约为0.06~0.08 mPa·s,黏度比接近60~80。高黏度比直接导致黏性指进和波及效率下降,这是CO2驱区别于水驱的核心矛盾。密度方面,CO2在超临界区密度可达0.6~0.8 g/cm³,与轻质原油接近,重力分异作用相对减弱,但在气顶或高渗条带中仍会形成上浮通道。

界面张力(IFT)随压力升高而下降,在最小混相压力(MMP)附近趋于零,这是混相驱替的物理基础。方案评价中必须明确:目标油藏压力是否高于MMP。若低于MMP,则属于非混相或近混相驱,采收率机理以溶胀、降黏和部分组分交换为主,不能按混相驱的驱替效率预期来设计注采比。

三、相态基础:组分模型与状态方程

CO2驱必须使用组分模型而非黑油模型。原因是CO2与原油之间发生多组分传质:CO2溶于油使油相体积膨胀(溶胀系数可达1.1~1.4),同时抽提原油中的轻质组分,形成随压力、组分变化的相行为。工程上通常采用Peng-Robinson状态方程描述气液相平衡,并通过实验(PVT分析、膨胀实验、多次接触实验)回归临界参数和二元交互系数。

相态描述的关键输出包括:饱和压力、溶胀因子、CO2在油相中的溶解度Rs、油气两相的相对体积和黏度。这些参数直接进入渗流方程中的相对渗透率和毛管压力曲线。需要提醒的是,状态方程对重质组分和沥青质沉淀的预测能力有限,若原油含蜡或沥青质较高,应补充沥青质沉淀实验,否则模型会高估长期注入能力。

四、达西定律与多相渗流方程的适用性

达西定律的原始形式 q = -(K·A/μ)·(dp/dx) 在CO2驱中仍作为基本骨架使用,但需注意三个适用性条件。第一,高速非达西效应:注气井近井地带气体流速高,惯性力不可忽略,需引入Forchheimer修正,即在压力梯度项中加入βρv²项,否则会高估注气能力。第二,多相渗流的相对渗透率:CO2驱中油气相对渗透率曲线受黏度比、界面张力和润湿性影响显著,低IFT条件下相对渗透率曲线趋于直线,需要在模型中采用随IFT变化的相对渗透率模型,而不是固定曲线。第三,毛管压力:在混相条件下毛管压力趋于零,但在非混相区仍是控制垂向和横向窜流的重要因素。

多相渗流方程(以油、气两相为例)可写为:

∂/∂t (φ·So·ρo) + ∇·(ρo·vo) = -qo
∂/∂t (φ·Sg·ρg) + ∇·(ρg·vg) = -qg
vo = -(K·Kro/μo)·(∇po - ρo·g·∇z)
vg = -(K·Krg/μg)·(∇pg - ρg·g·∇z)
pc = pg - po,  So + Sg = 1

在CO2驱中,上述方程还需补充组分守恒方程,因为CO2在油相和水相中均有溶解。工程上常采用“黑油+溶解气”的扩展形式或完整组分模型。方案评价阶段应至少验证:物质平衡闭合误差小于1%、注气量与产气量之差可由溶解和滞留解释、压力传播与实测井底压力趋势一致。

教学示例:某低渗油藏温度85 ℃、原始压力22 MPa、原油黏度3.5 mPa·s、MMP实测约24 MPa。由于原始压力低于MMP,方案初期应按非混相驱设计,注气压力需提升至24 MPa以上才能进入混相区。这一判断直接决定了注采井距和注入压力的取值,是方案评价的起点。
CO2驱油藏工程基础参数链与方程适用性 储层岩石物性 φ、K、Cr 有效应力敏感 K随压力变化 → 用K(p)乘子 流体物性 μo、μg、ρ、IFT 黏度比 60~80 溶胀 1.1~1.4 → 控制指进与波及 相态基础 PR状态方程 MMP、Rs、相组成 混相/非混相判别 → 组分模型必需 渗流方程 达西 + 多相 Forchheimer修正 Krg随IFT变化 → 组分守恒补充 方程适用性检查与修正路径 高速非达西 注气井近井 v 高 加 βρv² 惯性项 相对渗透率 低IFT曲线趋直 用IFT相关Krg模型 组分守恒 CO2溶于油与水 补组分方程 方案评价验证指标 物质平衡误差 < 1% | 注采气量差可由溶解与滞留解释 | 井底压力趋势与实测一致 | 压力高于MMP才按混相驱设计

五、从基础参数到方案评价的衔接

把上述基础串起来,方案评价阶段应形成一条可核查的参数链:岩心实验给出K(p)和相对渗透率曲线,PVT实验给出MMP、溶胀因子和溶解度,状态方程回归后进入组分模型,模型中的渗流方程按高速非达西和IFT相关相渗进行修正,最后用物质平衡和压力历史做闭合检验。任何一环缺失,都会使采收率预测偏离实际。例如,若忽略黏度比导致的指进,模型会高估波及体积;若忽略有效应力对K的影响,注气压力预测会系统性偏低。

因此,CO2驱的油藏工程基础回顾不是复习课,而是一次“适用性审计”:逐项确认每个基础参数和方程在当前油藏条件下是否仍然成立,不成立的部分就是方案评价中必须修正或补充实验的地方。这一审计结论,将直接决定后续注采参数设计和经济评价的可靠性边界。