采油工程原理与设计 第0绪论－1章油井流入动态与井筒多相流动计算.ppt

管英柱 男 吉林白山人 1992.9～1996.6：江汉石油学院开发系采油工程专业； 1996.7～2001.8：中石油吉林油田前大采油厂； 2001.9～2007.6：成都理工大学能源学院油气田开发工程专业（硕博）； 2007.7～ 长江大学石油工程学院 主讲：采油工程、储层改造 联系方式： Tel E-mail:guanyzh@ QQ 答疑地点：9教-100号 第一章 油井流入动态与井筒多相流计算 第一节 油井流入动态（IPR曲线） 油井流入动态： 油井产量(qo) 与井底流动压力(pwf) 的关系，反映了油藏向该井供油的能力。 图1-1 典型的流入动态曲线 单相流动时，油层物性及流体性质基本不随压力 变化。 图2-4 计算的溶解气驱油藏油井IPR曲线1-用测试点按直线外推；2-计算机计算值；3-用Vogel方程计算值 按Vogel方程计算的IPR曲线，最大误差出现在用小生产压差下的测试资料来预测最大产量时，但一般误差低于5％。虽然随着采出程度的增加到开采末期误差上升到20％左右，但其绝对值却很小。 如果用测试点的资料按直线外推，最大误差可达 70～80％，只是在开采末期约30%。 采出程度 Np 对油井流入动态影响大，而kh/μ、Bo、k、So等对其影响不大。 射孔完成——打开性质不完善； 未全部钻穿油层——打开程度不完善； 打开程度和打开性质双重不完善； standing方法计算不完善井IPR曲线的步骤： 图1-8 拟合的IPR曲线与实际曲线的对比 _____拟合的IPR曲线,……实际曲线 图1-9 参数v、n与采出程度之间的关系 (一) 采液指数计算(由测试点确定曲线) 第二节 井筒气液两相流基本概念 一、井筒气液两相流动的特性 奥齐思泽斯基（Orkiszewski）方法 第四节贝格斯-布里尔方法 图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； 倾斜度 -900～+900； 流 型 水平管流动的全部流型。 1）绘制IPR曲线 解： pwf 11.2 10.2 9.7 9.1 q 16.1 35.9 46.1 57.9 平均地层压力 2）计算采油指数 3）查表得 4）直线外推至q=0,求得 5）Pwf=8.8MPa时 Qo=20×(12-8.8)=64m3 已知D井平均地层压力为16MPa，Pb为13MPa,Pwf为8MPa时的产量为71.45m3/d，试计算Pwf为14MPa和7MPa时的产量并绘制该井的 IPR曲线。 例 题 a.计算