NTP 430: Gases licuados: evaporación de fugas y derrames

Liquefied gases: Evaporation of releases and spillages
Gaz liquéfiés: Evaporation de fuites et écoulements

Evaporación de gases licuados a presión

El estudio del fenómeno de la evaporación de un líquido o de un gas licuado, por fuga y derrame accidental, se puede realizar mediante modelos matemáticos para predecir la intensidad de evaporación en unas determinadas situaciones que conducen a la formación de una nube cuyos efectos dependerán de la magnitud de la evaporación. Dicho estudio del fenómeno requiere conocer el tipo de liquido derramado, las condiciones y tipo de emisión, el substrato y el área de confinamiento en que se derrama.

Esta situación se da cuando se produce una fuga y derrame por fallo del recipiente contenedor, accesorio, equipo o conducción de un gas licuado a presión. La temperatura inicial del líquido es la ambiental si el depósito está en equilibrio térmico con el medio ambiente. El gas licuado que escapa y se derrama normalmente está a una temperatura superior a su temperatura de ebullición a presión atmosférica, de ahí su nombre de líquido sobrecalentado.

Al producirse la pérdida de contención, el gas licuado sufre un descenso súbito de su presión de almacenamiento (la presión de vapor correspondiente a la temperatura de almacenamiento) hasta la presión atmosférica. Esto da lugar a una evaporación súbita, también llamada evaporación flash, que al mismo tiempo puede arrastrar una cantidad considerable de líquido en forma de aerosol (pequeñas gotas). Una parte de esas gotas puede volver a caer al suelo en forma de lluvia (rain out) debido al enfriamiento y condensación y otra parte se evaporará a causa de la absorción de calor procedente del aire ambiente que se mezcla con la emisión de fluido.

La parte de gas licuado restante forma un charco de líquido en el suelo que está a una temperatura igual a la de ebullición la cual es muy baja con respecto a la temperatura del suelo. Si el derrame es de poca magnitud, la evaporación del charco es muy rápida y todo el gas licuado se incorpora a la atmósfera en muy poco tiempo. En derrames importantes, tras la rápida evaporación inicial, el suelo (substrato) sufre un enfriamiento suficiente para que la subsiguiente evaporación del charco transcurra de forma menos rápida.

Cálculo de la evaporación rápida inicial

Esta evaporación ocurre por una expansión rápida que se puede considerar adiabática, es decir sin transferencia de calor con su medio circundante. El calor necesario para la evaporación (calor latente de evaporación) es extraído del exceso de energía del propio líquido sobrecalentado que sufre un descenso de temperatura, desde la existente en el momento del derrame T₀, hasta la de ebullición a presión atmosférica T_eb.

La masa inicial evaporada se puede calcular del balance de energía, en que se iguala la cantidad de calor necesario para la evaporación con la cantidad de calor proporcionado por la masa de gas licuado, obteniendo la expresión

La expresión anterior es válida para la evaporación de un solo componente. En mezclas multicomponentes, el cálculo manual es complicado y se pueden utilizar programas informáticos para tales cálculos.

El valor de la fracción vaporizada f_v = m_v/m_o suele ser muy inferior a los valores observados experimentalmente, lo cual se debe principalmente al arrastre de líquido en forma de gotas pequeñas (aerosol). La formación de este aerosol aumenta la masa total de gas licuado que se incorpora a la nube y afecta a las características de ésta, aumentando su densidad aparente y causando un enfriamiento adicional (por evaporación de las gotas de aerosol que roban calor del aire entremezclado, y que puede conducir a la condensación de la humedad del aire atmosférico). La predicción fiable del grado de formación de aerosol es difícil, requiriendo modelos complicados del proceso. Kletz propone una estimación simplificada, sugiriendo que la fracción procedente del aerosol iguala a la evaporación inicial, con lo que un valor aproximado de la fracción líquida que pasa a formar la nube inicial se considera igual a 2f_v. En la bibliografía especializada se citan trabajos experimentales en que no queda gas licuado remanente tras la evaporación inicial, a pesar de que la fracción teórica vaporizada era solamente el 20% (y en algunos casos el 10%). En ausencia de cálculos más sofisticados, el criterio más conservador sería admitir que toda la masa vertida se incorpora a la nube, excepto para los casos en que el cálculo de la fracción teórica vaporizada sea inferior al 10 -20%, en cuyo caso se podría aplicar la corrección sugerida por Kletz.

La figura 1 y la figura 2 muestran el valor de f_v en porcentaje para un grupo importante de gases licuados a distintas temperaturas de almacenamiento. La tabla 1, que se muestra más adelante, presenta valores de las temperaturas de ebullición T_eb, calores específicos c_p y calores latentes de vaporización h_v de los gases más comunes con los números de referencia para las dos figuras.

Tabla 1: Propiedades de los gases licuados a presión más comunes

Cálculo de la evaporación del gas licuado derramado sobre tierra

Esta situación se da en la evaporación de derrames de gases licuados por enfriamiento (refrigerados) que están a su temperatura de ebullición o próximos a ella y en la evaporación de derrames procedentes de gases licuados a presión después de haber ocurrido la evaporación rápida inicial y quedar un charco de gas licuado a su temperatura de ebullición.

Un aspecto importante para el cálculo de la evaporación es el tipo de substrato sobre el que se extiende el charco de gas licuado.

Evaporación de gases licuados en suelo no permeable

Cuando un gas licuado a su relativamente baja temperatura de ebullición entra en contacto con un terreno seco a mayor temperatura y que por ser impermeable no penetra en él de una manera significativa, se establece una transferencia de calor del medio a mayor temperatura (suelo) al de menor (masa de gas licuado).

La solución matemática a este fenómeno físico da el flujo de calor q_s (W/m²) transferido por conducción desde el terreno a la masa de gas licuado del charco

El flujo de evaporación m_e del gas licuado extendido formando un charco sobre un suelo no permeable y en función del tiempo será:

siendo e un coeficiente independiente del tiempo en unidades [kg/m² s^1/2] la expresión del flujo de evaporación en función del tiempo queda así:

Para la aplicación de la fórmula se necesitan los valores de las constantes de los gases licuados implicados y de los materiales del substrato que se indican en las tabla 1 y la tabla 2.

Tabla 2: Constantes del material de algunos suelos

Si el charco de gas licuado formado requiere un tiempo largo para su completa evaporación, el terreno ya no podrá aportar energía calorífica y a partir de un cierto momento debe tenerse en cuenta la influencia de la radiación solar y la evaporación ocasionada por la convección del aire circundante.

En cuanto a la radiación solar varía notablemente según la hora del día, la época del año, las condiciones ambientales y la latitud geográfica del lugar. Un valor promedio de la radiación solar es q_rs = 400 W/m² y en bibliografía especializada del tema se citan valores típicos de radiación solar máxima (a mediodía en días claros), oscilando entre 425 W/m² (enero) y 1170 W/m² (julio) para una latitud de 45° norte. Para una latitud de 35° norte los valores correspondientes serían 650 W/m² y 1140 W/m², respectivamente.

En un derrame que ocurra sobre una superficie aislada térmicamente, el efecto de la radiación solar o de la convección del aire circundante puede ser dominante al cabo de un período de tiempo corto.

Se recomienda comparar el valor calculado anteriormente con el que daría el flujo de evaporación de un líquido con temperatura de ebullición superiora la temperatura ambiente mediante la fórmula de Sutton y elegir el valor mayor de los dos. En general esta fórmula se emplea para líquidos cuya presión parcial de vapor es igual o menor que 2x10⁴ Pa. Esta fórmula tiene la siguiente expresión para una condición de estabilidad atmosférica neutra (ver NTP 329-1993).

m_e = 2 x 10^-3 · u₁₀^0,78 . r^-0,11 . M_m P_t / (RT). In [1+ (P_vs - P_{v air}) / (P_t -P_vs)]

m_e = 2 x 10^-3 · u₁₀^0,78 . L^-
0,11 . M_m P_t / (RT). In [1+(P_vs - P_{v air}) / (P_t -P_vs)]

Evaporación de gases licuados refrigerados en suelo permeable y seco

Las expresiones utilizadas para esta situación se basan en los ensayos experimentales realizados con un gas licuado refrigerado, en particular el gas natural licuado , evaporados en suelo arenoso y seco.

El flujo de evaporación producido por conducción térmica desde el suelo se ajusta a la ecuación

Si se utiliza la fórmula de evaporación de un gas licuado en suelo no permeable

me = e t^-1/2 en donde e = l_s. (T_s -T_eb) / [h_v (p a _s)^1/2] y se sustituyen los parámetros correspondientes al gas natural y a un suelo arenoso seco a 20°C se obtiene

Estableciendo la razón entre el valor experimental y el correspondiente a suelo no permeable se tiene

Esto indica que el flujo de evaporación calculado teóricamente es unas ocho veces menor que el determinado experimentalmente. La explicación a esta diferencia puede hacerse por la influencia de la penetración del líquido en el suelo. Según los datos experimentales observados, la fórmula empírica también se puede aplicar al propano y butano derramado sobre suelo arenoso y seco.

De todo lo razonado, para el flujo de evaporación de un gas licuado, derramado sobre un suelo permeable y seco y en función del tiempo, se puede aplicar la expresión siguiente

Igualmente se debe remarcar que para tiempos largos, se recomienda hacer el cálculo de la evaporación ocasionada por convección para líquidos (fórmula de Sutton) y adoptar el valor mayor obtenido.

Aplicación práctica

Calcular el flujo de evaporación en función del tiempo, de un derrame instantáneo ocasionado por la rotura de un recipiente soportado por pilares y que contiene 6000 kg de cloro licuado a 20°C sobre un cubeto circular de 5 m de diámetro cuyo suelo es de cemento.

Solución: La evaporación se debe separar entre la evaporación rápida inicial (flash) y la evaporación del charco que se extiende en las inmediaciones del recipiente contenedor y en este caso en la superficie del cubeto de retención.

Cálculo de la evaporación rápida inicial

Esta masa es el dato inicial de la siguiente fase a considerar, esto es la dispersión de una emisión instantánea, que será tratada en otra NTP.

La fracción evaporada inicialmente es f_v= m_v/ m₀ = 972,8/ 6000 = 0,162 = 16,2%.

Considerando la recomendación de T. Kletz se debería adoptar como valor más real de la masa de la nube formada, el doble de la masa teórica calculada, es decir 2 x 972,8 = 1 945,6 kg

La cantidad no evaporada inicialmente que llamaremos masa final m_f, será por tanto:

m_f = 6000 - 972,8 = 5 027,2 kg en cálculo teórico y m_f = 6 000 - 1945,6 = 4054, 4 kg en el más real.

Cálculo de la evaporación del charco formado en función del tiempo

La cantidad no evaporada inicialmente cae al suelo y se extiende en la superficie del cubeto.

El flujo de evaporación por unidad de superficie y en función del tiempo será:

m_e = 1, 1(293-239) / [2,9x10⁵(p. 10 x 10^-7)^1/2 [ t^-1/2 = 0,1155615 t^-1/2 [kg/(m². s)]

Considerando la superficie total de evaporación del cubeto A = p r² = p 2,5² = 19,63 m² se tiene

Esta ecuación nos daría los valores de la velocidad o tasa de evaporación para distintos tiempos.

Se puede determinar el tiempo a partir del cual empieza a ser significativo el calor recibido por radiación del sol, que varía entre 425 W / m² en enero y 1170 W / m² en julio, para una latitud de 45° norte y a mediodía en días claros. Para ello se toma la expresión que da el flujo calorífico procedente del terreno en función del tiempo

Se puede considerar significativo el efecto de la radiación solar cuando sea como mínimo el 10 % del flujo calorífico recibido del suelo. En enero se tendría 425 = 10%.4 250 W / m² para lo cual tendrían que transcurrir

En julio a la hora de mayor radiación se tendría el 10% cuando el flujo calorífico del suelo fuera 11 700 W / m² lo cual ocurriría al cabo de

Ahora se puede comparar con el resultado que proporcionaría la ecuación de Sutton aplicada en la evaporación de líquidos por convección, suponiendo que no fuera un gas licuado, en la superficie circular del cubeto, con una velocidad del viento de 2 m/s y considerando el límite máximo de presión parcial de vapor 2x10⁴ Pa y P_vair = 0 para la aplicación de la fórmula

m_e = 2 x 10^-3. 2^0.78 . 2,5^-0,11. 71.10⁵ / (8310.239) In [1 + (2 x 10⁴-0) / ( 10⁵-2 x 10⁴)] = 0,0024768 kg / (m². s)

Comparando con la tasa de evaporación calculada anteriormente se observa que el efecto de la transferencia de calor por convección en un líquido con una presión parcial de vapor de 2 x 10⁴ Pa tardaría unos 2400 s (40 min) en dar la tasa de evaporación del gas licuado considerado (cloro).

Los cálculos realizados teniendo en cuenta sólo la absorción de calor del suelo dan un tiempo de permanencia del gas licuado derramado superior al real ya que según se ha visto, el efecto de la radiación solar en verano empieza a los pocos segundos y en invierno al cabo de un minuto aproximadamente y se ha despreciado la evaporación por calor de convección, por lo que en cálculos más precisos se debería considerar el efecto combinado de las tres fuentes de calor. Por ejemplo, el flujo de evaporación m_rs debido a la radiación solar en el mes de julio es igual a

En los 19,63 m² del cubeto sería m_rs.19,63 = 0,0791952 kg/s »0,08 kg/s que podría ser añadido y dar un flujo total de evaporación en función del tiempo.

Si se desea conocer la masa total M_e (t) evaporada al cabo de un tiempo t, se calcularía

Teniendo en cuenta la influencia de la radiación solar se incluiría el flujo constante m_rs = 0,08 kg/s multiplicado por el tiempo transcurrido, con lo cual se tendría:

Se puede completar la tabla 3 con los valores anteriores y la masa total evaporada en función del tiempo.

Tabla 3: Flujos de evaporación y masa total evaporada en función del tiempo en el cubeto de retención del derrame de cloro

Bibliografía

(1) SANTAMARÍA, J.M. y BRAÑA, P.A.
Análisis y reducción de riesgos en la industria química.
Madrid, Editorial MAPFRE, 1994.

(2) TNO.
Methods for the calculation of physical effects resulting from releases of hazardous materials (liquids and gases). 2nd ed.
Voorburg, The Netherlands, TNO, 1992.


Fig. 1: Fracción teórica (calculada) vaporizada inicialmente en la expansión adiabática de un gas licuado a presión (número de referencia en tabla 1)	Fig. 2: Fracción teórica (calculada) vaporizada inicialmente en la expansión adiabática de un gas licuado a presión (número de referencia en tabla 1)