PROPIEDADES FISICAS

Propiedades físicas y termodinámicas

Maestría
Ingeniería de Procesos, Transporte, Almacenamiento y Distribución de Gas Natural

Jhordana Valeria Guzmán Fernández

Tarea # 2

Propiedades Físicas y Termodinámicas

Asignación # 1

Desarrollo:

a) Para calcular el peso molecular aparente, en Fig. 23.2 (Anexo 1), leer Peso Molecular de
cada componente presente y posteriormente multiplicarlo por la fracción molar,
obtenida de dividir el porcentaje molar entre 100. A continuación, sumar todos los
porcentajes molares obtenidos para obtener el peso molecular aparente de la mezcla.

Comp % molar yi Pmi Pma

C1 76.25 0.7625 16.042 12.232025

C2 7.42 0.0742 30.069 2.2311198
C3 5.5 0.055 44.096 2.42528
iC4 0.78 0.0078 58.122 0.4533516
nC4 0.94 0.0094 58.122 0.5463468
iC5 0.37 0.0037 72.149 0.2669513
nC5 0.41 0.0041 72.149 0.2958109
C6+ 0.67 0.0067 110 0.737
CO2 1.63 0.0163 44.01 0.717363
H2S 5.51 0.0551 34.082 1.8779182
N2 0.52 0.0052 28.0135 0.1456702
100 21.9288368

𝑷𝑴𝑨 = 𝟐𝟏. 𝟗𝟐𝟖𝟖

 b) Para el cálculo de Presión y Temperatura seudocríticas, en Fig. 23-2 (Anexo 1), leer la
Presión Critica (psia) y Temperatura Critica (ªF). Convertir la Temperatura Crítica en (ºR),
sumando 459.67. Multiplicar ambas por la fracción molar, obteniendo las columnas de
Presión y Temperatura seudocríticas. Para el cálculo de éstas propiedades de la mezcla,
sumar las mismas.

Pci Tci Tci Psc Tsc

Comp % molar yi Pmi Pma
(psia) (º F) (º R) (psia) (ºR)

C1 76.25 0.7625 16.042 12.232025 667 -116.66 343.01 508.5875 261.545125

C2 7.42 0.0742 30.069 2.2311198 706.6 89.92 549.59 52.42972 40.779578
C3 5.5 0.055 44.096 2.42528 615.5 205.92 665.59 33.8525 36.60745
iC4 0.78 0.0078 58.122 0.4533516 527.9 274.41 734.08 4.11762 5.725824
nC4 0.94 0.0094 58.122 0.5463468 550.9 305.55 765.22 5.17846 7.193068
iC5 0.37 0.0037 72.149 0.2669513 490.4 369 828.67 1.81448 3.066079
nC5 0.41 0.0041 72.149 0.2958109 488.8 385.8 845.47 2.00408 3.466427
C6+ 0.67 0.0067 110 0.737 418.3 576.6 1036.23 2.80295 6.942753
CO2 1.63 0.0163 44.01 0.717363 1070 87.76 547.43 17.441 8.923109
H2S 5.51 0.0551 34.082 1.8779182 1306.5 212.81 672.48 71.98815 37.053648
N2 0.52 0.0052 28.0135 0.1456702 492.5 -232.53 227.14 2.561 1.181128
100 21.9288368 702.777 412.484

En el caso de C6+, realizamos los cálculos con correlaciones numéricas con gravedad
especifica

𝑃𝑆𝐶 = 1188 − 431 log(𝑀𝐶7+ − 61.1) + [2319 − 852 × log(𝑀𝐶7+ − 53.71)](𝛾𝐶7+ − 0.8)

𝑇𝑆𝐶 = 608 + 364 log(𝑀𝐶7+ − 71.2) + [2450 log(𝑀𝐶7+ ) − 3800] × log⁡(𝛾𝐶7+ )

𝑷𝑺𝑪 ⁡(𝒑𝒔𝒊𝒂) = 𝟒𝟏𝟖. 𝟑

𝑻𝑺𝑪 ⁡(º𝑹) = 𝟏𝟎𝟑𝟔. 𝟐𝟑

c) Para el cálculo de Valor Calorífico Neto y Bruto, a partir de los componentes se leen los
mismos valores de la Fig. 23.2 (Anexo 3). Para componente C6+ se procede de igual
forma con la interpolación. Cuando tenemos los datos de todos, éstos se multiplican por
la fracción molar y se realiza la sumatoria.

PC NETO PC BRUTO PC NETO PC BRUTO

Comp yi Pmi
(Btu/ft3) (Btu/ft3) (Btu/ft3) (Btu/ft3)

C1 0.7625 16.042 909.4 1010 693.4175 770.125

C2 0.0742 30.069 1618.7 1769.7 120.10754 131.31174
C3 0.055 44.096 2314.9 251602 127.3195 13838.11
iC4 0.0078 58.122 3000.4 3252 23.40312 25.3656
 nC4 0.0094 58.122 3010.8 3262.4 28.30152 30.66656
iC5 0.0037 72.149 3699 4000.9 13.6863 14.80333
nC5 0.0041 72.149 3706.9 1008.7 15.19829 4.13567
C6+ 0.0067 110 5586.2 6023.9 37.42754 40.36013
CO2 0.0163 44.01 0 0 0 0
H2S 0.0551 34.082 586.8 637.1 32.33268 35.10421
N2 0.0052 28.0135 0 0 0 0
1091.19399 14889.9822
𝑷𝑪𝑵𝒆𝒕𝒐 ⁡(𝑩𝒕𝒖⁄ 𝟑 ) = 𝟏𝟎𝟗𝟏. 𝟏𝟗𝟒
𝒇𝒕

𝑷𝑪𝑩𝒓𝒖𝒕𝒐 ⁡(𝑩𝒕𝒖⁄ 𝟑 ) = 𝟏𝟒𝟖𝟖𝟗. 𝟗𝟖𝟐

𝒇𝒕

d) Para el cálculo de riqueza del gas, leemos densidad molar de componentes de Fig. 23.2
(Anexo 2). Posteriormente, dividimos peso molecular entre densidad para obtener
unidades de (gal/lbmol). Como último paso, multiplicamos 1000*fracción
molar*Densidad molar(gal/lbmol) y dividimos entre 379.6, para realizar la sumatoria
total.

Densidad Densidad
Comp yi Pmi Molar Molar GPM
(lb/gal) (gal/lbmol)
C1 0.7625 16.042 2.5 6.4168 12.8893836
C2 0.0742 30.069 2.9716 10.1187912 1.97790914
C3 0.055 44.096 4.2301 10.4243399 1.5103759
iC4 0.0078 58.122 4.6934 12.383773 0.25446109
nC4 0.0094 58.122 4.8696 11.9356826 0.29556221
iC5 0.0037 72.149 5.2074 13.8550908 0.13504699
nC5 0.0041 72.149 5.2618 13.7118477 0.14809951
C6+ 0.0067 110 5.7352 19.1798019 0.33852654
CO2 0.0163 44.01 6.8534 6.42163014 0.27574439
H2S 0.0551 34.082 6.6922 5.0927946 0.73923336
N2 0.0052 28.0135 6.727 4.16433774 0.05704572
18.6213884
𝑮𝑷𝑴 = 𝟏𝟖. 𝟔𝟐𝟏𝟒
Considerando:

𝐺𝑃𝑀 ≥ 3 → 𝐺𝑎𝑠⁡𝑅𝑖𝑐𝑜
𝐺𝑃𝑀 ≤ 3 → 𝐺𝑎𝑠⁡𝑃𝑜𝑏𝑟𝑒
El gas del que se realizan los cálculos, es un Gas Rico.
e) Factor de compresibilidad a P=2100 psia y T=150ºF

i) Método de Standing Katz

Para utilizar éste método se debe calcular
𝐴 = 𝑦𝐶𝑂2 + 𝑦𝐻2 𝑆 = 0.0714
𝐵 = 𝑦𝐻2 𝑆 = 0.0551
𝐹𝑠𝑘 = 120(𝐴0.9 − 𝐴1.6 ) + 15(√𝐵 − 𝐵4 ) = 12.9186

La temperatura y presión se ajustan por medio del factor de ajuste Fsk, de la

siguiente manera:

𝑇`𝑆𝐶 = 𝑇𝑆𝐶 − 𝐹𝑠𝑘

(𝑃𝑆𝐶 )(𝑇`𝑆𝐶 )
𝑃`𝑠𝑐 =
𝑇𝑆𝐶 + (𝐵)(1 − 𝐵)(𝐹𝑠𝑘)
Con Datos de Psc y Tsc obtenidos en el inciso b, realizamos la corrección.
Obteniendo los siguientes resultados:

𝑇`𝑆𝐶 = 399,5656
𝑃`𝑠𝑐 = 679.66
Calculamos propiedades seudoreducidas dividiendo la propiedad entre la misma
propiedad seudocrítica.

Con temperatura y presión corregidas, obtenemos propiedades seudoreducidas

y con las mimas leemos z en Fig. 23.4 (Anexo 4)
𝑇
𝑇𝑠𝑟 = = 1.53
𝑇𝑆𝐶
𝑃
𝑃𝑠𝑟 = = 3.09
𝑃𝑆𝐶

Con Tsr y Psr leer z para Tsr=1.5 y Tsr=1.6 e interpolar para Tsr=1.53

𝑧𝑇𝑠𝑟=1.5 = 0.775

𝑧𝑇𝑠𝑟=1.6 = 0.820

𝒛𝑻𝒔𝒓=𝟏.𝟓𝟑 = 𝟎. 𝟕𝟗𝟎

ii) Método de Papay

Con Temperatura y presión obtenidas en inciso (b), calculamos propiedades
reducidas
𝑇
𝑇𝑠𝑟 = = 1.48
𝑇𝑆𝐶
𝑃
𝑃𝑠𝑟 = = 2.99
𝑃𝑆𝐶
 Calculamos z en la siguiente ecuación:
3.52(𝑃𝑟 ) 0.274(𝑃𝑟2 )
𝑧=1− + = ⁡0.6951
10(0.9213(𝑇𝑟)) 10(0.8157(𝑇𝑟 ))

Calculamos error:
|𝑎𝑝𝑟𝑜𝑥 − 𝑒𝑥𝑎𝑐𝑡|
%𝑒𝑟𝑟𝑜𝑟 = × 100 = 12.0%
𝑒𝑥𝑎𝑐𝑡

iii) Método de Brill y Beggs

Con Temperatura y presión obtenidas en inciso (b), calculamos propiedades
reducidas
𝑇
𝑇𝑠𝑟 = = 1.48
𝑇𝑆𝐶
𝑃
𝑃𝑠𝑟 = = 2.99
𝑃𝑆𝐶
Calcular:
𝐴 = 1.39(𝑇𝑠𝑟 − 0.92)0.5 − 0.36𝑇𝑠𝑟 − 0.10 = 0.4063

0.066 2
0.32 6
𝐵 = (0.62 − 0.23𝑇𝑠𝑟 )𝑃𝑠𝑟 + [ − 0.037] 𝑃𝑠𝑟 + 9(𝑇 −1) 𝑃𝑠𝑟 = 1.471
𝑇𝑠𝑟 − 0.86 10 𝑠𝑟

𝐶 = 0.132 − 0.32 log 𝑇𝑠𝑟 = 0.0777

2
𝐷 = 𝑎𝑛𝑡𝑖𝑙𝑜𝑔(0.3106 − 0.49𝑇𝑠𝑟 + 0.1824𝑇𝑠𝑟 ) = 0.9657

1−𝐴 𝐷
𝑧=𝐴+ + 𝐶𝑃𝑠𝑟 = 0.77
exp(𝐵)

𝒛 = 𝟎. 𝟕𝟕

|𝑎𝑝𝑟𝑜𝑥 − 𝑒𝑥𝑎𝑐𝑡|
%𝑒𝑟𝑟𝑜𝑟 = × 100 = 3.0%
𝑒𝑥𝑎𝑐𝑡

iv) Método de Pitzer

Con Temperatura y presión obtenidas en inciso (b), calculamos propiedades
reducidas
𝑇
𝑇𝑠𝑟 = = 1.48
𝑇𝑆𝐶
𝑃
𝑃𝑠𝑟 = = 2.99
𝑃𝑆𝐶
Calculamos:
0.422
𝐵𝑂 = 0.083 − 1.6 = −0.1428
𝑇𝑠𝑟
0.422
𝐵` = 0.083 − 1.6 = 0.1057
𝑇𝑟
𝑅𝑇𝐶
𝐵 = (𝐵𝑂 + 𝑤𝐵` ) = −0.8720
𝑃𝐶
𝐵𝑃
𝑧 =1+ = 0.7200
𝑅𝑇
 𝒛 = 𝟎. 𝟕𝟐

|𝑎𝑝𝑟𝑜𝑥 − 𝑒𝑥𝑎𝑐𝑡|
%𝑒𝑟𝑟𝑜𝑟 = × 100 = 8.85%
𝑒𝑥𝑎𝑐𝑡

v) Método gráfico de la GPSA

Con Temperatura y presión obtenidas en inciso (b), calculamos propiedades
reducidas
𝑇
𝑇𝑠𝑟 = = 1.48
𝑇𝑆𝐶
𝑃
𝑃𝑠𝑟 = = 2.99
𝑃𝑆𝐶

Con Propiedades seudoreducidas, en Fig. 23.4, leemos z

𝑧𝑇𝑠𝑟=1.45 = 0.745

𝑧𝑇𝑠𝑟=1.5 = 0.78

𝒛𝑻𝒔𝒓=𝟏.𝟓𝟑 = 𝟎. 𝟕𝟔𝟔

|𝑎𝑝𝑟𝑜𝑥 − 𝑒𝑥𝑎𝑐𝑡|
%𝑒𝑟𝑟𝑜𝑟 = × 100 = 5.7%
𝑒𝑥𝑎𝑐𝑡

Desarrollo

𝑃𝑉 = 𝑛𝑅𝑇⁡⁡⁡(1)
𝑚⁡(𝑙𝑏)
𝑛= ⁡⁡⁡(2)
𝑙𝑏
𝑃𝑀⁡( )
𝑙𝑏𝑚𝑜𝑙

Reemp. (2) en(1)

 𝑚
𝑃𝑉 = 𝑅𝑇
𝑃𝑀
𝑃 × 𝑉 × 𝑃𝑀
𝑚= ⁡⁡⁡(3)
𝑅𝑇
Para reemplazar datos en ecuación primero calculamos Peso Molecular de la mezcla:

Comp yi Pmi Pma

H2S 0.02 34.082 0.68164
CO2 0.08 44.01 3.5208
C1 0.8 16.042 12.8336
C2 0.1 30.069 3.0069
1 20.04294

a) T=150ºF=609.67ºR, P=900psia

Calculamos V en ft3:

1𝑓𝑡 3
𝑉 = 20⁡𝑚3 × = 𝟕𝟎𝟔. 𝟑𝒇𝒕𝟑
1𝑚3
Reemplazando en ecuación 3:

𝒎 = 𝟏𝟗𝟒𝟕. 𝟓𝟗𝟒𝟕⁡𝒍𝒃
Reemplazando en ecuación 2:

𝒏 = 𝟗𝟕. 𝟏𝟕⁡𝒍𝒃⁡𝒎𝒐𝒍
b) n=60%*n inicial, T=130ºF=589.67ºR, V=706.3ft3

𝑛 = 0.6 ∗ 𝑛𝑖𝑛𝑖𝑐𝑖𝑎𝑙
𝒏 = 𝟓𝟖. 𝟑𝟎𝟐
En ecuación despejar P:
𝑛𝑅𝑇
𝑃=
𝑉
𝑷 = 𝟓𝟐𝟐. 𝟐𝟖⁡𝒑𝒔𝒊𝒂

Consideraciones de situación 1 y 2

P1=P2,
𝑛1 𝑅𝑇1
𝑃1 = ⁡⁡⁡(1)
𝑉1
 𝑛2 𝑅𝑇2
𝑃2 = ⁡⁡⁡(2)
𝑉2
Igualando ecuaciones (1) y (2)
𝑚1 𝑅𝑇1 𝑚2 𝑅𝑇2
=
𝑃𝑀1 𝑉1 𝑃𝑀2 𝑉2

Considerando que el PM de la mezcla seguirá siendo el mismo y el volumen del recipiente

cilíndrico que ocupa el gas también será igual, se anula el denominador, y como R se encuentra
en ambos miembros también se anula, quedando la ecuación de la siguiente manera:

𝑚1 𝑇1 = (𝑚1 − 34.73)𝑇2
𝑚1 𝑇1
𝑇2 = ⁡⁡⁡(3)
(𝑚1 − 34.73)
Para reemplazar en ecuación (3), necesitamos obtener el valor de m1
𝑅 × 𝑚1 𝑇1
𝑃1 = ⁡⁡⁡⁡(4)
𝑃𝑀 × 𝑉1
𝜌𝐻𝐶 = 𝜌𝐻2𝑂 × 𝛾𝐻𝐶

𝜌𝐻𝐶 = 40.56⁡ 𝑙𝑏⁄ 3

𝑓𝑡
𝑃𝑀 × 𝑃
𝜌𝐻𝐶 =
𝑅×𝑇
𝑃𝑀 × 𝑃 = 40.56 × 10.73 ∗ 559.67 = 243573.3091
Ordenando ecuación 4:
𝑃 ∗ 𝑃𝑀 ∗ 𝑉1
𝑚1 =
𝑅 × 𝑇1
243573.3091 ∗ (3.1416 × 32 × 9.8425)
𝑚1 =
10.73 × 559.67
𝑚1 = 11287.4477⁡𝑙𝑏

Reemplazando m1 en ecuación 3
11287.4477 × 559.67
𝑇2 =
(11287.4477 − 34.73)
𝑇2 = 561,4⁡º𝑅

Para estimar la Presión de vapor de propano utilizamos la Fig. 23.20 (Anexo 6). Con Temperatura
(ºF), subimos hasta Propano y leemos en lado izquierdo.
 𝑷𝑽 (𝒑𝒔𝒊𝒂) = 𝟑𝟎𝟎⁡𝒑𝒔𝒊𝒂

A partir del Peso Molecular 25 lbm/lbmol, estimamos Pc y Tc

Pc= 691.91 psia

Tc= 15.2373 ºF= 474.9073 ºR

Pr= 1.16

Tr= 1.39

En Fig 23.31 (Anexo 7), leemos conductividad termica (Ka) a P=1atm, posteriormente calculamos
el radio de conductividad termica (k/ka) a partir de Tr y Pr en Fig. 23.32 (Anexo 7). Despejamos
k y obtenemos conductividad termica a presión deseada.
ℎ𝑟 ∗ 𝑠𝑞𝑓𝑡 ∗ º𝐹
𝑘𝐴 = 0.020⁡𝐵𝑡𝑢/( )
𝑓𝑡
𝑘
= 1.2
𝑘𝐴
𝒉𝒓 ∗ 𝒔𝒒𝒇𝒕 ∗ º𝑭
𝒌 = 𝟎. 𝟎𝟐𝟒⁡𝑩𝒕𝒖/( )
𝒇𝒕

Para realizar las lecturas en gráficas, primero debemos calcular PM mezcla, Pr y Tr

Pc Tc Tc Psc Tsc
Comp yi Pmi Pma (Psia) (º F) (º R) (Psia) (º R)
CO2 0.11 44.01 4.8411 1070 87.76 547.43 117.7 60.2173
N2 0.07 28.0135 1.960945 492.5 -232.53 227.14 34.475 15.8998
C1 0.6 16.042 9.6252 667 -116.66 343.01 400.2 205.806
C2 0.11 30.069 3.30759 706.6 89.92 549.59 77.726 60.4549
C3 0.04 44.096 1.76384 615.5 205.92 665.59 24.62 26.6236
C4 0.03 58.122 1.74366 550.9 274.41 734.08 16.527 22.0224
C5 0.02 72.149 1.44298 488.8 305.55 765.22 9.776 15.3044
 8.754 24.203
C7+ 0.02 140 2.8 437.7 750.50 1210.17
689.78 430.53
1 27.485 5029.0 5042.2

En Fig 23.31 (Anexo 8), leemos conductividad termica (Ka) a P=1atm, posteriormente calculamos
el radio de conductividad termica (k/ka) a partir de Tr=1.35 y Pr=0.72 en Fig. 23.32 (Anexo 8).
Despejamos k y obtenemos conductividad termica a presión deseada.
ℎ𝑟 ∗ 𝑠𝑞𝑓𝑡 ∗ º𝐹
𝑘𝐴 = 0.0398⁡𝐵𝑡𝑢/( )
𝑓𝑡
𝑘
= 1.6
𝑘𝐴
𝒉𝒓 ∗ 𝒔𝒒𝒇𝒕 ∗ º𝑭
𝒌 = 𝟎. 𝟎𝟔𝟑𝟕⁡𝑩𝒕𝒖/( )
𝒇𝒕

Primero calcular MWa y z, para poder realizar todos los cálculos

Pci Tci Tci Psc Tsc

Comp yi Pmi Pma
(psia) (º F) (º R) (psia) (ºR)
H2S 0.05 34.082 1.7041 1306.5 212.81 672.48 65.325 33.624
CO2 0.1 44.01 4.401 1070 87.76 547.43 107 54.743
N2 0.08 28.0135 2.24108 492.5 -232.53 227.14 39.4 18.1712
C1 0.46 16.042 7.37932 667 -116.66 343.01 306.82 157.7846
C2 0.12 30.069 3.60828 706.6 89.92 549.59 84.792 65.9508
C3 0.06 44.096 2.64576 615.5 205.92 665.59 36.93 39.9354
iC4 0.05 58.122 2.9061 527.9 274.41 734.08 26.395 36.704
nC4 0.04 58.122 2.32488 550.9 305.55 765.22 22.036 30.6088
nC5 0.02 72.149 1.44298 488.8 385.8 845.47 9.776 16.9094
C7+ 0.02 130 2.6 395.7 658.7 1118.40 7.91455031 22.3679763
1 31.2535 706.38855 476.799176

Calcular Factor de Compresibilidad (Brill & Beggs)

Psc (psia)= 706.38855

Tsc (ºR)= 476.799176
Psr= 4.9547802
Tsr= 1.38353846
 A= 0.34828901
B= 5.35464561
C= 0.0868828
D= 0.95898617
z= 0.75450605

Realizar cálculos para obtener Viscosidad del Gas:

(9.4 + 0.02𝑀𝑊𝑎)𝑇 1.5

𝐾= = 116.1410
209 + 19𝑀𝑊𝑎 + 𝑇
986
𝑋 = 3.5 + + 0.01𝑀𝑊𝑎 = 5.3072
𝑇
𝑌 = 2.1 − 0.2𝑋 = 1.3385
𝑃 × 𝑀𝑊𝑎
𝜌𝑔 = 1.4935𝐸 − 3 × = 0.1084⁡𝑐𝑝
𝑧×𝑇

Primero debemos hallar la relación con el caudal de corrientes 1 y 2

Rel = 4.3/2.8= 1.54

Posteriormente hallamos composición de mezcla

Frac. Molar
Mol gas1+moles
Mol. Mezcla Mezcla
gas 2
Comp yi1 yi2 del Gas 1 (Mol
(yi2+Mol Mezcla
(Rel*yi1) G1+G2/Sum(mol
Gas 1)
G1+G2)

C1 0.805 0.91 1.23625 2.14625 0.85

C2 0.12 0.06 0.184285714 0.244285714 0.10
C3 0.035 0.005 0.05375 0.05875 0.02
iC4 0.01 0 0.015357143 0.015357143 0.01
CO2 0.03 0 0.046071429 0.046071429 0.02
N2 0 0.025 0 0.025 0.01
1 1 2.535714286 1

Para determinar poder calorífico, leer de tablas el PC neto de cada componente

Comp yimez PC neto

C1 0.846 909.4
C2 0.096 1618.7
C3 0.023 2314.9
iC4 0.006 3000.4
CO2 0.018 0
N2 0.010 0
1 7843.4

𝑷𝑪𝒏𝒆𝒕𝒐 = 𝟕𝟖𝟒𝟑. 𝟒⁡𝑩𝒕𝒖/𝒇𝒕𝟑