Capítulo 1
Capítulo 1
Capítulo 1
SEP
D.G.E.S.T.
INSTITUTO TECNOLGICO
DE CERRO AZUL
CAPTULO 1 PROPIEDADES
DE LOS FLUIDOS
CATEDRTICO:
ING. FRANCISCO BARRAGN OLIVEROS.
MATERIA:
FLUJO MULTIFSICO EN TUBERAS.
EQUIPO 2:
BAUTISTA CRUZ ADILENE.
FRANCISCO MORALES YENNI IBETH.
OLVERA ORTIZ JOS ALBERTO.
PREZ REYES DANIEL.
REYES SANTIAGO BLANCA XOCHITL.
SEMESTRE: 6
GRUPO: 1
1.1
CAPTULO 1
PROPIEDADES DE LOS FLUIDOS
INTRODUCCIN
ideal, es decir
compresibilidad.
siguiente manera:
n
gi q gi
g = i=1
qg
(1.1)
Donde:
n=
gi =
q gi =
i ( pies g3 a c . s ./d a ) .
( g =0.55 ) . Al
lb m
d a.
lbm
d a.
lb m
d a.
lb g
W g=R
)(
pie g a c . s .
bl a c . s
lb aire
qo o
x 0.0764
x g
bl o a c . s .
da
pie 3aire a c . s
pie g a c . s .
lb aire
pie
3
aire
a c.s.
W g=0.0764 R q o g (1.4)
lb gd
W gd=R s
0.0764
)(
pie gd a c . s .
bl a c . s
pie gd a c . s .
qo o
gd
x
bl o a c . s .
d a
lb aire
lb aire
pie aire a c . s
3
pie 3aire a c . s .
W gd=0.0764 R s qo gd (1.5)
lb gf
W gf =( RRs )
0.0764
pie 3gf a c . s .
bl a c . s
qo o
x gf
bl o a c . s .
da
lb aire
pie aire a c . s
3
)(
pie 3gf a c . s .
lb aire
pie
3
aire
ac. s.
gf
R g Rs gd
(1.8)
RR s
gd
R> Rs .
Por
Rs .
y sustituir
gd
por
si
g < gd
0.0764 gf
(1.10)
Bg
2.7044 p gf
(1.11)
Z (T + 460)
T +460
(1.16)
T pc
p pr =
p
(1.17)
p pc
r=0.27
p pr
(1.18)
Z T pr
Z =1+ ( A1 + A2 /T pr + A 3 /T pr
5
A5
2 A 5 A6 r
+ ( A7 2r /T 3pr )( 1+ A 8 2r ) exp ( A 8 2r ) (1.19)
) r + A 4 + T r + T
pr
pr
Donde:
A 1=0.31506
A 5=0.6123
A 2=1.0467
A 6=0.10489
A 3=0.5783
A 7=0.68157
A 4 =0.5353
A 8=0.68446
11
para
lb | |
2
pg
T pc T pc 3 (1.20)
p pc=
p pc T pc
(T pc + y H 2 S ( 1 y H 2 S ) 3 )
(1.21)
la siguiente manera:
1.6
0.5
4
3 ( R )=120 ( y CO20.9
H S y CO 2H S ) +15 ( y H S y H
2
)(1.22)
Donde:
y CO
H S =es la suma de las fracciones molares de CO 2 y H 2 S .
2
12
CO2 y H 2 S
de
Los valores de
T pc y p pc
T pc y p pc
((
g
g =K x 10 exp X
62.428
4
) )(1.23)
1.5
K=
209+550.4 gf + (T + 460 )
(1.24)
X =3.5+
986
+0.2897 gf (1.26)
( T +460 )
(1.27)
3
(1.28)
(1.29)
(1.309
13
Dnde:
g = es la viscosidad del gas natural calculada con la ec. (1.23)
C N 2 = es la correccin por presencia de
CCO 2 = es la correccin por
N2 .
C 02 .
N2 .
C 02 .
RS
p=18
10
gd
( )
(1.31)
14
RS
de la
RS = gd
((
))
(1.32)
B O=0.972+0.000147 (F)
(1.33)
Dnde:
F=RS
gd
+125 T
( )
(1.34)
RS
(1.35)
API
)
T + 460
(1.36)
API30
API>30
15
C1
0.0362
0.0178
C2
1.0937
1.1870
C3
25.724
23.931
API
)(C1 ,C 2 ,C 3 R S )
gs
(1.37)
COEFICIENTE
C1
API30
4
4.667* 10
API>30
4
4.67* 10
C2
1.751* 10
1.1* 10
C3
8
-1.811* 10
8
-1.337* 10
RS
1. Calcule p con:
log p=2.67364+2.35772 log p0.703988 log 2 p +0.098479 log 3 p
2. Calcule
RS
con:
16
(1.38)
RS = gd
p API 0.989
Ta
1
0.816
(1.39)
Dnde:
a= 0.130, para aceites voltiles.
a=0.172, para aceites negros.
Para lograr un mejor ajuste se puede variar el valor al exponente a.
3. Calcule
B O=R S
BO
con:
gd
0.526+ 0,968
O
( )
(1.40)
4. Determine
BO
con:
B O=1+10
(1.41)
Dnde:
a=0.58511+2.91329 log BO 0.27683 log 2 BO ( 1.42)
1.3.4 CORRELACIN DE J.A LASATER
La correlacin de Lasater se basa en 158 medidas experimentales de la
presin en el punto de burbujeo de 137 sistemas independientes,
producidos en Canad, en el centro y oeste de los estados unidos y
Amrica del sur. El error promedio en la representacin algebraica es del
3.8% y el mximo error encontrado es del 14.7%.
Las ecuaciones siguientes corresponden a la correlacin de Lasater para
un aceite saturado:
p=
pf (T + 460)
(1.4)
gd
17
Donde
pf
yg
R s /373.3
Rs
+ 350
379.3
(1.45)
El peso molecular del aceite en el tanque se correlaciono con los API del
aceite en el tanque de almacenamiento, a cuya curva se le ajustaron las
siguientes expresiones:
Si
Si
40 API <55, M O=
0.0996
( 1048.33
API )
RS
(1.46)
(1.47)
Y g O
( 1Y g ) Mo
(1.48)
pf
, se le ajusto la siguiente
ecuacin:
5
y g =419.54510 p f 591.42810
p f +334.51910 p f +169.87910
(1.49)
18
62.4 +0.01362 R S gd
Bo
(1.50)
(1.51)
(1.52)
(1.53)
om=10 x 1
(1.54)
x=YT 1.63
(1.55)
Y =10 Z
(1.56)
Z =3.03240.02023 API
(1.57)
(1.58)
19
(1.59)
(1.60)
( pbp )m
(1.61)
m=C 1 Pc 2 exp ( C 3+ C 4 p)
(1.62)
(1.63)
pb=18
Vzquez:
Oisten:
pb =
pb=10
(( )
R
(0.00091T0.0125 API )
0.83 10
gd
API
)
T +246
C 1 gd
R exp (c 3
(1.64)
1/ C 2
(1.65)
(1.7669+1.7447logPb0.30218 logPb)
(1.66)
20
Donde:
p
0.816
( R )
b=
gd
0.130
T
0.989
API
pb=
Lasater:
Donde
(1.67)
pf
Pf (t+ 460)
gd
(1.68)
(1.69)
62.43
Bw
(1.70)
w = A+
B
(1.71)
T
Dnde:
A= (-0.04518 + 0.0009313 (%NaCl) 0.000393 (%NaCl)2
(1.72)
B= 70.634 + 0.09576 (%NaCl)2
(1.73)
12
p (T 40)
(1.74)
280T
( w2 w 1 ) + w1 (1.75)
206
w 1=52.50.006 p(1.76)
w 2=76 exp (0.00025 p ) (1.77)
Donde:
w 1 Es la tensin superficial agua-gas a 280F.
w 2 Es la tensin superficial agua-gas a 74F.
22
con
las
siguientes
(1.78)
(1.79)
(1.80)
S= P (A+BT+CT+DT3)
Donde:
A=3.69051
(1.81)
B=0.08746 C=0.01129 D=-0.00647
Rsw
(1.84)
Donde:
B=-0.01052 + 4.77 X 10 -7p
(1.87)
23
24
Rs
B od
B
( of ),
B
( od )
por el
B
( odb) , Fig. 1.2.
26
B of
B ofb
(1.91)
B odb
( Rsf ) ,
en el tanque de
almacenamiento
y
en
el
separador, en funcin de la
presin del yacimiento. Como se
muestra en la fig. 1.3, los valores
de la relacin de solubilidad
varan para una presin dada,
dependiendo del proceso de
separacin
utilizado.
Estos
valores estn expresados en pies
cbicos a condiciones estndar por barril.
27
28
H 2 S =0.2327
API =30
R=750 (
pies 3
)
blo
lb | |
pg 2
ps=120
T s=75( F)
B ofb =1.455(b l o a c . b ./b l o a c . s .)
Rsfb =750( pies3 /b l o )
lb / p g 2||
pb=3,300
T b=220( F)
lb/ p g2||
p=2,000
T 200( F )
El primer paso, es seleccionar las correlaciones que mejor determinan
Pb
R
Bo
los valores de
, s
y
con sus respectivos factores de
correccin.
Para hacer las comparaciones, se tomaran los datos de laboratorio
,
B ofb
Pb
Rsfb
( 0.85
0.83
29
30
1
220+240 1.0937
=3,626
0.0362 0.8522
Pb=
Pb=
30
=610
220+ 460
(10.6850 ) 336.41
0.6850 0.8762
Rs =132,755
=752
30
30
( 0.85222
) (1.751 10
)+ 1.811 1 08 750=1.3729
Bolab
(1.94 )
B oc
B oc
(1.93), se tiene:
B =
B olab
B (1.95)
B oc oc
( )
Donde:
B oc =
B oc =
Es un valor cualquiera de
Es el valor de
Bo
Bo
Bo
el factor de correlacin.
B olab=
B =
Es el valor de
Es el valor de
Bo
corregido.
31
B olab
= Factor de correccin.
B oc
( )
Bo
, con la segunda.
gd = g '
por lo que
gd
Rs
, a las condiciones de
2. Se supone un valor de
gd .
32
Rs =456
Pf .
con la Ec.(1.49).
gd
Rs
gds = gd c
y se contina a partir
del punto 2.
p < Pb o Rs
R , el aceite
est saturado.
Aplicando las Ec.(1.8) y las Ecs. (1.12),(1.13),(1.16)-(1.22) por existir
gases amargos:
gf =
P pc=662.69
390.09
=609.29
419.13+ 0.2327 ( 10.2327 ) 29.04
33
T pr =
200+ 460
=1.6919
390.09
P pr=
2,000
=3.2828
609.24
z=0.86067
Bg y g
B g=
g=
0.0764 0.7962
=7.5813
0.0080236
960
+0.2897 0.7962=5.1852
(200+ 460)
Y =2.40.2 5.1852=1.3630
1.5
200+ 460
5.1852
7.5813
( 62.428
)=0.017145
1.363
CO2 y H 2 S
C N =0
2
34
g =0.018034
Clculo de las Propiedades del Aceite
El factor de volumen del aceite se obtiene con las Ecs.(1.40)-(1.42):
0.8847 0.526
+0.968 ( 200 ) =651.9215
0.8762
B o=456
B o=( 1+1 00.57895 ) C B =1.28363
0
Y =1 0
=266.38
200 1.163=0.56157
X =266.38
0.56157
om=1 0
1=2.6439
456+ 1500.338=0.6239
b=5.44
456+ 1000.515=0.4133
a=10.715
0.6239
2.6439
=0.7581
o=0.4133
Aplicando la Ec. (1.58), se obtiene la tensin superficial del aceite:
35
62.43
=61.795
1.0103
70.637
=0.30967
200
( 280200 )
( 46.09640.5 ) 40.5=42.67
206
=0.5847
( 2,000
2,276 )
36
5
T = ( 20032 )=93.333
9
'
T=
93.3390
=0.333
10
2
Ejemplo 1.2
Anotar en una tabla las variaciones numricas de
Rs
gd , gf
Bo
Rs
gd
, gf
3,300
750.26
0.8500
3000
701.58
0.8558
0.7660
2,500
580.19
0.8701
0.7813
2,000
456.24
0.8847
0.7961
1,500
331.02
0.8995
0.8109
1,000
209.60
0.9138
0.8253
500
99.65
0.9268
0.8382
250
50.60
0.9326
0.8440
100
23.28
0.9358
0.8473
50
14.53
0.9369
0.8483
25
10.22
0.9374
0.8488
15
8.51
0.9376
0.8490
Afectados por su correspondiente factor de correcin.
37
B o
1.4149
1.3504
1.2851
1.2199
1.1579
1.1034
1.0798
1.0669
1.0627
1.0607
1.0599
38
REFERENCIAS
1. Kermit Brown: Fluid Property Options, Technology of Artificial Lift
Methods.
2. Katz, D.L.: Prediction of the Shrinkage of Crude Oils, Drilling and
Prod. Prac. API. 1942.
3. Benedict, M., et. Al, An Empirical Equation for Thermodynamic
Properties of Light Hydrocarbons and Their Mixtures. J. Chem.
Phys. Vol. 8. 1940.
4. Standing, M.B. y Katz, D.L.: Density of Natural Gases, Trans. AIME
(1942) 140-149.
5. Wichert, E. y Aziz, K.: Calculate Zs for Sour Gases, Hydrocarbon
Processing, May, 1972.
6. Lee, A.L., et. Al.: The Viscosity of Hydrocarbon Gases Under
Pressure, J.P.T., October, 1954.
7. Carr, N.L., et. Al.: Viscosity of Hydrocarbon Gases Under
Pressure, J.P.T., October, 1954.
8. Standing, M.B.: A Pressure-Volume-Temperature Correlation for
Mixtures of California Oil and Gases, Drill, and Prod. Prac., API
(1947) 275-286.
9. Vzquez, M. y Beggs, M.D:: Correlations for Fluid Physical Property
Prediction, J.P.T. June, 1980.
10.
Oistein, Glaso.: Generalized Pressure- Volume-Temperature
Correlations, J.P.T. May, 1980.
11.
Lasater,
J.A.:
Bubble
Point
Pressure
Correlation,
Transactions AIME, Vol. 13, 1958.
12.
Beggs, H.D. y Robnson, J.R.: Estimating the Viscosity of
Crude Oil Systems, J.P.T. September, 1975.
13.
Baker, O.: Designing Pipelines for Simultaneous flow of Oil
and Gas, Pipeline Engineer, February, 1960.
14.
Dodson,
C.R.
y
Standing,
M.B.:
Pressure-VolumeTemperature and Solubility Relations for Natural-Gas-WaterMixtures, Drill, and Prod. Prac., API (1944) 173-179.
15.
Matthews, C.S. y Rusell, D.G.: Pressure Build-up and Flow
Test in Wells, Monograph Series SPE (1967).
39
16.
Philip L. Moses.: Engineering Applications of Phase Behavior
of Crude Oil and Condensate Systems, Journal of Petroleum
Technology, July 1986.
17. Rollins, J.B., et.al: Estimation of Solution GOR of Black Oils, J.P.T.
January, 1990.
40