📝 摘要
提高采收率（EOR）被认为是提高常规和非常规油藏最终采收率的关键技术。注气是应用最广泛的EOR方法之一，其中超临界CO₂因其易实现混相、与储层流体配伍性好、密度和粘度较高以及碳封存需求而备受青睐。然而，所有气体注入方法都不可避免地面临指进、窜流和重力超覆等挑战。
将注入的气相发泡是改善波及效率的成熟技术。泡沫是通过将气体与表面活性剂溶液共同注入而形成的。在现场应用泡沫之前，通常会进行岩心驱替实验以了解泡沫特性和流动特征，然后通过建模研究确定模型参数及泡沫降低气相活动性的能力。在储层尺度上模拟泡沫流变性具有挑战性，因为模型必须能够复现实验室小尺度数据，并预测泡沫在大油田尺度上的传播行为。本研究使用CMG STARS模拟器，研究了超临界CO₂泡沫提高采收率的多个方面，从机理泡沫模型与实验室岩心驱替数据的拟合，到泡沫EOR应用的现场尺度评价。
主要结论：与气相CO₂、N₂、烃气、烟道气等其他气体泡沫相比，超临界CO₂泡沫具有极低的启动压力梯度，可在低压力梯度环境下形成稳定泡沫并将其更深地注入储层。机理泡沫模型可为现场泡沫EOR设计提供指导，例如与CMG STARS模拟结合，可确定最佳的注入泡沫质量和注入速率。
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | CMG STARS（先进过程与热采模拟器） |
| 模型类型 | 等温、三相（水相、油相、气相）、多组分组分模拟器。 |
| 模拟对象 | 超临界CO₂泡沫提高采收率（EOR），包括： • 机理模型拟合：使用STARS（隐含纹理泡沫模型）进行重力超覆模拟（第3章） • 油田尺度评价：哥伦比亚Lisama油田Mugrosa组（五点井网）的CO₂泡沫EOR优化（第4章） |
| 应用方式 | • 历史拟合/模型校准：利用岩心驱替实验数据，通过机理泡沫模型（气泡种群平衡模型）确定泡沫参数（MRF、S_w*等） • 敏感性分析：使用STARS模拟泡沫注入条件下的重力超覆距离 • 油田尺度优化：使用STARS对五点井网区块进行模拟，结合机理模型确定的MRF值，绘制波及效率等值线图，确定最优注入参数 |
| 泡沫模拟方法 | • 机理模型（第2-3章）：基于气泡种群平衡，追踪泡沫生成、 coalescence、气相 trapped、非牛顿流变等 • 隐含纹理泡沫模型（第4章）：现场尺度模拟采用STARS的泡沫模块，通过指定MRF（STARS中称FM）来降低气相相对渗透率：k_rgᶠ = k_rg(S_w) × MRF⁻¹ |
| 主要研究参数 | 启动压力梯度(∇P₀)、极限水饱和度(S_w*)、活动性降低因子(MRF)、注入总速率(Q_t)、注入泡沫质量(f_g)、地层非均质性（页岩渗透率） |
| 几何模型 | • 第3章：假设的均质圆柱形油藏（半径1000 ft，厚度275 ft），中心一口注入井 • 第4章：哥伦比亚Lisama油田Mugrosa组区块，五点井网（1注4采），尺寸：30×30×31网格，单网格66.7 ft×66.7 ft×7-11 ft |
| 储层参数（第3章） | • 渗透率：450 mD • 孔隙度：22% • 厚度：275 ft • 温度：110 °F • 压力：1555 psia • 水/CO₂密度：62.6/36.5 lb/ft³ • 水/CO₂粘度：0.65/0.07 cp |
| 储层参数（第4章） | • 压力：2500 psia • 温度：140-160 °F • 初始含油饱和度：45% • 平均孔隙度：18.5% • 渗透率范围：20-200 mD（多层，A/B/C/D） • 油粘度：1.7 cp @ 2500 psia |
| 模拟时长 | 第3章：4000天（达到稳态）；第4章：20年 |
| 评价指标 | 波及效率、累计产油量、泡沫传播距离（R_csw， R_gs）、生产动态 |
文中明确指出：
“The field-scale simulation in this study applies CMG STARS module by incorporating MRF…” （第4.3.2节）
“CMG STARS simulation and the Stone and Jenkins model confirm that gravity segregation also limits foam propagation distance.” （第3.5节）
致谢：”I would like to acknowledge Computer Modeling Group (CMG) for generous donation of CMG-STARS for this research.”
🧪 模拟方案与主要结果
1. 机理泡沫模型拟合与传播距离预测（第2-3章）
-
模型校准：利用Yin (2007) 的超临界CO₂泡沫岩心驱替数据，校准了气泡种群平衡模型的参数（表2.4-2.6），成功匹配了“S”形曲线（三个泡沫状态）和两个强泡沫流态（高质量区/低质量区）（图2.8-2.10）。
-
启动压力梯度的影响：超临界CO₂泡沫的∇P₀（~1 psi/ft）远低于N₂泡沫（~30 psi/ft）。低∇P₀使得泡沫从弱泡沫到强泡沫的过渡平滑且单调，无中间态，因此能传播得更远（图3.3）。
-
传播距离受限于两种机制：
-
强泡沫→弱泡沫转换 (R_csw)：湿泡沫（f_g=60%）比干泡沫（f_g=90%）传播更远。低∇P₀ (1 psi/ft) 时，R_csw可达数百至数千英尺（图3.6-3.10）。
-
重力超覆 (R_gs)：通过CMG STARS模拟和Stone-Jenkins模型确认。结果与机理模型预测的MRF值一致（图3.11-3.14）。
-
组合机制：低∇P₀时，传播主要受重力超覆限制；高∇P₀时，主要受强→弱泡沫转换限制（图3.22-3.23）。
-
2. 油田尺度优化（第4章）- Lisama油田案例
-
MRF-f_g路径：利用机理模型生成不同注入速率(Q_t)下的MRF与注入泡沫质量(f_g)的关系曲线（图4.15）。
-
波及效率等值线图：通过CMG STARS对132个场景（3种Q_t × 11种f_g × 4种MRF）进行20年模拟，绘制了波及效率随MRF和f_g变化的等值线图（图4.16-4.18）。
-
最优条件确定：将MRF-f_g路径叠加在波及效率等值线图上，可直观确定最优f_g。结果表明，最优f_g随注入速率(Q_t)增加而增加（图4.19）：
-
低Q_t (23,358 ft³/d)：最优f_g ≈ 50%
-
中Q_t (46,717 ft³/d)：最优f_g ≈ 55%
-
高Q_t (70,075 ft³/d)：最优f_g ≈ 60%
-
-
敏感性分析：
-
页岩渗透率：页岩层具有渗透性会增加重力超覆，降低波及效率，并改变最优f_g（图4.20）。
-
经济性：若CO₂成本相对于表面活性剂有优势，最优f_g会向更干的方向偏移（图4.21）。
-
CO₂供应限制：若CO₂供应受限，可能无法达到最优f_g，从而降低波及效率（图4.22）。
-
✅ 主要结论
-
超临界CO₂泡沫的优势：其极低的启动压力梯度(∇P₀)是核心优势，使其能在低压力梯度下形成稳定泡沫，并比其它气体泡沫传播得远得多（1-2个数量级）。
-
机理模型的有效性：气泡种群平衡模型能够成功拟合岩心驱替实验数据，并可靠地预测泡沫在储层尺度上的传播距离（受限于泡沫转换和重力超覆）。
-
最优注入策略的确定方法：结合机理模型（提供MRF-f_g关系）和STARS场尺度模拟（提供波及效率等值线图），可形成一种稳健的图形化方法，用于确定特定油藏条件下最优的注入泡沫质量和注入速率。
-
重力与粘性力的竞争：在低注入速率下，重力（水底支撑）对提高采收率贡献更大，最优注入偏向湿泡沫；在高注入速率下，高MRF的粘性指进控制占主导，最优注入偏向更干的泡沫。
-
多因素影响：储层非均质性（页岩隔层渗透性）、经济性（CO₂与表面活性剂相对成本）、以及CO₂供应限制都会显著影响最优注入策略。
🏛️ 作者及单位信息
-
作者: Mohammad Izadi
-
单位: 路易斯安那州立大学（Louisiana State University），Craft & Hawkins石油工程系
-
导师: Seung Ihl Kam 博士
-
委员会成员: Mehdi Zeidouni 博士, Clinton Willson 博士, Jianwei Wang 博士
-
学位: 哲学博士 (Doctor of Philosophy)
-
提交年份: 2019年
致谢：
-
特别感谢Computer Modeling Group (CMG) 为本研究慷慨捐赠CMG-STARS模拟器。
-
感谢Dr. Seung Ihl Kam的悉心指导。
-
感谢Halliburton Scholars Research Experience program, LSU Alumni Association, 和哥伦比亚ANH的间接支持。
💡 补充说明
这篇博士论文代表了泡沫提高采收率数值模拟领域的前沿水平。其核心贡献在于建立了从微观机理（孔隙尺度气泡种群平衡） → 室内实验数据拟合（岩心尺度） → 现场应用优化（油田尺度） 的完整桥梁。研究创新地提出了“波及效率等值线图 + MRF-f_g路径”的图形化优化方法，为泡沫EOR的工程设计和决策提供了实用工具。论文清晰揭示了超临界CO₂泡沫区别于其他气体泡沫的物理机制，并证明了CMG STARS在模拟复杂泡沫流变性和EOR过程中的强大能力。