买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：西南石油大学

摘要：本发明公开了一种考虑油藏CO2气驱渗流时变的数值模拟方法，属于油藏注气开发技术领域。所述数值模拟方法包括以下步骤：S1：获取粘度与残余液相饱和度的函数关系一以及粘度与残余气相相对渗透率的函数关系二；S2：建立油藏CO2气驱数值模拟模型，并对油藏进行网格化以及赋予各个网格单元初始地层油粘度；S3：进行数值模拟，获得下一个时间步各网格单元的原油粘度；S4：确定所述下一个时间步各网格单元的残余液相饱和度和残余气相相对渗透率；S5：更新同时间步网格单元的残余液相饱和度和残余气相相对渗透率；S6：重复步骤S3‑S5，直至数值模拟时间结束，输出累计产油量。本发明能够考虑油藏CO2气驱渗流时变，使得数值模拟的结果更符合实际工况。

主权项：1.一种考虑油藏CO2气驱渗流时变的数值模拟方法，其特征在于，包括以下步骤：S1：获取粘度与残余液相饱和度的函数关系一以及粘度与残余气相相对渗透率的函数关系二；所述函数关系一和所述函数关系二为符合回归分析曲线走势的拟合公式；S2：建立油藏CO2气驱数值模拟模型，并对油藏进行网格化以及赋予各个网格单元初始地层油粘度；所述油藏CO2气驱数值模拟模型包括数值模拟运动方程和数值模拟渗流方程，所述数值模拟运动方程包括：油相数值模拟运动方程： （3）式中：vo为油相的渗流速度，cm3s；K为绝对渗透率，D；KroSg为含气饱和度Sg条件下的油相相对渗透率，无量纲；μo为油相粘度，mPa·s；po为油相压力，10-1MPa；x为坐标轴上的x方向；ρo为油相密度，gcm3；g为重力加速度，10-2cms2；α为流动方向与平面的夹角，°；气相数值模拟运动方程： （4）式中：vg为气相的渗流速度，cm3s；KrgSg为含气饱和度Sg条件下的气相相对渗透率，无量纲；μg为气相粘度，mPa·s；pg为气相压力，10-1MPa；ρg为气相密度，gcm3；残余气相数值模拟运动方程： （5）式中：vge为残余气相的渗流速度，cm3s；Krge为残余气相相对渗透率，无量纲；所述数值模拟渗流方程包括：油相数值模拟渗流方程： （6）式中：为拉普拉斯算子；KroSo为含油饱和度So条件下的油相相对渗透率，无量纲；μop为油藏压力p条件下的油相粘度，mPa·s；Bop为油藏压力p条件下的油相体积系数，m3m3；p为梯度压力，mPa；为地层孔隙度，%；So为含油饱和度，%；t为时间，d；气相数值模拟渗流方程： （7）式中：KrgSo为含油饱和度So条件下的气相相对渗透率，无量纲；μgp为油藏压力p条件下的气相粘度，mPa·s；Bgp为油藏压力p条件下的气相体积系数，m3m3；Swc为束缚水饱和度，%；Rsp为油藏压力p条件下的气体的溶解气油比，m3m3；残余气相数值模拟渗流方程： （8）式中：Sle为残余液相饱和度，%；S3：进行数值模拟，获得下一个时间步各网格单元的原油粘度；S4：根据所述函数关系一和所述函数关系二，确定所述下一个时间步各网格单元的残余液相饱和度和残余气相相对渗透率；S5：根据步骤S4的计算结果更新所述油藏CO2气驱数值模拟模型中同时间步网格单元的残余液相饱和度和残余气相相对渗透率；S6：重复步骤S3-S5，直至数值模拟时间结束，输出CO2气驱数值模拟的累计产油量。

全文数据：

权利要求：

百度查询： 西南石油大学 一种考虑油藏CO₂气驱渗流时变的数值模拟方法