抽油机井系统效率分析及提高措施1

1、薛建泉薛建泉电电 话：话：83937778393777234567%100PP入水抽油机输入功率抽油机有效功率8第二节第二节 系统效率计算与测试系统效率计算与测试 井地入光光水入水PPPPPP9第二节第二节 系统效率计算与测试系统效率计算与测试 321KP/P入光地10343232121PPPPPPP1电动机输入功率，即电动机输入功率，即P ，kW；P2电动机输出功率，电动机输出功率，kW；P3减速箱输出功率，减速箱输出功率，kW； P4光杆功率，即光杆功率，即P光光，kW。第二节第二节 系统效率计算与测试系统效率计算与测试 11第二节第二节 系统效率计算与测试系统效率计算与测试 7654P/

2、P光水井12第二节第二节 系统效率计算与测试系统效率计算与测试 787676565454PPPPPPPP13第二节第二节 系统效率计算与测试系统效率计算与测试 7654321K井地14 电动机的测试参量有：电动机的测试参量有：有功功率、无功功率、电流、电压、电有功功率、无功功率、电流、电压、电机转速、电机输出扭矩等。机转速、电机输出扭矩等。皮带、减速箱应当分别测试，从而分别计算它们的效率。但由皮带、减速箱应当分别测试，从而分别计算它们的效率。但由于抽油机的结构原因，单独测试困难，将它们作为一个组件于抽油机的结构原因，单独测试困难，将它们作为一个组件测试测试参量有：参量有：输出轴扭矩、输出轴转速

3、输出轴扭矩、输出轴转速第二节第二节 系统效率计算与测试系统效率计算与测试 15第二节第二节 系统效率计算与测试系统效率计算与测试 光杆的测试参数有：光杆的测试参数有：盘根盒的测试参数有：盘根盒的测试参数有： 抽油杆的测试参数为：抽油杆的测试参数为：16第二节第二节 系统效率计算与测试系统效率计算与测试 抽油泵的测试参数有：抽油泵的测试参数有： 井下管柱的测试参数有：井下管柱的测试参数有：17第二节第二节 系统效率计算与测试系统效率计算与测试 目前虽有一些计算有杆抽油系统各部分功耗的公式，但其一目前虽有一些计算有杆抽油系统各部分功耗的公式，但其一般误差都很大。只有采用一定的测试仪器，在现场进行实

4、测，般误差都很大。只有采用一定的测试仪器，在现场进行实测，才可能得到比较精确的耗功数据，进而计算系统效率。才可能得到比较精确的耗功数据，进而计算系统效率。 为了测试有杆抽油系统的分效率，可在整个系统布点为了测试有杆抽油系统的分效率，可在整个系统布点9 9处。处。18第二节第二节 系统效率计算与测试系统效率计算与测试 分分级级效效率率测测试试布布点点简简图图19第二节第二节 系统效率计算与测试系统效率计算与测试 位于电动机动力输入处，在该点测各种电参数位于电动机动力输入处，在该点测各种电参数用电功率表可以在用电功率表可以在1 1点测得电动机的输入功率。除此之外还要点测得电动机的输入功率。除此之外

5、还要测输入的电流、电压、无功功率（功率因数测输入的电流、电压、无功功率（功率因数coscos）等。用最）等。用最大电流比可以确定抽油机的平衡度；根据电压是否变化，可以大电流比可以确定抽油机的平衡度；根据电压是否变化，可以了解整个系统的平稳程度；无功功率或功率因数则说明所测抽了解整个系统的平稳程度；无功功率或功率因数则说明所测抽油机系统对整个供电系统的影响。油机系统对整个供电系统的影响。测点测点2 2位于电动机输出轴端，测点位于电动机输出轴端，测点3 3位于减速位于减速箱的曲柄轴输出端。用动态应变仪分别测箱的曲柄轴输出端。用动态应变仪分别测2 2、3 3点处的扭矩，同点处的扭矩，同时用转速计测该

6、两点处的转速；然后就可以分别求出电动机的时用转速计测该两点处的转速；然后就可以分别求出电动机的输出功率和减速箱输出功率。输出功率和减速箱输出功率。 20第二节第二节 系统效率计算与测试系统效率计算与测试 （3 3）测点测点4 4 位于悬绳器处。用动力仪可测得示功图；根据示位于悬绳器处。用动力仪可测得示功图；根据示功图面积及动力仪的比例系数可求得光杆功率。功图面积及动力仪的比例系数可求得光杆功率。（4 4）测点测点5 5 在井下抽油泵处。在泵的吸入口和排出口分别安在井下抽油泵处。在泵的吸入口和排出口分别安装振弦压力计，测试抽油泵吸入和排出压力。装振弦压力计，测试抽油泵吸入和排出压力。（5 5）测

7、点测点6 6、7 7、8 8、9 9 四点都在井口处。在四点都在井口处。在6 6点处用回声仪点处用回声仪测油井动液面深度，在测油井动液面深度，在7 7点处用压力表测套管压力，在点处用压力表测套管压力，在8 8点处用点处用压力表测油管压力，压力表测油管压力，9 9点为油管出口处，在该点用流量计测量点为油管出口处，在该点用流量计测量油井产量油井产量Q Q。21 根据根据5 5、6 6、7 7、8 8、9 9五点所测得数值，可分别求得抽油泵输五点所测得数值，可分别求得抽油泵输出功率和系统的有效功率，出功率和系统的有效功率，9 9个测试点可确定个测试点可确定6 6种功率（种功率（P P1 1、P P2

8、 2、P P3 3、P P4 4、P P5 5、P P6 6）。用这）。用这6 6种功率求全部的种功率求全部的7 7个分效率是不够的，个分效率是不够的，因为井口盘根盒的耗功和井下抽油杆耗功还是个未知数。因为井口盘根盒的耗功和井下抽油杆耗功还是个未知数。 要在直接测试要在直接测试P P5 5与与P P6 6有较大的困难，为此采用模拟测试方法。有较大的困难，为此采用模拟测试方法。 测试盘根盒耗功测试是室内完成的。用一根短管把两个盘根盒测试盘根盒耗功测试是室内完成的。用一根短管把两个盘根盒连接起来，管上安装有压力管通油泵，用油泵对盘根盒加压。光连接起来，管上安装有压力管通油泵，用油泵对盘根盒加压。光

9、杆上贴有电阻片测出光杆的摩擦力。见下图杆上贴有电阻片测出光杆的摩擦力。见下图第二节第二节 系统效率计算与测试系统效率计算与测试 22 盘根盒耗功测试装置简图盘根盒耗功测试装置简图1盘根盒；盘根盒；2连接管；连接管；3光杆；光杆；4加压管；加压管；5压力表压力表第二节第二节 系统效率计算与测试系统效率计算与测试 23 在研究分析抽油机系统各部分耗功时，进行分级效率测试在研究分析抽油机系统各部分耗功时，进行分级效率测试是必要的，因为只有这样才能深入地揭示系统中各部分的工况。是必要的，因为只有这样才能深入地揭示系统中各部分的工况。分级效率测试与计算的工作量很大，在日常管理工作中要对每一分级效率测试与

10、计算的工作量很大，在日常管理工作中要对每一口抽油机井都进行这样的测试工作是不可能的。为此，下面给出口抽油机井都进行这样的测试工作是不可能的。为此，下面给出一种简便而适应于日常管理工作的测试方法。一种简便而适应于日常管理工作的测试方法。（1 1）测输入功率）测输入功率P P及电流及电流I I用功率表测电动机输入功率。用钳形电流表测量抽油机上下冲用功率表测电动机输入功率。用钳形电流表测量抽油机上下冲程时最大的工作电流，从而确定抽油机平衡程度。程时最大的工作电流，从而确定抽油机平衡程度。第二节第二节 系统效率计算与测试系统效率计算与测试 24式中式中 I I下下maxmax抽油机下冲程的最大工作电流

11、，抽油机下冲程的最大工作电流，A A； I I上上maxmax抽油机上冲程的最大工作电流，抽油机上冲程的最大工作电流，A A； nn测量组数；测量组数； BB平衡度，平衡度，% %。%100nI/IBn1imaxmax上下第二节第二节 系统效率计算与测试系统效率计算与测试 25（2 2）测光杆功率）测光杆功率P P光光及水功率及水功率P P水水 通过动力仪测光杆功率；用回声仪测动液面深度；由压力表测通过动力仪测光杆功率；用回声仪测动液面深度；由压力表测套管压力、油管压力；在计量站测出油井产液量，从而求出抽油套管压力、油管压力；在计量站测出油井产液量，从而求出抽油机系统的有效功率。机系统的有效功

12、率。 求得求得P P、P P光光、P P水水后，可分别求得后，可分别求得、地地及及井井。这样的测。这样的测试工作，工作量不大，又可求得系统总效率以及地面与井下的分试工作，工作量不大，又可求得系统总效率以及地面与井下的分效率，对于生产井的管理工作有指导意义。效率，对于生产井的管理工作有指导意义。 第二节第二节 系统效率计算与测试系统效率计算与测试 26第二节第二节 系统效率计算与测试系统效率计算与测试 2728油井流入动态曲线（油井流入动态曲线（IPRIPR曲线）曲线）： 表示产量与井底流压关系的曲线，简称表示产量与井底流压关系的曲线，简称IPRIPR曲线。曲线。油井流入动态：油井流入动态： 油

13、井产量与井底流动压力的关系。油井产量与井底流动压力的关系。 它反映了油藏向井的供油能力，反映油藏压力、油层物性、它反映了油藏向井的供油能力，反映油藏压力、油层物性、流体物性、完井质量等对油层渗流规律的影响，是采油工程与油流体物性、完井质量等对油层渗流规律的影响，是采油工程与油藏工程的衔接点。通过油井流入动态研究为油藏工程提供检验资藏工程的衔接点。通过油井流入动态研究为油藏工程提供检验资料；为采油工程的下一步工作提供依据；检查钻井、固井、完井料；为采油工程的下一步工作提供依据；检查钻井、固井、完井和各项工艺措施等技术水平的优劣。和各项工艺措施等技术水平的优劣。第一节第一节 油井流入动态油井流入动

14、态291.1.单相液体流入动态单相液体流入动态asXBPPhkqoowfroo21ln)(2供给边缘压力不变圆形地层供给边缘压力不变圆形地层中心一口井的产量公式为：中心一口井的产量公式为：asXBPPhkqoowfroo43ln)(2圆形封闭油藏圆形封闭油藏，拟稳态条件下的油井产量公式为：，拟稳态条件下的油井产量公式为：第一节第一节 油井流入动态油井流入动态30 单相流动时，油层物性及流体性质基本不随压力变化，单相流动时，油层物性及流体性质基本不随压力变化，产量公式可表示为产量公式可表示为: :)(wfroPPJq采油指数采油指数J J的获得：的获得： 油藏参数计算油藏参数计算 试井资料：测得

15、试井资料：测得3 35 5个稳定工作制度下的产量及其流个稳定工作制度下的产量及其流压，便可绘制该井的实测压，便可绘制该井的实测IPRIPR曲线曲线第一节第一节 油井流入动态油井流入动态312.2.油气两相渗流时的流入动态油气两相渗流时的流入动态平面径向流，直井油气两相渗流时油井产量公式为：平面径向流，直井油气两相渗流时油井产量公式为：(1)(1)垂直井油气两相渗流时的流入动态垂直井油气两相渗流时的流入动态drdpBhrkqoooo2kkkorodpBKrrkhqewfPPooroweoln2第一节第一节 油井流入动态油井流入动态321)Vogel 1)Vogel 方法方法假设条件：假设条件：a

16、.a.圆形封闭油藏，油井位于中心；圆形封闭油藏，油井位于中心；b.b.均质油层，含水饱和度恒定；均质油层，含水饱和度恒定；c.c.忽略重力影响；忽略重力影响；d.d.忽略岩石和水的压缩性；忽略岩石和水的压缩性； e.e.油、气组成及平衡不变；油、气组成及平衡不变；f.f.油、气两相的压力相同；油、气两相的压力相同；g.g.拟稳态下流动状态。拟稳态下流动状态。2max8 . 02 . 01rwfrwfooPPPPqqVogelVogel方程方程第一节第一节 油井流入动态油井流入动态33dpBksrrkhqrwfPPooroweo43ln22 2）费特柯维奇方法）费特柯维奇方法溶解气驱油藏：溶解气

17、驱油藏：假设假设 与压力与压力 成直线关系，则：成直线关系，则：ooroBkprwfPPweocpdpsrrkhq43ln222243ln2wfrwePPcsrrkhcPBkooro第一节第一节 油井流入动态油井流入动态34rpoororKPc)(1式中：式中：rwfrPooroweoPPPKsrrkhqr243ln222则：则：第一节第一节 油井流入动态油井流入动态352PCPBAq (2)(2)斜井和水平井的斜井和水平井的IPRIPR曲线曲线 Cheng Cheng对溶解气驱油藏中斜井和水平井进行了数值模拟，对溶解气驱油藏中斜井和水平井进行了数值模拟，并用回归的方法得到了类似并用回归的方法

18、得到了类似VogelVogel方程的不同井斜角井的方程的不同井斜角井的IPRIPR回归方程：回归方程：第一节第一节 油井流入动态油井流入动态36 Bendakhlia Bendakhlia等用两种三维三相黑油模拟器研究了多种等用两种三维三相黑油模拟器研究了多种情况下溶解气驱油藏中水平井的流入动态关系。得到了不情况下溶解气驱油藏中水平井的流入动态关系。得到了不同条件下同条件下IPRIPR曲线。曲线。nrwfrwfooPPvPPvqq2max11BendakhliaBendakhlia用公式来拟合用公式来拟合IPRIPR曲线图版，发现吻合很好。曲线图版，发现吻合很好。第一节第一节 油井流入动态油井

19、流入动态37dpBKrrkhqrwfPPooroweo43ln2（1 1）基本公式）基本公式 当油藏压力高于饱和压力，而流动压力低于饱和压力当油藏压力高于饱和压力，而流动压力低于饱和压力时，油藏中将同时存在单相和两相流动，拟稳态条件下产时，油藏中将同时存在单相和两相流动，拟稳态条件下产量的一般表达式为量的一般表达式为: :3. 3. 时的流入动态时的流入动态wfbrPPP第一节第一节 油井流入动态油井流入动态38组合型IPR曲线（2）实用计算方法第一节第一节 油井流入动态油井流入动态39)(brbPPJq流压等于饱和压力时的产量为：流压等于饱和压力时的产量为： 当当 时，由于油藏中全部为单相液

20、体流动。时，由于油藏中全部为单相液体流动。bwfPP)(wfroPPJq流入动态公式为：流入动态公式为：第一节第一节 油井流入动态油井流入动态40)(8 . 02 . 01 2bwfbwfcboPPPPqqq) 1(8 . 18 . 1brbbcPPqJPq)(8 . 0)(2 . 01 8 . 12bwfbwfbbroPPPPPPPqJ 当当 后，油藏中出现两相流动。后，油藏中出现两相流动。bwfPP直线段直线段采油指数采油指数第一节第一节 油井流入动态油井流入动态413.3.油气水三相油气水三相IPR IPR 曲线（曲线（PetrobrasPetrobras方法）方法） 油气水三相IPR

21、曲线第一节第一节 油井流入动态油井流入动态42已知一个测试点已知一个测试点( ( 、 ) )和饱和压力和饱和压力 及油藏压力及油藏压力wftestpttestqbPrP 当当 时：时：bwftestPPwftestrttestlPPqJ 当当 时：时：)(8 . 1)1 (wftestrwbbrwttestlPPfAPPPfqJbwftestwfGPPP 当当 时：时：wfGwftestPP0bwrwbwftestwttestlpfpfppfqJ49849828 . 02 . 01bwftestbwftestPPPPA其中：其中：bwwfGpfP94直线段直线段采液指数采液指数第一节第一节 油

22、井流入动态油井流入动态43 ，则：，则：btqq 0ltrwfJqPP ，则按流压加权平均进行推导：，则按流压加权平均进行推导：maxotbqqqwfwaterwwfoilwwfPfPfP)1 (bobtbwfoilqqqqPPmax80811125. 0ltrwfwaterJqPp (2) (2)某一产量下的流压计算某一产量下的流压计算)(80811)1 (125. 0maxbobtbwltrwwfqqqqPfJqPfP所以：所以：第一节第一节 油井流入动态油井流入动态44lwqqtwfJfdqdPot98max因为：因为：lwotlorwwfJfqqJqPfP)98)()(maxmax所以

23、：所以：若若 ，则综合，则综合IPRIPR曲线的斜率可近似曲线的斜率可近似为常数。为常数。maxmaxttoqqq第一节第一节 油井流入动态油井流入动态45第一节第一节 油井流入动态油井流入动态46第二节第二节 已知动液面计算井底流压方法已知动液面计算井底流压方法 由动液面计算井底流压的方法是在套管不产气正常生产井由动液面计算井底流压的方法是在套管不产气正常生产井井筒中，将井筒流体划分为气柱、油柱及液柱三段，分别求出井筒中，将井筒流体划分为气柱、油柱及液柱三段，分别求出各段的压力降落，从而计算井底流压。其计算表达式如下：各段的压力降落，从而计算井底流压。其计算表达式如下： loglogcwfP

24、PPPPPPP 471.1.动液面处压力的计算动液面处压力的计算fgHTZcgePP03416. 02. 2. 油柱所产生压力的计算油柱所产生压力的计算 nifiiniioinioioHHhhgPP161110第二节第二节 已知动液面计算井底流压方法已知动液面计算井底流压方法 483.3.井底流压的计算井底流压的计算 当油管未下至油层中部时，在油管管鞋至当油管未下至油层中部时，在油管管鞋至油层中部深度处流体呈现多相混合流动。已知油层中部深度处流体呈现多相混合流动。已知油管管鞋处的压力、温度等，并以此作为起点，油管管鞋处的压力、温度等，并以此作为起点，利用多相流相关式计算井底流压。利用多相流相关

25、式计算井底流压。第二节第二节 已知动液面计算井底流压方法已知动液面计算井底流压方法 49 根据流体静力学的原理，关井后井底静压是井口根据流体静力学的原理，关井后井底静压是井口套压与各段流体压力降落之和，即套压与各段流体压力降落之和，即wogwogcsPPPPPPPP 第三节第三节 根据静液面计算静压方法根据静液面计算静压方法 501.1.井筒中油水界面深度的计算井筒中油水界面深度的计算 关井之后由地层续流到井筒中的流体含水率与关井之后由地层续流到井筒中的流体含水率与 油井正常生产时地面产液的含水率相同；油井正常生产时地面产液的含水率相同； 关井之后由抽油泵及油管漏失到井筒中的流体关井之后由抽油

26、泵及油管漏失到井筒中的流体的含水率与正常生产时地面产液的含水率相同。的含水率与正常生产时地面产液的含水率相同。 为了便于研究，特作如下假设：为了便于研究，特作如下假设：第三节第三节 根据静液面计算静压方法根据静液面计算静压方法 51因此，井筒中油水界面深度为：因此，井筒中油水界面深度为：wwowctosftctipiwimmbfVBVAAHHAAHHHHHHH21111第三节第三节 根据静液面计算静压方法根据静液面计算静压方法 522.2.计算气柱的压力降落计算气柱的压力降落03416. 0sgHTZcgePP3.3.计算油柱的压力降落计算油柱的压力降落 nisbiniioiHHooHHhhg

27、gdhPbs11661010第三节第三节 根据静液面计算静压方法根据静液面计算静压方法 534.4.计算水柱的压力降落计算水柱的压力降落 610bmwwHHgP第三节第三节 根据静液面计算静压方法根据静液面计算静压方法 54井筒多相流理论：井筒多相流理论：研究各种举升方式油井生产规律基本理论研究各种举升方式油井生产规律基本理论研究特点：研究特点：流动复杂性、无严格数学解流动复杂性、无严格数学解研究途径：研究途径：基本流动方程基本流动方程 实验资料相关因次分析实验资料相关因次分析 近似关系近似关系第四节第四节 井筒多相流理论与计算井筒多相流理论与计算55l 石油工业中的多相流技术研究始于石油工业

28、中的多相流技术研究始于19501950年左右年左右l 70 70年代，石油工业开始采用已在其它工业领域中使用的年代，石油工业开始采用已在其它工业领域中使用的 一些物理机理来预测多相流的流型一些物理机理来预测多相流的流型l 80 80年代初，计算机的引入极大地促进了多相流的发展年代初，计算机的引入极大地促进了多相流的发展l 80 80年代中期应用了核密度计、超声波、电容传感器、激光多年代中期应用了核密度计、超声波、电容传感器、激光多普勒仪和高速摄影机等设备对多相流机理进行研究普勒仪和高速摄影机等设备对多相流机理进行研究l 目前，双流体瞬态模拟方法和精确描述物理现象的稳态机理目前，双流体瞬态模拟方

29、法和精确描述物理现象的稳态机理模型是多相管流研究的主要方法模型是多相管流研究的主要方法1.1.多相流理论的研究简史多相流理论的研究简史第四节第四节 井筒多相流理论与计算井筒多相流理论与计算56 多相流问题未得到解析解多相流问题未得到解析解 复杂管道中的流动研究较少复杂管道中的流动研究较少 水平井段变质量流动研究较少水平井段变质量流动研究较少 缺乏专用研究仪器缺乏专用研究仪器 缺乏向下流动的综合机理模型缺乏向下流动的综合机理模型多相流理论的研究现存的问题多相流理论的研究现存的问题第四节第四节 井筒多相流理论与计算井筒多相流理论与计算57(1)(1)气液两相流动与单相液流的比较气液两相流动与单相液

30、流的比较第四节第四节 井筒多相流理论与计算井筒多相流理论与计算58流动型态（流动结构、流型）：流动过程中油、气的分布状态。流动型态（流动结构、流型）：流动过程中油、气的分布状态。(2)(2)气液混合物在垂直管中的流动结构变化气液混合物在垂直管中的流动结构变化影响流型的因素影响流型的因素各相介质的体积比例各相介质的体积比例介质的流速介质的流速各相的物理及化学性质各相的物理及化学性质( (密度、粘度界面张力等密度、粘度界面张力等) )流道的几何形状流道的几何形状壁面特性壁面特性管道的安装方式管道的安装方式第四节第四节 井筒多相流理论与计算井筒多相流理论与计算59流动型态的划分方法：两类流动型态的划

31、分方法：两类 第一类划分方法：根据两相介质分布的外形划分，包括泡状第一类划分方法：根据两相介质分布的外形划分，包括泡状流、弹状流或团状流、（层状流、波状流）、段塞流或冲击流、流、弹状流或团状流、（层状流、波状流）、段塞流或冲击流、环状流、雾状流环状流、雾状流 垂直气液两相流流型垂直气液两相流流型 水平气液两相流流型水平气液两相流流型第四节第四节 井筒多相流理论与计算井筒多相流理论与计算60第二类划分方法：按流动的数学模型或流体的分散程度划分，包括分第二类划分方法：按流动的数学模型或流体的分散程度划分，包括分散流、间歇流、分离流散流、间歇流、分离流 第一类划分方法第一类划分方法泡状流泡状流弹状流

32、或团状流弹状流或团状流层状流层状流波状流波状流段塞流或冲击流段塞流或冲击流环状流环状流雾状流雾状流第二类划分方法第二类划分方法分散流分散流间歇流间歇流分离流分离流分离流分离流间歇流间歇流分离流分离流分散流分散流两类划分结果两类划分结果的对应关系的对应关系第四节第四节 井筒多相流理论与计算井筒多相流理论与计算61两种分类方法比较两种分类方法比较 第一类划分方法较为直观第一类划分方法较为直观第二类划分方法便于进行数学处理第二类划分方法便于进行数学处理第四节第四节 井筒多相流理论与计算井筒多相流理论与计算62 油井生产中可能出现的流型自油井生产中可能出现的流型自下而上依次为：纯油下而上依次为：纯油(

33、 (液液) )流、泡流流、泡流、段塞流、环流和雾流。、段塞流、环流和雾流。 实际上，在同一口井内，一般实际上，在同一口井内，一般不会出现完整的流型变化。不会出现完整的流型变化。第四节第四节 井筒多相流理论与计算井筒多相流理论与计算油气沿井筒喷出时的流型变化示意图油气沿井筒喷出时的流型变化示意图纯油流；纯油流；泡流；泡流；段塞流；段塞流；环流；环流；雾流雾流63量从断面流出的流体能的能量在断面和之间耗失体额外所做的功在断面和之间对流进入断面的流体能量3.3.井筒气液两相流能量平衡方程井筒气液两相流能量平衡方程 两个流动断面间的能量平衡关系：两个流动断面间的能量平衡关系：222,VPU222,2m

34、ghmV121111,2,mghmVVPU倾斜管流能量平衡关系示意图第四节第四节 井筒多相流理论与计算井筒多相流理论与计算64222,VPU222,2mghmV121111,2,mghmVVPU倾斜管流能量平衡关系示意图倾斜管流能量平衡关系示意图2222221121112sin2sinVPmvmgZUqVPmvmgZU0)(sindqPVddZmgmvdvdU倾斜多相管流断面倾斜多相管流断面1 1和断面和断面2 2的的流体的能量平衡关系为：流体的能量平衡关系为：0sin1wdIdZgvdvdP第四节第四节 井筒多相流理论与计算井筒多相流理论与计算65适合于各种管流的通用压力梯度方程：适合于各种

35、管流的通用压力梯度方程：加速度摩擦举高)()()(dZdPdZdPdZdPdZdP则：则：令：令：2)()(sin)(2vdfdZdIdZdPdZdvvdZdPgdZdPw摩擦加速度举高2sin2mmmmmmmvdfdZdvvgdhdP第四节第四节 井筒多相流理论与计算井筒多相流理论与计算66说明：说明：a. a. 计算压力分布过程中，温度和压力是相关的；计算压力分布过程中，温度和压力是相关的；b. b. 流体物性参数计算至关重要，但目前方法精度差；流体物性参数计算至关重要，但目前方法精度差；c. c. 不同的多相流计算方法差别较大，因此在实际应用中不同的多相流计算方法差别较大，因此在实际应用

36、中有必要根据油井的实际情况筛选精度相对高的方法。有必要根据油井的实际情况筛选精度相对高的方法。第四节第四节 井筒多相流理论与计算井筒多相流理论与计算674.4.井筒气液两相流动的模型井筒气液两相流动的模型 气液两相流动的规律较单相流复杂得多气液两相流动的规律较单相流复杂得多, ,常采用简化的流动模型常采用简化的流动模型进行处理，以便探讨其流动规律。常用的模型有进行处理，以便探讨其流动规律。常用的模型有流动型态模型、均流动型态模型、均相流动模型、分相流动模型和漂移流动模型相流动模型、分相流动模型和漂移流动模型等。等。第四节第四节 井筒多相流理论与计算井筒多相流理论与计算68（1 1）流动型态模型

37、）流动型态模型定义：将两相流动分成几个典型的流动型态，然后按照不同的定义：将两相流动分成几个典型的流动型态，然后按照不同的流动型态分别研究其流动规律。流动型态分别研究其流动规律。特特点点针对性强，精确度高针对性强，精确度高数学处理复杂，计算量大数学处理复杂，计算量大流型界限确定困难流型界限确定困难流动型态模型是未来的研究方向流动型态模型是未来的研究方向第四节第四节 井筒多相流理论与计算井筒多相流理论与计算69（1 1）气相和液相的实际速度相等，即）气相和液相的实际速度相等，即vvvgl（2 2）两相介质已达到热力学平衡状态）两相介质已达到热力学平衡状态对于泡状流和雾状流，具有较高的精确性对于泡

38、状流和雾状流，具有较高的精确性对于弹状流和段塞流，需要进行时间平均修正对于弹状流和段塞流，需要进行时间平均修正对于层状流、波状流和环状流，则误差较大对于层状流、波状流和环状流，则误差较大两个假定两个假定特点特点定义：把气液两相混合物看成均匀介质，其物定义：把气液两相混合物看成均匀介质，其物 性参数取两相的均值而建立的模型性参数取两相的均值而建立的模型第四节第四节 井筒多相流理论与计算井筒多相流理论与计算70 定义：它是把气液两相流动看成为气、液相各自分开的流动，定义：它是把气液两相流动看成为气、液相各自分开的流动，每相介质都有其平均流速和独立的物性参数。因此需要建立每一相每相介质都有其平均流速

39、和独立的物性参数。因此需要建立每一相介质的流体动力特性方程式。这就要求预先确定每一相占有过流断介质的流体动力特性方程式。这就要求预先确定每一相占有过流断面的份额（即真实含气率）以及介质与管壁的摩擦力和两相介质之面的份额（即真实含气率）以及介质与管壁的摩擦力和两相介质之间的摩擦阻力，为了取得这些数据，目前主要是利用实验研究所得间的摩擦阻力，为了取得这些数据，目前主要是利用实验研究所得的经验关系式。的经验关系式。 近年来，随着计算流体力学的发展，有些数据已可以通过数学近年来，随着计算流体力学的发展，有些数据已可以通过数学模型用解析计算求得。模型用解析计算求得。 第四节第四节 井筒多相流理论与计算井

40、筒多相流理论与计算71分流模型的基本假设是：分流模型的基本假设是： (1) (1) 两相介质分别有各自的按所占断面积计算的断面平均流速；两相介质分别有各自的按所占断面积计算的断面平均流速； (2) (2) 虽然两相介质之间可能有质量交换，但两相之间处于热力学虽然两相介质之间可能有质量交换，但两相之间处于热力学平衡状态，压力和密度互为单值函数。平衡状态，压力和密度互为单值函数。分流模型适用于层状流、波状流和环状流。分流模型适用于层状流、波状流和环状流。第四节第四节 井筒多相流理论与计算井筒多相流理论与计算72 漂移流动模型简称漂移模型，它是由朱伯漂移流动模型简称漂移模型，它是由朱伯(Zuber)

41、(Zuber)和芬德和芬德莱莱(Findlay)(Findlay)针对均流模型、分流模型与实际的两相流动之间针对均流模型、分流模型与实际的两相流动之间存在的偏差而提出的特殊模型。存在的偏差而提出的特殊模型。 在均流模型中，没有考虑两相间的相互作用，而是用平均在均流模型中，没有考虑两相间的相互作用，而是用平均流动参数来模拟两相介质；分流模型中，尽管在流动特性方流动参数来模拟两相介质；分流模型中，尽管在流动特性方面分别考虑了每相介质以及两相界面上的作用力，但是每相面分别考虑了每相介质以及两相界面上的作用力，但是每相的流动特性仍然是孤立的；而在漂移流动模型中，既考虑了的流动特性仍然是孤立的；而在漂移

42、流动模型中，既考虑了气液两相之间的相对运动，又考虑了空隙率和流速沿过流断气液两相之间的相对运动，又考虑了空隙率和流速沿过流断面的分布规律。面的分布规律。第四节第四节 井筒多相流理论与计算井筒多相流理论与计算73 在实际研究过程中，不仅要根据所研究的问题选择合适的在实际研究过程中，不仅要根据所研究的问题选择合适的多相流计算模型，而且常常需将前面介绍的多种模型有机地结多相流计算模型，而且常常需将前面介绍的多种模型有机地结合起来使用，以各取其长，获得尽可能精确的计算结果。合起来使用，以各取其长，获得尽可能精确的计算结果。5.5.井筒气液两相流动的计算方法井筒气液两相流动的计算方法第四节第四节 井筒多

43、相流理论与计算井筒多相流理论与计算74 Poettmann-Carpenter Poettmann-Carpenter方法、方法、Fanch-BrownFanch-Brown相关式、相关式、Baxendell-Thomas Baxendell-Thomas 相关式、相关式、Hagedron-Brown Hagedron-Brown 关系式、关系式、Duns-RosDuns-Ros相关式、相关式、OrkiszewskiOrkiszewski相关式、相关式、Beggs-Brill Beggs-Brill 相关式、相关式、DuklerDukler相关式、相关式、Mukherjee-BrillMukh

44、erjee-Brill相关式、相关式、 Aziz Aziz相关式、相关式、EatonEaton相相关式、关式、AnsariAnsari相关式等。相关式等。 不同方法有其适用条件和精度，可根据具体油田实际选用。不同方法有其适用条件和精度，可根据具体油田实际选用。第四节第四节 井筒多相流理论与计算井筒多相流理论与计算75数学相关式大体分为三种类型：数学相关式大体分为三种类型：在计算井筒流体混合密度时不考虑液体滞留量的影响，而在计算井筒流体混合密度时不考虑液体滞留量的影响，而液体滞留量与管壁摩阻损失用一个经验摩阻系数来表达，不液体滞留量与管壁摩阻损失用一个经验摩阻系数来表达，不区分多相流体的流态分布

45、情况。区分多相流体的流态分布情况。在计算多相流体混合物密度时考虑液体滞留量的影响，而在计算多相流体混合物密度时考虑液体滞留量的影响，而摩阻系数要依据液体和气体的组成特征来确定。摩阻系数要依据液体和气体的组成特征来确定。考虑液体滞留量的影响，摩阻系数取决于多相流体中连续考虑液体滞留量的影响，摩阻系数取决于多相流体中连续相的特征。同时考虑了不同流态的影响，其中包括泡状流、相的特征。同时考虑了不同流态的影响，其中包括泡状流、段塞流、环流及雾状流。段塞流、环流及雾状流。第四节第四节 井筒多相流理论与计算井筒多相流理论与计算76 在进行气液两相流动的有关计算中，常需要确定工作条件下在进行气液两相流动的有

46、关计算中，常需要确定工作条件下原油、天然气、水及其混合物的物性参数。客观地讲，确定这些原油、天然气、水及其混合物的物性参数。客观地讲，确定这些物性参数最根本、最精确的方法是实验测定。然而，实际生产设物性参数最根本、最精确的方法是实验测定。然而，实际生产设计和计算中所遇到的原油、天然气及水的组成、工作温度和工作计和计算中所遇到的原油、天然气及水的组成、工作温度和工作压力等的范围都非常广泛，完全依赖实验方法测定各种工况条件压力等的范围都非常广泛，完全依赖实验方法测定各种工况条件下的油、气、水及其混合物的物性参数是很困难的。另外，过去下的油、气、水及其混合物的物性参数是很困难的。另外，过去曾一度建立

47、和使用的许多物性参数图版，也都难以适应目前广泛曾一度建立和使用的许多物性参数图版，也都难以适应目前广泛应用电子计算机进行工程计算的要求。应用电子计算机进行工程计算的要求。 因此，为了便于利用电子计算机进行气液两相流动的计算，因此，为了便于利用电子计算机进行气液两相流动的计算，建立原油、天然气、水及其混合物物性参数计算的相关公式是非建立原油、天然气、水及其混合物物性参数计算的相关公式是非常必要的。常必要的。第四节第四节 井筒多相流理论与计算井筒多相流理论与计算77 1. 1.假设条件假设条件l气体质量忽略不计气体质量忽略不计l井筒中液体流动为准稳定流，体积和流井筒中液体流动为准稳定流，体积和流态

48、变化的影响忽略不计态变化的影响忽略不计l流体对地层放热，其总传热系数为不变流体对地层放热，其总传热系数为不变l因天然气析出及膨胀吸热忽略不计因天然气析出及膨胀吸热忽略不计l油流在油管中流动时因摩擦而产生的热油流在油管中流动时因摩擦而产生的热量忽略不计量忽略不计第五节第五节 井筒多相流体温度分布计算井筒多相流体温度分布计算78 2. 2.公式的推导公式的推导传热公式传热公式DdhttKdqs油损失热量为油损失热量为GCdtdq根据热量守恒根据热量守恒GCdtDdhttKs地层温度与深度的关系地层温度与深度的关系httss1联立上述公式可得联立上述公式可得GCdtdhhttDKs1求解方程得求解方

49、程得GCDhKseDKGChtt11C C为原油为原油比热比热第五节第五节 井筒多相流体温度分布计算井筒多相流体温度分布计算793.3.公式中各项参数取值公式中各项参数取值参数取值参数取值(1) (1) 井底油温井底油温 井底油温也就是油层温度井底油温也就是油层温度(2) (2) 地温梯度地温梯度一般认为一般认为 0.030.03m m(3) (3) 井深井深h h取油层中部至井口的距离取油层中部至井口的距离(4(4）原油质量流量）原油质量流量油井生产时可实测，生产前由地油井生产时可实测，生产前由地质部门或作业者提供质部门或作业者提供(5(5）原油比热）原油比热C C一般计算时取一般计算时取C

50、 C21kJ21kJkg.kg.当原当原油含水时油含水时水水油油CGCGGC第五节第五节 井筒多相流体温度分布计算井筒多相流体温度分布计算80(6(6）总传热系数）总传热系数K K021111RKii(7) (7) 套管外径套管外径D D其数据在完井后即可提供其数据在完井后即可提供第五节第五节 井筒多相流体温度分布计算井筒多相流体温度分布计算814.4.假设条件的分析假设条件的分析(1)(1)气体质量流量气体质量流量(2)(2)流动型态流动型态 国内多数油田的气油比一般在国内多数油田的气油比一般在10108080之间，气量可忽略之间，气量可忽略不计，假设条件（不计，假设条件（1 1）可以成立。

51、但对气油比大于）可以成立。但对气油比大于100100的井，的井，可将气量折换成油量。可将气量折换成油量。 井筒内的流型较复杂，但考虑到只研究宏观整体，故可井筒内的流型较复杂，但考虑到只研究宏观整体，故可将流动视为准稳定流。将流动视为准稳定流。第五节第五节 井筒多相流体温度分布计算井筒多相流体温度分布计算82(3)(3)总传热系数总传热系数(4)(4)天然气析出和膨胀问题天然气析出和膨胀问题 井筒外部不同位置的岩层性质亦不同，井筒内的流体流型又井筒外部不同位置的岩层性质亦不同，井筒内的流体流型又有变化，所以严格地说，总传热系数值应该是一个变量，但其变有变化，所以严格地说，总传热系数值应该是一个变

52、量，但其变化与平均值之间的差值不大，取实测平均值即可。因此，计算时化与平均值之间的差值不大，取实测平均值即可。因此，计算时可按常数考虑。可按常数考虑。 油流沿井筒上升，当压力低于饱和压力时，天然气析出。析油流沿井筒上升，当压力低于饱和压力时，天然气析出。析出气体需要热量，已析出气体不断膨胀，又会吸收一部分热量，出气体需要热量，已析出气体不断膨胀，又会吸收一部分热量，这两部分热量的计算比较复杂。但在井底压力大于饱和压力且气这两部分热量的计算比较复杂。但在井底压力大于饱和压力且气油比小于油比小于100100时，气体析出膨胀而引起温降一般在时，气体析出膨胀而引起温降一般在22左右。而且左右。而且此温

53、降又被油气在油管中由于摩擦所产生的热量补偿掉一部分。此温降又被油气在油管中由于摩擦所产生的热量补偿掉一部分。因此，在一般工程计算中，可忽略此温降。因此，在一般工程计算中，可忽略此温降。第五节第五节 井筒多相流体温度分布计算井筒多相流体温度分布计算83第六节第六节 流体物性参数计算方法流体物性参数计算方法 1.1.原油密度原油密度 BRgso310206. 110002.2.原油的原油的APIAPI度度 5 .1315 .141oAPI3.3.油水混合液体的密度油水混合液体的密度 Lwwwff1844.4.原油体积系数原油体积系数 API当当 3030时时： gAPIabrrtRsRsB6032

54、8 . 1)108106. 110751. 1 (104.677+1=85-4gAPIabrrtRsRsB60328 . 1)10337. 1101 . 1 (104.677+1=95-4API当当 3030时时： 第六节第六节 流体物性参数计算方法流体物性参数计算方法 85当平均压力小于饱和压力时：当平均压力小于饱和压力时： bBB当平均压力大于饱和压力时：当平均压力大于饱和压力时： 51022.1461.121180328 . 12 .1751433ParrtRpCAPIga其中：其中：)22.13(abppCbeBB第六节第六节 流体物性参数计算方法流体物性参数计算方法 865.5.溶解

55、气油比溶解气油比 API当当 3030时时： 67.459328 . 1172.110937. 1s1064.2722.14= RaAPItragp67.459328 .1393.10187.11006.5622.14aAPItragspRAPI当当 3030时时： 第六节第六节 流体物性参数计算方法流体物性参数计算方法 876.6.液体粘度液体粘度 )ln()1ln(lnaoDtba 死油粘度：死油粘度：活油粘度：活油粘度：1000100BODouA水的粘度：水的粘度：10002528 . 13210982. 18 . 13210479. 1003. 1ttwe油水混合物粘度：油水混合物粘度

56、：wwwolfufuu1第六节第六节 流体物性参数计算方法流体物性参数计算方法 887.7.油、天然气的表面张力油、天然气的表面张力 ogAPIPte4240047 1832026710001015 107.8.8.水、天然气的表面张力水、天然气的表面张力 78.13778.13733.232068 . 1248twg9 9、油、水混合物和天然气的表面张力、油、水混合物和天然气的表面张力wwgwoglff1第六节第六节 流体物性参数计算方法流体物性参数计算方法 8910.10.天然气的粘度天然气的粘度 uKegxgy104103tMtMK8 . 1197018 . 149202. 04 . 9

57、5 . 1y=2.4-0.2xx=3.5+986/(492+1.8t)+0.01MM=28.96其中：其中：第六节第六节 流体物性参数计算方法流体物性参数计算方法 9011.11.天然气的密度天然气的密度 ggPZ t3484410273153.12.12.天然气的压缩因子天然气的压缩因子ZAATATAATATrrRrRRr11233452623/第六节第六节 流体物性参数计算方法流体物性参数计算方法 91第七节第七节 抽油机井生产节点系统分析抽油机井生产节点系统分析1.1.节点系统分析方法节点系统分析方法 节点系统分析法：节点系统分析法：应用系统工程原理，把整个油应用系统工程原理，把整个油井

58、生产系统分成若干子系统，研究各子系统间的相互井生产系统分成若干子系统，研究各子系统间的相互关系及其对整个系统工作的影响，为系统优化运行及关系及其对整个系统工作的影响，为系统优化运行及参数调控提供依据。参数调控提供依据。节点划分依据：节点划分依据：不同的流动规律相关式不同的流动规律相关式922.2.抽油机井生产系统的组成抽油机井生产系统的组成(1)(1)油气层子系统油气层子系统 (2)(2)井筒子系统井筒子系统 (3)(3)地面集输子系统地面集输子系统 (4)(4)采油设备子系统采油设备子系统 稳定工作条件：协调稳定工作条件：协调 第七节第七节 抽油机井生产节点系统分析抽油机井生产节点系统分析9

59、3节点系统分析对象：整个油井生产系统节点系统分析对象：整个油井生产系统节点系统分析实质：协调理论在采油应用方面的发展节点系统分析实质：协调理论在采油应用方面的发展协协调调条条件件质量守恒质量守恒能量能量( (压力压力) )守恒守恒热量守恒热量守恒第七节第七节 抽油机井生产节点系统分析抽油机井生产节点系统分析94求解点的选择：求解点的选择：主要取决于所要研究解决的问题。主要取决于所要研究解决的问题。求解点：求解点：为使问题获得解决的节点。为使问题获得解决的节点。0510152025010203040506070产 量压力节点流入曲线节点流入曲线节点流出曲线节点流出曲线协调点协调点第七节第七节 抽

60、油机井生产节点系统分析抽油机井生产节点系统分析950510152025010203040506070产 量压力节点流入曲线节点流入曲线节点流出曲线节点流出曲线协调点协调点第七节第七节 抽油机井生产节点系统分析抽油机井生产节点系统分析96节点系统分析思路：节点系统分析思路：以系统两端为起点分别计算不同流量下节点上、下游的压以系统两端为起点分别计算不同流量下节点上、下游的压力，并求得力，并求得节点压差流量曲线。节点压差流量曲线。根据描述节点设备根据描述节点设备( (油嘴、安全阀等油嘴、安全阀等) )的流量的流量压差相关式，压差相关式，求得求得设备工作曲线设备工作曲线。两条压差流量两条压差流量曲线的