石油学报 ›› 2015, Vol. 36 ›› Issue (3): 361-365.DOI: 10.7623/syxb201503011

• 油田开发 • 上一篇    下一篇

基于数字岩心技术的气体解析/扩散格子Boltzmann模拟

张磊, 姚军, 孙海, 孙致学, 高莹, 赵建林   

  1. 中国石油大学石油工程学院 山东青岛 266580
  • 收稿日期:2014-08-03 修回日期:2014-12-31 出版日期:2015-03-25 发布日期:2015-03-11
  • 通讯作者: 姚军,男,1964年3月生,1984年毕业于华东石油学院采油工程专业,2000年获石油大学(华东)油气田开发工程专业博士学位,现为中国石油大学(华东)教授、博士生导师,主要从事油气田开发工程的教学和科研工作。Email:rcogrf_upc@126.com
  • 作者简介:张磊,男,1984年9月生,2006年毕业于中国石油大学(华东)信息与计算科学专业,现为中国石油大学(华东)油气田开发工程专业博士研究生,主要从事多孔介质微观流动模拟方面的研究。Email:zhlei84@163.com
  • 基金资助:

    国家自然科学基金重点项目(No.51234007)、国家自然科学基金重大项目(No.51490654)、山东省自然科学基金项目(ZR2014EEP018,2013ZRE28068)和中国石油大学(华东)自主创新科研计划项目(14CX06090A)资助

Lattice Boltzmann simulation of gas desorption and diffusion based on digital core technology

Zhang Lei, Yao Jun, Sun Hai, Sun Zhixue, Gao Ying, Zhao Jianlin   

  1. School of Petroleum Engineering, China University of Petroleum, Shandong Qingdao 266580, China
  • Received:2014-08-03 Revised:2014-12-31 Online:2015-03-25 Published:2015-03-11

摘要:

在页岩气和煤层气等非常规储层中普遍存在着气体解析/扩散现象,为了描述微观尺度下的流动模拟,通过在格子Boltzmann方法的演化方程中添加源汇项,模拟了多孔介质中的解析/扩散现象。并以Langmuir等温吸附定律为例进行模拟,模型计算结果与宏观扩散方程的有限差分结果基本一致,但有限差分方法对时间步长要求较高,计算量比格子Boltzmann方法大,格子Boltzmann方法在迭代次数上具有明显的优势。最后基于CT扫描的二维数字岩心,模拟了气体在复杂结构多孔介质中气体解析/扩散过程。这种处理方式避免了复杂的边界条件,计算量不会明显增加,并可以用来模拟具有复杂孔隙结构的真实岩心中的解析/扩散现象,同时可以根据不同储层中气体满足的解析/吸附规律来配置源汇项。

关键词: 数字岩心, 解析/吸附, 扩散, 格子Boltzmann方法, 源汇项

Abstract:

Gas desorption and diffusion phenomenon widely exists in unconventional reservoirs such as shale gas and coalbed methane. In order to describe the flow simulation on a micro-scale, a source-sink term is added to evolution equation of the lattice Boltzmann method, so as to simulate gas desorption and diffusion phenomenon in porous media. Meanwhile, Langmuir isothermal adsorption is taken as a case for simulation. The model calculation results are basically consistent with the finite difference results of macro diffusion equation. However, the finite difference method has a higher requirement for time step with more calculated amount than the lattice Boltzmann method. The latter method has a significant advantage in iterations. Finally, 2D digital core based on CT scanning was used to simulate gas desorption and diffusion process in porous media with complex structure. This processing method is able to avoid complex boundary conditions, and the calculated amount will not increase significantly, but can be used to simulate gas desorption and diffusion phenomenon in real cores with complex porous structure. In the meantime, a source-sink term can be configured according to gas desorption and diffusion laws in different reservoirs.

Key words: digital core, desorption and adsorption, diffusion, lattice Boltzmann method, source/sink term

中图分类号: