一组两相管流水力学模型的评估与优选修改稿.docVIP

一组油气两相管流水力学模型的评估与优选 刘 武1，张鹏1 ，程富娟2，李闽1 （1.西南石油学院石工院,四川 南充 637001；2.大港油田设计院） 摘要：现有的两相管流水力学计算模型大多数是在特定的试验条件下，得到的经验或半经验公式，应用于油气田两相流生产管线，普遍适用性较差。本文对国内、外稳态两相流软件中常用的13种水力学组合模型，在持液率计算和流型判断方面进行了改进和补充，并用JZ20-2凝析气管道的生产数据对这些水力模型进行了评估，筛选出计算精度相对较高的EF、BBM、BBME三种模型。统计误差分析和灰色关联法用于水力学模型评估的可靠性在实例应用中得到了证明，为工程设计者选择合理的水力计算模型提供了一种有效的方法。 关键词：两相管流；水力学模型；评估；优选 中图分类号：TE 进行稳态或动态条件下油气混输管线长距离输送工艺计算，首先需要选择或建立描述多相流水力特征的数学模型，公开发表的两相管流水力计算公式很多，这些模型大都是在小管径，简单实验流体（比如空气和水）和低压试验基础上建立的半经验半理论模型，应用于油气混输工艺计算中会产生很大的误差，甚至会得出错误的结论。解决这一问题较为积极的方法是根据实际生产数据，优选现有模型和关系式，或根据各个计算式的特点，组合一种新的混合模型，但其计算结果也必须经过现场生产数据的检验和修正才有应用价值。 1两相管流水力计算方法 一个完整的两相管流水力学模型应包括流型判断、持液率和压降计算三部分。流型判别是进行两相管流水力计算的第一步，一般采用流型图或根据流型转换准则，利用流体流动参数来确定流型，应用较多的是Taitel Dukler[1]（1976）、Barnea[2]（1987）和XiaoBrill[3]（1990）流型判断方法，表1中XBKS模型根据原油/伴生气、湿天然气输送特点，将两相管流的流型划分为气泡流、分层流、间歇流、环状流和弥散流五种主要形式，以满足黑油和组分模型工艺计算的需要。截面持液率是气、液两相管流最重要的特征参数之一，常常采用经验式进行计算，如Armand、Martinelli Nelson 关联式[4]，因受实验数据源的影响，误差较大；另一类计算法是采用统一的理论模型，如Wallis[5] 的一维流动计算方法和Zuber Findlay[6]考虑径向速度分布的计算模型等，这些方法建立在某一理论和假设条件的基础之上，其经验系数中纳入了不同流型形式的影响，但在相同假设条件下，没有从根本上改变采用同一种机理存在的不足，因而使用范围和计算精度也是有限的；第三类方法是依据流型的结构特点，建立相应的数学物理模型进行持液率的计算，这已成为目前气液两相管流研究的重要内容之一，具有代表性的是1990年Xiao et al[3] 统一计算模型。两相管流的压降计算是工程上最关心的问题，由于气液相之间存在质量和能量交换，气相又具有可压缩性，事实上要准确计算压降是相当困难的，现已发表的压降计算方法大体上可分为以下四种[7]： ⑴ 均相流模型压降计算公式。把气液混合物看作一种均匀连续介质，相间没有相对速度，水力摩阻系数由试验或实测数据确定，压降按单相管路计算，该模型适用于分散气泡流和弥散流。 ⑵ 分相流模型压降计算公式。把气液两相作为完全分离的两种流体，存在着不同的特性和速度，用不同的计算公式计算压降，但不考虑气液相界面间的相互作用，该模型适用于分层流和环状流，较著名的有Lockhart-Martinelli[8]和Dukler[9]压降计算法。 ⑶ 流型模型压降计算法。这种方法首先确定流型，然后根据不同的流型选择不同的计算公式，由于不同流型的能量损失机理不同，压降计算公式也不一样。典型的计算公式有Mukherjee-Brill[10]、Beggs-Brill[11]、Oliemans[12]等。 ⑷ 组合压降计算法。实际应用中，较流行的做法是针对不同的计算对象，选择不同的公式分别计算摩阻、高程和加速产生的压降，力求取得较好的计算结果。例如用Dukler[9]公式计算摩阻压降损失，高程变化引起的压降由Flanigan[13]公式进行修正，加速压降损失则由Eaton[14]公式计算。国内外两相管流稳态计算软件中常采用的组合模型，见表1。 表1常用组合水力学模型 模型及代码 流型划分 截面含液率 摩阻压降 高程压降 加速压降 Dukler-Eaton-Flanigan(DEF) 无 Eaton Dukler Flanigan Eaton Dukler-Flanigan(DF) 无 Dukler Dukler Flanigan 无 Eaton Flanigan(EF) 无 Eaton Eaton Flanigan Eaton Lock