En este artículo se ha realizado una evaluación probabilista e integral del riesgo sísmico para las edificaciones de Lorca, la ciudad más afectada tras el sismo de mayo de 2011 en la Región de Murcia, España. La amenaza sísmica se ha representado mediante un conjunto de escenarios estocásticos que permite considerar el efecto de eventos pequeños, moderados y extremos en las posibles pérdidas futuras. Mediante funciones de transferencia espectrales, se ha considerado la respuesta dinámica del suelo. Se ha utilizado una base de datos de elementos expuestos con nivel de resolución edificio por edificio, la cual permite la desagregación de los resultados en diferentes categorías y, además, la obtención de mapas de riesgo sísmico para la visualización de la distribución geográfica de las pérdidas futuras. Para cada una de las tipologías estructurales identificadas se ha asignado una función de vulnerabilidad que permite obtener las pérdidas esperadas para diferentes niveles de aceleración. Los resultados de riesgo se han obtenido en forma de curva de excedencia de pérdidas a partir de la cual pueden derivarse otras métricas probabilistas de riesgo como la pérdida anual esperada y la pérdida máxima probable. Los resultados de riesgo resultan útiles para los tomadores de decisiones en ámbitos de planeación de emergencia, desarrollo de planes de refuerzo de edificaciones existentes y protección financiera mediante mecanismos tradicionales de seguro y reaseguro u otros instrumentos alternativos para la transferencia del riesgo.
A comprehensive and probabilistic seismic risk assessment has been performed for the buildings of Lorca, the most affected city in the region of Murcia, Spain, after the May 2011 earthquake. Seismic hazard is represented through a set of stochastic scenarios that allow accounting for small, moderate and extreme events in the future losses, also through spectral transfer functions the dynamic soil response has been considered. A building by building resolution level database has been used allowing the disaggregation of risk results in several categories besides allowing the generation of risk maps to visualize the geographical distribution of the future losses. For each of the identified building classes a vulnerability function has been assigned to determine the expected losses for different acceleration levels. Risk results have been obtained in terms of the loss exceedance curve from where other probabilistic risk metrics such as the average annual loss and the probable maximum loss can be derived. Risk results are useful for the decision-makers in the fields of emergency planning, existent building retrofitting schemes and financial protection through traditional insurance and reinsurance schemes or with alternative risk transfer instruments.
El 11 de mayo de 2011, un sismo de magnitud 5,1 (MW) sacudió la Región de Murcia y, en específico, la ciudad de Lorca fue la más afectada. El epicentro se localizó a poco más de 5km al norte del casco urbano de Lorca, lo que, combinado con la poca profundidad del evento y a pesar de la magnitud moderada del sismo, causó la muerte de 9 personas, heridas a más de 300, mientras que alrededor de 10.000 no pudieron regresar a sus hogares. Históricamente, varios eventos han sacudido la región, tales como aquellos ocurridos en marzo de 1911, enero de 1917, junio de 1948, agosto de 2002 y enero de 2005 [1], pero solo el de mayo de 2011 causó fallecidos, heridos, daños estructurales severos e interrupciones en el funcionamiento normal de la comunidad. Adicionalmente, debido a la inexperiencia en la aplicación de un plan de emergencias, los daños en la infraestructura desembocaron en una situación caótica, con el agravante que muchas de las acciones de respuesta tomaron más tiempo del esperado inicialmente por la comunidad [2].
En trabajos anteriores se han realizado comparaciones de los daños registrados y modelizados para escenarios sísmicos con características similares a las del sismo de Lorca [3,4], las cuales han permitido la validación de las metodologías e información base utilizada en dichos trabajos. Algunos estudios se han realizado incluso antes del sismo del 2011 [5]. En este artículo se cuantifica el riesgo sísmico no en forma de proporción medias de daño (mean damage ratio [MDR]), sino a partir de la curva de excedencia de pérdidas (CEP) obtenida tras la convolución de la amenaza y la vulnerabilidad de los elementos expuestos. A partir de esta se pueden derivar métricas probabilistas como la pérdida anual esperada (PAE) y la pérdida máxima probable (PMP), que permiten la caracterización con mayor detalle de los elementos expuestos y de la amenaza sísmica que las evaluaciones regionales previas [6].
Se ha ensamblado una base de datos de elementos expuestos con un nivel de resolución edificio por edificio, el más detallado para este tipo de evaluaciones, la cual captura información relevante en cuanto a parámetros y tipologías estructurales, los cuales han sido asignados a partir de índices basados en el más reciente censo de población y vivienda [7]. Es importante tener en cuenta que, cuando se realizan análisis probabilistas de riesgo sísmico, una de las suposiciones está relacionada con la ley de los grandes números, es decir, que se requiere un número considerable de elementos expuestos, con lo que aun cuando se produzcan errores por sobre o subestimación, el resultado final tiende al valor medio.
Para la obtención de los resultados de riesgo sísmico en términos de la CEP se requiere una representación probabilista de la amenaza sísmica a partir de escenarios estocásticos, así como de la definición de funciones de vulnerabilidad para las diferentes tipologías estructurales identificadas en el área bajo análisis. La evaluación probabilista del riesgo sísmico se ha realizado únicamente para las edificaciones, tanto públicas como privadas, localizadas dentro del área urbana de Lorca y considerando pérdidas económicas directas; sin embargo, los resultados presentados pueden servir como insumo para evaluaciones que abarcan otras dimensiones del riesgo, como las propuestas en las referencias [8–13].
Dentro de las diferentes herramientas para la evaluación del riesgo sísmico, se ha seleccionado la Plataforma CAPRA1[14–19], la cual consiste en diferentes módulos para la evaluación de la amenaza sísmica, la vulnerabilidad estructural y el riesgo físico. La plataforma realiza procesos de evaluación y simulación numérica por ordenador que abarca varios ámbitos, tal y como se presenta en la figura 1
2Evaluación probabilista de la amenaza sísmica
La evaluación de la amenaza sísmica para Lorca se ha realizado con un enfoque probabilista para dar cuenta de las incertidumbres asociadas al proceso de generación de terremotos tales como su localización, profundidad, magnitud, frecuencia de ocurrencia y movimiento del terreno [20]. Para ello, la amenaza en roca se expresa a partir de un conjunto de escenarios estocásticos y, además, debido a que existe información relacionada con la respuesta dinámica del suelo y se han identificado una serie de zonas de suelo homogéneo [21], se han definido funciones de transferencia espectrales para la obtención del movimiento del terreno en el suelo y no en el basamento rocoso. A continuación se presentan los detalles de la evaluación de la amenaza sísmica.
2.1Evaluación de la amenaza en rocaComo se ha mencionado anteriormente, la metodología de evaluación probabilista de riesgo aplicada en este estudio requiere un conjunto de escenarios estocásticos que representen de manera exhaustiva e integral la amenaza sísmica en Lorca. Para ello, se ha realizado una evaluación tanto probabilista como espectral de la amenaza sísmica que considera fuentes sismogénicas tanto en la península Ibérica como en zonas aledañas, donde cada una de dichas fuentes se ha caracterizado a partir de un modelo de sismicidad Gutenberg-Richter [22]. La zonificación tectónica es la propuesta en el marco del proyecto SHARE [23], más una fuente adicional en el norte de África, tal como se muestra en la figura 2.
Como no es posible predecir la ocurrencia de terremotos en el tiempo y la duración de una ventana de observación completa es desconocida, ha sido generado un conjunto con 50.982 escenarios estocásticos. Para caracterizar el proceso de ocurrencia de los sismos en cada una de las fuentes sismogénicas se ha asumido un modelo de Poisson donde la actividad sísmica está representada a partir de curvas de recurrencia de magnitudes, las cuales relacionan la frecuencia de ocurrencia de sismos con determinada magnitud en cada fuente. Consecuentemente, la sismicidad se define a partir de la siguiente ecuación:
que es una versión modificada del modelo Gutenberg-Richter [24], donde λ(M) es la tasa de recurrencia para eventos con magnitud igual o superior al valor umbral seleccionado M0 (3,5 en este caso); λ0, β y MU son los parámetros de sismicidad calculados para cada fuente sismogénica y M es la variable aleatoria que representa la magnitud de cada fuente. λ0 describe la tasa anual de ocurrencia con magnitud igual o superior a la umbral; β representa la pendiente inicial de la regresión logarítmica, lo que es igual a la relación entre sismos grandes y pequeños, y MU es la magnitud última asociada a cada fuente sismogénica. Se ha utilizado el catálogo sísmico publicado por el Instituto Geológico Nacional [25], donde se ha realizado una verificación de completitud para la magnitud umbral seleccionada, siguiendo la metodología propuesta por Tinti y Mulargia [26]. Los parámetros λ0 y β se calcularon utilizando un procedimiento de máxima verosimilitud [27], mientras que los valores esperados de MU se asignaron a partir de los publicados en el marco del proyecto SHARE [23].
Dentro de los diferentes modelos geométricos posibles, se ha seleccionado el de área para representar cada fuente sismogénica como un plano. Cada fuente principal se divide en subfuentes con la misma sismicidad por unidad de área así como con la misma probabilidad espacial de ocurrencia de sismos. Las subfuentes son polígonos triangulares donde, desde su centroide, se generan eventos con diferentes magnitudes y frecuencias, compatibles con los parámetros de sismicidad.
Se han utilizado relaciones de atenuación de movimiento fuerte calibradas para el entorno europeo [28], con lo que la amenaza sísmica se puede calcular en términos de curva de excedencia de aceleraciones (para diferentes ordenadas espectrales) mediante la siguiente ecuación:
donde ν(a) es la tasa de excedencia de la aceleración, N es el número de fuentes sismogénicas, Pr(A>a|M,Ri) es la probabilidad de que la aceleración sea excedida dada la magnitud y la distancia de un sismo entre la n-ésima fuente y el punto de interés, Ri. Asumiendo que, dada una magnitud y distancia, la aceleración sigue una distribución log-normal, la probabilidad Pr(A>a|M, Ri) puede calcularse como:donde Φ(·) es la distribución normal estándar, MED(A|M, Ri) es la mediana de la aceleración, dada por la relación de atenuación de movimiento fuerte asociada para un par magnitud-distancia conocido, y σLna es la desviación estándar del logaritmo natural de la aceleración.
Las curvas de excedencia de aceleraciones se calcularon para diferentes ordenadas espectrales, por lo que es posible obtener espectros de amenaza uniforme (EAU) para diferentes períodos de retorno. La evaluación de la amenaza sísmica se ha realizado utilizando el programa CRISIS 2014 V1.2 [29], que corresponde al módulo de amenaza sísmica de la Plataforma CAPRA. La figura 3 muestra las curvas de excedencia de aceleración para diferentes ordenadas espectrales, mientras que la figura 4 presenta los EAU para 225, 475 y 1.000 años de período de retorno. Estos resultados para Lorca tienen en cuenta la contribución de todas las fuentes sismogénicas incluidas en el modelo que estén dentro de la distancia de integración, es decir, 300km.
Las características dinámicas del suelo pueden modificar las ondas sísmicas tanto en su contenido frecuencial como en amplitud, por lo que, siempre que exista información para la determinación de zonas de suelo homogéneo (microzonificaciones sísmicas), es deseable incluirlo tanto en la evaluación de la amenaza como en la del riesgo sísmico.
Con base en el estudio de Navarro et al. [21] sobre la microzonificación sísmica del casco urbano de Lorca, para cada una de las zonas de suelo homogéneo se ha asignado un tipo de suelo de acuerdo al Eurocódigo-8. A partir de dicha información es posible calcular funciones de transferencia espectrales para cada una de las zonas de suelo homogéneo calculando el cociente entre la forma funcional para el tipo de suelo identificado y la forma funcional para roca. Es importante notar que la metodología empleada para la definición de las zonas homogéneas de suelo está basada en deformaciones pequeñas que no capturan completamente el comportamiento no lineal de este.
3Inventario de elementos expuestosLas evaluaciones de riesgo sísmico requieren la elaboración de bases de datos que contienen los elementos expuestos, en este caso las edificaciones del casco urbano de Lorca. Dichas bases de datos se pueden construir utilizando diferentes niveles de resolución, dentro de los cuales el de edificio por edificio es el más detallado. La selección del nivel de resolución depende no solamente de la información disponible sino del alcance de la evaluación; en el presente caso, dado que se ha realizado una evaluación de la amenaza sísmica detallada, incluyendo la microzonificación sísmica y, además, que la información catastral permite un alto nivel de detalle, se ha seleccionado esa escala. El proceso de elaboración de las bases de datos de elementos expuestos siempre ha significado un reto para la modelación dado que, en muchos casos, la información requerida no se encuentra disponible en una única fuente y que, en muchos casos, necesita ser generada o inferida a partir de índices obtenidos de múltiples fuentes. Para este estudio se requiere información relacionada con la ubicación geográfica y las características estructurales tales como edad, material, sistema estructural, número de plantas y tipología estructural y, por esto, dichos parámetros se han asignado a cada uno de los elementos expuestos mediante el procedimiento que se describe a continuación.
3.1Caracterización a partir de la información catastral disponibleSe dispone de la información catastral actualizada para Lorca [30] con un nivel de resolución de elemento por elemento. Tras la remoción de elementos diferentes a edificaciones, tales como terrazas, plazas y balcones, un total de 17.064 elementos quedaron incluidos en la base de datos inicial. La información catastral tiene datos relacionados con la ubicación geográfica y el número de plantas de cada elemento. La figura 5 presenta la distribución geográfica de alturas (en número de plantas) para cada elemento del casco histórico de Lorca. Puede observarse que la mayoría de edificaciones son de baja altura, entre 1 y 3 plantas.
3.2Información disponible para la asignación de la vulnerabilidad
Con base en el censo poblacional y de vivienda más reciente de España [7], fue posible definir la distribución de edades para las edificaciones en Lorca. Dichos parámetros fueron asignados a los elementos usando los datos presentados en la tabla 1.
Además, con base en estudios anteriores [5], y usando la información de edad presentada anteriormente, se realizó la clasificación de la vulnerabilidad basada en la escala EMS-98 [31] con los datos presentados en la tabla 2. De allí es posible observar que las estructuras se han clasificado en categorías entre A y D en dicha escala.
Clases de vulnerabilidad en la escala EMS-98 de acuerdo con la edad para las edificaciones de Lorca
| Clase de vulnerabilidad EMS-98 | A | B | C | D |
|---|---|---|---|---|
| Edad | ||||
| Antes de 1900 | 80% | 20% | – | – |
| 1900-1920 | 72% | 28% | – | – |
| 1921-1940 | 72% | 28% | – | – |
| 1941-1950 | 69% | 28% | 3% | – |
| 1951-1960 | 46% | 49% | 5% | – |
| 1961-1970 | 18% | 38% | 44% | – |
| 1971-1980 | 5% | 40% | 55% | – |
| 1981-1990 | – | 38% | 57% | 5% |
| 1991-2001 | – | 28% | 62% | 10% |
| 2002-2011 | – | 18% | 69% | 13% |
No se cuenta con el valor catastral de las edificaciones y, por dicha razón, se ha desarrollado un índice basado en el área total construida para inferir el valor de reposición de cada elemento. Dicho valor tiene el objetivo de capturar el valor de reparación o reposición de cada elemento para llevarlo a exactamente la misma condición con la que ha sido inicialmente caracterizado. El principal objetivo de esta valoración es el de establecer un orden de magnitud del valor de reposición de las edificaciones de Lorca como conjunto. En el presente trabajo, los valores de reposición no tienen en cuenta los valores históricos o patrimoniales de las estructuras.
Con base en INE [7] se ha definido un valor base de 1.247 euros por metro cuadrado construido pero, para considerar que no todos los elementos pueden tener la misma valoración, el parámetro de edad ha sido seleccionado como diferencial. Debido a que la reparación de mamposterías en piedra y bloques es más costosa que la de edificaciones en forjados de hormigón armado debido a la necesidad de mano de obra más especializada, ha sido definido un factor que aumenta con la edad, tal y como se presenta en la tabla 3. El valor total de reposición de los edificios de Lorca se ha establecido en 6.927 millones de euros.
Valores de reposición y factores por edad para Lorca
| Edad | Factor edad | Costo por m2 |
|---|---|---|
| Antes de 1900 | 2,00 | 2.494 € |
| 1900-1920 | 2,00 | 2.494 € |
| 1921-1940 | 1,75 | 2.182 € |
| 1941-1950 | 1,75 | 2.182 € |
| 1951-1960 | 1,50 | 1.871 € |
| 1961-1970 | 1,50 | 1.871 € |
| 1971-1980 | 1,50 | 1.871 € |
| 1981-1990 | 1,25 | 1.559 € |
| 1991-2001 | 1,25 | 1.559 € |
| 2002-2011 | 1,00 | 1.247 € |
Habiendo identificado la distribución de edades y de clases de vulnerabilidad, fueron identificadas una serie de tipologías estructurales a partir del trabajo de la referencia [7]. Fue asignada una clase de vulnerabilidad en la escala EMS-98 a cada tipología estructural. Las edificaciones en Lorca están construidas principalmente de diferentes tipos de mampostería (ladrillos y piedra) para las edificaciones de baja altura, mientras que para los edificios de más niveles son comunes los forjados en hormigón armado. En las zonas de uso industrial se encuentran pórticos en acero y estructuras prefabricadas en hormigón armado.
Combinando los parámetros para todos los elementos, ha sido asignada una única tipología estructural a cada edificación; en total se utilizaron en el análisis 10 tipologías estructurales (véase la tabla 4). En la segunda columna de la tabla se presenta el código asignado a cada una de ellas, mientras que en la tercera columna se muestra la clase de vulnerabilidad asociada. A partir de una visita de campo se dividió la ciudad en 11 zonas, para las cuales, en cada una de ellas, se realizó la inspección de validación de las tipologías estructurales inicialmente identificadas y su distribución relativa.
Tipologías estructurales, códigos de abreviación y clase de vulnerabilidad en la escala EMS-98
| Tipología estructural | Código | Clase de vulnerabilidad EMS-98 |
|---|---|---|
| Mampostería de piedra | M-PP | A |
| Tapial | M-TA | A |
| Encintado toledano | M-ET | B |
| Mampostería de ladrillo | M-L | B |
| Forjados de hormigón | M-H | C |
| Hormigón armado (pre 1995) | E-H | C |
| Hormigón armado (post 1995) | E-H2 | D |
| Hormigón con diagonales en acero | E-HX | D |
| Hormigón prefabricado | E-HF | C |
| Estructuras metálicas | E-MT | D |
La pérdida L se define como una variable aleatoria, por lo que la variación de los momentos estadísticos de ella para diferentes valores de demanda sísmica se describen mediante funciones de vulnerabilidad. La distribución de probabilidad de la pérdida pL|S(L) se asume como una función beta (ecuación 4) donde los momentos estadísticos corresponden a la media (conocida como la relación media de daño) y a la desviación estándar:
donde Γ es la función gamma y sus parámetros a y b son:E(L|S) es el valor esperado de la pérdida y c(L|S) es el coeficiente de variación de la pérdida dada una demanda sísmica S. Las funciones de vulnerabilidad, al contrario que las escalas cualitativas de daños, contienen toda la información necesaria para calcular la probabilidad de que se iguale o exceda un nivel de pérdida dada una demanda sísmica a partir de la siguiente ecuación:donde l es una pérdida dentro del dominio de definición de la variable aleatoria L, y S es de nuevo la demanda sísmica. El daño es cuantificado a partir de la relación media de daño, que se calcula como la relación entre el costo de reparación estimado y el valor total expuesto de cada elemento. Consecuentemente, la función de vulnerabilidad se define relacionando el daño medio y la intensidad sísmica expresada como aceleración máxima del terreno (para edificaciones de baja altura) o como aceleración espectral (para edificaciones de mayor altura). Dado un nivel de aceleración sísmica para cada tipología estructural, la relación medida de daño puede calcularse mediante la metodología propuesta por Miranda [32] y Ordaz [33]donde β es la pérdida, γ0 y γi son parámetros de vulnerabilidad estructural que dependen de la tipología estructural y de la edad de construcción, ¿ es la pendiente y E(·) es el valor esperado. Por definición, β es la relación entre el costo de reparación y el costo total de la edificación, la cual, dado que únicamente se están evaluando las pérdidas físicas directas, toma valores entre 0 y 1.
Adicional al valor esperado del daño, se calcula la dispersión de este para diferentes niveles de aceleración. Esta dispersión es igual a cero para valores extremos y toma su valor máximo cuando el daño medio es igual al 50%.
Se requiere una función de vulnerabilidad para cada tipología estructural identificada y, además, la diferencia en el comportamiento sísmico de las edificaciones puede considerarse a partir del período fundamental de vibración de cada una de ellas, es decir, cada función de vulnerabilidad tiene asociada una ordenada espectral que corresponde al período fundamental elástico típico de cada tipología cuyo daño está siendo caracterizado.
Se utilizaron en la simulación realizada un total de 22 funciones de vulnerabilidad, las cuales han sido desarrolladas en el marco del modelo global de riesgo [34,35] incluido en el Global Assessment Report on Disaster Risk Reduction de la Estrategia Internacional para la Reducción de Riesgo de Desastre de las Naciones Unidas [36]. La figura 6 muestra diferentes funciones de vulnerabilidad, de las cuales resulta claro que algunas tipologías estructurales, especialmente aquellas de mampostería simple, son mucho más vulnerables que las otras, teniendo, para un mismo nivel de aceleración, un nivel de daño esperado mayor. La nomenclatura de la tabla 4 ha clasificado el número de pisos en 3 categorías: bajo (L) para edificios entre 1 y 3 plantas, mediana (M) para edificios entre 4 y 7 plantas y de gran altura (H) para 8 y más niveles.
5Evaluación probabilista del riesgo sísmico5.1Metodología
Para la evaluación probabilista del riesgo sísmico se ha utilizado el procedimiento analítico propuesto por Ordaz [33], que es el que se encuentra implementado en la Plataforma CAPRA [14]. Así pues, se realiza el cálculo de la función de densidad de probabilidad para la pérdida l en el elemento expuesto j, condicionada a la ocurrencia del escenario i, mediante la siguiente ecuación:
Resulta evidente que esta distribución de probabilidad no puede calcularse de manera directa, por lo cual se realiza un encadenamiento de 2 distribuciones de probabilidad separadas, tal como puede verse en la siguiente ecuación:
Con esto, la función de densidad de probabilidad de la pérdida asociada a cada escenario es calculada mediante la agregación de las pérdidas de cada elemento expuesto. Dado que la pérdida se asume como una variable aleatoria, se requiere calcular el valor esperado de la pérdida y su varianza asociada para cada escenario para su correcta consideración:
NE es el número total de elementos expuestos, E(lj) es el valor esperado de la pérdida en el elemento expuesto j dada la ocurrencia del escenario i, σ2(lj) es la varianza de la pérdida en el elemento expuesto j dada la ocurrencia del escenario i, y cov(lk,lj) es la covarianza de la pérdida de 2 elementos expuestos separados. La covarianza se calcula asumiendo un coeficiente de correlación (σk,j) igual a 0,3 y considerando la desviación estándar para las pérdidas en los diferentes elementos expuestos:Una de las mejores maneras de comunicar los resultados de riesgo es a partir de la curva de excedencia de pérdida (CEP), que asocia las frecuencias con las que los eventos igualan o exceden un valor predeterminado de pérdida. Esta frecuencia, generalmente expresada en términos anuales, es conocida como la tasa de excedencia y se calcula mediante la siguiente ecuación, que es una de las formas del teorema de probabilidad total:donde v(l) es la tasa de excedencia de la pérdida l, N es el número de escenarios en el conjunto estocástico, Pr(L>l|Eventoi) es la probabilidad de exceder l dada la ocurrencia del evento i, mientras que FA(Eventoi) es la frecuencia anual de ocurrencia del evento i; la suma de la ecuación considera todos los eventos que generan cualquier nivel de pérdida. El valor inverso de v(l) es el período de retorno de la pérdida l.5.2Resultados
La tabla 5 presenta el resumen de los resultados de evaluación probabilista de riesgo sísmico en Lorca. Es posible observar que la pérdida anual esperada, PAE, es de aproximadamente 16,3 millones de euros, lo que, en términos relativos al valor total expuesto, representa el 2,36‰; dicho valor corresponde a un nivel alto de riesgo sísmico, el cual se encuentra principalmente dominado por la alta vulnerabilidad estructural de las edificaciones. Adicionalmente, se presentan los valores de pérdida máxima probable, PMP, para diferentes períodos de retorno donde, por ejemplo, para 500 años, se estima una pérdida de 690 millones de euros que corresponde al 10% del valor total expuesto. La estimación de las pérdidas aseguradas para el sismo de mayo de 2011 es del orden de 490 millones de euros [37], las cuales incluyen condiciones específicas de los contratos como deducibles, coaseguro y límites asegurables que no permiten una comparación directa con los valores calculados para las pérdidas físicas directas para la totalidad de edificaciones evaluadas.
Los resultados pueden desagregarse para cualquiera de los factores incluidos en la base de datos de elementos expuestos y se presentan en este caso de acuerdo a las tipologías estructurales para ver cuáles son las que concentran los mayores niveles de riesgo.
La tabla 6 muestra la desagregación de los resultados de riesgo sísmico por tipologías estructurales; de ella es evidente que las de construcciones en tierra, así como de mampostería, son las que mayores valores de riesgo tienen, esto en términos tanto absolutos como relativos. Lo anterior se encuentra asociado con la vulnerabilidad sísmica de dichas tipologías que tienen un pobre comportamiento ante la demanda sísmica debido a las prácticas constructivas deficientes, así como a la utilización de materiales no adecuados. Los edificios con forjados reticulares en hormigón también tienen una contribución importante en las pérdidas modeladas debido también a su alta vulnerabilidad física [38–42].
Desagregación de los resultados de riesgo por tipología estructural
| Tipología estructural | Número de elementos | Valor expuesto (millones de €) | % del valor expuesto | Pérdida anual esperada (millones de €) | Pérdida anual esperada (‰) | % de la pérdida |
|---|---|---|---|---|---|---|
| Mampostería de piedra | 1.838 | 848 | 12,2% | 3,26 | 3,84 | 19,94 |
| Tapial | 1.955 | 978 | 14,1% | 4,31 | 4,41 | 26,42% |
| Encintado toledano | 528 | 203 | 2,9% | 0,54 | 2,68 | 3,33% |
| Mampostería de ladrillo | 5.207 | 2.057 | 29,7% | 4,24 | 2,06 | 25,97% |
| Forjados de hormigón | 2.963 | 1.156 | 16,7% | 2,74 | 2,37 | 16,77% |
| Hormigón armado (pre 1995) | 3.432 | 1.293 | 18,7% | 0,66 | 0,51 | 4,04% |
| Hormigón armado (post 1995) | 485 | 161 | 2,3% | 0,12 | 0,75 | 0,74% |
| Hormigón con diagonales en acero | 35 | 8 | 0,1% | 0,01 | 0,89 | 0,04% |
| Hormigón prefabricado | 593 | 216 | 3,1% | 0,44 | 2,01 | 2,66% |
| Estructuras metálicas | 28 | 8 | 0,1% | 0,01 | 1,38 | 0,07% |
| Total | 17.064 | 6.927.7 | 100 | 16,3 | 2,36 | 100 |
Dado que la evaluación de riesgo sísmico se ha realizado en una base de datos georreferenciada, se han elaborado mapas de riesgo para observar la distribución geográfica de las pérdidas esperadas. Para una mejor comparación de los niveles de riesgo se ha seleccionado la PAE relativa, dado que constituye una métrica robusta y, en términos generales, insensible a las incertidumbres [43]. La figura 7 muestra dicha distribución.
6Conclusiones
Los modelos de riesgo sísmico a nivel urbano dan estimaciones que pueden resultar útiles para los tomadores de decisiones en términos de daños esperados y de recursos requeridos. La identificación del riesgo es la primera etapa en un esquema integral de gestión de riesgo [44] y sus resultados pueden ser utilizados en el diseño de planes de emergencia, regulaciones de planeación urbana, definición de planes para el reforzamiento de edificaciones existentes y estrategias de protección financiera [45,46].
La cuantificación del orden de magnitud de las pérdidas antes de la ocurrencia del desastre permite la evaluación de acciones ex-ante tales como bonos de catástrofe, créditos contingentes, fondos de reserva para desastre e incluso mecanismos tradicionales de seguro y reaseguro, en vez de las medidas ex-post tales como nuevos impuestos, reasignación presupuestal o ayuda internacional para la atención y recuperación de la zona afectada.
Se han obtenido estimaciones de las pérdidas directas para el portafolio de edificaciones en Lorca, las cuales pueden agruparse en diferentes categorías como edad, número de plantas, clase estructural y tipología estructural, entre otras. Lo anterior resulta útil para la determinación de dónde se encuentra concentrado el riesgo. Los resultados permiten la generación de mapas de riesgo con un nivel de resolución elemento por elemento mediante los cuales es posible visualizar la distribución geográfica de las pérdidas futuras; sin embargo, debe tenerse en cuenta que, en cualquier caso, el riesgo debe ser preferiblemente expresado en términos de tasas y probabilidades de excedencia de valores de pérdida. Estos resultados pueden actualizarse de manera periódica en la medida que mejor información exista relacionada con las características de las edificaciones y sus variaciones a lo largo del tiempo, así como de información de microzonificación sísmica más detallada.
Resultados de investigaciones futuras en los campos de la geotecnia, mediante la determinación de funciones de transferencia para diferentes niveles de excitación sísmica en el basamento rocoso y la ingeniería estructural mediante el desarrollo de funciones de vulnerabilidad más específicos para las edificaciones en Lorca, principalmente aquellas de forjados reticulares en hormigón armado, pueden utilizarse para la actualización de los resultados obtenidos.
Por último, estos resultados pueden usarse como base para la evaluación holística de riesgo sísmico en zonas urbanas. Este enfoque, además de las pérdidas directas, considera factores relacionados con la fragilidad social y la falta de resiliencia que agrava las condiciones del riesgo físico.
Los autores quieren agradecer el apoyo del Ministerio de Educación y Ciencia de España para el «Enfoque integral y probabilista para la evaluación del riesgo sísmico en España», CoPASRE (CGL2011-29063), y al Ministerio de Economía y Competitividad en el marco del programa de formación de investigadores (FPI).