采油工程课程设计

采油工程课程设计

课程设计姓名：孔令伟学号：201301509287中国石油大学（北京）石油工程学院2014年10月30日一、给定设计基础数据：TOC\o"1-5"\h\z\o"CurrentDocument"二、设计计算步骤 2\o"CurrentDocument"2.1油井流入动态计算 2\o"CurrentDocument"2.2井筒多相流的计算 4\o"CurrentDocument"悬点载荷和抽油杆柱设计计算 11\o"CurrentDocument"抽油机校核 15\o"CurrentDocument"2.5泵效计算 15\o"CurrentDocument"2.6举升效率计算 17\o"CurrentDocument"三、注水措施建议 18\o"CurrentDocument"四、设计计算总结果 20\o"CurrentDocument"五、课程设计总结 21一、给定设计基础数据：井深：2000+87X10=2870m套管内径：0.124m油层静压：2870/100X1.2=34.44MPa油层温度：90°C恒温层温度：16C地面脱气油粘度：30mPa.s油相对密度：0.84气相对密度：0.76水相对密度：1.0油饱和压力：10MPa含水率：0.4套压：0.5MPa油压：1MPa生产气油比：50m3/m3原产液量（测试点）：30t/d原井底流压(测试点)：16.35Mpa抽油机型号：CYJ10353HB电机额定功率：37kw配产量：50t/d泵径：56mm冲程：3m冲次：6rpm柱塞与衬套径向间隙:0.3mm沉没压力：3MPa二、设计计算步骤2.1油井流入动态计算油井流入动态是指油井产量与井底流动压力的关系，它反映了油藏向该井供油的能力。从单井来讲，IPR曲线表示了油层工作特性。因而，它既是确定油井合理工作方式的依据，也是分析油井动态的基础。本次设计油井流入动态计算采用Petrobras方法Petrobras方法计算综合IPR曲线的实质是按含水率取纯油IPR曲线和水IPR曲线的加权平均值。当已知测试点计算采液指数时，是按产量加权平均;预测产量时,按流压加权平均。采液指数计算已知一个测试点：Pwftest>qtxest和饱和压力及油藏压力。q txwstpP-P因为wftest，= =30/(34.44-12)=1.3/(d.Mpa)某一产量下的流压Pwf=j(p1-pb)=1.4x(34.44-10)=34.22t/djbqomzx=+1.8=34.44+1.4*10/1.8=42.22t/dQomzx-油IPR曲线的最大产油量。P-当0q时，令q=10t/d，则p=1qt-jP-当0q时，令q=10t/d，则p=1qt-j=15.754Mpa同理，q=20t/d,P=13.877Mpaq=30t/d，P=12.0Mpa当qqqomz时，令q=50t/d,则按流压加权平均进行推导得:（P-红） 初1-P=fw1j+0.125（1-f）P[-1+ qomzx-qb=8.166Mpa同理q=60t/d，P=5.860Mpa当qqt时，p_f（二q、（q-q ）（8f-9）p—J（p omzx） 1 omzx w wfw1omzx令q=71t/d，P=2.233Mpa综上，井底流压与产量的关系列表如下:Pwf/Mpa15.747 13.873 12.0 10.0 8.166 5.860 2.233Q/（t/d） 10 2030 40.653 50 60 71得到油井的流入动态曲线如下图:IPR曲线图IPR曲线图1油井IPR曲线2.2井筒多相流的计算井筒多相流压力梯度方程井筒多相管流的压力梯度包括：因举高液体而克服重力所需的压力势能、流体因加速而增加的动能和流体沿管路的摩阻损失，其数学表达式如下：dpdvv2= m+f ―m—dhPgsin9+PvdhmP/d*2式中P为多相混合物的密度；v为多相混合物的流速;f为多相混合物流动时的摩擦阻力系数;d为管径;p为压力;h为深度;g为重力加速度；9为井斜角的余角。井筒多相管流计算包括两部分：（1）由井底向上计算至泵入口处；（2）油管内由井口向下计算至泵出口处。1）由井底向上计算至泵入口处，计算下泵深度Lp。采用深度增量迭代方法,首先估算迭代深度。在本设计中为了减小工作量，采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用Orkiszewski方法判断流型，进行压力梯度计算，最后计算出深度增量和下泵深度Lp。按深度增量迭代的步骤:井底流压12Mpa,假设压力降为0.2Mpa；估计一个对应的深度增量=40m,即深度为1960m。由井温关系式可以计算得到该处的井温为：89.96°C。平均的压力和温度：=（90+89.96）/2=89.98C。平均压力=11.9Mpa。由平均压力和平均温度计算的得到流体的物性参数为：溶解油气比R=71.31；原油体积系数B=1.25原油密度P=739.00；油水混合液的密度P=843.40；死-4 -4油粘度口=6.537*10 ；活油粘度口O=3.318*10 ； 水的粘度口-4 -4=3.263*10；液体的粘度口=3.296*10；天然气的压缩因子Z=0.9567；天然气的密度卩g=90.70。以上单位均是标准单位。由以上的流体物性参数判断流型：不同流动型态下的和的计算方法不同,为此,计算中首先要判断流动形态。该方法的四种流动型态的划分界限如表1所示。

表1流型界限流动型态界限泡 流■qg<Lq bt段塞流q>L,v<Lq bgst过渡流L>v>LM g S雾 流v>Lg M其中=1.071-0.7277t2/D且>0.13(如果<0.13,则取=0.13)；vq—g=50+36qt；q0.75—g0.75=75+84(qt)由计算得到,由于该段的压力大于饱和压力的值,所以该段的流型为纯液流。dp计算该段的压力梯度dh。由压力梯度的计算公式：-dP=tdh+gpdh+pvdvfmmmm=843.40；=计算对应于的该段管长（深度差）从计。将第步计算得的从计与第②步估计的进行比较，两者之差超过允许范围，则以新的作为估算值，重复②〜⑤的计算，使计算的与估计的之差在允许范围内为止。该过程之中只迭代一次。2）由井口向下计算至泵出口处，计算泵排出口压力PZ。采用压力增量迭代方法，首先估算迭代压力。同样为了减小工作量，也采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用Orkiszewski方法判断流型，进行压力梯度计算，最后计算出压力增量和泵排出口压力PZ。按压力增量迭代的步骤

已知任一点（井底或井口）的压力，选取合适的深度间隔（可将管等分为n段）。估计一个对应于计算间隔的压力增量。计算该段的和，以及、下的流体性质参数。④计算该段压力梯度dP④计算该段压力梯度dP\~dh丿o⑤计算对应于的压力增量比较压力增量的估计量与计算值，若二者之差不在允许范围内，则以计算值作为新的估计值，重复第②〜⑤步，使两者之差在允许范围之内为止。计算该段下端对应的深度和压力P=P+卫APL—ixAh«o «i1以处的压力为起点压力重复第②〜⑦步，计算下一段的深度和压力，直到各段累加深度等于或大于管长时为止。2.2计算气-液两相垂直管流的Orkiszewski方法本设计井筒多相流计算采用Orkiszewski方法。Orkiszewski法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。如图1所示。在处理过渡性流型时，采用内插法。在计算段塞流压力梯度时要考虑气相与液体的分布关系。针对每种流动型态提出了存容比及摩擦损失的计算方法。，屯，屯图1气液混合物流动型态（Orkiszewski）1.压力降公式及流动型态划分界限

由前面垂直管流能量方程可知，其压力降是摩擦能量损失、势能变化和动能变化之和。由式(2-36)可直接写出多项垂直管流的压力降公式：—dP=tdh+gpdh+pvdvfmmmm(26)式中一压力，Pa；—摩擦损失梯度，Pa/m；一深度，m；一重力加速度，m/s2；一混合物密度，kg/m3；一混合物流速，m/s。动能项只是在雾流情况下才有明显的意义。出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：pvdvmmmWqpvdvmmmWqXPdp(27)式中一管子流通截面积，m2一流体总质量流量，kg/s；一气体体积流量，m3/s。将式(27)代入式(26),并取dh=-Ahk,dP=APk，Pm=Pm，P=P经过整理后可得:m1—Wqm1—Wqt看A2P]Ahk(28)式中一计算管段压力降，Pa；一计算管段的深度差，一计算管段的深度差，m；一计算管段的平均压力，Pa。不同流动型态下的和的计算方法不同,下面按流型分别介绍。泡流平均密度p=Hp+Hp-H卞+Hpm LL gg gL ggH+H=1Lg式中一气相存容比(含气率)，计算管段中气相体积与管段容积之比值一液相存容比(持液率)，计算管段中液相体积与管段容积之比值；

Pg、PjPm一在P、T下气、液和混合物的密度，kg/m3。可解出：气相存容比由滑脱速度来计算。滑脱速度定义为：气相流速与液相流速之差可解出：vv q q-qvs-Hg1-H—aHa(1-H)g pg p gH二2[1+VAsp式中一滑脱速度，由实验确定，m/s、一气相和液相的表观流速，m/s。泡流摩擦损失梯度按液相进行计算：vLHvLHA(1-H)pg式中一摩擦阻力系数；一液相真实流速，m/s。摩擦阻力系数可根据管壁相对粗造度£/D和液相雷诺数查图2。N=DVA液相雷诺数：Re巴式中一在P、T下的液体粘度，油、水混合物在未乳化的情况下可取其体积加权平均值，Pa.s。

an岂冠)，_.己芒丄兰岂.--.豊==亘.-:..11V葺-Q5

雷-:

甘.=

岂an岂冠)，_.己芒丄兰岂.--.豊==亘.-:..11V葺-Q5

雷-:

甘.=

岂.0L5-nsuw-0BgF档-H驚二百心^!；Zc.二HSpew=-<■!:5二2-5「2?「管=*-*驾，=

go

s

£wEwd

S.UT

-二

lnT=0段塞流混合物平均密度W+PVA。p—tl~sp+opmq+vA Ltsp式中—液体分布系数；一滑脱速度，m/s。滑脱速度可用Griffith和Wallis提出的公式计算:(35)v—CC^[gD(35)s 1 2过渡流过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。L—v

pL—v

p=〒 -pmL—LSLMSv—L

卡 pMSMi(36)(37)L一vv—LT —T+ ~T(37)tL—LSLL—LMiMS MS式中的、及、为分别按段塞流和雾流计算的混合物密度及摩擦梯度。雾流雾流混合物密度计算公式与泡流相同p—Hp+Hp—(1—H)p+HpmLLgg gLgg由于雾流的气液无相对运动速度，即滑脱速度接近于雾，基本上没有滑脱。所以H=-^t~gqL+忍 (38)摩擦梯度则按连续的气相进行计算，即T=fp_E2_Tf—f2D (39)式中—气体表观流速，vsg—qg/Ap，m/s。雾流摩擦系数可根据气体雷诺数Z和液膜相对粗糙度由图2查得。按不同流动型态计算压力梯度的步骤与前面介绍的用摩擦损失系数法基本相同，只是在计算混合物密度及摩擦之前需要根据流动型态界限确定其流动型态。图3为Orkiszewski方法的计算流程框图。

图3Orkiszewski方法计算流程框图2.3悬点载荷和抽油杆柱设计计算抽油杆柱设计的一般方法见《采油工程设计与原理》。之所以设计方法较复杂，原因之一是因为杆柱的最大、最小载荷与杆长不是线性关系。例如在考虑抽油杆弹性时的悬点载荷、在考虑杆柱摩擦时的悬点载荷公式与杆长不是线性关系。原因之二是因为杆、管环空中的压力分布取决于杆径，而杆柱的设计有用到杆、管环空中的压力分布。由于综合课程设计时间较少，所以这里提供一种简化杆柱设计方法。暂将杆、管环空中的压力分布给定（按油水两相、不考虑摩擦时的压力分布），杆柱的最大、最小载荷公式采用与杆长成线性关系的下面公式。它是针对液体粘度较低、直井、游梁抽油机的杆柱载荷公式。悬点最大、最小载荷的计算公式：

i SN2Pmax-孑“+WL）（1+鬲a）》W-》qLgj—1 rj j—1 rjrj （41）W'f（PZ-PN）（42）式中：——第i级杆每米杆在空气中的质量，Kg/m――第i级杆杆长，m；i――抽油杆级数，从下向上计数;PZ——泵排出口压力，Pa；PN——泵的沉没压力，Pa；N 冲次，rpm;S 光杆冲程，m；fP 活塞截面积，m2；g 重力加速度，m/s2；Pmin—》W，-丝工WPminj—1rj1790j—1rj（43）rjj—1—》W.-》p.x（f-f）rjj—1j—1 j—1 （44）式中：令式中：令fr0=0Pj――第j级抽油杆底部断面处压力,Pa：P—P+[Px（1-f）+pf]•g（L-艺L）jt0www t—1t（45）Pt——井口油压，Pa;p 地面油密度，kg/m3;fw——体积含水率，小数；应力范围比计算公式：cmax-cminPL—call-cmin（46）PmaxPminbmax= bmm=TOC\o"1-5"\h\z方 f （47）T ——ball=（—+0.5625bmin）SF抽油杆柱的许用最大应力的计算公式： 4式中：ball——抽油杆许用最大应力，Pa；T——抽油杆最小抗张强度，对C级杆，T=6.3*108Pa,对D级杆T=8.1*108Pa;bmin——抽油杆最小应力，Pa；——使用系数，考虑到流体腐蚀性等因素而附加的系数（小于或等于1.0）使用时可考表2来选值。表2抽油杆的使用系数使用介质APID级杆APIC级杆无腐蚀性1.001.00矿化水0.900.65含硫化氢0.700.50若抽油杆的应力范围比小于［］则认为抽油杆满足强度要求，此时杆组长度可根据［］直接推导出杆柱长度的显示公式。对于液体粘度低的油井可不考虑采用加重杆，抽油杆自下而上依次增粗，所以应先给定最小杆径（19mm）然后自下而上依次设计。有应力范围比的计算公式即给定的应力范围比（［］=0.85）计算第一级杆长L1,若L1大于等于泵深L,则抽油杆为单级杆，杆长为L,并计算相应的应力范围比，若L1小于泵深L,则由应力范围比的计算公式及给定的应力范围比计算第二级杆长L2,若L2大于等于（L-L1）,则第二级杆长为L2,并计算相应的应力范围比，若L2小于（L-L1）,则同理进行设计。在设计中若杆径为25mm仍不能满足强度要求，则需改变抽汲参数。在设计中若杆径小于或等于25mm并满足强度要求，则杆柱设计结束。此为杆柱非等强度设计方法。若采用等强度设计方法，则需降低［］重新设计杆的长度。在设计抽油杆的过程中油管直径一般取J'2（外径73mm,内径62mm）。

若泵径大于或等于70mm,则油管全用（外径89mm,内径76mm）,原因是作业时大柱塞不能下如小直径油管中；若采用25mm抽油杆，则相应油管直径应用，原因是25mm抽油杆节箍为55mm,与62mm油管间隙太小。当采用多级杆时油管长度比25mm杆长多10m。为了减小计算工作量，在本次课程设计中杆柱设计简化处理，采用单级杆设计（19mm）。设计内容如下：由于采用单级杆设计，且杆径为19mm，所以选用油管的直径为：62mm。计算内容和步骤：最大载荷：TOC\o"1-5"\h\z\o"CurrentDocument"W'=f(P-P) 6WLfp(ZN丿=0.0014999110=1499.9N=》qLg由于是单级的计算，所以简化为：rjrjrj由于是单级的计算，所以简化为：j=1 j=1WqLgr=rr=78509.81200192 x10-6一4=26174.24N\o"CurrentDocument"i SN2p=(Zw+w')(i+-)max rjL1790j=13x621790Pmax=(1499.9+26174.24)(1+ )=29343.66N2、最小载荷：P=Zw工WTOC\o"1-5"\h\zminrj1790rjj=1 j=1Yw、工W—工Px(f—f)rj rj jrj rj-1j=1 j=1 j=1式中：令fr0=0.P=P+[Px（1—f）+pf]•g（L—艺L）jt 0 w ww tt=1由于，在该设计过程之中，只有一级杆，所以公式变为P=P+[Px(1—f)+pf]•gL=1+10.631=11.631Mpat0www

=1+10.631=11.631MpaW'rj=Wrj-PX(fj-fr0)=26174.24N—11.631(0.0014999-0) 10=8728.9NP=W'-W x26174.24minrj1790rj=8728.9-1790 =7149.68N2.4抽油机校核1)最大扭矩计算公式Mmax二1800S+0,202S(Pmax-P斷)=18003+0.2023(29343.66-7149.68)=18849.55N.m2)电动机功率计算，1000M xn1000x18849.55x6N= max—t14388 = 14388 =7860.53W所以，可知电机的计算功率小于电机的额定功率，因而符合要求。2.5泵效计算1）泵效及其影响因素在抽油井生产过程中，实际产量Q一般都比理论产量Qt要低，两者的比值叫泵效，n叫泵效，n表示,50)2）产量计算根据影响泵效的三方面的因素，实际产量的计算公式为S0qX—PX -―leak51)51)ll式中：Q——实际产量，m3/d;Qt——理论产量，m3/d;Sp 柱塞冲程，m；S 光杆冲程，m；SpzS——抽油杆柱和油管柱弹性伸缩引起冲程损失系数;Bl——泵内液体的体积系数；

B――泵的充满系数；qleak 检泵初期的漏失量，m3/d;1）理论排量计算Qt二1440fpSN=14000.001499936=37.80m3/d2）冲程损失系数的计算根据静载荷和惯性载荷对光杆冲程的影响计算本设计按照油管未锚定计算。当油管未锚定时；S u2 WLLP=当油管未锚定时；S u2 WLLP=(1+ )- —(4+ 2S 2SEffr1 r2L L、+-^+)

f ftr3由于只有一级抽油杆柱，所以公式简化为：SP=(1+U2)-性(土+L(1+O」4782). 1499.9X( 1200 +1200)S2SEff=( 2 )3x2.06x1011(0.0014999 0.001521，r1t=1.018式中：u=sL/a=0.14783 曲柄角速度，rad/s；3=nN/30=n6/30=0.6283；a 声波在抽油杆柱中的传播速度，5100m/s；W1=（PZ-Pin）•fpXIXLfgfp=1100.0014999=1499.9NPZ——泵排出口压力，Pa；Pin——泵内压力，Pa；当液体粘度较低时，可忽略泵吸入口压力，故Pin-PN；PN——泵的沉没压力，Pa；fp、fr、ft 活塞、抽油杆及油管金属截面积，m2；L——抽油杆柱总长度，m；pl 液体密度，kg/m3;E——钢的弹性模数，2.06x1011Pa；Lf——动液面深度，m；L1、L2、L3——每级抽油杆的长度，m；fr1、fr2、fr3 每级抽油杆的截面积，m23）充满系数B的计算

1-KR1+R=0.4814式中：K——泵内余隙比；取0.1.R——泵内气液比；(R-R)(1-f)PTZPS W—in1-KR1+R=0.4814式中：K——泵内余隙比；取0.1.R——泵内气液比；(R-R)(1-f)PTZPS W—in-(P+105)T

in 0(50-10)x0.6x105x351.66x0.96(3x106+105)293=0.892=50，m3(标)/m3;=10m3(标)/m3;=3MPa;=0.4；P0=105Pa;T0=293K；=273+t=351.66；Z=0.964)泵内液体的体积系数BlB=B(1-f)+Bf沁B(1-f)+f-0462l0www0ww=1.04625)漏失量的计算检泵初期的漏失量为leak二21600(-De36卩21600x(兀x0.044x0.0000536x0.000531.5 0m3/d;D=0.044m；y=0.00053Pa•s；l=1.5m;AP~PZ一PN=10Pa；g=9.8m/s2；e=0.00005m;SN3x6V=—p30=30=0.6m/s;所以最终算出泵的效率：Q=(QxHx—-乩)/QtSBBt=(37.80X1.018x0.4814一1.0462)/37.80=4682%2.6举升效率计算光杆功率：P光=WrSN/60=8728.936/60=2618.67水力功率：P水力=Q实际(PZ—PN)/86.4=17.7110/86400=204.94井下效率：n井下=P水力/P光=0.0783地面效率：n地=卩光/P电机=0.3331系统效率：n总=卩地*P井下=0.0261三、注水措施建议1、水质要求：油田注水所要求的水源不仅量大，而且希望水源的水量和水质较为稳定。这样，在水源充足的地方，有个水源选择问题；水源缺乏的地方，需要寻找水源并进行选择。陆地水源包括地面的江、湖，泉水和地层水。海上包括海水和通过海底浅井抽取海水。水源选择要考虑到水质处理工艺要简便，还要满足油田注水设计的最大注水量。水源水量的估计以设计注水量为依据，如果采出的污水大部分回注的话，最终所需要的水量，大致为注水油层孔隙体积的150〜170%。（1） 地面水源——淡水河、湖、泉水已广泛用于注水。随着国家建设的发展，工农业对这种水源使用也愈来愈广，加上可能遇到自然干旱，对注水可能供不应求。所以使用这种水源一般要得到有关部门的批准。另外，这种地面水源，特别是小溪、泉的水量常是随着季节变化的，并且常常是高含氧，携带很多悬浮物和各种微生物，不同季节水质成分变化很大，从而给水质处理带来许多麻烦。胜利油田注水所用的黄河水属这种水源，特点是：有大量的泥沙和杂质，其含量在200ppm以上，并随季节变化；矿化度不高，一般在500〜600ppm；属于硫酸水型；含铁在0.5ppm左右。因此，黄河水要经过沉淀、过滤、杀菌和脱氧处理才能使用。（2） 来自河床等冲积层水源——淡水这种水是通过在河床打一些浅井到冲积层的顶部，从而使水质得到一定的改善。其特点是：水量稳定，水质变化不大，通常无腐蚀性；由于自然过滤，混浊度不受季节影响；水中含氧稳定便于处理，但由于硫酸还原菌深埋地下，这种水仍可能受到它的污染。因此，把井钻深一些，以便排除或减少这种细菌的影响。（3）地层水水源——淡水或盐水地层水源是根据地质资料，通过钻专门的水井而找到的来自地下的水源。找到高压、高产量的淡水层最好，盐水层也行，若找不到单一水层，多层水层也可以，但应注意，不同水层的水彼此不要产生化学反应而结垢。盐水也有它的好处，可以防止注水所引起的粘上膨胀。（4）海水——盐水近海和海上油田注水，一般用海水。因它既多又方便，但因高含氧和盐，腐蚀性强。悬浮的固体颗粒随季节变化较大，为改善这一点，通常钻一些浅井到海底，使其过滤从而减少水的机械杂质。2、水质处理：在水源确定的基础上，一般要进行水质处理。从防止设备腐蚀及地层堵塞的角度，对水质提出基本要求。主要是机械杂质堵塞、电化学腐蚀及生成物堵塞，细菌腐蚀及堵塞。机械杂质的含量根据地层性质来定。国外对低渗透的孔隙性地层，含量要求小于0.1〜0.5ppm,高渗透层或裂缝性地层含量可达lOppm。机械杂质颗粒大小，一般应小于岩石孔隙喉道的1／10。颗粒直径大于1／3喉道直径，对地层堵塞很快,但解诸也容易,固颗粒不能进入地层深部。国外用现场注入水通过孔径为0.45u的滤纸测定其半衰期（注水量降至1/2时，为半衰期），以半衰期的长短来对比水质，半衰期缩短了认为水质不合格。地面水和海水作为油田注入水时，一般都要经过除氧、杀菌处理，否则，会带来严重的后果。如加拿大帕宾那油田由于注入水未除氧、杀菌，结果天然气中由原