Introducción

A lo largo de los años, la ocurrencia de sismos en México, han dejado cuantiosos daños materiales y la irreparable pérdida de miles de vidas humanas. Son fenómenos naturales que nada podemos hacer para evitarlos y mucho menos predecirlos, por lo cual lo sociedad debe de habituarse a convivir con ellos y dentro de lo posible, estar preparados para mitigar sus potenciales daños. Precisamente uno de estos sismos, se presentó el 7 de septiembre de 2017 a las 23:49:17 horas del centro de México (04:49:17 UTC), con magnitud de 8.2, epicentro localizado a 140 km al suroeste de Pijijiapan, en el Estado Mexicano de Chiapas, a una latitud de 14.76° N y longitud -94.10° W, según lo reporto el Servicio Sismológico Nacional (SSN).

En el Estado Veracruz, los efectos dañinos generados por este evento sísmico, se apreciaron en los municipios sureños de Coatzacoalcos, Cosoleacaque, Hidalgotitlán, Jáltipan, Minatitlán y Texistepec, lo que obligó a la Secretaría de Gobernación ha emitir la declaratoria de desastre natural para estos 6 municipios (^{DOF, 20 de septiembre 2017}). De estos, el municipio con mayores reportes de daños fue Coatzacoalcos, donde viviendas de uno a dos niveles del sur de la ciudad, presentaron fallas asociadas al fenómeno de la licuación de suelos.

El objetivo principal de este trabajo es documentar los daños ocasionados a 21 viviendas por los asentamientos súbitos que se presentaron durante el sismo del 7 de septiembre de 2017 (M_w 8.2) y que son atribuibles al fenómeno de la licuación de los estratos arenosos que cubren zonas importantes de la ciudad de Coatzacoalcos, Veracruz, México, a partir de daños en viviendas y posterior evaluación del fenómeno mediante el criterio del National Center for Earthquake Engineering Research (NCEER) publicado por ^{Youd et al. (2001)} en donde evalúa el fenómeno mediante un Factor de Seguridad. Así mismo se espera que esta información coadyuve en la elaboración en un futuro no muy lejano del Reglamento de Construcción de la ciudad.

Como se puede apreciar en la Figura 1 el municipio de Coatzacoalcos se encuentra al sur del Estado de Veracruz, en México, a una latitud de 18.15° N y longitud -94.43º W, a una altura promedio de 10 metros sobre el nivel del mar (msnm). Limita al norte con el Golfo de México, al sur con los municipios de Nanchital, Cosoleacaque, Ixhuatlán del Sureste y Moloacán; al este con Agua Dulce y al oeste con Pajapan. La importancia del municipio, radica en que tiene una ubicación privilegiada con ventajas naturales y logísticas, lo cual permitió que se instalaran en décadas pasadas, una parte importante de la Petroquímica del país y recientemente para estimular su potencial de desarrollo el Gobierno de la República pretende construir el corredor interoceánico que conectaría el Golfo de México con el Océano Pacifico a partir de los puertos marítimos de Coatzacoalcos y Salina Cruz, este último ubicado en el estado de Oaxaca.

(M_w=8.2, Intensidad Instrumental ≈ IX)

Fuente:https://earthquake.usgs.gov/earthquakes/eventpage/us2000ahv0#map?ShakeMapStations=true&shakemapSource=us&shakemapCode=us2000ahv0

Figura 1 Ubicación de Coatzacoalcos, Veracruz, México, en el mapa de intensidades sísmicas de Mercalli Modificado, generado por el USGS del sismo del 7 de septiembre de 2017 

El sismo, que más daños ha generado a la infraestructura de la zona sur del estado de Veracruz, ocurrió el 26 de agosto de 1959, con una magnitud de 6.5, ubicado en el Golfo de México, con epicentro localizado a 35 km de la Ciudad de Coatzacoalcos a una latitud de 18.45° N y longitud -94.26° W, el cual generó rangos de aceleraciones superficiales de 0.07 a 0.2 de la gravedad (Cossio, 1960; ^{Marsal, 1961}). A este evento se le conoce como el sismo de Jáltipan, debido a los cuantiosos daños materiales y humanos generados a esta ciudad. Así mismo, en este evento, se documentó el primer caso de licuación de arenas en México (^{Marsal, 1961}). En la ciudad de Coatzacoalcos, este sismo generó daños en lo que hoy es el Astillero Número Tres de la Secretaría de Marina-Armada de México y en los muelles de lo que antes se llamaba Puertos Libres, hoy Administración Portuaria Integral (Figura 2). En el Astillero se presentaron daños en las cimentaciones de tres edificios, en los cuales se observaron asentamientos súbitos del orden de un metro en las columnas y en las paredes del taller mecánico. En la parte del muelle se observaron desplazamientos relativos horizontales de más de medio metro entre sus secciones. Estos daños fueron atribuidos al fenómeno de licuación parcial que se presentó en los mantos arenosos y limo -arenosos de esa zona (^{Marsal, 1961}).

Figura 2 Mapa de zonificación geotécnica de la zona urbana de la Ciudad de Coatzacoalcos, Veracruz, México, (^{Astudillo et al, 2010}), así como la vista general de zonas potencialmente licuables debido al sismo del 7 de septiembre de 2017 (M_w 8.2) 

Información geotécnica de la zona de estudio

En la última zonificación geotécnica propuesta por ^{Astudillo et al., (2010)}, se plantea que la zona urbana de la Ciudad de Coatzacoalcos está divida en cuatro zonas geotécnicas. En esta zonificación como se observa en la Figura 2, además de la zona geotécnica alta (ZGA) y la zona geotécnica baja (ZGB), como en un principio lo propusieron ^{Vieitez y Mosqueda (1970)}, incluye otras dos zonas; la zona geotécnica costera (ZGC) y la zona geotécnica de transición (ZGT).

Cada una de estas zonas, presentan características geotécnicas propias, como, por ejemplo; la zona geotécnica costera (ZGC) y la zona geotécnica alta (ZGA), están conformadas por un potente estrato de arenas de dunas y médanos de origen eólico, con espesores que van de los 15 a 30 metros de profundidad. Es escasa la información de los suelos que subyacen a las arenas de esta zona; sin embargo, se infiere que subyacen suelos arcillosos y limosos de origen marino o aluvial.

Figura 3 Ubicación de SPT, PVA y viviendas dañadas (VD) en las tres zonas potencialmente licuables 

La zona geotécnica de transición (ZGT), está formada por arenas de origen eólico igual que las dos zonas geotécnicas anteriores (ZGC y ZGA), diferenciándose, en el espesor inicial, el cual se reduce entre 4 y 12 metros. A partir de esta profundidad, se van alternando estratos cohesivos muy blandos o arenosos sueltos hasta llegar a los 20 metros de profundidad. Precisamente en estos estratos arenosos que se encuentran intercalados, se tienen características geotécnicas que pueden indicar vulnerabilidad a la licuación, así mismo empiezan a tomar relevancia los asentamientos generados por la consolidación de los estratos cohesivos blandos.

Por último, tenemos la zona geotécnica baja (ZGB), la cual tiene las mayores desventajas, debido a que está conformada por una secuencia de depósitos principalmente de origen fluvial (río Coatzacoalcos), constituidos en la parte superficial por un estrato de arenas de médanos (arenas limosas mal graduadas), de muy sueltas a media, con espesores que van de los 2 a 4 metros. A este suelo, le subyace un estrato de turba con espesor aproximado de 2 metros, de esa profundidad hasta los 38 metros, se localiza un estrato potente de limo orgánico de alta compresibilidad, consistencia blanda a media, plasticidad alta y olor fétido, intercaladas entre este gran estrato; bolsas de arena limosa y turba. A partir de esta profundidad y hasta los 45 metros aproximadamente, se localiza una arena limosa de compacidad densa a muy densa. En esta zona los asentamientos por consolidación son graves.

Metodología

Para la documentación de daños ocasionados a las viviendas por el sismo, se procedió a realizar en horas posteriores al evento, un censo, en donde se identificaron las viviendas que sufrieron daños tanto geotécnicos como estructurales, mediante su geolocalización, descripción de daños y reporte fotográfico, como se reporta en la Tabla 1 y Figuras 4 a la 8.

Figura 4 Vivienda dañada VD01 

Figura 5 Vivienda dañada VD05 

Figura 6 Vivienda dañada VD11 

Figura 7 Vivienda dañada VD12 

Figura 8 Vivienda dañada VD13 

Del censo realizado, se pudieron contabilizar 21 viviendas que resultaron seriamente dañadas, concentradas en tres zonas potencialmente licuables (Figura 2); en donde la mayoría de ellas presentaban un patrón común; durante el evento sísmico sufrieron asentamientos súbitos de gran magnitud en la cual emanó agua y arena del subsuelo (licuación de suelos). Este fenómeno se produce en medios granulares saturados y está asociado a un incremento en la presión de poros como resultado de la tendencia de las partículas a reacomodarse debido a las cargas cíclicas, lo que da como resultado una reducción en la resistencia al esfuerzo cortante, llegando a ser total, como fue en algunos de los casos que se reportan en este documento.

A continuación, se presenta el reporte fotográfico (casos representativos) de viviendas que sufrieron daños:

Para validar que los daños que se presentaron en las viviendas fueron debido al fenómeno descrito, fue necesario determinar la vulnerabilidad a la licuación de los estratos arenosos que subyacen en las zonas mencionadas, mediante el criterio del National Center for Earthquake Engineering Research (NCEER) publicado por ^{Youd et al. (2001)}, el cual compara la relación de resistencia cíclica (CRR) y la relación de esfuerzo cíclico (CSR). La primera correlaciona la resistencia a la licuación con las características del suelo expresadas a través del número de golpes de la prueba de penetración estándar (N_spt). Y la segunda está asociada a la aceleración máxima (a_máx.) que se presenta en el depósito de suelo durante el evento sísmico.

a) Sondeos geotécnicos.

Para contar con los perfiles estratigráficos de las zonas potencialmente licuables, se optó por ejecutar 5 sondeos de penetración estándar (SPT), ubicados como se muestra en la Figura 3a y 3b y con las características que se mencionan en la Tabla 2. Los perfiles estratigráficos obtenidos de los sondeos se presentan en las Figuras 16 al 20.

Tabla 2 Características de los sondeos de penetración estándar (SPT) e influencia de estos sobre las VD y PVA 

Sondeo	Coordenadas (UTM)		Prof. de Exploración	Influencia del sondeo sobre las VD	Influencia del sondeo sobre los PVA	Zona potencialmente licuable No.
	X	Y	(m)
SPT-01	342434.89	2005227.15	7.20	VD01 a la VD12	P23, P24	1
SPT-02	341972.00	2005257.00	9.60	VD14 y VD15	P25	1
SPT-03	344334.00	2004638.00	9.70	VD16 y VD17	P21	2
SPT-04	344200.00	2004741.00	9.70	VD18 y VD19	P21	2
SPT-05	342086.00	2005398.00	10.30	VD13	P25	1

UTM Universal Transversal de Mercator

VD Viviendas dañadas

PVA Puntos de vibración ambiental

Figura 9 Ejecución de los sondeos de penetración estándar 

b) Estimación de la aceleración máxima en superficie.

Las aceleraciones máximas en roca (a0r) estimadas con el programa PRODISIS (CFE-2015), tienen una diferencia mínima del 8% entre las ciudades de Minatitlán (mayores) y Coatzacoalcos (Figura 1); por ello se creyó apropiado emplear el acelerograma registrado en la estación MIHL (Figura 10) de Minatitlán, Ver., ubicada a 17 km de la zona en estudio y operada por el Instituto de Ingeniería de la UNAM¹. A este registro se le quitó el efecto de sitio para obtener el acelerograma en roca (AR), y posteriormente calcular los acelerogramas en superficie (ASC) en la zona de estudio, en la ciudad de Coatzacoalcos.

Figura 10 Registro del sismo de Tehuantepec (M_w 8.2) del 7 de septiembre de 2017 en estación MIHL¹ 

Es frecuente evaluar el efecto de sitio (^{Lermo et al., 1993}; ^{Chávez et al., 2014}), utilizando el cociente espectral de la componente horizontal con la vertical en el dominio de las frecuencias, para registros de sismos obtenidos en una estación sobre sedimentos, EHVSR, (Earthquake Horizontal to Vertical Spectral Ratio), cuando no se cuenta con una estación de referencia cercana, o que haya registrado los mismos eventos. El cociente espectral para el registro del sismo del 7 de septiembre de 2017 en la estación MIHL se determinó para cada dirección (N-S y E-W), y después se promediaron. Para estos cálculos se utilizó el programa DEGTRA (^{Ordaz et al., 2013}) y los resultados se muestran en la Figura 11a.

Figura 11 Determinación del acelerograma en roca (AR) 

Para eliminar el efecto de sitio al registro del sismo en superficie y determinar el acelerograma en roca (AR), se realizó un escalamiento entre la función de transferencia (EHVSR) determinada con el registro sísmico, y la componente horizontal del sismo que presentaba la mayor aceleración o energía. Primero, al acelerograma en superficie se le aplicó la transformada rápida de Fourier para obtener su respectivo espectro de amplitudes de Fourier (EAFS); de este último, sólo se considera el vector de amplitudes, pero no el de la fase (Figura 11b). Como la componente N-S del sismo presentó las mayores aceleraciones, se utilizó para determinar el acelerograma en roca (AR). Una vez que se tienen la EHVSR y el EAFS, en el dominio de las frecuencias, se calculó una “seudo” deconvolución (sólo se emplea el vector de amplitudes), dividiendo su respectivo efecto de sitio. Como segundo paso, al cociente se le aplica la transformada inversa de Fourier para obtener la señal en el dominio del tiempo, es decir el acelerograma en roca (AR) (Figura 12). El programa utilizado fue desarrollado en ambiente Matlab (J. Lermo Samaniego y D. Vega Rocha, comunicación personal, febrero 2011).

Figura 12 Acelerograma en roca (AR), amax = 33.8 cm/seg² 

c) Efecto de sitio en la zona de estudio

En la práctica, es común abordar los efectos de sitio en una de estas dos rutas: la primera consiste en estimar las frecuencias para obtener la amplificación y la amplitud de ésta, ya sea utilizando registros de sismos o de vibración ambiental. La segunda, aborda el problema determinando la geometría y las propiedades mecánicas de los materiales del subsuelo (^{Chávez-García et al, 2014}).

^{Lermo y Chávez-García (1993)} obtuvieron con ruido (vibración ambiental o microtremores), una buena estimación de la frecuencia de resonancia del suelo, usando la relación espectral entre los movimientos horizontales y verticales (H/V) de un mismo registro de ruido en un mismo sitio. Esto se ha estudiado ampliamente señalando sus ventajas y limitaciones, entre estas últimas se ha observado que la amplificación solo es una idea aproximada de la real. A la relación H/V se le han dado diversas explicaciones: una es suponer que los microtremores se componen principalmente de ondas superficiales, relacionándolas con las ondas Rayleigh, mientras que otros autores postulan que las ondas de cuerpo dan el pico máximo de la relación H/V.

Estudios recientes realizados por ^{Sánchez-Sesma et al., (2011b)}, proponen una teoría que sugiere que el cociente H/V (Técnica de Nakamura o MHVSR, Microtremor Horizontal to Vertical Spectral Ratio) es un estimador adecuado del nivel de amplificación del sitio, tal cociente, describe adecuadamente las principales propiedades dinámicas de suelo que interesan evaluar para determinar la respuesta del sitio (Efecto de sitio), y asume que el ruido es una onda de campo difuso que contiene todo tipo de ondas elásticas que dan una idea física de profundidad. Por ello, los cocientes espectrales H/V han mostrado ser útiles tanto para analizar registros de sismos como registros de vibración ambiental (^{Field y Jacob, 1995}). A partir de la publicación de ^{Lermo y Chávez-García (1994a)}, la técnica se popularizó y se han publicado una gran cantidad de artículos sobre el tema, lo que explica el éxito de una de las técnicas más utilizadas hoy en día.

Para este estudio se tomaron 13 registros de vibración ambiental en la zona de interés, distribuidos tomando en cuenta la ubicación de los daños en las viviendas. El criterio general de distribución fue, el de colocar una mayor densidad de puntos en los sitios que correspondan a zonas de depósitos sedimentarios con características de terrenos geotécnicamente blando (ZGB), infiriendo que en esta zona en particular se presentará una amplificación de ondas sísmicas, y así evaluar el efecto de sitio de los puntos P23, P24 y P25 (Figura 3a).

Los equipos de registro que se utilizaron para la campaña de vibración ambiental fueron dos sismómetros triaxiales digitales de banda ancha de la marca Guralp Systems CMG-6TD (6D37 y TD6618), los cuales cuentan con tres canales de registro de velocidades ortogonales, así como accesorios periféricos como GPS, batería externa de respaldo, cables y otros. Fueron orientados conforme al norte magnético y la duración de los registros fueron de una hora en cada sitio, con un muestreo de 100 muestras por segundo. Los registros de vibración ambiental fueron corregidos por línea base y filtrados con un filtro pasa bandas de 0.1 a 10 Hz y un método de filtrado de Butterworth de cuarto orden. Una vez acondicionada la señal se realizó el cociente espectral entre la raíz cuadrada del promedio de los espectros de amplitudes de Fourier de las señales horizontales y el espectro de amplitudes de Fourier de la señal vertical, al resultado se le aplicó un suavizado de Konno & Omachi con un factor de suavizado de 40; todos estos procesos se hicieron con el programa GEOPSY para así obtener los cocientes espectrales H/V que se muestran en la Figura 13a.

Figura. 13 Cocientes espectrales H/V 

En la zona en estudio se identificaron tres zonas potencialmente licuables y en la Tabla 3 se relacionan las zonas con los puntos donde se tomaron registros de vibración ambiental. En la Figura 13b, se indican los cocientes espectrales de los puntos incluidos en las zonas y el punto P16 donde hubo licuación en el sismo de Jáltipan de 1959.

Tabla 3 Vibración ambiental vs zonas potencialmente licuables 

Puntos	Frecuencia (Hz)			Periodo (seg)			Amplificación			Zona
P23, P24, P25	0.97	1.02	1.07	1.03	0.98	0.94	4.03	3.04	9.87	1
P21	0.94			1.06			5.63			2
P15	0.77			1.30			3.6			3

Se estimó el efecto de sitio en los lugares donde se tomaron registros de vibración ambiental, y con ellos se determinaron los acelerogramas en superficie. Para lo anterior, se realizó una “seudo” convolución entre el espectro de amplitudes de Fourier (EAFR) del AR, con los cocientes espectrales (H/V) de las señales de ruido sísmico registradas (HVNR), para ello se realizó el producto del EAFR con H/V y así considerar las amplificaciones por el efecto de sitio en la señal. Finalmente se aplicó una transformada inversa de Fourier para obtener la señal en el dominio del tiempo, es decir el acelerograma en superficie ASC (acelerogramas sintéticos).

En la Figura 14a se muestra el espectro de amplitudes del AR, línea azul, y los cocientes espectrales H/V del punto 25, línea roja, por ser el sitio en donde se presentan las mayores amplificaciones. Los acelerogramas del sitio se muestran en la Figura 14b, en donde en azul aparece el AR y en rojo el acelerograma en superficie ASC, apreciándose de esta manera la amplificación que sufre la señal al propagarse hasta llegar a la superficie.

Figura 14 Efecto de sitio en el punto P25 

Este mismo procedimiento se siguió para los demás puntos, y en la Tabla 4 su muestran las aceleraciones en superficie en los puntos incluidos en las zonas potencialmente licuables. También se muestran las aceleraciones en el punto P16, que corresponden a un punto cercano a la zona donde se presentó licuación de arenas durante el sismo de Jáltipan de 1959. El valor obtenido en este punto P16, corresponde a 0.064 de la gravedad (62.5 gales), valor muy cercano al mínimo dentro del rango (0.07 a 0.2 de la gravedad) que reportan Cossio, 1960; ^{Marsal, 1961} para el sismo de Jáltipan.

Tabla 4 Aceleraciones en superficie (gales) por zonas 

Zona	Sismo de Tehuantepec 7 septiembre 2017 (Mw 8.2)			Para diseño de acuerdo a PRODISIS 1				Punto P163
				S1	S2	S3	Promedio	f (hz)	T (seg)	A
1	P23	P24	P25	241.5	280.8	253.6	258.62	1.3	0.33	3.6
	67.6	55.4	141.0
2	97.1			182.1	214.6	184.2	193.6	62.5
3	79.4			123.7	115.1	146.7	128.5

¹ Considerando 3 acelerogramas sintéticos

² Se evalúo en el punto P25

³ Zona donde se presentó licuación en el sismo de Jáltipan 1959

Para tener una referencia del nivel de aceleraciones que se presentaron en la zona de estudio durante el sismo del 7 de septiembre de 2017, se hizo una comparación entre las aceleraciones en superficie que se alcanzaron durante este sismo, y las aceleraciones de diseño que se pueden presentar cuando se emplean las máximas aceleraciones probables, considerando el peligro sísmico de acuerdo con el PRODISIS, Tabla 4. Para ello se obtuvieron del PRODISIS tres acelerogramas sintéticos cuyos espectros de respuesta se ajustan al espectro de peligro uniforme (espectro de referencia, ER) para la zona de estudio, Figura 15. Como los acelerogramas sintéticos de la Figura 15b, representan al movimiento del suelo basal, es decir en roca, fue necesario realizar la “seudo” convolución (siguiendo el procedimiento descrito anteriormente) entre el espectro de amplitudes de Fourier de los acelerogramas sintéticos y el producto espectral HVNR, que se obtuvo de ruido sísmico, en los puntos de estudio.

Figura 15 Obtención de acelerogramas sintéticos empleando PRODISIS 

De la Tabla 4 podemos observar que, a pesar de que las aceleraciones en superficie del sismo fueron inferiores a las de diseño de acuerdo con PRODISIS, éstas provocaron que el suelo presentara fallas por licuación. Esto es una señal clara de la disposición que tiene el suelo de esta zona para presentar este tipo de falla, y el hecho de que las aceleraciones de diseño no se alcanzaran nos indican que, de hacerlo, la zona de daño por esta causa pudo haber sido mucho mayor.

Vulnerabilidad a la licuación

El método del National Center for Earthquake Engineering Research (NCEER) publicado por ^{Youd et al. (2001)}, evalúa la vulnerabilidad a la licuación de estratos arenosos, por medio de un factor de seguridad contra la licuación y la define como:

En donde CRR representa la capacidad de resistencia del suelo a licuarse y CSR representa la fuerza sísmica que actúa en el suelo. En teoría, cuando el valor del factor de seguridad (FS_lic) es menor que la unidad, se presenta el fenómeno a la profundidad analizada. En la práctica existe incertidumbre acerca del valor del FS_lic que debe utilizarse, ya que han existido casos en los cuales el suelo se ha licuado a pesar de tener un FS_lic>1 (^{Youd e Idriss, 2001}). Por lo anterior, el Building Seismic Safety Council (BSSC, 2003) recomienda que el FS_lic debe de estar sujeto a la importancia de la estructura y potencial de deformación del suelo. Martin y Lew (1999) sugirieron valores del FS_lic asociados al número de golpes de la prueba SPT, en donde para valores de N_spt<15, recomienda utilizar un FS_lic de 1.2. Esta última consideración fue adoptada para nuestro análisis.

Resultados

A continuación, se presentan los resultados obtenidos del análisis de la vulnerabilidad a licuación de cada uno de los sondeos;

1) Sondeo SPT-01: Como se puede observar en la Figura 16, durante el evento sísmico, se estimó que se licuó un espesor de arena de aproximadamente 1.3 m (profundidad -2.70 al -4.00). Esto se puede atribuir; a que en este nivel se tiene una arena limosa en condiciones saturadas, cuya compacidad transita de suelta a muy suelta y números de golpes N_SPT ≤ 7. Estas condiciones geotécnicas desfavorables predominan en esa parte de la ciudad, lo que puede explicar las amplificaciones del evento sísmico y los subsecuentes asentamientos súbitos que se presentaron en las viviendas VD01 a la VD12.

Figura 16 Perfil estratigráfico y resultados del análisis a la vulnerabilidad a la licuación del sondeo SPT-01 

2) Sondeo SPT-02: Como se puede observar en la Figura 17, durante el evento sísmico, se estimó que se licuaron dos estratos; el primero de un espesor aproximado de 1.8 m (profundidad -3.00 al -4.80), en donde se localiza una arena limosa saturada, cuya compacidad es suelta, con números de golpes N_SPT ≤ 6., el segundo estrato que se licuó, corresponde a un espesor de 1.2 m (profundidad -6.60 al -7.80), en la cual se tiene igualmente una arena limosa saturada, con una compacidad muy suelta y números de golpes N_SPT ≤ 4. En esta parte de la ciudad sufrieron asentamientos súbitos las viviendas VD14 y VD15.

Figura 17 Perfil estratigráfico y resultados del análisis a la vulnerabilidad a la licuación del sondeo SPT-02 

3) Sondeo SPT-03: Como se puede observar en la Figura 18, durante el evento sísmico, se estimó que se licuaron parcialmente dos estratos; en el primer estrato, que es una arena limosa mal graduada, la cual transita de una compacidad densa a suelta, se licuaron 0.4 m en su parte inferior (profundidad -2.60 al -3.00) y el segundo estrato que se licuó, corresponde a una arena limosa, cuya compacidad transita de media a muy suelta, con número de golpes N_SPT ≤ 15., el colapso se presentó en su parte superior (profundidad -6.10 al -7.30) como inferior (profundidad -7.90 al -9.30). En esta parte de la ciudad sufrieron asentamientos súbitos las viviendas VD16 y VD17.

Figura 18 Perfil estratigráfico y resultados del análisis a la vulnerabilidad a la licuación del sondeo SPT-03 

4) Sondeo SPT-04: Como se puede observar en la Figura 19, durante el evento sísmico se estimó que se licuaron dos estratos; en el primero estrato, se tiene una arena limosa mal graduada, con una compacidad muy suelta y número de golpes N_SPT ≤ 4, se licuó aproximadamente 2.7 m (profundidad -0.60 al -3.00). El segundo estrato que se licuó corresponde a una arena limosa, la cual transita de una compacidad muy suelta a media. En este estrato, el fenómeno se presentó tanto en la parte superior (profundidad -6.10 al -7.90) como inferior (profundidad -8.20 al -9.20). En esta parte de la ciudad sufrieron asentamientos súbitos las viviendas VD18 y VD19.

Figura 19 Perfil estratigráfico y resultados del análisis a la vulnerabilidad a la licuación del sondeo SPT-04 

5) Sondeo SPT-05: Como se puede observar en la Figura 20, durante el evento sísmico se estimó que se licuaron tres estratos; el primero de un espesor aproximado de 2.3 m (profundidad -1.30 al -3.60), en donde se localizó una arena limosa saturada transitando de una compacidad media a muy suelta en su parte inferior, el segundo estrato que se licuó, tiene un espesor de 1.3 m (profundidad -4.80 al -6.10), en el cual se localizó una arena limosa mal graduada saturada, con una compacidad suelta a muy suelta y por último, en el tercer estrato solo se licuó un espesor pequeño de 0.7 m (profundidad -9.00 al -9.70) en la cual igualmente se tiene una arena limosa, con una compacidad muy suelta. En los tres estratos en donde se presentó la licuación se tienen números de golpes N_SPT ≤ 7. En este punto se localizó la vivienda VD13, la cual fue la que más daños presentó durante el evento símico.

Figura 20 Perfil estratigráfico y resultados del análisis a la vulnerabilidad a la licuación del sondeo SPT-05 

De los 5 sondeos realizados para este trabajo, este sondeo en particular, presenta las condiciones geotécnicas más desfavorables, esto se atribuye a que, esta zona corresponde a terrenos ganados al pantano, lo que puede explicar las amplificaciones del movimiento sísmico (efecto de sitio) y el subsecuente hundimiento de la vivienda ya mencionada, que llego hacer de hasta 1.2 m, tal como se puede visualizar en la Figura 8.

Conclusiones

En este trabajo se presenta un análisis de daños ocasionados en viviendas ubicadas en la ciudad de Coatzacoalcos, Veracruz, México., por el sismo de magnitud 8.2 (M_w) que se suscitó el 7 de septiembre de 2017, en las costas de Pijijiapan, Chiapas. Se revisaron los daños en 21 viviendas, localizadas principalmente en la zona sur de la ciudad (ZGT y ZGB), las cuales presentaron asentamientos súbitos atribuibles al fenómeno de licuación parcial de los estratos arenosos que subyacen en esta zona. Para lo anterior, se realizó un censo de daños en horas y días posteriores al evento sísmico, se ejecutaron 5 sondeos de penetración estándar (SPT) y se llevó a cabo una campaña de vibración ambiental. Del análisis de resultados se concluye lo siguiente: