一、引言：一个”被判死刑”的油田
在加拿大萨斯喀彻温省的广袤平原之下，埋藏着一个让无数工程师头疼的油藏。这里的原油在常温下呈现出近乎凝固的状态——粘度高达上万厘泊，像蜂蜜一样粘稠。更糟糕的是，油藏底部还存在着活跃的水体，像一个永远关不紧的水龙头，随时准备淹没生产井。
在1987年以前，这个油田的产量低得令人沮丧。石油公司尝试过常规开采，但很快就放弃了——没有经济性的生产意味着没有投资的价值。在当时的技术条件下，蒸汽驱被视为不适当的选择，因为传统的直井无法避免底水锥进问题：一旦开井生产，底部的水就会像锥子一样刺穿油层，直抵井底，让油井迅速”水淹”。
然而，故事的转折点在1987年到来了。水平井技术的成熟，为这个”被判死刑”的油田带来了重生的希望。1991年10月，Jesperson和Fontaine在南萨斯喀彻温省CIM第四届石油会议上，介绍了这个成功的水平井热采先导项目。而我们现在要解读的这份CMG STARS模拟文件（hrw007.dat），正是这个项目的数字化重现。
| 核心主题：本文通过一个完整的CMG STARS数值模拟案例，深度解读如何利用水平井+蒸汽驱技术，实现高粘稠油藏的商业化开发。案例涵盖1982年至1995年共12年的生产历史拟合，包含26口井（含3口水平井）的完整井网配置。 |
二、油田概况：地质与流体的”双重诅咒”
2.1 地质构造与储层非均质性
这个油藏的地质条件可以说是一个”教科书级”的复杂案例。模拟模型采用35×45×5的网格系统，共7875个网格块，采用SI单位制。模型的顶部深度是变化的（Variable Top），从324米到354米不等，反映了构造的起伏特征。
模型的关键地质特征包括：
- 断层与岩性尖灭：在特定区域设置了传导率乘数（Transmissibility Multiplier），如I方向11:13、14:16等区域，这通常用于模拟断层或岩性变化对流体流动的遮挡作用
- 层状非均质性：垂向渗透率设置为水平渗透率的0.35倍，表明这是一个层理发育、垂向沟通受限的油藏
- 活跃的双重含水层：模型包含两个区域性热力含水层（2 Regional Thermal Aquifers），采用半解析法模拟。这意味着油藏的边部和底部存在活跃的水体，对压力保持和热损失都有重要影响
含水层的参数设置尤其值得关注：水层压缩系数为2.03×10⁻⁶ 1/kPa，岩石热容量为2.194×10⁶ J/(m³·°C)，热导率为6.6×10⁵ J/(m·day·°C)。这些参数直接决定了热量在水层中的传播速度和损失程度。
2.2 流体属性：万厘泊级稠油
这个油藏最致命的挑战来自其流体属性。模型包含三个组分：水（H2O）、重油（HOIL）和溶解气（SGAS）。从*VISCTABLE数据中可以清晰地看到原油粘度的温度敏感性：
| 温度 (°C) | 水粘度 (cP) | 油粘度 (cP) | 气粘度 (cP) |
|---|---|---|---|
| 20 | 1.005 | 11,313 | 108 |
| 50 | 0.552 | 756 | 51 |
| 100 | 0.283 | 52.8 | 20 |
| 150 | 0.185 | 11.8 | 10.3 |
| 200 | 0.139 | 3.87 | 5.76 |
| 250 | 0.112 | 2.33 | 4.27 |
| 300 | 0.091 | 1.39 | 3.10 |
从表中可以提取几个关键信息：
- 常温（20°C）下原油粘度高达11,313 cP，属于典型的特稠油范畴
- 温度升高到50°C时，粘度急剧下降至755.5 cP——这就是蒸汽驱的核心原理：通过加热让原油”流动起来”
- 在蒸汽温度（约240°C）下，粘度预计会降至1 cP以下，流动性极大改善
此外，模型的初始条件设定为：油藏温度20°C，参考压力4410 kPa（约44 bar），油水界面深度357.2米，气油界面深度324米。初始含油饱和度约78.72%，含气饱和度约21.28%。
三、技术破局：水平井+蒸汽驱的组合拳
3.1 为什么传统直井蒸汽驱会失败
在水平井技术出现之前，这个油藏的开发面临着两难困境：一方面，原油粘度极高，必须通过加热才能流动；另一方面，油藏底部存在活跃水体，直井生产 inevitably 会导致底水锥进。
具体来说，直井蒸汽驱面临以下技术瓶颈：
- 热损失过大——蒸汽注入后，热量会迅速向上下盖层和底水层散失，热效率极低
- 水锥进严重——直井的生产压差集中在井筒周围，容易将底水直接抽入井底
- 泄油面积有限——单井控制的泄油面积小，难以形成有效的蒸汽驱替前缘
- 经济性差——高投入的热采设施与低产量之间的矛盾，使项目不具备商业价值
3.2 水平井的核心优势
水平井技术的引入，从根本上改变了这个油藏的开发前景。在本案例中，3口水平井的加入带来了多重优势：
- 超大泄油面积：水平井的完井段长度可达数百米（本例中Well-1的完井段横跨I方向14-22网格，约400米），大幅增加了井筒与油藏的接触面积
- 低生产压差：相同的产量下，水平井可以在更低的压差下生产，有效抑制底水锥进
- 均匀驱替：配合蒸汽注入，水平井可以形成更均匀的热前缘推进，提高驱替效率
Jesperson和Fontaine在1991年的会议上报告了先导试验的成功结果：水平井+蒸汽驱的组合不仅显著提高了产量，还有效控制了水油比，使这个曾被认为”不可开发”的油藏具备了商业价值。
四、模拟实战：CMG STARS建模详解
4.1 网格系统与地质建模
本案例的网格系统采用变深度顶面设计（*GRID *VARI），35×45×5的网格规模对于油田级模拟而言属于中等精细度。网格在水平面（I、J方向）上的尺寸是不均匀的——井口附近网格更密（如I方向第10、22、23等位置的网格尺寸为24.86米，远小于默认的49.72米），这是为了在关键区域获得更高的计算分辨率。
模型的边界条件通过两个区域性含水层来实现：
- 含水层1：覆盖I方向19-34、J方向1的区域，沿K方向（垂向）延伸
- 含水层2：覆盖I方向35、J方向5-17的区域，同样沿K方向延伸
含水层采用半解析方法（*AQMETHOD *SEMI-ANALYTICAL）和矩形无限大几何模型（*AQGEOM *RECTANG *INFINITE），这相比于数值含水层计算效率更高，同时能够合理模拟水层的能量供给。
4.2 热力学参数设置要点
热采模拟的核心在于准确地描述热量的传递和损失。本案例的热力学参数设置如下：
| 参数 | 数值 | 单位 |
|---|---|---|
| 岩石热容量 (*ROCKCP) | 2.194×10⁶ | J/(m³·°C) |
| 岩石热导率 (*THCONR) | 6.6×10⁵ | J/(m·day·°C) |
| 水的热导率 (*THCONW) | 5.35×10⁴ | J/(m·day·°C) |
| 油热导率 (*THCONO) | 1.15×10⁴ | J/(m·day·°C) |
| 盖层热容量 (*HLOSSPROP) | 2.62×10⁶ | J/(m³·°C) |
| 盖层热导率 (*HLOSSPROP) | 1.87×10⁵ | J/(m·day·°C) |
| 岩石压缩系数 (*CPOR) | 2.9×10⁻⁶ | 1/kPa |
特别需要注意的是盖层和底层的热损失设置（*HLOSSPROP）。在热采模拟中，向上下围岩的热损失通常占总热损失的30%-50%，直接决定蒸汽驱的热效率和经济性。模型中设置盖层和底层的热容量为2.62×10⁶ J/(m³·°C)，热导率为1.87×10⁵ J/(m·day·°C)，并启用了温度依赖的热损失计算（*HLOSST 20）。
| 实操提示：在实际项目中，热损失参数通常需要通过历史拟合来校正。如果模拟的井下蒸汽干度或温度剖面与实际测井数据不符，首先应该调整*HLOSSPROP中的热导率参数。 |
4.3 井网配置与生产制度
本案例的井网设计体现了从冷采到热采、从直井到水平井的完整转型过程。全模型共26口井，操作方式包括：
| 井名 | 类型 | 操作特征 | 备注 |
|---|---|---|---|
| Well-1 | 水平生产 | I方向完井，长水平段 | 主力生产井 |
| Well-2~3,8~11,13~14,26 | 注汽井 | 温度239.7°C，干度75% | 蒸汽驱热源 |
| Well-4~7,15~25 | 生产井 | BHP下限100-200 kPa | 直井/水平井混合 |
| Well-12 | 水平生产 | I方向完井，13:23网格 | 后期主力井 |
注汽井的操作制度设置为：注入温度239.7°C（对应约3 MPa的饱和蒸汽压力），蒸汽干度0.75，最大井底流压限制14,000 kPa。生产井则设定最小井底压力100-200 kPa，防止出砂和水锥突进。
水平井的完井方式也有讲究。以Well-1为例，其采用I方向完井（*GEOMETRY *I），完井段横跨I方向14至22、J方向15、K方向4的网格，等效井筒半径0.086米， skin 因子0.2488。长水平段的设计使其拥有远超直井的泄油面积。
五、开发历程：12年历史拟合的四个阶段
这份.dat文件最珍贵的地方在于它完整记录了一个油田12年的开发历史。通过解析文件中频繁出现的*OPEN、*SHUTIN和*ALTER指令，我们可以清晰地还原出四个阶段的开发历程。
5.1 冷采试产期（1982-1986）
1982年12月，模拟启动。此时油田处于原始状态，温度20°C，所有井都处于关井状态。1983年1月，Well-6最早投产，但初始产量极低（约束为0.43 m³/d），反映了冷采条件下稠油的流动性之差。
1984年1月，Well-5、7、16、18、19、25陆续开井，但大多数井的产量约束仅为0.001 m³/d——几乎只是维持生产状态。由于油稠且底水活跃，部分井（如Well-8、25）在1986年被关停。这个阶段的惨淡经营，直观地展示了为什么仅凭冷采无法开发这个油藏。
5.2 转型阵痛期（1987-1988）
1987年7月是一个里程碑时刻——水平井Well-1正式投入生产。从数据中可以观察到其产量的剧烈变化：初始11.55 m³/d，到1988年1月飙升至124 m³/d，再到1988年7月达到324.5 m³/d。这种产量的大幅提升，正是水平井巨大泄油面积的直接体现。
然而，转型并非一帆风顺。产量的剧烈波动暗示着水锥进和系统压力调整的挑战——Well-5和Well-16在此期间被关停。这也提醒我们：水平井虽然强大，但在强底水油藏中仍需要精细的操作控制。
5.3 热采全面见效期（1989-1991）
1989年1月是另一个关键转折点——Well-2正式转为注汽井（*INJECTOR），标志着油田从冷采全面转向热采。注汽参数设定为：温度239.7°C，蒸汽干度75%，最大注入压力14,000 kPa。
随着蒸汽的注入，油层温度逐渐升高，原油粘度显著下降，产量开始全面提升。1990年7月Well-12开井后产量迅速攀升，1991年7月达到575.9 m³/d。同年，Well-10、11也转为注汽井，形成了较为完善的注采井网。
这个阶段的成功充分证明：蒸汽驱+水平井的组合拳能够有效攻克高粘稠油藏的开发难题。热量降低了原油粘度，水平井提供了充足的泄油通道，两者相辅相成。
5.4 加密与稳产期（1992-1995）
1992年至1993年，项目进入加密调整阶段。大量新井投入生产或注汽：Well-13、14转为注汽井，Well-15、24、26陆续开井。井网密度达到最高，全油田共有20余口井同时运行。
1994年进行了大规模的井网调整。3月部分老井（Well-1、8、11、14）短暂关停后重新开启，这通常对应着修井作业或井况转换。7月后全井网恢复运行，产量在高位保持稳定。
1995年1月，历史拟合阶段结束。此时全油田处于高密度井网运行状态，水平井+蒸汽驱的开发模式已经完全成熟。文件中预留了至2004年的*DATE卡片，只需取消注释即可进行未来产量预测。
六、操作指南：如何复现与扩展该模拟
6.1 快速上手步骤
对于想要学习和复现这个案例的工程师，建议按照以下步骤操作：
- 环境准备——确保已安装CMG Launcher，并拥有STARS模块的许可。将hrw007.dat文件及所有关联的.inc文件（如CMGTEMP_000.inc至CMGTEMP_003.inc）放在同一目录下
- 首次运行——在CMG Launcher中选择STARS模块，导入hrw007.dat文件。初次运行建议保持默认参数，观察12年历史拟合的运行过程和结果
- 结果分析——重点关注三个指标——全油田产油量、水油比、井底温度剖面。在历史拟合良好的情况下，模拟产量应与实际数据高度吻合
- 预测扩展——取消1995年之后的*DATE卡片注释（从”*DATE 1995 06 30″开始），运行预测模式，评估未来10年的开发效果
- 敏感性分析——尝试修改关键参数（如水平井长度、注汽速率、注汽干度），观察对产量的影响，加深对热采机理的理解
6.2 关键参数调试技巧
在实际应用中，以下几个参数是历史拟合的关键调控点：
| 参数 | 影响 | 调试建议 | ||
|---|---|---|---|---|
| *PERMI（渗透率） | 决定流体流动能力 | 先降低整体渗透率，再微调局部Multiplier | ||
| *HLOSSPROP（热损失） | 影响温度场分布 | 热损失大则前缘推进慢，产量偏低 | ||
| *KRTEMTAB（温度依赖Krg） | 决定热采效率 | 高温Krg增大→水相流动性增强→含水上升 | ||
| *TRANSI/J/K（传导率） | 控制流体路径 | 局部调高可模拟裂缝或高渗通道 | ||
| 井筒模型参数 | 影响单井产能 | skin因子、等效半径为主要调试对象 | ||
| 经验之谈：热采模拟的历史拟合通常遵循”先压力、后温度、最后饱和度”的顺序。首先确保全区压力趋势与实际吻合，然后调整热损失参数使温度剖面匹配，最后微调相对渗透率曲线来拟合含水上升规律。 | ||||
七、关键启示与技术要点总结
通过对这个经典案例的深度解读，我们可以提炼出以下对实际项目具有指导价值的关键启示：
1. 技术适用性判断
水平井+蒸汽驱的组合特别适用于以下类型的油藏：原油粘度在冷采条件下超过1000 cP、油层厚度适中（过薄则热损失过大）、底水活跃但并非极端强烈的油藏。如果底水过于活跃，即使水平井也难以控制水锥进。
2. 历史拟合的价值
本案例包含12年的完整生产数据，是极其宝贵的学习资源。历史拟合的过程就是不断校正模型参数、加深对油藏认识的过程。特别是在热采项目中，温度依赖的相对渗透率曲线（如*KRTEMTAB）往往难以通过实验直接获取，必须通过历史拟合来反演确定。
3. 动态调整的重要性
从文件中频繁的*OPEN、*SHUTIN和*ALTER指令可以看出，热采开发绝不是”设定好就一劳永逸”的事情。开发过程中需要根据生产动态不断调整井网、转换井别、优化注采参数。这种灵活性是热采项目成功的关键。
4. 经济性的终极考量
技术可行不等于经济可行。在评估类似项目时，必须综合考虑：蒸汽发生和注入成本、水平井钻井完井成本、油藏条件对热效率的影响、环保要求（特别是温室气体排放）等因素。CMG模型提供了产量预测基础，但最终的决策还需要经济评价模型的支撑。
| 维度 | 核心要点 |
|---|---|
| 地质认识 | 双重含水层+层状非均质性是压力与热动态的主控因素 |
| 技术选择 | 水平井解决水锥进问题，蒸汽驱解决粘度问题，两者缺一不可 |
| 模型构建 | 35×45×5网格+半解析水体，平衡了精度与计算效率 |
| 历史拟合 | 12年数据分四阶段拟合，优先匹配压力→温度→含水率 |
| 工程操作 | 频繁的开井/关井/调整反映了热采开发的动态管理需求 |
八、结语
这个来自CMG STARS的经典案例（hrw007.dat），远不止是一组冷冰冰的模拟参数。它是一个关于技术创新如何改变资源开发命运的真实故事——从”被判死刑”到”起死回生”，水平井+蒸汽驱的组合为高粘稠油藏的开发开辟了全新的道路。
对于从事油藏数值模拟工作的工程师来说，这个案例提供了极其丰富的学习素材：从地质建模的技巧，到热力学参数的设置；从井网设计的思路，到历史拟合的方法论。每一个*OPEN、*SHUTIN的背后，都蕴含着开发策略的深思熟虑。
更重要的是，这个案例传递了一种工程思维：面对看似无解的技术难题，不要急于否定，而应该积极探索新技术的组合应用。水平井技术的成熟，让曾经”不合适”的蒸汽驱方案重新焕发了生命力。这种”旧问题+新方法”的思维模式，或许正是工程创新的本质所在。
最后，希望每一位阅读本文的工程师，都能从这个案例中汲取到对自己工作有益的启示。毕竟，每一个.dat文件的背后，都藏着一个值得被讲述的故事。