NTP 385: Fugas en recipientes: emisión en fase gaseosa

NTP 385: Fugas en recipientes: emisión en fase gaseosa Fuites dans recipients: emmission en phase gazeuse Leakages in vessels: emission in gas phase Redactor: Emilio Turmo Sierra Ingeniero Industrial CENTRO NACIONAL DE CONDICIONES DE TRABAJO Este documento forma parte de una serie de esta misma colección dedicada al análisis de consecuencias originadas por la emisión o escape de sustancias peligrosas a la atmósfera, que pueden dar lugar a accidentes de diversa gravedad, en función de la cantidad emitida, de su peligrosidad y de la proximidad a personas o a instalaciones.
La emisión en fase gas se puede dar en muchos de los procesos existentes en una planta industrial y un aspecto previo a considerar serán las diferentes situaciones que pueden crear este peligro, antes de evaluar el propio fenómeno de la emisión.
Introducción La emisión de sustancias peligrosas es la primera fase a estudiar, una vez se tienen relacionadas las situaciones o escenarios de peligro mediante una cuidadosa auditoría de los procesos e instalaciones existentes.
En esta Nota Técnica se consideran las emisiones en fase gaseosa que pueden ocurrir en un recipiente de almacenamiento o proceso.
La relación entre presión, temperatura y volumen de un gas real es compleja, aunque a excepción de presiones elevadas, el comportamiento de los gases se aproxima al de un gas ideal, para el cual el volumen de una masa dada es inversamente proporcional a la presión y directamente proporcional a la temperatura absoluta.
En un proceso isotérmico (T = constante) se cumple: pV = nRT = constante p = presión absoluta del gas (Pa).
V = volumen ocupado por el gas (m3).
n = nº de kilomoles contenidos en el volumen V de gas (kmol).
R =constante universal de los gases (8314J /(kmol K)).
En un proceso adiabático (sin absorción o cesión de calor) se cumple de forma aproximada: pVk = constante Si la expansión se aproxima a un proceso adiabático reversible (isoentrópico), el coeficiente adiabático k se transforma en la razón de los calores específicos γ.
γ= Cp/Cv Cp = calor específico a presión constante (J/(kg K)).
Cv = calor específico a volumen constante (J/(kg K)).
Las capacidades caloríficas de los gases ideales dependen solamente de la temperatura y los gases reales se aproximan tanto más a un gas ideal cuanto menor sea la presión.
En expansiones no isoentrópicas el coeficiente adiabático k es menor que γ.
Las condiciones isotérmicas solo se pueden mantener si existe una buena transferencia de calor con el entorno y esta situación ocurre con caudales reducidos y tiempos prolongados de emisión en equipos pequeños.
Las condiciones más cercanas a un proceso adiabático se dan en situaciones contrarias, es decir grandes instalaciones con caudales elevados y tiempos reducidos.
Ejemplos típicos de emisión por fugas en fase gaseosa se tienen normalmente en las situaciones siguientes: ● Conexiones de tuberías y equipos.
● Roturas de tuberías y equipos.
● Orificios en equipos (conducciones, recipientes) que contienen gas a presión.
● Descarga en válvulas de alivio de seguridad (sólo para fase vapor).
Las NTP son guías de buenas prácticas.
Sus indicaciones no son obligatorias salvo que estén recogidas en una disposición normativa vigente.
A efectos de valorar la pertinencia de las recomendaciones contenidas en una NTP concreta es conveniente tener en cuenta su fecha de edición.
Año: 199 En lo referente a los posibles escenarios que dan origen a la emisión de sustancias peligrosas y aspectos a considerar en un escape, se remite al lector a la NTP362-1994 "Fugas en recipientes y conducciones.
Emisión en fase líquida".
El caudal de emisión de cualquier fluido y en nuestro caso de un gas o vapor por un orificio, se determina mediante las relaciones fundamentales de la mecánica de fluidos, que se pueden encontrar en los textos clásicos de esta materia.
Las expansiones que tienen lugar en las situaciones de emisión, son normalmente de uno de los dos tipos: expansión con estrangulamiento y expansión libre.
La expansión con estrangulamiento ocurre cuando el gas sale a través de una abertura o grieta muy estrecha.
Se tiene una gran resistencia de fricción y una energía cinética baja.
Cuando la abertura es mayor se alcanza una situación de expansión libre.
La resistencia por fricción es baja y la energía cinética elevada.
Las situaciones de mayor interés en nuestro estudio por el riesgo que implican, corresponden a una expansión libre, la cual se considera aproximadamente como adiabática.
En el estudio del flujo de un gas, caracterizado por una masa molecular M y una razón de calores específicos γ, por un orificio de características determinadas (diámetro y coeficiente de descarga Cd), en un recipiente a presión interna absoluta p y temperatura termodinámica T (absoluta), a su entorno ambiental exterior a presión atmosférica pa, se demuestra que mientras la razón de las presiones p/pa se mantiene superior a la razón crítica rcrít = [(γ+1)/2]γ/(γ-1) el caudal de emisión sólo depende de la presión p y es independiente de la presión pa.
En cambio, cuando la razón de las presiones p/pa es inferior a la rcrít, el caudal de emisión depende de p/pa.
Para seguir el fenómeno completo de la emisión, se puede considerar desde el instante inicial en que se vacía un recipiente a presión p0, hasta que finaliza cuando se iguala su presión interior con la atmosférica.
Al principio de ésta experiencia, que es lo que ocurre normalmente en todos los escapes de gas o vapor, se tiene un régimen de flujo sónico, hasta que se llega a la razón crítica de presiones tal que: p/pa = [(k+1)/2]k/(k-1) Si el proceso se considera adiabático, no isotérmico o intermedio entre adiabático e isotérmico se utiliza el valor k.
Si el proceso se considera isoentrópico (adiabático y reversible), el valor k se sustituye por γ.
Continuando con la experiencia del vaciado completo del recipiente, por debajo de la relación crítica se pasa a tener un flujo subsónico.
La razón de tal denominación viene de que en el flujo sónico la velocidad de escape del gas corresponde aproximadamente a la velocidad del sonido en el gas y en el flujo subsónico escapa a velocidades inferiores a la del sonido.
Existen escapes por toberas especiales en que la velocidad es supersónica, pero cuyo estudio trasciende el objeto de este documento.
Caudal de emisión de gas en un recipiente Caudal másico inicial (proceso adiabático) Para la deducción de la expresión del caudal másico inicial se considera un gas ideal y un proceso de expansión adiabático y reversible (isoentrópico), con lo cual se dispone de una fórmula aproximada para un caso real.
Una emisión real se aproxima a un proceso adiabático y a partir de ésta hipótesis se llega a una expresión de la siguiente forma: Qm (0) = Cd·A·Ψ[p0 ρ0γ [2/(γ+1)](γ+1)/(γ)-1]1/2 La equivalente, teniendo en cuenta la relación ρ0 = p0M/ (RT0) en las condiciones iniciales: Qm (0) = Ψ· Cd·A·p0 [ M/RT0) γ 2/(γ+1)](γ+1)/(γ)-1]1/2 Qm (0) = caudal másico inicial de gas (kg/s) Cd = coeficiente de descarga (sin dimensiones) A = área de la sección del orificio de salida (m2) Ψ= coeficiente del flujo de emisión o factor de flujo po = presión inicial del gas en el recipiente (Pa) ρ0 = densidad inicial del gas en el recipiente (kg/m3) γ = razón entre los calores específicos a presión y volumen constantes (Cp/Cv) M = Masa molecular (kg/kmol) R = constante universal de los gases (8314 J/ (kmol · K)) Para facilitar las operaciones se puede calcular previamente el término: ß = [ γ [ 2/(γ + 1)] (γ+1)/ (γ-1)] 1/2 También se emplea para el caudal másico en función del tiempo Qm(t), operación que se verá en la aplicación práctica.
El coeficiente de descarga Cd es un valor empírico inferior a 1 para tener en cuenta las pérdidas por la viscosidad y la turbulencia.
Este coeficiente es la relación entre el caudal real que pasa a través del orificio y el caudal ideal.
Un valor común es 0,85 para casos de rotura y es la estimación óptima para emisiones previstas tecnológicamente (válvulas de alivio, venteos).
Para un orificio de bordes afilados se suele tomar 0,61.
Para una boquilla bien redondeada en su unión con un recipiente, es aproximadamente la unidad.
El factor de flujo viene dado por la expresión: Ψ= { 2/(γ-1) [(γ+1/2] (γ+1) (γ-1) (pa/p) 2/γ [ 1-(pa/p) (γ-1)/γ]}1/2 y es función de la razón crítica de presiones [(γ+1)/2] γ/(γ-1).
Se aplica durante la emisión con flujo subsónico (subcrítico) en que p/pa < [(γ+1)/2] γ/(γ-1).
Si al comienzo con las condiciones iniciales se tiene po/ Pa ≥ [(γ+1)/2γ/(γ-1) dicho factor Ψo = 1 y el flujo se llama sónico o crítico.
La presión p es función del tiempo, y se calcula según se indica más adelante.
Dibujando un gráfico de y en función de p/pa y teniendo en cuenta que γ varía entre 1,1 y 1,4 se observa que Ψdepende muy poco del valor de γ y que Ψvale aproximadamente 1 cuando p ≥ 1,5 Pa (ver Fig.
1).
Fig.
1: Factor de flujo Ψen función de p/pa Caudal másico en función del tiempo (proceso adiabático) Haciendo las mismas suposiciones que en el caudal másico inicial se llega a la expresión siguiente: Qm (t) = Cd · A · ß (p0 ρo)1/2 [1 + 1/2 (γ-1) tr] -(γ+1)/γ-1) o en función del caudal másico inicial: Qm (t) = Cd · A · ß (poρo)1/2 [1 + 1/2 (γ 1) tr] -(γ+1) / (γ -1) En donde: Qm (t) = dm/dt = caudal másico (kg/s) en el tiempo t (s) ß = { γ [2/(γ+1)] γ+1)/(γ-1) } 1/2 tr se define como tiempo reducido (sin dimensiones) cuyo valor es: tr = Cd · A (po ρo) 1/2 mo -1 · ß t mo= masa inicial de gas en el interior del recipiente (kg) t = tiempo transcurrido desde el instante inicial (s) Los demás parámetros de estas expresiones ya se han indicado anteriormente.
Para gases ideales tr se puede expresar de esta forma: tr = Cd · A (RTo/M) 1/2 Vr -1 · ß t Vr = volumen del recipiente (m3) To = temperatura inicial absoluta (K) Si se conoce la presión p(t) en función del tiempo por las expresiones del apartado siguiente, se puede calcular el caudal másico por la ecuación: Qm (t) = Qm (0) · [p(t)/po] (1+γ)/2γ en la que Qm(0) es el caudal másico inicial, calculado por una de las dos ecuaciones propuestas.
Presión en función del tiempo (proceso adiabático) La presión del gas en el interior del recipiente en función del tiempo va disminuyendo a medida que se vacía su contenido y es interesante conocer su valor, el cual está expresado por la ecuación siguiente: p(t) = po {1 + (γ-1)/2 · ΨCdA/mo · (poρo)1/2 · [γ [2/(γ+1)] (γ+1)/(γ-1)]1/2 · t} -2γ/(γ-1) = po [ 1 + 1/2 (γ-1) tr] -2γ/(γ-1) en la cual se incluye tr definido anteriormente.
Otra expresión equivalente es: p (t) = po { 1 + (γ-1) / 2 · ΨCd A / Vr [γ [2/(γ+1)] (γ+1) / γ-1)]1/2 · [8314To / M] 1/2 t }-2γ/(γ-1) Para valores de γ 1 la fórmula expresada en función de tr pasa a ser: p (t) = po exp (-γ tr) Tiempo en que se alcanza la razón crítica de las presiones rcrít Se obtiene haciendo p(t)/pa = rcrít y despejando t.
Tener en cuenta que Ψ= 1.
Temperatura en función del tiempo Al haberse considerado el proceso adiabático y reversible, la temperatura interna del gas va disminuyendo progresivamente conforme se vacía el recipiente.
La expresión que nos determina la variación de esa temperatura en función del tiempo se puede calcular a partir de la presión interna existente en cada tiempo considerado.
Se desarrolla a partir de las dos fórmulas de un gas ideal considerando la expansión adiabática pV n RT y pVk = Cte, obteniéndose: T = To [p (t)/po] (k-1)/k y si es isoentrópica: T = To [p (t)/po] (k-1)/γ Aplicación práctica Determinar los valores correspondientes del caudal másico, presión y temperatura en función del tiempo, en una emisión de gas etileno a través de un orificio con un área de sección 0,003 m2 en un tanque de servicio para alimentación de un proceso.
La capacidad del tanque es 50 m3 de gas etileno a una temperatura ambiente de 17ºC y a una presión absoluta de 30 bar.
Se supone que la transferencia de calor con el entorno es despreciable, de forma que el proceso durante la emisión se puede considerar isoentrópico (adiabático y reversible).
Datos Coeficiente de descarga Cd = 0,61 Temperatura ebullición t= -104ºC Temperatura crítica tc = 9,7 ºC Presión crítica pc= 50,5 atm (5,12 MPa) Peso o masa molecular del etileno (eteno) C2H4 28,05 Razón Cp/Cv γ = 1,18 Solución El etileno es un producto obtenido del petróleo que normalmente se almacena en forma de gas licuado refrigerado, ya que su temperatura crítica 9,7 ºC es inferior a la temperatura ambiente y a esta temperatura sería imposible mantenerlo en estado líquido.
En este ejercicio se considera que el escape ocurre en un recipiente previo a un proceso y en el cual sólo existe fase gaseosa.
Las fórmulas empleadas se han deducido haciendo la suposición de gas ideal.
Para la aplicación de las fórmulas correspondientes se realizan unos cálculos previos que se indican paso a paso.
El caudal másico inicial (proceso adiabático) se calcula con la expresión: Qm (0) = Cd · A · Ψo { p0 ρ0 γ [2/ (γ+1)] γ+1) / (γ -1) } 1/2 = Cd · A · Ψ0 (p0 ρ0) 1/2 { γ [2/(γ+1)] (γ+1) / (γ-1) } 1/2 La densidad inicial ρo del gas se obtiene a partir de la conocida fórmula de los gases ideales pV = nRT junto con ρ= m/V y n = m/M. ρ0 = M p0 / (R/T0) = (28,05 · 30 · 105) / (8314 · 290) = 34,90 kg /m3 en donde el subíndice cero indica las condiciones iniciales.
M = 28,05 kg / kmol P0 = 30 · 105 Pa T0 = 290 K La razón tomada entre los calores específicos γ = Cp /Cv = 1,18 es un valor constante suponiendo un comportamiento del etileno como gas ideal.
Para su comportamiento como gas real se debería tomar un valor variable con la temperatura.
Para facilidad de cálculo se podría elegir un valor promedio en el campo de las variables del proceso.
El valor { γ [ 2/ (γ +1)] (γ+1) / (γ -1) } 1/2 que llamaremos ß se calcula y da: ß = { 1,18 [ 2/ (1,18 +1)] (1.
18+1) / (1,18-1)} 1/2 = 0,64446185 La relación crítica de presiones: rcrit = [(γ +1) / 2] γ / γ-1) = [(1,18 +1) / 2] 1,18 / (1,18-1) = 1,759 nos indica que mientras p/pa > 1,759 se tendrá flujo sónico y para p/pa < 1,759 se tendrá flujo subsónico.
Para el campo, de flujo sónico Ψ= 1 y sustituyendo valores en la expresión correspondiente se tiene el caudal másico inicial: Qm (0) = 0,61 · 0,003 · 1 [ 30 · 105 · 34,90] 1/2 · 0,6556185 = 12,07 kg/s.
La cantidad o masa inicial m0 de etileno es igual a: m0 = Vr · ρ0 = 50 · 34,9 = 1745 kg.
Vr = volumen del recipiente El caudal másico en el tiempo t se puede calcular según la 2ª de las expresiones propuestas, una vez ya se ha calculado el caudal másico inicial Qm (0).
Qm (t) = Qm (0) · [ 1 + (1/2) (γ-1) tr] (γ+1) / (γ-1) tr = Cd A (po ρo 1/2 mo-1 · ßt = 0,61 · 0,003 · (30 · 105 · 34,9) 1/2 · 0,6446185 · t / 1745 = 6,917 · 10-3 t Qm (t) = 12,07 [1 + (1/2) (1,18 -1) · 6,917 · 10-3t] -(1,18 +1) / (1,18-1) = 12,7 [1+6,225 · 10-4 t] -12,11 con la cual se tendrá el caudal másico en función del tiempo y se podrán obtener los valores para diferentes tiempos.
El tiempo en que se alcanza la razón crítica de las presiones se calcula de la siguiente forma: En ese instante se debe calcular el valor: p(t) = p0 { 1+(γ-1) / 2 · ΨCd A / Vr [ γ [2/(γ +1) ] γ+1) / (γ-1) ] ½ [8314T0 / M] ½ · t} -2γ / (γ-1) = 30 · 105 { 1 + (1,18-1) / 2 · º · 0,61 · 0,003/50 [1,18 [2/(1,18+1)] (1,18 + 1) / (1,181)] 1/2 [8314 · 290/28,05] ½ · t} -2.
1, 18/ (1,18-1) = 30 · 105 (1+6,225 · 10-4 · t) -13,11 p(t) / pa = rcrít = [(γ+1) / 2] γ/(γ-1) = 1,759 = 30 · 105 / (1,01325 · 105 · (1+6,225 · 104 · t) -13,11 y despejando t t = [ (1,759 / 29,607) 1/-13,11 -1] / (6,225 · 10-4) = 385,98 ≈386 s La presión p en el recipiente en función del tiempo t se calcula con la expresión: p (t) = po { 1+ (½) (γ -1) ΨCd A (po ρo) ½ mo -1 · [γ [2/(γ+1)] (γ+1) / γ-1{-2γ / (γ-1) = po [ 1 + (½) (γ-1) Ψtr ] -2γ/γ-1) = 30 · 105 [1+ (½) (1,18 -1) · 6,917 · 10-3 t ] -2.
1,18 / (1,18 -1) pudiéndose calcular los valores indicados en la tabla 1, al ir dando valores al tiempo t y considerando la fase sónica con Ψ= 1.
Para la fase subsónica como p(t) es función del tiempo y también del factor de flujo Ψy a su vez Ψes función de p, se puede recurrir al método de cálculo expuesto posteriormente.
Tabla 1: Variación de los parámetros de la emisión en fase gaseosa Para determinar los valores de la temperatura del gas en el interior del recipiente en función de su presión anteriormente calculada, se aplica la expresión: T (p) = To [p(t)/po ] (γ-1) / γ T (p) = To { [1 + (1/2) (1,18 1) · 6,917 · 10-3 t ] (γ-1) / γ la cual también a su vez se puede calcular en función del tiempo, sustituyendo los valores p (t) por su ecuación equivalente en función de t: T (t) = To {1 + (1/2) (1,18 1) · 6,917 · 10-3 t} -2.
1, 18 / (1,18 1) } (1,18 1) / 1,18 = 282,91 K para t = 20 segundos Dando valores al tiempo se tendrán las temperaturas correspondientes, conformando la tabla 1 con los diversos parámetros del fenómeno del escape de un gas contenido en un recipiente.
Para realizar los cálculos en la fase subsónica, se divide ésta en un cierto número de intervalos, yendo cada uno de ellos desde la presión p(t)n a la p(t)n+1 para valores sucesivos de n y en cada intervalo se considera un valor constante para el coeficiente o factor del flujo Ψel cual se supone igual al valor promedio de Ψn y Ψn+1.
Para limitar la inexactitud de la operación de promediado de Ψa un porcentaje bajo, se recomienda elegir los valores de pa/ p (t)n, n+1 que se indican en el párrafo siguiente para determinar el tamaño de esos intervalos.
Como la relación pa/p(t) en el punto crítico de cambio de fase sónica a subsónica vale aproximadamente entre 0,58 y 0,48 para valores corrientes de y entre 1,1 y 1,7 y en nuestro caso con γ = 1,18 y 1/ rcrit = [2 / (γ + 1)] γ / (γ 1) = 1: 1759 = 0,568 se adoptan normalmente los intervalos: 1/rcrit a 0,6 / 0,6 a 0,7 / 0,8 a 0,9 / 0,9 a 0,95 y sucesivas adiciones de la mitad del intervalo que queda 1 [pa / p (t)].
Así el siguiente intervalo iría de 0,95 a 0,95 + (1/2) (1 0,95) = 0,975 y así sucesivamente.
Con las expresiones incluidas en este documento y aunque sea laborioso se pueden calcular todos los valores presentados en la tabla 1.
Preparando un programa con ordenador, la confección de la tabla es más rápida y cómoda, permitiendo una presentación similar a la adjunta.
Como comprobación final se puede verificar que la cantidad inicial de gas menos laque queda en el recipiente debe ser igual a la cantidad fugada.
Cantidad inicial: 1745 kg Cantidad que queda: m = Vr · ρ172,97K ρ172,97K = M · p172,97K / (R.172,97) = 28,05 · 1,01325 · 105 / (8314 · 172,97) = 1,976 kg/cm3 m = 50.
1,976 = 98,82 kg 1745 98,82 = 1646,18 kg valor muy aproximado al reflejado al final de la última columna de la tabla que es 1646,65 kg, dando una pequeña diferencia de 0,47 kg, equivalente a un error del 0,03%.
Bibliografía (1) ABBOT, M.M. y H. C. VAN NESS Thermodynamics (2ª ed.
) Nueva York, Mc Graw-Hill, Inc.
, 1989, 362 págs.
(2) COULSON, J.M. y J.F. RICHARDSON Chemical Engineering (Vol.
1) Oxford, Pergamon Press, 1990, 6 vols.
(3) GILES, R.V. Mecánica de los Fluidos e Hidráulica Nueva York, Mc Graw-Hill, 1969, 273 págs.
(4) HOUGEN, O.A., K.W. WATSON y R.A. RAGATZ Termodinámica Barcelona, Ed.
Reverté, S.A., 1964 (5) LESS, F.P. Loss Prevention in the Process Industries Londres, Butterworth & Co.
(Publishers) Ltd.
, 1980, 2 vols.
, 1316 págs.
(6) PERRY, R.H., y D. GREEN Perry's Chemical Engineers' Handbook (6ª Edición) Nueva York, Mc Graw-Hill, 1984.
(7) TNO Methods for the calculation of the physical effects of the escape of dangerous material (Liquids and gases) (1ª Edición) Voorburg (Holanda), Directorate-General of Labour, Ministry of Social Affairs, 1979.
(8) TNO Methods for the calculation of the physical effects.
Resulting from releases of hazardous materials (Liquid and gases) (2ª Edición) Voorburg (Holanda), Directorate-General of Labour, Ministry of Social affairs and Employment, 1992.
© INSHT