多相管流-第三章垂直气液两相管流计算1

第三章 垂直气液两相管流压力梯度计算模型及方法

垂直气液两相管流压力梯度计算模型及方法•3.13.1流动模型

流动模型•3.23.2压力分布计算方法压力分布计算方法•

3.33.3垂直气液两相管流压力梯度计算模型及垂直气液两相管流压力梯度计算模型及方法

•3.43.4水平或倾斜管中气液两相流动计算模型水平或倾斜管中气液两相流动计算模型及方法

•3.53.5水平管中气体和非牛顿液体的两相流动水平管中气体和非牛顿液体的两相流动主要内容

第一节流动型态

流动型态的划分方法：两类

第一类划分方法：

根据两相介质分布的外形划分

泡状流、弹状流或团状流、（层状流、波状流）、段塞流或冲击流、环状流、雾状流

第二类划分方法：

按流动的数学模型或流体的分散程度划分

分散流、间歇流、分离流

分散流、间歇流、分离流

2010-3-26

3

垂直气液两相流流型

水平气液两相流流型

两种分类方法比较

第一类划分方法较为直观

第二类划分方法便于进行数学处理

第

一类划

分方法•泡状流

•弹状流或团状流

•层状流

•波状流

•段塞流或冲击流

•环状流

•雾状流

第

二

类

划

分

方

法

•分散流

•间歇流

•分离流

•分离流

•间歇流

•分离流

•分散流

两类划分结果

的对应关系

2010-3-264

垂直环空两相流型

第二节 压力分布计算方法

第二节

•由于多相管流中每相流体影响流动的物理参数（密度、粘度等）及混合物密度和流速都随压力和温度而变，沿程压力梯度并不是常数。因此，多相管流压降需要分段计算，并要预先求得相应段的流体性质参数。然而，这些参数又是压力和温度的函数，压力却又是计算中需要求得的未知数。所以，多相管流通常采用迭代法进行计算。

一、常用两相流压降计算方法

•早期均匀流方法（总摩阻系数法）

1952 Poettmann 1952 Poettmann—

—Carpenter 80 80’

’s 陈家琅 λ'~(N Re )2•经验相关式

1963 Duns--Ros 无因次化处理 N vL 、Nvg 、N D 、N L 1965 Hagedorm--Brow 现场实验 1967 Orkiiszewski 流型组合 1973 Beggs--Brill 倾斜管实验

1985 Mukherijee--Brill 改进实验条件•现代机理模型

SPE20630等考虑具体流型的物理现象

第二节第二节

压力分布计算方法

段塞流示意图环状流示意图2010-3-268

2010-3-26

当 单相液流， H L 、ρm 、f m 随两相流流型变化

b P P ≥

二、两相管流压降计算

根据地面条件应用关系式计算井底流压1 输入数据

油管数据：管长L 、管径D 、井斜角θ、粗糙度e 油气井产量：油气水日产量Q O 、Q SC 、Q W

或Q L 、f w 、GOR P (GLR P ) Q

Q W = f w Q L Q O = Q L - Q w

Q SC = GOR P Q O 或 Q SC = GLR P Q L 边界条件：井口压力P wh 、井口温度T wh 、地温梯度g t 考虑井温线性分布 T(Z)=T wh +g t Z 油气水相对密度γo 、γg、γw

第二节第二节

压力分布计算方法

2 输入数据单位处理

常用单位 统一单位

Q —m 3/d q —m 3/s μ—Pa.s P —MPa

p —Pa V Pa V—

—m/s D —mm d mm d——m T m T——℃ T T—

—K 第二节第二节

压力分布计算方法

3 输入流体物性资料

气：拟临界压力、温度 Pc ， Tc

偏差系数 Zg(Pr, Tr)

粘度μg

油：μo ， 溶解油气比 Rs

体积系数 Bo ， 油气界面张力σo

水：μw ，σw ，B w

第二节第二节

压力分布计算方法

)

,(2

23004hk p h z F k ++=第二节第二节

压力分布计算方法

龙格库塔数值解法

•压力梯度函数F(Z，P)计算步骤

(1) Z处流动温度 T(Z)=T0+g t Z

(2) 计算T、P条件下的有关物性

(3) 气液体积流量 q g，q L

(4) 气液表观流速V sg、V sL和V m

(5) 计算λL、μL、ρns、μns

(6) 无因次量N Rens、N L、N gV、N LV、N gvsM

(7) 计算H L、ρm

(8) 判别流型，计算f m

(9) 计算F(Z，P)

•2、迭代计算

第二节第二节

压力分布计算方法误差又能提高计算速度。

•已知井口（Z 0＝0）压力p 0=p wh 沿油管的压力分布计算步骤如下：

• 1. 记计算节点序号i = 1， 选取压力增量Δp 和对应的管长初值ΔZ 0；

• 2. 计算第i 节点位置Z i 及其温度

• Z i =Z i-1+ΔZ 0

•考虑流体温度沿井深线性变化，节点处的温度为• T i =T(Z i )=T 0+g T Z i

• 式中 g T ——温度梯度，K/m (℃/m)；• T 0——井口流动温度，K。

第二节第二节

压力分布计算方法