第0绪论-1章油井流入动态与井筒多相流动计算-(2).ppt

(三)多相垂直管流压力分布计算步骤 由于多相管流中每相流体影响流动的物理参数(密度、粘度等)及混合物密度和流速都随压力和温度而变，沿程压力梯度并不是常数，因此，多相管流需要分段计算；同时，要先求得相应段的流体性质参数，然而，这些参数又是压力和温度的函数，压力却又是计算中需要求得的未知数。所以，多相管流通常采用迭代法进行计算。 有两种不同的迭代途径：按深度增量迭代和按压力增量迭代。 ⑧以计算段下端压力为起点，重复②～⑦步，计算下一段的深 度和压力，直到各段的累加深度等于管长为止。 ⑥重复②～⑤的计算，直至 。 1)按深度增量迭代的步骤 ①已知任一点(井口或井底)的压力作为起点，任选一个合适 的压力降作为计算的压力间隔?p(0.5 ～1.0MPa)。 ②估计一个对应的深度增量?h估计，计算与之对应的温度 。 ③计算该管段的平均温度及平均压力，并确定流体性质参数。 ④计算该段的压力梯度dp/dh。 ⑤计算对应于的该段管长(深度差)?h计算。 ⑦计算该段下端对应的深度及压力。 2)按压力增量迭代的步骤（略） 思考题：根据上述步骤整理出计算压力分布的程序流程框图。 说明： a. 计算压力分布过程中，温度和压力是相关的； b. 流体物性参数计算至关重要，但目前方法精度差； c. 不同的多相流计算方法差别较大，因此在实际应用中有必要根据油井的实际情况筛选精度相对高的方法。 第三节 奥齐思泽斯基（Orkiszewski）方法 奥齐思泽斯基于1967年用三大类、148口井的实际资料对前人的研究 进行了评价，加上自己的研究，提出了此方法。其主要构成为： Orkiszewski流型图 他提出的四种流动型态是：泡流、段塞流、过渡流及环雾流 存容比（滞流率）：多相流动的某一管段中某相流体的体积与管段容积之比。（持液率） 出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为： 一、压力降公式及流动型态划分界限 由垂直管流能量方程可知，压力降是摩擦能量损失、势能变化和动能变化之和： 所以压降计算式为： 表1-3 Orkiszewski方法流型划分界限 不同流动型态下 和 的计算方法不同。 二、平均密度及摩擦损失梯度的计算 气相存容比(截面含气率、空隙率)Hg ：管段中气相体积与管段容积之比值，也等于fg /f 。 液相存容比(截面含液率、持液率)HL ：管段中液相体积与管段容积之比值，也等于fl /f 。 (1)泡流 平均密度: 滑脱速度：气相流速与液相流速之差。 则： 泡流摩擦损失梯度按液相进行计算： 图1-21 图1-21 摩擦阻力系数曲线 (2)段塞流 平均密度: 滑脱速度 vs 液体分布系数 滑脱速度 由泡雷诺数 查图确定 ～ 曲线 根据泡雷诺数 及雷诺数 查图确定 ～ 曲线 (1)假价设一个vs值，求得C1及C2 (2)用式 计算一个vs 值 (3)重复计算直到假设值与计算值接近为止 vs值的确定 迭代法 vs也可由公式进行计算 的计算 值需根据连续液相的类别 及气液总流速来选用计算公式 计算式详见教材公式(1-58) a ~e 时 时 计算得的 必须满足下面的条件： (3)过渡流 过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。 段塞流的摩擦梯度 式中的摩擦系数 f ， 根据管壁相对粗糙度和雷诺数由穆迪图查得。 雾流混合物平均密度计算公式与泡流相同： 由于雾流的气液无相对运动速度，即滑脱速度接近于零，基本上没有滑脱。 雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度查图得。 摩擦梯度: (4)雾流 所以： 图1-24 Orkiszewski方法计算流程框图 第四节贝格斯-布里尔方法 Beggs-Brill方法是可用于水平、垂直和任意倾斜气液两相管流动计算的方法，是1973年，Beggs和Brill根据在长15m，直径1inch和1.5inch聚炳烯管中，用空气和水进行实验的基础上提出的。 实验参数范围 气体流量 0～0.098 m3/s ； 液体流量 0～0.0019 m3/s ； 持液率 0～0.87 m3/m3； 系统压力 241～655 kpa(绝对压力)； 压力梯度 0～18 kPa/m； 倾斜度 -90