采油工程课程设计

采油工程课程设计课程设计姓名：孔令伟学号： 201301509287中国石油大学（北京）石油工程学院2014年 10月 30日一、给定设计基础数据： (2)二、设计计算步骤 (3)2.1油井流入动态计算 (3)2.2 井筒多相流的计算 (4)2.3 悬点载荷和抽油杆柱设计计算 (12)2.4 抽油机校核 (16)2.5 泵效计算 (16)2.6 举升效率计算 (19)三、设计计算总结果 (22)四、课程设计总结 (23)一、给定设计基础数据：井深： 2000+87×10=2870m套管内径： 0.124m油层静压： 2870/100×1.2 =34.44MPa油层温度： 90℃恒温层温度： 16℃地面脱气油粘度： 30mPa.s油相对密度： 0.84气相对密度： 0.76水相对密度： 1.0油饱和压力： 10MPa含水率： 0.4套压： 0.5MPa油压： 1 MPa生产气油比： 50m3/m3原产液量 (测试点 )：30t/d原井底流压 (测试点 )：16.35Mpa抽油机型号： CYJ10353HB电机额定功率： 37kw配产量： 50t/d泵径： 56mm冲程： 3m冲次： 6rpm柱塞与衬套径向间隙 :0.3mm沉没压力： 3MPa二、设计计算步骤2.1 油井流入动态计算油井流入动态是指油井产量与井底流动压力的关系， 它反映了油藏向该井供油的能力。

从单井来讲， IPR 曲线表示了油层工作特性。

因而，它既是确定油井合理工作方式的依据， 也是分析油井动态的基础。

本次设计油井流入动态计算采用 Petro bras 方法 Petro bras 方法计算综合 IPR 曲线的实质是按含水率取纯油 IPR 曲线和水 IPR 曲线的加权平均值。

当已知测试点计算采液指数时， 是按产量加权平均 ;预测产量时 ,按流压加权平均。

(1) 采液指数计算已知一个测试点：P wftest、qtxest和饱和压力 Pb及油藏压力 P 。

qtxwst因为 PwftestP=PPwfest=30/(34.44-12)= 1.3/( d.Mpa)b， j1(2) 某一产量qt下的流压 Pwfq b =j(P 1P b)=1.4 x （34.44-10）=34.22t/djP bqo m zx=q b+ 1.8=34.44+1.4*10/1.8=42.22t/dq omzx-油 IPR 曲线的最大产油量。

q t当 0q b 时，令 q，则P 1t t1wf 1jq=10 t/dp= =15.754 Mpa同理， q t 2 =20 t/d ， P wf 2 =13.877 Mpaq t 3 =30 t/d ， P wf 3 =12.0 Mpa当 q b q tq omzx时，令 q t 4 =50 t/d,则按流压加权平均进行推导得：w( P 1 q t )81 80(q tq b )]P wf 4 =f j +0.125(1-f w )P b [-1+qomzxq b=8.166Mpa同理 q t 5=60t/d， P wf 5 =5.860 Mpa当q omzxt p wf f w ( p1q omzx)( q t q omzx)(8 f w 9) q 时，J J令q t 6 =71t/d，P wf 6 =2.233 Mpa综上，井底流压与产量的关系列表如下：Pwf/Mpa15.74713.87312.010.08.166 5.860 2.233Q/(t/d)10203040.653506071得到油井的流入动态曲线如下图：IPR 曲线20)15apM(10fwP50102030405060708090Q（t/d ）图1 油井 IPR 曲线2.2井筒多相流的计算井筒多相流压力梯度方程井筒多相管流的压力梯度包括：因举高液体而克服重力所需的压力势能、流体因加速而增加的动能和流体沿管路的摩阻损失，其数学表达式如下：dp dv m v m2dhρ m gsinθ+ρm v m dhf mρm /d*2式中ρ m为多相混合物的密度;v m为多相混合物的流速 ;f m为多相混合物流动时的摩擦阻力系数 ;d 为管径 ;p 为压力 ;h 为深度 ;g 为重力加速度 ; θ为井斜角的余角。

井筒多相管流计算包括两部分：（1）由井底向上计算至泵入口处；（2）油管内由井口向下计算至泵出口处。

1）由井底向上计算至泵入口处，计算下泵深度 Lp 。

采用深度增量迭代方法，首先估算迭代深度。

在本设计中为了减小工作量，采用只迭代一次的方法。

计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用Orkiszewski 方法判断流型，进行压力梯度计算，最后计算出深度增量和下泵深度Lp。

按深度增量迭代的步骤 :井底流压 12Mpa，假设压力降为 0.2 Mpa；估计一个对应的深度增量h，=40m 即深度为 1960m 。

由井温关系式可以计算得到该处的井温为：89.96℃。

平均的压力和温度：T=（90+89.96）/2=89.98℃。

平均压力P=11.9 Mpa。

由平均压力和平均温度计算的得到流体的物性参数为：溶解油气比R S =71.31 ；原油体积系数 B0 =1.25原油密度P 0 =739.00；油水混合液的密度P z4O4=843.40；死油粘度μod =6.537*10 ；活油粘度μ =3.318*10；水的粘44度μw =3.263*10；液体的粘度μ= 3.296*10 ；天然气的压缩因子 Z=0.9567；天然气的密度g90.70。

以上单位均是标准单位。

由以上的流体物性参数判断流型：不同流动型态下的m和 f 的计算方法不同，为此，计算中首先要判断流动形态。

该方法的四种流动型态的划分界限如表 1 所示。

表 1 流型界限流动型态界限泡流q g L Bq t段塞流q gBg S q t L ,v L过渡流雾流L M v g L S v g L M其中LB =1.071-0.7277nt2/ D且LB >0.13（如果LB <0.13，则取LB =0.13）；q gL S=50+36v g q；tq g0.75L M v g q t=75+84（）。

由计算得到，由于该段的压力大于饱和压力的值，所以该段的流型为纯液流。

dP计算该段的压力梯度dh。

由压力梯度的计算公式：dP f dh g m dh m v m dv mm=843.40；f =计算对应于P的该段管长 (深度差 ) h计。

⑥将第步计算得的h计与第②步估计的h进行比较，两者之差超过允许范围，则以新的h作为估算值，重复②～⑤的计算，使计算的与估计的h 之差在允许范围内为止。

该过程之中只迭代一次。

2）由井口向下计算至泵出口处，计算泵排出口压力PZ。

采用压力增量迭代方法，首先估算迭代压力。

同样为了减小工作量，也采用只迭代一次的方法。

计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用 Orkiszewski 方法判断流型，进行压力梯度计算，最后计算出压力增量和泵排出口压力PZ。

按压力增量迭代的步骤①已知任一点 (井底或井口 )的压力P0，选取合适的深度间隔h(可将管L等分为 n 段 )。

②估计一个对应于计算间隔h的压力增量P 。

T PdP④计算该段压力梯度dhoP i h dP⑤计算对应于h的压力增量dho⑥比较压力增量的估计量P 与计算值P i，若二者之差不在允许范围内，则以计算值作为新的估计值，重复第②～⑤步，使两者之差在允许范围o 之内为止。

⑦计算该段下端对应的深度L和压力Pi iiL i i h P i P o P i1⑧以Li 处的压力Pi 为起点压力重复第②～⑦步，计算下一段的深度L i 1和压力Pi 1，直到各段累加深度等于或大于管长L 时为止。

2.2 计算气 -液两相垂直管流的 Orkiszewski 方法本设计井筒多相流计算采用Orkiszewski 方法。

Orkiszewski 法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。

如图1 所示。

在处理过渡性流型时，采用内插法。

在计算段塞流压力梯度时要考虑气相与液体的分布关系。

针对每种流动型态提出了存容比及摩擦损失的计算方法。

图1 气液混合物流动型态 (Orkiszewski)1.压力降公式及流动型态划分界限由前面垂直管流能量方程可知，其压力降是摩擦能量损失、势能变化和动能变化之和。

由式 (2-36)可直接写出多项垂直管流的压力降公式：dP f dhg m dh m v m dv m(26)式中P —压力，Pa；f—摩擦损失梯度， Pa/m；h—深度， m；g—重力加速度， m/s2；m—混合物密度， kg/m3；vm—混合物流速，m/s。

动能项只是在雾流情况下才有明显的意义。

出现雾流时，气体体积流量远大于液体体积流量。

根据气体定律，动能变化可表示为：m v m dv m W t q g2dpA P(27)p式中Ap—管子流通截面积， m2；tW—流体总质量流量， kg/s；qg —气体体积流量，m3/s。

将式 (27)代入式 (26)，并取dhh k， dP p k，m m， PP经过整理后可得：[m g f ]h k1W t q gP k A p2 P(28)式中Pk —计算管段压力降，Pa；hk —计算管段的深度差，m；P—计算管段的平均压力， Pa。

不同流动型态下的m 和 f 的计算方法不同,下面按流型分别介绍。

(1)泡流平均密度mH L L H g g 1 H g L H g g式中H g— 气相存容比 (含气率 )，计算管段中气相体积与管段容积之比值；HL—液相存容比 (持液率 )，计算管段中液相体积与管段容积之比值；g、Ｌ、m—在P、T下气、液和混合物的密度，kg/m3。

气相存容比由滑脱速度Vs来计算。

滑脱速度定义为： 气相流速与液相流速之差。

v sv sgv q gq t q gsLH g 1 H gA p H gA p (1 H g )1[1q t (1 q t)24q g可解出 H g： H g 2 v s A pv s A p v s A p式中v s—滑脱速度，由实验确定， m/s ；v sg、vsL—气相和液相的表观流速， m/s 。

泡流摩擦损失梯度按液相进行计算：Lv LH 2v L Hq L tf D 2A p (1 H g )式中f —摩擦阻力系数；vLH—液相真实流速， m/s 。

摩擦阻力系数f可根据管壁相对粗造度/ D 和液相雷诺数 NRe 查图 2。

N ReDv sL L液相雷诺数 :L式中L—在P、T下的液体粘度，油、水混合物在未乳化的情况下可取其体积加权平均值， Pa.s 。

图2(2)段塞流混合物平均密度W t L v s A pm L(34)q t v s A p式中—液体分布系数；vs —滑脱速度，m/s。