📝 摘要
深部咸水层的 CO₂ 封存被认为是减少温室气体的重要途径。然而，许多已发表的关于封存潜力的报告未能认识到将 CO₂ 封存在封闭系统中的必要性。通过解析模型，本研究表明所需的封存体积和注入井数量远超预期，这使得地质 CO₂ 封存成为一个极其不可行的减排选项，除非同时采出 brine 以创造空隙和压力释放。解析模型结果与 CMG 数值模拟器（包括 CO₂ 封存的多相流物理过程）吻合良好。
含水层压力升高威胁密封完整性并带来 CO₂ 泄漏风险。因此，长期监测 CO₂ 封存储层的完整性对于确保地质封存的安全性和可接受性至关重要。本研究讨论了 CO₂ 封存储层压力降落测试的重要方面，重点关注储层压力监测和泄漏检测，并探讨了在商业封存项目中定期进行压力降落测试以及如何帮助诊断含水层泄漏的重要性。
为消除自由气在顶部密封下方的积聚，本研究提出了一种工程化 CO₂ 注入系统，通过采出 brine 来控制含水层压力升高。该系统消除了浮力驱动的自由气积聚，并通过采出 brine 避免了含水层压力升高。对 30 年 CO₂ 注入及随后 1000 年自然梯度流的模拟表明，与纯 CO₂ 注入相比，工程化系统可以更安全地将 CO₂ 封存在相对较小的封闭含水层体积中，并具有更大的封存潜力。
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | • CMG GEM-GHG：温室气体封存模块，用于 CO₂ 在深部咸水层中封存的多相流、地球化学反应、相态行为和组分输运的数值模拟（论文第21页，第2章）。 • Kappa Engineering：用于压力瞬态分析（论文致谢部分）。 |
| 模型类型 | 三维均质/非均质模型，采用 Peng-Robinson 状态方程（PR-EOS）模拟 CO₂-brine 体系。 |
| 模拟对象 | 深部咸水层（均质和随机渗透率场），模拟 CO₂ 注入、brine 采出、brine 回注以及工程化系统（brine 注入井+CO₂ 注入井+brine 生产井）对 CO₂ 封存的影响。 |
| 应用方式 | • 相态与溶解度建模：使用 PR-EOS 模拟 CO₂-brine 体系；使用 Henry 定律（Harvey 相关性）计算 CO₂ 在 brine 中的溶解度，该溶解度是压力、温度和盐度的函数；使用 Harvey 盐析系数校正盐度对溶解度的影响。 • 残余气 trapping 建模：使用 Land 模型（Land，1968）模拟相对渗透率滞后和残余气 trapping，最大残余气饱和度由 Holtz（2002）相关性计算（孔隙度、渗透率的函数）。 • 多组分输运与反应：模拟 CO₂ 在气相和水相之间的分配、水蒸发、以及 CO₂ 溶解对 pH 的影响。 • 封闭/开放含水层对比：模拟封闭（无流边界）和开放（恒压边界，通过 VOLMOD 关键字在边界块施加巨大孔隙体积乘数实现）两种边界条件下 CO₂ 注入的压力响应。 • 工程化系统设计：设计了包含顶部 brine 注入井、中部 CO₂ 注入井和底部 brine 生产井的工程化系统（图6.1），模拟了 30 年注入期内 brine 循环对 CO₂ 溶解和 trapping 的加速作用，并模拟了 1000 年的自然梯度以评估长期封存安全性。 • 压力降落测试模拟：模拟了不同注入时间的压力降落测试，利用压力导数曲线识别干涸区、两相区和 brine 区，并用于检测含水层泄漏。 • 非均质性和倾角敏感性分析：使用随机渗透率场（序贯高斯模拟生成）和有倾角地层，测试工程化系统在更真实地质条件下的性能。 |
| 模型参数（基础案例） | • 网格：50×50×50（附录输入文件，第171页） • 模型尺寸：25,000 ft × 25,000 ft × 1,000 ft • 深度：6,000 ft • 初始压力：2,600 psi • 温度：150°F • 盐度：100,000 ppm • 孔隙度：0.25 • 渗透率：100 md（均质），kz/kx=0.1 • 注入量：52 MMscf/D（~100 万吨 CO₂/年），注入 30 年 • 束缚水饱和度：0.25 • 临界气饱和度：0.25 • 最大残余气饱和度：0.3（Land 模型） |
| 模拟方案 | • 边界条件对比：无流边界（封闭）vs. 恒压边界（开放） • 注入井数量：不同含水层尺寸（5-100 英里）下的 BHP 和注入能力对比 • 压力降落测试：不同注入时间（20 天、30 年）后的压力降落测试 • 泄漏检测：在干涸区、两相区和 brine 区设置泄漏，分析压力响应 • brine 采出：添加 brine 生产井控制压力升高 • 工程化系统：顶部 brine 注入井 + CO₂ 注入井 + 底部 brine 生产井（图6.1） • 非均质性：随机渗透率场（1-200 md） • 倾角：4° 倾角地层 |
| 主要结论 | • 封闭含水层的压力升高：与开放（恒压边界）含水层不同，封闭（无流边界）含水层中 CO₂ 注入会导致平均储层压力随时间线性升高。BHP 在早期因干涸区扩展而略有下降，但一旦压力波到达边界，BHP 开始上升。 • 存储潜力被高估：使用开放边界假设的模型严重高估了存储潜力。对于典型的 500 MW 燃煤电厂（3 MMT CO₂/年，30 年），所需的含水层面积可能高达 13,800 平方英里，或需要数百口注入井。 • 解析模型验证：开发的封闭含水层解析模型（考虑三区复合模型）与 CMG GEM 数值模拟结果吻合良好，可用于快速评估存储潜力和所需井数（图3.7、3.8）。 • 压力降落测试的价值：定期压力降落测试可以追踪干涸区和两相区半径的扩展、计算各区的流度，并通过 Horner 分析估算平均储层压力。压力降落测试还可以检测含水层泄漏的位置（干涸区、两相区或 brine 区）并估算泄漏量（图4.18-4.28）。 • brine 采出的效果：采出 brine 可显著降低压力影响范围（从 23 英里半径降至 0.5 英里），从而降低泄漏风险（图5.1、5.2）。 • 工程化系统的效果：设计的工程化系统（顶部 brine 注入+底部 brine 生产）使垂直气体速度降低 3 倍，水平速度提高 3 倍，提高了波及效率（图6.6、6.7）。注入结束后 20 年内，90% 的 CO₂ 被溶解或 trapping（图6.11）。工程化系统的存储潜力（8.5%）是纯 CO₂ 注入（~1%）的 4 倍。 • CMG 适用性：CMG GEM-GHG 能够有效模拟 CO₂ 封存中的多相流动、溶解度、残余气 trapping、压力响应和复杂工程方案，是 CO₂ 地质封存研究和风险评估的有力工具。 |
文中明确指出（Case140.pdf）：
“The results from our analytical model match well with a numerical reservoir simulator including the multiphase physics of CO₂ sequestration.” （第3页，Abstract部分）
“I would like to express my appreciation to the Computer Modeling Group (CMG) and Kappa Engineering for providing the uninterrupted license of the software packages. I am grateful to Mr. Bob Brugman from CMG for his help with the GEM-GHG (Green House Gas) module of CMG.” （第6页，Acknowledgements部分）
“In this section, some of the important features of CMG’s generalized EOS model simulator for green house gases (GEM GHG) are applied to set up the base case simulation parameters used throughout this study.” （第21页，第2章引言）
“The GEM simulator was used with the GHG option. Simulations for 30 years of CO₂ injection followed by 1,000 years of natural gradient flow were made.” （第26页，第2.4节）
🧪 模拟方案与主要结果
1. 基础模型与封闭/开放含水层对比（第2.4节、第3.1节，图3.1、3.2）
模型参数：25,000×25,000×1,000 ft，k=100 md，φ=0.25，注入 52 MMscf/D（~1 MMT/yr），30 年。
| 边界条件 | 实现方法 | BHP 行为 | 平均压力行为 | 存储潜力 |
|---|---|---|---|---|
| 开放（恒压） | 边界块使用大孔隙体积乘数（VOLMOD） | 随时间下降（干涸区扩展提高注入能力） | 几乎不变 | 被高估 |
| 封闭（无流） | 默认无流边界 | 早期下降，压力波到达边界后上升 | 线性上升 | 实际有限 |
关键结论：使用开放边界假设的模型严重高估了 CO₂ 封存潜力。封闭含水层中的压力升高会限制注入量，并可能超过裂缝压力。
2. 压力降落测试与泄漏检测（第4章，图4.5-4.28）
模拟方法：在不同注入时间后模拟压力降落测试，使用压力导数曲线识别三区复合模型。
| 区域 | 压力导数水平 | 对应物理区域 |
|---|---|---|
| 第一水平段 | 最低 | 干涸区（单相 CO₂，高流度） |
| 第二水平段 | 中等 | 两相区（CO₂+brine，中等流度） |
| 第三水平段 | 最高 | brine 区（单相水，低流度） |
压力降落测试应用：
-
追踪前缘：通过压力导数偏离时间计算干涸区半径和两相区前缘（图4.9）
-
计算流度：从压力导数水平值计算各区流度（图4.10-4.12）
-
平均压力估算：Horner 分析可估算平均储层压力，区分开放/封闭边界（图4.13-4.16）
-
泄漏检测：不同位置（干涸区、两相区、brine 区）的泄漏会导致 Horner 曲线特征变化，可检测泄漏并估算泄漏量（图4.18-4.28）
3. 工程化系统设计（第6章，图6.1-6.31）
系统设计（图6.1）：顶部 brine 注入井（I_B）+ 中部 CO₂ 注入井（I_CO₂）+ 底部 brine 生产井（P_B）
| 参数 | 均质基础案例 | 工程化系统 | 改善效果 |
|---|---|---|---|
| 垂直气体速度 | 0.03 ft/D | 0.01 ft/D | 降低 3 倍 |
| 水平气体速度 | 0.01 ft/D | 0.03 ft/D | 提高 3 倍 |
| 重力数（Ngv） | 1.25 | 0.5 | 降低，波 及效率提高 |
| 存储潜力 | ~1% | 8.5% | 提高 4 倍 |
| 1000 年后自由气 | 大量积聚在顶部密封下方（图6.8a） | 无自由气（图6.8b） | 消除泄漏风险 |
| 50 年后固定化 CO₂ | ~30% | 90% | 加速溶解和 trapping |
非均质性和倾角敏感性（第6.5节，图6.16-6.31）：
-
工程化系统在随机渗透率场（1-200 md）中仍然有效（图6.22、6.23）
-
在 4° 倾角地层中同样有效（图6.29、6.30）
-
垂直渗透率>25 md 时系统效率下降（图6.13）
-
水平渗透率影响较小，但需要调整井距以避免 CO₂ 突破（图6.14）
✅ 主要结论
-
封闭含水层压力升高是 CO₂ 封存的关键限制因素：与开放（恒压边界）含水层不同，封闭（无流边界）含水层中 CO₂ 注入会导致平均储层压力随时间线性升高。BHP 在早期因干涸区扩展而略有下降，但一旦压力波到达边界，BHP 开始上升。使用开放边界假设的模型严重高估了存储潜力。
-
解析模型与 CMG 模拟验证：开发的封闭含水层解析模型（考虑干涸区、两相区和 brine 区的三区复合模型）与 CMG GEM 数值模拟结果吻合良好，可用于快速评估存储潜力和所需井数。
-
存储潜力被严重高估：对于典型的 500 MW 燃煤电厂（3 MMT CO₂/年，30 年），所需的含水层面积可能高达 13,800 平方英里（比 Prudhoe Bay 油田面积大 40 倍），或需要数百口注入井。纯 CO₂ 注入的存储潜力仅约 1%。
-
压力降落测试是有效的监测工具：定期压力降落测试可以：
-
通过压力导数偏离时间追踪干涸区和两相区半径的扩展
-
从压力导数水平值计算各区的流度
-
通过 Horner 分析估算平均储层压力，区分开放/封闭边界
-
检测含水层泄漏的位置（干涸区、两相区或 brine 区）并估算泄漏量
-
-
brine 采出可控制压力升高：采出 brine（体积等于注入 CO₂）可将+50 psi 超压等值线范围从 23 英里半径缩小至 0.5 英里，显著降低泄漏风险。
-
工程化系统（brine 循环）大幅提高封存安全性和效率：
-
顶部 brine 注入井使垂直气体速度降低 3 倍，延迟 CO₂ 到达顶部密封的时间
-
底部 brine 生产井使水平速度提高 3 倍，提高波及效率
-
注入结束后 20 年内，90% 的 CO₂ 被溶解或 trapping（而纯 CO₂ 注入 1000 年后仍有大量自由气）
-
存储潜力从 1% 提高到 8.5%（提高 4 倍）
-
-
工程化系统的敏感性：
-
垂直渗透率 >25 md 时系统效率下降，需要更高的 brine 注入速率
-
水平渗透率对最优 CO₂ 注入速率影响较小，但需调整井距以避免突破
-
系统在随机渗透率场和 4° 倾角地层中仍然有效
-
-
CMG 软件的适用性：CMG GEM-GHG 能够有效模拟 CO₂ 封存中的多相流动（PR-EOS）、溶解度（Henry 定律，Harvey 相关性）、残余气 trapping（Land 模型）、压力响应和复杂工程方案，是 CO₂ 地质封存研究和风险评估的有力工具。
🏛️ 论文主要信息
| 项目 | 内容 |
|---|---|
| 论文标题 | Aquifer Management for CO₂ Sequestration |
| 作者 | Abhishek Anchliya |
| 学位 | 理学硕士（Master of Science） |
| 授予单位 | 德克萨斯农工大学（Texas A&M University） |
| 毕业年份 | 2009 年 12 月 |
| 所属院系 | 石油工程系（Department of Petroleum Engineering） |
| 研究领域 | CO₂ 地质封存、深部咸水层、压力管理、压力降落测试、泄漏检测、工程化注入系统 |
💡 补充说明
该硕士论文是一项关于 CO₂ 地质封存中压力管理和泄漏检测的系统性数值模拟研究，核心创新点在于：
-
封闭含水层压力行为的正确建模：首次明确指出大多数已发表研究错误地使用了开放（恒压）边界假设，严重高估了 CO₂ 封存潜力。通过 CMG 模拟和解析模型，证明了封闭含水层中压力会线性上升，这是 CO₂ 封存的关键限制因素。
-
解析模型的开发与验证：将 Burton 等人的开放系统解析模型扩展到了封闭系统，考虑了干涸区、两相区和 brine 区的三区复合流动。模型与 CMG GEM 模拟结果吻合良好，可用于快速评估存储潜力和所需井数。
-
压力降落测试在 CO₂ 封存中的应用：首次系统研究了压力降落测试在 CO₂ 封存监测中的应用，包括追踪前缘扩展、计算区流度、估算平均压力、检测泄漏位置和泄漏量。为 CO₂ 封存项目的长期监测提供了实用的技术方案。
-
工程化注入系统的创新设计：提出了“顶部 brine 注入 + 中部 CO₂ 注入 + 底部 brine 生产”的工程化系统，同时解决了两个核心风险：通过 brine 采出控制压力升高，通过顶部 brine 注入消除自由气在密封下方的积聚。模拟表明该系统可在注入结束后 20 年内固定 90% 的 CO₂。
-
加速溶解和 trapping 的量化：首次定量比较了工程化系统与纯 CO₂ 注入在溶解和 trapping 速率上的差异（90% vs 30% 在 50 年内），以及存储潜力的差异（8.5% vs 1%）。
-
非均质性和倾角敏感性分析：使用随机渗透率场（序贯高斯模拟生成）和有倾角地层，验证了工程化系统在更真实地质条件下的有效性，并确定了垂直渗透率的阈值（25 md）。
-
完整的 CMG 工作流：从基础相态建模（PR-EOS）、溶解度模型（Henry 定律）、残余气 trapping（Land 模型）到工程化系统设计和长期预测，完整展示了 CMG GEM-GHG 在 CO₂ 封存研究中的应用。
该研究对从事 CO₂ 地质封存、温室气体减排、油藏工程、压力瞬态分析和多相流模拟的工程师和科研人员具有重要的参考价值。