采油工程课程设计 YKK standardization office【 YKK5AB- YKK08- YKK2C- YKK18】
采油工程课程设计
课程设计
姓名:孔令伟
学号:
中国石油大学(北京)
石油工程学院
2014年10月30日
一、给定设计基础数据:
井深:2000+87×10=2870m
套管内径:
油层静压:2870/100× =
油层温度:90℃
恒温层温度:16℃
地面脱气油粘度:
油相对密度:
气相对密度:
水相对密度:
油饱和压力:10MPa
含水率:
套压:
油压:1 MPa
生产气油比:50m3/m3
原产液量(测试点):30t/d
原井底流压(测试点):
抽油机型号:CYJ10353HB 电机额定功率:37kw
配产量:50t/d
泵径:56mm
冲程:3m
冲次:6rpm
柱塞与衬套径向间隙:0.3mm 沉没压力:3MPa
二、设计计算步骤 油井流入动态计算
油井流入动态是指油井产量与井底流动压力的关系,它反映了油藏向该井供油的能力。从单井来讲,IPR 曲线表示了油层工作特性。因而,它既是确定油井合理工作方式的依据,也是分析油井动态的基础。本次设计油井流入动态计算采用Petro bras 方法Petro bras 方法计算综合IPR 曲线的实质是按含水率取纯油IPR 曲线和水IPR 曲线的加权平均值。当已知测试点计算采液指数时,是按产量加权平均;预测产量时,按流压加权平均。
(1) 采液指数计算 已知一个测试点:
wftest
P 、
txest
q 和饱和压力
b
P 及油藏压力P 。
因为wftest P ≥b P
,1j =
txwst
wfest
q P P -=30/= (
(2) 某一产量
t
q 下的流压Pwf
b
q =j(b P P
-1)= x ()=d m o zx q =b q +8
.1b
jP =+*10/=d omzx
q -油IPR 曲线的最大产油量。
当0?q t ?b q
时,令q 1t =10 t/d ,则p 1wf =
j q P t
-
1= Mpa
同理,q 2t =20 t/d ,P 2wf = Mpa q 3t =30 t/d ,P 3wf = Mpa 当q
b
?q t ?omzx q 时,令q 4t =50 t/d,则按流压加权平均进行推导得:
P 4wf =f )(1j q P t
w -+(1-f w )P b
= 同理q 5t =60t/d ,P 5wf = Mpa
当q
omzx
?q t
时,
1()(89)
()omzx t omzx w wf w q q q f p f p J J --=-
-
令q 6t =71t/d ,P 6wf = Mpa
综上,井底流压与产量的关系列表如下:
Pwf/Mpa Q/(t/d)
10
20
30
50
60
71
得到油井的流入动态曲线如下图:
图1 油井IPR 曲线
井筒多相流的计算
井筒多相流压力梯度方程
井筒多相管流的压力梯度包括:因举高液体而克服重力所需的压力势能、流体因加速而增加的动能和流体沿管路的摩阻损失,其数学表达式如下:
=
dh dp
ρm gsin θ+ρm v m m m f dh dv +ρm /d*22
m v
式中ρm 为多相混合物的密度;v m 为多相混合物的流速;f m 为多相混合物流动时的摩擦阻力系数;d 为管径;p 为压力;h 为深度;g 为重力加速度; θ为井斜角的余角。
井筒多相管流计算包括两部分:(1)由井底向上计算至泵入口处; (2)油管内由井口向下计算至泵出口处。
1)由井底向上计算至泵入口处,计算下泵深度Lp 。采用深度增量迭代方法,首先估算迭代深度。在本设计中为了减小工作量,采用只迭代一次的方法。计算井筒多相管流时,首先计算井筒温度场、流体物性参数,然后利用Orkiszewski 方法判断流型,进行压力梯度计算,最后计算出深度增量和下泵深度Lp 。
按深度增量迭代的步骤:
井底流压12Mpa ,假设压力降为 Mpa ;估计一个对应的深度增量h ?=40m ,即深度为1960m 。
由井温关系式可以计算得到该处的井温为:℃。
平均的压力和温度:T =(90+)/2=℃。平均压力P = Mpa 。由平均压力和平均温度计算的得到流体的物性参数为:溶解油气比R S = ; 原油体积系数B 0= 原油密度P 0=; 油水混合液的密
度P z =; 死油粘度μod =*104
-; 活油粘度μO
=*104-; 水的粘度μw =*104
-; 液体的粘度μ=
*104
-;天然气的压缩因子Z=; 天然气的密度
g ρ=
。以上单位均是标准单位。
由以上的流体物性参数判断流型:
不同流动型态下的m ρ和f τ的计算方法不同,为此,计算中首先要判断流动形态。该方法的四种流动型态的划分界限如表1所示。
表1 流型界限
其中
B
L =
2/t
D
且
B
L >(如果
B
L <,则取
B
L =);
S
L =50+36
g v g
t
q q ;
M
L =75+84 (g v
g
t q q )
0.75
。
由计算得到,由于该段的压力大于饱和压力的值,所以该段的流型为纯液流。
计算该段的压力梯度dh dP
。由压力梯度的计算公式:
m
m m m f dv v dh g dh dP ρ+ρ+τ=-
m ρ=;f τ=计算对应于P ?的该段管长(深度差)计h ?。
⑥ 将第步计算得的计h ?与第②步估计的h ?进行比较,两者之差超过允许范围,则以新的h ?作为估算值,重复②~⑤的计算,使计算的与估计的h ?之差在允许范围ε内为止。该过程之中只迭代一次。
2)由井口向下计算至泵出口处,计算泵排出口压力PZ 。采用压力增量迭代方法,首先估算迭代压力。同样为了减小工作量,也采用只迭代一次的方法。计算井筒多相管流时,首先计算井筒温度场、流体物性参数,然后利用Orkiszewski 方法判断流型,进行压力梯度计算,最后计算出压力增量和泵排出口压力PZ 。
按压力增量迭代的步骤
①已知任一点(井底或井口)的压力0P , 选取合适的深度间隔h ?(可将管L 等分为n 段)。 ②估计一个对应于计算间隔h ?的压力增量P ?。 ③计算该段的T 和P ,以及P 、T 下的流体性质参数。
④计算该段压力梯度o dh dP ?
?? ??
⑤计算对应于h ?的压力增量
o i dh dP h P ??? ???=? ⑥比较压力增量的估计量P ?与计算值i P ? ,若二者之差不在允许范围内,则以计算值作为新的估计值,重复第②~⑤步,使两者之差在允许范围o ε之内为止。
⑦计算该段下端对应的深度i L 和压力i P
⑧以i L 处的压力i P 为起点压力重复第②~⑦步,计算下一段的深度1+i L 和压力1+i P ,直到各段累加深度等于或大于管长L 时为止。 计算气-液两相垂直管流的Orkiszewski 方法
本设计井筒多相流计算采用Orkiszewski 方法。
Orkiszewski 法提出的四种流动型态是泡流、段塞流、过渡流及环雾流。如图1所示。在处理过渡性流型时,采用内插法。在计算段塞流压力梯度时要考虑气相与液体的分布关系。针对每种流动型态提出了存容比及摩擦损失的计算方法。
图1 气液混合物流动型态(Orkiszewski) 1.压力降公式及流动型态划分界限
由前面垂直管流能量方程可知,其压力降是摩擦能量损失、势能变化和动能变化之和。由式(2-36)可直接写出多项垂直管流的压力降公式:
m m m m f dv v dh g dh dP ρ+ρ+τ=- (26)
式中 P —压力,Pa ;
f τ—摩擦损失梯度,Pa/m ; h —深度,m ;
g —重力加速度,m/s2; m ρ—混合物密度,kg/m3; m v —混合物流速,m/s 。
动能项只是在雾流情况下才有明显的意义。出现雾流时,气体体积流量远大于液体体积流量。根据气体定律,动能变化可表示为:
dp P
A q W dv v p g
t m m m 2
-
=ρ (27)
式中 p A —管子流通截面积,m2; t W —流体总质量流量,kg/s ; g q —气体体积流量,m3/s 。
将式(27)代入式(26),并取k h dh ?-=,k p dP ?=,m m ρρ=,P P =经过整理后可得:
=
?k P k
p g
t f
m h P A q W g ?-+]1[
2τρ (28)
式中 k P ?—计算管段压力降,Pa ; k h ?—计算管段的深度差,m ; P —计算管段的平均压力,Pa 。
不同流动型态下的m ρ和f τ的计算方法不同,下面按流型分别介绍。 (1)泡流 平均密度
式中 g H —气相存容比(含气率),计算管段中气相体积与管段容积之比值; L H —液相存容比(持液率),计算管段中液相体积与管段容积之比值; m g ρρρ、、L—在T P 、下气、液和混合物的密度,kg/m3。
气相存容比由滑脱速度s V 来计算。滑脱速度定义为:气相流速与液相流速之差。
可解出g H :
H g =p
s g p s t p s t A v q A v q A v q 4)1(1[212-+-+
式中 s v —滑脱速度,由实验确定,m/s ; sg v 、sL v —气相和液相的表观流速,m/s 。 泡流摩擦损失梯度按液相进行计算: 式中 f —摩擦阻力系数; LH v —液相真实流速,m/s 。
摩擦阻力系数f 可根据管壁相对粗造度D /ε和液相雷诺数Re N 查图2。
液相雷诺数: L
L sL Dv N μρ=
Re
式中 L μ—在T P 、下的液体粘度,油、水混合物在未乳化的情况下可取其体积加权平均值,。
图 2
(2)段塞流混合物平均密度
L
p
s t p
s L t m A v q A v W δρ++ρ+=
ρ (34)
式中 δ—液体分布系数; s v —滑脱速度,m/s 。
滑脱速度可用Griffith 和Wallis 提出的公式计算:
gD
C C v s 21= (35)
(3)过渡流
过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算,然后用内插方法来确定相应的数值。
Mi
S
M s g SL S M g M m L L L
v L L v L ρ--+ρ--=ρ (36)
Mi
S
M g g SL S M g M t L L L
v L L v L τ--+τ--=
τ (37)
式中的SL ρ、SL τ及Mi ρ、Mi τ为分别按段塞流和雾流计算的混合物密度及摩擦梯度。 (4)雾流
雾流混合物密度计算公式与泡流相同:
由于雾流的气液无相对运动速度,即滑脱速度接近于雾,基本上没有滑脱。所以
g
L g g q q q H +=
(38)
摩擦梯度则按连续的气相进行计算,即
D
v f
sg
g f 22ρ=τ (39)
式中 sg v —气体表观流速,p g sg A q v /= ,m/s 。
雾流摩擦系数可根据气体雷诺数()g N Re 和液膜相对粗糙度由图2查得。
按不同流动型态计算压力梯度的步骤与前面介绍的用摩擦损失系数法基本相同,只是在计算混合物密度及摩擦之前需要根据流动型态界限确定其流动型态。图3为Orkiszewski 方法的计算流程框图。
图3 Orkiszewski 方法计算流程框图
悬点载荷和抽油杆柱设计计算
抽油杆柱设计的一般方法见《采油工程设计与原理》。之所以设计方法较复杂,原因之一是因为杆柱的最大、最小载荷与杆长不是线性关系。例如在考虑抽油杆弹性时的悬点载荷、在考虑杆柱摩擦时的悬点载荷公式与杆长不是线性关系。原因之二是因为杆、管环空中的压力分布取决于杆径,而杆柱的设计有用到杆、管环空中的压力分布。
由于综合课程设计时间较少,所以这里提供一种简化杆柱设计方法。暂将杆、管环空中的压力分布给定(按油水两相、不考虑摩擦时的压力分布),杆柱的最大、最小载荷公式采用与杆长成线性关系的下面公式。它是针对液体粘度较低、直井、游梁抽油机的杆柱载荷公式。
悬点最大、最小载荷的计算公式:
)
17901)((2
1max SN W W P L i
j rj +'+=∑= (40) g
L q W
rj i
j rj i
j rj
∑∑===1
1
(41)
)
(N Z p L P P f W -=' (42)
式中:i
r q ——第i 级杆每米杆在空气中的质量,Kg/m
ri
L ——第i 级杆杆长,m ;
i —— 抽油杆级数,从下向上计数; PZ ——泵排出口压力,Pa ; PN ——泵的沉没压力,Pa ; N ——冲次,rpm; S ——光杆冲程,m ; fP ——活塞截面积,m2; g ——重力加速度,m/s2;
∑∑==-'=i j rj i
j rj W SN W P 121min
1790 (43)
)
(111
1
1
-===-?-='∑∑∑j r rj i
j j i
j rj
i
j rj
f f P W
W (44)
式中:令fr0=0
Pj ——第j 级抽油杆底部断面处压力,Pa :
)
(])1([1
10∑-=-?+-?+=j t t w w w t j L L g f f P P ρρ (45)
Pt ——井口油压,Pa; ρ0——地面油密度,kg/m3; fw ——体积含水率,小数; 应力范围比pL 计算公式:
min min
max σσσσ--=
all PL (46)
fr P max max =
σ fr P min
min =σ (47)
抽油杆柱的许用最大应力的计算公式:
SF
T
all min)5625.04
(σσ+=
式中:all σ——抽油杆许用最大应力,Pa ;
T ——抽油杆最小抗张强度,对C 级杆,T=*108Pa,对D 级杆T=*108Pa; min σ——抽油杆最小应力,Pa ;
SF ——使用系数,考虑到流体腐蚀性等因素而附加的系数(小于或等于),使用时可考表2来选值。
表2 抽油杆的使用系数
若抽油杆的应力范围比小于[pL ]则认为抽油杆满足强度要求,此时杆组长度可根据[pL ]直接推导出杆柱长度的显示公式。
对于液体粘度低的油井可不考虑采用加重杆,抽油杆自下而上依次增粗,所以应先给定最小杆径(19mm )然后自下而上依次设计。有应力范围比的计算公式即给定的应力范围比([pL ]=)计算第一级杆长L1,若L1大于等于泵深L ,则抽油杆为单级杆,杆长为L ,并计算相应的应力范围比,若L1小于泵深L,则由应力范围比的计算公式及给定的应力范围比计算第二级杆长L2,若L2大于等于(L-L1),则第二级杆长为L2,并计算相应的应力范围比,若L2小于(L-L1),则同理进行设计。在设计中若杆径为25mm 仍不能满足强度要求,则需改变抽汲参数。在设计中若杆径小于或等于25mm 并满足强度要求,则杆柱设计结束。此为杆柱非等强度设计方法。若采用等强度设计方法,则需降低[pL ]重新设计杆的长度。
在设计抽油杆的过程中油管直径一般取
"
21
2(外径73mm ,内径62mm )。若泵径大于或等
于70mm ,则油管全用3''(外径89mm, 内径76mm ),原因是作业时大柱塞不能下如小直径油管中;若采用25mm 抽油杆,则相应油管直径应用3'',原因是25mm 抽油杆节箍为55mm ,与62mm 油管间隙太小。当采用多级杆时3''油管长度比25mm 杆长多10m 。
为了减小计算工作量,在本次课程设计中杆柱设计简化处理,采用单级杆设计(19mm )。 设计内容如下:
由于采用单级杆设计,且杆径为19mm ,所以选用油管的直径为:62mm 。 计算内容和步骤: 最大载荷:
)
(N Z p L P P f W -=' =??6
=
1
1
i
i
rj
rj rj j j W
q L g
===∑∑;由于是单级的计算,所以简化为:
r W =r q r L g =???26
19104π-??÷=
max P =
(+)?(1+
2
361790
?) =
2、 最小载荷: 式中:令fr0=0.
由于,在该设计过程之中,只有一级杆,所以公式变为:
0[(1)]t w w w P P f f gL
ρρ=+?-+? =1+= Mpa
'0()
rj rj rj r W W P f f =-?- = – ?()? 106
=
2'min
1790rj
rj SN P W W =-= - 2
3626174.24
1790??=
抽油机校核
1)最大扭矩计算公式
m in)
(202.01800max max P P S S M -+= =1800 ? 3+ ?3 ?( - )=
2)电动机功率计算,
143881000max n M N t ?=
=100018849.55614388??=
所以,可知电机的计算功率小于电机的额定功率,因而符合要求。
泵效计算
(1)泵效及其影响因素
在抽油井生产过程中,实际产量Q 一般都比理论产量Qt 要低,两者的比值叫泵效,η表示,
t Q Q
=
η (50)
(2)产量计算
根据影响泵效的三方面的因素,实际产量的计算公式为
l
leak
l P t B q B S S Q Q -??
=β (51)
式中:Q ——实际产量,m3/d; Qt ——理论产量,m3/d; Sp ——柱塞冲程,m ; S ——光杆冲程,m ;
S
S p ——抽油杆柱和油管柱弹性伸缩引起冲程损失系数;
Bl ——泵内液体的体积系数; β——泵的充满系数;
qleak ——检泵初期的漏失量,m3/d; 1)理论排量计算
SN
f Q p t 1440==1400 ? ??=37.80 m3/d
2)冲程损失系数S
S p 的计算
根据静载荷和惯性载荷对光杆冲程的影响计算 本设计按照油管未锚定计算。
当油管未锚定时;)
()21(3322112t r r r l P f L
f L f L f L SE W u S
S +++'
-+= 由于只有一级抽油杆柱,所以公式简化为:
211(1)()2l P r t W S L u L S SE f f '=+-+=2110.14781499.912001200(1)()
23 2.06100.00149990.001521+-?+??=
式中:u =ωL/a=
ω——曲柄角速度,rad/s ;ω=πN/30=π?6/30=; a ——声波在抽油杆柱中的传播速度,5100m/s ;
p
f l p in Z l gf L f P P W ?≈?-='
ρ)( =1? 106
? =
PZ ——泵排出口压力,Pa ;
Pin ——泵内压力,Pa ;当液体粘度较低时,可忽略泵吸入口压力,故Pin≈PN ; PN ——泵的沉没压力,Pa ;
fp 、fr 、ft ——活塞、抽油杆及油管金属截面积,m2; L ——抽油杆柱总长度,m ; ρl ——液体密度,kg/m3; E ——钢的弹性模数,×1011Pa ; Lf ——动液面深度,m ;
L1、L2、L3——每级抽油杆的长度,m ; fr1、fr2 、fr3——每级抽油杆的截面积,m2 3) 充满系数β的计算
R KR
+-=
11β =
式中:K ——泵内余隙比;取. R ——泵内气液比;
050)10()1)((T P Z T P f R R R in in W S P +--= = 565
(5010)0.610351.660.96(31010)293-?????+=
P R =50,m3(标)/m3;
S
R =10m3(标)/m3;
in P =3M Pa;
W
f =;P0=105Pa; T0=293K ;
in
T =273+t=;Z=
4) 泵内液体的体积系数Bl
00(1)(1)l w w w w w
B B f B f B f f =-+≈-+ =
5)漏失量的计算 检泵初期的漏失量为
)
6(216003p leak
DeV l P
De q πμπ-??==
36
0.0440.000051021600(
0.0440.000050.6)
60.00053
1.5
ππ????-????≈0m3/d;
D=;μ=·s ;l= 1.5m;
ΔP≈PZ —PN=106
Pa ;g=s2;e= ;
30SN V p =
=36
30?=0.6m/s;
所以最终算出泵的效率:
()/leak
P t t l l
q S Q Q Q S B B β=?
?-=(37.80 1.0180.4814 1.0462)/37.80??÷=﹪
举升效率计算
光杆功率:P 光= 'r W SN/60 = ??60=
水力功率:P 水力=Q 实际(PZ —PN )/=?6
86400= 井下效率:η井下=P 水力 / P 光 = 地面效率:η地=P 光/ P 电机 =
系统效率:η总=P 地* P 井下 =
三、注水措施建议
1、水质要求:油田注水所要求的水源不仅量大,而且希望水源的水量和水质较为稳定。这样,在水源充足的地方,有个水源选择问题;水源缺乏的地方,需要寻找水源并进行选择。陆地水源包括地面的江、湖,泉水和地层水。海上包括海水和通过海底浅井抽取海水。水源选择要考虑到水质处理工艺要简便,还要满足油田注水设计的最大注水量。水源水量的估计以设计注水量为依据,如果采出的污水大部分回注的话,最终所需要的水量,大致为注水油层孔隙体积的150~170%。
(1)地面水源——淡水
河、湖、泉水已广泛用于注水。随着国家建设的发展,工农业对这种水源使用也愈来愈广,加上可能遇到自然干旱,对注水可能供不应求。所以使用这种水源一般要得到有关部门的批准。
另外,这种地面水源,特别是小溪、泉的水量常是随着季节变化的,并且常常是高含氧,携带很多悬浮物和各种微生物,不同季节水质成分变化很大,从而给水质处理带来许多麻烦。
胜利油田注水所用的黄河水属这种水源,特点是:有大量的泥沙和杂质,其含量在200ppm 以上,并随季节变化;矿化度不高,一般在500~600ppm;属于硫酸水型;含铁在左右。因此,黄河水要经过沉淀、过滤、杀菌和脱氧处理才能使用。
(2)来自河床等冲积层水源——淡水
这种水是通过在河床打一些浅井到冲积层的顶部,从而使水质得到一定的改善。其特点是:水量稳定,水质变化不大,通常无腐蚀性;由于自然过滤,混浊度不受季节影响;水中含氧稳定便于处理,但由于硫酸还原菌深埋地下,这种水仍可能受到它的污染。因此,把井钻深一些,以便排除或减少这种细菌的影响。
(3)地层水水源——淡水或盐水