UNIVERSIDAD NACIONAL AUTNOMA DE MXICO.

FACULTAD DE QUMICA.

PROYECTO INGENIERA DE CALOR

ALUMNOS:
Garca Ochoa Jess Christian
Hernndez Salazar Emilio
Isidro Guadalupe Ana Lucero
Mendoza Meja Diana Itzel
Pea Reyes Vctor Leonel
Ruiz Prez Brenda

Ciudad Universitaria, Mxico, DF

CARTA DE TRANSMISIN

Proyecto Ingeniera de calor

Mxico Distrito Federal, 13 de Mayo del 2014
Para el Ingeniero Qumico Carlos lvarez Maciel.
A travs de la presente, nuestra compaa SOLIQ (Soluciones en Ingeniera Qumica)
nos permitimos presentar el proyecto Tren de precalentamiento para una torre de
fraccionamiento de petrleo crudo despuntado.
Para que un flujo de 105 000 BPD de petrleo crudo despuntado a 90F salga de un
tren de precalentamiento a 500F, y que despus entre al calentador a fuego directo
(horno) previo a la torre de fraccionamiento de cual se requiere que el flujo de
petrleo crudo despuntado salga a 660F, se hacen las siguientes recomendaciones:
El tren de precalentamiento estar compuesto de 17 intercambiadores de haz y
envolvente; se emplearan las corrientes provenientes del fraccionamiento primario del
crudo para calentar el petrleo crudo; la topologa de las corrientes para todos los
intercambiadores ser que el petrleo estar en tubos y los fraccionamientos
primarios del crudo (nafa, turbosina, querosina, disel, gasleo pesado y residuo) por
la coraza.
Para este proceso de precalentamiento se consider el uso de una bomba para
impulsar al petrleo crudo despuntado hacia el tren de precalentamiento, y un
calentador a fuego directo (horno) para aumentar finalmente la temperatura del
petrleo crudo a 660F, que despus entrara a la torre de fraccionamiento.
El Costo estimado para la inversin inicial en equipos de Intercambiadores de Calor y
el Horno, Agua y Vapor de Proceso es de: 9, 825, 760 US$
Sin ms por el momento nos despedimos desendole xito en todas sus labores.

Cordialmente SOLIQ

_____________
Isidro Guadalupe Ana Lucero

_____________

____________

Hernndez Salazar Emilio

Garca Ochoa Jess Christian

_____________
______________
Mendoza Meja Diana Itzel

_____________
Pea Reyes Vctor Leonel

PLANTEAMIENTO DE LA RESOLUCIN DEL PROBLEMA

Pgina 2

Ruiz Prez Brenda

Proyecto Ingeniera de calor

El petrleo es un recurso de alto valor industrial y econmico, pues de l se

extraen agentes qumicos de gran utilidad, adems de diversos tipos de
combustibles: gasolina, disel, combustleo, entre otros. La separacin del
petrleo crudo en sus fracciones principales se lleva a cabo en una primer torre
de destilacin (fraccionamiento primario), la cual consume una gran cantidad
de energa. Para aprovechar al mximo los recursos energticos, se aprovecha
la energa trmica de las fracciones de la destilacin primaria del petrleo para
calentar al crudo previo a su entrada a la torre de destilacin. Esto se efecta
mediante un tren de precalentamiento consistente en una serie de
intercambiadores de calor, como se muestra en la figura siguiente:

A las diferentes fracciones de petrleo empleadas en los intercambiadores de

calor del tren de precalentamiento, se les somete a un refinamiento posterior
para transformarlas en productos de un alto valor agregado.
PROBLEMA:
Una corriente de petrleo crudo de 35API y con un flujo de 105,000 BPD,
entra a una serie de intercambiadores de calor para incrementar su
temperatura de 90F hasta 500F, empleando las corrientes provenientes del
fraccionamiento primario del crudo. Si los intercambiadores de calor son de haz
y envolvente, disea el tren de precalentamiento.

Pgina 3

Proyecto Ingeniera de calor

Se trabaj con las siguientes propiedades, las cuales fueron calculadas en

base a la temperatura promedio (tm o Tm) de entrada y salida para cada
sustancia:

Topologa de las corrientes:

Se plantea que por dentro del tubo vaya el petrleo debido a su alta viscosidad
Ejemplo Querosina
Algoritmo de clculo:
vq=4700

Barriles
ft
=0.31
Da
s

W =vqdensidad querosina
W =0.31

f t3
lb
lb
54.3123 3 =16.8368
s
s
ft

Petrleo
vp=105000 BPD=6.83

f t3
s

w=vpdensidad petroleo
lb
3600 s
ft
lb
s
w=6.83
59.3065 3 =405.06
=1458228.222 lb /h
s
1h
ft
3

Pgina 4

Proyecto Ingeniera de calor

Primero se plantea calcular la carga trmica en base a la temperatura de salida

del intercambiador anterior, que es t 1 , o en el caso de la nafta, la
temperatura del petrleo es 90F, t 2 es una temperatura supuesta a la que
deseamos salga el petrleo.
En base a esto se plantea lo siguiente, calculando la carga trmica:
q=wC p petrleo(t 2t 1)
Con lo que se obtiene el calor que se transferira.
Despejando la temperatura de salida de la querosina

T1

q
+T
wC p destilado 2

T1=

Calculando el coeficiente de transferencia de calor del lado de los tubos

i 0.14 di
( )
w
d0
d i vi i 0.8 Cp i i 13

i
ki
hio di
=0.027
ki

( )

Primero debemos conocer los dimetros, se propone trabajar con un dimetro

calculado mediante la ecuacin de Barlow
d 0=1, di=0.87
Se puede conocer el nmero de Reynolds con las propiedades del petrleo
=

di v i i d iG
=
i
i

Pr=

Cp
k

0.8

1
3

( Pr ) (

ki
)
d0

hio
=0.027
t
Con lo que podemos calcular

hio BTU
(
)
t h f t 2 F

Podemoscalcular el a ft
Pgina 5

Proyecto Ingeniera de calor

a ft =

Wt
(f t 2)
Gt

GT =flujo volumtrico petrleo (

a ft =

lb
)
2
f t s

W
(f t 2)
GT

2 2
a f /tubo = d i (f t )
4
W
(f t 2 )
a
G
N tubos= ft = T
=Nmero de tubos
a f
2 2
di (f t )
tubo
4
Pt=1.25 d 0
Pt=c d 0 c =Ptd 0 ()
Coeficiente de transferencia del lado de la coraza:
0 0.14

w
d 0 v 0 0 0.8 Cp0 0 13

0
k0
ho de
=0.027
k0

d e =Dimetro equivalente , considerandoun arreglo tringular

2

1
1 d
4 ( P t ( 0.86 ) )
2
2 4
de=
()
1
d
2
Con las propiedades del petrleo a las temperaturas medias se puede obtener
0.8

1
3

( Pr ) (

ks
)
de

ho
=0.027
s
Una vez que se calcula ese trmino se puede sustituir en la siguiente ecuacin
considerando nuevamente que el petrleo va por dentro y el fluido caliente va
por fuera:
Pgina 6

Proyecto Ingeniera de calor

t w =t m+

1
hio
t

()
( )( )
1
1
+
ho
hio
t
s

(T m t m)

En caso de que el fluido caliente vaya por fuera, que es lo que ocurre en el
proyecto:

T w =T m

1
ho
t

()
( )( )
1
1
+
ho
hio
t
s

(T mt m )

Una vez realizado esto podemos obtener las viscosidades acorde a la

temperatura de pared y calcular t y s
t 0.14

wt
t =
s 0.14

ws
s=

Posteriormente procedemos a calcular el coeficiente de transferencia de calor

pero antes de eso necesitamos obtener hio y h o
hio =

hio
BTU
t (
)
t
h f t2 F

( )

ho
BTU
s (
)
2
s
hft F

( )

ho =

Una vez calculado esto podemos obtener el coeficiente global de transferencia

de calor limpio:
U c=

BTU
)
1 1 h f t 2 F
+
hio ho
(

Pgina 7

Proyecto Ingeniera de calor

Se puede calcular el coeficiente de calor sucio considerando el

ambos datos del Kern.
U D=

1
1
+R +R
U c di do

Rdi y Rdo

BTU
)
h f t2 F

Una vez obtenido este valor se puede proceder a calcular el rea de

transferencia de calor:
q=U DA TLMTDF t (
Ft

BTU
)
h f t 2 F

es un valor que se puede obtener del Kern, es un factor de correccin

Posteriormente podemos calcular el rea de transferencia de calor y el LMTD

con las temperaturas establecidas previamente.
AT =

q
(f t 2)
U DLMTD

LMTD=

( T 1t 2) (T 2t 1)
( F )
( T 1t 2 )
ln (

(T 2t 1 )

Posteriormente se procede a calcular el nmero total de tubos:

N TotalTubos =

AT
A T (f t 2 )
=
A tubos d 0L(f t 2 )

Con este valor se puede calcular el nmero de pasos, recordando que se tiene
que redondear a un valor entero.
Con el nmero total de tubos se puede obtener el nmero de pasos
N pasos t =

N TotalTubos
N Tubos / paso

N pasos s=

N pasos t
2

Una vez obtenida la configuracin del intercambiador podemos obtener la

correccin del rea de transferencia de calor A T
R=

T 1T 2
t 1t 2
Pgina 8

Proyecto Ingeniera de calor

S=

t 2t 1
T 1t 1

Estos valores se pueden correlacionar en las grficas de las pginas 933-938

del Kern, de acuerdo a la configuracin que se tenga.
Obtenemos con ello un valor de correccin Ft , el cul podremos emplear
para obtener el rea de transferencia de calor con una modificacin, donde
A T es el rea de transferencia de calor previamente obtenido, el cual es
equivalente a tener un Ft =1 .
A Tc=

AT
( f t2 )
Ft

Podemos obtener de nueva cuenta el nmero total de tubos con esta rea
corregida.
N TotalTubos =

A Tc
d 0L

Una vez obtenido esto se puede proceder a calcular el dimetro de la coraza,

para esto se puede considerar un arreglo triangular
Ds=1.15Pt N TotalTubos ()
Arreglo de dimetro equivalente para el clculo del nmero de Reynolds
4 ( Pt 2
Arreglo cuadrado de=

d2
)
4

4 Pt 2 d 2
d

3.44 P t 2 d 2
Arreglo tringular de=
d
Cadas de presin
Por el lado de los tubos
Pt = P tubos+ P retornos
Para calcular la cada de presin en los tubos utilizamos:
P tubos=

2 f ' N +1
5.22 X 1010 di Sg t

Para poder obtener el factor de friccin en los tubos, calculamos el nmero de

Reynolds en los mismos.

Pgina 9

Proyecto Ingeniera de calor

Ret=

diGt
t
N +1=12

Donde

L
B

Para el clculo de la cada de presin en los retornos

2

P retornos=

( )

4n
( psi)
Sg 2 gc

Del lado de la coraza utilizamos:

2

Ps=

f ( Gs ) Ds( N +1)
10

5.22 X 10 deSgs

Tambin para poder obtener el factor de friccin , calculamos el nmero de

Reynolds
Res=

deGs
s

Para realizar la correccin del ho se realiza la siguiente correccin mediante el

siguiente mtodo.
Mtodo de Bell
hcoraza =hidealJ c J lJ bJ sJ r
Dnde:
J c =Factor de correccn por laconfiguracin de lamampara
J l=Factor de correcin para losefectos de fugasla mampara
J b=Factor de correcin para los efectos de derrames en el haz de tubos
J s=Factor de correcin para el espaciamiento desigual
de las mamparas enla entrada y /o slida
J r =Factor de correccin para la acumulacin de gradientes adversos
con nmeros de Reynolds menores a 100

Para poder calcular estos factores se plantea calcular ciertos parmetros que
intervienen, los cuales se muestran a continuacin.

Pgina 10

Proyecto Ingeniera de calor

Espaciamiento de los tubos paralelos al flujo

pp

Para arreglos triangulares:

p p= psen

( 3 )

Para arreglos cuadrados

p p= p
Para espaciamiento de los tubos en direccin normal al flujo
Para arreglos triangulares:

pn= pcos( )
3
Para arreglos en cuadrados:
pn= p
Nmero de hileras de tubos cruzados en una seccin de flujo transversal
D s [12(
N c=

lc
)]
Ds

pp

Dnde:
l c = Alturade lamampara
lc
=Corte del deflector
Ds
Podemos calcular la fraccin de tubos
Ds=Dametrod e lacoraza

2D s2 l c
1
Ds2lc
Ds2 lc
Fc = [ +
sen cos1
2 cos1 (
)]

D otl
Dotl
Dotl

Dotl =

Ds1
=Dmetrolmite exterior de tubos de la carcaza
12

Tambin podemos calcular el nmero de hileras eficaces de flujo transversal en

cada ventana N cw
N cw=

0.8 l c
pp
Pgina 11

Proyecto Ingeniera de calor

El rea de flujo transversal en la lnea central o cerca de ella para una seccin
de flujo transversal S m se puede obtener de la siguiente manera:
Para configuraciones cuadradas:
S m=B [D sD otl +

Dotl d 0
( pd 0 ) ]
pn

Dnde:
d 0=Dimetro exterior de los tubos
Para configuraciones triangulares:
S m=B [D sD otl +

Dotl d 0
( pd 0 ) ]
p

Otro parmetro es la fraccin de rea de flujo transversal disponible para el

flujo de desviacin Fbp
Fbp =

( d iDotl )B
Sm

rea de fuga entre tubos y mamparas:

2

S tb =bd 0 N T(1+ Fc)

Dnde b es un factor
b=1.701104
Se pueda calcular el rea de fuga entre la carcasa y la mampara:
S sb =

D s sb
2l c
1
[cos 1
]
2
Ds

sb =Espaciolibre de flanqueo diametral entre lacarcasa y el deflector

En la siguiente tabla se muestran valores tpicos de

sb

para tubo comercial

sb

Di (ft)

(ft)

0.67 - 1.16

0.0083

1.17 - 1.49

0.0104

1.5 - 2

0.0125

Pgina 12

Proyecto Ingeniera de calor

Clculos para obtener el rea de flujo a travs de la ventana:

cos1 1

) (

2l c
2l c
2l c
1
1 1
Ds
Ds
Ds

)(

S wg=

D s2

Dnde:
S wt =

NT
( 1F c )d 02
8

Sumando ambos trminos obtenemos:

S w =S wg + Swt
Retomando:
S wg = reatotal de laventana
S wt = reaocupada por lostubos en la ventana
S w = rea de flujo a travs de la ventana

Podemos calcular el ngulo de corte de la mamparas ( b en radiantes

mediante la siguiente ecuacin:
b=2cos1 (1

2l c
)
Ds

Calculando el nmero de mamparas

N b=

12L
1
B

Una vez calculados estos parmetros se puede proceder a obtener los factores
de correccin mediante grficos proporcionados.
Correccin de las cadas de presion
2

f w N
Pbk=b . k 2 . m

Sm

0.14

( )

b=1.5432e-07 para sistema internacional (adim.)

Pgina 13

Proyecto Ingeniera de calor

b=3.3282e-11 para sistema ingls (adim)

Cada de presin para una seccin una ventana ideal ( P k w )
si 100
W 2 (2+0.6 N cw )
PWk =b .
Sw . Sm.
b=1.5432E-07 para sistema internacional (adim)
b=3.3252E-11 para sistema ingls (adim)

posteriormente buscamos en graficas los factores de correccin

buscamos el factor de correccin por el efecto de la fuga en los deflectores (Ri)
tambin el, factor de correccin para la desviacin del haz de tubos (Rb)
calculamos la cada de presin del lado de la carcasa

( Nb1 ) . Pbk
Ncw
MEMORIA DE CLCULO
1+
nc
Ps=[ . Rb+ Nb. PmK ] Rs+ 2 Pbk . Rb .
En las siguientes tablas se muestran los resultados para cada arreglo con el
algoritmo antes especificado.
Arregl 1 (nafta-petrleo)

Gt (lb/h*ft^2)
Ret
Prt
hio/t
(btu/h*ft^2*F)
aft por tubo (ft^2)
numero de tubos
por paso

Nmero total de
tubos
Nmero de pasos
lado tubo
nmero de pasos
shell

500000
3242.073
75
97.57524
91
49.92520
87
0.008659
01

R
S

149

Ft

pt (in)

1.725

Ds (in)

C (in)

0.225

Ds (ft)
Pgina 14

1132
8
4
4.440284
84
0.214973
92
0.98
72.56263
48
6.046886
24

METODO BELL

Proyecto Ingeniera de calor

lc
1.8141
Swg
29.9409
Ssb+stb/sm
0.0022
pp
0.1245
Sw
33.4008
ssb/ssb+stb
0.8870
Gspn
(lb/h*ft^2)
100000
afs
(ft^2)
0.0719
Nb
13.1852
jl
0.9800 1.196376
68
Nc
19.4291
p wk
0.0603
nc+ncw
31.0866
Res
8583.768
B (ft)
1.516850
Dotl
5.9636
Fk43
0.1500
jr*
0.4400
3
Fc Prs0.5019
p bk
0.2223
jr
3.999606
LADO
Vt (ft/s)1.0000 2.649540
Ncw 11.6575
Nss
Jc
0.8900
65
TUBOS
04
Smho/s
17.6492
Rb
1.0000
jb
0.9000
99.80634

0.00038
(btu/h*ft^2*F)
Fbp
0.5035
Ri85
0.9000
JS
1.0200
Tw (F)
132.2700
0.84
stb
0.0044
ps
3.8727
multiplicacionSg 0.8007
41
ssb
0.0348
Sb
0.0058
t
1.108994
pt (psi)
2.976973
Swg 29.9409
71
88
a s 0.4000
0.985357
p retorno(psi)
0.032
Swt
3.4599
78
55.36679
3.200973
hio (btu/h*ft^2*F)
24
p total tub(psi)
88
78.74280
LADO
ho (btu/h*ft^2*F)
5
CORAZA
Sg
0.65
32.50875
Uc (btu/h*ft^2*F)
87

0.0022
28.76792
0.529908
Ud (btu/h*ft^2*F)
11
Vs (ft/s)
01
28.55484
3.977968
LMTD
42
ps tot (psi)
39

Arreglo 2 (turbosina-petrleo)

Gt (lb/h*ft^2)

240000

Ret

2606.627297

Prt
hio/t (btu/h*ft^2*F)

61.52150581
35.49359148

Nmero total de
tubos
Nmero de pasos
lado tubo
nmero de pasos
shell
R

aft por tubo (ft^2)

numero de tubos por
paso
pt (in)

0.008659015

2
0
0.170367
57

303
1.725

Ft
Ds (in)

0.98
69.67671

Pgina 15

1234
4

Proyecto Ingeniera de calor

C (in)

0.225

Ds (ft)

Gs (lb/h*ft^2)
Res

40.0106398
10076.59772

afs (ft^2)
B (ft)

Prs
ho/s (btu/h*ft^2*F)
Tw (F)

3.8533336
69.67671127
213.5762174

t
s

1.108994712
0.985357783

hio

(btu/h*ft^2*F)

LADO
TUBOS

pt (psi)
p retorno(psi)

55.36679243

ho (btu/h*ft^2*F)
Uc (btu/h*ft^2*F)

78.74280505
32.50875869

Ud (btu/h*ft^2*F)
LMTD

28.76792107
28.55484417

Vt (ft/s)

Sg

p total tub(psi)
LADO
CORAZA

Sg

Vs (ft/s)
ps tot (psi)

13
5.806392
61
0.883433
22
0
0.225
0.00033
0.82
0.300142
3
0.0085
0.325642
3
0.63
0.0022
0.542852
8
2.521822

metodo de bell
lc
pp
pn
Nc
Dotl
Fc
Ncw
Sm
Fbp
stb
ssb
Swg
a
Swt

1.7419
0.1245
0.0719
18.6564
5.7231
0.5022
11.1938
13.0113
0.5037
0.0048
0.0334
27.6067
0.4000
3.7695

Swg
Sw
Nb
p wk
Fk
p bk
Nss
Rb
R1
ps
Sb

27.6067
31.3762
18.1457
3.2674
0.0200
34.7040
1.0000
0.7500
3.7310
0.0058

Ssb+stb/sm
ssb/ssb+stb
jl
nc+ncw
jr*
jr
Jc
jb
JS
multiplicacion

0.0029
0.8737
0.9800
29.8502
0.4900
1.0000
0.8800
0.8400
1.0200
0.7389

Arreglo 3 (querosina-petroleo)
Gt (lb/h*ft^2)

400000

Ret
Prt

4965.004374
55.65267795
Pgina 16

Nmero total de
tubos
Nmero de pasos
lado tubo
nmero de pasos

587
2
1

Proyecto Ingeniera de calor

shell
hio/t (btu/h*ft^2*F)

56.73247636

aft por tubo (ft^2)

numero de tubos por
paso

0.008659015

359

Ft

pt (in)

1.725

Ds (in)

C (in)

0.225

Ds (ft)

Gs (lb/h*ft^2)

64.49756844

afs (ft^2)

Res

6621.764216

B (ft)

Prs
ho/s (btu/h*ft^2*F)
Tw (F)

6.585332069
48.02453155
247.8960785

t
s

1.136872081
0.989926

pt (psi)
p retorno(psi)

(btu/h*ft^2*F)

64.49756844

p total tub(psi)

ho (btu/h*ft^2*F)
Uc (btu/h*ft^2*F)

61.47316551
31.47453043

Ud (btu/h*ft^2*F)

27.87711244

Vs (ft/s)

LMTD

67.35212921

ps tot (psi)

hio

LAD0
TUBOS

Vt (ft/s)

Sg

LADO
CORAZA

Sg

metodo de bell
lc
pp
pn
Nc
Dotl
Fc
Ncw
Sm
Fbp
stb
ssb

1.2006 Swg
0.1245
Sw
0.0719
Nb
12.8589 p wk
3.9187
Fk
0.5053 p bk
7.7153
Nss
6.7484
Rb
0.5054
R1
0.0023 ps
0.0230 Sb

13.1149
14.8968
23.3633
6.9678
0.0200
26.2434
1.0000
0.9000
3.6053
0.0058

Ssb+stb/sm
ssb/ssb+stb
jl
nc+ncw
jr*
jr
Jc
jb
JS
multiplicacion

Pgina 17

0.0037
0.9092
0.9900
20.5742
0.5200
1.0000
0.9000
0.8400
1.0200
0.7634

21.63629
56
0.045542
09
0.98
48.02453
16
4.002044
3
0.466842
45
0.894324
29
2.198304
67
0.00033
0.81
0.418444
14
0.02
0.438444
14
0.68
0.0021
0.549576
17
1.226963
44

Proyecto Ingeniera de calor

Swg
a
Swt

13.1149
0.4000
1.7819

Arreglo 4 (diesel-petroleo)

Gt (lb/h*ft^2)
Ret
Prt
hio/t
(btu/h*ft^2*F)
aft por tubo (ft^2)
numero de tubos
por paso

Nmero total de
tubos
Nmero de pasos
lado tubo
nmero de pasos
shell

600000
7571.631
67
33.70962
21
100.9163
37
0.002922
47

R
S

140

Ft

pt (in)

1.15

Ds (in)

C (in)

0.15

Ds (ft)

Gs (lb/h*ft^2)

90000
4213.849
Res
96
8.482225
Prs
3
ho/s
73.74604
(btu/h*ft^2*F)
38
257.5484
Tw (F)
45
1.087299
t
59
s
0.989926
109.7262
hio (btu/h*ft^2*F)
92
54.52807
ho (btu/h*ft^2*F)
9
36.42620
Uc (btu/h*ft^2*F)
82
31.69294
Ud (btu/h*ft^2*F)
73
45.34758
LMTD
48

0.00026

Sg

p total tub(psi)

0.69
6.268984
26
0.02
6.408984
26

Sg

0.8

0.0022
0.500801
28
4.083749
47

pt (psi)
p retorno(psi)

Vs (ft/s)
ps tot (psi)

Pgina 18

4
5.205465
04
0.191671
44

Vt (ft/s)

B (ft)

LADO
CORAZA

0.98
44.22183
24
3.685152
7
0.414373
85
0.862071
79
3.338675
21

afs (ft^2)

LADO
TUBOS

1119

Proyecto Ingeniera de calor

metodo bell
lc
pp
pn
Nc
Dotl
Fc
Ncw
Sm
Fbp
stb
ssb
Swg
a
Swt

1.1055
0.0830
0.0479
17.7610
3.6018
0.5061
10.6566
6.0664
0.5032
0.0019
0.0212
11.1202
0.4000
1.5071

Swg
Sw
Nb
p wk
Fk
p bk
Nss
Rb
R1
ps
Sb

11.1202
12.6272
24.1999
8.8000
0.0200
35.6461
1.0000
0.9000
4.8461
0.0058

Ssb+stb/sm
ssb/ssb+stb
jl
nc+ncw
jr*
jr
Jc
jb
JS
multiplicacion

0.0038
0.9157
0.9500
28.4177
0.6500
1.0000
0.9000
0.8400
1.0400
0.7469

ARREGLO 5 (gasleo pesado-petrleo)

Gt (lb/h*ft^2)
Ret
Prt
hio/t
(btu/h*ft^2*F)
aft por tubo (ft^2)
numero de tubos
por paso

Nmero total de
tubos
Nmero de pasos
lado tubo
nmero de pasos
shell

800000
24825.02
19
25.26483
17
155.9309
04
0.008659
01

R
S

162

Ft

pt (in)

1.725

Ds (in)

C (in)

0.225

Ds (ft)

Gs (lb/h*ft^2)

100000
17827.82
67
2.987467
03
106.5867
77
307.9760
29
1.072356
71

afs (ft^2)

Res
Prs
ho/s
(btu/h*ft^2*F)
Tw (F)
t

3
11.49651
35
0.086357
68

Vt (ft/s)

0.00021

Sg

0.78
3.517927
09

pt (psi)

Pgina 19

0.98
58.79232
46
4.899360
38
0.940309
89
1.471425
24
4.878643
74

B (ft)
LADO
TUBOS

879

Proyecto Ingeniera de calor

0.989926
167.2135
hio (btu/h*ft^2*F)
52
75.30616
ho (btu/h*ft^2*F)
29
51.92242
Uc (btu/h*ft^2*F)
2
42.80913
Ud (btu/h*ft^2*F)
15
61.24386
LMTD
66

p retorno(psi)
p total tub(psi)

0.085
3.942927
09

Sg

0.88

0.0018
0.609830
47
1.454125
56

LADO
CORAZA

Vs (ft/s)
ps tot (psi)

metodo de bell
lc
pp
pn
Nc
Dotl
Fc
Ncw
Sm
Fbp
stb
ssb
Swg
a
Swt

1.4698 Swg
0.1245
Sw
0.0719
Nb
15.7420 p wk
4.8160
Fk
0.5034 p bk
9.4452
Nss
13.7437
Rb
0.5044
R1
0.0034 ps
0.0282 Sb
19.6553
0.4000
2.6782

19.6553
22.3336
14.5923
4.5759
0.0200
14.2371
1.0000
0.9000
1.8213
0.0058

Ssb+stb/sm
ssb/ssb+stb
jl
nc+ncw
jr*
jr
Jc
jb
JS
multiplicacion

0.0023
0.8912
0.9800
25.1873
0.5200
1.0000
0.8500
0.8400
1.0200
0.7137

Arreglo 6 (residuos-petrleo)

Gt (lb/h*ft^2)

600000

Ret

16108.68086

Prt

15.00154945

Nmero total de
tubos
Nmero de pasos
lado tubo
nmero de pasos
shell

hio/t (btu/h*ft^2*F)

179.2730916

aft por tubo (ft^2)

numero de tubos por
paso

0.001860172

170

Ft

pt (in)

0.8625

Ds (in)
Pgina 20

1358
8
4
2.999351
49
0.332500
24
0.98
39.72119
4

Proyecto Ingeniera de calor

C (in)

0.1125

Ds (ft)

Gs (lb/h*ft^2)

172.1962482

afs (ft^2)

Res

18540.9398

B (ft)

Prs
ho/s (btu/h*ft^2*F)
Tw (F)

2.500258241
39.72119403
397.2992019

t
s

0.960524788
0.989926

hio

(btu/h*ft^2*F)

LADO
TUBOS

Vt (ft/s)

Sg
pt (psi)
p retorno(psi)

172.1962482

p total tub(psi)
LADO
CORAZA

ho (btu/h*ft^2*F)
Uc (btu/h*ft^2*F)

147.4470804
79.4317659

Sg

Ud (btu/h*ft^2*F)

59.91819345

Vs (ft/s)

LMTD

48.19535274

ps tot (psi)

metodo de bell
lc
pp
pn
Nc
Dotl
Fc
Ncw
Sm
Fbp
stb
ssb
Swg
a
Swt

0.9930 Swg
0.0719 Sw
0.0719
Nb
18.4214 p wk
3.2268
Fk
0.5074 p bk
11.0529 Nss
2.9127
Rb
0.9978
R1
0.0013 ps
0.0190 Sb
8.9719
0.4000
1.0262

8.9719
9.9981
22.8707
2.2732
0.0200
116.8046
1.0000
0.9000
11.8432
0.0058

Ssb+stb/sm
ssb/ssb+stb
jl
nc+ncw
jr*
jr
Jc
jb
JS
multiplicacion

Horno
Carga trmica:
q=W petrleoCp m petroleo(T 2T 1)
q=1196376.678

lb
BTU
0.63
( 660393 ) F
h
lb F
Pgina 21

0.0070
0.9346
0.9800
29.4743
0.4900
1.0000
0.9000
0.8400
1.0200
0.7557

3.310099
5
0.394822
49
0.914465
7
3.603603
6
0.00021
0.72
6.884915
28
0.0067
6.931815
28
0.865
0.0018
0.960960
96
10.72194
26

Proyecto Ingeniera de calor

q=201242521.1

BTU
h

Se tiene que considerar un poder calorfico inferior de 6,324 MJ/bbl

J
1 bbl
bbl
BTU
3
BTU
6 5.615 f t
6,32410
=1067502.889
1055.05 J
f t3
BTU
f t 3
3
ft
BTU
1067502.889
=17099.75581
62.427961lb
lb
Poder calorifico=17099.75581

BTU
lb

Para obtener el flujo de combustleo:

W combustleo =

q
lb
( )
PClT h

BTU
h
lb
W combustleo =
( )
BTU
h
17099.75581
0.87
lb
201242521.1

W combustleo =13527.28338

W
ncombustleo = combustleo =
PM

lb
h

13527.28338
226

lb
h

lb
lbmol

ncombustleo =59.85523621 lbmol

Combustleo y aire:
Combustleo 25 C
Combustin
Aire +

O2 +

Reaccin qumica:
Pgina 22

C O2 +

O2 +
H2O + N2

Proyecto Ingeniera de calor

49
O =16 CO 2+ 17 H 2 O
2 2

C16 H 34+

Balance de

O2 (10% de exceso)

n Alimentados nreaccinan=n salida

1613.098616lbmol 1466.453287lbmol =nsalida
n salida=146.6453287 lbmol
Flujo de aire a la entrada
W O =1613.098616lbmol32
2

lb
lbmol

W aire =46926.50519 O2 +176533.0433 N 2

Balance de

CO2

n formadosCO =nsalida CO =957.6837793 lbmol

2

Balance de

H2O

n formados H O =nsalida H O =1017.539016lbmol

2

Balance de

N2

nentrada =nsalida =6068.323364 lbmol

Flujo msico de gases de combustin
W gases por hora=957.6837 lbmol44
W gases=235140.8088

lb
h

H f Reactivos = H f O + H f N + H f C
2

Hf C

16

H 34

lb
lb
CO2 ) +146.6453 lbmol32
(
( O ) +1017.539 lbmol18 lb
lbmol
lbmol 2

16

H 34

( para N 2 y O2 es cero)

=11690556.23 BTU

H f Reactivos = H f C

16

H 34

H f productos= H f C O + H f H
2

H f C O =162023955.6 BTU
2

Pgina 23

Proyecto Ingeniera de calor

H f H O=125044433.4 BTU
2

H 1= H f Productos H f Reactivos
H 1=275377832.8 BTU
Tad

H 2=

Productos

298.15

i=1

ni Cpi dT

q= H 1 + H 2
Resolviendo el balance de energa
Tad

Q H f Productos + H f Reactivos =+

298.15

Despejando

Productos

i=1

n i Cp i dT

T ad :

q= H 1 + H 2
q=286803703.4+134933.8725T ad 40230534.09
T ad=

201242521.1+286803703.4 +40230534.09
134933.8725

T ad=3976.96 F=2191,64 C
Transferencia de calor hacia el crudo y clculo del rea de transferencia de
calor.
Para este clculo consideraremos que ah transferencia de calor por
conveccin y radiacin.

rea de transferencia de calor:

q
=( h conv +h rad )(T G T s)
AT
Dnde:
T G =3976.96 F=4436.63 R
T s=2800 F=3259.67 R
Clculo del coeficiente de transferencia de calor por radiacin:

Pgina 24

Proyecto Ingeniera de calor

TG 4 T s
+
100
100

0.173 s
hrad =

( )( )

Obteniendo un coeficiente de transferencia por radiacin de:

hrad =20.177

BTU
h f t 2 F

Clculo del coeficiente de transferencia de calor por conveccin:

k 3 G3CpG BT
hconv =116 G
d 0 G

0.25

Dnde:
Cp G=capacidad calorifica de los gases
k G =conductividad trmica de los gases
G =densidad de los gases
B=

1
TG

T =T GT s
G=viscosidad de los gases
d 0=dimetro comercial propuesto(aceroinoxidable de 15Cd 40)
hconv =0.171739894
Se puede notar que el coeficiente de radiacin es mucho mayor comparado
con el de conveccin, siendo muy considerable para el diseo de un horno.
Retomando la ecuacin del flux de calor.
q
=( h conv +h rad )(T G T s)
AT
q
BTU
=23950.47987 2
AT
f t h
BTU
, sugeridos en el documento, esto puede
2
f t h
deberse al calor necesario para que el crudo alcance los 660 F y a la
Siendo mayor que 15000

Pgina 25

Proyecto Ingeniera de calor

composicin del combustleo sugerida, sin embargo se aproxima nuestro

diseo.
Clculo del rea de transferencia:
Para este clculo consideraremos el flux de calor de:
AT =

q
BTU
=12000 2
AT
f t h

q
12000

BTU
f t 2h

BTU
h
BTU
12000 2
f t h

201242521.1
AT =

A T =16770.21009 f t 2

Temperatura de los gases de combustin despus de haber cedido calor al

crudo:
T ag=Temperatura de los gases despes de haber cedido calor al crudo
q=W gasesCpm gases( T ad T ag)
T ag=T ad

q
W gasesCp m gases

BTU
h
T ag=3976.96 F
lb
BTU
235140.8088 0.4
h
lb F
201242521.1

T ag=1837.3642 F
Para obtener el flujo de agua con cambio de estado a vapor saturado, se
propone el siguiente balance trmico:
q w =w s vaporizacin + ws Cp l ( t saturacint 1 )
Incgnita q w
w s=

y ws
qw

vaporizaci n +Cp l ( t saturaci nt 1 )

Dnde:
Pgina 26

Proyecto Ingeniera de calor

BTU
lb

vaporizaci n =923.13

t saturaci n=282.27 F t 1=77 F

Cp lquido=1.0021

BTU
@ tm=180 F
lb F

Para calcular el calor emitido por los gases se puede hacer uso de la siguiente
ecuacin:
q w =Cp gases decombustinW gases(T agT sg )
Dnde:
T ag=Temperatura de los gases despes de ceder calor al crudo
q w =Cp gases decombustinW gases(T agT sg )
T sg =t saturaci n+ 10 C=332.27 F

q w = 0.42

BTU
lb
235140.8088 (1837.3642 F332.27 F)
lb F
h

q w =148641817.8
w agua=ws =

BTU
h

qw
vaporizaci n +Cpl ( t saturaci n t 1 )
148641817.8

w s=

923.13

BTU
h

BTU
BTU
+1.0021
( 282.27 F77 F )
lb
lb F
lb
vapor saturado
h

w agua=782.783384

Bomba
POTENCIA REQUERIDA EN LA BOMBA

W
WHP=

( lbs )W ( lbflbft ) =HP

f

550

Trabajo de flecha que debe desarrollar la bomba

Pgina 27

Proyecto Ingeniera de calor

12
1 ft

(426.6814.7)
W f=

lb

PD Ps
=

Flujo msico de crudo

W =1288 504

lb 1 h
lb
=358
h 3600 s
s

WHP=736 HP
Sabiendo

WHP
=MB
EHP

MB= MB =0.81
Obteniendo la potencia elctrica

EHP=WHP MB=596.2 HP
Watts consumidos

EHP=

Watts
745.7

Watts=EHP745.7=444 586.3

Watts
h

Costos
La estimacin de costos se hace con referencia a las grficas del libro Chemical Enginering Plant Cost
Index de Peters & Timmerhaus, tomado como parmetros importantes el rea de transferencia de calor,
dimetro y longitud de tubos y presin de operacin, tomando como material de tubos acero al carbn.
En primer instancia, se toma un precio de referencia con base en el rea de transferencia de calor para
cada intercambiador, Para el caso de la corriente de Nafta con una rea de 8 716.5 ft 2, se obtiene un
costo de 40, 000.00 dlares.

Pgina 28

Proyecto Ingeniera de calor

Despus se aplican a tal precio, factores que corrigen el costo segn el dimetro y longitud de tubos y
presin de operacin. Estos factores son tomados de grficas.
En el caso de la corriente de Nafta se tienen dos intercambiadores, por lo que el precio se duplica.
Por ltimo, se actualizan los precios del ao de 1990 al 2014, mediante el ndice CEPCI.

actual
( CEPCICEPCI
ao del mtodo )

Costo actual=Costo ao delmtodo

=192,610.24
( 637.3
368.3 )

Costo actual=$ 111,310.77

El mtodo anterior se realiza para todos los intercambiadores de las corrientes provenientes de la torre.
Los costos que se obtienen se muestran a continuacin:

Para el caso de horno y bomba, no se tomaron factores de ajuste, ms que la actualizacin de precio en
el caso de la bomba

Pgina 29

Proyecto Ingeniera de calor

Para los costos fijos se toman en cuenta todos los costos generados por la compra de los
intercambiadores, horno y bomba. Obteniendo as

Para el clculo de los costos variables se toma en cuenta el mantenimiento anual de los equipos, la
energa electrica consumida y costo por la generacin de vapor de agua, adems de la utilidad
correspondiente a la generacin de vapor (lo que representa un costo negativo, pues es ganancia).

Consideraciones:

El mantenimiento anual se toma como el 10% de la inversin inicial en equipos.

Electricidad

444. 588 kWh consumidos

1.22
( $kWh
)=$ 542.4=41.34 US $

Vapor de agua

Costo vapor agua=W combustible . gen .vapor Preciocombustible +W aguaPrecioagua

Costo vapor agua = $48 384 = 3, 687.8 US$

La utilidad se considera como el 30% del costo de generacin de vapor y se marca negativa por
representar una ganancia y no un costo.

Evaluando el costo total del proyecto a 10 aos incluyendo los gastos variables generados en ese tiempo,
se tiene:

Pgina 30

Proyecto Ingeniera de calor

RESULTADOS Y CONCLUSIONES

Pgina 31

Proyecto Ingeniera de calor

RESULTADOS

Pgina 32

Proyecto Ingeniera de calor

Temperatura del crudo al final de cada etapa

12
10
8

Temperatura (F)

T de salida

4
2
0

Corriente

La temperatura mxima alcanzada

intercambiadores fue de 393.61F

en

el

ltimo

arreglo

de

CONCLUSIONES
La eficiencia ptima de un intercambiador de calor requiere, de un buen
diseo y una apropiada seleccin de los materiales para su fabricacin.
Para lograr esto, es conveniente contar con una metodologa, puesto
que la investigacin cientfica se define como la serie de pasos que
conducen a la bsqueda de conocimientos o desarrollo de proyectos
cientficos y tecnolgicos aplicando mtodos y tcnicas. El mtodo para
el desarrollo del conocimiento cientfico es un procedimiento riguroso,
de orden lgico, con el propsito de demostrar una teora o resolver un
problema prctico.
Este trabajo muestra que desarrollando una metodologa es posible
realizar proyectos de desarrollo tecnolgico, que tanto se necesitan en
nuestro pas, obteniendo resultados ms eficientes y confiables, la
metodologa es una herramienta que relaciona una propuesta de
investigacin o un problema a resolver con los resultados.
CONSIDERACIONES
Topologa:
Tubos: Petrleo debido a su mayor viscosidad y para facilitar limpieza
del equipo
Coraza: Corrientes de destilacin
Arreglo triangular: Provoca ms turbulencia y por lo tanto una mayor
trasferencia
Pgina 33

Proyecto Ingeniera de calor

Recomendaciones TEMA:
Ds
B Ds
5
50000 Gs 300000
50000 > 100000

v2
Lquidos en general

500

Mximo tolerable

4000

Pala lquidos abrasivos

1500

Presiones
pt < p max = 7 psia
ps < p max = 10 psia
Cabezal fijo: Mas barato, por lo tanto el interior de los tubos puede ser
limpiado usando un cepillo rotatorio o un cepillo de alambre
Material: Acero al carbn debido a que es ms econmico

Pgina 34