《采油工程课程设计沐风书屋》由会员分享，可在线阅读，更多相关《采油工程课程设计沐风书屋（24页珍藏版）》请在金锄头文库上搜索。
1、采油工程课程设计课程设计姓名：孔令伟 学号：201301509287中国石油大学（北京）石油工程学院2014年10月30日一、给定设计基础数据：1二、设计计算步骤22.1 油井流入动态计算22.2 井筒多相流的计算32.3悬点载荷和抽油杆柱设计计算112.4 抽油机校核152.5 泵效计算152.6举升效率计算18三、设计计算总结果19四、课程设计总结20一、给定设计基础数据： 井深：2000+8710=2870m 套管内径：0.124m 油层静压：2870/1001.2 =34.44MPa 油层温度：90 恒温层温度：16 地面脱气油粘度：30mPa.s 油相对密度：0.84 气相对密度：0
2、.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 曲线表示了油层工作特性。因而，它既是确定油井合理工作方式的依据，也是分析油井动态的基
3、础。本次设计油井流入动态计算采用Petro bras方法Petro bras方法计算综合IPR 曲线的实质是按含水率取纯油IPR曲线和水IPR曲线的加权平均值。当已知测试点计算采液指数时，是按产量加权平均;预测产量时,按流压加权平均。(1) 采液指数计算已知一个测试点：、和饱和压力及油藏压力。因为，=30/(34.44-12)= 1.3/( d.Mpa)(2) 某一产量下的流压Pwf =j()=1.4 x（34.44-10）=34.22t/d=+=34.44+1.4*10/1.8=42.22t/d-油IPR曲线的最大产油量。当0q时，令q=10 t/d，则p=15.754 Mpa同理，q=20
4、 t/d，P=13.877 Mpa q=30 t/d，P=12.0 Mpa当qq时，令q=50 t/d,则按流压加权平均进行推导得：P=f+0.125(1-f)P-1+=8.166Mpa 同理q=60t/d，P=5.860 Mpa当qq时，令q=71t/d，P=2.233 Mpa综上，井底流压与产量的关系列表如下：Pwf/Mpa15.74713.87312.010.08.1665.8602.233Q/(t/d)10203040.653506071 得到油井的流入动态曲线如下图：图1 油井IPR曲线2.2 井筒多相流的计算井筒多相流压力梯度方程井筒多相管流的压力梯度包括：因举高液体而克服重力所需
5、的压力势能、流体因加速而增加的动能和流体沿管路的摩阻损失，其数学表达式如下：gsin+v/d*式中为多相混合物的密度;v为多相混合物的流速;f为多相混合物流动时的摩擦阻力系数;d为管径;p为压力;h为深度;g为重力加速度; 为井斜角的余角。井筒多相管流计算包括两部分：（1）由井底向上计算至泵入口处；（2）油管内由井口向下计算至泵出口处。1）由井底向上计算至泵入口处，计算下泵深度Lp。采用深度增量迭代方法，首先估算迭代深度。在本设计中为了减小工作量，采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用Orkiszewski方法判断流型，进行压力梯度计算，最后计算出
6、深度增量和下泵深度Lp。按深度增量迭代的步骤:井底流压12Mpa，假设压力降为0.2 Mpa；估计一个对应的深度增量=40m，即深度为1960m 。由井温关系式可以计算得到该处的井温为：89.96。平均的压力和温度：=（90+89.96）/2=89.98。平均压力=11.9 Mpa。由平均压力和平均温度计算的得到流体的物性参数为：溶解油气比R=71.31 ； 原油体积系数B=1.25 原油密度P=739.00； 油水混合液的密度P =843.40； 死油粘度=6.537*10； 活油粘度=3.318*10； 水的粘度=3.263*10； 液体的粘度= 3.296*10；天然气的压缩因子Z=0.
7、9567； 天然气的密度90.70。以上单位均是标准单位。由以上的流体物性参数判断流型：不同流动型态下的和的计算方法不同，为此，计算中首先要判断流动形态。该方法的四种流动型态的划分界限如表1所示。表1 流型界限流动型态界 限泡 流段 塞 流过 渡 流雾 流其中=1.071-0.7277且0.13（如果0.13，则取=0.13）； =50+36 ； =75+84 （）。由计算得到，由于该段的压力大于饱和压力的值，所以该段的流型为纯液流。计算该段的压力梯度。由压力梯度的计算公式： =843.40；=计算对应于的该段管长(深度差)。 将第步计算得的与第步估计的进行比较，两者之差超过允许范围，则以新的
8、作为估算值，重复的计算，使计算的与估计的之差在允许范围内为止。该过程之中只迭代一次。2）由井口向下计算至泵出口处，计算泵排出口压力PZ。采用压力增量迭代方法，首先估算迭代压力。同样为了减小工作量，也采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用Orkiszewski方法判断流型，进行压力梯度计算，最后计算出压力增量和泵排出口压力PZ。按压力增量迭代的步骤已知任一点(井底或井口)的压力， 选取合适的深度间隔(可将管等分为n段)。估计一个对应于计算间隔的压力增量。计算该段的和 ，以及、下的流体性质参数。计算该段压力梯度计算对应于的压力增量比较压力增量的估计量与
9、计算值 ，若二者之差不在允许范围内，则以计算值作为新的估计值，重复第步，使两者之差在允许范围之内为止。计算该段下端对应的深度和压力 以处的压力为起点压力重复第步，计算下一段的深度和压力 ，直到各段累加深度等于或大于管长时为止。2.2计算气-液两相垂直管流的Orkiszewski方法本设计井筒多相流计算采用Orkiszewski方法。Orkiszewski法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。如图1所示。在处理过渡性流型时，采用内插法。在计算段塞流压力梯度时要考虑气相与液体的分布关系。针对每种流动型态提出了存容比及摩擦损失的计算方法。图1 气液混合物流动型态(Orkiszewski
10、) 1.压力降公式及流动型态划分界限由前面垂直管流能量方程可知，其压力降是摩擦能量损失、势能变化和动能变化之和。由式(2-36)可直接写出多项垂直管流的压力降公式： (26)式中 压力，Pa； 摩擦损失梯度，Pa/m； 深度，m； 重力加速度，m/s2； 混合物密度，kg/m3； 混合物流速，m/s。动能项只是在雾流情况下才有明显的意义。出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为： (27)式中 管子流通截面积，m2； 流体总质量流量，kg/s； 气体体积流量，m3/s。将式(27)代入式(26)，并取，经过整理后可得： (28)式中 计算管段压力降，Pa； 计
11、算管段的深度差，m； 计算管段的平均压力，Pa。不同流动型态下的和的计算方法不同,下面按流型分别介绍。(1)泡流平均密度 式中 气相存容比(含气率)，计算管段中气相体积与管段容积之比值； 液相存容比(持液率)，计算管段中液相体积与管段容积之比值； 在下气、液和混合物的密度，kg/m3。气相存容比由滑脱速度来计算。滑脱速度定义为：气相流速与液相流速之差。 可解出： 式中 滑脱速度，由实验确定，m/s； 、气相和液相的表观流速，m/s。泡流摩擦损失梯度按液相进行计算： 式中 摩擦阻力系数； 液相真实流速，m/s。摩擦阻力系数可根据管壁相对粗造度和液相雷诺数查图2。液相雷诺数: 式中 在下的液体粘度
12、，油、水混合物在未乳化的情况下可取其体积加权平均值，Pa.s。图 2(2)段塞流混合物平均密度 (34)式中 液体分布系数； 滑脱速度，m/s。滑脱速度可用Griffith和Wallis提出的公式计算： (35) (3)过渡流过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。 (36) (37)式中的、及、为分别按段塞流和雾流计算的混合物密度及摩擦梯度。(4)雾流雾流混合物密度计算公式与泡流相同： 由于雾流的气液无相对运动速度，即滑脱速度接近于雾，基本上没有滑脱。所以 (38)摩擦梯度则按连续的气相进行计算，即 (39)式中 气体表观流速， ，m/s。雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度由图2查得。按不同流动型态计算压力梯度的步骤与前面介绍的用摩擦损失系数法基本相同，只是在计算混合物密度及摩擦之前需要根据流动型态界限确定其流动型态。图3为Orkiszewski方法的计算流程框图。图3 Orkiszewski方法计算流程框图2.3悬
