采油工程工课程设计

1、采油工程课程设计采油工程课程设计课程设计课程设计姓名：孔令伟姓名：孔令伟 学号：学号：9287中国石油大学（北京）中国石油大学（北京）石油工程学院石油工程学院20142014 年年 1010 月月 3030 日日1一、给定设计基础数据：一、给定设计基础数据：.1 1二、设计计算步骤二、设计计算步骤.2 2 油井流入动态计算 .2 井筒多相流的计算 .3悬点载荷和抽油杆柱设计计算.11 抽油机校核 .15 泵效计算 .15举升效率计算.18三、设计计算总结果三、设计计算总结果.1919四、课程设计总结四、课程设计总结.20202一、给定设计基础数据：一、给定设计基础数据： 井深：2000+871

2、0=2870m 套管内径： 油层静压：2870/100 = 油层温度：90 恒温层温度：16 地面脱气油粘度： 油相对密度： 气相对密度： 水相对密度： 油饱和压力：10MPa 含水率： 套压： 油压：1 MPa 生产气油比：50m3/m3 原产液量(测试点)：30t/d 原井底流压(测试点)： 抽油机型号：CYJ10353HB 电机额定功率：37kw 配产量：50t/d 泵径：56mm 冲程：3m 冲次：6rpm 柱塞与衬套径向间隙:0.3mm 沉没压力：3MPa 3二、设计计算步骤二、设计计算步骤 油井流入动态计算油井流入动态计算油井流入动态是指油井产量与井底流动压力的关系，它反映了油藏向

3、该井供油的能力。从单井来讲，IPR 曲线表示了油层工作特性。因而，它既是确定油井合理工作方式的依据，也是分析油井动态的基础。本次设计油井流入动态计算采用 Petro bras 方法 Petro bras 方法计算综合 IPR 曲线的实质是按含水率取纯油 IPR 曲线和水 IPR 曲线的加权平均值。当已知测试点计算采液指数时，是按产量加权平均;预测产量时,按流压加权平均。(1) 采液指数计算已知一个测试点：wftestP、txestq和饱和压力bP及油藏压力P。因为wftestPbP，1j=txwstwfestqPP=30/= ( (2) 某一产量tq下的流压 Pwf bq=j(bPP 1)=

4、x（）=dmozxq=bq+8 . 1bjP=+*10/=domzxq-油 IPR 曲线的最大产油量。当 0qtbq时，令 q1t=10 t/d，则 p1wf=jqPt1= Mpa同理，q2t=20 t/d，P2wf= Mpa q3t=30 t/d，P3wf= Mpa当 qbqtomzxq时，令 q4t=50 t/d,则按流压加权平均进行推导得：P4wf=f)(1jqPtw+(1-fw)Pb-1+8180()tbomzxbqqqq=4 同理 q5t=60t/d，P5wf= Mpa当 qomzxqt时，1()(89)()omzxtomzxwwfwqqqfpfpJJ令 q6t=71t/d，P6wf

5、= Mpa综上，井底流压与产量的关系列表如下：Pwf/MpaQ/(t/d)102030506071 得到油井的流入动态曲线如下图： IPR 曲线051015200102030405060708090Q（t/d）Pwf(Mpa)图 1 油井 IPR 曲线 井筒多相流的计算井筒多相流的计算井筒多相流压力梯度方程井筒多相管流的压力梯度包括：因举高液体而克服重力所需的压力势能、流体因加速而增加的动能和流体沿管路的摩阻损失，其数学表达式如下：dhdpmgsin+mvmmmfdhdvm/d*22mv式中 m为多相混合物的密度;vm为多相混合物的流速;fm为多相混合物流5动时的摩擦阻力系数;d 为管径;p

6、为压力;h 为深度;g 为重力加速度; 为井斜角的余角。井筒多相管流计算包括两部分：（1）由井底向上计算至泵入口处；（2）油管内由井口向下计算至泵出口处。1）由井底向上计算至泵入口处，计算下泵深度 Lp。采用深度增量迭代方法，首先估算迭代深度。在本设计中为了减小工作量，采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用Orkiszewski 方法判断流型，进行压力梯度计算，最后计算出深度增量和下泵深度 Lp。按深度增量迭代的步骤:井底流压 12Mpa，假设压力降为 Mpa；估计一个对应的深度增量h=40m，即深度为 1960m 。由井温关系式可以计算得到该处的

7、井温为：。平均的压力和温度：T=（90+）/2=。平均压力P= Mpa。由平均压力和平均温度计算的得到流体的物性参数为：溶解油气比 RS= ； 原油体积系数 B0= 原油密度 P0=； 油水混合液的密度 Pz =； 死油粘度 od=*104； 活油粘度 O=*104； 水的粘度 w=*104； 液体的粘度 = *104；天然气的压缩因子 Z=； 天然气的密度g。以上单位均是标准单位。由以上的流体物性参数判断流型：不同流动型态下的m和f的计算方法不同，为此，计算中首先要判断流动形态。该方法的四种流动型态的划分界限如表 1 所示。表 1 流型界限流动型态界 限泡 流BtgLqq段 塞 流SgBtg

8、LvLqq,过 渡 流SgMLvL6雾 流MgLv 其中BL=2/tD且BL（如果BL，则取BL=） ； SL=50+36 gv gtqq； ML=75+84 （gvgtqq）0.75。由计算得到，由于该段的压力大于饱和压力的值，所以该段的流型为纯液流。计算该段的压力梯度dhdP。由压力梯度的计算公式：mmmmfdvvdhgdhdP m=；f=计算对应于P的该段管长(深度差)、h。 将第步计算得的、h与第步估计的h进行比较，两者之差超过允许范围，则以新的h作为估算值，重复的计算，使计算的与估计的h之差在允许范围内为止。该过程之中只迭代一次。2）由井口向下计算至泵出口处，计算泵排出口压力 PZ。

9、采用压力增量迭代方法，首先估算迭代压力。同样为了减小工作量，也采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用Orkiszewski 方法判断流型，进行压力梯度计算，最后计算出压力增量和泵排出口压力 PZ。按压力增量迭代的步骤已知任一点(井底或井口)的压力0P， 选取合适的深度间隔h(可将管L等分为 n 段)。估计一个对应于计算间隔h的压力增量P。计算该段的T和P ，以及P、T下的流体性质参数。7计算该段压力梯度odhdP计算对应于h的压力增量oidhdPhP比较压力增量的估计量P与计算值iP ，若二者之差不在允许范围内，则以计算值作为新的估计值，重复第步，

10、使两者之差在允许范围o之内为止。计算该段下端对应的深度iL和压力iP hiLi iioiPPP1以iL处的压力iP为起点压力重复第步，计算下一段的深度1iL和压力1iP ，直到各段累加深度等于或大于管长L时为止。计算气-液两相垂直管流的 Orkiszewski 方法本设计井筒多相流计算采用 Orkiszewski 方法。Orkiszewski 法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。如图 1 所示。在处理过渡性流型时，采用内插法。在计算段塞流压力梯度时要考虑气相与液体的分布关系。针对每种流动型态提出了存容比及摩擦损失的计算方法。图 1 气液混合物流动型态(Orkiszewski)

11、1.压力降公式及流动型态划分界限由前面垂直管流能量方程可知，其压力降是摩擦能量损失、势能变化和动能变化之和。由式(2-36)可直接写出多项垂直管流的压力降公式：mmmmfdvvdhgdhdP (26)式中 P压力，Pa；8 f摩擦损失梯度，Pa/m； h深度，m； g重力加速度，m/s2； m混合物密度，kg/m3； mv混合物流速，m/s。动能项只是在雾流情况下才有明显的意义。出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：dpPAqWdvvpgtmmm2 (27)式中 pA管子流通截面积，m2； tW流体总质量流量，kg/s； gq气体体积流量，m3/s。将式(

12、27)代入式(26)，并取khdh，kpdP ，mm，PP经过整理后可得：kPkpgtfmhPAqWg12 (28)式中 kP计算管段压力降，Pa； kh计算管段的深度差，m； P计算管段的平均压力，Pa。不同流动型态下的m和f的计算方法不同,下面按流型分别介绍。(1)泡流平均密度 ggLLmHHggLgHH 1 1gLHH 式中 gH气相存容比(含气率)，计算管段中气相体积与管段容积之比值； LH液相存容比(持液率)，计算管段中液相体积与管段容积之比值；9 mg、在TP、下气、液和混合物的密度，kg/m3。气相存容比由滑脱速度sV来计算。滑脱速度定义为：气相流速与液相流速之差。)1 (1gp

13、gtgpggsLgsgsHAqqHAqHvHvv 可解出gH：HgpsgpstpstAvqAvqAvq4)1 (1 212 式中 sv滑脱速度，由实验确定，m/s； sgv、sLv气相和液相的表观流速，m/s。泡流摩擦损失梯度按液相进行计算： 22LHLtvDf )1 (gpLLHHAqv 式中 f摩擦阻力系数； LHv液相真实流速，m/s。摩擦阻力系数f可根据管壁相对粗造度D/和液相雷诺数ReN查图 2。液相雷诺数: LLsLDvNRe 式中 L在TP、下的液体粘度，油、水混合物在未乳化的情况下可取其体积加权平均值， 。10图 2(2)段塞流混合物平均密度LpstpsLtmAvqAvW (3

14、4)式中 液体分布系数； sv滑脱速度，m/s。滑脱速度可用 Griffith 和 Wallis 提出的公式计算：gDCCvs21 (35) (3)过渡流过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，11然后用内插方法来确定相应的数值。MiSMsgSLSMgMmLLLvLLvL (36)MiSMggSLSMgMtLLLvLLvL (37)式中的SL、SL及Mi、Mi为分别按段塞流和雾流计算的混合物密度及摩擦梯度。(4)雾流雾流混合物密度计算公式与泡流相同： ggLgggLLmHHHH)1 (由于雾流的气液无相对运动速度，即滑脱速度接近于雾，基本上没有滑脱。所以gLggqqqH

15、 (38)摩擦梯度则按连续的气相进行计算，即Dvfsggf22 (39)式中 sgv气体表观流速，pgsgAqv/ ，m/s。雾流摩擦系数可根据气体雷诺数gNRe和液膜相对粗糙度由图 2 查得。按不同流动型态计算压力梯度的步骤与前面介绍的用摩擦损失系数法基本相同，只是在计算混合物密度及摩擦之前需要根据流动型态界限确定其流动型态。图 3 为 Orkiszewski 方法的计算流程框图。12以井口油压或井底流压为起点，选择合适的压力间隔P，假设h计算平均 P 和 T，并求得在此 P 和 T 下的流体性质参数和流动参数，以及相应的流动型态界限LB、Lg和 LM确定流动型态雾流计算气相存容比、平均密度

16、及摩擦梯度过渡流分别按段塞流和雾流计算平均密度及摩擦梯度，并进行内插段塞流计算滑脱速度、液体分布系数、平均密度和摩擦梯度泡流计算气相存容比、平均密度和摩擦梯度计算并比较h，重复上述计算使h的计算值与假设值相等或在允许的误差范围内重复上述步骤，直到h的等于或大于油层深度为止图 3 Orkiszewski 方法计算流程框图悬点载荷和抽油杆柱设计计算悬点载荷和抽油杆柱设计计算 抽油杆柱设计的一般方法见采油工程设计与原理 。之所以设计方法较复杂，原因之一是因为杆柱的最大、最小载荷与杆长不是线性关系。例如在考虑抽油杆弹性时的悬点载荷、在考虑杆柱摩擦时的悬点载荷公式与杆长不是线性关系。原因之二是因为杆、管

17、环空中的压力分布取决于杆径，而杆柱的设计有用到杆、管环空中的压力分布。 由于综合课程设计时间较少，所以这里提供一种简化杆柱设计方法。暂将杆、管环空中的压力分布给定（按油水两相、不考虑摩擦时的压力分布） ，杆柱的最大、最小载荷公式采用与杆长成线性关系的下面公式。它是针对液体粘度较低、直井、游梁抽油机的杆柱载荷公式。 悬点最大、最小载荷的计算公式：13 )17901)(21maxSNWWPLijrj （40） gLqWrjijrjijrj11 （41） )(NZpLPPfW （42） 式中：irq第 i 级杆每米杆在空气中的质量，Kg/m riL第 i 级杆杆长，m； i 抽油杆级数，从下向上计数

18、； PZ泵排出口压力，Pa；PN泵的沉没压力，Pa；N冲次，rpm;S光杆冲程，m；fP活塞截面积，m2；g重力加速度，m/s2；ijrjijrjWSNWP121min1790 （43）)(11111jrrjijjijrjijrjffPWW （44） 式中：令 fr0=0 Pj第 j 级抽油杆底部断面处压力,Pa： )()1 (110jttwwwtjLLgffPP （45） Pt井口油压，Pa; 0地面油密度，kg/m3;14 fw体积含水率，小数；应力范围比pL计算公式：minminmaxallPL （46） frPmaxmax frPminmin （47） 抽油杆柱的许用最大应力的计算公式

19、：SFTallmin)5625. 04(式中：all抽油杆许用最大应力，Pa； T抽油杆最小抗张强度，对 C 级杆，T=*108Pa,对 D 级杆T=*108Pa; min抽油杆最小应力，Pa； SF使用系数，考虑到流体腐蚀性等因素而附加的系数（小于或等于） ，使用时可考表 2 来选值。表 2 抽油杆的使用系数使用介质API D 级杆API C 级杆无腐蚀性矿化水含硫化氢若抽油杆的应力范围比小于pL则认为抽油杆满足强度要求，此时杆组长度可根据pL直接推导出杆柱长度的显示公式。对于液体粘度低的油井可不考虑采用加重杆，抽油杆自下而上依次增粗，所以应先给定最小杆径（19mm）然后自下而上依次设计。有

20、应力范围比的计算公式即给定的应力范围比（pL）计算第一级杆长 L1，若 L1 大于等于泵深L，则抽油杆为单级杆，杆长为 L，并计算相应的应力范围比，若 L1 小于泵深 L,则由应力范围比的计算公式及给定的应力范围比计算第二级杆长 L2，若 L2 大于15等于（L-L1）,则第二级杆长为 L2，并计算相应的应力范围比，若 L2 小于（L-L1）,则同理进行设计。在设计中若杆径为 25mm 仍不能满足强度要求，则需改变抽汲参数。在设计中若杆径小于或等于 25mm 并满足强度要求，则杆柱设计结束。此为杆柱非等强度设计方法。若采用等强度设计方法，则需降低pL重新设计杆的长度。在设计抽油杆的过程中油管直

21、径一般取212（外径 73mm，内径 62mm） 。若泵径大于或等于 70mm，则油管全用3 （外径 89mm, 内径 76mm） ，原因是作业时大柱塞不能下如小直径油管中；若采用 25mm 抽油杆，则相应油管直径应用3 ，原因是 25mm 抽油杆节箍为 55mm，与 62mm 油管间隙太小。当采用多级杆时3 油管长度比 25mm 杆长多 10m。为了减小计算工作量，在本次课程设计中杆柱设计简化处理，采用单级杆设计（19mm） 。设计内容如下：由于采用单级杆设计，且杆径为 19mm，所以选用油管的直径为：62mm。计算内容和步骤：最大载荷： )(NZpLPPfW = 6=11iirjrjrjj

22、jWq L g；由于是单级的计算，所以简化为： rW=rqrL g= 2619104=)17901)(21maxSNWWPLijrj maxP （+）（1+23 61790） = 2、 最小载荷：ijrjijrjWSNWP121min1790 161111()iiirjrjjrjrjjjjWWPff 式中：令 fr0=0. )()1 (110jttwwwtjLLgffPP 由于，在该设计过程之中，只有一级杆，所以公式变为： 0(1)twwwPPffgL =1+= Mpa0()rjrjrjrWWPff = （） 106= 2min1790rjrjSNPWW= - 23 626174.241790

23、= 抽油机校核抽油机校核1）最大扭矩计算公式 min)(202. 01800maxmaxPPSSM =1800 3+ 3 （ - ）= 2）电动机功率计算，143881000maxnMNt =1000 18849.55 614388= 所以，可知电机的计算功率小于电机的额定功率，因而符合要求。 泵效计算泵效计算 （1）泵效及其影响因素在抽油井生产过程中，实际产量 Q 一般都比理论产量 Qt 要低，两者的比值叫泵效， 表示，tQQ （50）（2）产量计算 根据影响泵效的三方面的因素，实际产量的计算公式为 lleaklPtBqBSSQQ （51）17式中：Q实际产量，m3/d; Qt理论产量，m3

24、/d; Sp柱塞冲程，m； S光杆冲程，m； SSp抽油杆柱和油管柱弹性伸缩引起冲程损失系数； Bl泵内液体的体积系数； 泵的充满系数； qleak检泵初期的漏失量，m3/d;1）理论排量计算 SNfQpt1440=1400 =37.80 m3/d 2）冲程损失系数SSp的计算根据静载荷和惯性载荷对光杆冲程的影响计算本设计按照油管未锚定计算。当油管未锚定时；)()21 (3322112trrrlPfLfLfLfLSEWuSS 由于只有一级抽油杆柱，所以公式简化为： 211(1)()2lPrtWSLuLSSEff=2110.14781499.912001200(1)()23 2.06 100.0

25、0149990.001521=式中：uL/a= 曲柄角速度，rad/s；N/30=6/30=； a声波在抽油杆柱中的传播速度，5100m/s； pflpinZlgfLfPPW)( =1 106 =PZ泵排出口压力，Pa；Pin泵内压力，Pa；当液体粘度较低时，可忽略泵吸入口压力，故PinPN；PN泵的沉没压力，Pa； fp、fr、ft活塞、抽油杆及油管金属截面积，m2；18 L抽油杆柱总长度，m； l液体密度，kg/m3; E钢的弹性模数，1011Pa； Lf动液面深度，m； L1、L2、L3每级抽油杆的长度，m； fr1、fr2 、fr3每级抽油杆的截面积，m23) 充满系数 的计算 RKR

26、11 = 式中：K泵内余隙比；取. R泵内气液比； 050)10()1)(TPZTPfRRRininWSP = 565(50 10) 0.6 10351.66 0.96(3 1010 )293= PR=50，m3(标)/m3;SR=10m3(标)/m3;inP=3M Pa;Wf=；P0=105Pa; T0=293K；inT273t=；Z= 4） 泵内液体的体积系数 Bl 00(1)(1)lwwwwwBBfB fBff = 5）漏失量的计算 检泵初期的漏失量为)6(216003pleakDeVlPDeq=360.044 0.000051021600 (0.044 0.00005 0.6)6 0.

27、000531.50m3/d; D=；=s；l= 1.5m;PPZPN=106Pa；g=s2；e= ;30SNVp=3 630=0.6m/s;19所以最终算出泵的效率：()/leakPttllqSQQQSBB=(37.80 1.018 0.4814 1.0462)/37.80= 举升效率计算举升效率计算 光杆功率：P 光= rWSN/60 = 60=水力功率：P 水力Q 实际（PZPN）/=686400= 井下效率： 井下P 水力 / P 光 = 地面效率： 地P 光/ P 电机 = 系统效率： 总P 地* P 井下 = 20三、注水措施建议三、注水措施建议1、水质要求：油田注水所要求的水源不仅

28、量大，而且希望水源的水量和水质较为稳定。这样，在水源充足的地方，有个水源选择问题；水源缺乏的地方，需要寻找水源并进行选择。陆地水源包括地面的江、湖，泉水和地层水。海上包括海水和通过海底浅井抽取海水。水源选择要考虑到水质处理工艺要简便，还要满足油田注水设计的最大注水量。水源水量的估计以设计注水量为依据，如果采出的污水大部分回注的话，最终所需要的水量，大致为注水油层孔隙体积的 150170%。（1）地面水源淡水河、湖、泉水已广泛用于注水。随着国家建设的发展，工农业对这种水源使用也愈来愈广，加上可能遇到自然干旱，对注水可能供不应求。所以使用这种水源一般要得到有关部门的批准。 另外，这种地面水源，特别

29、是小溪、泉的水量常是随着季节变化的，并且常常是高含氧，携带很多悬浮物和各种微生物，不同季节水质成分变化很大，从而给水质处理带来许多麻烦。 胜利油田注水所用的黄河水属这种水源，特点是：有大量的泥沙和杂质，其含量在 200ppm 以上，并随季节变化；矿化度不高，一般在500600ppm；属于硫酸水型；含铁在左右。因此，黄河水要经过沉淀、过滤、杀菌和脱氧处理才能使用。（2）来自河床等冲积层水源淡水这种水是通过在河床打一些浅井到冲积层的顶部，从而使水质得到一定的改善。其特点是：水量稳定，水质变化不大，通常无腐蚀性；由于自然过滤，混浊度不受季节影响；水中含氧稳定便于处理，但由于硫酸还原菌深埋地下，这种水

30、仍可能受到它的污染。因此，把井钻深一些，以便排除或减少这种细菌的影响。（3）地层水水源淡水或盐水地层水源是根据地质资料，通过钻专门的水井而找到的来自地下的水源。找到高压、高产量的淡水层最好，盐水层也行，若找不到单一水层，多层水层也可以，但应注意，不同水层的水彼此不要产生化学反应而结垢。盐水也有它21的好处，可以防止注水所引起的粘上膨胀。（4）海水盐水近海和海上油田注水，一般用海水。因它既多又方便，但因高含氧和盐，腐蚀性强。悬浮的固体颗粒随季节变化较大，为改善这一点，通常钻一些浅井到海底，使其过滤从而减少水的机械杂质。2、水质处理：在水源确定的基础上，一般要进行水质处理。从防止设备腐蚀及地层堵塞的角度，对水质提出基本要求。主要是机械杂质堵塞、电化学腐蚀及生成物堵塞，细菌腐蚀及堵塞。 机械杂质的含量根据地层性质来定。国外对低渗透的孔隙性地层，含量要求小于，高渗透层或裂缝性地层含量可达 10ppm。机械杂质颗粒大小，一般