石油工程1009高起专采油工程设计(最新整理by阿拉蕾)

石油工程1009高起专采油工程设计(最新整理by阿拉蕾)内容摘要：

???? iioi PPP 1 ⑧ 以 iL 处的压力 iP 为起点压力重复第 ② ～ ⑦ 步，计算下一段的 11 深度 1?iL 和压力 1?iP ，直到各段累加深度等于或大于管长 L 时为止。 计算气 液两相垂直管流的 Orkiszewski 方法 本设计井筒多相流计算采用 Orkiszewski 方法。 Orkiszewski 法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。 如图 1 所示。 在处理过渡性流型时，采用内插法。 在计算段塞流压力梯度时要考虑气相与液体的分布关系。 针对每种流动型态提出了存容比及摩擦损失的计算方法。 图 1 气液混合物流动型态 (Orkiszewski) 由前面垂直管流能量方程可知，其压力降是摩擦能量损失、势能变化和动能变化之和。 由式 (236)可直接写出多项垂直管流的压力降公式： mmmmf dvvdhgdhdP ??????? (26) 式中 P—压力， Pa； f? —摩擦损失梯度， Pa/m； h —深度， m； g —重力加速度， m/s2； 12 m?—混合物密度， kg/m3； mv —混合物流速， m/s。 动能项只是在雾流情况下才有明显的意义。 出现雾流时，气体体积流量远大于液体体积流量。 根据气体定律，动能变化可表示为： dpPA qWdvv p gtmmm 2??? (27) 式中 pA —管子流通截面积， m2； tW —流体总质量流量， kg/s； gq —气体体积流量， m3/s。 将式 (27)代入式 (26)，并取 khdh ??? ， kpdP?? ， mm ??? ， PP? 经过整理后可得： ??kPkpgtfm hPAqWg ??? ]1[ 2?? (28) 式中 kP? —计算管段压力降， Pa； kh? —计算管段 的深度差， m； P —计算管段的平均压力， Pa。 不同流动型态下的 m? 和 f? 的计算方法不同 ,下面按流型分别介绍。 (1)泡流 平均密度 ggLLm HH ????? ? ? ggLg HH ????? 1 1?? gL HH 式中 gH —气相存容比 (含气率 )，计算管段中气相体积与管段容积之比值； 13 LH —液相存容比 (持液率 )，计算管段中液相体积与管段容积之比值； mg ??? 、 Ｌ —在 TP、 下气、液和混合物的密度， kg/m3。 气相存容比由滑脱速度 sV 来计算。 滑脱速度定义为：气相流速与液相流速之差。 )1(1 gp gtgp ggsLgsgs HA qqHA qHvHvv ??????? 可解出 gH ： Hg?psgpstpst Av qAv qAv q 4)1(1[21 2 ???? 式中 sv—滑脱速度，由实验确定， m/s； sgv 、 sLv —气相和液相的表观流速， m/s。 泡流摩擦损失梯度按液相进行计算： 22LHLt vDf ??? )1( gp LLH HA qv ?? 式中 f —摩擦阻力系数 ； LHv —液相真实流速， m/s。 摩擦阻力系数 f 可根据管壁相对粗造度 D/? 和液相雷诺数 ReN查图 2。 液相雷诺数 : L LsLDvN ???Re 式中 L? —在 TP、 下的液体粘度，油、水混合物在未乳化的情况下可取其体积加权平均值，。 14 图 2 (2)段塞流混合物平均密度 Lpst psLtm Avq AvW ???????? (34) 式中 ?—液体分布系数； sv —滑脱速度， m/s。 滑脱速度可用 Griffith 和 Wallis 提出的公式计算： gDCCvs 21? (35) 15 (3)过渡流 过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。 MiSM sgSLSM gMm LL LvLL vL ????????? (36) MiSM ggSLSM gMt LL LvLL vL ????????? (37) 式中的 SL? 、 SL? 及 Mi? 、 Mi? 为分别按段塞流和雾流计算的混合物密度及摩擦梯度。 (4)雾流 雾流混合物密度计算公式与泡流相同： ggLgggLLm HHHH ?????????? )1( 由于雾 流的气液无相对运动速度，即滑脱速度接近于雾，基本上没有滑脱。 所以 gLgg qqqH ?? (38) 摩擦梯度则按连续的气相进行计算，即 Dvf sggf 2 2??? (39) 式中 sgv —气体表观流速， pgsg Aqv /? ， m/s。 雾流摩擦系数可根据气体雷诺数 ? ?gNRe和液膜相对粗糙度由图 2查得。 按不同流动型态计算压力梯度的步骤与前面介绍的用摩擦损失系数法基本相同，只是在计算混合物密度及摩擦之前需要根据流动型态界限确定其流动型态。 图 3 为 Orkiszewski 方法的计算流程框图。 16 以井口油压或井底流压为起点，选择合适的压力间隔P?，假设 h?计算平均 P 和 T ，并求得在此 P 和 T 下的流体性质参数和流动参数，以及相应的流动型态界限 LB、 Lg和 LM确定流动型态雾流计算气相存容比、平均密度及摩擦梯度过渡流分别按段塞流和雾流计算平均密度及摩擦梯度，并进行内插段塞流计算滑脱速度、液体分布系数、平均密度和摩擦梯度泡流计算气相存容比、平均密度和摩擦梯度计算并比较h?，重复上述计算使h?的计算值与假设值相等或在允许的误差范围。