Segundo coeficiente virial y temperatura de boyle usando

una ecuación de estado aplicada a yacimientos

 

Leonardo Cupil Jiménez

leolipuc@gmail.com

sarai.alejandro@ujat.mx

david.guerrero@ujat.mx

Universidad Juárez Autónoma de Tabasco

Jalpa de Méndez, Tabasco, México

 

RESUMEN

Las ecuaciones de estado (EdE) son modelos matemáticos que permiten calcular las propiedades termodinámicas de las sustancias. En este trabajo se llevó a cabo la estimación del segundo coeficiente virial y la temperatura de Boyle usando una EdE cúbica de tres parámetros desarrollada para aplicaciones en yacimientos de petróleo y gas natural, para determinar su capacidad de predicción a condiciones supercríticas. Los resultados se compararon con valores recomendados en la literatura y se calcularon el error absoluto medio, el sesgo y el error medio cuadrático para determinar la precisión y exactitud de las predicciones, además, estos estadísticos se correlacionaron con el factor acéntrico. Se encontró que la EdE cúbica de tres parámetros predice el segundo coeficiente virial con un error absoluto medio de 0.06, semejante a las EdE cúbicas de dos parámetros; en todos los casos se observa una disminución de la precisión en la predicción cuando las moléculas se alejan del modelo de esfera. Para la temperatura de Boyle se encontró que la EdE cúbica de tres parámetros presenta un comportamiento parabólico con un valor mínimo para un factor acéntrico de 1.80, mientras que las EdE cúbicas de dos parámetros presentan un comportamiento asintótico.

Palabras clave: ecuaciones cúbicas; peng-robinson; soave; esmaeilzadeh-roshanfekr.

Second virial coefficient and Boyle temperature using an

equation of state for reservoirs

 

ABSTRACT

The Equations of State (EoS) are mathematical models to compute thermodynamic properties for substances. In this work, we calculate the second virial coefficient and Boyle’s temperature using a three-parameter cubic EoS developed for oil and natural gas reservoirs application to establish its ability to calculate thermodynamic properties in supercritical regions. The results were compared with recommended values from the literature, and the average absolute deviation, BIAS, and Root-Square Mean Deviation were computed to find the accuracy of the EoS; furthermore, these statistical values were correlated with the acentric factor. For the three-parameter cubic EoS, the second virial coefficient is predicted with an average absolute deviation of 0.06, like the value obtained for two parameters cubic EoS. The three EoS studied indicate a drop of accuracy when the molecule goes away from the spherical model. For Boyle’s temperature, we found that three-parameter cubic EoS shows a parabolic behavior with a minimum value for acentric factor equal to 1.80. In contrast, two-parameter cubic EoS shows an asymptotic behavior.

Keywords: cubic equations; peng-robinson; soave; esmaeilzadeh-roshanfekr.

 

Artículo recibido: 10. Junio. 2021

Aceptado para publicación: 16. Julio. 2021

Correspondencia: david.guerrero@ujat.mx

    Conflictos de Interés: Ninguna que declarar

 

 

 

 

 

 

 

 

 

1. INTRODUCCIÓN

Desde la postulación de la tesis doctoral de Van der Waals (VdW) (Van der Waals, 1873) se ha encontrado una diversidad de aplicaciones para las EdE, principalmente las cúbicas, caracterizadas por ser modelos matemáticos sencillos que permiten predecir el comportamiento PVT de las sustancias y estimar las propiedades termodinámicas de componentes puros y mezclas en fase homogénea y heterogénea, tales como el Equilibrio Líquido – Vapor (ELV). En la literatura existen una infinidad de EdE, que varían en grado de complejidad dependiendo de su aplicación, siendo las más utilizadas en la industria de los hidrocarburos y la simulación de procesos químicos las cúbicas con dos parámetros, tales como la propuesta por Redlich y Kwon (RK) (Redlich & Kwong, 1949), Peng y Robinson (PR) (Peng & Robinson, 1976) o la modificación hecha por Soave a la EdE de RK (SRK) (Soave, 1972). De forma general, la familia de las EdE de VdW pueden ser representadas por

	(1)
	(2)
	(3)

donde , ,  y  son la presión absoluta, la constante universal del gas ideal, la temperatura absoluta y el volumen molar, respectivamente. Los parámetros  y  representan las fuerzas moleculares atractivas y repulsivas; el subíndice  indica la propiedad en el punto crítico, los parámetros  y  son específicos para la EdE y  es una función de la temperatura;  y  dependen de la forma de la EdE, si es de dos parámetros, estarán relacionados con ; en caso de modificaciones para el desarrollo de una EdE de 3 parámetros,  y  pueden estar relacionadas con  u otras propiedades como el factor acéntrico o incluso un nuevo parámetro  (Esmaeilzadeh & Roshanfekr, 2006). La elección del parámetro dependerá de la modificación que se pretenda llevar a cabo en la EdE; algunos autores han intentado mejorar la predicción de las propiedades de la fase líquida, principalmente en el ELV, tal es el caso de la EdE propuesta por Esmaeilzadeh y Roshanfekr (2006) (ER), en la que se propone que el factor de compresibilidad ( ) sea diferente para cada sustancia, mejorando la predicción de las propiedades termodinámicas del líquido saturado, favoreciendo así su aplicación a los fluidos en los yacimientos petroleros. La EdE ER se muestra a continuación:

(4)

sujeta a las siguientes restricciones:

	(5)
	(6)
	(7)

donde el parámetro  se denomina factor de compresibilidad crítico empírico y puede ser un valor constante para cada sustancia o estimarse como una función del factor acéntrico y la temperatura; los parámetros  y  para la Ec. 4 se obtiene a partir de las Ecs. 2 y 3, mientras que el parámetro  se calcula mediante:

(8)

los parámetros ,  y  se calculan mediante el procedimiento descrito en el apartado Estrategias Metodológicas.

La EdE PR se validó para conocer sus capacidades en la predicción de la densidad del líquido y el vapor saturado, la presión de vapor, así como la entalpía y entropía de vaporización, para sustancias puras presentes en los yacimientos petroleros (Esmaeilzadeh & Roshanfekr, 2006). Además, en otros trabajos se han llevado a cabo algunas modificaciones de la ecuación ER, tales como: obtener los parámetros de solubilidad de Hildebrand (Eslamimanesh & Esmaeilzadeh, 2010), modificar la función  para mejorar la predicción del comportamiento  (Bonyadi & Esmaeilzadeh, 2008) y del ELV (Esmaeilzadeh & Samadi, 2008) para las sustancias puras y las mezclas, estimación del coeficiente de Joule-Thompson (Akhouri & Kaur, 2019), predicción de la formación de hidratos de gas (Khan, Kumari, Dixit, Majumder, & Arora, 2020; Talaghat, 2009), aplicaciones en la industria del petróleo y gas como la simulación de yacimientos de gas seco (Ejraei Bakyani, Taghizadeh, Nematollahi Sarvestani, Esmaeilzadeh, & Mowla, 2018) o la precipitación de asfaltenos (Nourbakhsh, Yazdizadeh, & Esmaeilzadeh, 2011), así como la predicción de las densidades en líquidos iónicos (Farzi & Esmaeilzadeh, 2016a) y la tensión superficial (Farzi & Esmaeilzadeh, 2016b). Sin embargo, existen otras propiedades termodinámicas que no han sido verificadas con esta EdE, como el segundo coeficiente virial ( ) y la temperatura de Boyle ( , que se utilizan como punto de validación para la predicción de vapor sobrecalentado y a condiciones supercríticas. En este trabajo se llevó a cabo la predicción del segundo coeficiente virial y la temperatura de Boyle para 16 sustancias, que representan compuestos que se encuentran comúnmente en la industria del petróleo y el gas natural, utilizando la EdE propuesta por Esmaeilzadeh y Roshanfekr, considerando el factor de compresibilidad empírico independiente (ER1) y dependiente (ER2) de la temperatura. Se llevaron a cabo comparaciones con las EdE PR y SRK, tomándolas como referencia debido a su amplio uso en la industria del petróleo y el gas natural.

2. ESTRATEGIAS METODOLÓGICAS

El segundo coeficiente virial ( ) se puede determinar a partir de las EdE, de acuerdo con la definición

(9)

donde  es la densidad molar. Para poder obtener una expresión que permita calcular  a partir de la Ec. 1, es necesario convertir el volumen molar a densidad molar, posteriormente usar la definición , derivar con respecto a la densidad molar a temperatura constante y aplicar el límite en el que , obteniendo:

(10)

La Ec. 10 puede expresarse en su forma reducida, si se sustituyen  y  por las expresiones en las Ecs. 2 y 3 y se define el segundo coeficiente virial reducido como

(11)

De esta forma, la Ec. 11 es función del factor acéntrico y la temperatura reducida. Para la ecuación ER1 y ER2, es necesario llevar a cabo el siguiente procedimiento: primero se encuentra el valor de  tomando la raíz más pequeña obtenida al resolver la función no lineal

(12)

En este trabajo se usó el método de Tartaglia-Cardano (Soto & Mosquera López, 2018) para resolver la Ec. 12, ya que el discriminante permite distinguir si la ecuación no lineal cuenta con una, dos o tres raíces reales y así determinar cuál es la menor. El factor de compresibilidad crítico empírico puede obtenerse como un valor constante para cada sustancia, dependiente del factor acéntrico ( )

(13)

o puede ser considerado dependiente de la temperatura al superar cierto valor para la temperatura reducida ( ), de acuerdo con las siguientes condiciones: si ,  se obtiene de la Ec. 13; si , utilizar la Ec. 14

(14)

para conocer la temperatura reducida límite en el que  es independiente de la temperatura ( ), se usa la Ec. 15

(15)

En caso de no conocer el factor de compresibilidad crítico de la sustancia, se puede usar la correlación propuesta por Nath (Nath, 1982)

(16)

Una vez conocidos los valores de ,  o , se pueden evaluar  y  mediante

	(17)
	(18)

finalmente, el parámetro  se evalúa usando

	(19)
	(20)
	(21)

Para las EdE PR y SRK, los parámetros ,  y  se encuentran en la Tabla 1.

Tabla 1. Parámetros para la ecuación de estado de PR (Peng & Robinson, 1976) y SRK (Soave, 1972).

EdE
PR	0.45724	0.07780
SRK	0.42747	0.08664

Por otra parte, la temperatura de Boyle es aquella en la cual el segundo coeficiente virial se iguala con cero (Mamedov & Somuncu, 2020); a partir de la Ec. 11 se obtiene

(22)

donde  es la temperatura de Boyle reducida. La raíz para la Ec. 22 se obtuvo utilizando el método Gradiente Reducido Generalizado (GRG, por sus siglas en inglés), usando derivadas numéricas centrales y una tolerancia de , tomando como variable independiente el factor acéntrico.

3. RESULTADOS Y DISCUSIÓN

En la Tabla 2 se proporciona la lista de sustancias empleadas para llevar a cabo la predicción del segundo coeficiente virial; una vez obtenidos los valores para esta propiedad, se evaluaron el Error Absoluto Medio (EAM) y el Sesgo (S), de acuerdo con:

	(23)
	(24)

donde  se refiere al número de datos experimentales utilizados en la comparación, para cada sustancia;  es el segundo coeficiente virial reducido, de acuerdo con la Ec. 11, finalmente, los superíndices  y  hacen referencia a los datos calculados y experimentales recomendados (Dymond, Marsh, Wilhoit, & Wong, 2002), respectivamente. Los errores no se normalizaron, ya que al aproximarse a la temperatura de Boyle el segundo coeficiente virial tiende a cero, ocasionando que el error porcentual tienda a infinito. El  proporciona información de la desviación promedio entre los valores calculados con las EdE y los determinados experimentalmente, entre más se aproxima a cero las predicciones son más precisas. Por otro lado, el  es un indicador de la sobrepredicción (si  es positivo) o subpredicción (si  es negativo) de la EdE.

Se observa que las cuatro EdE tienen un  y un  en el mismo orden de magnitud, en promedio; sin embargo, la EdE SRK presenta un  superior para los hidrocarburos con más de cinco carbonos.

Tabla 2. Sustancias y factor acéntrico ( ) utilizados para la estimación del segundo coeficiente virial y la temperatura de Boyle. Resultados para el Error Absoluto Medio (EAM) y el Sesgo (S), para cada sustancia y promedio para cada ecuación de estado.

Sustancia		ER1		ER2		SRK		PR
		EAM	S	EAM	S	EAM	S	EAM	S
Argón (Ar)	0.0000a	0.03	0.03	0.03	0.03	0.04	0.03	0.04	0.00
Monóxido de carbono (CO)	0.0660b	0.01	0.00	0.01	0.00	0.01	0.01	0.02	-0.02
Dióxido de carbono (CO2)	0.2280b	0.02	0.00	0.02	0.00	0.02	0.02	0.02	-0.01
Metano (C1)	0.0120b	0.02	0.00	0.02	0.01	0.02	0.01	0.03	-0.02
Etano (C2)	0.1000b	0.02	0.00	0.02	0.00	0.02	0.02	0.03