Pérdidas de Calor en superficie-Tubería desnuda. 1. Calculo Ts, Tsurf, Tp y rins. (Ts=115,1Ps0,225) Calculo de las pérdidas de calor Q.

Ts−Ta Volumen especifico

1. Temperatura de saturación (Tsurf=Ts) a. Tsurf = =℉
2 Vws=VsX+(1-X)Vw=X(Vs-Vw)+Vw
Ts=115,1 Ps0,225=℉ Tsurf −Ta Vws=XVs= pie3/lb
2. Propiedades físicas del aire a Tp. Tp= =℉
2 Vs=363,9(Ps-0,9588)=pie3/lb
Tsurf −Ta rins=rto+rins (rins espesor del aislante) Rango de presiones 15 a 1000 lpca
Tp= =℉
2 de+ 2 Vw=0,01602+0,000023(-6,6+3,74075Ps0,225+8,70394257P0,45)
de= =pie
Kha=0,01328+ 2,471 x 10−5 T −4,247 x 10−9 T 2=BTU /hpie ℉ 12 Calor sensible del agua Hw.
μa=0,04+6,155 x 10−5 T −1,22 x 10−8 T 2=lb/ pieh de Hw=Cw (Ts-32)=BTU/lb
rins=
BTU 2 Cw: calor especifico del agua 1BTU/lbF
Cpa=0,2382+1,39 x 10−5 T +1,027 x 10−8 T 2= ℉
lb Propiedades físicas del aire (Kha, a, Cpa, a, a) Hw=91,0 (P 0,2574)
−5 −8 2 3
ρa=0,771−8,848 x 10 T −3,744 x 10 T =lb / pie hr= BTU/hpie2F Rango de presiones 15 a 1000 lpca
1 2F
hc=BTU/hpie Calor latente de vaporización.
βa=0,0024−0,757 x 10−5 T + 0,169 x 10−7 T 2 −0,148 x 10−10 T 3=
℉ de transferencia de calor total Utins.
Coeficiente Lv=1318(Ps-0,08774)=BTU/lbF
−1
3. Coeficiente de transferencia de calor por Radiación hr rins Rango de presiones 15 a 1000 lpca
(Temperaturas en ºR)
2
1

σ =0,1713 x 10−8
2
hr =σε (T +T )( T 1 + T 2 )=BTU /h pie ℉
2
BTU
2 Utins= [
rins ∙ ln

Khins
rto

Pérdidas de Calor Q.
+
( ) 1
(
hc +hr ) ] =¿
BTU/hpie2F Calor total o entalpia de vapor seco y saturado.
Hs=Hw+Lv=BTU/lb
Hs=1119(Ps 0,01267)=BTU/lb
h pie 2 ℉ 4 Rango de presiones 15 a 1000 lpca
Q=2 πrins ∙Utins ( Ts−Ta ) ∆ L=BTU /h
4. Coeficiente de transferencia de calor por convección Calidad del vapor y vapor húmedo
Calculo de Tsurf=ºF
libre o natural hc=BTU/hpie2F Hws=Hw+xLv=BTU/lb
a) Se calcula U*
Calidad de vapor.
0,53 Kha gβa Cpa ∙ μa
−0,25 Método de separador
Khins
hc=
de [
de 3 (Tsurf −ta ) 2
va Kha ] U ¿=
rins
=BTU /h pie 2 ℉ X=mv/(mv+ml)
rins∙ ln ( )
rto
Método de los cloruros
8
g=4,17 x 10 pie /h 2 Cl−¿en el agua de alimentacion
X =100−% ¿
A=2 π ∙ rins ∆ L= pie 2 % Cl−¿en la parte liquida ∙ 100 ¿
5. Coeficiente de transferencia de calor total Uto. Método de la conductividad eléctrica
−1
Q %σ en el agua de alimentacion
Uto= ( hc 1+hr ) =BTU /h pie 2 ℉ Tsurf ¿ cal=Ts−
A U¿
=℉ X =100−
%σ en la parte liquida
∙ 100

6. Perdidas de Calor Q. =conductividad eléctrica

Q=2 πrt oUto ( Ts−Ta ) ∆ L=BTU /h Comparar Tsurf calculado con el supuesto. Si no son iguales Calidad del vapor en el cabezal del pozo Xst,wh.
dentro de una aprox 0,1º repetir el paso 1 usando el calculado QL
como el nuevo valor supuesto. (Iterar)
Xst , wh=Xst , gen−
WLv
Q=BTU/h, L=pies, W=lb/h, Lv=BTU/lb
Pérdidas de Calor-Tubería aislada. Iter. Tsurf Tsurf)cal hr hc Uto
 Marbellsy Nazaretth Sanchez Salazar Pr=
Cpan∙ μan K insf (t )
Tco=Te+ ( Ts−Tins )
Explotación de Crudos Pesados Khan rins
hc=
0,049 Khan ( GrPr )0,333 Pr 0,074 KheLn ( )
rto
rci 4. (Kha, a, Cpa, a, a) a Tp
Caída de presión en líneas de superficie.
fL V ws W 2
rtoLn
rto ( ) 5. hr donde T1=Tins y T2=Tci
∆ P=1,35731 x 10
−11
=lpc 6. Calculo de Gr y Pr para hc
d 5ti 7. Calculo de Uto. rto es rins y To es Tins
Factor de fricción rh −1

1 2

( 1
rtoLn
rco( ) )
f=

( −4 log { ε
3,7065
5,0452
− ℜ log
ε
2,8257
7,149
[
+ ℜ
1,1098

(
0,8981

) ]} )
Uto=
hc +hr
+

8. Calculo de Th.
Khcem 7. Calculo de Uto.

rins −1

( )
=1,5x10-4pies tuberías de acero ordinario
=5,0x10-6pies tuberías lisas o suaves
ℜ=0,5264166
w Th=
KheTe
Tsf ( t ) +
rtoUto
Khe
Uto= ( rtoLn
rto
Khins
+
rto
rins(hc+ hr) )
μws ∙ dti f ( t )+ 8. Calculo de Tins
rtoUto
( rins
Perdidas de calor en el pozo-Tubería desnuda.
rto )
1. Calculo de f(t) para tiempo=Dias 9. Calculo de Tco. rtoUtoLn
KheTe+ rtoUtof ( t ) Ts
f ( t )=ln ( ) 2 √ αt
rh
−0,29 rtoUtoLn ( rcorh ) ( Ts−Th ) Tins=Ts−
Khins (Ts−
Khe+rtoUtof (t ) )
Tco=Th+
Kh Khcem
α = = pie 2 /h Comparar Tco calculado y Tci supuesto si difieren de una
Comparar Tins) Supuesto con el Tins calculado si difieren de
M una aprox de 0,1º repetir desde el paso 3 usando el
M= (1-)rCr+(Sooco+SwwCw+SggCg) aprox de 0,1º repetir desde el paso 3 usando el calculado
calculado como el nuevo supuesto (Iterar)
Kh=BTU/hpieF, M=BTU/pie3F como el nuevo supuesto (Iterar)
Iter. Tins Th Tco hr hc Uto
2. Calculo de  Iter. Tci Th Tco hr hc Uto
9. Perdidas de calor Q.
−1 10. Tasa de pérdida de calor Q.
1 rto 1 Q=2 π ∙rtoUto ( Ts−Th ) ∆ L
ε=
[+
εto rci εci
−1 ( )] Q desde el interior de la tub de iny hasta cem-form
Q=2 π ∙rtoUto ( Ts−Th ) ∆ L Q=
2 πKhe (Ts−Tco)∆ L
f (t )
3. Tci)supuesto Q desde la interfase cem-form y formación
Ts+Tc 2 πKhe (Th−Te )∆ L 2 πKhins (Ts−Tins) ∆ L
Tci= =℉ Q= Q=
2 f (t) rins

Ts+Tci
Q a través del cemento
ln
rto ( )
Tp= =℉ 2 πKhcem (Tco−Th)∆ L Convección forzada  10 millas/h
2 Q=
4. Propiedades físicas del aire (Kha, a, Cpa, a, a)
5. hr donde T1=Ts y T2=Tci
ln
rh
( )
rco
log ( hcde
kha )
=−0,0757+0,3082 ( logRe )+ 0 ,=379(logRe ) 2

6. Calculo de Gr y Pr para hc Perdidas de calor en el pozo-Tubería aislada. 10008.800deva 50.000

1. Calculo de f(t) para tiempo=Dias (rh=rco) ( ρa∙ va)0,6
( rci−rto )3 g ρan2 βan(Ts−Tci) hc=19,3 Cpa
Gr= ( μan2 ) 2. Calculo de  rto es rins y to es ins
3. Caculo Tco con Tins)supuesto
de 0,4
8000 deva50.000
 Kha 5.280 ρa∙ va∙ de 0,805
b 13. Calculo de Hw,s
hc=0,0239
de ( μa ) ln ( μ1 ) =ln ( a ) +
T 1 + 460 H w =C w ( T ¿avg ¿−32 ) =BTU /lb
1metro-3,2808pies b H r =C w ( T r−32 ) =BTU /lb
ln ( μ2 ) =ln ( a ) +
Te= (GGL)+tsup T 2+ 460
H w , s=5,615 ρW ( RW ( H W −H R ) + R S LV ) =BTU /BN
(W= 62,4 lb/pie3 a condiciones normales)
4

Marbellsy Nazaretth Sanchez Salazar μo=5,664 x 10 −9 ( 1,60929 x 10 )

e T +460 =cp 14. Calculo de 
Explotación de Crudos Pesados H f =q oh ( H o , g + H w , s ) =BTU / D
8. Tasa de flujo de petróleo antes de la estimulación
2 π (1,127) hk ( P ¿ ¿ e−Pw ) 1 Hf
q oc = =B N D ¿ δ= 2
ℜ 2 π r h h́ M s (Ts−Tr )
μoc ln
rw 15. Temperatura promedio del yacimiento calculada
Inyección cíclica de vapor
Marx y Langenheim T avg ¿ c al¿=T r + ( T s−T r ) ( v r v z (1−δ )−δ )=° F
1. Calculo de las propiedades del vapor
9. Tasa de petróleo estimulada Comparar Tavg Supuesta con el Tavg calculada si difieren de
Ts=115,1 ( Ps 0,225 ) =° F 2 π (1,127) hk ( P ¿ ¿ e−P w ) una aproximación de 0,1 repetir desde el paso 6 usando el
Lv=1.318( Ps−0,08774 )=BTU /lb° F q oh= =B N D ¿ calculado como el nuevo supuesto (Iterar)
ℜ rh
2. Calculo de la tasa de inyección de calor. μoc ln + μoh ln Iter. Tavg)sup Tavg)cal qoh oh Ho,g Hw,s 
rh rw
Qi= ( 350
24 ) ist ( Cw (Ts−Tr ) + Xst Lv )=
BTU
h
10.Calculo de Vr y Vz
hi t D
Inyección continúa de vapor
Mand y Volek
3. Calculo del tiempo adimensional. (t en horas)
h́= =pies
F1 1. Calculo de la razón B y el factor F2c
4 K hob ∙ M ob ∙ t X st Lv X st Lv
t D= K hob t 4 K hob t
M h 2
s
2
t
x=log
( M ob r h 2
)
; y=log (
M ob h́2 ) B= =
H w −H r Cw ( T s−T r )
4. Calculo de F1 1
tD v r=0,180304−0,41269 x +0,18217 x 2+0,149516 x 3 +0,024183
F 2c = x 4
1+ B
F 1= v z=0,474884−0,56832 y −0,239719 y 2 −0,035737 y 3 tDc se obtiene por la tabla 5.1
1+0,85 √ t D
11. Calculo de Ho,g 2. Calculo del tiempo critico tc
5. Área de la zona de vapor calentada y Radio calentado
H o , g=( 5,615 ρ o C o + R g C g ) ( T avg−T R )=BTU /BN t M s 2 h t 2t Dc
Qi M s h t F 1 c=¿ =hrs ¿
A s= = pie 2 (o=60 lb/pie3 Co=0,5BTU/lbF a condiciones estándar)
4 K hob M ob
4 K hob M ob(Ts−Tr ) 3. Calculo del volumen de la zona e vapor
r 12. Cálculo de Rs (BN/BN) Si no se conoce Ps (tiny  tc) Determinar tD y F1 para tiempo de producción
A
h=¿
√ π
s
= pie ¿
Ps = (Ts
)
1
0,225
=lpca V s=
Qi M s h t2 F 1
=
pie 3
=acre− pie
Para calcular el Vs (volumen de la zona de vapor) ht2 115,1 4 K h ob M ob(Ts−Tr ) 43.560
6. Temperatura promedio del yacimiento supuesta Si Ps>Pw toda el agua es producida como vapor Rs=Rw (tiny > tc) Vs en función de F3 que se obtiene de la Fig 5.5
Ts+Tr Si Ps<Pw Rs se calcula (Rs no debe exceder Rw) 4. Calculo de la recuperación acumulada de petróleo
T ¿avg ¿= =° F
2 ∅ V s ( S oi−S orst ) hn
7. Calculo oh a Tavg y oc a Tr R s=
1 ρ g R g Ps
5,615 Pw Pw −Ps
−4
≈ 1,6 x 10 Rr
Ps
Pw −Ps ( ) N p=
5,615
Ec
ht
=BN( )
Constantes a y b – Ec de Andrade
 5. Calculo de la tasa de producción de petróleo
24 Q i ∅(S oi −S orst ) hn
q o=
5,615 M s (Ts−Tr )
Ec ( )
F =B N D
ht 2
6. Calculo de la razón petróleo/vapor acumulada
Np
OSR=
i st t
Razón petróleo/vapor instantánea
qot
OSR(t )=
i st
Myhill y Stegemeier
Razón petróleo/vapor
∅ (S oi −Sorst ) ( ρ w C w )
F os= ( 1+B ) F 4
Ms
Tasa de utilización de calor
Q s =ht M s ( Ts−Tr )=BTU /hr pie 2
Tasa de pérdidas de calor
K h ob (Ts−Tr ) BTU
Qc = =
√ π αob t hr pie 2
Calor ac umulado perdido Q ob=1−Es
F
Eficiencia E s= 1
tD
Marbellsy Nazaretth Sánchez Salazar
Explotación de Crudos Pesados