西南石油大学学报(自然科学版)

• 石油与天然气工程 • 上一篇    下一篇

裂缝性油藏二氧化碳驱的多组分数值计算模型

邸元1 *,吴玉树2,鞠斌山3   

  1. 1. 北京大学工学院,北京海淀100871;2. 美国科罗拉多矿业大学石油工程系,科罗拉多高登80401;
    3. 中国地质大学(北京)能源学院,北京海淀100083
  • 出版日期:2015-04-01 发布日期:2015-04-01
  • 通讯作者: 邸元,E-mail:diyuan@mech.pku.edu.cn
  • 基金资助:

    国家科技重大专项(2011ZX05014);973 计划(2011CB201004);国家自然科学基金(10932001)。

Model Formulation for Compositional Numerical Simulation of CO2
Flooding in Fractured Reservoirs

Di Yuan1*, Wu Yushu2, Ju Binshan3   

  1. 1. College of Engineering,Peking University,Haidian,Beijing 100871,China
    2. Department of Petroleum Engineering,Colorado School of Mines,Golden,Colorado 80401,USA
    3. School of Energy Resources,China University of Geosciences(Beijing),Haidian,Beijing 100083,China
  • Online:2015-04-01 Published:2015-04-01

摘要:

对于预测油藏二氧化碳驱的效果,黑油模型不再适用,必须采用组分模型。研究了裂缝性油藏二氧化碳驱的
三相多组分数值计算模型,给出了相关的计算方法和公式。采用多重连续介质理论来模拟裂缝性油藏中的输运过程,
每个组分的质量传输包括传导和扩散过程。多相流体的流动采用Darcy 定律来描述;组分的扩散采用Fick 定律来描
述。二氧化碳在岩石颗粒上的吸附及油相重组分在岩石骨架上的沉淀采用非线性等温吸附来描述。用串状毛管模型
来考虑沥青质沉积对储层渗透率的影响。采用积分差分法对各组分的质量守恒方程进行空间离散,并按全隐式格式
通过Newton 迭代法进行求解。根据Soave-Redlich-Kwong 状态方程计算流体的物性。对于油–水–气(二氧化碳)三
相的相平衡计算,分别给出了平衡常数法、闪蒸计算法和最小Gibbs 自由能法的计算公式。

关键词: 数值模拟, 组分模型, 二氧化碳驱, 裂缝性油藏, 相平衡

Abstract:

Compositional simulation is required in evaluating CO2 flooding in EOR operations when black-oil simulation is
no longer adequate. A three-phase compositional modeling approach for CO2 flooding in fractured reservoirs is presented. The
mathematical models and formulations are given. A generalized multi-continuum approach is adopted to handle transport in
naturally fractured reservoirs and the mass flux of each mass component is contributed by advection and diffusion processes.
The multiphase fluid flow is described by Darcy equation,and the diffusion of multi-component is described by Fick′s law.
The precipitation of heavy oil components and absorption of CO2 on the solid grains are modeled based on reversible nonlinear
isotherm. The permeability reduction effects due to asphaltene deposition are considered by a tubes-in-series model. The
governing partial differential equations for conservation of each component are discretized using an integral finite differences
method and solved fully implicitly by Newton iteration. The equation of state (EOS) by Soave-Redlich-Kwong is used to
calculate the physical properties of fluids. Three methods of three-phase equilibrium calculation for hydrocarbon-water-CO2
mixtures,K-value approach,flash calculation and Gibbs energy minimization are also presented.

Key words: numerical simulation, compositional model, carbon dioxide flooding, fractured reservoirs, phase equilibrium