UNIVERSIDAD DE LOS ANDES

DEPARTAMENTO DE INGENIERÍA QUÍMICA

Marco de referencia para el análisis del riesgo asociado a eventos

Natech provocados por inundaciones.

Presentado por: Nicolás Villalba Hernández

Director
Felipe Muñoz Giraldo, PhD

Jurados
Felipe Salcedo, PhD
Mauricio Sánchez Silva, PhD

1
Resumen

En los últimos años, se han evidenciado cambios sustanciales en la variabilidad climática del planeta,
afectando significativamente los patrones de temperatura y precipitación, aumentando la
frecuencia y la severidad de eventos naturales de tipo meteorológico como son las inundaciones.
Este tipo de eventos naturales, suelen provocar accidentes mayores sobre la población que habita
zonas urbanas, así como a instalaciones industriales, áreas de almacenamiento de sustancias
químicas o cualquier infraestructura que albergue materiales peligrosos. Generando la pérdida de
contención de dichos materiales, la cual se puede transformar en incendios, explosiones y
dispersiones toxicas. Este tipo de eventos accidentales se conocen como Natech (Natural Hazard
Triggering Technological accidents). Los eventos Natech pueden acarrear pérdidas sociales,
ambientales y económicas. Por lo tanto, la necesidad de analizar los riesgos asociados al impacto de
una inundación sobre una actividad industrial, ha ganado gran importancia. Este trabajo presenta
una metodología para el análisis cuantitativo del riesgo asociado a eventos Natech debido al
impacto de inundaciones sobre tanques de almacenamiento verticales. A través de la cual es posible
establecer los escenarios accidentales finales, estimar la probabilidad de daño debido al pandeo de
las paredes del tanque, integrando el tratamiento de la incertidumbre asociada al modelo y evaluar
las consecuencias de una dispersión tóxica de sustancias peligrosas en fase líquida, utilizando una
aproximación económica. La metodología propuesta es aplicada a un caso de estudio, permitiendo
la estimación del riesgo asociado al impacto de diferentes tipos de inundaciones, sobre diferentes
configuraciones estructurales de tanques de almacenamiento. Esta aplicación permitió verificar la
aplicabilidad y repetitividad de la metodología como herramienta para el análisis cuantitativo del
riesgo. Adicional a esto, por medio del cálculo del riesgo asociado a los eventos Natech considerados
en este trabajo, se rectificó que la aproximación económica propuesta, busca representar mediante
valores monetarios, la magnitud de las afectaciones más no una pérdida económica y así relacionar
la frecuencia de los escenarios accidentales con la magnitud de los efectos asociados a la dispersión
de sustancias peligrosas en fase líquida durante un evento de inundación.

2
Contenido
1. Introducción ................................................................................................................................ 5
2. ¿Cómo se ha abordado la problemática de los eventos Natech? ............................................... 6
3. Metodología para el análisis de riesgos asociados a eventos Natech generados por
inundaciones ....................................................................................................................................... 9
3.1 Definición de los posibles escenarios accidentales ............................................................. 9
3.2 Estimación de frecuencias de los posibles escenarios finales........................................... 11
3.2.1 Cálculo de frecuencia: ocurrencia de una inundación .............................................. 12
3.2.2 Cálculo de la probabilidad de daño del tanque: Pandeo de las paredes debido al
impacto de una inundación....................................................................................................... 12
3.2.3 Modelo Probit para la estimación de la probabilidad de daño de tanques debido a la
deformación de sus paredes por el impacto de una inundación. ............................................. 14
3.2.4 Determinación de la probabilidad de ocurrencia de fallas estructurales debido al
impacto de una inundación....................................................................................................... 24
3.3 Estimación de los efectos de una pérdida de contención durante un evento de
inundación ..................................................................................................................................... 26
3.4 Evaluación de las consecuencias: Costos asociados a la pérdida de contención de
material peligroso debido al impacto de una inundación ............................................................ 28
3.5 Estimación del riesgo asociado a eventos Natech generado por inundaciones ............... 29
3.6 Evaluación del riesgo a partir de su calificación................................................................ 30
3.6.1 Matriz de riesgo para eventos Natech generados por inundaciones ....................... 32
4. Caso de estudio ......................................................................................................................... 34
4.1.1 Descripción del caso de estudio ................................................................................ 34
4.2 Resultados ......................................................................................................................... 35
5. Conclusiones.............................................................................................................................. 39
6. Referencias ................................................................................................................................ 41

3
Objetivos
Objetivo General

Evaluar el riesgo de actividades industriales que puedan ser objeto de eventos Natech por
escenarios de inundación.

Objetivos Específicos

• Plantear una aproximación cuantitativa que permita evaluar el riesgo en una instalación
industrial susceptible a eventos Natech por inundación.

• Identificar y aplicar métodos para la reducción de incertidumbre en el cálculo del riesgo

asociados a eventos Naetch por inundación.

• Realizar la aproximación propuesta sobre un caso de estudio.

4
1. Introducción
La variabilidad climática, representa cambios sustanciales en el clima de la tierra. Usualmente estos
cambios se presentan durante periodos de millones de años. A pesar de esto, últimamente se ha
presenciado que dichos cambios se presentan en intervalos de tiempo más cortos (cada 10 o 20
años), lo cual se considera una variabilidad climática anormal[1, 2]. Esta variabilidad anormal, se ve
representada en el aumento de la temperatura atmosférica, lo cual trae como consecuencia
cambios en los ecosistemas del planeta y aumentos en la severidad de los eventos de origen natural
[3]. Un evento de origen natural, se puede definir como: cualquier fenómeno natural que pueda
causar, pérdidas sociales y económicas, lesiones personales, daños a la propiedad o al medio
ambiente [4]. Dentro de este grupo de eventos se encuentran: tormentas tropicales, huracanes,
sequias, fuertes lluvias, deslizamientos, inundaciones, entre otros.

Cuando un evento de origen natural se materializa, pude transformase en un desastre natural,

interrumpiendo el funcionamiento de la sociedad generando pérdidas humanas, económicas y
ambientales. Las cuales, exceden las habilidades de la comunidad para responder a la situación[4].
Durante un desastre natural, pueden ocurrir desastres secundarios como pérdidas de contención
de material peligroso, consecuencia del impacto del evento natural sobre: instalaciones industriales,
áreas de almacenamiento de sustancias químicas o cualquier infraestructura que albergue
materiales peligrosos, ha aumentado [5]. Este tipo de eventos se conocen como Natech (por sus
siglas en inglés Natural Hazard Triggering Technological accidents)[6]. Diferentes estudios han
encontrado que, los eventos naturales que más se han registrado como causa de un accidente
industrial son sismos e inundaciones, seguidos por deslizamientos de suelos, huracanes y descargas
eléctricas, por último se encuentran las sequias[7, 8]. Además se ha encontrado que, las pérdidas
de material peligroso durante estos eventos, pueden transformase en incendios, explosiones o
dispersiones tóxicas sobre cuerpos de agua o el ambiente [9]. Dentro de las afectaciones producto
de un evento Natech, se pueden encontrar: afectaciones sociales, como daño a viviendas familiares,
pérdidas económicas, debido a los daños a propiedad ajena y los costos de remediación, así como
afectaciones ambientales, que involucran ecosistemas, cuerpos de agua, entre otros[10-12]. Vale la
pena recalcar que, las afectaciones generalizadas sobre el territorio, dificultan la atención y control
de los incidentes industriales, producidos por el evento de origen natural.

En los últimos años la ocurrencia de desastres naturales ha aumentado. Por ejemplo, en la década
de los noventas, se presentaron 1022 inundaciones en todo el mundo. Mientras que, en la última
década la ocurrencia de este tipo de eventos naturales, ha aumentado en un 74%[13]. Debido a esto
y las afectaciones asociadas a eventos Natech, se le ha prestado mayor atención al estudio de los
riesgos asociados a este tipo de eventos. Por lo tanto, el presente estudio propone una metodología
para el análisis cuantitativo del riesgo asociado a eventos Natech, debido al impacto de
inundaciones sobre tanques verticales de almacenamiento. A través de la cual es posible establecer
los escenarios accidentales finales, estimar la probabilidad de daño debido al pandeo de las paredes
del tanque, integrando el tratamiento de la incertidumbre asociada al modelo y evaluar las
consecuencias de una dispersión tóxica de sustancias peligrosas en fase líquida, utilizando una
aproximación económica. Por último, se presenta una matriz de riesgo, la cual se utiliza como
herramienta para la evaluación de los riesgos asociados a eventos Natech generados por
inundaciones.

5
2. ¿Cómo se ha abordado la problemática de los eventos Natech?
Debido al aumento en la ocurrencia de eventos Natech en los últimos años, diferentes autores han
desarrollado conocimiento para el análisis y gestión de los riesgos asociados a este tipo de eventos.
El estudio de eventos Natech y sus riesgos, puede realizarse a través dos aproximaciones: análisis
aposteriori y a priori. El primero, consiste en identificar información que permita caracterizar el
evento a partir del análisis de datos históricos y reportes de eventos pasados. Este tipo de análisis,
usualmente alimenta los métodos basados en la aproximación a priori. Por otro lado, el análisis a
priori de eventos accidentales, consiste en identificar los posibles escenarios accidentales y analizar
los riesgos que estos representan. Esta aproximación hace uso de metodologías sistemáticas, los
cuales permiten estimar los riesgos asociados a un incidente. Cada una de las etapas del marco de
referencia se debe alimentar con información proveniente de modelos que buscan representar las
situaciones que componen el escenario accidental.
Tabla 1. Estado del arte en el análisis a posteriori de eventos Natech.

Autor Año Tipo de desastre Natural Contribución Ref.

Identifica los eventos naturales más relevantes,
Fenómenos atmosféricos,
K. Rasmussen 1995 como causa de pérdida de contención de material [14]
geológicos y acuáticos.
peligroso.
Terremotos, Identifica de forma cualitativa las características
S. Young et al. 2002 inundaciones, erupciones de impacto y los efectos a la salud de diferentes [8]
volcánicas y sequias tipos de desastres naturales.
Presenta las lecciones aprendidas de uno de los
A. M. Cruz Huracanes, tormentas
2009 eventos Natech de mayor impacto en la historia, [10]
E. Krausmann tropicales e inundaciones.
provocados por el huracán Katrina.
Presenta un análisis de datos históricos, donde
caracteriza las inundaciones relacionadas con
V. Cozzani et accidentes tecnológicos. Identifica los equipos de
2010 Inundación. [15]
al. proceso involucrados en este tipo de eventos, así
como los posibles tipos de daño estructural
debido al impacto de la inundación.
Presenta un análisis de datos históricos,
identificando los equipos de proceso
E. Renni et al. 2010 Relámpagos. involucrados, las dinámicas de falla y los modos de [16]
daño, así como las consecuencias del evento final,
debido al impacto de relámpagos.
Presenta un análisis de datos históricos, donde
establece la severidad del terremoto en función
de los parámetros del pico de movimiento del
suelo. Adicional a esto se identificaron los equipos
E. Krausmann Terremotos, inundaciones
2011 más propensos a fallar en este tipo de situaciones, [17]
et al. y relámpagos.
así como los modos de daño que se pueden
presentar. Por último, se identificaron los
escenarios finales más frecuentes debido al
impacto de un terremoto.
Presenta un análisis de datos históricos dirigido a
accidentes que involucran líneas de transporte.
S. Girgin Meteorológicos,
2016 Identifica los costos asociados a un evento [18]
E. Krausmann geológicos e hidrológicos.
Natech, así como el volumen de material
derramado, con el objetivo de presentar la

6
 relevancia de los eventos Natech en la gestión del
riesgo.

Se realiza una revisión bibliográfica de los diferentes trabajos desarrollados entorno al análisis del
riesgo de eventos Natech. De los trabajos más relevantes, se realizó una categorización con base en
las aproximaciones presentadas anteriormente. De esta forma, la Tabla 1, presenta los trabajos
desarrollados en la última década, que utilizan la aproximación a posteriori. Mientras que en la Tabla
2, se encuentran las metodologías y modelos propuestos en la literatura para el análisis de eventos
Natech causados por diferentes tipos de eventos de origen natural.
Tabla 2. Estado del arte en el análisis a priori de eventos Natech

Autor Año Tipo de desastre Natural Contribución Ref.

Modelos para la estimación de daños
Presentan las funciones Probit para el
cálculo de la probabilidad de daño de
E. Salzano et al. 2003 Sismo [19]
tanques atmosféricos debido al impacto de
un sismo.
Proponen un modelo para el cálculo de la
probabilidad de daño, debido al pandeo de
G. Landucci et al. 2012 Inundación [20]
la pared de tanques verticales, a partir de
correlaciones empíricas.
Proponen un modelo para el cálculo de la
probabilidad de daño, debido al
G. Landucci et al. 2014 Inundación [21]
desplazamiento de vasijas horizontales, a
partir de correlaciones empíricas.
Presentan un conjunto de árboles de falla y
de eventos genéricos, para la identificación
C. El Hajj et al. 2015 Inundación [22]
de escenarios accidentales que involucran
el impacto de inundaciones
Metodología
Proponen una metodología para el cálculo
de riesgo asociado a eventos Natech,
G. Antononi et al. 2007 Sismo generados por sismos. Los equipos de [23]
interés en el estudio son tanques de
almacenamiento verticales y horizontales.
Proponen un marco de referencia para el
análisis del riesgo asociado a eventos
Natech generados por cualquier tipo de
v. Cozzani et al. 2014 Sismos e inundaciones [24]
desastre natural. Presenta elementos para
el estudio del impacto de sismos e
inundaciones.
Proponen una metodología para evaluar el
riesgo asociado a eventos Natech debido al
G. Antononi et al. 2015 Inundaciones impacto de inundaciones sobre tanques de [25]
almacenamiento. Únicamente evalúan los
escenarios finales de incendio y explosión.

7
 Proponen una metodología para evaluar el
riesgo asociado a eventos Natech debido al
impacto de relámpagos sobre tanques de
A. Necci et al. 2016 Descargas eléctricas almacenamiento. Presentan modelos para [26]
el cálculo de la probabilidad de daño de
tanques de almacenamiento debido al
impacto de un relámpago.

Una vez realizada la revisión bibliográfica, se identificó que existe una incertidumbre asociada a los
parámetros de entrada de los modelos propuestos para la estimación de la probabilidad de daño,
de los tanques debido al impacto de una inundación, lo cual indica que la representación del modelo
puede no ser la más acertada. En cuanto a las metodologías para el análisis del riesgo de eventos
Natech generados por inundaciones, se encontró que no hacen uso explícito de árboles de eventos
para la identificación de escenario finales. Lo cual no permite estructurar un escenario accidental
final completo, donde ciertos eventos secundarios entre el evento inicial y el final, pueden pasar
desapercibidos y generar una estimación de frecuencias alejada de la realidad. Por último, se
evidencia que, en las metodologías presentadas en la literatura, los análisis de riesgo se llevan a
cabo únicamente para dos escenarios finales: incendio y explosión, dejando de lado la dispersión
tóxica de sustancias en fase líquida, la cual resulta ser la consecuencia más recurrente, en eventos
Natech provocados por inundaciones.

Teniendo en cuenta los elementos enunciados anteriormente, a continuación, se presenta un marco

de referencia para el análisis del riesgo asociado a eventos Natech generados por inundaciones, por
medio del cual, se complementarán las metodologías propuestas por los autores presentados en la
Tabla 2. Por lo tanto, se realiza un tratamiento de la incertidumbre asociado a las variables de
entrada del modelo para cálculo de probabilidad de daño, debido al pandeo de las paredes del
tanque, a través de un marco de propagación puramente probabilístico, además, se presentan
arboles de eventos para identificar los posibles daños, fallas y escenarios finales debido al impacto
de la inundación sobre tanques de almacenamiento verticales. Por último, la metodología, se
enfocará en el escenario final de dispersión tóxica de sustancias en fase líquida.

8
3. Metodología para el análisis de riesgos asociados a eventos Natech generados por
inundaciones
La Figura 1, presenta la metodología propuesta para el análisis y evaluación de riesgos asociados a
accidentes tecnológicos provocados por un evento de inundación, la cual está compuesta por seis
(6) pasos, presentados a continuación: i) Definir los posibles escenarios accidentales debido al
impacto de la inundación, ii) Estimar la frecuencia del posible escenario final, iii) Evaluar las
consecuencias del evento final, iv) Estimar los efectos de las consecuencias identificadas, v) Estimar
el riesgo asociado a eventos Natech generados por inundaciones y vi) Evaluar el riesgo.

Figura 1. Metodología para el análisis y evaluación del riesgo asociado a eventos Natech generados por inundaciones

Cada uno de los pasos presentados en la Figura 1, estan compuestos por procedimientos, modelos
y herramientas, los cuales en conjunto con información que debe ser alimentada, permiten llevar a
cabo el análisis y la evaluación de riesgos. A continuación, se describirán cada uno de los pasos, con
el objetivo de establecer los elementos necesarios para desarrollar cada una de las etapas que
componen la metodología propuesta.

3.1 Definición de los posibles escenarios accidentales

Para poder definir los posibles escenarios accidentales finales (fase 1, Figura 1), que se pueden
desencadenar del al impacto de una inundación sobre un tanque de almacenamiento vertical, se
hace uso del método de árbol de eventos. La construcción de un árbol de eventos puede realizarse
siguiendo guías, como la propuesta por ARAMIS [27], adaptándola para abordar eventos Natech
provocados por inundaciones. En principio debe determinarse un evento iniciador, el cual suele
representar la falla de un componente o una falla externa. Posterior a esto, se define la secuencia
de eventos, condiciones o funciones de seguridad que siguen al evento iniciador y que deben
sobrepasarse para obtener cierto resultado. Este tipo de eventos se denominan eventos
secundarios (ES), mientras que los eventos obtenidos al final de cada rama se conocen como eventos

9
accidentales finales (EF)[28]. La Figura 2 muestra una representación de un árbol de eventos y sus
elementos.

Figura 2. Representación de un árbol de eventos[28].

Considerando lo anterior, el evento iniciador seleccionado es el impacto frontal de la inundación

sobre un tanque de almacenamiento. Por lo cual, se seleccionaron la altura y la velocidad de la
inundación como los parámetros que permitirán caracterizar los impactos que pueden generar las
inundaciones [7, 15, 20]. Debido a esto, se han establecido tres (3) tipos de frentes de inundación,
los cuales son: 1. Velocidad baja, este tipo de inundaciones se caracterizan por tener una velocidad
del frente de agua de magnitudes cercanas a 0 m/s. 2. Velocidad media, este tipo comprende las
inundaciones con velocidades iguales o menores a 1 m/s. 3. Velocidad alta, el cual enmarca las
inundaciones cuya velocidad es mayor a 1 m/s [9]. En cuanto a los eventos secundarios (ES), su
selección, depende de las consecuencias que acarrea el impacto de la inundación sobre un tanque
vertical. Cozzani et al.[15] propone un análisis de datos histórico donde identificó diferentes tipos
de daños estructurales que pueden sufrir los tanques de almacenamiento durante un evento de
inundación. Con base en este análisis se identificaron tres tipos de eventos secundarios, los modos
de daño, los tipos de falla debido al daño y los tipos de pérdida de contención. Para el primero
tenemos que, el tanque puede sufrir daños por: desplazamiento del equipo o pandeo se sus
paredes. Como consecuencia del daño pueden presentarse cinco (5) tipos diferentes de fallas, las
cuales son: 1. Colapso de la estructura, 2. Falla total de la conexión, 3. Falla parcial de la conexión,
4. Falla del techo del tanque y 5. Ruptura de la coraza. Una vez se presenta alguna de las posibles
fallas, ocurre la pérdida de contención de material peligroso, para este evento secundario se
establecieron tres (3) modos de liberación, los cuales harán parte de la fase de estimación de
volumen derramado y dependen de la tipología de falla que se presenta. Estos modos de liberación
se presentarán de forma detallada en la Sección 3.4. En la Figura 3, se presenta el árbol de eventos
genérico para la identificación de escenarios accidentales finales (sucesos) asociados al impacto del
frente de inundación, que involucran los eventos presentados anteriormente. En esta misma Figura,
se indican con las líneas en negrilla, los escenarios de interés para este estudio, en un trabajo futuro
se desarrollarán los demás escenarios.

10
 Figura 3. Árbol de eventos para la secuencia de eventos debido al impacto de diferentes tipos de inundación a tanques
de almacenamiento verticales.

Una vez identificado el escenario final, es necesario estimar la frecuencia de ocurrencia del suceso,
para continuar con el proceso de análisis del riesgo.

3.2 Estimación de frecuencias de los posibles escenarios finales

Esta sección, presenta el procedimiento para la estimación de frecuencias de los posibles escenarios
finales (Figura 1, fase 2). El cálculo de la frecuencia de eventos finales, hace uso de dos operaciones,
la intersección de eventos independientes y la unión de los eventos en las diferentes ramas que
terminan en un mismo evento accidental final.

11
Cálculo de la intersección de eventos (Ecuación 1)
(1)
𝑝𝑗 = ∏ 𝑝𝑖
𝑖∈𝑗

Cálculo de la unión de eventos (Ecuación 2)

(2)
𝑝𝐸𝐹 = ∑ 𝑝𝑗
𝑗=𝐸𝐹

Donde i se define como cada uno de los eventos secundarios como y j se define como cada una de
las ramas del árbol. Con base en esto, la frecuencia de cada uno de los eventos finales (EF) será la
multiplicación de la frecuencia del evento iniciador 𝑓𝐸𝐼 por la probabilidad del evento final 𝑝𝐸𝐹 .
Como se puede ver en la Figura 2 y en las ecuaciones presentadas, la estimación de la frecuencia
del escenario final, depende de la frecuencia del evento iniciador y la probabilidad de ocurrencia de
cada uno de los ES.

La determinación de estas frecuencias y probabilidades puede realizarse a través de registros

históricos, datos operacionales o modelos que representen el fenómeno que se está estudiando
[29]. Para los eventos señalados en la Figura 3 con la línea en negrilla, este trabajo, propone
diferentes aproximaciones para el cálculo de la frecuencia o probabilidad de ocurrencia
dependiendo del evento, las cuales serán presentadas a continuación.

3.2.1 Cálculo de frecuencia: ocurrencia de una inundación

Uno de los parámetros más utilizados al momento de estudiar la ocurrencia de una inundación
dentro de la gestión del riesgo de desastres naturales, es el periodo de retorno del evento natural
(𝑡𝑟 ). Los valores para este parámetro, usualmente se encuentra disponibles en la literatura. Esto,
debido a que algunas entidades gubernamentales han calculado la tasa de retorno para
inundaciones, que se pueden presentar en zonas de alto riesgo [30, 31], en caso de no estar
disponible, su estimación es posible por medio de modelos hidrológicos[32, 33]. La frecuencia de
ocurrencia de una inundación en función de su periodo de retorno se calcula a partir de la Ecuación
3.
1 (3)
𝑓=
𝑡𝑟
Una vez establecido el cálculo de la frecuencia del evento iniciador, presentaremos el modelo
propuesto para la estimación de la probabilidad de ocurrencia de uno de los modos de daño
presentados en la Figura 3, daño del tanque debido al pandeo o deformación de sus paredes.

3.2.2 Cálculo de la probabilidad de daño del tanque: Pandeo de las paredes debido al impacto de
una inundación
El pandeo o deformación, es una de las posibles consecuencias producto del impacto de una
inundación sobre tanques de almacenamiento verticales (Figura 4). La determinación de la
probabilidad de ocurrencia de este evento, es necesaria para el cálculo de la frecuencia del
escenario final. Para estimar la probabilidad de daño debido al pandeo de tanques verticales,
Landucci et al. [20] proponen un conjunto de relaciones matemáticas, derivadas de un modelo

12
mecánico, que describe el sistema presentado en la Figura 4. Durante el desarrollo de las
correlaciones Landucci et al. asumen valores constantes para ciertos parámetros de entrada, siendo
que su comportamiento natural presenta variabilidad. Esto genera una incertidumbre asociada a
dichos parámetros, lo que se traduce en una representación de la probabilidad de daño poco
representativa de la realidad. La incertidumbre presente en los modelos, se encuentra asociada a la
falta de conocimiento que existe del comportamiento de un sistema. Lo que resulta problemático,
cuando la información proporcionada por el modelo debe utilizarse en un proceso de toma de
decisiones [34].

Figura 4. Representación gráfica del daño por paneo de un tanque de almacenamiento, debido al impacto de un frente
de inundación.

Para poder considerar la incertidumbre en un modelo, es necesario conocer qué tipo de

incertidumbre se está presentando. Parry[35], en su estudio, identifica dos tipos de incertidumbres:
incertidumbre aleatoria e incertidumbre epistémica. La primera, es aquella que se presenta cuando
se realiza el modelaje de sistemas cuyas variables poseen una variación natural, por lo tanto, el
conocimiento de su comportamiento es inexacto[36] y suele estar ligada a las variables de entrada
del modelo. Por otro lado, la incertidumbre epistémica, depende de que tan bien el modelo
representa el sistema que se encuentra describiendo, sabiendo que el conocimiento relacionado
con las variables y parámetros que lo constituyen, es pobre[35]. Una vez establecido el tipo de
incertidumbre en el modelo, esta se debe tratar con el fin de conocer su efecto en la respuesta, lo
cual permitirá obtener resultados más cercanos a la realidad. Para esto existen diferentes marcos
de propagación, de los cuales solo se mencionarán dos, el marco puramente probabilístico y el
marco puramente posibilístico[37]. El primer marco considera que las variables de entrada poseen
incertidumbre aleatoria, la cual puede ser representada por medio de funciones de probabilidad,
los métodos utilizados para la propagación de esta incertidumbre son: expansiones en series de
Taylor de la función que describe la respuesta del sistema y métodos de muestreo como la
simulación de Monte Carlo[34]. Mientras que el marco posibilístico, considera la incertidumbre
epistémica presente en las variables de entrada del modelo y lleva a cabo la representación de estas,
mediante distribuciones de posibilidad, esta propagación se desarrolla a través del principio de
extensión de la teoría de conjuntos difusos[37, 38].

Este trabajo propone un modelo Probit para la estimación de la probabilidad de daño de un tanque
debido al pandeo de las paredes, consecuencia del impacto de una inundación, el cual integra la

13
incertidumbre aleatoria asociada a las variables de entrada del modelo, para el tratamiento de la
incertidumbre se hará uso del método de simulación de Mote Carlo.

3.2.3 Modelo Probit para la estimación de la probabilidad de daño de tanques debido a la

deformación de sus paredes por el impacto de una inundación.
Este numeral, presenta el proceso para la construcción del modelo Probit para el cálculo de la
probabilidad de daño de tanques de almacenamiento, debido al impacto de una inundación,
específicamente para el modo de daño de pandeo de las paredes del tanque. La probabilidad estará
determinada por la magnitud de la inundación en términos de la altura y la velocidad, la
configuración estructural del reservorio, la condición de operación (nivel de llenado) en el momento
de la inundación y, por último, la densidad del fluido que almacena.

Para la construcción de las funciones Probit, se propone el siguiente algoritmo, el cual está
compuesta por cuatro (4) pasos secuenciales. En el primer paso, se lleva a cabo una revisión
bibliográfica, en busca de un modelo que representara el sistema de interés. Paso seguido, a partir
del modelo seleccionado, se establece la función de estado límite, la cual debe representar la
respuesta del sistema bajo una solicitación. En el tercer paso, se calculará la probabilidad de daño
para el sistema de interés, integrando la incertidumbre asociada a los parámetros de entrada
utilizando simulaciones de Monte Carlo. El último paso del algoritmo consiste en ajustar las
probabilidades calculadas, a una función Probit logarítmica, la cual permite estimar la probabilidad
de daño debido al pandeo de las paredes del tanque. El diagrama que explica el algoritmo se
presenta en la Figura 5.

Figura 5. Metodología para la construcción del modelo Probit para la estimación de la probabilidad de daño de tanques
debido a la deformación de sus paredes por el impacto de una inundación.

14
En los siguientes numerales, se presentan de forma detalla los procedimientos desarrollados en
cada uno de los pasos que componen el algoritmo para la construcción de funciones Probit, con el
fin de presentar las funciones para el cálculo de la probabilidad de daño mencionadas
anteriormente.

[Link] Modelo para la determinación de la deformación de un tanque de almacenamiento

debido al impacto de una inundación
El modelo que se presenta a continuación, permite estimar, si el equipo sufrirá daño por pandeo o
deformación de sus paredes Figura 4. El modelo propuesto por Landucci et al. [20], relaciona las
presiones que actúan sobre la pared de un tanque de almacenamiento vertical, al ser impactado por
na inundación con una velocidad de frente de agua (𝑣𝑤 ) y altura de inundación (ℎ𝑤 ), con la
resistencia que opone el tanque debido a su estructura y fluido almacenado.

Figura [Link] de las presiones que actúan sobre la pared de un tanque de almacenamiento debido al impacto de un
frente de inundación [19].

Como se puede ver en la Figura 6, la inundación genera una presión externa sobre el tanque
designada como carga externa (𝑃𝑤 ), la cual está compuesta por: la presión estática máxima (𝑃𝑤𝑠 ) y
la presión dinámica máxima (𝑃𝑤𝑑 ). La carga externa se puede calcular de la siguiente manera:

𝑃𝑤 = 𝑃𝑤𝑠 + 𝑃𝑤𝑑 (4)

Donde 𝑃𝑤𝑠 y 𝑃𝑤𝑑 se definen como:

𝑃𝑤𝑠 = 𝜌𝑤 𝑔ℎ𝑤 (5)

1 (6)
𝑃𝑤𝑑 = 𝜌𝑤 𝑘𝑤 𝑣𝑤2
2
Por otro lado, existen presiones que producen resistencia a los efectos de la carga externa, como se
puede ver en la Figura 6, la presión del fluido (𝑃𝑓 , Ecuación 7) cumple con este objetivo y se calcula
en función de la densidad del fluido almacenado, la altura del tanque y el nivel de llenado. Siendo

15
el modo de daño de interés el pandeo de las paredes del tanque, no se puede pasar por alto el hecho
de que, la estructura está diseñada para resistir una carga mínima, antes de presentar una
deformación[39]. Dicha carga, se conoce como presión critica (𝑃𝑐𝑟 , Ecuación 8), la cual se puede
estimar a partir de las propiedades mecánicas del material del tanque. Al igual que la carga externa,
la combinación de la presión crítica y la presión del fluido, componen la presión de resistencia del
tanque (𝑃𝑟 ), las ecuaciones para el cálculo de las presiones se presentan a continuación:

𝑃𝑓 = 𝜌𝑓 𝑔𝐻𝜙 (7)

2𝐸𝑡 1 𝑡2 2
2𝑛2 − 1 − 𝜈
𝑃𝑐𝑟 = + (𝑛 − 1 + ) (8)
𝐷
2𝑛𝐻 2
3𝐷 2 (1 − 𝜈 2 ) 2𝑛𝐻 2
( 𝑛2 − 1) (1 + ( ) 1 + ( 𝜋𝐷 )
𝜋𝐷 )
( )
𝑃𝑟 = 𝑃𝑓 + 𝑃𝑐𝑟 (9)

Con este modelo, es posible establecer una relación entre las cargas que actúan sobre el tanque al
momento de ser afectado por una inundación, para determinar si el equipo sufrirá daño por el
pandeo o deformación de sus paredes. La determinación de dicha relación se aborda en la siguiente
sección.

[Link] Determinación de la ecuación de estado límite

En esta sección, se establece una relación para definir si el tanque resistirá la carga externa,
producto del impacto del frente de inundación o si presentará daño debido al pandeo. Para esto, se
utiliza el problema básico de confiabilidad, el cual permite definir un criterio de daño y por
consiguiente una probabilidad de daño. El problema básico de confiabilidad considera la relación
entre la resistencia del sistema R y la solicitación S a la que este se encuentra sometido (Ecuación
10). Debido a la naturaleza aleatoria de las variables, la probabilidad de daño puede definirse a
través de la Ecuación 11.

𝑝𝑑 = 𝑝(𝑅 − 𝑆 < 0) (10)

𝑝𝑑 = 𝑝(𝑔(𝑅, 𝑆) ≤ 0) (11)

Donde la función 𝑔(𝑅, 𝑆) se conoce como función de estado límite del sistema[34]. Ahora,
considerando el sistema de interés para este trabajo, se puede definir la carga externa generada por
el impacto de la inundación como la solicitación y la presión de resistencia (𝑃𝑟 ) como la resistencia
del sistema. Así posible establecer la función de estado límite (Ecuación 12) la cual permite conocer
la probabilidad de daño por pandeo, en términos de las cargas (𝑃𝑤 ) que actúan sobre el tanque de
almacenamiento impactado por un frente de inundación y las presiones que ejercen resistencia (𝑃𝑟 ).

𝑔(𝑃𝑟 , 𝑃𝑤 ) = 𝑃𝑟 − 𝑃𝑤 (12)

Por lo tanto, el criterio para el cálculo de la probabilidad de daño debido a la deformación de las
paredes del tanque, consecuencia de la carga dinámica y estática generadas por el impacto del
frente de inundación, tendrá la siguiente forma

𝑝𝑑 = 𝑝(𝑃𝑟 − 𝑃𝑤 ≤ 0) (13)

16
Una vez establecida la función de estado límite, es posible aplicar el método de simulación de Monte
Carlo. La aplicación de este método, se presenta a continuación.

[Link] Tratamiento de la incertidumbre dentro del cálculo de la probabilidad de daño debido al

pandeo de tanques verticales y construcción de datos sintéticos
En esta etapa (paso 3, Figura 5), se calcula la probabilidad de deformación o paneo de un tanque de
almacenamiento vertical, consecuencia del impacto de una inundación. Dicha estimación se realiza,
haciendo uso de la simulación de Monte Carlo con el fin de tratar la incertidumbre asociada a los
parámetros del modelo.

En la Figura 7, se presenta el diagrama que explica la metodología para el cálculo de la probabilidad

de daño debido al pandeo del tanque, integrando el tratamiento de la incertidumbre asociada a los
parámetros de entrada del modelo.

Figura 7. Metodología para el cálculo de la probabilidad de daño, debido al pandeo de las paredes del tanque con
tratamiento de la incertidumbre dentro de un marco de propagación puramente probabilístico.

17
Para iniciar el proceso presentado en la Figura 7, se deben establecer los parámetros que presentan
variabilidad. Puntualmente, Landucci et al. [20] en su estudio establece que para el desarrollo del
modelo se estimaron valores puntuales para el coeficiente hidrodinámico del agua (𝑘𝑤 ) y la
densidad del agua (𝜌𝑤 ), conociendo que estos parámetros presentan una variabilidad en su
naturaleza. Con base en esto, dichos parámetros se seleccionaron como los parámetros aleatorios
del modelo. En la se presentan los parámetros estadísticos para 𝑘𝑤 y 𝜌𝑤 .
Tabla 3. Parámetros estadísticos de las distribuciones asociadas a los parámetros aleatorios seleccionados[40, 41].

Variable Unidad Tipo de distribución Media (𝜇) Coeficiente de variación (𝐶𝑉)

𝑘𝑤 - Normal 1,8 1,1%
𝜌𝑤 𝑘𝑔/𝑚3 Normal 1100 9,1%

Una vez establecidos los parámetros aleatorios, se deben seleccionar los valores que se le asignaran
a las diferentes variables que componen el modelo. Como se presentó en la Sección [Link], el
modelo está especificado en función del vector de impacto, el nivel de llenado, la geometría del
tanque, las propiedades mecánicas del material del tanque y la densidad del fluido almacenado. Se
seleccionaron diferentes valores para algunos de estos parámetros, así como diferentes
configuraciones estructurales de reservorios, para asegurar una estimación completa de la
probabilidad de daño, ya que esto, permite que se integran diferentes escenarios de inundación y
de operación durante el cálculo de la probabilidad. Para la configuración de los tanques se utilizó
una base de datos, la cual contiene 116 diferentes geometrías verticales atmosféricas (Tabla 9 y 11
en Anexos), en cuanto a las demás variables, los rangos y valores se presentan en la Tabla 4.
Tabla 4. Valores y rangos para el desarrollo de la simulación de Monte Carlo

Variable Unidad Valores Referencia

𝑣𝑤 m/s [0,5 – 4]
[20, 24]
ℎ𝑤 m [0,5 – 4]
𝜙 % [10 – 70] [20, 21]
𝐸 Pa 196x109 [42]
𝜈 - 0,294 [42]
𝜌𝑓 kg/m3 760 - 1100 [43]

Para establecer los diferentes escenarios de inundación, sobre los que se quiere estimar la
probabilidad de daño debido al pandeo del tanque, el rango de velocidad de inundación se divide
en rangos de menor tamaño, con el fin de facilitar el tratamiento de los resultados. Los rangos se
asignaron de la siguiente manera, [0,5-1], [1-1,5], [1,5-2], [2-2,5], [2,5-3], [3-3,5] y [3,5-4] m/s. Para
cada uno de estos intervalos, junto con el de altura de la inundación, se generan 100 elementos,
mientras que el nivel de llenado se dividió en 7 elementos para llevar a cabo la simulación. En cuanto
al material de construcción del tanque se seleccionó acero inoxidable 316 y de allí se derivan sus
parámetros mecánicos. En cuanto a la densidad del fluido almacenado, se seleccionan dos valores,
relacionados con las sustancias que más han estado involucradas en eventos Natech, como: crudos
pesados y gasolina[15]. A continuación, se presenta la ecuación para el cálculo de la probabilidad
de daño por pandeo de las paredes de un tanque

18
 𝑁𝑠𝑖𝑚 (14)
𝑔(𝑖) 1 𝑖𝑓 𝑃𝑟 − 𝑃𝑤 ≤ 0
𝑝𝑑 (ℎ𝑤 , 𝑣𝑤 , 𝜙) = ∑ 𝑦 𝑔(𝑖) = {
𝑁𝑠𝑖𝑚 0 𝑑. 𝑙. 𝑐
𝑖=1

El cálculo de la probabilidad de daño debido al pandeo del tanque provocado por el impacto de una
inundación, se llevó a cabo siguiendo la metodología para el tratamiento de incertidumbre
presentada en la Figura 7. La programación del algoritmo se realizó en Matlab® con 100000
iteraciones dentro de la simulación de Monte Carlo. El tiempo promedio de simulación para un valor
de nivel de llenado, un rango de velocidad de inundación y las 116 configuraciones es de 7,3 horas
en un servidor con dos procesadores Intel® Xeon® de 3,47GHz cada uno.

Después de realizada la simulación de Monte Carlo, se obtuvieron las probabilidades de daño debido
a la deformación de las paredes del tanque, para cada una de las configuraciones estructurales
seleccionadas, cada uno de los niveles de llenado y cada combinación del vector de impacto de la
inundación. La Figura 8, presenta la distribución de probabilidad obtenida de la simulación de Monte
Carlo, los perfiles que se presentan, hacen referencia a la probabilidad de daño por pandeo de un
tanque de geometría ID 113, nivel de llenado 50%, densidad del fluido almacenado 760 kg/m 3. La
probabilidad de daño se presenta en función de un rango de altura y cuatro valores de velocidad de
la inundación. De esta figura, se puede observar como a medida que la inundación se hace más
profunda, la probabilidad de que el tanque sufra daño por pandeo de sus paredes aumenta.

Figura 8. Comparación de la distribución de probabilidad de daño debido al pandeo del tanque, consecuencia del
impacto de una inundación, para diferentes velocidades de inundación.

19
Ahora, considerando los cambios en la velocidad de la inundación, se puede ver como un aumento
en esta variable desplaza el perfil de probabilidad hacia alturas de inundación menores. Esto
significa que un aumento en la velocidad de la inundación, a un valor fijo de altura, aumentará la
probabilidad de daño, lo cual es consistente con el sistema ya que se espera que a medida que el
impacto de la inundación sea mayor, más probable será que el tanque se deforme.

Figura 9. Comparación de la distribución de probabilidad de daño debido al pandeo del tanque, impactado por una
inundación de alta velocidad.

En la Figura 9, se presenta una comparación de las diferentes distribuciones de probabilidad

acumuladas, calculadas en función de diferentes condiciones de operación del tanque. En esta figura
Se puede apreciar como a medida que se aumenta el nivel de llenado del tanque, la curva de
probabilidad acumulada se desplaza hacia valores más altos de altura de la inundación. Esto nos
permite interpretar que un tanque con mayor cantidad de material almacenado, va a tener más
oportunidad de sobrellevar el posible daño generado por inundaciones de bajas alturas, contrario a
los recipientes con menor inventario, cuyo perfil presenta una mayor vulnerabilidad al impacto de
inundaciones no tan profundas.

Las curvas de probabilidad acumulada obtenidas con la simulación de Monte Carlo, permiten
estimar un valor puntual de probabilidad de daño debido a la deformación del tanque, para los
valores presentados en la Tabla 4. A pesar de esto, dicha información (la cual llamaremos datos
sintéticos), no es de uso versátil y está restringida a los valores seleccionados para su desarrollo. Por
lo tanto, es necesario valorizarla y para esto se construyeron un conjunto de funciones Probit con
base en los datos sintéticos calculados. El desarrollo de dichas funciones se presenta en la siguiente
sección.

20
[Link] Funciones Probit para la estimación de la probabilidad de daño de tanques de
almacenamiento debido al pandeo de sus paredes
Esta sección, presenta el último paso de la metodología presentada para la construcción de las
funciones Probit para la estimación de la probabilidad de daño debido a la deformación del tanque.
Como se puede observar en las figuras 7 y 8, el comportamiento que presentan las funciones de
probabilidad acumulada, se asemeja a la de las curvas de dosis respuesta. Las curvas de dosis
respuesta, son funciones que buscan estimar la probabilidad de afectación de un ser vivo que ha
sido expuesta a un dosis[44]. A pesar de su uso común en la toxicología, este mismo tipo de curvas
ha sido aplicado a otro tipo de exposiciones como la radiación de calor o sobrepresiones, donde se
busca estimar la probabilidad de que existan daños materiales (entiéndase daños materiales como,
daños a edificios, casas, autos, etc) debido a este tipo de exposiciones[45].

Las curvas de dosis respuesta, suelen ser funciones de probabilidad acumulada con un
comportamiento de tipo sigmoidea. Con el objetivo de establecer una función matemática que
permita estimar la probabilidad de afectación dependiendo de la dosis, se han desarrollado diversos
métodos para la representación de dichas curvas, la mayoría busca establecer una función de línea
recta[44]. De los métodos más utilizados, se encuentra el método Probit. Este método consiste en
utilizar la relación entre la probabilidad de afectación y los puntos Probit (Ecuación 15), para
establecer una función lineal que permita la estimación de los puntos Probit en función de la dosis.
Para determinar dicha función, se hace uso de técnicas estándar para el ajuste de curvas.

𝑌−5 |𝑌 − 5| (15)
𝑝̂ = 0.5 (1 + ) erf ( )
|𝑌 − 5| √2
Con base en lo anterior, se puede establecer que las funciones de probabilidad de daño calculadas
en la sección anterior, tienen el mismo comportamiento de una curva dosis respuesta. Siendo el
vector de impacto de la inundación la dosis y la respuesta, que tan probable es que el tanque sufra
daño por el pandeo de sus paredes, consecuencia del impacto de dicha inundación. Por lo tanto, se
aplicará el método Probit, para obtener un conjunto de expresiones matemáticas que permitan
estimar la probabilidad de daño debido a la deformación del tanque, en función de la configuración
estructural del tanque, la velocidad y altura de la inundación.

Para la aplicación del método Probit, se propone una función logarítmica la cual relaciona los puntos
Probit (𝑌) con el vector de impacto de la inundación (𝑣𝑤 y ℎ𝑤 ), la cual se presenta a continuación

𝑌 = 𝑘1 + 𝑘2 ln(𝑣𝑤 ) + 𝑘3 ln(ℎ𝑤 ) (16)

Donde 𝑘1 , 𝑘2 y 𝑘3 son las constantes del modelo. Cabe recordar que las constantes Probit del
modelo, se estiman para fluidos almacenados con densidad en fase líquida cercanos a los 760 kg/m3
y a los 1100 kg/m3. Para poder estimar las constantes del modelo que mejor ajusten los datos a una
línea recta, se utilizó el método de ajuste por mínimos cuadrados, el cual busca que las
probabilidades estimadas con la Ecuación 16, sean iguales a los datos sintéticos, maximizando el
coeficiente de determinación de la siguiente manera
𝐵 (17)
𝑅𝑖2
𝐹. 𝑂 ≡ max {∑ }
𝐵
𝑖=1

21
Donde 𝐵 es el tamaño de cada intervalo de velocidad de inundación 𝑣𝑤 y 𝑅 2es el coeficiente de
determinación, el cual se calcula como
𝑁 2 (18)
𝑝
∑𝑗=1 (𝑝𝑓𝑗 − 𝑝̂𝑗 )
2
𝑅 =1− 2
𝑁
𝑝
∑𝑗=1 (𝑝𝑓𝑗 − ̅̅̅)
𝑝𝑓

Donde 𝑁𝑝 es el número de datos sintéticos generados para cada uno de los intervalos de velocidad
de inundación. Así se obtendrá un conjunto de constantes Probit, para cada curva de probabilidad
acumulada, resultado de la simulación de Monte Carlo.

Una vez realizado el ajuste de los datos sintéticos, se calcularon los coeficientes de correlación entre
dichos datos y el modelo Probit propuesto, encontrando que la correlación entre estos, registra
valores alrededor del 98%. Lo cual se puede evidenciar en la Figura 10, la cual presenta las curvas
de Probabilidad de daño en función del vector de impacto, para los datos sintéticos y la función
Probit propuesta. Se puede observar que las curvas se superponen, lo que confirma los niveles de
ajuste mencionados anteriormente.

Figura 10. Comparación datos sintéticos contra probabilidad de daño por pandeo, estimadas mediante el modelo Probit
propuesto, para un tanque de ID 93, velocidad de inundación 2 m/s y densidad de fluido almacenado de 760 kg/m3.

Realizado el ajuste y determinar los coeficientes para la función Probit, es estableció que no era
suficiente calcular las constantes Probit para cada uno de los rangos de velocidad y nivel de llenado.
Esto, debido a que cada configuración estructural tendrá una gran cantidad de conjuntos de
constantes Probit. Para reducir la cantidad de tripletas de constantes, se planteó realizar otro ajuste
de datos, esta vez en busca de relacionar las constantes de la función Probit, con la geometría del
tanque.

Lo anterior no solo permitirá reducir la cantidad de información para el uso de las funciones, sino
que permitirá que el modelo Probit propuesto se aplique a cualquier tanque de almacenamiento

22
vertical. Por lo tanto, se propone una función cuadrática (Ecuación 19), la cual calcula cada uno de
los parámetros Probit, en función de la geometría del tanque.

𝑘𝑖 = 𝑎𝑖 + 𝑏𝑖 𝑥 + 𝑐𝑖 𝑥 2 | 𝑖 = 1, 2, 3 (19)

Donde 𝑥 se define como el factor geométrico y se expresa de la siguiente forma

𝐻𝑡 (20)
𝑥=
𝐷
El factor geométrico, presenta una relación entre los parámetros de la configuración estructural del
tanque. La cual, permitirá llevar a cabo el ajuste entre la geometría del reservorio y las constantes
Probit. Para el ajuste se hace uso de los conjuntos de las tripletas de constantes calculadas y la base
de datos de configuración de tanques presentada en la Tabla 8. El método de mínimos cuadrados
se aplica nuevamente para estimar las constantes 𝑎, 𝑏 y 𝑐.

La Figura 11, presenta los resultados obtenidos para el ajuste del factor geométrico presentado en
la Ecuación. 20 y las constantes de la función Probit (𝑘1 , 𝑘2 y 𝑘3 ), para un tipo de configuración ID
93. En la figura se presentan los coeficientes de correlación para el ajuste, se puede ver que dichos
coeficientes están alrededor del 94%, lo cual indica que la función propuesta para relacionar la
geometría del tanque con cada una de las constantes Probit, ajusta correctamente los datos.

Figura 11. Ajuste de las constantes Probit (𝑘1 , 𝑘2 y 𝑘3 ) con el factor geométrico.

23
Retomando lo expuesto en esta sección, se realizan los ajustes a cada una de las funciones
propuestas para obtener las constantes que permitirán el cálculo de la probabilidad de daño debido
al pandeo de un tanque, consecuencia del impacto de una inundación. En Anexos, se presentan las
constantes para la estimación de la probabilidad de daño por pandeo, para sustancias almacenadas
con densidad de 760 kg/m3 (Anexos Tabla 9) y 1100 kg/m3 (Anexos Tabla 10).

Como se evidencia durante el desarrollo de las funciones Probit, se llevó a cabo un proceso de
tratamiento de la incertidumbre asociada a los parámetros de entrada del modelo, que permite la
estimación de daño por el pandeo de las paredes del tanque, debido al impacto de una inundación,
con el fin de lograr la estimación de la probabilidad de daño, integrando el comportamiento natural
de dichos parámetros. Lo anterior, mejorará los resultados del análisis del riesgo asociado a eventos
Natech, provocados por inundaciones, debido a que los valores calculados, son información de
entrada para los pasos siguientes del marco propuesto en la Figura 1. Por último, cabe recalcar que
las funciones propuestas, son una herramienta que puede ser aplicada dentro del marco para el
análisis de riesgos propuesto en este trabajo. Otras metodologías pueden hacer uso de la
herramienta, siempre y cuando involucre las variables consideradas en las funciones Probit para el
cálculo de la probabilidad de daño.

A continuación, se presentará la determinación de la probabilidad de ocurrencia de las posibles

fallas, producto del daño que puede sufrir el tanque al ser impactado por una inundación.

3.2.4 Determinación de la probabilidad de ocurrencia de fallas estructurales debido al impacto

de una inundación
La determinación de la probabilidad de ocurrencia de las fallas debido al impacto de una inundación
sobre un tanque de almacenamiento vertical, se desarrolla a partir de la aproximación
determinística, asignando valores a los eventos que se identificaron como posibles fallas del
sistema. Esta asignación, se realizó con base en un análisis de accidentes industriales provocados
por inundaciones. Donde se establece que, dentro de las posibles fallas, la que presenta mayores
registros es falla en las conexiones, seguido por la falla del techo del tanque y el colapso de la
estructura. Con esta información y los elementos presentados en la Figura 2, se asignan las
probabilidades de las diferentes fallas que se presentan en el árbol de eventos para el impacto de
una inundación (Figura 12). Como se puede ver en la Figura 12, la probabilidad de los eventos falla
de una misma rama suman 1, lo cual es consistente con las características de un árbol de eventos,
presentadas en la Figura 2. En trabajo futuro, se buscará tratar la incertidumbre asociada a los
eventos de falla, consecuencia del impacto de una inundación sobre tanques de almacenamiento
verticales.

Con esto, se finaliza la etapa de estimación de frecuencias de los eventos finales y se procede a
evaluar las consecuencias del escenario accidental final.

24
 25
Figura 12. Arboles de eventos para la identificación de la secuencia de eventos debido al impacto de una inundación en función de su velocidad. Las
probabilidades de ocurrencia de los tipos de falla se presentan en cada una de sus ramas.
3.3 Estimación de los efectos de una pérdida de contención durante un evento de inundación
A partir de una pérdida de contención de sustancias peligrosas, se pueden generar diferentes
escenarios como: explosión, incendio y dispersión tóxica. Los cuales tendrán un impacto específico,
de llegar a materializarse[29]. Dentro el análisis de riesgo convencional, se han desarrollado
diferentes modelos para el cálculo de efectos de las posibles consecuencias derivadas de una
pérdida de contención. Los efectos suelen estar asociados a sobrepresiones consecuencia de
explosiones, a la radiación de calor producto de los diferentes tipos de incendio y fuegos y, por
último, niveles de toxicidad los cuales son generados a partir de la dispersión de sustancias
tóxicas[45].

En cuanto la estimación de los efectos desarrollados en metodologías para el análisis de riesgo de

eventos Natech, suele hacerse uso de los modelos convencionales, debido al enfoque que se le da
al análisis sobre los escenarios de incendio y explosión[25, 46]. Este trabajo, por el contrario, se
enfoca en analizar los riesgos relacionados con el escenario de dispersión tóxica de sustancias en
estado líquido, sobre cuerpos de agua y suelo. Durante la ocurrencia de este escenario en un evento
de inundación, las sustancias tóxicas derramadas suelen propagarse a través de las aguas de
inundación, contaminándolas y por consiguiente contaminando las zonas que entran en contacto
con estas. Por lo tanto, sus efectos están directamente relacionados con la contaminación de los
medios. Algunos de estos medios son, casas, edificios, hospitales, estaciones de bomberos, vías
públicas, aguas subterráneas y ecosistemas[12]. Considerando lo anterior, el volumen de material
derramado, será el parámetro que representará los efectos de contaminación. El cual, en fases
posteriores permitirá estimar las consecuencias asociadas a dichos efectos. A partir del análisis de
eventos Natech generados por inundaciones, se ha identificado que las sustancias que han estado
involucradas en la mayoría de este tipo de eventos son crudos y gasolina, por lo cual se hará foco
en la estimación de los efectos causados por la pérdida de contención de este tipo de materiales[25].

Para estimar la cantidad de material vertido después de que se presenta la falla en el sistema, es
necesario establecer los diferentes modos de liberación de material. Diferentes estudios, han
establecido 3 modos de liberación asociados a la pérdida de contención de material peligroso,
debido a la falla de un tanque, consecuencia del impacto de una inundación[15, 24, 25]. Los modos
de liberación son los siguientes

 Modo de liberación 1: Liberación instantánea del inventario completo.

 Modo de liberación 2: Liberación del inventario completo, después de un tiempo
determinado.
 Modo de liberación 3: Liberaciones menores a través de orificios de diámetros pequeños.

Una vez identificados los modos de liberación, se procede a relacionarlos con los posibles escenarios
accidentales que se pueden presentar si un tanque es impactado por una inundación. En la Figura
3, el último evento del árbol, indica que modo de liberación se presentará para cada una de las
fallas, lo cual permitirá estimar el volumen de material derramado.

Una vez establecida la modalidad de liberación, se hace uso de modelos fuente o modelos de
descarga, los cuales permiten describir el proceso de liberación de material, que se caracteriza a
través de los siguientes conceptos: Tasa de liberación, cantidad total liberada y fase de liberación
(sólido, líquido, gaseoso etc.). Existen un gran número de modelos fuente para estimar la tasa de

26
liberación y tiempo de duración de esta, cada uno de los modelos es útil para situaciones específicas.
Por lo tanto, es necesario establecer ciertos criterios en relación con el tipo de pérdida de
contención establecido en la sección anterior[29, 47].

Para el modo de liberación 2, se debe establecer que este tipo de falla provoca la pérdida de
contención, ya sea por la ruptura total de las conexiones o por la ruptura de la carcasa del tanque.
La libración generad por ambas fallas, se puede modelar como un flujo a través de un orificio en un
tanque en estado dinámico. La diferencia radica en que, para la ruptura de la conexión, se conoce
el diámetro del orificio mientras que para la fractura de la carcasa debe estimarse. En cuanto al
modo de liberación 3, se puede decir que estará relacionado con la ruptura parcial de conexiones,
la cual podrá modelarse de la misma forma que una ruptura total, debido a la flexibilidad del modelo
con respecto al tamaño del orificio. Una vez establecidas las características de la liberación y que
modelo las representa de mejor forma, procederemos a realizar una breve descripción del modelo
de flujo a través de un orificio en un tanque.

Para estimar el flujo a través de un orificio en un tanque, ubicado a una altura ℎ𝐿 por debjo del nivel
del fluido, se debe realizar un modelamiento con base a las ecuaciones de conservación de masa y
la ecuación de Bernoulli[44]. Adicional a esto, el modelo se rige por las siguientes suposiciones.

 Flujo incompresible
 Presión constante en el tiempo
 Trabajo de eje despreciable
 Cambio de velocidad despreciable dentro de la unidad de proceso

Por lo tanto, para el cálculo del caudal de fuga en cualquier tiempo se utiliza la siguiente ecuación
(21)
𝑔𝑐 𝑃𝑔 𝜌𝑔𝐶𝑜2 𝐴2
𝑄𝑚 = 𝜌𝐶𝑜 𝐴√2( + 𝑔ℎ𝐿0 ) − 𝑡1
𝜌 𝐴𝑡

Donde ℎ𝐿 se cálcula como

(22)
𝐶𝑜 𝐴 𝑔𝑐 𝑃𝑔 𝑔 𝐶𝑜 𝐴 2
ℎ𝐿 = ℎ𝐿𝑜 − √2 ( + 𝑔ℎ𝐿0 𝑡) + ( 𝑡 )
𝐴𝑡 𝜌 2 𝐴𝑡 1

Para la estimación del volumen derramado se hace uso de la siguiente expresión

𝑡𝑓 (23)
𝑉𝐿 = ∫ 𝑄𝑚 𝑑𝑡
0

Con las ecuaciones presentadas enteriormente, es posible estimar el volumen de material

derramado. El modelo presentado, es de las mejores aproximación para la estimación del volumen
derramado de material a traves de un orificio, en un tanque de almacenamiento. Una vez, cálculada
la canntidad de sustancia liberada, se procede a estimar los efectos que tendrán las posibles
consecuancias asociadas a dicha pérdida.

27
3.4 Evaluación de las consecuencias: Costos asociados a la pérdida de contención de material
peligroso debido al impacto de una inundación
Las consecuencias de los accidentes tecnológicos, suelen estar asociadas a lesiones o fatalidades de
humanos cerca de la zona del accidente, así como pérdidas monetarias y daños al medio ambiente.
Diferentes relaciones matemáticas, permiten estimar dichas consecuencias en función a los efectos
físicos del escenario accidental [45].

Contrario a lo presentado anteriormente, para eventos Natech provocados por inundaciones, no se

presentan modelos para la estimación de consecuencias. Esto, debido a que las afectaciones de los
accidentes tecnológicos generados por inundaciones, suelen ser impactos sociales y ambientales,
los cuales no son sencillos de representar matemáticamente, ya que darle valor a un individuo a un
organismo es de las tareas más trabajosas en el ámbito del análisis del riesgo.

Antes de presentar la propuesta para la estimación de los efectos de una perdida de contención
durante un evento de inundación, es importante recalcar que, en este tipo de eventos naturales, las
consecuencias que se han identificado debido a la pérdida de contención de sustancias peligrosas
en fase líquida son: daños a propiedad privada debido a la contaminación de casas, edificios,
automóviles, entre otros, daños a la propiedad pública debido a la contaminación de las vías públicas
y la afectación a la fauna y flora, debido al retorno de las aguas de inundación contaminadas a los
canales[12].

Teniendo en cuenta lo anterior y a través de una búsqueda bibliográfica y revisión de bases de datos
de accidentes, se establece una posible medición económica de los efectos, debido a la
contaminación del medio por la pérdida de contención de sustancias peligrosas. La medición de los
impactos consistirá en establecer una relación entre el volumen derramado de material con los
diferentes costos asociados al accidente. Perez et al. [48] proponen un modelo matemático, el cual
relaciona el volumen derramado de material con las consecuencias económicas producto de la
pérdida de contención. Para definir dichas consecuencias, se hizo uso de la base de reportes
accidentales PHMSA, la cual registra información relevante sobre los accidentes, que involucran
pérdida de sustancias peligrosas en fase líquida. Dentro de la información disponible se encuentran
los costos asociados a los derrames. Los diferentes costos que registra la base de datos son: (i) Costo
estimado asociado a daños y reparación de propiedad privada y pública, (ii) costo estimado por la
pérdida de producto, (iii) costo estimado de daños y reparación de la propiedad del operador, (iv)
costo estimado asociado a respuesta a emergencias, (v) costo estimado, asociado a la remediación
del medio ambiente y (vi) otros costos[49].

El modelo para la estimación de costos debido a una pérdida de contención que se utilizará, hace
uso también de características del accidente como:

 Tipo de sustancia involucrada.

 Densidad poblacional en el área del accidente.
 Tipo de incidente.

Para la primera característica se consideran tres tipos de sustancias, refinados del petróleo, crudo y
líquidos altamente volátiles. En cuanto a la densidad poblacional se consideran dos tipos de área
con alta densidad y con baja densidad. Por último, el tipo de incidente recae en dos categorías,

28
incidentes subterráneos y no subterráneos. Lo cual permitirá estimar los costos para diferentes
combinaciones de escenarios.

Para la construcción del modelo, los autores realizaron un análisis sobre cada uno de los sub costos
que reportan y presentan la siguiente expresión para el cálculo del costo en función del volumen
derramado

𝐶𝑖 = 𝛽𝑜𝑖 + 𝛽1𝑖 𝑉 + 𝛽2𝑖 𝑉 2 (24)

Donde i representa cada uno de los sub costos que se quieran integrar a la estimación y la variable
𝛽 son los coeficientes para cada combinación de escenario. La Ecuación 25 presenta la relación para
el cálculo del costo total.
𝑖 (25)
𝐶𝑇 = ∑ 𝐶𝑛
𝑛=1

De esta forma es posible estimar los efectos de un derrame de material peligroso, en función de los
costos asociados al accidente. Como se puede ver, los sub costos que componen el modelo,
consideran los daños que puede generar el derrame a la propiedad privada, la propiedad pública y
el medio ambiente. Lo cual representa los posibles efectos que se pueden presentar durante un
accidente tecnológico provocado por una inundación. Las constantes para el cálculo de los sub
costos se presentan en Anexos, Tabla 11.

Hasta este punto se han presentado los elementos para el cálculo de frecuencias de los posibles
eventos finales y la estimación de las consecuencias de dichos eventos. Por lo tanto, en la siguiente
sección, se presentará la propuesta para el cálculo del riesgo asociado e eventos Natech
consecuencia del impacto de una inundación.

3.5 Estimación del riesgo asociado a eventos Natech generado por inundaciones
El riesgo es un concepto inherente a nuestras vidas, cada día personalmente, tomamos decisiones,
realizamos un sin número de acciones o evitamos hacer otras, las cuales tienen asociado cierto
grado de posibilidad de generar algún daño a nosotros mismos, a quienes están en nuestro
alrededor, nuestra propiedad o la de otros. Esa probabilidad de sufrir algún daño es lo que suele
llamarse riesgo[50]. En la ingeniería, la interpretación más común del riesgo se define como el valor
esperado de las pérdidas. El cual se estima a través de la probabilidad de que se presente un nivel
de pérdida especifico debido a la ocurrencia de un escenario accidental[34]. Este concepto puede
extrapolarse a diferentes áreas como son, el campo financiero y de inversión, la medicina y los
procesos industriales, entre otros; siendo esta última de interés para este trabajo. El Centro para la
Seguridad de Procesos Químicos, CCPS (por sus siglas en inglés, Center for Chemical Process Safety)
definió el concepto de riesgo para accidentes industriales como, la medida de la pérdida económica,
lesiones personales o daño al medio ambiente, en términos de la probabilidad de ocurrencia o
frecuencia del incidente y la magnitud de las pérdidas o lesiones[29]. Es importante resaltar que, el
riesgo no solo es una medida, es un concepto que permite evaluar la ocurrencia de un evento no
deseado para apoyar la toma de decisiones en función de reducir dicha ocurrencia[34].

29
Considerando la definición del riesgo propuesta por la CCPS, el cálculo del riesgo se puede llevar a
cabo combinando, la frecuencia de ocurrencia del evento accidental con las consecuencias de dicho
evento. La Ecuación 26 presenta la expresión algebraica para la estimación del riesgo

𝑅 =𝑓x𝐶 (26)

Haciendo uso del concepto de riesgo, una expresión para el cálculo del riesgo asociado a eventos
Natech provocados por inundaciones, se construyó en función de las frecuencias de los eventos
accidentales que pueden generar este tipo de eventos naturales y el costo asociado a los diferentes
efectos que produce una pérdida de contención producto del impacto de una inundación. A
continuación, se presenta la expresión para el cálculo de este riesgo

𝑅𝑁 = 𝑓𝑒 x 𝐶𝑇 (27)

Donde 𝑓𝑒 es la frecuencia del evento final y 𝐶𝑇 son los costos totales asociados al derrame de
material almacenado. Se puede decir que el riesgo calculado es un riesgo económico, debido al uso
de una variable monetaria.

El proceso de análisis del riesgo para eventos Natech generados por inundaciones, finaliza con la
estimación del riesgo. En los análisis convencionales la representación del riesgo individual y social
se desarrolla mediante herramientas gráficas, con el fin de comunicar de forma efectiva la
información obtenida mediante el análisis. Debido a que este marco de referencia utiliza una
aproximación económica, para hacer una representación del riesgo o poder categorizarlo dentro del
contexto, se presentará una matriz de riesgo, la cual permitirá no solo asignar un nivel a cada uno
de los riesgos sino también evaluar si es necesario tomar medidas para tratarlo.

3.6 Evaluación del riesgo a partir de su calificación

La evaluación de los riesgos estimados a través de los QRA convencionales o las metodologías
alternas, como la propuesta en este trabajo, se desarrolla con el objetivo de alimentar los procesos
de toma de decisiones. Una herramienta muy utilizada son las matrices de riesgo. Las matrices de
riesgo, son diagramas simples que permiten categorizar y priorizar riesgos asociados a diferentes
tipos de eventos, a través de un método gráfico[51]. Las cuales, buscan representar las frecuencias
de ocurrencia del evento junto con los niveles de impacto de las consecuencias. Para establecer los
rangos de frecuencia y consecuencia, suele utilizarse categorías cualitativas, más que valores
cuantitativos, lo cual amplia la incertidumbre al momento de priorizar el riesgo.

30
 Clasificación de consecuencias
C1 C2 C3 C4

F3 R9 R10 R11 R12

Clasificación de
frecuencias F1 R5 R6 R7 R8

F1 R1 R2 R3 R4

Figura [Link] de riesgo de 3x4 categorías y 12 niveles de riesgo [51].

Usualmente, cuando se hace uso de matrices de riesgo para evaluar riesgos, tres niveles son
establecidos y se les han asignado diferentes colores para establecer su diferenciación. Los
diferentes niveles son: eventos con riesgo inaceptable (nivel alto), los cuales hacen uso del color
rojo para diferenciar, riesgos de nivel intermedio, los cuales suelen indicarse con el color amarillo y
cumplen con el criterio ALARP (por sus siglas en inglés As Low As Reasonably Practicable) establecido
para ese tipo de eventos y por último los riesgos de bajo nivel, los cuales suelen ser aceptables y se
les indica con el color verde[51]. Vale la pena aclarar que establecer una región ALARP permite
tomar medidas para reducir los riesgos hasta un nivel razonable dentro de lo factible, esto quiere
decir que existe un nivel durante la reducción donde seguir aplicando las medidas hace que el
proceso sea inoficioso y poco rentable para quien lo aplica. Por lo tanto, debe establecerse un límite
el cual indica el momento de detener las tareas de reducción del riesgo para llegar a un nivel
tolerable y mantener la rentabilidad del proceso[52].

A pesar de que las matrices de riesgo suelen ser clasificadas como herramientas cualitativas, varios
autores han desarrollado estudios donde presentan un análisis, en el cual se establece que las
matrices de riesgo muestran una compatibilidad lógica con los riesgos cuantitativos como los
representados por la Ecuación 27. Por lo tanto, dicha expresión es el principio que gobierna el diseño
del diagrama, debido a esto, la representación del riesgo en la matriz debe replicar lo presentado
en dicha definición matemática. Dicho esto, se establece el uso de matrices de riesgo para evaluar
riesgos cuantitativos[53].

Una vez establecidos los criterios básicos de una matriz de riesgo, se puede proceder a la
construcción de la matriz. Para esto, se presentan tres axiomas que se deben cumplir para la
construcción de una matriz bien definida[53, 54]. Los axiomas se presentan a continuación.

 Poca coherencia: Este axioma consiste en que cada riesgo que este posicionado en la
categoría de color rojo, por inferencia, debe ser un riesgo mayor que los que se encuentran
localizados en la zona verde.
 Intersección de regiones: Siempre se debe cumplir que, si se dibuja una línea recta con
pendiente positiva, que inicia en la zona verde y termina en la zona roja, debe pasar a través
de la categoría de color amarillo.
 Coloración coherente: Los riesgos que tengan un mismo valor o aproximado, deben estar
asignados a una categoría del mismo color.

31
Con base en la teoría presentada anteriormente, a continuación, se presenta la construcción de la
matriz de riesgo para evaluar eventos Natech provocados por inundaciones.

Para la construcción de la matriz de riesgo, se utilizarán los componentes de la definición del riesgo
de eventos Natech generados por inundaciones, presentados anteriormente, para ser asignados
como los ejes de la matriz, de esta forma tenemos que la evaluación se realizará en base a la
frecuencia del evento final y los costos asociados al derrame del material.

3.6.1 Matriz de riesgo para eventos Natech generados por inundaciones

Antes de presentar la matriz de riesgo para eventos Natech provocados por inundaciones, es
importante mencionar que los rangos seleccionados para los ejes pueden distribuirse en su escala
de forma logarítmica. Esto es una ventaja ya que permite realizar una mejor asignación de nivel de
riesgo y la comparación entre estos. Contrario a las matrices escaladas de forma uniforme, las
matrices con escala logarítmica, no permiten que se presenten problemas compresión de rango,
esto quiere decir que riesgos con valores diferentes suelen ser agrupados en una misma categoría,
llevando a desinformar al usuario de la matriz. Adiciona a esto, las matrices escaladas
uniformemente, solo se apoyan en la categorización por zonas de color mientras que las matrices
logarítmicas permiten una categorización del riesgo más amplia, ya que se puede utilizar la
asignación de zonas por colores y la categorización por cada una de las celdas que compone a matriz
(Figura 14)[55].

Figura [Link] de riesgo para eventos Natech provocados por inundaciones.

En la Figura 14, se presenta la matriz de riego para evaluar los riesgos asociados a eventos Natech
generados por inundaciones. La matriz, muestra las diferentes categorías de riesgo que se pueden
evaluar, identificados con las letras en mayúscula. En su trabajo, Levine[55] presenta el uso de letras
para categorizar los riesgos y asegura que el uso de zonas de colores es innecesaria. a pesar de esto,
el autor considera que el uso de colores permite una mejor comunicación de los niveles de riesgo
además de considerar que el uso de colores es una herramienta muy común en el diseño de matrices
de riesgo. Adicional a esto, los analistas ya se encuentran acostumbrados a ello y un cambio podría

32
prestar una mala lectura de la matriz, generando toma de decisiones desviadas. Es por esto que se
decide implementar tanto la categorización por celdas de la matriz y por colorimetría. Por último,
cabe resaltar que la matriz de riesgo es una herramienta versátil ya que permitirá establecer las
regiones de riesgo bajo, alto o ALARP en función de las necesidades del analista. En la Figura 15 se
presenta la forma en que se estableció la región ALARP de la matriz de riesgo, se estableció que los
valores de riesgo que estén entre el rango de 10-1 y 10, serán riesgos tolerables. Dichos rangos de
tolerabilidad se pueden construir con respecto a valores umbral, establecidos por parte de
entidades gubernamentales o industriales.

Figura 15. Matriz de riesgo con los rangos de tolerabilidad indicados.

Con la presentación de la herramienta para la evaluación de los riesgos asociados a eventos Natech
provocados por inundaciones, se finaliza la presentación de la metodología propuesta para el
análisis y evaluación de este tipo de riesgos. Como se puede apreciar, se desarrolló un marco para
el análisis del riesgo de accidentes tecnológicos causado por el impacto de una inundación. La
metodología presentada, integra en cada uno de sus pasos los elementos relacionados con el evento
natural de inundación, permitiendo estudiar los sucesos y eventos que componen un escenario
accidental provocado por una inundación. La metodología permite realizar la estimación de la
probabilidad de daño por la deformación del tanque, para diferentes escenarios de inundación y
diferentes estructuras de almacenamiento, bajo un modelo que integra el tratamiento de la
incertidumbre asociada a los parámetros de entrada. Así mismo, es posible estimar desde un
enfoque económico los efectos de la pérdida de contención de material, así como el cálculo del
riesgo asociado a la dispersión de material peligroso durante un evento de inundación. A
continuación, se desarrollará la aplicación de la metodología propuesta para el análisis del riesgo
asociado a eventos Natech consecuencia del impacto de una inundación, sobre un caso de estudio.

33
4. Caso de estudio
En esta sección se presenta un caso de estudio, sobre el cual se aplicará el marco de referencia
propuesto (Figura 1) para el análisis del riesgo asociado a eventos Natech provocados por
inundaciones.

4.1.1 Descripción del caso de estudio

Con el objetivo de demostrar la aplicación de la metodología propuesta, un análisis del riesgo
asociado a eventos Natech generados por inundaciones se llevará a cabo. En la Figura 16 se presenta
la distribución del parque de almacenamiento, el cual será afectado por una inundación. El estudio
se desarrollará sobre ocho tanques de almacenamiento verticales que operan a presión atmosférica.
El parque de almacenamiento cumple con la función de albergar gasolina producto del proceso de
refinación, la cual está lista para su distribución. En la Tabla 6 se presentan las configuraciones
estructurales de los tanques de almacenamiento presentes en el área del accidente, junto con las
propiedades de la sustancia almacenada.
Tabla 5. Características de los escenarios de inundación considerados en el estudio.

Escenario de Periodo de Frecuencia de Altura de la Velocidad de la

inundación retorno ocurrencia inundación inundación
[años] [1/año] [m] [m/s]
Caso 1 1000 1x10-3 2 0,5
-3
Caso2 200 5x10 0,5 2
Caso 3 30 3,3x10-2 1 1

Figura 16. Distribución de los tanques seleccionados para el estudio. Tomado de [25]

En cuanto al escenario de inundación, se aplicará el marco de referencia sobre tres posibles

escenarios de inundación. El primer caso considera una inundación del tipo inmersión lenta y de alta

34
profundidad, mientras que el segundo caso involucrará una inundación de alta velocidad y baja
profundidad. Por otro lado, el tercer escenario manejará un vector de impacto no muy fuerte,
presentando una velocidad moderada y una profundidad media. En la Tabla 5, se presentan los
valores para velocidad y altura de inundación de cada uno de los casos, junto con sus respectivos
periodos de retorno, sobre los cuales se manejarán 3 frecuencias diferentes, alta frecuencia de
ocurrencia, media y de baja.

Tabla 6. Características de los tanques considerados en el estudio.

Características del tanque T1-T2 T3-T4 T5-T8

Diámetro [m] 4,5 12 24
Altura [m] 9 14,4 14,4
Espesor de la coraza [mm] 5,3 12,5 12,5
Nivel de llenado de operación [%] 50 10 70
Substancia almacenada Gasolina Gasolina Gasolina
Presión [atm] 1 1 1
Densidad en fase líquida [kg/m3] 760 760 760
Inventario [kg] 32175 122100 3418275

En estudios anteriores se han analizado los riesgos asociados a escenarios finales como explosión e
incendios. Por lo tanto, cabe aclarar que en este caso de estudio se analizaran los riesgos asociados
a la perdida de contención de materiales peligrosos (en este caso gasolina), producto de la falla de
tanques de almacenamiento consecuencia del impacto de una inundación. El modo de daño que se
considerará en la implementación del caso de estudio, será el pandeo de la coraza del tanque. El
tipo de falla que se seleccionará para la implementación de la metodología propuesta sobre cada
uno de los casos de inundación, serán los que presenten mayor probabilidad de ocurrencia (ver
Figura 3). Por último, se asume que no existen sistemas de protección para contrarrestar la pérdida
de contención debido a la falla completa de la conexión, por lo tanto, para los escenarios que
consideren esta falla se establece el derrame del inventario total que este almacenado en el tanque.

4.2 Resultados
La metodología propuesta para el análisis del riesgo de eventos Natech, fue aplicada a 3 escenarios
de inundación diferentes que podrían impactar sobre un conjunto de tanques de almacenamiento.
Los resultados del análisis sobre los casos se presentan la Tabla 7. En principio, podemos observar
como las probabilidades de daño, debido al pandeo de las paredes de tanque, estimadas a través
de las funciones Probit, son diferentes para cada conjunto de tanques, en cada uno de los casos. Se
observa que, para el primer caso, el cual considera una velocidad baja y una altura alta de
inundación. Presenta valores de probabilidad cercanos a 1, por lo tanto, si dicha inundación llegase
a ocurrir, los tanques T3-T4 sufrirían una deformación de sus paredes, mientras que los demás
tanques tienen una probabilidad baja de sobrellevar el daño. Para el caso 2, el cual considera una
inundación de baja altura y alta velocidad, se puede observar que la probabilidad de daño calculada

35
es baja para dos conjuntos de tanques, mientras que los tanques T5-T8, no presentarán daño por la
deformación de sus paredes. Esto se debe, a que el componente estático de la carga externa
aplicada sobre el tanque, tiene mayor efecto que el componente dinámico en la deformación por
pandeo de las paredes del tanque. Recordemos que, el componente estático lo proporciona la altura
de la inundación. Observando los resultados para el tercer caso, se evidencia que para los tanques
T1-T2 y T5-T8, presentan valores bajos de probabilidad de daño, contrario a los tanques T3-T4, que
presentan una probabilidad alta. Esto evidencia, la influencia que tiene el nivel de llenado del tanque
en la resistencia del sistema al daño por pandeo, ya que los tanques T3-T4 están operando con el
nivel mínimo de almacenamiento, por lo que el aporte a la resistencia por parte de la presión que
genera el fluido almacenado es muy bajo. Mientras que los otros dos conjuntos de tanques, están
operando con más de la mitad de su inventario total.
Tabla 7. Valores de frecuencia del escenario final, probabilidad de daño debido al pandeo del tanque y costos asociados
al derrame para los escenarios de inundación considerados en el caso de estudio

Cálculo del riesgo T1-T2 T3-T4 T5-T8

Probabilidad de daño 0,9 1 0,84
Falla total de la Falla total de la Falla total de la
Tipo de falla
conexión conexión conexión
Probabilidad de falla 0,4 0,4 0,4
Caso 1 Frecuencia del escenario 4x10-4 4x10-4 3x10-4
final [1/año]
Volumen derramado
360 1366 38236
[bbl]
Costo asociado al
4,6x105 1,08x106 1x107
derrame [$USD]
Probabilidad de daño 0,0017 0,38 0
Falla total de la Colapso de la
Tipo de falla No hay falla
conexión estructura
Probabilidad de falla 0,35 0,15 0
Caso 2 Frecuencia del escenario 3x10-6 3x10-4 0
final [1/año]
Volumen derramado
360 1366 No hay pérdida
[bbl]
Costo asociado al
4,6x105 1,08x106 0
derrame [$USD]
Probabilidad de daño 0,003 0,95 5x10-5
Falla total de la Falla total de la Falla total de la
Tipo de falla
conexión conexión conexión
Probabilidad de falla 0,5 0,5 0,5
Frecuencia del escenario
Caso 3 4x10-5 0,01 8x10-7
final [1/año]
Volumen derramado
360 1366 38236
[bbl]
Costo asociado al
4,6x105 1,08x106 1x107
derrame [$USD]

36
Como se estableció en la descripción del caso de estudio, se asume que el volumen de material
derramado es la totalidad del inventario que almacena el tanque durante su operación. De esta
forma, se estiman los costos asociados al derrame para cada uno de los tanques, los cuales se
presentan en la Tabla 7. Como se puede ver, los tanques del conjunto T5-T8 son los de mayor
capacidad y los que operan con un nivel de llenado mayor, por lo que la pérdida total de su
inventario, generará unos costos mayores comparados con los demás conjuntos de tanques.

Una vez estimada la frecuencia del escenario final accidental y las consecuencias asociadas a este,
se procede a realizar la evaluación de los riesgos asociados a los tanques, presentados en el caso de
estudio. Antes de entrar a discutir sobre los riesgos estimados, es importante presentar la función
que tiene el riesgo dentro de este marco de referencia. Debido a la falta de modelos para estimar
los efectos de una dispersión de sustancias toxicas en fase líquida, se estableció un riesgo económico
el cual no busca representar pérdidas monetarias, sino establecer la magnitud de la afectación en
función del dinero necesario para remediar las afectaciones. Por lo tanto, el riesgo de eventos
Natech generados por inundaciones, indicará, la cantidad de dinero necesaria para remediar las
consecuencias si dicho evento llegase a ocurrir, en un periodo de tiempo establecido.

La Figura 17, presenta la evaluación de los riesgos calculados para el caso 1. De la figura se puede,
observar como el riesgo asociado a los tanques T1-T2 cae dentro de la zona de riesgo intermedio, el
cual es un riesgo tolerable. A pesar de esto, debido a que, en el caso de estudio no se presentan
barreras para mitigar o prevenir el suceso, se puede establecer que existen medidas que, de ser
tomadas disminuirían la categoría del riesgo asociado a estos tanques. En cuanto a los tanques T3-
T4 y T5-T8, son categorizados como riesgos de alto nivel, debido al alto impacto que tendría el
evento si llegase a suceder.

Figura [Link]ón de los riesgos calculados en el caso 1.

37
La Figura 18, presenta la categorización de los riesgos estimados para el segundo caso de estudio.
Ya que el impacto de la inundación sobre los tanques T5-T8 no genera ningún daño por pandeo, no
se presentará pérdida de contención de material peligroso. Por lo tanto, no se desarrolló la
estimación del riesgo. Volviendo a los otros dos conjuntos de tanques, podemos ver como el primer
grupo se mantuvo dentro de la región de riesgo intermedio, al igual que los tanques T3-T4, los cuales
se mantienen en un nivel de riesgo alto. Vale la pena comparar las matrices de riesgo para los casos
1 y 2, ya que la diferencia entre estos es que sus vectores de impacto tienen las magnitudes
invertidas. Como se mencionó anteriormente las altas velocidades de inundación no son las
principales causantes del daño por pandeo a los tanques y esto se ve reflejado, al ver como una
disminución en la altura de la inundación, disminuye la probabilidad de daño del equipo de forma
considerable, obteniendo para el caso 1 una categoría de riesgo F y para el caso 2 una categoría D.
Lo cual nos indica que el riesgo asociado a los tanques T1-T2 es mayor para el caso 1 que para el
caso 2. Además, esta comparación deja ver la consistencia de las matrices propuestas para evaluar
los riesgos generados por eventos Natech, ya que representa de forma concreta lo que sucede con
el riesgo cuando hay cambios en el evento iniciador.

Figura [Link]ón de los riesgos calculados en el caso 2.

La evaluación de los riesgos calculados para el último de los casos de estudio, se presenta en la
Figura 19. Para este caso, se evidencia que el riesgo para los taques T3-T4 ha sido categorizado en
la letra I, la cual asigna la categoría más alta de riesgo, debido a sus altas consecuencias y alta
posibilidad de ocurrencia. Esto debido a que la probabilidad de sufrir el pandeo de sus paredes es
alta y al combinarla con la frecuencia de ocurrencia de la inundación, se obtiene un escenario final
accidental con altas chances de ocurrir. En cuanto a los otros dos conjuntos de tanques, podemos
ver como su categoría se mantiene, ya que los tanques T1-T2 presentan una mayor ocurrencia que
los tanques T5-T8, pero en cuanto a sus consecuencias ocurre lo contrario.

38
 10-1
E F G H I
T3-T4
10-2
D E F G H
Frecuencia [año-1]

10-3
C D E F G

10-4
B C D E F
T1-T2
10-5
A B C D E
T5-T8
10-6
102 103 104 105 106 107
Costo del derrame [$USD]

Figura [Link]ón de los riesgos calculados en el caso 3.

5. Conclusiones

Para el análisis y evaluación de los riesgos asociados a eventos Natech provocados por inundaciones,
se propone una metodología sistemática y repetible, que integra información cualitativa y
cuantitativa relacionada con las causas y consecuencias de accidentes industriales, causados por el
impacto de un frente de inundación. El foco de este trabajo es, el de analizar los eventos indeseados,
que se pueden presentar si un tanque de almacenamiento vertical, llegase a fallar al ser afectado
por una inundación. Para la identificación de los posibles escenarios finales accidentales, se
construyeron arboles de eventos, los cuales permiten establecer a partir de un escenario de
inundación la secuencia de eventos que pueden ocurrir, para que suceda una pérdida de contención
del material almacenado en el tanque. Se encontró que pueden existir dos modalidades de daño y
que la materialización de estas puede generar cinco tipos de falla, los cuales se traducirán en
derrames de material peligroso.

Dentro de la estimación de frecuencias de los escenarios finales accidentales, se desarrolla una

herramienta para el cálculo de la probabilidad de daño debido al pandeo de las paredes del tanque,
consecuencia del impacto de un frente de inundación. Esta herramienta, son funciones Probit las
cuales, se construyeron con base a datos sintéticos obtenidos a partir del tratamiento de la
incertidumbre asociada a los parámetros de entrada del modelo, para la estimación de daño por la
deformación del tanque. El tratamiento se realiza utilizando la simulación de Monte Carlo, la cual
permite realizar el cálculo de la probabilidad de daño para diferentes configuraciones estructurales,
en función de un vector de impacto de la inundación y así obtener resultados que representen de
forma más acertada la realidad. Adicional al cálculo de la probabilidad, una gran cantidad de
información es valorizada a través de la aplicación del método Probit, generando un conjunto de

39
ecuaciones, las cuales permiten estimar la probabilidad de daño debido al pandeo de la estructura,
para cualquier geometría de tanque vertical afectada por una inundación con una velocidad y altura
determinados.

Uno de los fuertes de este trabajo, es el enfoque en la estimación de las consecuencias de la

dispersión de crudo y gasolina sobre las aguas de inundación. La estimación de las consecuencias de
este escenario accidental, se realizó a partir de una aproximación económica, donde se seleccionó
un modelo para el cálculo de los costos asociados a un derrame de las sustancias mencionadas. La
decisión de representar las consecuencias a través de un concepto monetario, se tomó con el fin de
poder establecer una magnitud para las consecuencias, más no para establecer una pérdida
económica.

Se presenta un concepto tradicional del riesgo para eventos Natech provocados por inundaciones.
El cual considera la combinación de la frecuencia del escenario accidental final con las consecuencias
asociadas a la materialización de dicho escenario. En cuanto a la fase de evaluación, se elabora una
matriz de riesgo en escala logarítmica, para la categorización y comparación de los riesgos
considerados en este trabajo.

Por último, se aplica la metodología propuesta a un caso de estudio, donde se consideran diferentes
escenarios de inundación, las cuales impactan, diferentes configuraciones estructurales de los
tanques. Se evidencia, que la aplicación de la metodología permite estimar el riesgo asociado a
diferentes escenarios accidentales. Así como determinar, que el nivel de llenado del tanque es un
parámetro que afecta el resultado análisis, ya que genera un aporte importante a la resistencia de
la estructura al ser impactada por la inundación. Adicional a esto, mediante el uso de la herramienta
de evaluación del riesgo, se encuentra que las categorizaciones son consecuentes con el escenario
accidental que se está analizando.

A pesar de proponer diferentes métodos para el desarrollo de las fases de la metodología propuesta,
el trabajo futuro en el área del análisis de riesgos, asociados a eventos Natech es grande,
considerando las afectaciones que puede generar la ocurrencia de estos. Este trabajo se centra en
el desarrollo de herramientas para la estimación de la probabilidad de ocurrencia de uno de los
modos de daño, se plantea aplicar el algoritmo para la construcción de funciones Probit, enfocado
al modo de daño de desplazamiento del tanque. En cuanto a la estimación de consecuencias, es
necesario invertir mayores esfuerzos en poder estimar las afectaciones, de forma que permitan
establecer riesgos individuales y sociales, los cuales se comunican de una mejor manera y facilitan
la toma de decisiones con base en riesgo.

40
6. Referencias
[1] J. T. Hardy, Climate change: causes, effects, and solutions. John Wiley & Sons, 2003.
[2] J. E. Montealegre, "Estudio de la variabilidad climática de la precipitación en Colombia
asociada a procesos oceánicos y atmosféricos de meso y gran escala," Disponible desde
Internet en: [Link] ideam. gov. co/jsp/812 (con acceso 01/11/2014).[Links],
2009.
[3] S. Banholzer, J. Kossin, and S. Donner, "The Impact of Climate Change on Natural
Disasters," in Reducing Disaster: Early Warning Systems For Climate Change, A. Singh and
Z. Zommers, Eds. Dordrecht: Springer Netherlands, 2014, pp. 21-49.
[4] M. UNISDR, "UNISDR Terminology for Disaster Risk Redution," United Nations
International Strategy for Disaster Reduction (UNISDR) Geneva, Switzerland, 2009.
[5] A. M. Cruz and N. Okada, "Methodology for preliminary assessment of Natech risk in
urban areas," Natural Hazards, vol. 46, no. 2, pp. 199-220, 2008// 2008.
[6] A. M. Cruz and N. Okada, "Consideration of natural hazards in the design and risk
management of industrial facilities," Natural Hazards, journal article vol. 44, no. 2, pp.
213-227, 2008.
[7] M. Campedel, "Analysis of major industrial accidents triggered by natural events reported
in the principal available chemical accident databases," Report EUR, vol. 23391, 2008.
[8] S. Young, L. Balluz, and J. Malilay, "Natural and technologic hazardous material releases
during and after natural disasters: a review," Science of The Total Environment, vol. 322,
no. 1–3, pp. 3-20, 4/25/ 2004.
[9] E. Krausmann, A. Cruz, and E. Salzano, Natech risk assessment and management : reducing
the risk of natural-hazard impact on hazardous installations. Boston, MA: Elsevier, 2016, p.
pages cm.
[10] A. M. Cruz and E. Krausmann, "Hazardous-materials releases from offshore oil and gas
facilities and emergency response following Hurricanes Katrina and Rita," Journal of Loss
Prevention in the Process Industries, vol. 22, no. 1, pp. 59-65, 1// 2009.
[11] E. Krausmann, A. M. Cruz, and B. Affeltranger, "The impact of the 12 May 2008 Wenchuan
earthquake on industrial facilities," Journal of Loss Prevention in the Process Industries, vol.
23, no. 2, pp. 242-248, 3// 2010.
[12] D. Davis, "The Aftermath of Hurricanes Katrina and Rita on South Louisiana.," L. A. O. S. R.
a. D. Program, Ed., ed. Baton Rouge, LA: United States Enviromental Protection Agency.
[13] E. CRED, "EM-SAT: The OFDA/CRED International Disaster Database, Université Catholique
de Louvain, Brussels (Belgium)," ed, 2010.
[14] K. Rasmussen, "Natural events and accidents with hazardous materials," Journal of
Hazardous Materials, vol. 40, no. 1, pp. 43-54, 1// 1995.
[15] V. Cozzani, M. Campedel, E. Renni, and E. Krausmann, "Industrial accidents triggered by
flood events: Analysis of past accidents," Journal of Hazardous Materials, vol. 175, no. 1–3,
pp. 501-509, 3/15/ 2010.
[16] E. Renni, E. Krausmann, and V. Cozzani, "Industrial accidents triggered by lightning,"
Journal of Hazardous Materials, vol. 184, no. 1–3, pp. 42-48, 12/15/ 2010.
[17] E. Krausmann, E. Renni, M. Campedel, and V. Cozzani, "Industrial accidents triggered by
earthquakes, floods and lightning: lessons learned from a database analysis," Natural
Hazards, vol. 59, no. 1, pp. 285-300, 2011// 2011.
[18] S. Girgin and E. Krausmann, "Historical analysis of U.S. onshore hazardous liquid pipeline
accidents triggered by natural hazards," Journal of Loss Prevention in the Process
Industries, vol. 40, pp. 578-590, 3// 2016.

41
[19] E. Salzano, I. Iervolino, and G. Fabbrocino, "Seismic risk of atmospheric storage tanks in
the framework of quantitative risk analysis," Journal of Loss Prevention in the Process
Industries, vol. 16, no. 5, pp. 403-409, 9// 2003.
[20] G. Landucci, G. Antonioni, A. Tugnoli, and V. Cozzani, "Release of hazardous substances in
flood events: Damage model for atmospheric storage tanks," Reliability Engineering &
System Safety, vol. 106, pp. 200-216, 10// 2012.
[21] G. Landucci, A. Necci, G. Antonioni, A. Tugnoli, and V. Cozzani, "Release of hazardous
substances in flood events: Damage model for horizontal cylindrical vessels," Reliability
Engineering & System Safety, vol. 132, pp. 125-145, 12// 2014.
[22] C. El Hajj, E. Piatyszek, A. Tardy, and V. Laforest, "Development of generic bow-tie
diagrams of accidental scenarios triggered by flooding of industrial facilities (Natech),"
Journal of Loss Prevention in the Process Industries, vol. 36, pp. 72-83, 7// 2015.
[23] G. Antonioni, G. Spadoni, and V. Cozzani, "A methodology for the quantitative risk
assessment of major accidents triggered by seismic events," Journal of Hazardous
Materials, vol. 147, no. 1–2, pp. 48-59, 8/17/ 2007.
[24] V. Cozzani, G. Antonioni, G. Landucci, A. Tugnoli, S. Bonvicini, and G. Spadoni,
"Quantitative assessment of domino and NaTech scenarios in complex industrial areas,"
Journal of Loss Prevention in the Process Industries, vol. 28, pp. 10-22, 4/11/ 2014.
[25] G. Antonioni, G. Landucci, A. Necci, D. Gheorghiu, and V. Cozzani, "Quantitative
assessment of risk due to NaTech scenarios caused by floods," Reliability Engineering &
System Safety, vol. 142, pp. 334-345, 10// 2015.
[26] A. Necci, G. Antonioni, S. Bonvicini, and V. Cozzani, "Quantitative assessment of risk due to
major accidents triggered by lightning," Reliability Engineering & System Safety, vol. 154,
pp. 60-72, 10// 2016.
[27] O. Salvi et al., "ARAMIS: Accidental Risk Assessment Methodology for Industries in the
framework of SEVESO II directive," in Colloque" Urbanisation et Prévention des Accidents
Majeurs", 2002.
[28] O. Fabio, "Marco para el tratamiento de incertidumbre en el análisis de riesgo cuantitativo
en transporte de material peligroso a través de tuberías.," 2016.
[29] C. f. C. P. Safety, Guidelines for chemical process quantitative risk analysis. Center for
Chemical Process Safety/AIChE, 2000.
[30] C. Castaño-Uribe, R. Carrillo, and F. Salazar, "Sistema de Información Ambiental de
Colombia Tomo III. Perfil del estado de los recursos naturales y del medio ambiente en
Colombia 2001," IDEAM. Min. Medio Ambiente. Bogotá. 595p.[Links], 2002.
[31] B. Mundial, "Análisis de la gestión del riesgo de desastres en Colombia: un aporte para la
construcción de políticas públicas," 2012.
[32] M. T. Anees et al., "Numerical modeling techniques for flood analysis," Journal of African
Earth Sciences, vol. 124, pp. 478-486, 12// 2016.
[33] I. Portugués-Mollá, X. Bonache-Felici, J. F. Mateu-Bellés, and J. B. Marco-Segura, "A GIS-
Based Model for the analysis of an urban flash flood and its hydro-geomorphic response.
The Valencia event of 1957," Journal of Hydrology, vol. 541, Part A, pp. 582-596, 10//
2016.
[34] S. Mauricio, "Introducción a la confiabilidad y evaluación de riesgos," Facultad de
Ingeniería, Universidad de los Andes, Ediciones Uniandes. Bogotá, Colombia, 2004.
[35] G. W. Parry, "The characterization of uncertainty in probabilistic risk assessments of
complex systems," Reliability Engineering & System Safety, vol. 54, no. 2, pp. 119-126,
1996.

42
[36] J. C. Helton, J. D. Johnson, and W. L. Oberkampf, "An exploration of alternative approaches
to the representation of uncertainty in model predictions," Reliability Engineering &
System Safety, vol. 85, no. 1, pp. 39-71, 2004.
[37] T. Aven, E. Zio, P. Baraldi, and R. Flage, Uncertainty in risk assessment: the representation
and treatment of uncertainties by probabilistic and non-probabilistic methods. John Wiley
& Sons, 2013.
[38] C. Baudrit, D. Dubois, and D. Guyonnet, "Joint propagation and exploitation of probabilistic
and possibilistic information in risk assessment," IEEE transactions on fuzzy systems, vol.
14, no. 5, pp. 593-608, 2006.
[39] S. P. Timoshenko and J. M. Gere, Theory of elastic stability. Courier Corporation, 2009.
[40] J. Tilton, "Perry’s Chemical Engineers’ Handbook Section 6: Fluid and Particle Dynamics,"
McGraw-Hill Professional, vol. 4, pp. 12-13, 2007.
[41] O. T. Gudmestad and G. Moe, "Hydrodynamic coefficients for calculation of hydrodynamic
loads on offshore truss structures," Marine Structures, vol. 9, no. 8, pp. 745-758, 1996.
[42] A. S. o. M. Engineers, Boiler and Pressure Vessel code, Section VIII, Div. 1. New York:
American Society of Mechanical Engineers, 1989.
[43] B. Poling, G. Thomson, D. Friend, R. Rowley, and W. Wilding, "Physical and chemical data,"
Perry’s chemical engineers’ handbook, vol. 8, pp. 2-446, 2008.
[44] D. A. Crowl and J. F. Louvar, Chemical process safety: fundamentals with applications.
Pearson Education, 2001.
[45] T. G. Book, "Methods for the determination of possible damage," Rijswijk, The
Netherlands: TNO, 1989.
[46] G. Antonioni, S. Bonvicini, G. Spadoni, and V. Cozzani, "Development of a framework for
the risk assessment of Na-Tech accidental events," Reliability Engineering & System Safety,
vol. 94, no. 9, pp. 1442-1450, 9// 2009.
[47] P. U. De Haag and B. Ale, "Guidelines for quantitative risk assessment (purple book),"
Committee for the Prevention of Disasters, The Hague (NL), 1999.
[48] S. C. Andrés Perez, Felipe Muñoz, "Model for Optimal Sectioning of Liquid Hidrocabon
Pipeline Transport," 2016.
[49] P. a. H. M. S. Administration. Accident Reort Hazardous Liquid Pipeline Systems [Online].
Available: [Link]
transmission-and-gathering-lng-and-liquid-accident-and-incident-data
[50] American Institute of Chemical Engineers. Center for Chemical Process Safety., Guidelines
for developing quantitative safety risk criteria. New York
Hoboken, NJ: Center for Chemical Process Safety ;
Wiley, 2009, pp. xxxv, 211 p.
[51] N. J. Duijm, "Recommendations on the use and design of risk matrices," Safety Science,
vol. 76, pp. 21-31, 7// 2015.
[52] "Chapter 9 - Hazard Assessment A2 - Mannan, Sam," in Lees' Loss Prevention in the Process
Industries (Fourth Edition)Oxford: Butterworth-Heinemann, 2012, pp. 284-404.
[53] L. Anthony Cox, "What's Wrong with Risk Matrices?," Risk Analysis, vol. 28, no. 2, pp. 497--
512, 2008.
[54] R. Flage and W. Røed, "A reflection on some practices in the use of risk matrices," in 11th
International Probabilistic Safety Assessment and Management Conference and the Annual
European Safety and Reliability Conference, 2012, pp. 881-891.
[55] E. S. Levine, "Improving risk matrices: the advantages of logarithmically scaled axes,"
Journal of Risk Research, vol. 15, no. 2, pp. 209-222, 2012/02/01 2012.

43
Nomenclatura
Letras latinas

𝐴 Área del orificio (m2)

𝐴𝑡 Área transversal del tanque (m2)
𝐶𝑖 Sub costo asociado al derrame ($USD)
𝐶𝑜 Coeficiente de descarga (-)
𝐶𝑇 Costo total del derrame ($USD)
𝐷 Diámetro del tanque (m)
𝐸 Módulo de elasticidad (Pa)
𝑓 Frecuencia de ocurrencia de una inundación (año-1)
𝑓𝑒 Frecuencia del escenario final accidental (año-1)
𝑔 Gravedad de la tierra (m/s2)
𝐻 Altura del tanque de almacenamiento (m)
ℎ𝐿 Altura del orificio (m)
ℎ𝐿𝑜 Nivel inicial del fluido sobre el orificio de la fuga (m)
ℎ𝑤 Altura de la inundación (m)
𝑘𝑤 Coeficiente hidrodinámico del agua (-)
𝑛 Parámetro para minimizar la presión crítica (-)
𝑃𝑐𝑟 Presión crítica del tanque (Pa)
𝑃𝑓 Presión generada por el fluido almacenado (Pa)
𝑃𝑔 Presión manométrica (Pa)
𝑃𝑤 Carga externa aplicada por un frente de inundación (Pa)
𝑃𝑤𝑑 Presión dinámica máxima generada por la inundación (Pa)
𝑃𝑤𝑠 Presión estática máxima generada por la inundación (Pa)
𝑃𝑟 Presión de resistencia del tanque (Pa)
𝑝𝑑 Probabilidad de daño debido al pandeo del tanque (-)
𝑄𝑚 Tasa de liberación (kg/s)
𝑅𝑁 Riesgo asociado a eventos Natech provocados por inundaciones($USD/año)
𝑡 Espesor e la coraza del tanque (mm)
𝑡1 Tiempo de descarga (s)
𝑡𝑓 Tiempo final d la fuga (s)
𝑡𝑟 Periodo de retorno de una inundación (año)
𝑉𝐿 Volumen derramado (m3)
𝑣𝑤 Velocidad de la inundación (m/s)
𝑌 Puntos Probit (-)
𝑥 Factor geométrico (-)

Letras griegas

𝛽 Coeficiente para la estimación de los costos asociados al derrame (-)

𝜈 Coeficiente de Poisson (-)
𝜌 Densidad del fluido derramado (kg/m3)
𝜌𝑤 Densidad del agua (kg/m3)
𝜙 Nivel de llenado del tanque (-)

44
Anexos
Tabla 8. Configuración estructural de los tanques utilizados para el cálculo de la probabilidad de daño debido al pandeo
del tanque y la generación de datos sintéticos

ID Diámetro [m] Altura [m] Espesor [mm]

1 3 10,8 5
2 4,5 5,4 5
3 4,5 7,2 5
4 4,5 9 5,3
5 4,5 10,8 5,8
6 6 3,6 5
7 6 5,4 5
8 6 7,2 5,5
9 6 9 6,1
10 6 10,8 6,8
11 6 12,6 7,4
12 6 14,4 8
13 7,5 3,6 5
14 7,5 5,4 5,3
15 7,5 7,2 6,1
16 7,5 9 6,9
17 7,5 10,8 7,7
18 7,5 12,6 8,5
19 7,5 14,4 9,3
20 7,5 16,2 10,1
21 7,5 18 10,9
22 9 3,6 5
23 9 5,4 5,7
24 9 7,2 6,7
25 9 9 7,7
26 9 10,8 8,6
27 9 12,6 9,6
28 9 14,4 10,6
29 9 16,2 11,5
30 9 18 12,5
31 10,5 3,6 5,1
32 10,5 5,4 6,2
33 10,5 7,2 7,3
34 10,5 9 8,4
35 10,5 10,8 9,6
36 10,5 12,6 10,7
37 10,5 14,4 11,8

45
38 10,5 16,2 12,5
39 10,5 18 12,5
40 12 3,6 5,4
41 12 5,4 6,6
42 12 7,2 7,9
43 12 9 9,2
44 12 10,8 10,5
45 12 12,6 11,8
46 12 14,4 12,5
47 12 16,2 12,5
48 12 18 12,5
49 13,5 3,6 5,7
50 13,5 5,4 7,1
51 13,5 7,2 8,6
52 13,5 9 10
53 13,5 10,8 11,5
54 13,5 12,6 12,5
55 13,5 14,4 12,5
56 13,5 16,2 12,5
57 13,5 18 12,5
58 15 3,6 6
59 15 5,4 7,6
60 15 7,2 9,2
61 15 9 10,8
62 15 10,8 12,4
63 15 12,6 12,5
64 15 14,4 12,5
65 15 16,2 12,5
66 15 18 12,5
67 18 3,6 6,5
68 18 5,4 8,5
69 18 7,2 10,4
70 18 9 12,3
71 18 10,8 12,5
72 18 12,6 12,5
73 18 14,4 12,5
74 18 16,2 12,5
75 18 18 12,5
76 21 3,6 7,1
77 21 5,4 9,4
78 21 7,2 11,6

46
79 21 9 12,5
80 21 10,8 12,5
81 21 12,6 12,5
82 21 14,4 12,5
83 21 16,2 12,5
84 24 3,6 7,7
85 24 5,4 10,3
86 24 7,2 12,5
87 24 9 12,5
88 24 10,8 12,5
89 24 12,6 12,5
90 24 14,4 12,5
91 26 7,2 12,5
92 27 3,6 8,3
93 27 5,4 11,2
94 27 7,2 12,5
95 27 9 12,5
96 27 10,8 12,5
97 27 12,6 12,5
98 30 5,4 12,1
99 30 7,2 12,5
100 30 9 12,5
101 30 10,8 12,5
102 36 3,6 10,1
103 36 5,4 12,5
104 36 9 12,5
105 42 3,6 11,3
106 42 5,4 12,5
107 42 7,2 12,5
108 48 3,6 12,4
109 48 5,4 12,5
110 48 7,2 12,5
111 54 3,6 12,5
112 54 5,4 12,5
113 60 3,6 12,5
114 60 5,4 12,5
115 66 3,6 12,5
116 66 5,4 15

47
Tabla 9. Constantes Probit para el cálculo de la probabilidad de daño por pandeo de un tanque debido al impacto de una inundación, para sustancias almacenadas con valor de
densidad líquida de 760 kg/m3.

Nivel de Velocidad de la K1 K2 K3
llenado inundación a b c a b c a b c
0.5-1 6,7917 -1,0309 0,0232 1,1378 -0,0772 0,0021 10,2749 0,0457 -0,0013
1-1.5 6,6784 -0,9661 0,0212 3,2038 -0,2194 0,0060 9,0785 0,1355 -0,0038
1.5-2 4,9345 -0,7313 0,0140 6,0002 -0,3866 0,0102 7,1003 0,2876 -0,0080
0.1 2-2.5 2,3500 -0,5002 0,0082 8,3384 -0,4078 0,0092 4,9499 0,3899 -0,0100
2.5-3 0,4793 -0,5296 0,0124 9,1029 -0,1486 -0,0012 3,3718 0,3424 -0,0067
3-3.5 -1,1643 -0,6099 0,0168 9,1361 0,1323 -0,0102 2,5937 0,1824 -0,0005
3.5-4 -1.8052 -0.7326 -0.0198 8.6711 0.3246 0.014 2.2588 0.0445 0.0030
0.5-1 3,9059 -1,1916 0,0269 0,8757 -0,0673 0,0020 10,4782 0,0090 0,0010
1-1.5 4,0753 -1,1601 0,0261 2,4768 -0,1907 0,0054 9,5662 0,0865 -0,0015
1.5-2 3,2642 -1,0292 0,0224 4,7487 -0,3560 0,0100 8,1271 0,2125 -0,0055
0.2 2-2.5 1,2418 -0,7758 0,0151 7,4211 -0,5301 0,0146 6,1046 0,3902 -0,0110
2.5-3 -1,3605 -0,5155 0,0082 9,6107 -0,5754 0,0146 4,0632 0,5092 -0,0138
3-3.5 -2,9521 -0,5199 0,0111 10,3228 -0,3619 0,0058 2,6215 0,4856 -0,0116
3.5-4 -4,3794 -0,5831 0,0155 10,5050 -0,1085 -0,0032 1,8280 0,3550 -0,0060
0.5-1 1,8586 -1,3219 0,0253 0,7388 -0,0706 0,0026 10,6345 -0,0689 0,0101
1-1.5 2,2091 -1,3659 0,0295 2,1006 -0,1929 0,0068 9,9044 0,0393 0,0042
1.5-2 1,4157 -1,1146 0,0183 3,9550 -0,3352 0,0109 8,7777 0,1248 0,0012
0.3 2-2.5 0,1041 -0,9280 0,0143 6,2849 -0,5149 0,0162 7,1635 0,2743 -0,0047
2.5-3 -2,0539 -0,6744 0,0083 8,6977 -0,6628 0,0200 5,1325 0,4626 -0,0116
3-3.5 -4,4144 -0,4536 0,0036 10,5420 -0,6817 0,0190 3,3105 0,5706 -0,0148
3.5-4 -5,5501 -0,5098 0,0091 10,9781 -0,4559 0,0095 2,0422 0,5530 -0,0134
0.5-1 1,2512 -1,9678 0,0569 0,6290 -0,0648 0,0029 10,2098 0,1438 0,0040
0.4
1-1.5 1,4125 -1,8833 0,0532 1,7877 -0,1835 0,0081 9,6717 0,1797 0,0011

48
 1.5-2 1,0318 -1,7183 0,0456 3,4970 -0,3625 0,0159 8,7448 0,2550 -0,0025
2-2.5 -0,0766 -1,4682 0,0347 5,6443 -0,5754 0,0247 7,3518 0,3801 -0,0080
2.5-3 -1,9967 -1,0935 0,0199 7,9666 -0,7858 0,0327 5,6966 0,5340 -0,0151
3-3.5 -4,5941 -0,6541 0,0033 9,8932 -0,8256 0,0303 3,9055 0,6435 -0,0184
3.5-4 -6,7363 -0,4271 -0,0017 11,1967 -0,7349 0,0232 2,6106 0,6230 -0,0158
0.5-1 0,4326 -2,5276 0,1003 0,5794 -0,0677 0,0036 10,2101 0,2040 0,0058
1-1.5 0,3032 -2,2447 0,0760 1,5580 -0,1633 0,0083 9,6315 0,2905 -0,0023
1.5-2 0,2920 -2,2142 0,0784 3,0457 -0,3147 0,0154 8,7360 0,4114 -0,0120
0.5 2-2.5 -0,6104 -1,9561 0,0672 4,9806 -0,5323 0,0259 7,6715 0,4827 -0,0165
2.5-3 -1,9576 -1,6654 0,0496 7,1267 -0,7404 0,0361 6,1250 0,6351 -0,0218
3-3.5 -4,1821 -1,2042 0,0275 9,1560 -0,8571 0,0382 4,4619 0,7488 -0,0261
3.5-4 -6,7726 -0,7508 0,0063 11,0642 -0,9520 0,0406 2,8231 0,8464 -0,0294
0.5-1 -0,2144 -3,1870 0,1722 0,5243 -0,0826 0,0069 10,1416 0,3407 -0,0015
1-1.5 0,3767 -3,2958 0,1838 1,4948 -0,1952 0,0143 9,4450 0,4921 -0,0154
1.5-2 0,0966 -3,0908 0,1740 2,9635 -0,4403 0,0334 8,8306 0,5322 -0,0222
0.6 2-2.5 -1,4367 -2,1900 0,0769 4,6079 -0,5841 0,0381 7,9411 0,4983 -0,0127
2.5-3 -2,6242 -1,9233 0,0616 6,3090 -0,6127 0,0317 6,6333 0,6347 -0,0214
3-3.5 -4,3687 -1,5339 0,0429 8,3772 -0,8152 0,0412 5,0617 0,7722 -0,0291
3.5-4 -6,3125 -1,2284 0,0282 10,2411 -0,9271 0,0460 3,3106 0,9351 -0,0368
0.5-1 -0,4584 -4,1776 0,3147 0,4774 -0,0717 0,0070 9,8012 0,6968 -0,0362
1-1.5 -0,6999 -3,5814 0,2151 1,3695 -0,2136 0,0211 9,5012 0,6073 -0,0220
1.5-2 -1,0399 -3,2489 0,1920 2,6860 -0,4159 0,0378 8,9759 0,5986 -0,0305
0.7 2-2.5 -1,6741 -2,8236 0,1299 4,2660 -0,6075 0,0523 8,0298 0,6202 -0,0172
2.5-3 -2,4476 -2,7210 0,1409 6,1796 -0,8505 0,0671 6,6974 0,8483 -0,0446
3-3.5 -4,6054 -1,8786 0,0661 7,9326 -0,9048 0,0586 5,3982 0,8910 -0,0438
3.5-4 -6,2269 -1,5963 0,0467 9,3307 -0,7887 0,0417 4,0291 0,8969 -0,0364

49
Tabla 10. Constantes Probit para el cálculo de la probabilidad de daño de un tanque debido al impacto de una inundación, para sustancias almacenadas con valor de densidad
líquida de 1000 kg/m3

Nivel de Velocidad de K1 K2 K3
Llenado inundación a b c a b c a b c
0.5-1 5,838513 -1,11305 0,025807 1,040983 -0,07257 0,00193 10,34158 0,044337 -0,00122
1-1.5 5,811082 -1,05264 0,023822 2,937562 -0,21041 0,005795 9,225261 0,131275 -0,00363
1.5-2 4,447625 -0,85335 0,017737 5,551433 -0,38175 0,010324 7,448689 0,273111 -0,00773
0.1 2-2.5 1,866797 -0,56783 0,009582 8,232122 -0,49713 0,012714 5,291865 0,412733 -0,01129
2.5-3 -0,25047 -0,49305 0,009813 9,460581 -0,33189 0,005331 3,493001 0,432 -0,01032
3-3.5 -2,44896 -0,46757 0,011462 10,08199 -0,11017 -0,00278 2,516978 0,298433 -0,00447
3.5-4 -3,07196 -0,64081 0,017837 9,61297 0,157325 -0,01032 2,073009 0,131317 0,001133
0.5-1 2,330081 -1,20201 0,022921 0,76858 -0,0601 0,001768 10,53159 -0,03257 0,005416
1-1.5 2,561654 -1,17297 0,022235 2,185205 -0,18356 0,005733 9,829809 0,025955 0,003245
1.5-2 2,052234 -1,08079 0,020586 4,213046 -0,3457 0,010569 8,571853 0,146695 -0,00137
0.2 2-2.5 0,532707 -0,88216 0,01536 6,678239 -0,5242 0,015594 6,818981 0,312325 -0,00703
2.5-3 -1,8762 -0,60902 0,008539 9,133331 -0,65582 0,018622 4,739208 0,490977 -0,01308
3-3.5 -4,02084 -0,4611 0,006412 10,59688 -0,59093 0,014703 3,015709 0,560908 -0,01465
3.5-4 -4,58669 -0,6498 0,016029 10,51948 -0,26589 0,002093 1,879045 0,507785 -0,01166
0.5-1 1,557428 -2,0751 0,068184 0,646657 -0,06888 0,00319 10,22758 0,133766 0,003163
1-1.5 1,37285 -1,81244 0,04983 1,810752 -0,18722 0,008219 9,741373 0,133577 0,00382
1.5-2 1,018285 -1,6817 0,043025 3,548129 -0,36982 0,016188 8,810267 0,233209 -0,00105
0.3 2-2.5 -0,08738 -1,4164 0,031346 5,687818 -0,58355 0,025085 7,44085 0,35213 -0,00629
2.5-3 -2,26633 -0,97245 0,0119 7,933837 -0,72935 0,028387 5,830244 0,475159 -0,01095
3-3.5 -4,7904 -0,59996 0,000498 10,11444 -0,85391 0,031332 3,960442 0,624691 -0,01709
3.5-4 -6,81632 -0,41444 -0,00174 11,24805 -0,72791 0,02271 2,595161 0,62129 -0,01573
0.5-1 0,211178 -2,62968 0,100416 0,553536 -0,06911 0,004499 10,13542 0,261359 0,00579
0.4
1-1.5 0,441011 -2,59516 0,105497 1,588833 -0,20206 0,012726 9,700897 0,306466 -0,00199

50
 1.5-2 0,011751 -2,30458 0,083768 2,937309 -0,29064 0,014282 8,771498 0,447733 -0,01365
2-2.5 -0,84693 -2,00641 0,066913 4,792897 -0,49385 0,024409 7,770265 0,468443 -0,01297
2.5-3 -2,11932 -1,76201 0,056072 6,902924 -0,70519 0,034718 6,298128 0,639245 -0,02369
3-3.5 -4,10056 -1,35918 0,036827 8,841457 -0,79336 0,034528 4,624041 0,765326 -0,02861
3.5-4 -6,67784 -0,84854 0,010687 10,81169 -0,94351 0,041475 2,984648 0,853928 -0,03056
0.5-1 -0,89804 -3,21915 0,157076 0,487423 -0,06744 0,005978 9,961712 0,482602 -0,00347
1-1.5 -0,67199 -3,167 0,167336 1,38772 -0,18961 0,015457 9,609285 0,492974 -0,01401
1.5-2 -0,521 -3,1864 0,1719 2,737829 -0,38837 0,031692 8,850795 0,587819 -0,02164
0.5 2-2.5 -0,49971 -2,00206 0,067361 4,959088 -0,51678 0,02541 7,653487 0,4962 -0,01648
2.5-3 -2,40155 -2,49158 0,129355 6,52246 -0,93307 0,069347 6,530148 0,868173 -0,05122
3-3.5 -4,28012 -1,83894 0,056363 7,81236 -0,70374 0,037836 5,316455 0,799285 -0,03004
3.5-4 -6,42704 -1,37555 0,03994 9,766175 -0,87329 0,043574 3,757452 0,907801 -0,03737
0.5-1 -2,24669 -3,5642 0,220872 0,433014 -0,06173 0,006876 10,11888 0,48764 0,002021
1-1.5 -1,06867 -4,29817 0,31 1,210967 -0,15856 0,017363 9,205545 0,984133 -0,05593
1.5-2 0,05467 -3,06091 0,170378 2,963053 -0,44325 0,0339 8,862526 0,510965 -0,021
0.6 2-2.5 -1,33521 -2,27385 0,087761 4,583891 -0,5606 0,034948 7,882988 0,546201 -0,01883
2.5-3 -3,30786 -2,58022 0,076847 5,614073 -0,74283 0,07113 7,049998 0,74356 -0,01083
3-3.5 -4,53514 -2,35874 0,098604 7,570427 -0,98983 0,083599 5,637215 0,97496 -0,0499
3.5-4 -6,14017 -2,03016 0,093097 9,504219 -1,22011 0,095815 3,986016 1,222358 -0,07681
0.5-1 -3,75512 -3,99594 0,383794 0,399915 -0,05811 0,008202 10,08041 0,906948 -0,09095
1-1.5 -2,95906 -4,26894 0,388848 1,116929 -0,1534 0,020596 9,561575 0,986072 -0,07622
1.5-2 -2,14442 -4,37245 0,260413 2,145925 -0,30058 0,049308 8,760226 1,014948 -0,01272
0.7 2-2.5 -2,29342 -4,25891 0,294341 3,485966 -0,47537 0,067418 7,93563 1,116309 -0,0601
2.5-3 -2,25917 -2,81989 0,1531 6,145192 -0,84403 0,066725 6,542661 0,932208 -0,05465
3-3.5 -4,58352 -3,31029 0,240097 6,970201 -0,93357 0,094554 5,98984 1,162491 -0,08942
3.5-4 -6,69944 -2,19022 0,078712 8,80758 -1,20118 0,120274 4,647619 1,210237 -0,08052

51
 Tabla 11. Coeficientes para la estimación de costos asociados al derrame de material peligroso

Sub costos
Costo asociado a la
Costo por la pérdida Costo asociado a Otros
Coeficientes remediación del
de producto respuesta a emergencias costos
medio ambiente
𝛽𝑜 78,14 284,3 1,099 18,75
𝛽1 -0,0206 -0,112 -0,358 -0,0122
𝛽2 - 11,837 - 1,127

52