📝 摘要
通过蒸汽注入的热力采收被认为是实现高最终采收率的有前景方法。在诸如向油藏注入蒸汽等提高采收率项目中，可能会出现复合油藏。热力压降试井分析提供了一种快速估算波及体积和蒸汽区性质的方法。大多数用于分析的模型假设两区复合油藏具有不同但均匀的性质，并由尖锐的垂直界面分隔作为不渗透边界。因此，波及区表现为封闭油藏，压力响应以拟稳态（PSS）行为为特征。由于试井分析方法较水平井更简单，大多数研究考虑的是垂直井。然而，使用水平井对的蒸汽辅助重力泄油（SAGD）工艺是一种有前景的采收技术。对垂直井和水平井（SAGD井对）的蒸汽注入进行了数值模拟研究，以评估热力试井分析方法的适用性和准确性，以及多个参数对结果的影响。初步结果表明获得了相当合理的估算。然而，压力图上观察到的一些趋势无法用现有模型解释，且体积估算存在误差，这可能与常规模型的简化假设有关。
🖥️ CMG软件应用情况总结
| 项目 | 内容 |
|---|---|
| 所用软件 | CMG STARS（2015版本）：先进过程与热采模拟器，用于模拟蒸汽注入和热力压降试井 |
| 模型类型 | • 热采数值模拟模型（考虑热传导、热对流、相变） • 垂直井模型：径向/笛卡尔网格，注入井位于油藏中心 • 水平井模型（SAGD）：笛卡尔网格，注入井位于生产井上方（垂直间距15 ft） |
| 模拟对象 | 加拿大Athabasca油砂区稠油油藏，模拟蒸汽注入后的压降试井过程 |
| 应用方式 | • 压降试井模拟：注入蒸汽直至波及显著岩石体积，关井后读取井筒网格块压力随时间变化 • 历史拟合/输入参数：使用Athabasca油样实测数据（粘度、摩尔质量、密度）作为模拟输入 • 加热器控制：对于水平井SAGD过程，使用STARS的加热器功能预热油藏，建立注入井与生产井之间的连通 • 参数敏感性分析：分析注入时间、注入速率、蒸汽干度、井距、垂向/水平距离等对估算精度的影响 |
| Athabasca油样参数（图1） | • 摩尔质量：534 kg/kgmole • 油密度（标准条件）：1.01286 g/cm³ • 粘度-温度关系：如图1所示 |
| 垂直井模型参数（表1a） | • 油藏面积：~150,000 ft²，厚度：40 ft • 初始压力：700 psi，初始温度：93°F • 孔隙度：0.35，初始水饱和度：21% PV • 水平渗透率：700 md，垂向渗透率：70 md • 注入蒸汽温度：580°F，干度：0.8 • 注入速率：500 STB/D，注入时间：30天，关井时间：24小时 |
| 水平井模型参数（表1b） | • 油藏尺寸：200 ft × 100 ft × 50 ft • 注入井位于生产井上方15 ft • 孔隙度：0.3，初始水饱和度：21% PV • 水平渗透率：7000 md，垂向渗透率：2100 md • 注入蒸汽温度：400°F，干度：0.85 • 注入速率：200 STB/D，注入时间：20天，关井时间：50小时 • 地层体积热容：113 BTU/(ft³-°F)，热导率：78 BTU/(ft-D-°F) |
| 相对渗透率 | 使用Coates等人（1973）的水/油和气/油相对渗透率 |
| 主要结论 | • 垂直井：渗透率估算高估10-20%，波及体积估算误差<15% • 水平井：渗透率估算误差<10%（多数情况），波及体积估算误差<15% • 注入时间过短或过长均影响估算精度（过长导致不规则波及区形状和早期突破） • 注入速率越高，体积估算精度越差，但径向流特征更清晰 • 蒸汽干度增加可改善波及体积估算（但非整个范围） • 常规PSS模型因忽略重力导致的倾斜前缘、热损失和渐变性质，存在局限性 |
文中明确指出（Case125.pdf）：
“The thermal simulator STARS (CMG 2015) was used to simulate the steam falloff tests. Steam is injected into the reservoir models until appreciable rock volumes are swept. Pressure falloff tests are then simulated by shutting the injection well and reading the wellbore gridblock pressures as a function of time.” （第III节，第3页）
“Because of the high viscosity of Athabasca heavy oil and in order to have some initial flow, the reservoir is first heated up to initiate the communication between the horizontal injector and producer. This is done by using the STARS’s heater control.” （第III节，第3-4页）
致谢部分（第6页）：“Authors would like to gratefully thank the Department of Petroleum Engineering and Applied Geophysics at NTNU (Trondheim) for all the support for doing this research.”
🧪 模拟方案与主要结果
1. 试井分析方法（第II节）
垂直井分析（Satman等，1980模型）：
-
早期径向流（半对数直线）→ 计算蒸汽有效渗透率和表皮因子（式1）
-
拟稳态（PSS）流（笛卡尔直线）→ 计算波及体积（式3）
-
两相压缩性（式4-5）：考虑相变引起的体积变化，忽略单相压缩性
水平井分析：
-
可能流动形态：井筒存储（单位斜率）、早期径向流（零斜率）、早期线性流（半斜率）、拟径向流（零斜率）、晚期线性流（半斜率）、PSS流（单位斜率）
-
早期径向流可用于估算渗透率（式7）
-
后期PSS直线用于估算波及体积（同式3）
2. 模拟研究（第III节）
垂直井基案例：
-
注入速率：500 STB/D，注入时间：30天，关井24小时
-
蒸汽温度：580°F，干度：0.8
水平井基案例（SAGD）：
-
注入速率：200 STB/D，注入时间：20天，关井50小时
-
蒸汽温度：400°F，干度：0.85
-
生产井位于底部，注入井在其上方15 ft
-
使用STARS加热器功能预热油藏以建立井间连通
3. 结果与讨论（第IV节）
渗透率估算：
| 井型 | 估算误差 | 说明 |
|---|---|---|
| 垂直井 | 高估10-20%（多数情况） | 可能反映井筒附近高蒸汽饱和度区渗透率 |
| 垂直井（部分情况） | 高度高估 | 井筒存储效应掩盖径向流 |
| 水平井（多数情况） | 误差<10% | — |
波及体积估算：
-
大多数情况下误差<15%
-
极端情况下方法敏感（不规则波及区形状）
-
PSS方法给出的是蒸汽腔+热水区总体积？本研究中蒸汽前缘处流度对比足够高，视为封闭边界，估算体积仅包括蒸汽区
注入时间影响：
-
注入时间过短：体积高估（可能原因）
-
注入时间过长：不规则波及区形状 + 早期突破 → 估算精度下降
-
水平井长注入时间后观察到线性流（图2）
注入速率影响（图3-4）：
-
速率越高：体积估算精度越差
-
但径向流特征更清晰，有时可观察到代表未波及区的第二径向流
-
垂直井和水平井均观察到此行为
其他影响因素：
-
蒸汽干度增加可改善估算（但非整个范围，尤其极端值）
-
生产井与注入井的垂向/水平距离影响估算
-
重力导致的倾斜前缘、热损失、渐变性质未被常规模型考虑
✅ 主要结论
-
CMG STARS的有效性：STARS能够有效模拟蒸汽注入和压降试井过程，包括热传导、热对流、相变和SAGD预热，为热力试井分析方法验证提供了可靠的数值平台。
-
渗透率估算精度：垂直井渗透率估算高估10-20%，水平井误差<10%。高估可能源于井筒存储效应掩盖部分径向流，或估算值反映的是井筒附近高蒸汽饱和度区渗透率。
-
波及体积估算精度：PSS方法在大多数情况下估算误差<15%。但方法对极端情况敏感（不规则波及区形状）。
-
注入时间的关键作用：注入时间过短或过长均影响估算精度。过长注入时间导致不规则波及区形状和早期突破，反而损害估算准确性。
-
注入速率的影响：高注入速率下体积估算精度变差，但径向流特征更清晰，有时可观察到代表未波及区的第二径向流。此行为在垂直井和水平井中均观察到。
-
常规模型的局限性：
-
忽略重力导致的倾斜前缘
-
忽略热损失
-
假设尖锐饱和度梯度（实际为渐变）
-
假设封闭边界（实际边界并非完全封闭）
-
仅考虑PSS直线上数据，忽略低饱和度过渡区
-
-
SAGD试井的特殊考虑：预热是建立井间连通的必要步骤；垂向井距增加可降低早期突破风险，提高估算精度。
🏛️ 作者及单位信息
| 作者 | 单位 |
|---|---|
| Ashkan Jahanbani G. | 挪威科技大学（NTNU），石油工程与应用地球科学系，特隆赫姆，挪威 |
| Tom A. Jelmer | 挪威科技大学（NTNU） |
致谢：感谢NTNU石油工程与应用地球科学系对本研究的支持，感谢Jon Kleppe教授的宝贵意见。
💡 补充说明
该论文是热力试井分析方法在蒸汽驱和SAGD过程中应用的系统性数值模拟研究，核心创新点在于：
-
垂直井与水平井的系统对比：首次在同一框架下系统对比了垂直井和水平井（SAGD）热力压降试井分析的适用性、准确性和误差来源。
-
多参数敏感性分析：系统研究了注入时间、注入速率、蒸汽干度、井距、垂向/水平距离等参数对渗透率和波及体积估算精度的影响，为现场试井设计提供了定量依据。
-
常规模型局限性的识别：明确指出PSS模型因忽略重力导致的倾斜前缘、热损失、渐变性质和实际非封闭边界而存在系统误差，为改进模型提供了方向。
-
SAGD预热模拟：使用STARS的加热器功能模拟SAGD预热阶段，实现了从预热到蒸汽注入再到压降试井的全过程数值模拟。
-
压力导数曲线的深入分析：利用压力导数曲线（log-log图）识别不同流动形态（图2-4），包括井筒存储、径向流、线性流、PSS流等，提高了流动形态识别的准确性。
-
两相压缩性的正确处理：采用Grant和Sorey（1979）的两相压缩性模型（式4-5），考虑相变引起的体积变化，而非忽略相变贡献。
-
实际油样数据的应用：使用Athabasca油样的实测粘度、摩尔质量和密度数据作为模拟输入，增强了研究的实际意义和结果的可信度。
该研究对从事稠油热采、蒸汽驱、SAGD工艺、热力试井设计与解释、以及热采数值模拟的工程师和科研人员具有重要的参考价值。