CMG STARS油藏数值模拟入门：蒸汽吞吐建模与收敛问题排查

简介这份PPT是CMG数模软件STARS模块的入门教程面向油气藏数值模拟初学者及相关工程技术人员内容源自中国石油大学北京的培训课程。STARS模块是热采、泡沫及复杂流体模拟的三维数值模拟器教程围绕模拟流程展开重点讲解数据准备与处理涵盖网格与油藏特性定义、流体组分设置、初始条件、井定义与生产动态数据、数值方法控制等核心环节并演示热采模型建立、运行与结果后处理以及氮气及氮气泡沫压水锥数值模拟方法。资源为单个PPTX演示文稿容量5.52MB便于按目录分段学习。目前已有196人浏览学习适合希望系统掌握STARS模块基础操作、为后续复杂油藏数值模拟打下基础的读者参考。1. 先说清楚STARS到底能干什么收到一位做稠油热采的工程师留言说最近在整理课件时翻到一份《CMG数模软件STARS模块使用入门教程.pptx》打开发现里面有些内容丢失还想给文件加个密码问我有没有遇到过类似问题。这让我想起早几年刚接触数值模拟时对着STARS界面一脸懵的日子。先说结论CMG是目前油藏数值模拟领域用得比较多的软件之一而STARS模块是CMG里专门处理热采、化学驱、蒸汽吞吐、SAGD这类复杂开发方式的模拟器。区别于传统黑油模型STARS能模拟温度场变化、多组分相态变化、化学反应、非牛顿流体等简单说黑油模型只关心油、气、水三组分在地下怎么流动而STARS还要额外考虑“温度”“组分”“反应”三个维度。那么这篇内容适合谁看如果你刚接触CMG正在用STARS做蒸汽吞吐、蒸汽驱或者化学驱方案设计或者你只是想看懂别人建的STARS模型到底在算什么这篇内容就是给你准备的。我会把从建模到运行全流程的关键环节拆开讲也会把你可能在现场踩的坑提前告诉你怎么绕开。2. 建模前必须想清楚的三件事很多新手一打开STARS就急着画网格、布井结果跑到一半模型不收敛再回头检查发现问题全出在最前面的设置上。这里的经验是建模之前先花三个小时想清楚问题能省下后面三周的调试时间。2.1 这个模型是干什么用的先说一个最容易犯的错目的不清。同一个油藏你拿来做蒸汽吞吐方案对比和拿来做采收率标定模型的精细程度、网格规模、流体物性处理方式完全不同。举个实际例子如果是做单井蒸汽吞吐历史拟合模型不需要太大通常用径向模型或者小范围扇形模型就够用网格可以适当加密到几米甚至更细但如果是做整个区块的蒸汽驱方案优选模型就得扩大到全区范围网格尺寸自然放宽否则计算量会爆炸。这个道理跟拍照一样拍一张证件照和拍一幅全景图用的镜头和曝光参数能一样吗所以在打开STARS之前先写清楚这几行字模拟目标单井分析 / 井组对比 / 全区方案流体类型稠油 / 超稠油 / 化学驱体系开发方式蒸汽吞吐 / 蒸汽驱 / SAGD / 化学驱时间尺度短期注入阶段 / 长期生产动态这个清单不需要多专业但能逼着你想清楚模型怎么建避免后期返工。2.2 数据准备其实占了七成工作量建模本身的点鼠标操作反而是最简单的真正费时间的是数据准备。STARS需要的数据大概分四类油藏静态数据、流体PVT数据、相对渗透率与毛管力数据、以及动态生产数据。静态数据包括构造顶深、有效厚度、孔隙度、渗透率、含油饱和度、初始压力、初始温度这些通常来自地质建模或测井解释。流体PVT数据是决定模拟成败的关键STARS支持组分模型也就是说你需要定义油、气、水各组分在温度和压力下的物性参数比如粘度、密度、体积系数、溶解气油比等。稠油热采最核心的一个参数是粘温关系STARS内置了一些关联式但最好还是输入实验实测的粘温数据模拟结果才能贴近实际。相对渗透率曲线和毛管压力曲线尽量使用相渗实验和压汞实验的结果。如果没有实验数据也可以用STARS里面自带的相对渗透率模型生成但要意识到这是模型外推结果可信度要打折。最后是动态数据主要是投产以来的产油、产水、注汽量和井底流压等。历史拟合阶段需要用这些数据来反推模型参数数据质量直接影响拟合效果这就好比做菜食材不新鲜厨艺再好也白搭。2.3 组分选择不是越复杂越好STARS属于多组分热采模拟器组分的定义直接影响计算速度和稳定性。有的工程师喜欢把油、气、水定义得很细甚至把C1到C10分开列组分结果模型运行时间翻了几倍还不收敛。这里的原则是能少则少。常规热采模拟中最简单的组分组合是水组分、重油组分、溶解气组分。如果要做SAGD或者蒸汽驱可以考虑加入一个轻烃组分来模拟蒸馏效应。做化学驱再加表面活性剂或聚合物组分。判断标准很简单如果某一组分对你要回答的问题没有直接贡献就别加。组分越多方程组规模越大非线性越强收敛难度呈指数上升。这不是拼配置的游戏目的是把问题算清楚不是把模型弄得很唬人。3. 核心实操从零搭一个蒸汽吞吐模型接下来进入正题我用一个单井蒸汽吞吐模型作为案例把STARS操作流程的关键步骤拆开讲。这个流程是从建网格到出结果的全过程中间会注明每个环节容易出问题的地方。3.1 新建模型与网格设置打开CMG的Builder界面选择STARS模块后第一步是定义网格。蒸汽吞吐单井分析一般用径向柱坐标网格x方向代表径向距离y方向代表角度z方向代表深度。井底附近网格必须加密因为蒸汽注入后井筒附近温度、压力变化最剧烈网格太粗会严重高估热损失。我一般把第一层网格尺寸设为0.5米左右然后按1.2到1.5倍的比例向外扩展这个做法跟有限元分析里网格加密的道理一样物理量变化最剧烈的地方网格分辨率必须跟上。z方向根据油层厚度划分一般每层1到2米油层上下加上盖层和底层用于模拟热损失。盖层和底层的孔隙度、渗透率可以设得很低但热容参数必须给真值因为热损失全靠它们控制。3.2 流体物性与相对渗透率参数网格定义完成后进入组分和流体物性定义环节。以稠油为例核心参数包括油的分子量通常在300到600之间油的临界压力、临界温度、临界体积油在不同温度下的粘度这是最关键的一组数据溶解气油比如果考虑溶解气驱作用岩石和流体的热容、导热系数在STARS里粘度跟温度的关系通常用Andrade或Walther关联式拟合。如果只有几个温度点的实验粘度可以用软件自带工具拟合出完整曲线。注意稠油粘度对温度极其敏感40摄氏度下可能上万厘泊到了250摄氏度可能只有几十厘泊这里差两三个数量级都是正常的拟合结果里如果粘温曲线过于平缓说明实验数据或者拟合参数有误趁早回去核对。相对渗透率方面热采模拟里常用的有三组相渗油水相渗、油气相渗、液气相渗。蒸汽吞吐过程中近井区域存在“热水驱油”和“蒸汽驱油”两种机制两组不同温度下的相渗曲线往往能得到更好的拟合效果。3.3 初始化、射孔与井定义初始化条件包括油藏初始压力梯度、初始温度、初始饱和度分布。在蒸汽吞吐模型里初始温度一般取原始油藏温度初始压力根据深度和流体梯度计算。井的定义是很多新手容易出错的环节。STARS里定义注汽井时除了井轨迹和射孔段还要定义井底约束条件。蒸汽吞吐分三个阶段注汽阶段焖井前注入井定注入速率和注入温度一般注入温度为250到300摄氏度注入速率按井口蒸汽干度折算焖井阶段关井不注不采模拟热量向油层扩散的过程生产阶段把注入井切换为生产井设定井底流压或者最大产液量约束STARS的时间步控制策略是注入阶段用较小的初始步长然后逐步放大焖井阶段可以放宽生产阶段注意观察含水率变化含水率出现剧烈波动时时间步会自动缩短。这些自动控制参数在Numerical Control面板里可以调整。3.4 运行控制与输出的设置运行控制里有几个参数值得专门说Max Timestep Size最大时间步一般初始设为0.01天后期可以逐步放宽Max Timestep Cut最大时间步切割次数默认是10次如果频繁达到这个上限说明模型存在严重的收敛问题Iteration Parameters是指定迭代方法和收敛准则默认的Newton迭代配上相对误差标准在多数情况下够用不建议新手乱改输出设置上我习惯保存每步的网格温度场、压力场、饱和度场这样后续画图分析能看出蒸汽腔的扩展形态排查问题也更直观。4. 模型不收敛手把手帮你排查运行报错是STARS学习中绕不开的坎也是很多入门者最头疼的部分。这里集中把我遇到过的高频问题和排查思路整理成一张速查表方便你对照检查。4.1 高频报错与排查对照现象可能原因检查思路时间步不断缩减计算始终进行不下去数值控制参数过于严苛或物性突变严重先检查有没有组分物性在温度压力范围内出现不合理值提示“Negative composition”或“Saturation out of range”组分设置不合理某组分在计算中消失或为负检查初始饱和度是否合理相渗端点值是否给得太低压力场出现大幅振荡初始压力梯度设置错误或网格物性突变查看初始化压力分布云图确认是否设置了正确的参考深度压力井底压力异常偏低或偏高井底流压约束或表皮因子设置不当检查井的处理方式中是否有非达西流参数运行时间特别长但结果不合理网格数太多或组分太多先跑一个小规模的粗网格模型确认逻辑没问题再加密4.2 收敛问题的三个关键瓶颈第一个瓶颈是物性参数跨数量级突变。稠油热采里粘度随温度变化跨度达四五个数量级如果网格温度梯度太大相邻网格的粘度相差悬殊迭代矩阵刚性很强数值求解器很难收敛。解决办法是网格加密 限制温度变化速率。第二个瓶颈是相对渗透率曲线不光滑。有些实验曲线毛刺多斜率突变点正好落在工作区间内这会让相对渗透率的偏导数值剧烈跳动。遇到这种情况可以用STARS自带的曲线光滑工具或者手动修改端点和曲率让曲线变化更平缓。第三个瓶颈是初始时间步设置太大。默认初始步长如果是0.1天而井底压力从初始值瞬间降到流压约束值这一步的气相饱和度变化量可能超过允许范围直接触发步长切割。我通常先把初始时间步设为0.001天跑过开始几个步再允许逐步放大用一点点计算时间换稳定性非常划算。4.3 一个蒸汽吞吐模型的典型报错实录这里说一个之前帮同事排查过的案例。模型在注汽阶段运行正常一切入生产阶段就疯狂切割时间步屏幕上全是步长收敛警告。排查过程是这样的先圈定问题时间点发现报错集中在注汽转生产的瞬间再看该时间点的输出场发现井底附近的含油饱和度在突然下降且伴随气相饱和度剧烈上升。原因很快清楚切换生产阶段时井底压力瞬间从注入高压降到采出低压近井地带发生了闪蒸轻质组分突然气化导致饱和度场剧烈变化。解决方案是在注入阶段末期增加一个短时间的关井焖井阶段也就是实际工艺中本来就有的焖井过程让压力有一个缓慢释放的过程同时把初始生产阶段的时间步设得更小。改完这两个地方模型一口气跑完再没报错过。5. 结果后处理与模型可信度判断模拟跑完不代表工作结束后处理阶段做得好不好直接影响你对模拟结果的信心。5.1 先看基础指标再看云图打开Results页面我习惯先看几个整体指标累计注汽量、累计产油量、累计产水量、油汽比OSR。油汽比是热采方案里最重要的经济指标之一蒸汽吞吐通常第一周期油汽比最高后面周期逐步下降如果你的模型算出来第二周期油汽比高于第一周期那大概率什么地方出了问题。看完整体指标再看剖面图。温度场分布要看蒸汽腔是否合理扩展饱和度场要不要看剩余油分布压力场看压力传播范围。如果温度场出现不切实际的局部高温团块可能局部网格物性给错了。5.2 历史拟合的常见套路如果是做实际区块的历史拟合STARS模型的拟合参数优先级一般是先调相对渗透率曲线再调粘温关系再调整体渗透率和孔隙度分布最后再考虑局部alteration zones。调参的时候每次只动一组参数观察指标变化不要一次同时改五六个地方否则根本不知道是哪个参数起的作用。5.3 结果输出图和表都要规范CMG的Results可以导出多种格式。我习惯用Report模块生成标准的开发指标报表配合温度场和饱和度场截图。导出图件时注意统一色标范围不同方案对比时尤其重要否则两个方案用了不同色标范围对比就失去了意义。6. 关于课件PPT的两句闲话最后回应一下开头提到的那位工程师。他问的是《CMG数模软件STARS模块使用入门教程.pptx》打开时提示“发现不可读取的内容”以及想给课件加个密码这里一并说下我的处理习惯。PowerPoint提示文件有不可读取内容时绝大多数情况是因为文件里某些格式元素与本机Office版本不完全兼容。最稳妥的处理方式是把文件另存为.pptx副本每次另存后用“文件—信息—检查问题—检查兼容性”跑一遍能提前发现不兼容项。另一个补充建议是养成定期另存为的习惯每改完一轮重要内容就存一个新文件命名加上版本号这套习惯在频繁改课件的场景里非常省心。如果课件涉及甲方数据或者内部研究成果需要设置打开密码在“文件—信息—保护演示文稿—用密码进行加密”里设置即可。密码要单独记好遗忘后处理起来比较麻烦我不建议依赖第三方破解工具合规途径是找回密码或重建文件。7. 学STARS这条路建议这样走如果你刚入门我给的建议是先别碰复杂模型找一个简单的蒸汽吞吐算例完整跑通一遍流程看明白每个模块输入什么、输出什么。第二遍再尝试改参数感受不同参数对结果的影响。第三遍再尝试用实际油藏数据建模。我个人的体会是数值模拟软件只是个计算工具真正决定模型质量的是你对油藏工程原理的理解深度。工具给的是答案而判断答案是否合理靠的是对物理过程的直觉。这种直觉怎么来就是把每一个网格的压力、温度、饱和度变化都能在脑子里还原成地下流体实际流动的画面。等你哪天看到模拟结果能下意识感觉到“这里温度场不对蒸汽可能窜了”说明你才真正上手了。本文还有配套的精品资源点击获取