垂直管流中气—液两相的研究现状.doc

垂直管流中气—液两相的研究现状 摘要：深井气举采油由于流程长、产量变化范围大而有着一系列的特点。因流程长，使得气液两相混合物在流动过程中压力变化很大，不可避免地要产生两种甚至多种流动型态（Flow Pattern），而不同流动型态对压降、持液率（空隙率）、混合物密度影响非常大，影响着气举井的设计和运行。有人认为，在深井、超深井气举中，有1000m左右的管段可形成段塞流。另外产量变化也极大地影响着压力梯度和持液率，而压力梯度是气举井设计的最基本依据（以此确定注气点位置、气举阀分布及数量、注气压力），因此，要使气举井工作最优，必须研究铅直气液两相管流问题。 关键词：垂直管流 流型 两相 １ 研究现状 对于气－ 液两相管流的研究， 一般多是从能量平衡和物质守恒关系出发， 来计算气液混合物在管中的平均流速、密度、水头损失、压力梯度等有关参数问题。早在 1914 年Wisconsin大学的Davis-Weider便发表了在直径为31.75mm短玻璃管内以空气举升水的大量实验数据。他们把管内摩擦因数与滑移关联起来， 试图得到以DVP为自变量，与单相流摩擦因数f相一致的曲线，由于当时并未考虑混合物总的流动密度，只是使用了水的密度，而没有达到预期的效果，但却为后来的两相流研究奠定了基础。经过近一个世纪国内外学者的大量研究工作，已经预测流动模型形成和发展了各种各样的力学模型。研究方法上可概括为不考虑流动型态和考虑流动型态两类。 1.1不考虑流动型态 方法是将流型、滑脱损失及加速度影响计入两相阻力系数中，根据现场和试验数据得出压力梯度相关式。在早期Poettmann-Carpenter（1952）和M.R.Tek(1961)的摩擦损失系数法的基础之上，Baxendell、Thomax(1961)、Brown(1963)、Hagedorn(1965)及Govier、Aziz(1972)和大庆石油学院陈家琅(1979)等国内外学者先后针对井内气液比、密度和粘度等参数随井深变化、油井总流量、以及油管直径的影响等，在计算方法上作了进一步的修正， 为推动后来的各种经验和半经验压降预测模型的建立和确保模型的精确度， 起了重要的指导作用。 1.2考虑流动型态 这种方法的共同特点是以某一标准划分流动型态， 然后根据不同流动型态确定气液混合物的密度、空隙率及压力梯度。Duns-Ros两位学者于1961年至1963年先后通过数千次的实验，从10个无因次群中确定出了4个有意义的无因次群，首次得到了真正意义上的流动形态分布图；后来，针对前人提出的流态分布图，不同的学者通过总结前人的不足，形成了各自关于流态划分、流态判别以及流态转换的解释模型，先后发展起来了 Orkiszewaki方法、Aziz-Gorier-Fogaras方法、Taite-Bornea-Dukler方法、Hasan-Kabir方法、Ansari等的方法， 以及现在的 Kaya方法、H.ShiJ.A.Holmes等的井筒三相稳定流漂移流动模型和E.A.Osman等提出的基于反向传播学习算法的人工神经网络(ANN)方法(表1) 表 １ 气－ 液两相管流研究方法 方 式 模型（ 方法） 计算条件 评价及应用 不 考 虑 流 动 型 态 Poettmann-Carpenter摩擦损失系数法 气液两相混合物视为均匀介质 20 世纪 60 年代以前最经典垂直多相管流模型， 适用于高流量低气液比油井 Baxende及M.R.Tek和 Thomax等摩擦损失系数法 考虑液体粘度和气液质量比的气液两相雷诺数 是对模型在考虑气液比、密度及粘度等因素影响之后的修正， 在低产量低密度油井范围内仍存在误差 Hagedorn- Brown方法 以气液两相滑脱 为 基 础（ 滑移模型） 通过计算持液率来修正摩擦因数，从而求解由滑移增大的静压梯度，在水平多相管流中得到了广泛应用 考 虑 流 动 型 态 Duns-Ros压降计算方法 利用因次分析，从 10 个无因次群中确定出 4 个有意义的无因次群 提出的无因次群对油井多相垂直管流方法的发展起了重要作用，在工程上可以达到很好的精度，它更适用于较短的管段， 而对深度或压差很大的井， 需进行一连串的分段计算 Orkiszewaki方法 将压力梯度按流动型态划分，采用不同的方法进行计算 率先对每个流型单独进行计算并定义了液体分布系数，使Griffith-Walks的段塞流计算延伸到高流量范围，是业界公认的可靠方法之一 Aziz-Gorier-Fogaras方法 在密度和摩擦损失项中，通过气液两相分离作用，引入当地气相体积因素 这种分布图流型转变界线明确，有表达式，计算机处理方便，是业界普遍认可的方法之一 Taite-Bornea-Dukl