El yacimiento de Gas seco (G-01) presenta el siguiente análisis cromatográfico:

Componentes % Ni Yi Mi Pc Tc Psc Tsc

C1 90,32 0,903 16,043 667,8 343,37 603,023 310,063
C2 3,94 0,039 30,070 707,8 550,09 27,604 21,454
C3 2,6 0,026 44,097 616,3 666,01 16,024 17,316
n-C4 1,52 0,015 58,124 550,7 765,65 8,261 11,485
n-C5 0,52 0,005 72,151 488,6 845,70 2,443 4,229
C7+ 0,84 0,008 110 371,522 1008,706 2,972 8,070
H2S 0,11 0,001 34,076 1306 672,70 1,306 0,673
CO2 0,15 0,002 44,010 1071 547,90 2,142 1,096
∑ 100 1 663,775 374,386

%Ni
Yi=
100

Psc=Yi∗Pc

Tsc=Yi∗Tc

Adicional el Yacimiento presenta los siguientes datos de producción Ty=720 oR, γ g =0,65

Py(Lpca) 2900 2800 2700 2650 2600 2500 2000 1500 1000 500

 Evaluar el Factor Volumétrico del Gas (Bg) en función de la presión para este
Gas seco.

( )
0,02829∗Z∗T y Pcy
Bg =
Py Pcn

Corrección de impurezas por Carr, Kobayashi y Burrows

'
Ps c =Psc−166 Y N 2+ 440 Y CO 2+ 600Y H 2 S
'
Ps c =663,775+ 440 ( 0,002 ) +600 ( 0,001 )
'
Ps c =665,255

'
Ts c =Tsc−250 Y N 2−83,3 Y CO2 +130 Y H 2 S
'
Ts c =374,386−83,3 ( 0,002 ) +130 ( 0,001 )
'
Ts c =374,349
La Temperatura permanece constante, de tal manera que solo evaluaremos el Bg en
función de los cambios en la presión.

Ty 720 R
Tsr= Tsr= Tsr=1,92
Tsc
'
374,349

Bg a 2900 Lpca
Py 2900 Lpca
Psr= Psr = Psr =4,359
Ps c
'
665,255

Tsr=1,92 Psr =4,359

Z=0,92

0,02829∗0,92∗720 R PCY
Bg = =0,006
2900 Lpca PCN

Bg a 2800 Lpca

Tsr=1,92 Psr =4,209

Z=0,925

0,02829∗0,925∗720 R PCY
Bg = =0,007
2800 Lpca PCN

Bg a 2700 Lpca

Tsr=1,92 Psr =4,059

Z=0,922
0,02829∗0,922∗720 R PCY
Bg = =0,0070
2700 Lpca PCN

Bg a 2650 Lpca

Tsr=1,92 Psr =3.983

Z=0,918
0,02829∗0,918∗720 R PCY
Bg = =0,0071
2650 Lpca PCN

Bg a 2600 Lpca

Tsr=1,92 Psr =3,908

Z=0,92
 0,02829∗0,92∗720 R PCY
Bg = =0,0072
2600 Lpca PCN

Bg a 2500 Lpca

Tsr=1,92 Psr =3,75

Z=0,93

0,02829∗0,93∗720 R PCY
Bg = =0,008
2500 Lpca PCN

Bg a 2000 Lpca

Tsr=1,92 Psr =3,006

Z=0,93

0,02829∗0,93∗720 R PCY
Bg = =0,009
2000 Lpca PCN

Bg a 1500 Lpca

Tsr=1,92 Psr =2,25

Z=0,94

0,02829∗0,94∗720 R PCY
Bg = =0,013
1500 Lpca PCN

Bg a 1000 Lpca

Tsr=1,92 Psr =1,50

Z=0,95

0,02829∗0,95∗720 R PCY
Bg = =0,019
1000 Lpca PCN

Bg a 500 Lpca

Tsr=1,92 Psr =0,75

Z=0,97

0,02829∗0,97∗720 R PCY
Bg = =0,040
500 Lpca PCN

P(Lpca) 2900 2800 2700 2650 2600 2500 2000 1500 1000 500
Bg(PCY/PCN) 0,00 0,007 0,0070 0,0071 0,0072 0,008 0,00 0,013 0,019 0,040
6 9

Análisis del comportamiento de la Gráfica presentada

En la gráfica podemos notar que a medida que la presión del yacimiento va

declinando debido a la producción, el factor volumétrico del gas va aumentando de
manera exponencial. Podemos notar que los valores de los factores volumétricos del
gas eran muy similares al inicio, posterior a esto ocurre una rápida declinación de
presión ocasionando mayor diferencia en los valores del factor volumétrico. A altas
presiones se requiere de un valor muy bajo del factor volumétrico del gas a
condiciones de yacimiento para producir un Ft3 de gas a condiciones normales, esto
se debe a que el gas se expande ocupando un volumen mayor al volumen ocupado
en el yacimiento. A medida que la presión va declinando requerimos de un valor
mayor del factor volumétrico del gas a condiciones de yacimiento para producir un
Ft3 de gas a condiciones normales.

 Determinar la compresibilidad del gas a través de dos métodos

Método gráfico (Trube)

 Csr
C g=
Psc '

Dónde:

C g=Compresibilidad del Gas Csr=Compresibilidad Seudorreducida

Psc=Presion Seudocrítica

La Compresibilidad Seudorreducida para gases naturales, la podemos obtener

mediante la gráfica de Trube.

A 2900 Lpca

Tsr=1,92 Psr =¿ 4,359

Csr=0,21

0,21 −4 −1
C g= =3,157 x 10 Lpc
665,255

A 2800 Lpca

Tsr=1,92 Psr =4,209

Csr=0,22

0,22 −4 −1
C g= =3,307 x 10 Lpc
665,255

A 2700 Lpca

Tsr=1,92 Psr =4,059

Csr=0,23

0,23 −4 −1
C g= =3,457 x 10 Lpc
665,255

A 2650 Lpca

Tsr=1,92 Psr =3.983

Csr=0,24

0,24 −4 −1
C g= =3,608 x 10 Lpc
665,255
A 2600 Lpca

Tsr=1,92 Psr =3,908

Csr=0,245

0,245 −4 −1
C g= =3,683 x 10 Lpc
665,255

A 2500 Lpca

Tsr=1,92 Psr =3,75

Csr=0,26

0,26 −4 −1
C g= =3,913 x 10 Lpc
665,255

A 2000 Lpca

Tsr=1,92 Psr =3,006

Csr=0,35

0,35 −4 −1
C g= =5,261 x 10 Lpc
665,255

A 1500 Lpca

Tsr=1,92 Psr =2,25

Csr=0,44

0,44
C g= =6,614 x 10−4 Lpc−1
665,255

A 1000 Lpca

Tsr=1,92 Psr =1,50

Csr=0,76

0,76
C g= =11,424 x 10−4 Lpc−1
665,255

A 500 Lpca

Tsr=1,92 Psr =0,75

Csr=1,06
 1,06 −4 −1
C g= =15,934 x 10 Lpc
665,255

Método De Matar, Brar y Aziz

Psc=663,775 Tsc=374,386
Yi CO 2=0,15Yi H 2 S=0,11

Corrección de impurezas por Whichert y Aziz

'
Ts c =Tsc−Fsk

' Psc∗Tsc '

Ps c =
Tsc +B ( 1−B )∗Fsk

Fsk=120 ( A 0,9− A1,8 )+ 15(√ B−B4 )

Donde :
A=Suma de las Fracciones molares del CO 2 y H 2 S
B=Fraccion molar del H 2 S

Corrigiendo Impurezas

Fsk=120 ( 0,26 −0,26 ) +15 ( √ 0,11−0,11 )

0,9 1,8 4

Fsk=30,052
Ts c ' =374,386−30,052
'
Ts c =344,334
' 663,775∗344,334
Ps c =
374,386 +0,11 ( 1−0,11 )∗30,052
'
Ps c =605,734
La Compresibilidad del gas podemos obtenerla mediante la ecuación:

Csr
C g=
Psc '

A diferencia del método de Trube el cálculo del parámetro Csr se obtiene mediante
la siguiente ecuación:
Csr=
1 0,27
− 2
Psr Z ∗Tsr
[ ] ( ∂∂ Zρ )
1+ (
Z ∂ρ)
ρ ∂Z
r
Tsr

Tsr

El parámetro ( ∂∂Zρ )
Tsr
se obtiene:

( ) ∂Z 2

=T 1+ 2T 2 ρ r +5T 3 ρ4r +2 T 4 ρr (1+ A 8 ρ2r − A 28 ρ4r )e (−A 8 ρr )

∂ρ Tsr

Dónde:

A8=0,68446549

T1 = -0,249

T2=0,243

T3=0,031

T4=0,075

ρr =Densidad seudorreducida

0,27 Psr
ρr =
Z Tsr

Una vez planteadas las ecuaciones y corregidas las impurezas debemos calcular las
presiones seudorreducidas Psr para cada variación de presión que tengamos, Con
las condiciones seudorreducidas determinados el factor de compresibilidad Z, de
igual manera las densidades seudorreducidas ρr y las compresibilidades
seudorreducidas Csr variaran con los cambios en las presiones seudorreducidas.

La temperatura seudorreducida será constante

Ty 720
Tsr= Tsr= Tsr=2,091
Tsc
'
344,334

A continuación una muestra de cálculos

 Cuando la P es 2900 lpca

Py 2900
Psr= Psr = Psr=4,788
Ps c
'
605,734

Z=0,95

0,27∗4,788
ρr = ρ =0,651
0,95∗2,091 r

( ∂∂ Zρ )
Tsr
=−0,249+ 2∗0,243∗0,651+5∗0,031∗0,651 + 2∗0,075∗0,651 ( 1+ 0,68446549∗0,651 −0,68446549 ∗
4 2 2

( ∂∂ Zρ )
Tsr
=¿ 0,647

Csr=

Csr=0,145
1
−
0,27
Psr Z 2∗Tsr
[ ]
( ∂∂ Zρ )
1+ (
Z ∂ρ)
ρ ∂Z
r
Tsr

Tsr

Repitiendo el mismo procedimiento a las condiciones de Presión variable y

Temperatura constante del yacimiento tenemos:

P(Lpca) Psr Z ρr
( ∂∂Zρ )
Tsr
Csr

2900 4,788 0,95 0,651 0,647 0,145

2800 4,632 0,96 0,621 0,645 0,152
2700 4,457 0,94 0,612 0,639 0,158
2650 4,375 0,94 0,601 0,632 0,163
2600 4,292 0,95 0,583 0,619 0,169
2500 4,127 0,93 0,573 0,613 0,176
 2000 3,302 0,945 0,451 0,534 0,241
1500 2,476 0,96 0,333 0,458 0,349
1000 1,651 0,955 0,223 0,390 0,353
500 0,875 0,97 0,110 0,319 0,362

Finalmente con los valores de la Csr podemos calcular la compresibilidad del gas

A 2900 Lpca
Csr=0,145

0,145 −4 −1
C g= =2,0 x 10 Lpc
605,734

A 2800 Lpca

Csr=0,152

0,152
C g= =3,0 x 10−4 Lpc−1
605,734

A 2700 Lpca

Csr=0,158

0,158 −4 −1
C g= =3,1 x 10 Lpc
605,734

A 2650 Lpca

Csr=0,163

0,163
C g= =3,9 x 10−4 Lpc−1
605,734
A 2600 Lpca

Csr=0,169

0,169 −4 −1
C g= =3,8 x 10 Lpc
605,734

A 2500 Lpca
Csr=0,176

0,176 −4 −1
C g= =3,9 x 10 Lpc
605,734

A 2000 Lpca

Csr=0,241

0,241
C g= =4,1 x 10−4 Lpc−1
605,734

A 1500 Lpca
Csr=0,349

0,349 −4 −1
C g= =5,8 x 10 Lpc
605,734

A 1000 Lpca

Csr=0,353

0,353 −4 −1
C g= =5,83 x 10 Lpc
605,734

A 500 Lpca
Csr=0,362

0,362
C g= =5,9 x 10−4 Lpc−1
605,734
 REPÚBLICA BOLIVARIANA DE VENEZUELA

UNIVERSIDAD DE ORIENTE

NÚCLEO DE ANZOÁTEGUI

ESCUELA DE INGENIERÍA Y CIENCIAS APLICADAS

DEPARTAMENTO DE PETRÓLEO
 Asignación N0 1. Evaluación del Factor Volumétrico y
Compresibilidad del Gas.

Bachilleres: Profesor:

Barrios José Luis C.I 24225692 Ing. Jairo Uricare

Gómez Kendis C.I 25426734

Salazar Noriurka C.I 25844025

Barcelona, Julio 2016