CO2驱水二阶PDE两相流模拟：Comsol超负驱替过程与油气藏CCUS研究

张开发
• 2026/4/11 22:58:01 • 15 分钟阅读

最新文章

推荐文章

相关文章

分享文章

CO2驱水二阶PDE两相流模拟：Comsol超负驱替过程与油气藏CCUS研究
comsol co2驱水的二阶PDE两相流 可模拟超负驱替过程。 以定注入速度向地层注CO2采油也适用于油气藏ccus相关研究。CO2驱水模拟在油气田开发中越来越火特别是CCUS碳捕集与利用封存技术搞地下封存的时候。今天咱们来唠唠怎么用COMSOL整活这个两相流模型手把手教你用二阶PDE模拟CO₂把水顶出去的魔性过程。先看核心机制CO₂相和水相在孔隙介质里互相挤兑。这里要玩转两个关键方程——达西定律和质量守恒。举个栗子CO₂的动量守恒方程长这样% 达西速度定义 u_co2 -(k_rel_co2*k_abs/mu_co2)*grad(p) % 质量守恒方程 phi*rho_co2*ds_co2/dt div(rho_co2*u_co2) Q_co2这里的krelco2是相对渗透率函数得用Corey模型或者Brooks-Corey模型来定义。注意孔隙度phi这个参数实测中发现把它的空间导数考虑进去也就是用二阶项能更好地模拟超临界CO₂的骚走位。边界条件设置是个技术活。注入口直接给速度边界// 注入边界 u_n u_inj * (t t_inj) // 时间相关注入速度生产端建议用压力边界避免数值震荡。特别提醒当CO₂饱和度超过临界值时记得触发相对渗透率突变处理不然计算分分钟崩给你看。comsol co2驱水的二阶PDE两相流 可模拟超负驱替过程。 以定注入速度向地层注CO2采油也适用于油气藏ccus相关研究。后处理阶段重点关注两相界面变化。用这个代码抓取饱和度锋面% 提取CO₂饱和度锋面位置 s_co2_threshold 0.15; front_position find(s_co2 s_co2_threshold, 1, last);实测案例显示当注入速度超过临界值时会触发指进现象。这时候网格必须足够细建议用自适应网格配合二阶单元既能保证精度又不至于算到地老天荒。模型验证有个骚操作对比解析解和数值解的饱和度分布曲线。注意看二阶项带来的曲率变化——这货就是模拟超负压驱替的关键所在。当毛细管压力项和二阶导数项开始battle时真正的物理过程才被还原出来。最后给个实用技巧在study步骤里把瞬态求解器的相对容差调到1e-4绝对容差设1e-6时间步长用自动调整。亲测这么搞既能保证收敛又不浪费算力。模型跑完后别急着关记得导出饱和度场数据做CCUS封存量估算这才是甲方爸爸最关心的硬指标。

更多文章