Correlaciones PVT-Carlos Banzer
Correlaciones PVT-Carlos Banzer
Correlaciones PVT-Carlos Banzer
11111116111, l
1
}
i
'' . .. ,
,,
....... ,
• ·:1:-,
u.
1
1
-UCA
·-
u,J-1
1:
... , 1998
1rre ac 1111 111 ricas
111.11r1111ttnz1r s.
Elllcllll Pr1ll•l11r
•
CORRELACIONES NUMÉRICAS P.V.T.
11radec1n11ento
El autor expresa su sincero agradecimiento al Instituto de
Investigaciones Petroleras de la Universidad del Zulia
(INPELUZ), y a la Fundación Laboratorio de Servicios
Técnicos Petroleros (FLSTP), a nombre de sus Directores:
Prof. Renato Acosta (Director-Presidente INPELUZ-
FLSTP) e lng. Carmelo Urdaneta Galué (Director-Gerente
FLSTP), por haberme facilitado los materiales y equipos
necesarios para la elaboración de este trabajo.
Especial agradecimiento al Prof. Douglas A. Alvarado,
quien de manera desinteresada y generosa me proporcionó
parte del material original que se presenta.
Finalmente, especial agradecimiento a todas aquellas
personas que puedan apreciar la modesta importancia de
este esfuerzo.
I
-~ ('1 ,\ ,. ..___...~
Ing, Carlos Bánzer S.
Maracaibo, enero de 1996
11
CORRELACIONES NUMÉRICAS P.V.T.
ontenldo
Página
Agradecimiento-· 11
Nomenclatura VI
Introducción 1
Correlaciones P.V.T. 1
111
CORRELACIONES NUMÉRICAS P.V.T.
----------------------------------------------------------------
Página
Propiedades Físicas del Petróleo 48
Gravedad Específica del Petróleo, r O
48
Presión de Burbujeo, pb 48
Correlación de Standing. M.B. 49
Correlación de Lasater, J.A. 51
Correlación de Vázquez, M.E. y Beggs, H.D. 54
Correlación de Glaso, O. 55
Correlación de la TOTAL, C.F.P. 57
Correlación de AI-Marhoun, M.A. 57
Correlación de Dokla, M.E. y Osman, M.E. 58
Correlación de Petrosky, G.E., Jr. y Farshad, F.F. 59
Correlación de Kartoatmodjo, T. y Schmidt, Z. 59
Razón Gas Disuelto-Petróleo, R5 63
Correlación de Standing, M.B. 65
Correlación de Lasater, J.A. 65
Correlación de Vázquez, M.E. y Beggs, H.D. 66
Correlación de Glaso, O. 66
Correlación de la TOTAL, C.F.P. 66
Correlación de AI-Marhoun, M.A. 67
Correlación de Dokla, M.E. y Osman, M.E. 67
Correlación de Petrosky, G.E., Jr. y Farshad, F.F. 67
Correlación de Kartoatmodjo, T. y Schmidt, Z. 68
Factor Volumétrico del Petróleo, 80 72
Correlación de Standing, M.B. 72
Correlación de Vázquez, M.E. y Beggs, H.D. 74
Correlación de Glaso, O. 75
Correlación de la TOTAL, C.F.P. 75
Correlación de AI-Marhoun, M.A. 75
Correlación de Dokla, M.E. y Osman, M.E. 76
Correlación de Petrosky, G.E., Jr. y Farshad, F.F. 76
Correlación de Kartoatmodjo, T. y Schmidt, Z. 76
Factor Volumétrico Total, 81 80
Correlación de Glaso, O. 82
Correlación de AI-Marhoun, M.A. 82
Compresibilidad del Petróleo, c0 83
Correlación de Vázquez, M.E. y Beggs, H.D. 84
Correlación de Petrosky, G.E., Jr. y Farshad, F.F. 85
Correlación de Kartoatmodjo, T. y Schmidt, Z. 85
Correlación de McCain, W.D., Jr., Rollins, J.B. y Villena-Lanzi, A.J. 86
Viscosidad del Petróleo, µ0 93
Petróleo Muerto
Correlación de Beal, C. 98
Correlación de Beggs. H.D. y Robinson, J.R. 99
~----------------------------------------------------------~-
IV
CORRELACIONES NUMÉRICAS P.V.T.
Página
Correlación de Glaso, O. 100
Correlación de Egbogah, E.O. 100
Correlación de Kartoatmodjo, T. y Schmidt, Z. 100
Petróleo Saturado
Correlación de Chew, J.N. y Connally, C.A., Jr. 101
Correlación de Beggs, H.D. y Robinson, J.R. 103
Correlación de Kartoatmodjo, T. y Schmidt, Z. 103
Petróleo Subsaturado
Correlación de Beal, C. 103
Correlación de Vázquez, M.E. y Beggs, H.D. 105
Correlación de Kartoatmodjo. T. y Schmidt, Z. 105
Densidad del Petróleo, Po 108
Tensión lnterfacial Gas-Petróleo, ergo 110
Correlación de Baker, O. y Swerdloff, W. 110
Referencias 111
V
CORRELACIONES NUMÉRICAS P.V.T.
-----------------------------------
Nomenclatura
Simbolo Descripción Unidades
~-------------------------------------------·------~-
VI
CORRELACIONES NUMÉRICAS P.V.T.
VII
CORRELACIONES NUMÉRICAS P.V.T.
11uoduccló1
En el análisis del comportamiento de yacimientos, cálculo de reservas y diseño de
equipos, se requiere el conocimiento de las propiedades fisicas de los fluidos. Estas
propiedades normalmente se determinan en el laboratorio mediante el análisis de muestras de
fluidos tomadas en el fondo de los pozos o mediante una apropiada recombinación de muestras
tomadas en superficie.
El conjunto de pruebas necesarias para determinar estas propiedades se denomina
análisis Presión-Volumen-Temperatura, P.V.T., como comúnmente se le llama, y consiste en
determinar las relaciones entre presión, volumen y temperatura para una mezcla de
hidrocarburos (líquido y gas) en particular. ·
Muchas veces sin embargo, no se dispone de información experimental, debido a que no
se pueden obtener muestras representativas de los fluidos o porque el horizonte productor no
garantiza el gasto en realizar un análisis P.V.T. de los fluidos del yacimiento. En estos casos,
las propiedades fisicas de los fluidos deben ser determinadas por analogía o mediante el uso de
correlaciones empíricas.
En el pasado, las correlaciones P.V.T. fueron presentadas en forma tabular y/o gráfica;
sin embargo, con la aparición de las calculadoras manuales programables y las computadoras
personales, tales correlaciones han sido reducidas a simples ecuaciones numéricas o
expresiones analíticas con el propósito de utilizarlas en programas de computación.
En este sentido, el objetivo de este trabajo es presentar una revisión de las propiedades
fisicas de los fluidos y las correlaciones que se han desarrollado y reducido a la forma
anteriormente mencionada para su determinación. Adicionalmente, se presentan ejemplos
numéricos con el propósito de ilustrar el uso de las correlaciones sin intentar establecer
criterios comparativos entre las mismas.
Las correlaciones P.V.T. son desarrolladas a partir de datos de laboratorio y/o de campo
y formuladas de manera que puedan ser utilizadas con datos obtenidos sin mucha pérdida de
tiempo y/o inversión de esfuerzo. Estos datos son:
a) Propiedades de los fluidos
• Gravedad específica del petróleo, r O
• Gravedad específica del gas, r g_
• Razón gas producido-petróleo, R
b) Propiedades del yacimiento
• Presión,p
• Temperatura, T
Antes de determinar una propiedad, se debe asegurar que la aplicación de interés se
encuentra dentro del rango de los datos para la cual la correlación fue desarrollada. Una vez
hecho esto, la correlación utilizada será adecuada y se obtendrán resultados confiables.
Carlos Bánzer S.
1
PROPIEDADES FISICAS DEL GAS
-------------------------------------------------------------------
El gas natural puede ocurrir como tal, en yacimientos de gas libre o asociado con
yacimientos de petróleo y de condensado. En vista de que el gas natural normalmente se
encuentra presente en los yacimientos de hidrocarburos, es de gran importancia el
conocimiento de ciertas propiedades fisicas del mismo que son fundamentales en el análisis del
comportamiento de yacimientos de petróleo, gas y condensado.
En esta sección se presenta una revisión de las propiedades fisicas de gases naturales
hidrocarburos y las correlaciones numéricas y analíticas que se han desarrollado para s11
determinación.
Las relaciones que describen el comportamiento Presión-Volumen-Temperatura, P.V.T.
de gases son llamadas ecuaciones de estado. La ecuación de estado más simple, es llamada 1~,
ley de los gases ideales••:
p V = nRT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ( 1.1 )
donde V es el volumen de n moles de gas a la presión p y temperatura T, y R la constante
universal de los gases.
Esta ecuación es el resultado de los esfuerzos combinados de Boyle, Charles, Avogadro
y Gay Lussac, y es solo aplicable a presiones cercanas a la presión atmosférica para. la cual fue
experimentalmente obtenida y a las cuales los gases se comportan como ideales.
Debido a que n es igual al peso total de gas, m, dividido por el peso molecular, M, la Ec.
( 1.1) puede escribirse como:
pV = .:. RT ( 1.2)
Las Ecs. (1.1) y (1.2) son formas de la ecuación general de los gases ideales.
El valor de la constante R depende de las unidades empleadas de presión, volumen )'
temperatura. Por ejemplo, considérese p en libras por pulgada cuadrada absoluta (lpca), Ve11
pies cúbicos, Ten grados Rankine (ºR) y los moles de gas en lbs-mol. La hipótesis de:
Avogrado establece que 1 lb-mol (n = 1) de cualquier gas ideal ocupa un volumen de 379.4
pies cúbicos a 14.7 lpca. y 60 ºF. Por lo tanto:
R = PV = (14.7 lpca)(379.4 p.cu.) = 10.73 lpca - p.cu.
nT (1 lb- mol)(520 ºR) lb- mol ºR
Para otras unidades de p, V, Ty n, la Tabla 1.1 presenta diferentes valores de R.
Gravedad Especifica del Gas, y9:
La gravedad específica del gas se define como la razón de la densidad del gas a la
densidad del aire, ambas medidas a las mismas condiciones de presión y temperatura. Como la
gravedad del gas se mide generalmente a condiciones normales (14.7 lpca y 60 ºF), tanto el gas
como el aire se comportan como gases ideales y puede escribirse:
Pg
rr = p •••••••••••••••••••••••••••••••••••••••••••••••••••••••••••• ( 1.3)
O
Carlos Bánzer S.
*Referencias al final de la sección 2
CORRELACIONES NUMÉRICAS P.V.T.
p V T n R
La densidad de un gas se define como la masa por unidad de volumen, así, de la Ec.(1.2)
se tiene:
m pM
p8 =- = , ( 1.4)
V RT
Reemplazando la Ec. (1.4) en la Ec. (1.3) para el gas y el aire a las mísmas condiciones
de presión y temperatura:
pMg
r = RT = Ms = M .......................... ,, ( 1.5)
g P~ª- Mª 28.96
RT
donde;
Mg : Peso molecular del gas, lbs/lb-mol o grs/gr-mol
M0 Peso molecular del aire, ( = 28.96 lbs/lb-mol)
:
Si no se conoce el peso molecular del gas, pero se conocen los pesos moleculares de las
fracciones molares de la mezcla, M¡, se determina el peso molecular aparente de la misma
mediante:
,,
M = LY; M; ( 1.6)
;.J
donde;
M; : Peso molecular del eomponente i en la mezcla, lbs/lb-mol o grs/gr-mol
Y;: Componente i en la mezcla, fracción molar
n : Número de componentes en la mezcla
Así, una vez obtenido el peso molecular aparente de la mezcla mediante la Ec. ( 1.6), se
determina la gravedad específica de la misma mediante la Ec, ( 1.5). La Tabla 1.2 presenta
algunas propiedades flsicas de gases hidrocarburos párafínicos puros y de otros componentes
comúnmente encontrados en el gas natural.
Presión y Temperatura Seudocriticas de Gases Naturales, Psc y T.c:
Con el propósito de estudiar el comportamiento de mezclas gaseosas, se ha introducido
el concepto de presión y temperatura seudocríticas, Psc y Tsc·
Carlos Bánzer S.
3
PROPIEDADES FISICAS DEL GAS
-------------------------------------------------------------------------------
1 ... .........
11
(!i &, • 1 ...
f8M~§ ...... •an to ~ :.. •
il00 § (O ... i
&t)
......
... ...
(')
w .......
U)
•
co • • •
~
• • •
~ ~~
• ao• .....•
~
•
o o o • • • •
ºººº
•
o •
00
ººº
!o E:.... ¡o ~.....
• • • •
i~ . . . B:t
g~~~~
• • • •
~~
• •
o o o o
•
C') 00
~-.! ...
co• CIO• t') CO ~ O> CIO CO O O O> CD
.... 1~
o• o•
....
,~
I'- ..-
:§~-!~í~~~!!
• "1111"
is
•
CD (W)
CD
...... co
~
a.
(W) ....
~ ~~
.... "1111" ~
~
·e -
o ax
' • o
~._.
~
...
lt)
00
•
o
....
•
~i
! x• ~N x<» :é'
O O N N
:é' :é'lt> :é' :é'<O (.):é' :é'CIO :é'O, ~.- <;¡N
.. <O m CIO ~ N
oN d" ! u.
i "':e z •
O O O ~
'llllt •
C")
º~~~ lt)
e <i 'E ~ Z :C O
N NON
2
:sI!
~
e
o
o e
.1
e
!
!5 t :g
.e
m
::::,
foo
G) ~
o -8 g o
"O e
1 g o o
o
2 ~ e o e o
e ca e e o
e o e §, e :2
fi ~ e ! )(
J! ca
g g:eg~- ... ~ 1c.
¡; J!
2w
8" ~
ñ: ~ .
G) Q)
Q...
·..&.
Q)
~!
t~ é ~~ :! JJ ! d -ca
~ ~ ~
:'-
Cario• Banzer S.
4
CORRELACIONES NUMÉRICAS P.V.T.
Una serie de métodos o reglas de combinación han sido propuestos para determinar estas
propiedades seudocríticas cuando se conoce la composición de la mezcla gaseosa. Estos
métodos o reglas de combinación suponen que cada componente contribuye a la presión y
temperatura seudocríticas en proporción a su porcentaje de volumen en el gas y a la presión y
temperatura críticas, respectivamente, de dicho componente.
Entre los métodos o reglas de combinación más conocidos se tienen los siguientes:
Método de Kay, W.B.3:
La presión y temperatura seudocríticas están dadas por:
11
donde;
PscM : Presión seudocrítica de la mezcla, lpca.
TscM : Temperatura seudocrítica de la mezcla, ºR
Pci : Presión crítica del componente i, )pea.
Te; : Temperatura crítica del componente i, ºR
y; : Componente i en la mezcla, fracción molar
n : Número de componentes en la mezcla
Método de Stewart, W.F., Burkhardt, S.F. y Voo, D.':
Las ecuaciones de este método son:
2
1 n 2 n
J=- LY; (~/Pe);+- LY; (J;. / Pc);112 ............................ ( 1.9)
3 i=l 3 i=I
a.t»,
11
1/2
K= LY; )¡ (1.10)
donde;
K2
T_,cM = J (J.11)
Sutton5, modificó el método de Stewart y Cols. para aplicarlo a gases naturales con alto
peso molecular (y g > 1.25) . Los términos J y K fueron ajustados mediante los siguientes
factores:
J'=J-eJ (1.13)
K'=K-eK (1.14)
donde;
2 1
s J = 0.608 lFJ + l.1325FJ - l 4.004FJ Y<·1+ + 64.434.F:, Yc·1+ •••••••••••••• ( 1.15)
Carlos Bánzer S.
5
PROPIEDADES FISICAS DEL GAS
------------------------------------------------------------------
y
1 2 ,, 2
LY ( LY ( ~ / Pe)
ti
112
F¡ = - t; /Pe) +- .....................( 1.17)
3 , •• c·1+ 3 /:al c.·1+
•. -
,
' '
• ,
•
'
~
'I
'
'
500 1 'J
_í
'
.
1
.
. -
....
_,;._
~·
• ' 1
450
"~~
• • - '
~~
.
1
- • •
G~ .. - -
. 1 l
. 1 •
.
--~;~
:
'. ¿
1 1 1
•
:,,.. r;J- ~tlO '
~ G~S j
' •
• .
•' ' '
'1 1
•
. ' •
300 l l
'
0.5 0.6 0.7 o.a 0.9 1.0 1.1 1.2
GRAVEDAD DEL GAS , AIRE • 1
Fig. 1.1.- Presión y temperatura seudocrfticas de gases naturales (Ref. 6).
Carlos Bánzer S.
6
CORRELACIONES NUMÉRICAS P.V.T.
Las curvas en esta figura representan mezclas gaseosas con cantidades despreciables de
N2, C02 y H2S. Si la mezcla gaseosa contiene cantidades apreciables de componentes no
hidrocarburos (mayor de 5% por volumen), la gravedad específica de la mezcla debe ser
corregida por estos componentes. Así, si y gM es la gravedad específica de la mezcla total de
gas, la gravedad específica de la porción de gas hidrocarburo, y a« , está dada por7:
y gM - 0.967 YN 2 - 1.52Yco 2 - l.18yH2s
Y gHC = l ( 1.18)
- YN2 - Yco2 - YH2s
donde;
r gHC Gravedad
: específica de la porción de gas hidrocarburo, (aire= 1)
r gM : Gravedad específica de la mezcla total de gas, (aire= 1) ·
Y N2 : Contenido de N2, fracción molar
Y C02 : Contenido de C02, fracción molar
y82s : Contenido de H2S, fracción molar
En la Ec. ( 1.18), r gHC ~ O.SS (gravedad específica del metano).
7:
Un ajuste de las curvas de la Fig. 1.1, está dado mediante las siguientes ecuaciones
Gas Natural:
PscHC = 677 + 15y gH(~ -37.Sy gH(:2 ••••••••.•••••••••••• ,, •••••••••••••••• (1.19)
2
T.rcH(: = 168 + 325y gff(,' - 12.Sy gH(.' ······:············,, •••••••••••••••• ( 1.20)
Gas Condensado:
2
P.,·cff(,' = 706 - 51. 7 y gH(.' - 11. ly gH(~ ••••• : ( 1.21)
2
T:·cH(: = 169 .2 + 349 .Sy gH(~ - 74y gHc· ••••••••••••••.•••••••••••••••• ( 1.24)
1'.fcM = (1- YN2 - Yco2 - YH2S) J:cHc.· + 227YN2 + 548Yco2 + 672yH2s ·········· (1.26)
En estas ecuaciones, la presión y temperatura seudocríticas de la porción de gas
hidrocarburo PscHC y TscHc, pueden ser representadas por las ecuaciones anteriores dependiendo
del tipo de gas hidrocarburo considerado.
Presión y Temperatura Seudocritlcas del Componente Pesado, Paocr+ y TacC7+:
En el análisis cromatográfico del gas natural y del gas condensado, solo se analizan los
componentes hidrocarburos del metano al hexano, C1a C6 y los componentes no hidrocarburos
Carlos Bánzer S.
7
PROPIEDADES FISICAS DEL GAS
-----------------------------------------------------------------·
N2, C02, H2S, 02. El resto, o sea los elementos más pesados que aparecen en pequeñas
proporciones como el C6+ (hexanos y componentes más pesados) o C1+ (heptanos y
componentes más pesados) se expresan como un solo componente del sistema y generalmente
son caracterizados por su peso molecular y gravedad específica.
Las siguientes correlaciones permiten determinar la presión y temperatura seudocríticas
del C7+:
Correlación de Mathews, T.A., Roland, C.H. y Katz, D.L.1:
Mathews y Cols. en base a una serie de micro-destilaciones de fracciones de C7+ de
condensados obtuvieron la correlación gráfica de la Fig. 1.2, por medio de la cual puede
determinarse la presión y temperatura seudocrfticas del C7+ en función de la gravedad
específica, y c1+ (agua = 1) y del peso molecular, Me,+· Para fracciones de pentanos y
compuestos más pesados, C5+, y hexanos y compuestos más pesados, C6+, la correlación
también puede ser utilizada.
1400
f ... ~
5 1300 o.
-
~
....a:
F= Gt, E•P·
1200 - -_......•
-
o.a
º·'º
~
:;11111 ... IIIIC
en ,.,,. -
.,,,,... ~--== o'15
.
! 1100 --"" ~
l,,illl"'
<
o::
~
- t:=:-"'
or
-
~
LaJ ~
e, 1000 -
~
LIJ ~
1-
900
100 120 140 160 180 200 220 240
PESO MOLECULAR
600
s
-e..
A.
--
u
1-
500
o:
~
400
::>
l&J
CI)
zo
-en
l¡.J
300
o::
a.
200
100 120 140 160 180 200 220 240
PESO MOLE CU LAR
Fig. 1.2.-Presión y temperatura seudocrfticas del C7+ (Ref. 8).
Carlos Bánzer S.
8
CORRELACIONES NUMÉRICAS P.V.T.
Un ajuste de las curvas de la Fig. 1.2, está dado por las siguientes ecuaciones 7:
P~,cc·1+ = 1188 - 431 log( Mc.1+ - 61.1) + [2319 - 852 log( Mc,1+ · - 53. 71 )](y e·~¡+ - 0.8). ( 1.27)
r:.ce 7+ = 608 + 364 log( Mc.·1+ - 71.2) + (2450 logMc.~7+ - 3800) logy < ?+ •••••••.••••••••••••••• ( 1.28)
donde;
PscC7+ : Presión seudocrítica del C7+, lpca.
TscC1+ : Temperatura seudocrítica del C7+, ºR
Me7+ : Peso molecular del C7+, lbs/lb-mol
r ('?+ : Gravedad específica del C7+, (agua= 1)
La presición de estas ecuaciones en el ajuste de más de 20 valores ele las curvas de
Mathews y Cols. es la siguiente':
PscC1+ TscC1+
Número de puntos 25 28
Error relativo promedio -0.51% -0.29o/o
Desviación estandar 1.61% 0.97%
'rtc{'1+ = 341. 7 + 81 ly ('7+ + (0.4244 + O. l l 74y ('7+) T,> + (0.4669 - 3 .2623y ('7 +) X
l05 .( 1.30)
T,,
donde;
PscCi+ : Presión seudocrítica del C7+, lpca.
TscC7+ : Temperatura seudocrítica del 'C7+, ºR
y ('1+ : Gravedad específica del C7+, (agua= 1)
Tb : Temperatura de ebullición, ºR
La temcPeratura de ebullición se determina mediante la siguiente ecuación desarrollada
1
por Whitson :
3
T,, = [ 4. 5 5 79 Me· 7 + o .1 5 l 7 8
y (. 7 + o. l 54 2 7 ]
•••••••••••.••.••..••••••.•••••.•••. ( 1 .3 1 )
Carlos Bánzer S.
9
PROPIEDADES FÍSICAS DEL GAS
---------------------------·---- ..
Ejemplo 1.1:
Utilizando los métodos de Kay, W.B. y de Stewart, W.F., Burkhardt, S.F. y V<>o, D. y las
correlaciones de Brown, G.G., Katz, D.L., Oberfell, G.G. y Alden, R.C. y de Sutton, R.P.
determinar la presión y temperatura seudocríticas para el gas cuya composición se presenta en
la Tabla 1.3.
Solución:
La composición se presenta en la Col. ( 2 ), y las propiedades físicas, Cols. ( 3 )-( 5 ), se
tomaron de la Tabla 1.2.
Correlación de Mathews, T.A., Roland, C.H. y Katz, D.L.
De las Ecs. ( 1.27) y ( 1.28):
Pj·c<'?+ =1188-4311og(142-61.1)+[2319-852 log(l42-53.71)](0.807-0.8)
Psce?+ = 370.3 lpca.
I:d:1+ = 608 + 364 log( 142 - 71.2) + (2450 log(l 42) - 3800) log(0.807)
(3.2623)(0.807)] los
818.87
I:·c<:?+ = 1156.80 ºR
Método de Kay, W.B.
Las sumas totales de las Cols. ( 7 ) y ( 8 ) corresponden a la presión y temperatura seudo-
críticas de la mezcla:
Psc:M = 709.2 lpca.
T.fctvi = 505.8 ºR
Método de Stewart, W.F., Burkhardt, S.F. y Voo, D.
De las sumas totales de las Cols. ( 11 ), ( 12 ) y ( 13 ) y Ecs. ( 1.9) y ( 1.1 O):
Carlos Bánzer S.
10
CORRELACIONES NUMÉRICAS P.V.T.
o
--
......
)(
~
-
~
•
N o
N
- '°" . .- ("')
~ººººº 8 º"""""
8 ..-
IO CO ("') (") ("')..-
Ml'-IOl'-1'-MMO..-MOCO
C.O O
-
N ("') I'- C.O V ..- N ......
N V
.......
00• • • • • • • • • • • • •
ºººººººººººº o
-..... -
.....
O> (") ~
LO
"IQ' ~
8
N
N M ..- N (") (") O
..-
CJ> ("') 1'- (") N O CO
- --
(O
X N ll>v ..... N..- ..... MLO
00 NOOOOOOON
o•
• • • ... $ :.. • • • • • "
N 000000000000
-
<01"'-
CIO C.O O M O C.O ..... N 1"'- ee ("') O)
11 ,r & • • • • • • • • •
...- ..-..-..-NNMMM("')•IO
-
- --
LO
--
N
-
'
CX) )( O')
o
,..._
N
- -
s, . . . .. . . . .
OMN..-r--M..-0'>1.00NIO . CX)
--
. .
-
' ' •
)( M O ..... CX) O'> O I"- "IQ' U) CD M I.O 11)
N
("') U) O') (") (")
......
..... ..... a,
s
- --
("')
- N
)(
-- L()
("')
•
-
-º-
'C
'
-~-
• c.. o::
o
t,,;.
- .
-~ ·-
~
o
~ :,.,
o
8
o •
- ~
-.....
<X)
-
N•
LL.
- -
.o ~
-
N I"'-
-
...........
~ ~
-
...... ll
,: rJ
~º.UJ
.
* *
Carlos Bánzer S.
11
PROPIEDADES FISICAS DEL GAS
------------------------------
T. = K2 = (20.2972)2 = 529.9 ºR
jcM J 0.77739
p scM =
T:·cM
529.9
= J
J
- = 68 ) . 7 pea.
0.77739 .
Aplicando la modificación de Sutton, se tiene:
De la Ec. ( 1 .1 7):
Carlos Bánzer S.
13
PROPIEDADES FISICAS DEL GAS.
----------------------------------------------------------·--------
e,,·- 1.,0
4
Q-- TEMP. 2SOºF
C(
o 1.00 ~--i1----t-----+---...---GR. ESP.=0.66,
--
.J
m
TEMP. ISOºF
......._ GR. ESP.•0.66~
PRESION, lpca
Fig. 1.3.- Efectos de presión. temperatura y composición sobre
el factor de desviación del gas.
Este factor se puede obtener experimentalmente dividiendo el volumen real de n moles.
de un gas a p y T por el volumen ideal ocupado por la misma masa de gas a iguales condiciones
dep y T.
El valor de z para diferentes gases ha sido determinado en base al Teorema de los
Estados Correspondientes, el cual establece que: a las mismas condiciones de presión }'
temperatura seudorreducidas, todos los gases tienen el mismo factor de compresibilidad z. La
presión y temperatura seudorreducidas P.rr y T.. r, están definidas por:
p
Psr = ( 1.34)
P.fcM
T
T.-;r = •••••••••••••••••••••••••••••••••••••••••••••••••,.······· (1.35)
T.1cM
donde p y T son la presión y temperatura absolutas a la cual se desea determinar z y P.,·c:M y 1;.cM
la presión y temperatura seudocríticas, tal como se discutió anteriormente.
Correlación Gráfica de Standing, M.B. y Katz, D.L.11:
Basados en el Teorema de los Estados Correspondientes, Standing y Katz desarrollaron
la correlación gráfica mostrada en la Fig. 1.4, la cual puede utilizarse para determinar el factor
de compresibilidad de un gas natural a partir de la presión y temperatura seudorreducidas,
definidas anteriormente. Las Figs. 1.5 y 1.6 presentan factores de compresibilidad de gases
naturales a bajas presiones seudorreducidas y a presiones cercanas a la presión atmosférica.
Las Figs. 1.4, 1.5 y 1.6 fueron desarrolladas independientemente, cada una basada sobre
un grupo de datos diferentes. Por lo tanto, factores de compresibilidad obtenidos de estas
figuras a las mismas condiciones seudorreducidas pueden ser ligeramente diferentes.
El método gráfico de Standing y Katz ha tenido bastante aceptación en la industria del
petróleo por dos razones básicas: Exactitud dentro de un 3% en relación a los valores experi-
14
CORRELACIONES NUMÉRICAS P.V.T.
ftMPIIIATUM
t.t
~~~ ...... ++++++ ............ NUDORRIOUCIDA t+H+H-t++t+++~io+!
1.0 t.O
&01
u
0.95
N
&Q
. N
•
g 0.7 t.6 ~
:1 :::1
-
ii5
fl) -w
15
O)
f 0.6
:E
1.5 f
:E
o o
u o
w
Q
w
Q
0.4
0.3~t:m:ttttm ~
0.25 .................................
a 9 10 11 12 13 14
PRESION SEUDORREDUCIDA , Par
•
<D
•
•
•
..., . ds:
en
d
-
•
l.t)
' .. . •
...............
•
·-
u.
C)
Carlos s,nzer S.
16
CORRELACIONES NUMÉRICAS P.V.T.
mentales de z y facilidad en los cálculos. Para tener buenos resultados con este método se debe
tener en cuenta las siguientes limitaciones del mismo:
1. El gas debe ser rico en metano ( C1 > 90% ).
2. El gas no debe tener hidrocarburos aromáticos.
3. El gas no debe tener impurezas. El contenido de 20% de N2 produce un error del 4%.
La presencia de C02 produce un error en el cálculo de z igual al valor del porcentaje
de C02 en la mezcla.
4. No presenta buenos resultados a presiones y temperaturas cercanas a la crítica.
5. No se recomienda su uso en el cálculo de z a presiones mayores de 1 ().000 lpca.
Cuando el gas natural contiene significantes fracciones de dióxido de carbono (C02) y/o
sulfuro de hidrógeno (H2S), se recomienda utilizar la corrección hecha por Wichert y Aziz•2.13
al método de Standing y Katz. Esta conección consiste en calcular la presión y temperatura
seudocríticas utilizando las siguientes ecuaciones:
s = 120 [(Yco2 + YH2s )º·9 -(Yco2 + YH2s) 1.6 1 + 15 (YH2s o.s - YH1s
4)
( 1.36)
n
r;cM = L yi ~; - e . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . < I .3 7)
l=l
,,
LY;Pci .T;cM
i=I
p ;cM = -,,-, --------- ,, ( 1.3 8)
LY;~; + YH2s (1- YH2s).e
i•i
donde;
e : Factor de ajuste, ºR
Pe; : Presión crítica del componente i, lpca.
7;; : Temperatura crítica del componente i., ºR
Yco2 : Contenido de C02, fracción molar
y82s : Contenido de H2S, fracción molar
Y; : Componente i en la mezcla, fracción molar
P!cM : Presión seudocrítica de la mezcla corregida por C02 y/o H2S, lpca,
i;~M : Temperatura seudocrítica de la mezcla corregida por C02 y/o H2S, ºR
Wichert y Aziz determinaron que sus ecuaciones presentaron un error absoluto de 0.97%
con un error máximo de 6.59% sobre el siguiente rango de datos: 154 < p(lpca) < 7026 y 40 <
~ºF) < 300, y para concentraciones de C02 < 54.46% (% molar) y H2S < 73.85% (% molar).
El factor de ajuste, e, también puede ser obtenido de la Fig. 1.7.
Luego, teniendo p;cM y T:cM corregidas por C02 y/o H2S, se calcula r; y t; con las
cuales se obtiene z de la Fig. 1.4, 1.5 ó 1.6.
Cuando no se conoce la composición del gas, se deben reemplazar las sumatorias de las
Ecs. ( 1.3 7) y ( 1.3 8) por P"cM y r;cM calculadas en base a la gravedad específica del gas.
Ajustes de las Curvas de Standing, M.B. y Katz, D.L.:
El método gráfico de Standing y Katz ha sido utilizado con mucha confidencia por la
industria del petróleo durante los últimos 40 años, y por ello muchos autores han tratado de
ajustar estas curvas por métodos numéricos con el fin de obtener valores de z por medio del
computador. Entre los métodos o ajustes más conocidos se tienen los siguiente·s:
Carlos Bánzer S.
17
PROPIEDADES FfSICAS DEL GAS
----------------------------------------------------------------·
30
o 10 20 30 •o no •o 10 10
PORCENTAJE DE H2S
Fig. 1. 7.- Factor de ajuste por componentes no-hidrocarburos para la presión
y temperatura seudocríticas de gases naturales (Ref. 12).
Método de Sarem, A.M.14:
Este método para determinar z se basa en los Polinomios de Legendre de grado O a 5. La
ecuación básica de ajuste es:
s s
z = L L Aii P¡ ( x) J'_¡ (y) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ( 1. 3 9)
/•Oj•O
donde;
Carlos Bánzer S.
18
CORRELACIONES NUMÉRICAS P.V.T.
Los Polinomios de Legendre de grado O- 5, !'¡ y P_1 tienen las siguientes formas:
Po(a) = 0.7071068
P¡(a) = l.224745a
_fi(a)=0.7905695 (3a2-1)
~(a)= 0.9354145 (5a3 -3a)
~(a)=0.265165 (35a4 -30a2 +3)
f's(a)=0.293151 (63a5-70a3+15a)
donde a se reemplaza por "x" y "y" al efectuar la sumatoria de la Eo. ( 1.39). La Tabla 1.4
presenta los valores de los coeficientes Ay.
I
•
i= o j= 1 j= 2 j= 3 }=4 i= 5
o 2.1433504 0.0831762 -0.0214670 -0.0008714 0.0042846 -0.0016595
1 0.3312352 -0.1340361 0.0668810 -0.0271743 0.0088512 -0.002152
2 0.1057287 -0.0503937 0.0050925 0.0105513 -0.0073182 0.0026960
3 0.0521840 0.0443121 -0.0193294 0.0058973 0.0015367 -0.0028327
4 0.0197040 -0.0263834 0.019262 -0.0115354 0.0042910 -0.0081303
5 0.0053096 0.0089178 -0.0108948 0.0095594 -0.0060114 0.0031175
Para los intervalos 0.1 s P.\·r s 14. 9 y 1.05 s Y:,, s 2.95, el error del método con respecto a
los valores leídos en las curvas de Standing y Katz fue menor de 0.4o/o. El método de Sarem
puede ser utilizado cuando el gas natural contiene impurezas corrigiendo P.~\:·M y T,. cM por el
método de Wichert y Aziz.
Método de Papay, J.15,16:
La ecuación para el cálculo de z es la siguiente:
= l_ 3.52 P.,·r 0.274 P,\·r-'
Z 098131' + 08 571' ....•....•...•••.••..............•.....• (1.42)
10· sr 10·1 ,\Y'
Takacs16, determinó que el error promedio de este método con respecto a los valores
leídos en las curvas de Standing y Katz fue de -4.873o/o para presiones y temperaturas
seudorreducidas en el rango de: O. 2 s r; s 15. O y 1. 2 :::; 7;., s 3. O. El método de Papay puede
ser utilizado cuando el gas natural contiene impurezas corrigiendo P.'í<:M y 'r,.,M i,or el método de
Wichert y Aziz.
Método de Hall, K.R. y Yarborough, L.11,1&:
Carlos Bánzer S.
19
PROPIEDADES FISICAS DEL GAS
------------------------------------------------------------------
2 3 4
y+y+y-y 2 /) '
-APsr+ 3 -By +Cy =0 (l.44J
(1- y)
donde;
A= 0.06125 t exp[-1.2 (1-t)2]
B= 14.76 t-9.76 t2 +4.58 t3
C=90.7 t-242.2 t2 +42.4 t3
D = 2.18 + 2. 82 t
Takacs16, determinó que el error promedio en el cálculo de z por este método fue de
-0.51 So/o en comparación con el obtenido por el método de Standing y Katz, Para obtener
· buenos resultados con este método, se recomienda que P.,r y Ts, estén entre los siguientes
rangos: 0.1Spsr~24.0 y l.2~7;.r~3.0. Para gases naturales con impurezas se recomienda
corregir Pj·cM y T.-;·cM por el método de Wichert y Aziz.
Debido a que la Ec. ( 1.44) es no lineal, se requiere una solución de ensayo y error para
resolverla. Un método de ensayo y error frecuentemente utilizado, es el método de Newton-
Raphson 19 el cual utiliza el siguiente procedimiento iterativo (aplicado a este caso):
l. Suponer un valor inicial de y1 y calcular F(y1), donde F(y) es el término de la
derecha de la Ec. ( 1.44 ).
=
2. Si F (y1) O se encuentra dentro de una tolerancia especificada ( ± l 0-4· ), se puede:
ó
0.066 1 2 0.32
J
B = (0.62 - 0.23 T_,.,) Psr + t; - 0.86 - 0.037 p; + ] 09(7;,-J) Psr 6
C = 0.132-0.32 logi:r
D = antilog(0.3106-0.49 7;, + 0.1824 7;;,2)
Carlos Bánzer S.
20
CORRELACIONES NUMÉRICAS P.V.T.
La exactitud entre los valores calculados y leídos del gráfico de Standing y Katz fue
verificado para valores de P.,, de O, 1, 2, 3, 4, 5, 7, 9, 11 y 13, y valores de t; de 1.2, 1.3, 1.4,
1.5, 1.7, 2.0 y 2.4, los resultados fueron los siguientes7:
Rango de Psr , l'sr
O.O< Psr < 5.0 O.O<Psr < 13
1. 3 < 7;, < 1. 7 1.2 < Tsr < 2.4
Número de puntos 24 70
Error relativo promedio 0.19% 0.02%
Desviación estandar 0.96% 1.18%
El método no es válido para temperaturas seudorreducidas fuera del rango de 1.2 a 2 . 4.
Para gases naturales con impurezas se recomienda corregir PscM y Tsckt por el método de
Wichert y Aziz.
Método de Dranchuk, P.M., Purvis, R.A. y Robinson, D.e.21:
Este método es el resultado de un ajuste realizado a la ecuación de estado de Benedict,
Webb y Rubin22 (BWR), la cual escribieron en la siguiente forma:
z=l+ exp(-Agp;)
.......... (1.47)
La densidad reducida, p, , se calcula mediante la siguiente expresión (obtenida de la ley
de los gases):
pf 2 l
F=z- l+ 3 exp(-A8p, >J = O
t;
.......... ( 1.49)
Carlos Bánzer S.
21
PROPIEDADES FISICAS DEL GAS
------------------------------------------------------------------·
y
=l+ P, +2
z
2A p2
;
3'
[l + ÁsP~ -(A8p~ )2] exp(-A8p~) (1.50)
z sr
2
r-.«: 1+ P, ·-
2
A1 A8
p~+ A10 (1 + A11p;) Pr3 exp(-A11p;) =O (1.52)
r. +
t;
2
T:r
y
=1+ Pr +2 -
z z
A1 As
T_._,
+
r. 2
Carlos Bánzer S.
22
CORRELACIONES NUMÉRICAS P.V.T.
0.2 Y 1.2 1.05 Y 1.2 Ps,( 1.6643 t; - 2.2114 )- 0.3647 t; + 1.4385 ••••••••••• (1.54)
1.2+ Y 1.4 r; ( 0.0522 t; - 0.85]} )- 0.0364 t; + J .0490 ••••••••••• (} .55)
1.4+ y 2.0 s; ( 0.1391 T:, - 0.2988 )+ 0.0007 7;;, + 0.996'9 ( 1.56)
2.0+ y 3.0 P.,, ( 0.0295 T:, - 0.0825 )+ 0.0009 t; + 0.9967 (1.57)
1.2+ y 2.8 1.05 y 1.2 P.,, (- 1.3570 t; + 1.4942 ) + 4.6315 T:,. - 4. 7009 ( 1.58)
1.2+ y 1.4 P.,r ( 0.1717 T:, - 0.3232 ) + 0.5869 T:,. + 0.1229 (1.59
1.4+ y 2.0 P.,, ( 0.0984 T.,, - 0.2053 ) + 0.0621 T:, + 0.858·0 ( 1.60)
2.0+ y 3.0 P.,, ( 0.0211 T:, - 0.0527 ) + O.O 127 4, + 0.954:9 ··~······· ( 1.61)
2.8+ Y 5.4 1.05 Y 1.2 P.rr (- 0.3278 i:, + 0.4752 ) + 1.8223 T:, - 1.903·6 ••••••••••• ( 1.62)
1.2+ y 1.4 P.,, (- 0.2521 7;,. + 0.3871 ) + 1.6087 T:, - 1.6635 ( 1.63)
J.4+y2.0 p1,.(-0.0284 T,, + 0.0625 )+ 0.4714 i:r- O.OOJI ••••••••••• (1.64)
2.0+ y 3.0 P.,, ( 0.0041 r: + 0.0039 ) + 0.0607 7;;, + O. 792.7 .....•..... ( 1.65)
5.4+ y 15 1.05 y 3.0 Ps,( 0.711 + 3.66 J:.r)-1'4667 -1.637 / (0.319 t; . ._
0.522) + 2.071 ( 1.66)
El método de Gopal puede ser utilizado cuando el gas natural contiene impurezas
corrigiendo PscM y TscM por el método de Wichert y Aziz. .
Ejemplo 1.2:
Utilizando los métodos de Standing, M.B. y Katz, D.L., Sarem, A.M., Papay, J., Hall,
K.R. y Yarborough, L., Brill, J.P. y Beggs, H.D., Dranchuk, P.M., Purvis, R.,A. y Robinson,
D.B., Dranchuk, P.M. y Abou-Kassem, J.H. y de Gopal, V.N., determinar el factor de
compresibilidad z para el gas cuya composición se presenta en la Tabla 1.3 a una presión de
381 O )pea. y temperatura de 194 ºF.
Solución:
Corrección de PscM y TscM por C02 y H2S <P.reM = 709.2 lpca. y T.vcM= 505.8 ºR, Ejemplo 1.1 ).
De la Ec. ( 1.36):
s= 120 [(0.1101+0.0491)º·9 -(0.1101+0.0491)1·6]+ 15 (0.0491°·5 -0.04914)= 19.94 ºR
De las Ecs. ( 1.3 7) y ( 1.38):
7;..'cM = 505.8-19.94 = 485.9 ºR
, _ (709.2)(485.9) _
680·0
PscM - 505.8 + (0.0491)(1-0.0491)(19.94) - lpca.
De las Ecs. (1.34) y (1.35):
3810
Pvr = 680.0 = 5.60
194 + 460
T:r =
485.9
= 135
Carlos Sánzer S.
23
____________________ PROPIEDADES FISICAS DEL GAS ......¡..._ _
Carlos Bánzer S.
24
CORRELACIONES NUMERICAS P.V.T.
luego:
z =0.77343
Método de Papay, J.
De la Ec. (1.42):
_ (3.52)(5.60) (0.274 )(5.60)2
z - 1- l ()(0.9813)(1.35) + l o<0.8157)(1.lS)
z = 0.74746
Método de Hall, K.R. y Yarborough, L.
485.9
1 = 194 + 460 = O. 74
A =(0.06125)(0.743) exp[-1.2 (1-0.743)2]=0.0420
B = (14. 76)(0. 743)-(9. 76)(0. 743)2 + ( 4.58)(0. 743)3 = 7.4573
e= (90. 7)(0. 743)-(242.2)(0. 743)2 + ( 42.4 )(O. 743)3 = -48.9249
D = 2.18 + (2.82)(0. 743) = 4.2753
Sustituyendo estos valores en la Ec, ( 1.44 ), se tiene:
-0.2352 +y+ r. ~ ~ 3- l_.4 - 7.4573y2 -48.9249y42753 = 0
Suponiendo
., un. valor .de y1 = 0.001 para iniciar los cálculos, y sustituyendo este valor en la
ecuacion anterior, se tiene:
F(y1) = -0.23420:;: O (y1 = 0.001 no es la solución)
Derivando la ecuación anterior con respecto a y se tiene:
luego;
dFJ;) = 0.99312
Y1
Carlos Bánzer S.
25
PROPIEDADES FISICAS DEL GAS
Carlos Bánzer S.
26
CORRELACIONES NUMÉRICAS P.V.T.
exp(-(0.68446549)(1.86667)2)
dF(zi) = 8.08878
dz
Calculando un nuevo valor de z 1:
(-0.76596)
Z2 = 0.6 - = 0.69469
8.08878
Nuevo valor de z1 a probar: z1 = 0.69469
La siguiente tabla muestra los tanteos realizados para obtener z:
Pr
0.60000 -0.76596 8.08878 1.86667 0.69469
0.69469 -0.67126 7.1225 l 1.61222 0.78894
0.78894 -0.15844 3.28282 1.41963 0.83720
0.83720 0.05089 2.01272 1.33779 0.81192
0.81192 0.06795 1.72989 1.37945 0.77264
0.77264 0.00869 1.82378 1.44958 0.76812
0.76812 -0.03704 2.25004 1.45810 0.78458
0.78458 -0.02631 2.26722 1.42751 0.79619
0.79619 0.00510 2.10067 1.40670 0.79376
0.79376 0.01499 2.01130 1.411 O 1 0.78630
0.78630 0.00506 2.03973 1.42439 0.78382
0.78382 -0.00535 2.10351 1.42890 0.78637
0.78637 -0.00557 2.12088 1.42427 0.78899
0.78899 -0.00011 +- Se puede considerar como solución
luego;
z= 0.78899
Método de Dranchuk, P.M. y Abou-Kassem, J.H.
Suponiendo un valor de z1 ·= 0.6 para iniciar los cálculos:
De la Ec. ( 1.48):
(0.27)(5.60)
Pr = (0.6)(1.35) = 1.86667
De la Ec. ( 1.52):
1.07 0.5339 o.o 1569 0.05165
F(z1) = 0.6- l + 0.3265- - 3 + 4 - 5 (1.86667) +
1.35 1.35 1.35 1.35
0.7361 0.1844 0.7361 0.1844
0.5475- + 2 (1.86667)2-(0.1056 - + 2 (l.86667'15+
1.35 1.35 1.35 1.35 "
Carlos Bánzer S.
27
PROPIEDADES FfSICAS DEL GAS
dFJzz = 8.25445
1)
Z¡ F(:¡) Pr
0.60000 -0.78279 8.25445 1.86667 0.69483
0.69483 -0.68796 7.26435 1.61190 0.78954
0.78954 -0.16334 3.34348 1.41856 0.83839
0.83839 0.05137 2.03268 1.33589 0.81312
0.81312 0.06966 1.73816 1.37742 0.77304
0.77304 0.00936 1.93708 1.44883 0.76821
0.76821 -0.03830 2.27943 1.45794 0.78501
0.78501 -0.02778 2.29994 1.42674 0.79709
0.79709 0.00501 2.12444 1.40511 0.79473
0.79473 0.01573 2.02857 1.40928 0.78698
0.78698 0.00553 2.05745 1.42316 0.78429
0.78429 -0.00556 2.12568 1.42804 0.78691
0.78691 -0.00600 2.14528 1.42329 0.78970 ''
luego;
z= 0.78970
Método de Gopal, V.N.
De la Ec, (1.66):
-1.4667 1.637
z = (5.60)(0.711 + 3.66 X 1.35) - (OJ l 9)(l.35) + 0.522 + 2.07]
luego;
z=0.79414
Carlos Bánzer S.
28
CORRELACIONES NUMÉRICAS P.V.T.
donde, para la misma masa de gas nR se cancelan, y dado que p1 y T1 son condiciones normales
=
( 14. 7 lpca. y 60 ºF), luego z 1 1.0 y la Ec. ( 1.67) se reduce a:
zp1T (14.7)zT zT PCY
Bg = z1pfa = (1.0)(60 + 460)p = 0.02827 p ' PCN .... ················ (1.68)
Si se desea expresar Bg en BY/PCN, se divide por 5.615:
zT BY
Bg = 0.00503 p , PCN ( 1.69)
donde;
Bg : Factor volumétrico del gas, PCY/PCN o BY/PCN
z : Factor de compresibilidad del gas, adim.
p : Presión, lpca.
T : Temperatura, ºR ( = ºF + 460 )
En algunos casos se utiliza el inverso de Bg (factor de expansión del gas, E8)27:
1 p PCN
Eg = Bg = 35.37 zT ' PCY (l.70)
p PCN
= 198.8 zT , BY (1.71)
z
U Q024r-r-----,------,r------~------r------...-----.......------.. 500 ~
~
m
- ~
a
~
m:E i
,-
::,
.... G)
~
0.012 ......--- ~---r-------1-------+-----~1-------1-~----1200 .. t
a:
e
~
o 500 1000 1500 2000 2500
PRESION, lpca
Fig. 1.8.- Comportamiento típico de 89 y E9 vs. presión a temperatura
constante.
Para calcular Bg o E1 a una determinada presión y temperatura, el valor del factor de
compresibilidad, z, a esas condiciones debe ser conocido. Por lo tanto, si no se dispone de un
valor experimental de z, entonces es necesario recurrir a los métodos que se presentaron
anteriormente para su determinación.
Ejemplo 1.3:
Determinar los factores volumétricos del gas cuya composición se presenta en la Tabla
1.3 a una presión de 3810 lpca. y temperatura de 194 ºF para un factor de desviación de 0.78.
Solución:
De la Ec. ( 1.68):
(0.78)(654) PCY
B1 =0.02827 =0.00379 PCN
3810
De la Ec. (1.70):
3810 PCN
E, =3537 (0.78)(654) =264.17 PCY
Estos factores volumétricos significan que un pie cúbico de gas a condiciones normales
(14.7 lpca. y 60 ºF) ocupará un espacio de 0.00379 pie' en el yacimiento a 3810 lpca. y 194 ºF;
o recíprocamente, un pie cúbico de espacio poroso del yacimiento de gas contiene 264.17 PCN.
Compresibilidad del Gas, e,,:
La compresibilidad isotérmica de un gas, cg en lpc", se define como el cambio
fracciona) en volumen cuando la presión es cambiada a temperatura constante, esto es:
1 8V
e ==--
V ................................................. ( 1 .72)
g
q., 1'
=V! z oz --
1
.......................................... (1.75)
qJ 1' p
donde p en lpca. Así, un gas ideal a 1000 lpca. tendrá una compresibilidad de 1/1000 ó 1000 x
1 o-6 lpc-1• La Fig. 1.9 presenta el comportamiento típico de c8 vs, presión a temperatura
constante.
K.
--
<,
O.O!,
i
1 1
•o 1
11 --- CCJ: p
-; 0.04 .. C4 = .L, .!. (.ll.)
'\\
P Z ¡p T
en
el:
<.!)
..J \ ..
.
~
• Q 0.03
o \
s- \
\
-
~
..J
-
~~
CD 0.02
....
~----- -· . _ ----· ·----- -·------
~~
CI)
L&J
o: ~
CL
~ 0.01 .,
de donde z y su derivada ( & / ép.,·r)1. pueden obtenerse de las Figs. 1.4, 1.5 ó 1.6 y por lo tanto
determinarse e, y desde luego e g • SI'
1 0.27 op, t;
e= --- ................................. (1.80)
r P» z2 T';
1+ P~
Z Opr t;
La derivada ( & / t9p, )1.. se obtiene diferenciando la Ec. ( 1.4 7), la ecuación de estado
desarrollada por Dranchuk, Purvis y Robinson21, esto es:
az 4
ºP, P, + 5AsA6 i~ +
2A,Pr 2 2 2 2
--3- [1 + A8p r - ( A8p,) ] exp( -A8p r) ( 1.81)
r;
Utilizando las Ecs. ( 1.80) y ( 1.81) y la definición de compresibilidad seudorreducida, la
compresibilidad de un gas puede ser determinada siempre que la presión y temperatura del gas
se encuentren dentro de los rangos especificados por I~ ecuación de Dranchuk y Cols.
y
Mattar Cols. utilizando las Ecs. (1.80) y ( 1.81) desarrollaron las Figs. 1.1 O y 1.11. En
estas figuras, el producto de e,.~-, está graficado como una función de la presión y temperatura
seudorreducidas, e; Y T.r,·
Debe notarse, que los puntos mínimos de las isoterrnas del factor z de la Fig, 1.4 son
bastantes pronunciados a valores de temperaturas seudorreducidas por debajo de 1.4. Las
pendientes de estas lineas cambian abruptamente de positivo a negativo, por lo tanto, la Ec.
(1.80) no predice exactamente estas pendientes cerca de estos mínimos.
Así, las Ecs. ( 1.80) y ( 1.81) no deben utilizarse a valores de 7;;, menores de 1.4 y P.,·,
entre 0.4 y 3.0. Esta dificultad fue reconocida en la preparación de las Figs. 1.1 O y 1.11. En los
mínimos, los gradientes ( & I éps, )1. son cero, por lo tanto, debe aplicarse la Ec. (l. 77). Las
isotermas de las Figs. 1.1 O y 1.11 fuéron· ajustadas para pasar a través de estos puntos.
; Carlos Bánzer S.
32
CORRELACIONES NUMÉRICAS P.V.T.
e!
::,
......
~ --··-·
,iiÑN fl
-N-c-.-.-·
N' Ñ - .. _.. ,. •
~
:;?¿
~
~~~ E
::::: ~~~~~
~. _.... ~ -;,,,i ¡;,_.- V
JB
~ V'
-~
¡;;. ~ ¡:;..
~
.J=::::: 'º
-
l.-
·-
-
~
-- i.-- :: L--- V e,
~
--e:
-""'
- - .. "'
o
s- c.e: O')N
~
•
__,. ~ ~~
o Q)
.._ts o:::
~ ~
~
~
l&J
o::
Q:
.._,
o º~~
~
~
o ·-o
:::> QJ
w
en
-o ::s
e -o
z ~Q ~
- - o o o
s:;
r~~
r,-/
~
-~
"'
,, U>
uJ
cr
Q.
·~ "C
>$
m :::,
~
1
l. %a ...
1
•
~
... ~
~
~-~ ~ •
r... •- r>
e- ,, "":.
~ ~
•
·-
LL
C7)
o-
-
o - - od
Jll J:)
~
-
8
~
:J
........
~
Q)
o.
h E
Q)
+J
0~ -o
~
e:
- • ·-
~
o
l\~,oJt..z::: 1"' .... ~~
_..,
--
e,
e(
...
"'c.. O>
e?~
-
_
-f.\)~ ~· ~ ~ rlf!!. -
o·
"~"\)1'~.40 ...:; .-C3
~ ~~
N
~ ~
. u e .
1€ '~ ---,~ .. ~ ~ _.,,,. - - p :) 0-...
o (1)
o
....-.:~ ~
~-- \J
,.,o . ~ . . /
~~ ......
l&J
~
\,,. et:
._e,).._,
o'~
~
o
o Q) ·-
=., -o o::s
,-
,,,, e -e
LaJ
<
z: - z ·-
o
·O e
- o o
~ ~
_.
en -~ "'C
w (U ::,
a:
e, >~
1
o
~
•
•
~
en
~ •
·-
LL
o-
-
o - o-
. -o
o
Carlos Bánzer S.
33
PROPIEDADES FÍSICAS DEL GAS
A Ap ; 0'
-- -- . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ( 1 . 84)
ipsr 7·sr y y2 ipsr T,,
donde;
(] - y)4
=A ........................ (1.85)
l + 4y + 4y2 - 4y3 + y4 - (1 - y)4 (2By - CDy<1J-1>)
_
- --::-------------1_- __A ~--- + CDpg (J)-t>
ip_\·r 1_~, 0.132 1.92 5'
(0.62 - 0.23 T.vr) + ---·- - · - - 0.074 Psr + 1 o 9(.1'sr -1) Psr . exp( B)
T:·r - 0.86
.......... ( 1. 86 )•
. Carlos Bánzer S.
34
CORRELACIONES NUMÉRICAS P.V.T.
Gopal, V.N.
Ecuaciones para
Rango de P.vr Rango de 7;, ips, r."
0.2 y 1.2 1.05 y 1.2 1.6643 t; -2.2114 (1.87)
1.2+ y 1.4 0.0522 r;.., - 0.8511 , ( 1.88)
1.4+ y 2.0 0.1391 T.vr - 0.2988 (1.89)
2.0+ y 3.0 0.0295 T.fr - 0.0825 , (] .90)
1.2+ y 2.8 1 .05 y 1.2 -1.3570 T;..r + 1.4942 (1.91)
1.2+ y 1.4 0.1717 r:., -
0.3232 (1.92)
1.4+ y 2.0 0.0984 T;..r - 0.2053 (1.93)
2.0+ y 3.0 0.0211 T.. r - 0.0527 ······························ ( 1.94)
2.8+ y 5.4 1.05 y 1.2 -O. 3 2 7 8 J:,. + O. 4 7 5 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ( 1 . 9 S)
1.2+ y 1.4 -0.2521 T.. ,+ 0.3871 (1.96)
1.4+ y 2.0 -0. 0284 r:.,
+ O. 0625 ( 1 . 97)
2.0+ y 3.0 0.0041 r:.,
+ 0.0039 , (1.98)
5.4+ y 15 1.05 y 3.0 (0.711 + 3.66 ~",)-1'4667 ••••••••••••••••••••••.•••••••••••••••• (1.99)
Las exactitud de las ecuaciones anteriores no ha sido verificada, sin embargo; se cree que
éstas puedan ser utilizadas dentro de los rangos de presión y temperatura especificados por las
ecuaciones de las cuales fueron obtenidas.
Ejemplo 1.4:
Utilizando el método de Mattar,L., Brar, G.S. y Aziz, K. y las ecuaciones anteriormente
obtenidas, determinar la compresibilidad del gas cuya composición se presenta en la Tabla 1.3.
Solución:
Método de Mattar, L., Brar, G.S. y Azlz, K.
De la Ec. ( 1.48), para z = O. 78970 (Ejemplo 1.2):
p,. = 1.41826
De la Ec. ( 1.81) para 7;, = 1.35 y e; = 5.60 (Ejemplo 1.2):
= 0.58744
8p, 1:,.
De la Ec. (1.80):
Carlos Bánzer S.
35
PROPIEDADES FISICAS DEL GAS
1 027 0.58744
e= ------ ------- = 0.08694
r 5.60 (1.35)(0.78970)2 1.41826
1 + 0:-78970 (0.58744)
Finalmente;
= 0.08690 = 127 S 10-6 l -1
es 680 . x pe
Sarem, A.M.
De las Ecs. (1.40) y (1.41):
x=-0.257
y=-0.684
Polinomios de Legendre de grado O - 5, y argumentos x e y:
P0(x)= O.O P0(y)= O. 7071068
P1(x)= 0.16551 Pi(y) = -0.83773
P2(X) = -0.16474 P2 (y) = 0.31905
P3(X) = -0.25399 P3(y) = 0.42274
P4 (x) = 0.46738 P4 (y)= -0. 89481
P5(x) = 0.04735 P, (y)= o. 79402
Aplicando la Ec. ( 1.82):
Ec. (1.82) = (2.1433504)(0)(0.7071068) + (0.0831762)(0)(-0.83773) +(-0.0214670)(0)(0.31905)
'
+ (-o.ooos114)(0)(0.42274) + (0.0042846)(0)(-0.8948t) + (-o.0016s,s )(0)(0.79402)
+ (0.3312352)(0.16551)(0.7071068) + (-0.1340361)(0.16551)(-0.83773)
+ (0.0668810)(0.16551)(0.31905) + (-0.0271743)(0.16551)(0.42274)
+ (0.0088512)(0.16551)(-0.89481) + (-0.0021521)(0.16551)(0.79402)
+ (0.1057287)(-0.16474)(0.707 l 068) + (-0.0503937)(-0.16474)(-0.83773)
+ (0.0050925)(-0.16474)(0.31905) + (0.0105513)(-0.16474)(0.42274):
+ (-0.0073182)(-0.16474)(-0.89481) + (0.0026960)(-0.16474)(0.794(1)2)
+ (-0.0521840)(-0.25399)(0.7071068) + (0.044312 l)(-0.25399)(-0.83773)
+ (-0.0193294)(-0.25399)(0.31905) + (0.0058973)(-0.25399)(0.4227~)
'
+ (0.00 l 5367)(-0.25399)(-0.89481) + (-0.0028327)(-025399)(0.79402)
+ (O.O 197040)(0.46738)(0.707 l 068) + (-0.0263834 )(0.46738)(-0.83 773)
+ (0.0192621)(0.46738)(0.31905) + (-O.Ol 15354)(0.46738)(0.42274)
+ (0.0042910)(0.46738)(-0.89481) + (-0.0081303)(0.46738)(0.79402)
+ (0.0053096)(0.04735)(0.7071068) + (0.0089178)(0.04735)(-0.83773)
+ (-0.0108948)(0.04735)(0.31905) + (0.0095594)(0.04735)(0.42274)
+ (-0.0060114)(0.04735)(-0.8948 l) + (0.0031175)(0.04735)(0.79402)
luego;
Carlos Bánzer S.
36
CORRELACIONES NUMÉRICAS P.V.T.
ª
itJsr T..'tf'
= 0.06918
= 0.0247
De la Ec. ( 1.84):
a = 0.0420 _ (0.0420)(5.6)(0;02411) =
0_07547
itJsr T.
,\'r
0.2992 (0.2992)
Carlos Bánzer S.
37
PROPIEDADES F(SICAS DEL GAS
0.08266 -6 -1
c8 = = 121.6 X 10 lpc
680
Brlll, J.P. y Beggs, H.D.
De la Ec. (1.86):
ª
qJ
=~---------------------------------------------+
1-0.3245
sr t; 0.132 (1.92)(5.60)5
(0.62 - 0.23 X 1.35) + l.35 _ 0_86 - 0.074 (5.60) + Ú)9(i~3S-ij___ exp(l 1.7837)
= 0.08054
q-}sr T,,.
Carlos Bánzer S.
38
CORRELACIONES NUMÉRICAS P.V.T.
Q.
u
CI) -- 0.04
3
_,
~
o 0.03
s-
CI)
o
u
->
(1) Q02
200ºF
0.01 IOO F
O IOOO 2000 3000 4000 6000 7000
PRESION, lpca
Fig. 1.12.- Comportamiento tfpico de µ8 vs. presión a tres temperaturas
diferentes.
Los siguientes factores afectan la viscosidad de un gas natural:
a. A bajas presiones ( < 1000 - 1500 lpca) a medida que aumenta la temperatura
aumenta la viscosidad del gas debido al incremento de la energía cinética de las
moléculas que producen gran número de choques intermoleculares.
b. A elevadas presiones(> 1000-1500 lpca) a medida que aumenta la temperatura
disminuye la viscosidad del gas debido a la expansión térmica de las moléculas. A
elevadas presiones las distancias intermoleculares de los gases son pequeñas y un gas
tiende a comportarse como un líquido.
c. A cualquier temperatura, la viscosidad de un gas aumenta con el incremento de
presión debido a la disminución de las distancias intermoleculares.
d. A medida que un gas es más pesado, sus moléculas serán más grandes y por tanto su
viscosidad será mayor.
La viscosidad de un gas natural puede ser determinada experimentalmente o por medio
de ecuaciones. La determinación de µ~ en el laboratorio es dificil debido a que su valor es
muy pequeño ( 1::1 O. 2 cp.) para ser medido con exactitud. Por esta razón, se prefiere utilizar
correlaciones gráficas o numéricas para su determinación.
Correlación de Carr, N.L., Kobayashi, R. y Burrows, D.B.3º:
Esta correlación se presenta en la Fig. 1.13, y permite determinar 'la viscosidad del gas a
presión atmosférica y temperatura del yacimiento, µg1, a partir· de su peso molecular o
gravedad específica.
Carlos Bánzer S.
39
PROPIEDADES FISICAS DEL GAS
co,
0..0015 P't"P'l
1~
.......... "l""l""'P'l.......... ~I
.........,~1o++++++o11-Mo+~..........,..., 0.00IOl""t1"1~
o
s •o o
Mol. 4% Nz
20 30 40 so 60 70 80 . 90 100
PESO MOLECULAR
Fig. 1.13.- Viscosidad de gases naturales a 1 atm. y temperatura del ya~imiento (Ref. 30).
Los cuadros insertos en la figura permiten realizar las correcciones por contenido de
H2S, C02 y/o N2• Así, si el gas contiene impurezas, µg1 se corrige de la siguiente manera:
µ g 1 e = µ R 1 + CCO 2 + CH 2 S + CN 2 • • • • • • • • • • • • • • • • • • 1 • • • • • • • • •• • • • • • • • ( } .1 00)
donde;
: Viscosidad del gas a 1 atm. y T, corregida por impurezas, cp.
: Viscosidad del gas a 1 atm. y T, sin corregir (leída de la Fig. 1.13 }, cp.
: Correcciones por presencia de H2S, C02 y/o N2, cp. Estos vale. t.: . se
obtienen de las lineas que aparecen en los cuadros insertos en la Fig. 1.13.
Un ajuste de las curvas de la Fig. 1.13 esta dado por7:
µ,, = (1.709 x 10-s -2.062 x 10-6y ,) T+8.188 x 10-3 -6.15 x 10-311ogy 8
'
•••• (1.101)
El error de ajuste de la Ec. (1.101) a las curvas de la Fig. 1.13 fue de 0.38% con una
desviación estandar de 0.46% para 30 valores en los rangos de: 0.55 < r8 < 1.55 y 100 < 7tºF)
<300.
La viscosidad del gas a la presión requerida se obtiene mediante la determinación del
cociente µ8 I µ11 utilizando las Figs. 1.14 ó 1.1 S. Estas figuras se basan en el Teorema de los
Estados Correspondientes, el cual establece que: a las mismas condiciones de presión y
temperatura seudorreducidas todos los gases naturales tienen el mismo cociente de viscosidad
µg I µ81, siendo µs y µg1 las viscosidades a p y T y a l atm. y Ten cp., respectivamente.
La presión y temperatura seudorreducidas, se determinan en base al procedimiento
anteriormente descrito. En este caso no se debe hacer corrección adicional a Psr y 7s,. por
componentes no hidrocarburos. Finalmente, µ8 está dada por:
µg
.µglc (1.105)
•••••••••••••••••••••••••••••••·•••••••••••••·
µgl
.. 1,
e '
! 4 ~-+--+--+-+-+-~~~~\11-+..+ ........
~
,¡;.
U,t---+--+--+-+-+-~~I-M-M~·..............w.-1....... w... .............~~..._.......___......
1A / f
I
.-·
-
•• ... , ..... . ,.........l. -.~,-"J.t--1
.
.... l .
•.o--
..... r""
'-- ...,__
¡ --
r""
:::. 1
~
- - 1:.111 1
= =:111
1
;:
1.1 t.4
TEMPERATURA SEUDORAEDUCIDA. Ta,
•• ' o a.a
donde;
a0 = -2.46211820 E+ 00 a8 = -7.93385684 E-01
a1 = 2.97054714 E+ 00 a9 = 1.39643306 ·E + 00
a2 = -2.86264054 E-01 a10 = -1.49144925 E- 01
a3 =8.05420522 E - 03 a11 = 4.41015512 E-03
a4 = 2.80860949 E+ 00 a12 = 8.39387178 E - 02
a5 = -3.49803305 E+ 00 a13 =-1.86408848 E-01
a6 = 3.60373020 E - 01 a14 = 2.03367881 E - 02
a1 = -1.04432413 E-02 a1s = -6.09579263 E - 04
( = -0.000609579263 )
El ajuste de Dempsey presenta buenos resultados cuando P» y T.-., están entre los
siguientes intervalos: 1.0 S »; S 20 y 1.2 ~ T.,, s 3.0.
· Carlos Bánzer S.
42
CORRELACIONES NUMÉRICAS P.V.T.
µg
In 7s, = 1.37480
µgl
luego;
µ8 ~, = exp(l.37480) = 3.95429
Pg1
Finalmente;
_ (3.95429)( O.O 1148) _ O
µg - - .03363 cp.
. 1.29 ·
Correlación de Lee, A.L., Gonzélez, M.H. y Eakin, B.E.
Del Ejemplo 1.2, para z = O. 78 (Método de Standing, M.B. y Katz, D.L.):
-3 (3810)(35.969)
P« = l.4935 x 10 (0.78)(654) =0.40122
= (9.4+0.02 X 35.969)(654)1.S =
109.44492
K 209 + (19)(35.969) + 654
· 986
X=3.5+ +(0.01)(35.969)=5.36734
654
Y= 2.4-(0.2)(5.37) = 1.326
De la Ec. ( 1.107), la viscosidad del gas a p y T:
109.44492 exp(S.36734 x 0.401221•326)
µg = 104
Finalmente;
µg = 0.05475 cp.
Un resumen de los resultados de µ1 es como sigue:
Correlación µ9 , cp.
Carlos Bánzer S.
44
CORRELACIONES NUMÉRICAS P.V.T.
Carlos Bánzer S.
45
PROPIEDADES FISICAS DEL GAS
3. Kay, W.B.: ''Density of Hydrocarbon Gases and Vapors at High Temperature and
Pressure," Ind. Eng. Chem. (Sept. 1936) 28, 1014-1019.
4. Stewart, W.F., Burkhardt, S.F., and Voo, D.: "Predictíon of Pseudocritical Parameters
for Mixtures," paper presented at the AIChE Meeting, Kansas City, MO, May 18,
1959.
5. Sutton, R.P.: ''Compressibility Factors for High-Molecular-Weigbt Reservoir Gases,''
paper SPE 14265 presented at the SPE 60th Annual Technical Conference and
Exhibition, Las Vegas, NV, Sept. 22-25, 1985.
6. Brown, G.G., Katz, D.L., Oberfell, G.G., and Alden, R.C.: Natural Gasoline and the
Volatile Hydrocarbons, Natural Gasoline Ass. of America, Tulsa (11948), 24-32.
7. Standing, M.B.: Vo/umetric and Phase Behavior o/ Oil Fie/d Hydrocarbon Systems,
SPE, Dallas ( 1977).
8. Mathews, T.A., Roland, C.H., and Katz, D.L.: ''High Pressure Gas Measurement,"
Proc., NGAA (1942) 41.
9. Kessler, M.G. and Lee, 8.1.: ''Improve Prediction of Enthalpy of Fractions," Hyd
Proc. (March 1976) 55, 153-158.
10. Whitson, C.H.: ''Effect of Physical Properties Estimation on Equation-of-State
Predictions,'' paper SPE 11200 presented at the 57th Annual Fall Technical
Conference and Exhibition, New Orleans, LA, Sept. 26-29, 1982. :
11. Standing, M.B. and Katz, D.L.: ''Density of Natural Gases," Trans., AIME (1942)
146, 140-149.
12. Wichert, E. and Aziz, K.: ''Compressibility Factor of Sour Natural Gases," Can. J.
Chem. Eng. (Apr. 1971) 49, 267-273.
13. Wichert, E. and Aziz, K.: ''Calculate Z's for Sour Gases,'' Hyd. Proc. (May 1972) 51, ·
119-122.
14. Sarem, A.M.: ''Z-Factor Equation Developed for Use in Digital Computers," Oil and
Gas J. (Sept. 18, 1961) 118.
15. Papay, J.: ''A Termeléstchnológiai Paraméterek Változása a Gástelepek Müvelése
Során," OGIL Musz, Tud. Kuzl., Budapest (1968) 267-273.
16. Takacs, G.: ''Comparisons Made for Computer Z-Factor Calculations," Oil and Gas
J. (Dec. 20, 1976) 64-66. .
•
17. Hall, K.R. and Yarborough, L.: ''A New Equation ofState for Z-F~ctor Calculations,"
Oil and Gas J. (June 18, 1973) 82-85, 90-92. '
18. Yarborough, L. and Hall, K.R.: ''How to Solve Equation of State for Z-Factors," Oil
andGasJ. (Feb. 18, 1974) 86-88.
19. Camahan, B., Luther, H.A., and Wilkes, J.O.: Applied Numerieal Methods, John
Wiley & Sons Inc., New York City (1969).
20. Brill, J.P. and Beggs, H.D., University of Tulsa: ''Two-Phasq Flow in Pipes,"
INTERCOMP Course, The Hague, 1974.
21. Dranchuk, P.M., Purvis, R.A., and Robinson, D.B.: "Computer Calculations of
Natural Gas Compressibility Factors Using the Standing and Katz Correlation,"
Institute of Petroleum Technical Series, No. IP 74-008 (1974) 1-13~
Carlos Bánzer S.
46
CORRELACIONES NUMÉRICAS P.V.T.
22. Benedict, M., Webb, G.B., and Rubio, L.C.: ''An Empirical Equation for
Thermodynamics Properties of Light Hydrocarbons and Their Mixtures, l. Methane,
Ethane, Propane and n-Butane," J. Chem. Phys. (April 1940) 8, 334-345.
23. Katz, D.L., Cornell, D., Kobayashi, R., Poettman, F.H., Vary, J.A., Elenbaas, J.R.,
and Weinaug, C.F.: Handbook o/ Natural Gas Engineertng, McGraw-Hill Book Co.,
Inc., New York City ( 1959).
24. Dranchuk, P.M. and Abou-Kassem, J.H.: ''Calculation of Z-Factors for Natural Gases
Using Equations of State," J. Cdn. Pet. Tech. (July-Sept., 1975) 34-36.
25. Starling, K.E.: Fluid Thermodynamic Properties for Light Petroleum Systems, Gulf
Publishing Co., Houston ( 1973 ).
26. Gopal, V.N.: "Gas Z-Factor Equations Developed for Computer," Oil and Gas J.
(Aug. 8, 1977) 58-60.
27. Dake, L.P.: Fundamentals of Reservoir Engineering, Elsevier, Amsterdam, 1978.
28. Trube, A.S.: ''Compressibility ofNatural Gases," Trans., AIME (1957) 210, 355-357.
29. Mattar, L., Brar, G.S., and Aziz, K.: ''Compressibility of Natural Gases," J. Cdn. Pet.
Tech. (Oct-Dec., 1975) 77-80.
30. Carr, N.L., Kobayashi, R., and Burrows, D.B.: "Viscosity of Hydrocarbon Gases
Under Pressure," Trans., AIME (1954) 201, 264-272.
31. Dempsey, J.R.: ''Computer Routine Treats Gas Viscosity as a Variable," Oíl and Gas
J. (Aug. 16, 1965) 141.
32. Lee, A.L., González, M.H., and Eakin, B.E.: ''The Viscosity of Natural Gases," J.
Pet. Tech. (Aug. 1966), 997-1000.
carios Bánzer S.
47
PROPIEDADES FISICAS DEL PETRÓLEO
En esta sección se presenta una revisión de las propiedades fisicas del petróleo (petróleos .
negros) y las correlaciones que se han desarrollado para su determinación,
Las correlaciones que se presentarán en esta sección han sido desarrolladas utilizando ·
sistemas de crudos de varias regiones productoras. Estos crudos exhiben una tendencia
regional en su composición química que los categoriza como: parofinicos, nofténicos o
, .
aromáticos.
Debido a estas diferencias en composición, las correlaciones desarrolladas a partir de
muestras regionales de una base química predominante pueden generar resultados erróneos
cuando se aplican a crudos de otras regiones. Por lo tanto, el uso efectivo de las correlaciones
radica en el entendimiento de su desarrollo y en el conocimiento de sus limitaciones.
Dado que las correlaciones que se presentarán en esta sección, utilizan Ja gravedad
específica, y0, o la gravedad API del petróleo, yAPJ, es necesario definir estas propiedades.
Gravedad Específica del Petróleo, rO :
La gravedad específica del petróleo, y"' se define como la razón de la densidad del
petróleo a la densidad del agua, ambas medidas a la misma presión y temperatura":
r o = Pº (2.1)
Pw
La gravedad específica (adimensional) normalmente es reportada como 6(')º/60°, esto
significa que las densidades del petróleo y del agua fueron medidas a 60 ºF~ y presión
atmosférica.
Para el petróleo es de uso común la gravedad API (American Petroleum lnstitute),
definida mediante la siguiente relación:
141.S
Y API = - }31.5 (2.2)
r ,,
Donde, r APJ , es la gravedad API del petróleo, º API, y y" la gravedad específica del
petróleo a 60°/60°. Obsérvese que petróleos más livianos que el agua y < 1.0, tienen r APJ
0
mayor de 1 O, y más pesados que el agua tienen y APJ menor de 1 O. Al agua, r w == 1.0, le
corresponde una y APJ de 1 O.
Presión de Burbujeo, pb:
La presión de burbujeo o presión del punto de burbujeo como también se le llama, se
define como la presión a la cual se forma la primera burbuja de gas al pasar un sistema del
estado líquido al estado de dos fases, donde la fase líquida está en equilibrio con una cantidad
infinitesimal de gas libre.
A presiones por debajo de la presión de burbujeo se forman dos fases en equilibrio:
líquida (crudo con gas en solución) y gaseosa (gas natural). La fase líquida está constituída por
crudo saturado con gas natural, por esto a estas presiones se les llama de saturación, Ps·
Comúnmente, la presión de saturación.p; se utiliza como sinónimo de presión de burbujeo.jz,
Carlos Bánzer S.
*Referencias al final de la sección 48
CORRELACIONES NUMÉRICAS P.V.T.
6000 p~
,~ ~
1
5000 f
-vft t
400o • •
'
1
.
~ 3000 1
~
..
O 2000
..,
w
::,
o'
;a ~
m
o:
4~ J
:::, i4
)
m 1000 . !
w
e e00
z
ºffi
~
V
00 ... , ..
-z
6
50 o '
o:
a. 4oc -;1- '
,..
•
300
~-
~
..,
'
200 / -
V
100
10 20 30 40 50 60 80 200 300 400 'ºº 600 1000
FACTOR DE CORRELACION , F
Fig. 2.1.- Correlación de Standing para la presión de burbujeo (Ref. 2).
La siguiente ecuación es un ajuste de la correlación gráfica de la Fig. 2.13:
Ps = 18.2 [F - 1.4] (2.3.)
donde;
0.83
F= Rsh 10(0.0009J 7'-0.0J 25y .IP/) 0 •••••• •••••••••• • • • • •. •••• • • • • ••••••• • • (2.4)
Yg
Pb : Presión de burbujeo, lpca.
R.th : Razón gas disuelto-petróleo a p ~ PCN/BN r..
T : Temperatura del yacimiento, ºF
y8 : Gravedad específica del gas, ( aire = 1 )
r APJ : Gravedad del petróleo, º API
Además, para facilitar los cálculos la ecuación fue llevada a un nomograma el cual se
presenta en la Fig. 2.2.
Carlos Bánzer S.
49
PROPIEDADES F(SICAS DEL PETRÓLEO
EJEMPLO
DETERMINAR:
u presi6n ele ._.... a JOO •f • un pet,MIO dt JO
• API con 11111 ,uón P• dhuello,peb"IO de JSO PCNIIN, y
11t1 wavtc11d ttptcHka del ps de 0.75.
PROCEDIMIENTO: ~"+
Tal CGIDO lo indica ta llecha, ¡ putlt dll valor dt 350 se ~
tityf horizonc1l1111te a la lfnu de 0.75 de ••cltd del ~'
pt; vtrtic:llmentt 111 linea dt 30 • API, hori101Wwt1 1 ~ ·
la linea de 200 •f y vtllkalmente se lee ... valor • pre1ldft ~
• b ..... • ltJO lpa. CJ'l, ·
s.
~.,.
et.~
~
Cario• Bánzer S.
50
CORRELACIONES NUMÉRICAS P.V.T.
' ta
• 4.8
~
)C
.._ 4.4
.Q
e,
-o
1 1
4.0
..,::::,
UJ
m 3.6
a:
:::,
m 3.2
UJ
Q
z
-
o
,n
2.8
~
D. 2.4
=5
w
Q 2.0
a::
o
P-
1.6
u.
~
1.2
Debido a que el peso molecular efectivo del petróleo, ~' es una cantidad desconocida,
Lasater relacionó el peso molecular efectivo a la gravedad API del petróleo. Esta relación se
presenta en la Fig. 2.4.
1: 5S
fta ~
' '~
200 300 400 500
L--&.---"'-"--1~--l----+---~--+----f---ri-~,.-t,,"7"'1.000
li---+--~--+--+--+--tGAA~ D1L GAi 1,...-l.....-'.....,.q.+,,,,., _,._,,. 5000
.Alll• 1.0
e
~--"----l----+----+---+-~t-;-~~~,,Y--t--11000 ~
a
CIMWDADOIL NihOl.aO O
-l---·~M~~~·•J:...._. .....L...¡_ ~.,.q.~~~Cifr~+---t---12000 W
m
>--4--1--H'f~f--t---+-· -1--...C~-,,c..i~+,1,.....,._-t--+t--+---tr--t lOOO ;
en
f
100
~---+-~~""Hlll---i--...,_-+-i.,&-,~~~~+-~~-t+--t---'1-4IOO
A--~~u..-~:.uu~---l."-~·~;..,.q.~"4-----+--+---t---M----t----r--t400
§
'9 IIO
~'~~+-~~'-,L,l,~'-+-+4--4---+--+--+----I--~~~~
;~~~~~~~---h'--+---+--+--+--+--+--7'~~~
flEMPLO
DETERMINAR:
•A,I
La PNll6n de
CCN1
a 200 •f
111 ,, ualto r,nM10
* IIR p:11••
1111
100
• JO
PCN$N, y
tllll pawedad ..., .. nea ....... o.to.
PROCEDIMIENTO:
Tal como lo illdlca la •ectaa. a partk dea valor • seo •
.... hDJIIOfttzl ¡¿ •••• 1 11 ....... 10 ·APt1 "'1kal ,11nte
a II cuna• co,11lacNIII htalaom e Ir 1nte a la ffw ~00
••11ca
clll ps , "'''°
.,, V'lftlclb, 1 nlt .... ,.,......•• , ... d ,,
,_....,,,1111111•
ele bultu.,eo • 2W
•
IN ... Yllor efe ptl•l6n
fpa ..
Carlos Bánzer S.
53
PROPIEDADES FISICAS DEL PETRÓLEO
Las siguientes ecuaciones desarrolladas por Beggs5 pueden ser utilizadas en lugar de las
Figs. 2.3 y 2.4.
Flg. 2.3:
Para y 8 ~ 0.60:
PbY
T
g
= 0.679 exp(2.786y 8)
- 0.323 (2.6)
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (2 .1 O)
donde;
Ps : Presión de burbujeo, lpca.
Rsh : Razón gas disuelto-petróleo a p ~ ph, PCN/BN
T : Temperatura del yacimiento, ºF
r ge : Gravedad específica del gas corregida a 100 lpcm., ( aire = 1)
r API : Gravedad del petróleo, º API
carios Bánzer s.
54
CORRELACIONES NUMÉRICAS P.V.T.
Debido a que la gravedad del gas depende de las condiciones bajo las cuales el gas es
separado del petróleo, Vázquez y Beggs desarrollaron una correlación para normalizar la
gravedad del gas a una presión de separación de 100 lpcm. ( 114.7 lpca), Esta presión fue
tomada como representativa de las condiciones promedias de separación en el campo. La
gravedad específica del gas obtenida de la correlación debe ser utilizada en las correlaciones
desarrolladas por Vázquez y Beggs. La correlación es· ta siguiente:
donde;
y ge : Gravedad específica del gas corregida a l 00 lpcm., (aire= l)
y8 : Gravedad específica del gas sin corregir, (aire= 1)
r AJ'J : Gravedad del petróleo, º API
Psp : Presión del separador, lpca.
Tsp Temperatura del separador, ºF
:
donde;
Carlos Bénzer S.
55
PROPIEDADES FISICAS DEL PETRÓLEO
osre To.112
'
............................................... (2 .13)
r API 0.989
Un total de 45 muestras de crudos la mayoría de la región del Mar del Norte fueron
utilizadas en el desarrollo de esta correlación. La Tabla 2.6 presenta el rango de los datos
utilizados en su desarrollo. Glaso reporta un error promedio de 1.28% con una desviación
estandar de 6.98% para la correlación.
e H2S = 1.0 - (0.903 5 + 0.00 l 5y API ) YH¡S + o.o 19 (45 - r Al'/) YH2S
2
(2.16)
donde, Ten ºF y YN , Yco y YH s las fracciones molares de N2, C02 y H2S en los gases de
supe rflci
cíe, l 2 2
Carlos Bánzer S.
56
CORRELACIONES NUMÉRICAS P.V.T.
donde;
Ps : Presión de burbujeo, lpca.
R_..h : Razón gas disuelto-petróleo a p ~ p., PCN/BN
T : Temperatura del yacimiento, ºF
yK : Gravedad específica de I gas, ( aire = 1 )
r APJ : Gravedad del petróleo,
0
API
Esta correlación fue obtenida a partir de 336 análisis PVT de crudos del Oriente de
Venezuela. El análisis estadístico mostró que el 86.5% de 289 valores de Pb calculados con la
correlación, presentaron un error menor del 20% en comparación c:on los valores
experimentales. La Tabla 2.7 presenta los valores de las constantes C1, C2, C3 y C4 de acuerdo
a la gravedad API del crudo.
C1
Carlos Sánzer S.
58
CORRELACIONES NUMÉRICAS P.V.T.
P» = ·········································· (~-~:3)
donde;
Pb : Presión de burbujeo, lpca.
s; : Razón gas disuelto-petróleo a p ~ Pb, PCN/BN
T : Temperatura del yacimiento, ºF
: Gravedad específica del gas corregida a 100 lpcm., (aire= 1)
r se : Gravedad del petróleo, 0 API
Y API
Un total de 5392 puntos de· datos obtenidos de 740 muestras diferentes de crudos del
Sud-Este de Asia (Indonesia), Norte América, Medio Oriente y América Latina fueron
Carlos Bánzer S.
59
PROPIEDADES FISICAS DEL PETRÓLEO
utilizados en el desarrollo de esta correlación. Los datos fueron separados en dos grupos
debido a la volatilidad de los crudos. El primer grupo contenía crudos con gravedades s 30
O
API, y el segundo crudos con gravedades> 30 °API. La Tabla 2.11 presenta los valores de las
constantes C1, C2, C3 y C4 de acuerdo a la gravedad API del crudo. La Tabla 2.12 presenta el
rango de los datos utilizados en el desarrollo de la correlación, Kartoatmodjo y Schmidt
reportan un error promedio relativo de 3.34032% con una desviación estandar de 26:3.08% para
la correlación. ·
•
Y API S 30 0.05958 0.7972 13.1405 0.9986
Y API > 30 0.03150 0.7587 11.2895 0.9143
r ec = y 1. 1 + 0.1595y
~·
.,,º· º ª T:.n -o.
;u-
4 7
»r
2466 log P.,p
114 • 7
(2.24)
donde;
y,,. : Gravedad específica del gas corregida a 100 lpcm., (aire= 1)
r, : Gravedad específica del0 gas sin corregir, (aire= 1)
r APJ : Gravedad del petróleo, API
Psp : Presión del separador, lpca.
Tsp : Temperatura del separador, ºF
Cario• Bánzer S.
60
CORRELACIONES NUMÉRICAS P.V.T.
Ejemplo 2.1 :
En una prueba de producción, un pozo produjo 700 BID con un estrangulador de 1/8'' y
una razón gas-petróleo de 675 PCN/BN con un separador operando a 100 lpca. y 85 ºF. La
gravedad del petróleo es 31 °API, la gravedad específica del gas producido del separador es
0.95 y la temperatura de fondo, determinada de registros eléctricos es 180 ºF. El gas liberado
del separador contenía 20% (% molar) de C02 y 10% (% molar) de H2S.
Existe una pequefia capa de gas pero el pozo está produciendo de la zona de petróleo.
Determinar la presión de saturación del petróleo utilizando las correlaciones de Standing, M.B.,
Lasater, J.A., Vázquez, M.E. y Beggs, H.D., Glaso, O., TOTAL. C.F.P., Al-Marhoun, M.A.,
Dokla, M.E. y Osman, M.E., Petrosky, G.E., Jr. y Farshad, F.F., y de Kartoatmodjo, T. y
Schmidt, Z.
Solución:
Correcciones por C02 y H2S
De las Ecs. (2.15) y (2.16):
Cco2 = l.O - ( 693.8)(0.20Xl 80)-1.ss3 = 0.956
CH2S = 1.0 -[0.9035 + (0.0015X31)](0.1 O)+ (0.019)( 45 - 31)(0.10)2 = 0.902
Correlación de Standing, M.B.
De las Ecs. (2.3) y (2.4 ):
675 0.83
F= 1 o<o.00091 x iso-o.eus x 31> = 139_03097
0.95
P.,. = 18.2 [139.03097 - 1.4]
P.,· = 2505 lpca.
Corrección de Ps por C02 y H2S:
P«: = (2505)(0.956)(0.902) = 2160 lpca.
Correlación de Lasater, J.A.
De la Fig. 2.4 6 Ec. (2.8):
M0 = 630 - (10)(31) = 320
De la Ec. (2.5):
Yg = 0.651 > 0.60
De la Fig. 2.3 ó Ec. (2. 7):
PY g = (826)(0.651)3's6 + l.95 = 3.742
T
Finalmente;
. = (3.742) 180 + 460
Ps 0.95
P.-. = 2521 lpca.
P« = (2521)(0.956)(0.902) = 2174 lpca.
Carlos Bánzer S.
61
PROPIEDADES FISICAS DEL PETRÓLEO
1.187
675
P.f =
23.931 X 3)
(O.O 178)(0.941) exp --- ···-·-
180 + 460
Ps = 2855 lpca.
P« = (2855)(0.956)(0.902) = 2462 lpca,
Correlación de Glaso, O.
De las Ecs. (2.12) y (2.13):
675 0.816 100.172
F= =17.37159
0.95 31 0989
·
P.,· =} o(l.7669+1.7447 log(l7.371S9)-0.30218 (log(l?.37159))2)
Carlos .Sánzer S.
1
62
CORRELACIONES NUMÉRICAS P.V.T.
ps = 2060 lpca.
P«: = (2060)(0.956)(0.902) = 1776 lpca.
Correlación de Petrosky, G.E., Jr. Farshad, F.F.
De las Ecs. (2.21) y (2.22):
0.S774
F- 675 10(4.561 x 10-' x 1801·1911-7.916x 10-' x3J1·541º) _
- 0.950.8439 - 36.11787
P, = 112.727[36.11787 -12.34]
P . . = 2680 lpca.
P« = (2680)(0.956)(0.902) = 2311 lpca.
Correlación de Kartoatmodjo, T. y Schmidt, Z.
De la Ec. (2.24):
1
r ge = (0.95) 1 + (0.1595)(31)º·4º73 (85)-0.2466 log OO = 0.938 (aire= 1)
114.7
De la Ec. (2.23):
De la Tabla 2.11: C1 = 0.0315, C2 = 0.7587, C3 = 11 ~2895, C4 = 0.9143
0.9143
675
P.¡=
l l .289S X 31 )
(
(0.03) 5)(0.938)º·75B7 ¡ 0 IS0+4óO
Se dice que un crudo está saturado con gas cuando tiene en solución la máxima cantidad
de gas que el admite, y en este caso, al reducir un poco la presión ocurre liberación de gas. Si
el crudo tiene una cantidad de gas inferior a la que puede admitir a ciertas condiciones de p y T
se dice que el crudo está subsaturado (no saturado), o sea que tiene una deficiencia de gas.
En un crudo, si se dispone de suficiente gas y se puede someter el petróleo a altas
presiones y temperaturas, puede decirse que el gas es infinitamente soluble en el petróleo. Lo
que limita la solubilidad del gas en el petróleo en el yacimiento, son las condiciones de presión
y temperatura que allí existan y por supuesto la cantidad de gas disponible. Por esta razón, se
encuentran yacimientos con capa de gas inicial (exceso de gas a las condiciones de p y T
existentes) y yacimientos subsaturados (no saturados).
La Fig. 2.6 ilustra el comportamiento típico de R, vs. presión a temperatura constante.
z
CD
600--------------.....------------~---------------...-----~--
S67 PCN¡ IN
'~ 500.__ __ _.___ __. ....___ --+- __
2 3
t:; A 1200 1 pea
e
a.
~----....,_ Rs=337
~----,.__ .._ e
~
.._ l
~- 200t----
i
..J
f
~ ~
:
~ ;.
o .... ..... ..__ _,_ _._
~
._. ~
O 1000 1!500 ~ 2SOO !SOO
PRESION, lpca
Fig. 2.6.- Comportamiento tipico de R, vs. presión a temperatura
constante.
Obsérvese que R, aumenta con presión hasta llegar a la presión de burbujeo, Pb, a partir
de la cual se mantiene constante. La razón de esto, es que al llegar a) punto de burbujeo no
existe más gas disponible para entrar en solución con el petróleo·. Así, por arriba de la presión
de burbujeo el petróleo está subsaturado ya que una disminución de presión no causa
liberación de gas iR, es constante); sin embargo, por debajo de la presión de burbujeo el
petróleo está saturado, ya que una disminución de presión origina una liberación de gas debido
a que R, disminuye.
Los siguientes factores afectan la solubilidad del gas en el petróleo:
Presión: Al aumentar la presión aumenta R;
Temperatura: Al aumentar la temperatura disminuye R;
Gravedad del petróleo: Al aumentar la gravedad API aumenta R;
Gravedad del gas: Al aumentar la gravedad específica del gas aumenta R;
Tipo de liberación: La manera como se separa el gas del petróleo produce diferentes
valores de R; Existen dos tipos de liberación: a) Instantánea y b) diferenclal 1•
Carlos Bánzer S.
64
CORRELACIONES NUMÉRICAS P.V.T.
La exactitud de los resultados obtenidos con la correlación presentaron un error del 10%
respecto a los valores determinados experimentalmente. Si se desea se puede utilizar el
nomograma presentado en la Fig, 2.2.
Correlación de Lasater, J.A.4:
Esta correlación se presenta en las Figs. 2.3 y 2.4. Si se desea, se puede utilizar el
nomograma presentado en la Fig. 2.5. El procedimiento para obtener Rsa partir de las Figs. 2.3
y 2.4 es el siguiente:
1. Calcular el ''factor de la presión de burbujeo", p 1 = (Phr s I 1''), donde Pb es
cualquier presión de saturación.
2. De la Fig. 2.3 obtener Yg·
3. De la Fig. 2.4 obtener M0 a partir de la gravedad API del petróleo.
4. Determinar R, mediante la siguiente ecuación:
132755y ,)Y g
· R.f = (2.26)
M0(l-yg)
Las ecuaciones anteriormente presentadas por Beggs', también pueden ser utilizadas
para determinar R, mediante la Ec. (2.26). Así, M0 se obtiene de las Ecs. (2.8) o (2.9), y y8 de
las Ecs. (2.6) o (2. 7) dependiendo de si el valor del factor de la presión es menor o mayor que
3.29, esto es:
Para PY 8 < 3.29:
T
}.473py g
Yg =0.359
In ---+0.476 ..................................... (2.2'7)
T
Carlos Bánzer S.
65
PROPIEDADES FISICAS DEL PETRÓLEO
py
Para g ~ 3 .29:
T
0.281
0.121py g
y8 = -0.236 ........................................ (2.~~)
T
donde Ten ºR.
Correlación de Vázquez, M.E. y Beggs, H.D.1:
La correlación es la siguiente:
e C3y API
R, = C1 y ge p ·2 exp (2.29)
T+460
donde;
R, : Razón gas disuelto-petróleo a p < Pb , PCN/BN
p : Presión de interés, lpca.
T : Temperatura del yacimiento, ºF
r ge : Gravedad específica del gas corregida a 100 lpcm., (aire= 1)
r APJ : Gravedad del petróleo, 0 API
R, = Y g • F r AP1
To.112 ··········································· (~.3())
donde;
F = l 012.8869-(14.1811-3.3093 1011P>º·'1 (2.31)
R, : Razón gas disuelto-petróleo a p < Pb , PCN/BN
p : Presión de interés, lpca.
T : Temperatura del yacimiento, ºF
yg : Gravedad específica del gas, (aire= 1)
r APJ : Gravedad del petróleo, 0 API
donde;
R, : Razón gas disuelto-petróleo a p < Ph , PCN/BN
p : Presión de interés, lpca,
T : Temperatura del yacimiento, ºF
y8 : Gravedad específica del gas, (aire= 1)
r APJ : Gravedad del petróleo, 0 API
¡Carlos ,Sánzer S.
66
CORRELACIONES NUMÉRICAS P.V.T.
La Tabla 2.14 presenta los valores de las constantes C1, C2, C3 y C4 de acuerdo a la
gravedad API del crudo.
R., = [ 185 .84321 PY g i.s11s4 y" -3.t437 r-1.326s1 r.3984 ················ (2.33)
donde;
R, : Razón gas disuelto-petróleo a p < Ph , PCN/BN
p : Presión de interés, lpca.
T : Temperatura del yacimiento, ºR
rg : Gravedad específica del gas, (aire= 1)
r O
: Gravedad específica del petróleo, (agua= l)
Correlación de Dokla, M.E. y Osman, M.E.12:
La correlación es la siguiente:
R. . = [ 0.11956 x 10-3 rr s 1.01049 r "-0.101991 Tº·9s2ss4 ]
l.3811
..................... (2.34)
donde;
R, : Razón gas disuelto-petróleo a p < Pb , PCN/BN
p : Presión de interés, lpca.
T : Temperatura del yacimiento, ºR
yg : Gravedad específica del gas, (aire= 1)
r" : Gravedad específica del petróleo, (agua= 1)
Correlación de Petrosky, G.E., Jr. y Farshad, F.F.13:
La correlación es la siguiente:
- 1.73184
(7.916 10-"r.-,,,,1.s.aio_4_561 10,-57'1·''111)
R.,· = y V 0.8439 p + 12.34 10 X X
..... (2.3 S)
~ 112.727
donde;
Rs : Razón gas disuelto-petróleo a p < Ps , PCN/BN
p : Presión de interés, lpca.
T : Temperatura del yacimiento, ºF
Yg : Gravedad específica del gas, (aire= 1)
y AJJJ : Gravedad del petróleo, API
0
Petrosky y Farshad reportan un error promedio relativo de -0.0So/o con una desviación
estandar de 4. 79% para la correlación. La Tabla 2.1 O presenta el rango de los datos utilizados
en su desarrollo.
Carlos Bánzer S.
67
PROPIEDADES FiSICAS DEL PETRÓLEO
s, = s; = 560 PCN/BN
Para p = 2000 lpca. (p < p,, ).
F = l o(2.ss69-(t4.1s11-J.3093 to1<2000>>º·'1 = 12.08324
1.2255
310.989
R.,· =(0.95) (12.08324) --
1800.172
s .,b • 15.0057
Rs = 615 PCN/BN
Correlación de Dokla, M.E. y Osman, M.E.
De la Ec. (2.34 ):
Para p = 4000 lpca. (p > Pb ), entonces R.-; = R.,·h
1.3811
R.,· = R.fb = [ co.11956 x 10-3 )(2500)(0.95) 1.01049 co.s11)-0.1o,991 < 640) o.9s2ss4 ]
.
iCsrlos Bánzer S.
70
CORRELACIONES NUMÉRICAS P.V.T.
Rs, PCN/BN
Carlos Sánzer S.
71
.
PROPIEDADES FISICAS DEL PETRÓLEO
lo A 1200
lpca = 1.210 o
zo
~
--
.;..,
l~:::, e
..... -u
z
-
!!
zm
o...,
!-
en
a;c
l&I
~
-"'f
L
PRESION, lpca
Fig. 2.7.- Comportamiento típico de 80 vs. presión a temperatura
constante.
Debido a que no se libera gas de solución cuando la presión disminuye desde su valor
inicial de 3500 lpca. hasta la presión al punto de burbujeo a 2500 lpca., el fluido de yacimiento
permanece en estado monofásico (petróleo subsaturado), sin embargo; como los liquidas son
ligeramente compresibles, el volumen aumenta de 1.31 O BY/BN a la presión inicial a 1.333
BY /BN a la presión del punto de burbujeo.
Por debajo de la presión de burbujeo, la expansión líquida continúa, pero su efecto sobre
el aumento en volumen es insignificante comparado con un efecto mucho mayor: la reducción
en el volumen líquido debido a la liberación del gas en solución. A 1200 lpca. el volumen
disminuye a 1.210 BY/BN, y a la presión atmosférica y 160 ºF disminuye a ,J.040 BY/BN. La
disminución de 1.040 BY/BN a 1.000 BY/BN a presión atmosférica y 60 ºF, se debe a la
reducción o me1111a por efectos de temperatura, ya que la presión se mantuvo.constante. Siendo
B O proporcional a Rs, los factores que afectan a R." afectan a B" de la misma manera.
Las siguientes correlaciones permiten determinar B,, a la presión del punto de burbujeo
(p = Pb).
Correlación de Standing, M.B.2,3:
Standing encontró una correlación entre el denominado factor de correlación, F, y el
factor volumétrico del petróleo a la presión del punto de burbujeo, B"h· La Fig. 2.8 presenta la
correlación de estos factores.
Carlos Bánzer S.
72
CORRELACIONES NUMÉRICAS P.V.T.
...o
~"
"
'/
-
~
·1 -
1.6 ~
~.f ~~
1. 5
,r
i 1.4
•
• ·~
,_
- ~ UNEADEa
DE DESVIACION
ll!_ 1.3
2 ¿. 1'
1,' e e ~'
..•
,
i
...1
t.2
,/
- •
tia. - " • '
....w / 11
o 4 I
O 1.10 r,/_
~ f
o 1•
...
ii: 1.08
w ;l. o ,,•
11 ;
~ 1.06 .,¡ 1
5 1.05
,
•
j '
"
•
~ 1.04
e~ 1.03 I
J
,• '
u.. I
"'
V
1.02 j
' 4
I
(
I
1.01
100 200 300 400 600 1000 2000 4000 6000 10, 000
FACTOR DE CORRELACION • F
Fig. 2.8.- Correlación de Standing para el factor volumétrico del petróleo (Ref. 2).
Un ajuste de la correlación gráfica de la Fig. 2.8 está dado por la.siguiente ecuación':
B,,h = 0.97 59 + 12 X 1 o-s F 1·2
•••••••••••••••••••••••.••••••••••••••••• (2 .3 7)
donde;
r,
- + 125 T (2.3·8)
ro
B0h Factor volumétrico del petróleo a p,,, BY/BN
:
R.,h : Razón gas disuelto-petróleo a p 2! p,,, PCN/BN
T : Temperatura del yacimiento, ºF
r 8 : Gravedad especifica del gas, ( aire = 1)
r Gravedad específica del petróleo, (agua= 1)
O
:
Standing reporta un error promedio de O.O 1 % con una desviación estándar de diferencias
de 0.34% para la correlación (ocho puntos en el rango de l 50 < F < 188f)). La Tabla 2.1
presenta el rango de los datos utilizados en su desarrollo. Además, para facilitar los cálculos la
correlación fue llevada a un nomograma el cual se presenta en la Fig. 2.9.
Carlos Bánzer S.
73
PROPIEDADES FISICAS DEL PETRÓLEO
EJEMPlO
DETERMINAR:
R factor volum4trico • 200 •f. U1t pellMeo • 30 •API
coa una ruclll IU dltu1ho patr41eo • 350 PCNIIN, y w
pavedad ttpecfflc• clel III de 0.75.
PROCEDIMIENTO:
Tal como lo Indica la flecha. • ,a,tir del walo, de JSO •
Fig. 2.9.- Determinación del factor volumétrico del petróleo, 80, por la correlación
de Standing (Ref. 2).
Correlación de Vázquez, M.E. y Beggs, H.D.ª:
La correlación es la siguiente:
r ge:
donde;
B0h : Factor volumétrico del petróleo a ph, BY/BN
Rs1, : Razón gas disuelto-petróleo a p 2= pb, PCN/BN
T : Temperatura del yacimiento, ºF
r g,· : Gravedad específica del gas corregida a 100 lpcm., (aire= l)
r APJ : Gravedad del petróleo, 0 API
Vázquez y Beggs reportan un error promedio de 4.7% para la correlación. La Tabla 2.4
presenta el rango de los datos utilizados en su desarrollo. La Tabla 2.15 presenta los valores de
las constantes C1, C2 y C3 de acuerdo a la gravedad API del crudo.
e,
Y API S 30 4.677 X 10-4 1.751 X 10-5 -1.8106 X 1018
Y API > 30 4.670 X 10-4 1.100 X 10-S 1.3370 X 10~9
Carlos s•nzer S.
74
CORRELACIONES NUMÉRICAS P.V.T.
B0b = 1.022+4.857 10-4 R,.,. -2.009 I0-6(T-60) Y API + 17.569 X 10-9 R (T-60) Y API
X
Jl1
X
r g sh r g
.......... (2.42)
donde;
B0h : Factor volumétrico del petróleo a Pb, BY/BN
Rsh : Razón gas disuelto-petróleo a p ~ Pb, PCN/BN
T : Temperatura del yacimiento, ºF
rg : Gravedad específica del gas, (aire= 1)
r APJ : Gravedad del petróleo, ºAPI
El estudio estadístico realizado por la TOTAL determinó que el 97% de 289 valores de
B0h calculados con la correlación, presentó un error del 10% en comparación con los valores
experimentales.
Correlación de AI-Marhoun, M.A.11:
La correlación es la siguiente:
Bob = 0.497069 + 0.862963 x 10-3 T + 0.182594 x 10-2 F + 0.318099 x J 0-5 F2 ••• (2.43)
donde;
F = Rsh 0.14239 r 8 o.323294 y" -1.20204 ••••••••••••••••••••••.•••••••••••••••• (2.44)
Al-Marhoun reporta un error promedio relativo de -0.01 % con una desviación estandar
de 1.18% respecto a los valores esperimentales determinados con la correlación. La Tabla 2.8
presenta el rango de los datos utilizados en el desarrollo de la correlación,
carios Bánzer S.
75
PROPIEDADES F(SICAS DEL PETRÓLEO
Dokla y Osman reportan un error promedio relativo de 0.023% ccn una desviación
estandar de 1.681 % respecto a los valores experimentales determinados con la correlación. La
Tabla 2.9 presenta el rango de los datos utilizados en el desarrollo de la correlación,
Correlación de Petrosky, G.E., Jr. y Farshad, F.f .13:
La correlación es la siguiente:
B0b = 1.0113 + 7.2046 X 10-s F3·0936 •••••••••••••••••••••••••••••••••••• (2.47)
donde;
Las correlaciones anteriores pueden ser utilizadas para estimar valores de BO a presiones
por debajo de la presión de burbujeo. En este caso, se debe utilizar R, a la presión deseada en
vez de Rsb • Además, se debe considerar la posible variación de las gravedades del petróleo de
tanque, y API , y del gas en solución, y g, con presión.
Los valores B0 por debajo del punto de burbujeo están afectados tanto por la solubilidad
como por la compresibilidad, mientras que por arriba del punto de burbujeo la solubilidad es
constante y por tanto solo influye la compresibilidad. Por lo tanto, si se conoce la compresi-
bilidad del petróleo puede determinarse B0 a presiones mayores que la presión del punto de
burbujeo (p > ph) mediante la siguiente ecuación':
B" =Bh 0 exp]c, .(ph - p)] (2.51)
donde;
B0 : Factor volumétrico del petróleo a p > ph, BY/BN
B0h : Factor volumétrico del petróleo a ph, BY/BN
Pb : Presión de burbujeo, lpca.
p : Presión de interés, lpca.
1
e" : Compresibilidad isotérmica del petróleo, lpc -
Ejemplo 2.3:
Utilizando las correlaciones de Standing, M.B., Vázquez, M.E. y Beggs, H.D., Glaso, O;
TOTAL, C.F.P., Al-Marhoun, M.A., Dokla, M.E. y Osman, M.E., Petrosky, G.J~., Jr. y Farshad,
F.F., y de Kartoatmodjo, T. y Schmidt, Z., determinar el factor volumétrico del petróleo a
presiones de 2000, 2500 y 3000 lpca. dado los siguientes datos:
Pb = 2500 lpca.
T = 180 ºF
y API = 31 ºAPI
r g = 0.95 (aire= 1)
c0 = 9.61 X 10-6 lpc"'
Solución:
Correlación de Standing, M.B.
De las Ecs. (2.37) y (2.38):
Para p = 2500 lpca. (p = Pb ), entonces B,, = B,>h
0.95
F = (673) O + (1.25)(180) = 927.85832
.8 71
B0 = B0h = 0.9759+(12 X 10-5)(927.85832)1'2
B, = Bob = 1.41258 BY/BN
Para p = 2000 lpca. (p < Pb ).
0.95
F = (516) + (1.25)(180) = 763.89286
o.8 71
B0 = 0.9759 + (12 X 10-5)(763.89286)1·2
B0 = 1.32170 BY/BN
De la Ec. (2.51 ):
Parap = 3000 lpca. (p > P» ).
Carlos Bánzer S.
77
PROPIEDADES FISICAS DEL PETRÓLEO
B0 = 1.25031 BY /BN
Parap = 3000 lpca. (p > Ph ).
B0 = 1.31755 BY /BN
Correlación de la TOTAL, C.F.P.
De la Ec. (2.42):
Para p = 2500 lpca. (p = Ph ), entonces B" = B oh
31
B0 = B h = 1.022 + (4.857 X 10-4)(689)-(2.009
0 X 10-6)(180-60) + (17.569 X 10-~)(689)(180-60) 31
095 0.95
B, = Bob = 1.39618 BY/BN
Para p = 2000 lpca. (p < Pb ).
~ariosB~nzerS.
78
CORRELACIONES NUMÉRICAS P.V.T.
31 31
B0 = 1.022 + (4.857 X 10-4)(539)-(2.009 X 10-6)(180-60)
0.95
+ (17.569X 10-6)(539)(180-60)
0.95
B0 = 1.31301 BY/BN
Para p = 3000 lpca. ( p > Pb ).
B0 = 1.38949 BY/BN
Correlación de AI-Marhoun, M.A.
De las Ecs. (2.43) y (2.44):
Parap = 2500 lpca. (p= Pb), entonces B0 = B0b
F = (840)º·'4239 (0.95)0.J23294 (0.871)-1.20204 = 172.12916
B0 = Bob = 0.497069 + (0.862963 x 10-3)(640) + (0.182594 x 10-2)(172.12916) + (0.318099 x l10-s)(172.12916)2
B0 = 1.35804 BY/BN
Para p = 3000 lpca. ( p > Pb ).
B0 = 1.4.5092 BY/BN
Correlación de Dokla, M.E. y Osman, M.E.
De las Ecs. (2.45) y (2.46):
Para p = 2500 lpca. (p = Pb ), entonces B0 = B0b
F = (882)º·773572 (0.95)º·40402 (0.871)-o.ss2605 = 210.12232
B0 = Bob = 0.431935 X 10-I + (0.156667 X 10-2)(640) + (0.139775 x 10-2 )(210.12232) + (0.380525 X 10-S)(210.12232)2
B" = B b = 150757
0 BY/BN
Para p = 2000 lpca. (p < Pb ).
F = (648) 0.773S72 (0.95) 0.40402 (0.871)-0.882605 = 165.53 730
B0 =0.431935x 10-1 +(0.156667 x I0-2)(640)+(0.139775x 10-2)(16S.53730)+(0.380S25x I0-5)(16553730)2
B0 = 1.38152 BY/BN
Parap = 3000 lpca. (p > Ph ).
B0 = 1.50034 BY /BN
Correlación de Petrosky, G.E., Jr. y Farshad, F.F.
De las Ecs. (2.4 7) y (2.48):
Para p = 2500 lpca. (p = P» ), entonces B" = B,>h
o 950.29)4
F = (624)º·3738 •
0.8710.6265
+ (0.24626)(180)º·53'1 = 15.91633
B0 = B0b = 1.0113 + (7.2046 x 10-5)(15.91633)3·0936
Carlos Bánzer S.
79
PROPIEDADES FISICAS DEL PETRÓLEO
B0 = 1.26787 BY/BN
Parap = 3000 lpca. (p > Ps ).
B0 = l.31463BY/BN
Un resumen de los resultados de B0 es como sigue:
80, BY/BN
En otras palabras, este factor incluye el volumen líquido, B0, más el volumen de la
diferencia entre la razón inicial gas disuelto-petróleo, Rs;, y la razón gas disuelto-petróleo a las
condiciones actuales del yacimiento, R; Si el factor volumétrico del gas es Bg en BY/PCN,
para ~l. 8ias en solución, el factor volumétrico total puede determinarse mediante la siguiente
ecuacion :
B, = B,, + ( R.,i - R, ). B R ••••••••••••••••••••••••••••••••••••••••••••• (2.52)
donde;
B, : Factor volumétrico total, BY/BN
B0 : Factor volumétrico del petróleo, BY/BN
Bg : Factor volumétrico del gas a p y T, BY/PCN
Rs; : Razón gas disuelto-petroleo a p;, PCN/BN
R4 : Razón gas disuelto-petróleo a p, PCN/BN
La Fig. 2.1 O ilustra el comportamiento típico de B0 y B, vs. presión a temperatura
constante.
i 2.40
~
CD
•
.J 220 -----t------+-------+------t------__.
~
~ 2.00
IIJ
.-
L&J
1.ao---- ................. -
a.
~ 1.60----- (Rsi-Rs).BcJ_
e,~'º
8- a ~
o::
1-
L&J 1.40-----..--
i
.......--~--- ........-----
2
::, •!
....._ ....,.___ ..,....__ z __ __.., ~
120
. o
ü5
w
~
a:
Q.
f
~
1.00
o 2000
1000 3000 4000 · 5000
PRESION, lpca
Fig. 2.10.- Comportamiento tfpico dé B0y 81 vs. presión
a temperatura constante.
Por arriba de la presión de burbujeo, Rs; = R,, y el .factor volumétrico total es igual al
factor volumétrico del petróleo, B, = B0• Sin embargo; por debajo de la presión de burbujeo a
medida que la presión disminuye el B0 disminuye pero el B, aumenta debido a que el término
(Rs1 - Rs).B, aumenta puesto que B1aumenta y R., disminuye mientras que R.,1 es un valor
constante.
·carios Sánzer S.
<
81
PROPIEDADES FISICAS DEL PETRÓLEO
82
CORRELACIONES NUMÉRICAS P.V.T.
= 2000 =3.04
e; 657.41 De la Fig. 1.4: z = 0.68
T.. = 180 + 460 = 1.37
sr 465.47
De la Ec. ( 1.69):
B, = 1.17366 BY/BN
Correlación de AI-Marhoun, M.A.
De las Ecs. (2.55) y (2.56), para s, = 61 S PCN/BN (Ejemplo 2.2):
F = (61 S) o.644SJ6 (0.9S)-1.01934 (O.S7 l) o.124s14 (2000)-o.16191 ( 640)2.00621 = 7,8115.59
B, = 125992 BY /BN
Un resumen de los resultados de B, es como sigue:
Carlos Sánzer S.
83
PROPIEDADES FfSICAS DEL PETRÓLEO
o
..,"'::,
ID
60 a:
::,
CD
!
z
-
o
en
l&I
o:
A. :
1
!
Carlos Bánzer S.
84
CORRELACIONES NUMÉRICAS P.V.T.
donde;
c0 Compresibilidad isotérmica del petróleo, lpc'"
:
· Carlos Bánzer S.
85
PROPIEDADES FISICAS DEL PETRÓLEO
Correlación
Notesé que la Ec. (2.62) se transforma en la Ec. (2.58) para presiones mayores que la
presión de burbujeo debido a que R, es constante con respecto a la presión. Especial cuidado
debe tenerse en la evaluación de la Ec. (2.62) debido a que Bg se calcula en una variedad de
unidades.
La Fig .2.12 presenta el gráfico completo de c0 vs. presión a temperatura constante.
Como puede observarse, existe una discontinuidad a la presión del punto de burbujeo debido a
que la liberación de la primera burbuja de gas causa un cambio drástico en el valor de la
compresibilidad.
Cario• Bánzer S.
86
CORRELACIONES NUMÉRICAS P.V.T.
----+---~ti----......_ PRESION
..,, DE BURBUJEO--.....,_-----t
1 s
1
Carlos Bánzer S.
87
PROPIEDADES FISICAS DEL PETRÓLEO
iJR.s
ép T
= C4 Rs
p
(2.74)
AI-Marhoun, M.A.
éJB" 1.35556 X 10-3 F + 4.723 X l 0-6 F2
= (2.75)
iJR..,. T s;
iJR.s = 1.3984 R.. . ,. (2. 76)
ép 1' p
Carlos ,Bánzer S.
88
CORRELACIONES NUMÉRICAS P.V.T.
Carlos Bánzer S.
89
PROPIEDADES FISICAS DEL PETRÓLEO
= 2000 = 3.18
r; 629.1 O
De la Fig, 1.4: z = O. 762
T = 180+460 = 1.47
sr 434.44
De la Ec. ( 1.69):
B8 = 0.00503 ( O.
762 )( 640) = 0.00123 BY /PCN
2000
Standing, M.B.
De las Ecs. (2.67) y (2.68):
C0 =- l [(567.3 X 10-6)(0.30693)-(0.00123)(0.30693)]
1.32170
Finalmente;
C0 = 153.9 X ] 0-6 lpc "
Vézquez, M.E. y Beggs, H.D.
De las Ecs. (2.69) y (2.70):
31
iJBO =4.67 x 10-4 +(1.337 x 10-9)(180-60) =4722 x 10-6 BY,f>CN
iJR..,· T 0.95
447
iJRs = 1.187 = 0.26529 PCN/BN-lpc
ép 1' 2000
1Carlos Bánzer S.
90
CORRELACIONES NUMÉRICAS P.V.T.
De la Ec. (2.66):
1
Co = - 1.25416 [(4722 X 10-6)(026529)-(0.00123)(0.26529)]
Finalmente;
c0 = 160.3 x 10-6 lpc-1
Glaso, O.
De las Ecs. (2.41 ), (2. 71) y (2. 72):
0.526
F = (433) 0.95 + (0.968)(180) = 627.47242
0.871
/JBO = [2.913 29 - 05 5 3 66 )og( 62 7, 4 724 2)] 1 0(-6.SSS 11 + 2.91329 log(627.47242)-0.27683 (log( 627 .4 7242))2]
iJR.s T
O.S26
l
627.47242
º·0.871
95 =569.7 x 10-6 BY/PCN
433
8R.s = 2.02777 = 0.24326 PCN/BN-lpc
qJ 1. (2000)(14.1811 - 3.3093 log(2000)]º·5
De la Ec. (2.66):
C0 =- l [(569.7 X 10-6)(0.24326)-(0.00123)(0.24326)]
1.25031
Finalmente;
e; = 128.5 x 10-6 lpc -t
TOTAL, C.F .P.
De las Ecs. (2.73), (2.74):
31
IJBO =4.851 x 10-4 +(17.569 x 10-9)(180-60) =554.5 x 10-6 BY/PCN
iJR.s T 0.91
8R.s
qJ T
= 1.095 ;::o = 0.2951 O PCN/BN-lpc
De la Ec. (2.66):
C0 =- l [(5545 X 10-6)(0.29510)-(0.00123)(0.29510)]
1.31301
Finalmente;
e; = 151.8 x 10-6 lpc "
AI-Marhoun, M.A.
De las Ecs. (2.43), (2.75) y (2.76):
Carlos Bánzer S.
91
PROPIEDADES FISICAS DEL PETRÓLEO
= 1.3984 = 0.43001
IJR.s
qJ T 2 ~:o PCN/BN-lpc
De la Ec. (2.66):
C0 =- l [(444.2 X 10-6)(0.43001)-(0.00123)(0.43001)]
1.35804
Finalmente;
c =248.8 x 10-6 lpc-1
0
C0 =- l [(5252 X 10-6)(0.44748)-(0.00123)(0.44748)]
1.38152
Finalmente;
c0 =2283 x 10-6 lpc"
Petrosky, G.E., Jr. y Farshad, F.F.
De las Ecs. (2.48), (2. 79) y (2.80):
o 950.2914
F = (492)º·3738 • + (0.24626)(180)º·5371 = 14.90308
0.8710.6265
éJB o 950.2914
-º = (83.313 X 10-6)(14.90308)2·0936(492)-0·6262 • = 527.8 X 10-6 BY/PCN
éJR.s T 0.87 J 0.626S
492
IJR.s = 1.73184 = 025127 PCN/BN-lpc
qJ 1· 2000 + 139.05
De la Ec. (2.66):
e =- l [(527.8 X 10-6)(025127)-(0.00123)(025127)]
º 1.31838
Carlos Bánzer S.
92
CORRELACIONES NUMÉRICAS P.V.T.
Finalmente;
e" = 13.3.8 x 10-6 lpc "
Kartoatmodjo, T. y Schmldt, Z.
De las Ecs. (2.50), (2.81) y (2.82):
F = ( 435)º·755 (0.95)º·25 (0.871)-l.S + (0.45)(180) = 200.25128
IJBO = 1.1325 x 10-4 (200.25128)1.s -(0.45)(180X200.25128)0.S =4393 x 10-6 BY/PCN
8R.s T 435
435
8R.s = l = 0.23789 PCN/BN-lpc
ép T 0.9143 2000
De la Ec. (2.66):
C0 =- l [(4393 X 10-6X0.23789)-(0.00123X0.23789)]
1.26787
Finalmente;
C0 = 148.4 X 10-6 lpc-l
Un resumen de los resultados de c0 es como sigue:
Correlación C0, 1 pe -1
McCain, W.D., Jr. y Cols.
Ec. (2.63) 143.6 X 10-6
Ec. (2.64) 146.7 X ]0-6
Ec. (2.65) 162.4 X 10-6
Standing, M.B. 153.9 X 10-6
Vázquez, M.E. y Beggs, H.D. 160.3 X 10-6
Glaso, O. 1285 X 10-6
TOTAL, C.F.P. 151.8 X )0-6
AI-Marhoun, M.A. 248.8 X 10-6
Dokla, M.E. y Osman, M.E. 228.J X 10-6
Petrosky, O.E., Jr. y Farshad, F.F. 133.8 X 10-6
Kartoatmodjo, T. y Schmidt, Z. )48.4 X 10-6
Carlosl Bánzer S.
93
PROPIEDADES FISICAS DEL PETRÓLEO
2.80
Q.
u
...
o
_.
l&.I 2.40
o
et::
t-
l&J 2.00
a.
...J
LIJ
o
o 1.60
s-
U)
o 1.20
u
->
CI) PRESION DE SATURACION
--
l-
oao
Carlos .Bánzer S.
94
CORRELACIONES NUMÉRICAS P.V.T.
ao.ooo,-
D.OOOI-.....P....~~+---+--t---+--41---1---+-+-_,___,__...,._...._,1--......._,-""-I--'-....._.
,.coa,--.~.-....,;: ...--.- -1--t--t--+--f--1-......_4-,o.,......_.._,J---1---1--+--+-...._..
1.oac,1....,___~
ICOOt---+--
ao,....,_~...,..-~....__
U
4
. IOt--,,.,i
••t---f--lloi
~~~_,p ........
.....;
-
u
t¡ 101- -~ ~~-IJ.ó111,p.t,,.-.;
2
w
z
o
o
.,,-
~
8!! a,q--+o--+--+--+---t~+--~w
>
9CJ too IIO 140 180 ,., IDO 240 NO JIO NO «,o 440 410 B10
TEMPERATURA, GRADOS FAHRENHEIT
LEYENDA
I e • C,U,O O[ CAUl'ONN1A (9.,VI) I O .. CRUOO ,r400 DO. NOHII Dt ~-
2 6 • CRUDO COLOM•ANO (11.~API)
flO.ftPI) 9 • • CRUDO ROJO Ott. COHTffCNTE MEOtO
l D • ll!SttlJO OEL CCJNTINENTE MIOfO
4 V • CRUDO DE CALIFORNIA
114.nA)
(114.#tll
10 * · CftUOO DEL SUR DE TEXAS
11 • • CRUDO DE LA COSTA DEL GOl.10
C24.0'APtJ
(l?.8W'f)
tzo.!1-..)
5 O• CRUDO OIL tONTINfNT! MEDIO flMA 12 A • FETflOI.IO LIVfANO flAIWPINtCO
ACEITI PARA OLINOROS (22... Afll) DEL CONTINENTE MtOIO
t e • AC!fTE PESADO DI\. CONTINENTE (2t.71AP4)
MEDfO PMA MOIONIS (llraPIJ
13 • - DESTILADO AL YAaO DE fllNSILVAMA
14 e • CRUDO DI WYOMING
csz.1•••>
(2t.4•Afll)
1 • • OUX>De LA ms,a IM!l GOUO m IIECICO Cts.1•APIJ IS Y• OISTft.ADO D! ,1LT11-o, ..,ICENSA
DEL CXNTININTE MEDIO CZO! A,t)
Fig. 2.14.- Viscosidad de petróleos libres de gas vs. temperatura (Ref. 19).
Como la densidad del petróleo no depende fuertemente de la temperatura, frecuente-
mente se hacen gráficos como el de la Fig. 2.14 utilizando la viscosidad dinámicaµ vs.
temperatura. Tales gráficos son útiles (los errores generalmente son menores del 25% sobre
los rangos prácticos de temperatura) en proporcionar buenos estimados de viscosidades de
crudos a otras temperaturas diferentes de las medidas.
Las cartas ASTM D-341 se pueden obtener en varios tamaños, unidades e intervalos
de viscosidad y temperatura. La Fig. 2.1 S presenta un ejemplo de esta carta.
Ejemplo 2.7:
La viscosidad de un crudo de 26.7 ºAPI es de 5.2 cts. a 140 ºF y de 1.5 cts, a 280 ºF.
Determinar la viscosidad a 160 ºF.
Solución:
Determinación de la constante B (Ec. 2.86):
-B = log[log( 1.5 + O. 7)] - log[log(S.2 + O. 7)]
log( 140 + 460) - log(280 + 460)
-B = 3.869
Carlos Bánzer S.
96
CORRELACIONES NUMÉRICAS P.V.T.
--
• -
1-
....
es ... ·: 1
i: :1
t:
1 -1
..
- -
- ~
•
TIMPIAATIHII , DKIU!ES FAKAINHl!n
±
;~}
log[log( v.
I
+ O. 7)] = log[log(51 + 0.7)] + log[log(l.5 + 0.7)] = -0.578
±(log
i:)
T)
I
2 = (log(140 + 460))2 + (log(280 + 460))2 = 15.948
2
±
'=]
log TI
= (log(l40 + 460) + log(280 + 460))2 = 31.889
De la Ec. (2.89):
1
-1.648 - i (5.647)(-0.578)
B= l =-4.0
15.948- 2 (31.889)
Carlos Bánzer S.
97
PROPIEDADES FfSICAS DEL PETRÓLEO
De la Ec. (2.88):
-0.578-(-4.0)(5.647) .
A=
2
= 11.005
Sustituyendo A y Ben la Ec. (2.85), se tiene:
log[log( v + 0.7)] = 11.005- 4.0 log T
La viscosidad a 160 ºFes:
log[log( v + 0.7)] = 11.005 - 4.0 log(l 60 + 460) = -0.165
luego;
V+ 0.7 =4.84
Finalmente;
v = 4.14 cts.
Las siguientes correlaciones permiten determinar la viscosidad de crudos sin gas en
solución (muertos), µ0d en cp., a presión atmosférica y temperatura del yacimiento,
Correlación de Beal, C.20:
Esta correlación se presenta en la Fig, 2.16 donde la viscosidad del petróleo libre de
gas (muerto) está en función de la gravedad API del petróleo y a varias temperaturas.
En el desarrollo de esta correlación se utilizó un total de 655 datos de viscosidades de
crudos sin gas en solución a 100 ºF y 98 a temperaturas mayores de 100 ºF obtenidos de
492 campos diferentes, 358 de los cuales correspondían a los Estados Unidos. La Tabla 2.18
presenta el rango de estos datos. Beal reportó un error promedio de 24.2o/o entre los valores
leídos del gráfico y los valores experimentales.
Beal, C. Glaso, O.
Temperatura, ºF 98 - 250 50 - 300
Gravedad del petróleo, ºAPI 10 - 52.5 20.1 - 48.1
Viscosidad, cp. 0.865 - 1.55 0.616 - 39.1
Un ajuste de la correlación gráfica de la Fig. 2.16 está dado por la siguiente ecuación":
7 a
1.8 X 10 360
µod = 0.32 + 4.53 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ( 2. 90)
Y Af'I T+200
donde;
8.33
a=antilo 0.43+ (2.91)
Y API
µ0d : Viscosidad del petróleo libre de gas a 1 atm. y T, cp.
y AJlJ : Gravedad del petróleo, º API
T : Temperatura del yacimiento, ºF
Para 25 puntos leídos de la Fig, 2.16, la ecuación anterior presentó un error promedio
de -l.58o/o con una desviación estandar de diferencias de 9.4%3•
Carlos Bánzer S.
98
CORRELACIONES NUMÉRICAS P.V.T •
•. • •. •
•
. •
.
.......
' .
¡ • •
e
iI!
¡ ,o•: . . .
. .
o '
• .. .' '.
.. .
' ~~--.•. - ~. -.- ----.._
•
.
... - ..._.. ---"·
'
__• .
• • • < • '
......_.. : • 1 • • • ~· •
.•
a•~ •;
.. .
..
• • • . .. . ... .
' > •• -
-
~ • • $ ~ •
• ....•• • •
-·-
...
•• • • - .. • 1 • •< • i
.
.
• .
·-
•
....
•
.
•
. -~·"·· ~-·· ..
.......
4._ ••••.• ~:
•••• .. • , 1
• • •
. ..... ·~··-
• ••• 1 •.• __.. •
......... t " • • •
... • • • •• • :i,
'º ao •o •o • • 70
GRAVEDAD DEL PETROLEO. •API A 80 •f Y PRl!SION ATMOSFERICA
Carlos Binzer S.
99
PROPIEDADES FISICAS DEL PETRÓLEO
Un total de 460 datos de viscosidades de crudos sin gas en solución obtenidos de 600
muestras diferentes de crudos fueron utilizados en el desarrollo de esta correlación. La Tabla
2.19 presenta el rango de estos datos. Beggs y Robinson reportan un error promedio de
-0.64%con una desviación estandar de 13.53% entre las viscosidades medidas y calculadas
con la correlación.
µod
3
= · 141 X
10 10 T-3.444 (1og Y API )(10.313 log T-36.447) ••• ·• • •• • •• • •• • • • • .... ••
(2 93)
•
donde;
µod : Viscosidad del petróleo libre de gas a l atm. y T, cp.
r APJ : Gravedad del petróleo,
0
API
T : Temperatura del yacimiento, ºF
Esta correlación fue desarrollada a partir de 26 muestras diferentes de crudos
cubriendo el rango de datos mostrados en la Tabla 2.18.
Correlación de Egbogad, E.0.22:
La correlación es la siguiente:
log[log(µ,"' + 1)]= l.8653-0.025086y APJ -0.5644 logT , (2.94)
donde;
µ0J : Viscosidad del petróleo libre de gas a I atm. y T, cp.
r APJ : Gravedad del petróleo,
0
API
T : Temperatura del yacimiento, ºF
Egbogad reporta un error promedio de 6.6% para la correlación. El rango de los datos
utilizados en el desarrollo de la correlación presentaban las siguientes características:
59 < T(º F) < t 76 y s < r APJ <º API) < ss.
Correlación de Kartoatmodjo, T. y Schmldt, z.14:
La correlación es la siguiente:
µod = 16 •O X 10 8 T-2.8177 () ogy API )(S.7S26 logl'-26.9718)
••••••••••••• • • •• •.....
(2 • 95)
donde;
µoJ : Viscosidad del petróleo libre de gas a 1 atm. y T, cp.
r AJ>J : Gravedad del petróleo, 0 API
T : Temperatura del yacimiento, ºF
En el desarrollo de esta correlación se utilizó un total de 661 muestras diferentes de
crudos, y reporta un error promedio de -13.158%. La Tabla 2.20 presenta el rango de los
datos utilizados en su desarrollo.
Carlos Bánzer S.
100
CORRELACIONES NUMÉRICAS P.V.T.
Temperatura, ºF 75-320
Gravedad del petróleo. eAPI 14.4 - 68.9
Viscosidad, cp. 0.5-682
cñ~
3 ~ 10
~i :~=+=~~=~~~~==~=~=t4==~~~±j=s.'-'
2 ?t---t--t-+-+-+-l-t----f--l--+-+--1--1--1-
c ~ •r-.......-t-1-1-t-rt----+--+--l-+--1--'-
i~
g¡~ '
4t--t--t-~i-+-t-+----+--l-~l--
o f l,-.-t--t--t--1-+-t-+----+--1--"-
> ltAION GAi i*SUl'LfO•l'lnlOUO.
~ ~ lr--t--t-+-+-t-t-+---+-- lfCN/8N ~·,.. ~........iM
~s
!o ~!u 1.0.-- .......................
•t--+-+-t-+-+
e ..J •
~ ~ .?t---t--t-+-
0
u
en A:
~ ·---~
> ¡:
- .... 4
---
e
EJEMPL
1.0 IO IOO
VISCOSIDAD DEL PETRÓLEO MUERTO, cp
( A TEMP. DEL YACIMIENTO Y PRESION ATMOSFERICA)
Fig. 2.17.-Viscosidad del petróleo saturado de gas a presión atmosférica
y temperatura del yacimiento (Ref. 23).
Carlos Bánzer S.
101
PROPIEDADES FISICAS DEL PETRÓLEO
11.1 ~
o. PENbltNTt t. ,.,.
-
Q.
o ~
1) -
~ as
Q. ~
-·-
~,~Jj
~ 0.4
u .
«
~ o.s
~
--......!~
.....
A
-
-
z .....
-
0.1
-o 1 1 t 1 1 1 1 1 1 • 1 1 • t100 1 1 1
-
IOO IOO ,00 400 IOO IOO 100 IOO 'IOO IODO IIOO ll(líiIU 1400 l900 •oo
RAZON GAS DISUELTO- PETROLEO, PCN/BN
Fig. 2.18.- Valores del intercepto A y pendiente b utilizados en la correlación
de Chew y Connally (Ref. 23).
Carlos Bánzer S.
102
CORRELACIONES NUMÉRICAS P.V.T.
correlaciones.
Correlación de Beal, c.20:
'
Esta correlación se presenta en la Fig. 2.19 y permite determinar el incremento en
viscosidad, ~µ0, por cada 1.000 lpc de aumento de presión por arriba de p11•
En el desarrollo de esta correlación se utilizaron 52 datos de viscosidades tomados de
26 muestras de crudos de 20 campos diferentes, 11 de los cuales correspondían a California.
La mitad de los datos fueron tomados de crudos a la presión de burbujeo, mientras los datos
restantes fueron obtenidos de crudos a presiones mayores que la presión de burbujeo. La
Tabla 2.23 presenta el rango de los datos utilizados en el desarrollo: de la correlación, Beal
reporta un error promedio de 2.7o/o para la correlación.
Carlos Bánzer S.
103
PROPIEDADES F(SICAS DEL PETRÓLEO
6
4
.
2 ~
-en
l&J 10
a: .
CL .r
<( 6 lt.
.11
y
_J
l&Jo 4
o g_
en -
~s
º-
2 ,
<t, . 1
e o.
CI)u
i
=
8
->
(1) 6
4 .
:5
~ 2
z
~
L&J
2 .
LAJ
a: 6
u
-
z 4
-· ·-
• 4 6 1 2 4 6
' 10 2 4 101 •L,.• 3
Carlos Bánzer S.
104
CORRELACIONES NUMÉRICAS P.V.T.
Carlos Bánzer S.
105
PROPIEDADES FISICAS DEL PETRÓLEO
8.33
a= antilo 0.43 +
31
= 4.997
4.997
1.8 X 107 •
360
µod = 0.32 + 4 53
31 · 180+ 200
luego;
µ"J = 2.65 cp.
Correlación de Beggs, H.D. y Robinson, J.R.
De la Ec. (2.92):
z = 3.0324 - (0.02023)(31) = 2.405
y= 102·4º5 = 254.097
X= (254.097)( 180)-l.163 = 0.606
luego;
- 100.606
µod - -
1
Finalmente;
µ"d = 3.04 cp.
Correlación de Glaso, O.
De la Ec. (2.93):
µ od = (3 .141 X 1 OJO )(180)-3.444 (log(3 l ))
(10.313 log(JS0)-36.447)
luego;
µod = 2.76 cp.
Correlación de Egbogad, E.O.
De la Ec. (2.94 ):
log[log( µ0d + 1)] = 1.8653 - (0.025086)(31 )- (0.5644) log(l 80)
log[log( µ"d + 1)] = -0.18524
luego;
µ0d + 1 = 4.49539
Finalmente;
µod = 3.50 cp.
Correlación de Kartoatmodjo, T. y Schmidt, Z.
De la Ec. (2.95):
«: = ( 16.0 X l 08 )(l S0)-2.8177 (Jog(J) )) (5.7526 log(180)-26.9718)
luego;
µod = 2.63 cp.
Carlos Bánzer S.
106
CORRELACIONES NUMÉRICAS P.V.T.
Crudo Saturado:
Correlación de Chew, J.N. y Connally, C.A., Jr.
De las Ecs. (2.97) y (2.98):
A= antilog[(675)((2.2 x io:' )(675)- 7.4 x 1 o-4 )] = 0.399
0.68 0.25 0.062
b = 1 o(S.62 X
s
10- )(67S)
+ 1 o(l.l -' . \
X 10. )(675)
+ 1 o<J.74 X
-3
10 )(67S)
= 0.640
De la Fig. 2.17 ó Ec. (2.96) para µ,,J = 2.65 cp. (Beal, C.):
µ0h = (0.399)(2.65)0·640
luego;
µoh = 0.74 cp.
Correlación de Beggs, H.D. y Robinson, J.R.
De las Ecs. (2.100) y (2.1 O 1 ):
a= 10.715 (675 +) 00)-0.SIS = 0.348
b = 5.44 (675 + ISO)-O.J3S = 0.562
De la Ec. (2.99) para µ d = 2.65 cp, (Beal, C.):
0
µ0¡, = (0.348)(2.65)0'562
luego;
µoh = 0.60 cp.
Correlación de Kartoatmodjo, T. y Schmldt, Z.
De las Ecs. (2.103) y (2.104) para µ,,J = 2.65 cp. (Beal, C.):
b = I o<-o.ooos1)(6?s> = o.28396
A= co.200 t + o.8428 x t o<-0.ooos4s)c61s> )(2.6s)<º.43+(0.s16s)(o.2s396)) = o.74859
De la Ec. (2.102);
µ"b = -0.0682) + (0.9824)(0.74859) + (40.34 X ) o-s )(0.74859)2
luego;
µoh = 0.67 cp.
Crudo Subsaturado:
Correlación de Beal, C.
De la Fig. 2.19 ó Ec, (2.105) para µ"h = 0.74 cp. (Chew, J.N. y Connally, C.A . , Jr.):
µ" - 0.74 t.6 O.S6
(0.001)( 4000- 2500) = (0.024)(0.74) + (0.038)(0.74)
luego;
µo =0.8lcp.
Correlación de Vázquez, M.E. y Beggs, H.D.
Carlos Bánzer S.
107
PROPIEDADES FISICAS DEL PETRÓLEO
De las Ecs. (2.106) y (2.107) para µ0,, = 0.74 cp. (Chew, J.N. y Connally, C.A., Jr.):
m = (2.6)( 4000) i.1s, exp(-11.513 - (8.98 x l o-s )( 4000)) = 0.342
4000 0.342
µo = (0.74 2500
luego;
µo= 0.87 cp.
Correlación de Kartoatmodjo, T. y Schmldt, Z.
De la Ec. (2.108) para µ0h = 0.74 cp. (Chew, J.N. y Connally, C.A., Jr.):
µoh = (1.0008] )(0.74) + (1.127 X 10-3)(4000- 2500)((-65.17 X } 0-4 )(0.74)1·8148 +
(0.038)(0. 74)1'59]
luego;
u; = 0.77 cp,
Un resumen de los resultados de µ,, es como sigue:
•
¡
o•
...J
w
Q
e
<
e
UJ
ie
o
••
RAZON GAS DISUELTO-PETROLEO , PCN/BN;
IOO IOOO 1400
La tensión interfacial gas-petróleo, O'go en dinas/cm, se define como la fuerza por unidad
de longitud en la interfase entre dos fluidos inmiscibles. La tensión interfacial entre un gas y
un hidrocarburo líquido varía aproximadamente entre 35 dinas/cm a bajas presiones y
gravedades API a O dinas/cm a la presión crítica cuando ocurre solubilidad completa, Esta
propiedad es requerida para estimar fuerzas de presión capilar en cálculos de ingeniería de
yacimientos y es un parámetro utilizado en algunas correlaciones en el cálculo de flujo
multifásico en tuberías.
La siguiente correlación permite determinar la tensión interfacial gas-petróleo.
Correlación de Baker, O. y Swerdloff, W.28:
Esta correlación se presenta en las Figs, 2.21 y 2.22. La Fig. 2.21 permite determinar la
tensión interfacial del petróleo libre de gas (muerto) a presión atmosférica, O'ad.
42....-----...------...-----...-----.....---- ......---- .....---- .....----....,...----....
40
TENSION SUPERFICIAL DEL PETROLEO
38
A PRESION ATMOSFERICA
36
eo
'oe
tn 32
·-
o
--
J 28
26
24
22
20 ..... .._ ...._ _._ ..... ~ .......... ~
O 10 20 30 40 50 60 70 80 90
GRAVEDAD API A 60ºF
Fig. 2.21.- Tensión superficial del petróleo a presión atmosférica (Ref. 26).
Las siguientes ecuaciones son un ajuste de la Fig, 2.215:
0"68 = 39- 0.257]y API ••••••••••••••••••••••••••••••••••••••.•••••• (2.112)
0'100 = 37.5- Q.2571y AJ>/ ••••••••••••••••••••••••••••••••••••••••• (2.} }3)
donde;
o-68 : Tensión interfacial a 68 ºF, dinas/cm
o-100 : Tensión interfacial a 100 ºF, dinas/cm
r API : Gravedad del petróleo, 0 API
Debido a que el efecto de la temperatura sobre u d es desconocida, la extrapolación más
0
allá del rango de temperatura de 100 < T(º F) < 68 no es recomendado. Por lo tanto, se sugiere
que si la temperatura es mayor de 100 ºF, el valor a 100 ºF debe ser utilizado. De igual
Carlos Bánzer S.
110
CORRELACIONES NUMÉRICAS P.V.T.
manera, si T < 68 ºF, utilizar el valor calculado a T= 68 ºF. Para temperaturas intermedias,
utilizar la siguiente interpolación lineal entre los valores obtenidos a 68 y 100 °F5:
(T- 68Xa6s - ª100}
Ur = (768 - 32 ., (2.114)
IOOO 2000
Carlos Sánzer S.
111
PROPIEDADES FISICAS DEL PETRÓLEO
3. Standing, M.B.: Volumetric and Phase Behavior o/ Oíl Field Hydrocarbon Systems,
SPE, Dallas ( 1977).
4. Lasater, J.A.: ''Bubble Point Pressure Correlation," Trans., AIME ( 1958) 213, 379-
381.
S. Beggs, H.D.: "Oil System Correlations," Petroleum Engineering Handbook, H.C.
Bradley (ed.), SPE, Richardson, TX (1987) 1, Chap. 22.
6. Vázquez, M.E. and Beggs, H.D.: "Correlations for Fluid Physical Property
Prediction," J. Pet. Tech. (June 1980), 968-970.
7. Glaso, O.: "Generalized Pressure-Volume-Temperature Correlations," J. Pet. Tech.
(May 1980), 785- 795.
8. Jacobsen, H.A.: ''The Effect of Nitrogen on Reservoir Fluid Saturation Pressure, '' J.
Cdn. Pet. Tech. (July-Sept., 1967) 101-105.
9. Ostermann, R.O., Ehlig-Economides, C.A., and Owolabi, O.O.: "Correlatíons for the
Reservoir Fluid Properties of Alaskan Crudes," paper SPE 11703 presented at the
1983 California Regional Meeting held in Ventura, California, March 23-25, 1983.
1 O. TOTAL Compagnie Francaise Des Petroles: "Proyectos de Inyección de Fluidos -
Correlaciones P.V.T. para Crudos del Oriente de Venezuela," S.A. MENEVEN, Sept.
de 1983.
11. AI-Marhoun, M.A.: ''PVT Correlation for Middle East Crude Oils," J. Pet. Tech.
(May 1988), 650-666.
12. Dokla, M.E. y Osman, M.E.: "Correlation of PVT Properties for UA:E Crudes,"
Trans.; AIME (1992) 293, 41-46.
13. Petrosky, G.E., Jr. and Farshad, F.F.: "Pressure-Volume-Temperature Correlations
for Gulf of Mexico Crude Oils," paper SPE 26644 presented at the 68th Annual
Technical Conference and Exhibition, Houston, Texas, Oct. 3-6, 1993.
14. _ Kartoatmcdjo, T. and Schmidt, Z.: ''Large Data Bank lmprove Crude Physical
__ : - Property Correlations," Oil and Gas J. (July 4, 1994) 51-55.
IS. Martin, J.C.: "Simplified Equations of Flow in Gas Orive Reservoirs and the
Theoretical Foundation of Multiphase Pressure Buildup Analyses," J. Pet. Tech. (Oct.
1959), 309-311.
16. McCain, W.D., Jr., Rollins, J.B., and Villena-Lanzi, A.J.: ''The Coefficient of
Isothermal Compressibility of Black Oils at Pressures Below the Bubblepoint,''
Trans., AIME (1988) 2~5, 659-662.
17. Gottfried, B.S.: ''A Mathematical Model of Thermal Oil Recovery in Linear
Systerns," Trans., AIME (1964) 234, 196-210.
18. ASTM Standard Viscosity-Temperature Chart 0361-39 Nat. Bur. Stds. Circular
C461, 1947.
19. Braden, W.B.: ''A Viscosity-Temperature Correlation at Atmospheric Pressure for
Gas-Free Oils," Trans., AIME (1964) 237, 1487-1490.
20. Beal, C.: ''The Viscosity of Air, Water, Natural Gas, Crude Oil and Its Associated
Gases at Oil-Field Temperatures and Pressures," Trans., AIME (1946) 165, 94-115.
21. Beggs, H.D. and Robinson, J.R.: "Estimating the Viscosity of Crude Oil Systems," J.
Pet. Tech. Forum (Sept. 1975), 1140-1141.
Cario, Bánzer S.
112
CORRELACIONES NUMÉRICAS P.V.T.
22. Egbogah, E.O.: ''An Improved Temperature-Viscosity Correlation for Crude Oil
Systems,'' paper 83-34-32 presented at the 1983 Annual Technical Meeting of the
Petroleum Society of CIM, Banff, Alberta, May 10-13, 1983.
23. Chew, J.N. and Connally, C.A., Jr.: ''A Viscosity Correlation for Gas-Saturated
Crude Oils,'' Trans., AIME (1959) 216, 23-25.
24. Katz, D.L.: "Prediction of the Shrinkage of Crude Oils," Drlll. and Prod. Prac., API
(1942), 137.
25. Brill, J.P. and Beggs, H.D.: ''Two-Phase Flow in Pipes,'', University of Tulsa, Tulsa,
OK, 1988.
26. Baker, O. and Swerdloff, W.: ''Finding Surface Tension of Hydrocarbon Liquids,"
Oil and Gas J. (Jan. 2, 1956) 125.
Car/os Sánzer S.
113
PROPIEDADES FlslCAS DEL AGUA
Carlos Sánzer S.
*Referencias al final de la sección 114
CORRELACIONES NUMÉRICAS P.V.T.
-
"O
(1)
-e
(U
-o
-·-.o:::,
·-
o
-
~
(/)
... o
V,
o ..m
•
~
~M
L&J
::, "§ l,
-
(/)
Q
-o~
(O
,:,
·-e
-
V,
o ·-
o
----+-----+--------+-----------..__ __.º....
-o
.j
mi
(O
en
_J
iJ
.9 e
i!
o j
w
KJ
O>
~
1
•
N•
ti)
•
·-
O)
8- o
o, i ed ¡
o
¡ 2d 1nº
-
LL
d d ó d o
\f~nc1 'ln!>y N3
syg 1¡a ara111en1os
\f~3nW1\tS N3 S\19 130 a\10111en1os
§ 1
i
_____ ..,..___ -+----+----++-_._.-+-...,__._~ ---+--+-- ~ti--+-1
-
•
<"')
•
C')
010.1Ph{.... -- ·-
...... ...._ ..._ _._ ...._ _..~----. . . .----------~os LL
: ~ ~ N ~ m
N8/N~d ''lijíld 130 O\t01118íl10S
\tíl9\t N3 ONVJ.3W
Carlos Bánzer S.
115
PROPIEDADES FiSICAS DEL AGUA
donde;
A = 2.12 + 3.45 x 10-3 T - 3.59 x 1 o-s r2
B = o.o 101 - 5.26 x 1 o-s T + t.48 x 10-1 r2
e = -8. 7 s x 10-1 + 3.9 x 10-9 T - 1.02 x 10-11 r2
p : Presión, lpca.
T : Temperatura, ºF
y el factor de corrección por salinidad está dado por:
RSW -4 .
=1-(0.0753-l.73xl0 T).S (3.4)
R.swp
donde, Ten ºF y S (salinidad) en porcentaje por peso de sólidos disueltos (1% = 10.000 ppm).
Ejemplo 3.1 :
Utilizando las correlaciones de Culberson, 0.L. y McKetta, J.J., Jr., y de McCoy, R.L.
determinar la solubilidad de un gas natural en una salmuera conteniendo 20.000 ppm. de
sólidos disueltos a una presión de 5000 lpca. y temperatura de 200 ºF.
Sólución:
Correlación de Culberson, O.L. y McKetta, J.J., Jr.
De la Fig. 3 .1 ó Ec. (3 .1 ), la solubilidad del gas natural en agua pura a 5000 lpca. y 200 ºF es:
Rswp = 1.84729 + (5.08 X l 0-3 )(5000)- (2.] 7767 X 10-7 )(5000)2
R.-.wp = 21.8 PCN/BN
De la Fig, 3.2 ó Ec. (3.2), el factor de corrección para una salinidad de 20.000 ppm. a 200 ºFes:
Rsw = 1 o<-0.08406SS)(2)(200)-0.285854 =0.92
Rj'WP
Carlos Bánzer S.
116
CORRELACIONES NUMÉRICAS P.V.T.
s
8
-f: I.OOL----1------..---+-----+----+--
IIJ
2
3
g 0.,01-----1----...+----t---+-----+---- _,.__ .....
a:
PRES~ DE BURBUJEO
u
~~--~1....---- ......-----
~
La figura muestra que cuando la presión del yacimiento es reducida desde su presión
inicial a la presión del punto de burbujeo, el factor volumétrico del agua aumenta debido a la
expansión del agua en el yacimiento. A presiones por debajo de la presión de burbujeo, gas es
liberado, pero el factor volumétrico continúa aumentando debido a que la disminución (merina)
en el volumen de agua, resultante de la liberación de gas, es insuficiente para contrarrestar la
expansión del líquido. Este efecto se debe a la baja solubilidad del gas natural en el agua.
Las siguientes correlaciones permiten determinar el factor volumétrico del agua a
presiones menores o iguales a la presión del punto de burbujeo (p S Pb>·
Correlación de McCaln, W.D., Jr.2:
Esta correlación se presenta en las Figs. 3.4 y 3.5, donde el cambio en volumen durante
la reducción de presión es representado por ~ Vw, , y el cambio en volumen debido a la
reducción en temperatura es representado por aVwr como funciones de presión y temperatura.
El factor volumétrico del agua puede determinarse utilizando estos valores mediante la
.siguiente
. ecuacion:.,
Bw = ( 1 + ~ Vwp )( 1 + ~ VwT) ••••••••••...••••••••...•••••..•• , ...•••... (3 .5)
Un ajuste de las curvas de las Figs. 3.4 y 3.5 está dado mediante las siguientes
•
ecuaciones:
~Vw1' = -1.0001 X 10-2 + 1.33391 X 10-4 T + 5.50654 X 10-7 T2 •••••••••••••••••••••••••••· ••••••• (3.6)
dVwp = -1.95301 X 10-9 pT-1.72834 X 10-13 p2 T- 358922 X 10-7 p-2.25341 X 10-1º p2
.......... (3. 7)
donde, p en lpca. y Ten ºF.
Esta correlación es válida para aguas de yacimientos con concentraciones de sal
ampliamente variables. Así, un aumento en la concentración de sal causa un ligero aumento en
el coeficiente de expansión térmica del agua, lo cual produce un ligero aumento en &V...r- Un
aumento en la concentración de sal causa una disminución en la solubilidad del gas en el agua,
causando una ligera disminución en ó Vwp. Estos cambios en 8 V.,r y ó V.,P están dentro del 1 %
para concentraciones de sal en el rango de O a 30%.
Correlación de McCoy, R.L.4:
La correlación es la siguiente:
B wp = A + Bp + Cp2 •••••••••••••••••••••••••••••••••••••••••••••••••• (3 .8)
donde para agua pura libre de gas,
A = 0.9947 + 5.8 X 10-6 T + 1.02 X 10-6 T2
B.= -4.228xl·0-6+1.8376xI0-8T-6.77xl0-11T2
C = 1.3 X 1 o-to - 1.3855 X 1 o-tl T + 4.285 X 1 o-IS T2
y para agua pura saturada con gas,
A = 0.9911 + 6.35 x 10-s T + 85 x 10-7 T2
B =
-1.093 X 10-6 -3.497 X 10-9 T+4.51 X 10-12 T2
C = -5.0xl0-11+6.429xl0-13T-l.43xl0-1sT2
donde, p en lpca. y Ten ºF.
Cario• Bánzer S.
118
CORRELACIONES NUMÉRICAS P.V.T.
' •
•
.,
'
'
o
..1 1
... I a
o-
u
A -...
8
i•
~
8
-
o
N ¡;
l
~
<
8- 1•
lt)
•
('I')
•
·-
D>
1.1..
o N
o
o 8d §
d
-8
N
d
• • • •
--------..-------..-------....--------.--------,------.--.§
-CI)
""'!
~ E
o
• J!!
.!
C(
----+-----+------+---·---t 8 cti:: Q)
te
"' ::, -e
e .-..-:..
~N
o:: ~ .
~ al
a..... ª"
8 •
1
•
~
•
~
•
·-
u.
C)
.,
o
..q N
8
-8 o
d o
Carlos Sénzer S.
119
PROPIEDADES FISICAS DEL AGUA
......... (3.9)
donde, p en lpca., Ten ºF y S (salinidad) en porcentaje por peso de sólidos disueltos (1 % =
10.000 ppm).
El factor volumétrico del agua a presiones mayores que la presión del punto de burbujeo
(p > Ph ), se determina de la misma manera que para un sistema gas-petróleo, esto es:
B"' = Bwb expjc; .(p,, - p)] (3.1 O)
donde;
B; : Factor volumétrico del agua a p > p¡,, BY/BN
B..,b : Factor volumétrico del agua a Pb, BY/BN
Ps : Presión de burbujeo, lpca.
p : Presión de interés, lpca.
cw : Compresibilidad isotérmica del agua, lpc "
Ejemplo 3.2:
Utilizando las correlaciones de McCain, W.D., Jr. y de McCoy, R.L., determinar el
factor volumétrico de una salmuera conteniendo 20.000 ppm. de sólidos disueltos a una presión
de 5000 lpca. y temperatura de 200 ºF.
Solución:
Correlación de McCain, W.D., Jr.
De las Figs. 3.4 y 3.5 ó Ecs. (3.6) y (3.7):
dVwT = 0.03870
~vwp = -0.0102
De la Ec. (3 .5):
B; = (1 + 0.03870)(1-0.01025) = 1.02806 BY/BN
Correlación de McCoy, R.L.
De la Ec. (3.8), considerando agua pura saturada con gas:
Bwp = 1.03780- (1.61 X l 0-6 )(5000) + (2.138 X l o-t t )(5000)2
luego;
Bwp = 1.03028 BY/BN
De la Ec. (3.9), el factor de corrección por salinidad es:
Finalmente,
B..,= (l.03028)(1.00068) = 1.03098 BY/BN
;Carlos Bánzer S.
120
CORRELACIONES NUMÉRICAS P.V.T.
e=--
l oVW =-
l ºPw =--
1
(3 .11)
V q., s;
e •• f O •• •. O• 0 •• O ••• 0 O
w
w qJ T e; T
Esta propiedad es afectada por presión, temperatura y solubilidad del gas en el agua, la
cual a su vez es afectada por la salinidad. Así, un aumento en presión causa tina reducción en
la compresibilidad, mientras que un aumento en temperatura produce un aumento en esta
propiedad. Por otro lado, a una presión y temperatura dada, el efecto del gas disuelto en el
agua es aumentar la compresibilidad en relación al agua pura a las mismas condiciones de p y
T, mientras que un aumento en la salinidad del agua disminuye la compresibilidad.
Las siguientes correlaciones permiten determinar la compresibilidad del agua a presiones
mayores que la presión del punto de burbujeo (p > Ph ).
Correlación de Dodson, C.R. y Standing, M.B.6:
Esta correlación se presenta en la Fig. 3.6 (parte superior) donde la compresibilidad del
agua pura libre de gas, cwp , está en función de presión y temperatura. La parte inferior de esta
figura presenta la corrección que debe efectuarse a cwp por efectos del gas disuelto.
La siguiente ecuación es un ajuste de las curvas de la Fig. 3.6 (parte superior)':
A+BT+CT2
Cwp = 106 ,. (3.12)
donde;
A = 3.8546- t.34 x 10-4 p
B ,= -0.01052 + 4.77 x 10-7 p
C = 3.9267 X 10-S -8.8 X 10-JO p
p : Presión, lpca .
T : Temperatura, ºF
y el factor de corrección de cw, por solubilidad del gas Fig. 3.6 (parte inferior), está dado por7:
ew = 1 + 8.9 X 10-3 Rl'W (3 .13)
cwp
Carlos Bánzer S.
121
PROPIEDADES FÍSICAS DEL AGUA
; 40 ..- ~
-.
Jt
Cw : - _J_ _lh_
Bw o"p
s 3.6
~
IC
-!
c.,
e
..J 3.2 ....... _......,.__ ~~--+-
!
io
.,-¿ 2.8.......,;
en
LIJ
.r
CL
!u 2A,~60......._...._._._..._._.. ...................~..._ ........................ _..
100 140 180 220 · 260
TEMPERATURA. °F
5 10 '' 20 25
RAZON GAS DISUELTO-AGUA. PCN/BN
20.000, 200 < T(ºF) < 270 y O < S(grs/litro) < 200. Osif no presenta estimados de la exactitud
de sus resultados o de la ecuación.
Ejemplo 3.3:
Utilizando las correlaciones de Dodson, C.R. y Standing, M.B., y de Osif, T.L.,
determinar la compresibilidad de una salmuera conteniendo 20.000 ppm. de sólidos disueltos a
una presión de 5000 lpca. y temperatura de 200 ºF. La presión de saturación del crudo es de
3500 lpca.
Carlos Bánzer S.
122
CORRELACIONES NUMÉRICAS P.V.T.
Solución:
Correlación de Dodson, C.R. y Standing, M.B.
De la Fig. 3.6 ó Ec. (3.12):
3.1846- (0.00814)(200) + (3.5 X ] 0-s )(200)2 _6 _1
c.,P = = 2.96 x 10 lpc
106
De la Fig. 3.6 6 Ec. (3.13), el factor de corrección por solubilidad (Rsw = 17.8 PCN/BN,
Ejemplo 3.1) es:
e
_w =1+(8.9xl0-3)(17.8)=1.16
Cwp
luego;
Cw = (2.96 X 10-6Xl.)6) = 3.43 X 10-6 Jpc-l
De la Ec. (3 .14 ), el factor de corrección por salinidad es:
c., =1 +2º·7[-52 X 10-2 +5.4 X 10-2 -456 X 10-2 +8.97 X 10-3] =0.94
cwp
Finalmente;
Cw = (3.43 X 10-6 )(0.94) = 3.22 X 10-6 lpc -I
Correlación de Oslf, T.L.
De la Ec. (3.15):
1
e=~~------------------------------------~
" (7. 033)( 5000) + ( 541.5)(20000 / 58443) - ( 53 7)(200) + 403300
Cw = 3.02 X 10-6 lpc-1
Un resumen de los valores de cw es como sigue:
Correlación Cw, 1 pe -1
Dodson, C.R. y Standing, M.B. 3.22 X 10-6
Osif, T.L. 3.02 X ] 0-6
Para presiones menores que la presión del punto de burbujeo, la compresibilidad del
agua está definida por":
1 8R..,·w
e=-- ................................... (3.16)
w B
w qJ T
Carlos Bánzer S.
123
PROPIEDADES FISICAS DEL AGUA
.....
m
'z
~
~
EL r1tlfOlf PlfOIIEDIO NltA '' ., A ,so•, ES
1.5 % • EL EltltOlt M•XIMO E :t 15 %
o __ ......
O 1000 1000 4000 50-,00.-
P RESION, lpca
Fig. 3.7.- Cambio del gas en solución en agua de formación
con presión vs. presión (Ref. 9).
Un ajuste de la curva de la Fig. 3.7 está dado por2:
éJR.f••
q, T = B + 2Cp (3 .17)
Esta ecuación se obtuvo diferenciando la Ec. (3.1). Es de notar que esta ecuación no
ajusta muy bien la curva de la Fig. 3.7. Una expresión similar para (8R.."w / 4')7. se obtiene
diferenciando la Ec. (3.3).
Ejemplo 3.4:
Utilizando la correlación de Ramey, H.J., Jr. determinar la compresibilidad de una
salmuera conteniendo 20.000 ppm. de sólidos disueltos a una presión de 3000 lpca. y
temperatura de 200 ºF.
Solución:
De la Fig. 3 .6 6 Ec. (3 .12); la compresibilidad del agua pura a 3000 lpca. y 200 ºF e·s:
3.4526- (0.00909)(200) + (3.6627 X } 0-s )(200)2
c..,P = 106
c.,, = 3.10 X 10-6 lpc ..:.i
Cario• Bánzer S.
124
CORRELACIONES NUMÉRICAS P.V.T.
De la Fig. 3 .1 ó Ec. (3 .1 ); la solubilidad del gas natural en agua pura a 3000 lpca. y 200 ºF
es:
Rswp = 1.84729 + (5.08 X 10-3 )(3000) - (2.17767 X 10-7 )(3000)2
R.,wp = 15.13 PCN/BN
.
De la Fig. 3.2 ó Ec. (3.2), el factor de corrección para una salinidad de 20.000 ppm. a 200 ºF
es:
'
= 0.92
«;
Rsw
luego;
Cw = (3.10 X 10-6)(1.]2) = 3.47 X 10-6 lpc "
De la Ec. (3 .14 ), el factor de corrección por salinidad es:
Cw = 0.94
cwp
Finalmente,
Cw = (3.4 7 X 10-6 )(0.94) = 3.26 X 10-6 lpc -I
De las Figs. 3 .4 y 3 .5 ó Ecs. (3 .6) y (3. 7):
~Vw1· = 0.03870
~ vwp = -0. 00459
De la Ec. (3.5):
B; = (1 + 0.03870)(1-0.00459) = 1.03393 BY/BN a 3000 lpca. y 200 ºF
= 3000 = 4.47
r; 671.57
De la Fig. 1.4: z = 0.89
T = 200+460 = l.?9
sr 367.79
De la Ec. ( 1.69):
Carlos Bánzer S.
125
PROPIEDADES FiSICAS DEL AGUA
(0.89)(660)
B g = 0.00503 = 0.00098 BY /PCN
3000 •
De la Ec. (3.16):
•
A.
u
.,¡ 0.75
:::,
c.,
e ---
~
LIJ
o 0.50
o
et
-8
o
CI)
0.25
,,,,
-> .
•
o 2!500 3500
o 500 - 1000 1500 2000
PRESION, lpca
Fig. 3.8.- Comportamiento tfpico de µw vs. presión a temperatura
constante.
Carlos Bánzer S.
126
CORRELACIONES NUMÉRICAS P.V.T.
1.8
\
Salmuero, 60000 ppm
\ -<>--
o •a
o 147 lpca
a 14. 2 lpco
1.6 \ o 7100 lpca
\ 6
\ a la presión de vapor
~ ...
~
1.4
atJ.
.
u
C(
2-
o 1.2 1
e:(
z- ~
\
o \
o 1.0 fj ..,
<(
e
-
(1)
o
'o
u
0.8
• A\
-
CI)
> ~~~~
• \
0.6 .\
~
0.4 r-----'1-----+---~· -
02r-------t------+------+------+------+----- --.
O~----~----~----~~--.....a.-------'------~-----~
O 50 100 150 200 250 350 30(),
TEMPERATURA, ºF
Fig. 3.9.- Viscosidad del agua en función de temperatura para
varias presiones (Ref. 10).
La siguiente ecuación es una aproximación promedio a las curvas de la Fig. 3.911•
µw = exp(l.003 - 1.4 79 x 10-2 T + 1.982 x 1 o' T2) •••••••••••••••••••••• (3 .18)
donde, µw en cp. y Ten ºF.
Carlos Bánzer S.
127
PROPIEDADES FISICAS DEL AGUA
Temperatura, ºF µw f
40 - 120 1
-5
1.8
120 - 212 s 5 ~ 1.10
212 - 400 10 ~
~ l08 .........._...__.......
a.
-«
u
el
1.6 8 l(~ t---+---kl--+·~-1
"'eno
~
la. o: 1.04 t----+-_,_,t.....,'--1--
~
~ 1.4 ~ 1.02 .,.._-,..~..-=.-,
~ ~
0.4
TEMPERATURA, ºF
Fig. 3.10.- Viscosidad del agua a una atmósfera en función de tempe·ratura y salinidad,
para otras presiones seguir la corrección indicada (Ref. 12).
La siguiente ecuación es un ajuste de las curvas de la Fig, 3 .1013:
B
µw1 = A+ T (3.19)
donde;
Carlos Sánzer S.
128
CORRELACIONES NUMÉRICAS P.V.T.
A =
-0.04518+0.009313 S-0.000393 S2
B = 70.634 + 0.09576 S2
T : Temperatura, ºF
S : Salinidad,% por peso de sólidos disueltos (1% = 10.000 ppm).
y el factor de corrección por presión, está dado por!':
f = 1 + 3.5 X 10-12 p2 (T-40) (3.20)
donde, p en lpca. y Ten ºF. Estas correlaciones fueron desarrolladas para ell siguiente rango
de datos: p(lpca.) s 10.000, 60 ~ T(ºF) S400 y S(%) ~ 26.
Correlación de McCain, W.D., Jr.2:
Esta correlación es también un ajuste de las curvas de la Fig. 3.1 O, y está dada por:
B
µwi=AT (3.21)
donde;
A= 109.574-8.40564 s+0.313314 s2 +8.12213 x 10-3s3
B=-1.12166+2.63951 x 10-2s-6.7946I x 10-4s2-s.41119 x 10-5s3 +
isssas x 1 o-6 S4
T : Temperatura, ºF
S : Salinidad,% por peso de sólidos disueltos (1% = 10.000 ppm).
Esta correlación ajusta las curvas de la Fig. 3 .1 O dentro de un 5% para temperaturas
entre 100 y 400 ºF y salinidades hasta 26%. La Fig. 3 .11 presenta el factor de corrección por
presión 14•
1.0
o 2000 4000 6000 8000 10000 12000
PRES ION, lpca
Fig. 3.11.- Efecto de la presión sobre la viscosidad del agua (Ref. 14).
Un ajuste de la curva de la Fig. 3 .11, está dado por2:
Carlos Bánzer S.
129
PROPIEDADES FISICAS DEL AGUA
........................................... (3.23)
µw = 1-1.87 X 10-3 s" + 2.18 X 10-4 s2·5 + (Tº·5 -1.35 X 10-2 T)(2.76 X 10-3 S-3.44 X 10-4 S1'5)
µ~ .
. (3.24)
donde, Ten ºF y S (salinidad) en porcentaje por peso de sólidos disueltos (lo/o = 10.000
ppm).
Ejemplo 3.5:
Utilizando las correlaciones de Van Wingen, N., Matthews, C.S. y Russel, D.G.,
McCain, W.D., Jr., y de McCoy, R.L. determinar la viscosidad de una salmuera conteniendo
20.000 ppm. de sólidos disueltos a una presión de 5000 lpca. y temperatura de 200 ºF.
Solución:
Correlación de Van Wlngen, N.
De la Fig. 3.9 ó Ec. (3.18):
µw =exp[l.003-(1.479 X 10-2)(200)+(1.982 x I0-5)(200)2]
luego;
µw=0.313cp.
Correlación de Matthews, C.S. y Russel, D.G.
De la Fig, 3 .1 O ó Ec. (3 .19); la viscosidad de una salmuera conteniendo 2% de sólidos
disueltos a 1 atm. y 200 ºFes:
µ...,1 = 0.327 cp.
Del cuadro inserto en la Fig. 3.10 ó Ec. (3.20), la corrección de µ..,1 por presión es:
j = 1 + (3.5 X 10-12 )(5000)(200- 40) = 1.0140
Finalmente;
µw = (0.327)(1.0140) = 0.332 cp
Correlación de McCain, W.D., Jr.
De la Fig. 3.10 6 Ec. (3.21); la viscosidad de una salmuera conteniendo 2% de sólidos
disueltos a 1 atm. y 200 ºFes:
µwl = (94.086)(200)-1·º12 = 0.32 cp .
•
Carlos Bánzer S.
130
CORRELACIONES NUMÉRICAS P.V.T.
De la Ec. (3.24), la corrección por salinidad para una salmuera conteniendo 2% de sólidos
disueltos a 200 ºF es:
µw = 1.0506
µwp
Finalmente;
µw = (0.299)(1.0506) = 0.314 cp
Un resumen de los resultados deµ..., es como sigue:
Correlación
donde;
p; : Densidad del agua a p y T, lbs/pie3
Pwi Densidad del agua a condiciones normales ( 14. 7 lpca. y 60 ºF), lbs/pie'
:
rw Gravedad específica del agua, adim,
:
B; : Factor volumétrico del agua ap y T, BY/BN
La gravedad específica del agua de formación, r w , puede ser estimada mediante la
siguiente ecuación si se conoce la concentración de sólidos disueltos":
y w = 1.0 + 0.695 X ) 0-6 S (3 .26)
donde, Ses la concentración de sólidos disueltos (salinidad) en mg/litro.
Carlos Bánzer S.
131
PROPIEDADES FÍSICAS DEL AGUA
.-
~
72
~
<[..,
o:: uJ
1&.1- 70
::> a..
2~
~m 68
e:( ..J
(/) ..
et u,
..J o 16
o
~ co
e 64
s-
V) 62
z
L&J
Q
60
o 10 15 20 25 30
TOTAL SOLIDOS DISUELTOS1 º/o
Fig. 3.12.- Densidad de la salmuera en función del porcentaje total
de sólidos disueltos (Ref. 2).
La siguiente ecuación ajusta la curva de la Fig. 3 .12 dentro de un 0.1 %2:
p wt = 62.368 + 0.438603 S + ).60074 X 10-3 S2 •••••••••••••••••••••••••• (3 .27)
donde, Ses la salinidad en porcentaje por peso de sólidos disueltos ( 1 % = 10.000 ppm ).
Ejemplo 3.6:
Utilizando la correlación de McCain, W.D., Jr. determinar la densidad de una
salmuera conteniendo 20.000 ppm. de sólidos disueltos a una presión de 50f)0 lpca, y
temperatura de 200 ºF.
Solución:
De la Fig. 3.12 6 Ec. (3.27); la densidad de una salmuera conteniendo 20.00() ppm. de
sólidos disueltos a 14.7 lpca. y 60 ºF, es:
p wl = 62.368 + (0.438603)(2) + (].60074 X 10-3 )(2)2
Carlos Bánzer S.
132
CORRELACIONES NUMÉRICAS P.V.T.
eu re
'g.. 65
·-
.
"O
....
~
~
<[
60
,.~
t
"A)~,,...
~ !5!5 ~
(I~
~, o,-
-'
~ 50
c.) ~
74
~
et:
k.l 45
t-
z-
!-
U)
40
212
•
z ......
uJ ~
1- 35 ~
-
lo...
3 O
30 'et
O 2000 4000 6000 8000 10000 12000
PRESION, lpca
Fig. 3.13 .- Tensión interfacial gas-agua (Ref. 15).
Un ajuste de las curvas de la Fig. 3.13, está dado por2:
ugw = A + Bp + Cp 2
(3 .2·8)
donde;
A = 79.1618-0.118978 T
B = -5.28473 x 10-3 +9.87913 x 10-6r
C = (2.33814-4.57194 x 10-4T-7.52678 x 10-6T2) x 10-1
p : Presión, lpca.
T : Temperatura, ºF
Es recomendable no utilizar esta ecuación a presiones mayores de 8000 lpca., y
temperaturas mayores de 350 ºF2.
Carlos Bánzer S.
133
PROPIEDADES
•
FÍSICAS DEL AGUA
Ejemplo 3.7:
Utilizando la correlación de Jennings, H.Y., Jr. y Newman, G.H. determinar la tensión
interfacial gas-agua a una presión de 5000 lpca. y temperatura de 200 ºF.
Solución:
De la Fig. 3 .13 ó Ec. (3 .28):
O"gw= 55.3662-(0.00331)(5000) + (1.94563 X 10-7 )(.5000)2
ugw = 43.68 dinas/cm.
Referencias:
1. Culberson, O.L. and McKetta, J.J., Jr.: ''Phase Equilibria in Hydrocarbon-Water
Systems 111-The Solubility of Methane in Water at Pressures to 10,000 psia,"
Trans., AIME (1951) 192, 223-226.
2. McCain, W.D., Jr.: The Properties o/ Petroleum Fluids, 2nd. ed. Tulsa, OK.: PennWell
Books, 1990.
3. McKetta, J.J. and Wehe, A.H.: "Hydrocarbon-Water and Formaiion Water
Correlations," Petroleum Production Handbook, Vol. 11, T.C. Frick and R.W.
Taylor (eds.), SPE, Dallas (1962) 22-13.
4. McCoy, R.L.: Microcomputer Programs for Petroleum Engineers: Vol. 1,
Reservoir Engineering and Formation Evaluation, Gulf Publishing Co., Houston
(1983).
5. Numbere, D., Brigham, W.E., and Standing, M.B.: Correlations for Physical
Properties of Petroleum Reservoir Brines, Stanford University Petroleum Research
Institute, SUPRI-TR-1, Nov. 1977.
6. Dodson, C.R. and Standing, M.B.: "Pressure-Volume-Ternperature and Solubility
Relations for Natural Gas-Water-Mixtures," Drill. and Prod. Prac., API (1944)
173-179. ·
7. Meehan, D.N.: ''A Correlation for Water Compressibility," Petroleum Engineer
(Nov. 1980) 125-126.
8. Osif, T.L.: ''The Effects of Salt, Gas, Temperature, and Pressure on the
Compressibility of Water," SPE Res. Eng. (Feb. 1988) 3, No. r, 175-181.
9. Ramey, H.J., Jr.: ''Rapid Methods for Estimating Reservoir Compressibilities,"
Trans., AIME (1964) 231, 447-454.
10. Van Wingen, N.: "Viscosity of Air, Water, Natural Gas, and Crude Oil at Varying
Pressure and Temperatures," Secondary Recovery of Oil in the United States, API
(1950) 127.
11. Brill, J.P. and Beggs, H.D.: ''Two-Phase Flow in Pipes," University of Tulsa,
Tulsa, OK, 1988.
12. Matthews, C.S. and Russel, D.G.: Pressure Buildup and Flow Test in Wells,
Monograph Series, Society of Petroleum Engineers of AIME, Dallas ( 1967) 1,
Appendix G.
13. Meehan, D.N.: ''A Correlation for Water Viscosity," Petroleum Engineer (July
1980) 117-118.
Carlos Bánzer S.
134
CORRELACIONES NUMéRICAS P.V.T.
Carlos Bánzer S.
135
UNIDADES Y FACTORES DE· CONVERSIÓN
Compresibilidad lpc-1 X --
1.450 377 E - 01 kPa-1
Densidad lbs/pie3 X 1.601 846 E+ 01 - kg/m3
Factor volumétrico del gas BY/PCN X 5.551 931 E+ 00 - m3/m3
Factor volumétrico del petróleo BY/BN X 1.0E + 00 ª -- m3/m3
Gravedad del petróleo ºAPI 141.5/(131.S+y API) -b grs/cm3
Presión lpca. X 6.894 757 E+ 00 - kPa
Solúbilidad PCN/BN X 1.801 175 E - 01 - m3/m3
Temperatura (ambiente) ºF (ºF - 32)/1.8 b - ºC
Temperatura (absoluta) ºR ºR/1.8 b - K
Tensión interfacial dinas/cm X 1.0E + 00 ª - mN/m
Viscosidad dinámica cp. X 1.0E-03 ª - Pa.s
Viscosidad cinemática cts. X 1.0E-06ª -- m2/seg
Volumen pie3 X 2.831 685 E- 02 -- m3
bis X 1.589 873 E - 01 - m3
gal. X 3.785 412E- 03 -- m3
"lndica que el factor de conversión es exacto.
'No es un factor de conversión (efectuar la operación indicada).
* Adaptada de: Prats, M.: Thermal Recovery, Monograph Series, SPE of AIME, Dallas ( 1982) 7.
Ejemplo 4.1 :
La correlación de Standing, M.B. para estimar la presión del punto de burbujeo está dada
por:
0.83
p,, = 18.2 Rsh 1 o<0.000917'-0.012Sr .41'/) - 1.4
Yg
donde, Pb en lpca., Rsb en PCN/BN, Ten ºF, r APJ en º API y y g (adim.). Calcular· la presión
de burbujeo en kPa para los siguientes datos en unidades SI, primero convirtiendo las unidades
SI a unidades de campo y luego convirtiendo el resultado a unidades SI. Posteriormente,
derivar las constantes de la correlación para unidades SI.
s; = 122 m3/m3 = 675 PCN/BN
T = 356 K = 180 ºF
y,) = 0.871 grs/cm3 = 31 ºAPI
rg = 0.95 (aire= 1)
Carlos Bánzer S.
136
CORRELACIONES NUMéRICAS P.V.T.
Solución:
Sustituyendo los datos en unidades de campo en la correlación, se tiene:
0.83
675 } Q(0.00091 X 180-0.012S X 31) _ ) .4
Pb = 18.2
0.95
luego;
P» = 2505 lpca. = 17271 kPa
Derivando las constantes de la correlación para unidades SI, se tiene:
0~ 3 3
Rs,, (PCN I BN) (m /m)
( 1 ) (kPa) - 18 2
P» pea. x 6.894757 (lpca.) - · X-----------
Yg ).801175 x l~)-t (PCN/BN)
luego;
0.83
Rsb to<0.00164T-l.76875/y u+l.2251S) - 1.4
Pb = (6 ·894757)(18 ·2) 4 · 1 S
r,
Finalmente;
0.83
Rsb 1o<t.22SIS+0.001647'-1.7687S/y,,) - 0.34
Pb = 520 · 76
r,
donde.p, en kPa., R31, en m3/m3, Ten K, r O en grs/cm3 y r 1 (adim.).
Sustituyendo los datos en unidades SI, se tiene:
0.83 .
122 ) Q(l.22SIS+0.00164 x 3S6-l.7687S/0.871) _ 0.34
Pb =520.76
0.95
luego;
Pb = 17406 kPa
Carlos Sánzer S.
137
GLOSARIO DE TÉRMINOS Y DEFINICIONES
Carlos Sánzer S.
138
CORRELACIONES NUMÉRICAS P.V.T.
Carlos Bánzer S.
139
GLOSARIO DE TÉRMINOS Y DEFINICIONES