没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
采用格子Boltzmann方法(LBM)研究了水平通道内非混相驱替过程中的黏性指进现象.考察流动方向与重力方向垂直时,驱替和被驱替流体间黏度比及壁面可润湿性的影响,比较了不同条件下的指进前缘形状,并利用突破时间和面积扫掠效率分析了驱替效率.研究表明,在重力的影响下,指进前缘形状关于水平通道的中心线是非对称的.随着非混相流体间黏度比的增大,指进现象也随之增强,导致驱替效率降低.当驱替流体不易润湿通道壁面时,指进现象被强化;当驱替流体易于润湿通道壁面时,指进现象则被抑制.
资源推荐
资源详情
资源评论
第52卷第3期
2 0 1 2 年 5 月
大 连 理 工 大 学 学 报
Journal of Dalian University of Technology
Vol .52 , No .3
May 2 0 1 2
文章编号 :1000‐8608(2012)03‐0343‐07
水 平 通 道 内 黏 性 指 进 现 象 LB M 模 拟
李 维 仲
倡
, 董 波 , 宋 永 臣
( 大连理工大学 海洋能源利用与节能教育部重点实验室 ,辽宁 大连 116024 )
摘要 :
采用格子 Boltzmann 方法(LBM )研究了水平通道内非混相驱替过程中的黏性指进现
象 .考察流动方向与重力方向垂直时 ,驱替和被驱替流体间黏度比及壁面可润湿性的影响 ,比
较了不同条件下的指进前缘形状 ,并利用突破时间和面积扫掠效率分析了驱替效率 .研究表
明 ,在重力的影响下 ,指进前缘形状关于水平通道的中心线是非对称的 .随着非混相流体间黏
度比的增大 ,指进现象也随之增强 ,导致驱替效率降低 .当驱替流体不易润湿通道壁面时 ,指
进现象被强化 ;当驱替流体易于润湿通道壁面时 ,指进现象则被抑制 .
关键词 :黏性指进 ;非混相流体 ;表面可润湿性 ;黏度比
中图分类号 :TE319 文献标志码 :A
收稿日期 :2010‐10‐11 ; 修回日期 :2012‐03‐25 .
基金项目 :国家自然科学基金资助项目 (50976017 ;50736001) .
作者简介 :李维仲
倡
(1956‐) ,男 ,教授 ,博士生导师 ,E‐mail :wzhongli@ dlut .edu .cn .
0 引 言
黏性指进现象广泛存在于油气开采过程中 ,
是多相渗流的重要特征之一 .指进现象的出现会
导致驱替流体的过早突破 ,使得波及体积减小 ,原
油采收率降低 ,因此研究者对这一现象进行了大
量的研究 .线性稳定性分析方法被广泛用来分析
毛细扩散 、黏度比以及岩石可润湿性对指进现象
的影响 ,但至今 ,还没有一套基于连续介质假设的
相对成熟的理论来描述指进现象的演化过程和定
量地刻画其形态特征
[1 、2]
.此外基于分形理论的
扩散限制凝聚模型
[3]
也被用来模拟多孔介质内非
混相驱替过程中的指进现象 ,虽然取得了一定进
展 ,但是还需进一步探讨分形理论的物理意义与
黏性指进的物理机制的内在联系
[4]
.
格子 Boltzmann 方法 (LBM )是一种介于微
观的流体分子动力学模型和宏观的连续介质模型
的介观方法 .该方法算法简单 ,易于程序化处理 ,
具有良好的并行特性 ,近十几年来发展非常快 .
LBM 已经被广泛地应用于多相多组分流体系统
及工程传热传质问题等多个领域 ,包括黏性指进
等复杂现象 .Chin 等
[5]
研究了二维通道内的两组
分流体驱替问题 ,分析了不同黏度比对驱替过程
的影响 .Grosfils 等
[6]
通过分析密度等值线讨论
了 Hele‐Shaw 槽 的 指 进 不 稳 定 性 问 题 .Kang
等
[7]
考察了二维通道内黏度比 、毛细管数 、表面可
润湿性对指进现象的影响 .Huang 等
[8]
研究了二
维多孔介质内的驱替过程中 ,黏性比以及壁面可
润湿性对黏性耦合问题的影响 .Boek
[9]
模拟了二
维通道内的黏性指进现象 ,以及多孔介质内的两
相驱替问题 .这些研究讨论了黏度比的影响 ,但是
大多忽略了重力作用 .Chau 等
[10]
考察了二维多
孔介质内排驱过程中 ,流体流动方向和重力方向
相同时的情况 .然而 ,流体流动方向与重力方向垂
直情况下的研究尚未见报道 .本文着重讨论水平
通道内非混相驱替过程中黏性指进现象 ,考察当
流动方向与重力方向垂直时 ,黏性比以及壁面可
润湿性对指进形成 、前缘形状以及驱替效率的影
响 .
1 LBM 简介
本文采用 Shan 和 Chen 于 1993 年
[11]
提出的
多组分模型(伪势模型)来进行两相流体的驱替过
程模拟 .该模型中 ,对于模拟非混相流体系统而
言 ,流体间界面可以自动形成 ,而不需要像传统方
法那样进行复杂的界面重构或者界面捕捉 .不同
的流体组分采用不同的流体粒子分布函数 ,每类
流体的粒子分布函数的演化方程为
f
σ
α
(x
+
e
α
δ
t
,t
+
δ
t
) =
f
σ
α
(x ,t) -
1
τ
σ
[
f
σ
α
(x ,t) -
f
σ
(eq)
α
(x ,t)] (1)
式中 :
f
σ
α
(x ,t) 是第
σ
类流体在 t 时刻 、x 处的粒子
分布函数 ,
f
σ
(eq)
α
(x ,t) 为第
σ
类流体对应的平衡态
粒子分布函数 ;
τ
σ
是第
σ
类流体的松弛时间 ;e
α
是
粒子离散速度向量 ,
α
为离散速度方向 .
每类流体粒子的平衡态粒子分布函数为
f
σ
(eq)
α
=
w
α
ρ
σ
1
+
e
α
·
u
eq
σ
c
2
s
+
(e
α
·
u
eq
σ
)
2
2c
4
s
-
u
eq
σ
·
u
eq
σ
2c
2
s
(2)
式中 :c
s
=
c/ 3 ,是格子声速 .w
α
为权系数 ,对于
D2Q9 模型 ,当
α
=
0 时 ,w
α
=
4/9 ;当
α
=
1
~
4
时 ,w
α
=
1/9 ;当
α
=
5
~
8 时 ,w
α
=
1/36 .
该模型中 ,假设流体粒子之间存在非局部的
相互作用势 ,如果只考虑邻近和次邻近粒子之间
的相互作用力 ,第
σ
类流体粒子受到的作用力
F
σ
f
(x) 可表示为
[12 ]
F
σ
f
(x) = -
ψ
σ
(x)G
σσ
∑
x
′
∑
σ
ψ
σ
(x
′
)e
α
(3)
式中 :x
′
=
x
+
e
α
,为 x位置在第
α
方向上的邻近点
位置 ;而
ψ
σ
(x) 称为第
σ
组分的“有效密度” ,是局
部密度的函数 ,这里采用
ψ
σ
(x) =
ρ
σ
.参数 G
σσ
用
来控制相互作用力的强度 ,对于 D2Q9 模型 ,当
e
α
=
0 时 ,G
σσ
=
0 ;当 e
α
=
1 时 ,G
σσ
=
G ;当
e
α
=
2 时 ,G
σσ
=
G/4 .对于非混相流体系统 ,
当
σ
=
σ
时 ,G
=
0 ;而当
σ
≠
σ
时 ,G 为正值 .通过
调节 G 的取值大小 ,可以调整流体之间的界面张
力的大小 .
当考虑固体表面可润湿性对流体流动的影响
时 ,流体和表面相互作用力 F
σ
w
(x) 可表示为
[13]
F
σ
w
(x) = -
ρ
σ
(x)
∑
α
G
σ
w
s(x
+
e
α
)e
α
(4)
其中 s 的取值表征固相的存在 ,对于固相 ,s
=
1 ,
对于非固相 ,s
=
0 .类似地 ,G
σ
w
是控制流体和表面
相互作用力强度 的 参 数 ,对 于 D2Q9 模 型 ,当
e
α
=
0 时 ,G
σ
w
=
0 ;当 e
α
=
1 时 ,G
σ
w
=
G
w
;当
e
α
=
2 时 ,G
σ
w
=
G
w
/4 .通常情况下 ,当流体为
非润湿流体时 ,G
σ
w
取正值 ;当流体为润湿性流体
时 ,G
σ
w
取负值 .通过调节 G
σ
w
的取值 ,可以得到不
同的表面润湿性条件 .
粒子受到外力 F
σ
e
(x) 表示为
F
σ
e
(x) =
ρ
σ
g
(5)
其中
g
为单位质量力 ,如果是考虑重力的影响 ,则
g
为重力加速度 .
作用在流体粒子上作用力 F
σ
(x) 包括流体间
的相互作用力 F
σ
f
(x) 、流体和表面间的作用力
F
σ
w
(x) ,以及外力 F
σ
e
(x) .通过改变式 (2) 中平衡
态速度 u
eq
σ
来体现 F
σ
(x) 的影响 ,即
u
eq
σ
=
u
′
+
τ
σ
F
σ
/
ρ
σ
(6)
其中 u
′
定义为 u
′
=
∑
ρ
σ
u
σ
τ
σ
∑
ρ
σ
τ
σ
,这样可保证
在没有作用力时碰撞过程动量守恒
[12]
.
对于 D2Q9 模型 ,混合流体的压力为
p
(x) =
c
2
s
∑
σ
ρ
σ
(x) +
3
2
∑
σ
,
σ
G
σσ
ψ
σ
(x)
ψ
σ
(x) ,而流体的速度
定义为 u
=
∑
σ
ρ
σ
u
σ
ρ
+
∑
σ
F
σ
2
ρ
.
2 数值模拟
2 .1 研究体系及边界条件
本文的研究对象为二维水平通道内的非混相
驱替过程 .通道的长宽比为 5 ∶ 1 ,采用笛卡儿直
角坐标系 ,坐标原点位于通道的左下角 ,如图 1 所
示 .研究区域划分为 400 × 80 个格子 ,所有物理量
均为格子单位 ,与实际物理单位可以通过归一化
处理进行相互转换 .初始时 ,对于 x
≥
A ,通道内
充满 了 静 止 的 流 体 2(被 驱 替 流 体 ) ;而 对 于
0
≤
x
<
A ,通道内则为初始速度为 u 的流体
1(驱替流体) .对于通道的上表面和下表面 ,采用
半步长反弹格式以实现无滑移的边界条件 .入口
条件为恒定速度边界条件 ,出口条件为外推边界
条件 .
2 .2 界面张力的确定
伪势模型中 ,界面张力是根据 Laplace 定律
确定的 .根据 Laplace 定律 ,气泡或液滴的曲率半
径和其内外的压力差满足 Δ
p
=
γ
/R ,其中参数
γ
为界面张力 .在 100
×
100 的格子系统内 ,最初 ,正
方形的流体 2 被流体 1 所环绕 ,不考虑外力作用 ,
系统最终演化到稳定状态 .由于界面张力的作用 ,
流体 2 会呈现一个圆形液滴 .图 2 给出了不同的
黏度比 M 下的压力差 Δ
p
和曲率 1/R 的关系如图
443
大 连 理 工 大 学 学 报
第 52 卷
剩余6页未读,继续阅读
资源评论
weixin_38612648
- 粉丝: 12
- 资源: 920
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功