RGAC
Número 70(1) • Julio-diciembre 2023
ISSN: 1011-484X • e-ISSN 2215-2563
Doi: https://dx.doi.org/10.15359/rgac.72-2.12
Recibido: 01/03/2023 • Aceptado: 30/08/23
URL: www.revistas.una.ac.cr/index.php/geografica/
Licencia (CC BY-NC-SA 4.0)

Análisis geocientífico de colapsos estructurales en la ciudad de Juan C. Bonilla, Puebla, México (2021)

Geoscientific Analysis of Structural Collapses in the City of Juan C. Bonilla, Puebla, Mexico (2021)

Análise Geocientífica do Colapso Estrutural na Cidade De Juan C. Bonilla, Puebla, México (2021)

Rogelio Ramos Aguilarn1

Benemérita Universidad Autónoma de Puebla, Puebla,México

Patricia Máximo Romero2

Benemérita Universidad Autónoma de Puebla, Puebla, México

Máximo Ávila Cruz3

Benemérita Universidad Autónoma de Puebla, Puebla, México

Víctor Galindo López4

Benemérita Universidad Autónoma de Puebla, Puebla, México

María de los Dolores Ramírez Álvarez5

Benemérita Universidad Autónoma de Puebla, Puebla, México

Lorena Cárdenas López6

Benemérita Universidad Autónoma de Puebla, Puebla, México

Argenis Eduardo Morales González7

Benemérita Universidad Autónoma de Puebla, Puebla, México

María de la Cruz Vázquez García8

Imágenes Geográficas S.A. México

Victoriano Covarrubias Salvatori9

Consejo de Ciencia y Tecnología del Estado de Puebla, México

Resumen

El trabajo se basa en un análisis general desde las geociencias del colapso en los campos de cultivo de Santa María Zacatepec, de la ciudad Juan C. Bonilla. En el desarrollo se aplica la técnica de Percepción Remota, permitiendo interpretar combinaciones de bandas y determinar la geometría, geología, geomorfología y cartografía del sitio. Se presenta un análisis resumido del trabajo que se realizó, además de mostrar modelos de simulación del gasto hidrológico y la aceleración del agua en las pendientes estructurales de escorrentía, datos que se usaron para la cartografía y modelado gráfico finales. Se realizó un estudio geomorfológico, geológico y topográfico del lugar, con técnicas telemétricas satelitales y aéreas; también se utilizó software especializado que se aplicó para levantamientos geodésicos, realizando cartografía en modelos 2D y 3D del área de análisis.

Las prospecciones del lugar se realizaron con el objetivo de conocer la forma de la estructura e identificar áreas potencialmente riesgosas para la población cercana al colapso. Se obtuvieron ortofotos para calcular y determinar la pendiente por distancia-elevación, la aceleración-altura, el gasto y la fuerza del agua; que sirvió para concluir las posibles causas que ocasionaron el colapso, y contar con datos técnicos suficientes que se compartieron con las autoridades de protección civil.

Palabras clave: colapso, geomorfología, geoestadística, aceleración, modelos.

Abstract

This paper presents a general analysis from the geosciences of structural collapses. The analysis was conducted in the crop fields of Santa María Zacatepec, in the city of Juan C. Bonilla, Mexico. The analysis implemented the Remote Sensing technique, facilitating the interpretation of band combinations and determining the geometry, geology, geomorphology and cartography of the site under study. A summary analysis of the work performed is presented, in addition to showing simulation models of the hydrological outflow and water acceleration in the runoff structural slopes, data that were used for the final mapping and graphic modeling. A geomorphological, geological, and topographical study of the site was conducted, using satellite and aerial telemetric techniques. Specialized software was also used for geodetic surveys, mapping in 2D and 3D models of the analysis area. Surveys of the site were carried out to determine the shape of the structure and identify potentially risk areas for the population near the collapse events. Orthophotos were obtained to calculate and determine the slope by distance-elevation, acceleration-height, water flow and force, which served to conclude the possible causes that led to the collapses and have sufficient technical data that were shared with the civil protection authorities.

Keywords: collapse, geomorphology, geostatistics, acceleration, models

Resumo

O estudo fundamenta-se em uma análise geral das geociências do colapso nos campos de cultivo de Santa María Zacatepec, na cidade de Juan C. Bonilla, México. A técnica de sensoriamento remoto é aplicada no desenvolvimento, permitindo a interpretação de combinações de bandas e a determinação da geometria, geologia, geomorfologia e cartografia do local. É apresentada uma análise resumida do trabalho realizado, além de modelos de simulação do gasto hidrológico e da aceleração da água nos taludes estruturais de escoamento. Esses dados foram utilizados para o mapeamento final e a modelagem gráfica. Foi realizado um estudo geomorfológico, geológico e topográfico do local, utilizando técnicas telemétricas aéreas e de satélite. Também foi utilizado um software especializado para levantamentos geodésicos e mapeamento em modelos 2D e 3D da área de análise. Os levantamentos do local foram realizados com o objetivo de conhecer o formato da estrutura e identificar áreas de risco potencial para a população próxima ao colapso. Foram obtidas ortofotos para calcular e determinar a inclinação por distância-elevação, aceleração-altura, fluxo de água e força, o que serviu para concluir as possíveis causas que levaram ao colapso e para ter dados técnicos suficientes que foram compartilhados com as autoridades de proteção civil.

Palavras-chave: colapso, geomorfologia, geoestatística, aceleração, modelagem

Agradecimientos

Se agradece el apoyo del Consejo de Ciencia y Tecnología del estado de Puebla para llevar a cabo el estudio, bajo convenio de colaboración CONCYTEP/214/2021.

Introducción

Los hundimientos geológicos existen en todo nuestro planeta, se originan en algunas ocasiones por excavación, filtración de agua o falta de ella, por un peso extra en la superficie, presentándose principalmente en zonas kársticas. La investigación se centra en el colapso estructural que se generó el 29 de mayo de 2021, en los campos de cultivo ubicados en la localidad de Santa María Zacatepec, perteneciente al municipio de Juan C. Bonilla, Puebla, México, dentro de la cuenca hidrológica del Río Alto Atoyac región hidrológica RH18 (INEGI, 2020).

Se realizaron estudios: geomorfológicos, geohidrológicos y geoestadísticos, con el objetivo de contar con datos e información estructural relacionada con la topografía, geología, geomorfología, geohidrología, cartografía y percepción remota del lugar; para comprender las posibles causas que originaron el colapso estructural (Máximo et al, 2017).

El trabajo se realizó con técnicas de prospección de campo y telemétricas, el cual se complicó, debido a que se desarrolló durante la temporada de lluvia, ya que el procesamiento de algunas de las imágenes satelitales, contaban con nubosidad arriba del 30%, lo que hacía imposible dar un seguimiento (Ramos et al, 2013).

Se realizaron cálculos de la prospección de campo y se compararon con los datos telemétricos para determinar si la precipitación de lluvia acumulada pudiera haber influido en el colapso.

Se encontró que la estructura colapsada pone en riesgo a la población que se encuentra cerca, ya que los resultados determinaron que probablemente exista un cambio paulatino de deformación, también se han generado pequeñas fisuras superficiales junto a la zona de estudio que suponen un mecanismo de falla continúa.

Los datos obtenidos, concluyen que la geodinámica terrestre y atmosférica de la zona, influyeron en el hundimiento geológico.

Antecedentes

El 29 de mayo de 2021 se presenta un colapso estructural en la localidad de Santa María Zacatepec del municipio de Juan C. Bonilla, ubicada a aproximadamente a 50 km de la capital del estado de Puebla. El trabajo de campo y gabinete cuenta con datos y salidas a la zona que abarcan hasta junio 26 de 2021; se determina el área de la superficie del hundimiento para realizar el desarrollo del cálculo del líquido que se desplaza dentro del mismo, se generaron modelos topográficos y geométricos y se obtuvo la media de la zona estudiada de 11,281.46 m2, cuyo dato sirvió para calcular el área interior del colapso en un modelo 3D, dando un valor promedio de 32,100.79 m2.

Para el análisis, se relaciona el modelo del geoide terrestre con el elipsoide terrestre, para poder desarrollar los modelos geoestadísticos y de riesgo geohidrológicos y geomorfológicos y graficar e interpretar los resultados.

Se utilizaron datos telemétricos de plataformas satelitales como RapidEye, SkySat y Landsat, con software especializado para procesar la información y obtener modelados cartográficos y gráficos en 2D y 3D; las coordenadas utilizadas son geográficas y UTM. Se localizaron 18 puntos como referencia para generar el estudio gravimétrico, geoestadístico y de gasto.

Los datos de longitud, latitud y altura elipsoidal se convirtieron a una altura ortométrica que es la distancia entre la superficie física de la Tierra y del geoide en forma vertical, se obtuvo con el software libre del Instituto Nacional de Estadística y Geografía (INEGI) la curva tangente a la dirección de la gravedad en cualquier punto (Ramos et al, 2012).

La localidad de Santa María Zacatepec se ubica en el municipio de Juan C. Bonilla, en el estado de Puebla, se localiza en la parte centro-oeste del estado. Sus coordenadas geográficas son: 19º 05’ 30’’ y 19 º 08’ 36’’ de latitud norte y 98º 18’ 24’’ y 98º 25’ 36’’ de longitud oeste (Figura 1).

Figura 1. Ubicación de localidad de Santa María Zacatepec, Juan C. Bonilla, Puebla

Elaboración propia

Metodología

Se aplica la aceleración gravitacional normal para obtener la aceleración observada en el punto, el modelo utilizado es el EIGEN-6C4, a continuación, se calcularon las reducciones de Bouguer y aire libre (incluyendo la eliminación de la masa entre el punto y el geoide e igualmente la reducción de aire libre) para calcular sus respectivas gravedades y anomalías.

En el proceso de imágenes satelitales se combinaron los canales y bandas de los sensores satelitales para determinar algunas características de la región; se utilizaron datos multitemporales de las plataformas SkySat, Landsat y RapidEye para analizar la geología, topografía y uso de suelo.

El procesamiento de imágenes reveló algunas de las características que rodean a la zona colapsada, identificándose varias zonas de cultivo con condiciones de suelo blando y geología de roca volcánica.

Cálculo gravimétrico de las anomalías

Durante el desarrollo gravimétrico (Ramos et al, 2021) se realizó un análisis de los datos obtenidos para determinar la diferencia entre los campos gravitatorios terrestres reales y normales del área de estudio, la anomalía o perturbación potenciales (T), diferencia entre potenciales de gravedad (W) y energía de gravedad normal (U), el modelo matemático se muestra a continuación (1).

T(x,y,z)=W(x,y,z)-U(x,y,z) (1)

Para el desarrollo de los cálculos se utilizó el modelo de referencia geodésico GRS80, que es un sistema formado por el elipsoide de referencia global y del campo gravitatorio (el geoide), en la práctica, los elipsoides de referencia GRS80 y WGS84 son los mismos; se aplicó la Fórmula Internacional de Gravedad Normal adoptada por la Sociedad Internacional de Geodesia (2).

γ = 9.780318 (1 + 0.0053024 sen2 φ – 0.0000058 sen22 φ) (2)

Donde φ es la latitud y γ la gravedad normal dada en m/s2.

Se aplicaron para el estudio modelos prácticos matemáticos (3), donde F es la reducción por aire libre en mgal (milésima de gal, que equivale a 1 cm/s2) y H es la altura.

F=0.000003086*H (3)

La anomalía de aire libre (ΔgAL), utiliza la gravedad observada (gobs), la gravedad normal (γ) y la reducción de aire libre (F) (4).

ΔgAL=gobs+F-γ (4)

En este análisis se utilizó el valor de ρ=2.67 gr/cm3, modelo matemático (5).

AB=0.1119*H (5)

La anomalía de gravedad de Bouguer representa los cambios de grosor en la corteza terrestre y la densidad a nivel cortical, modelo (6).

ΔgB=gobs-γ-F+A+B (6)

Los modelos matemáticos generales para calcular las anomalías son (7,8, 9, 10, 11, 12):

● Atracción Gravitatoria de la capa de Bouguer:

AB=2πGρh=0.1119 h mgales (7)

● Gravedad reducida por Bouguer:

gB=gAl+At-AB (8)

● Anomalía de Bouguer:

ΔgtB ΔgAl+At-AB (9)

ΔgB ΔgAl+At-AB (10)

ΔgP=g+δgF+δgP-γ (11)

δgF+δgP=0.1967*h (12)

Con los datos de la Tabla 1 se aplicaron los modelos matemáticos anteriores (tabla utilizada en el reporte presentado al Consejo de Ciencia y Tecnología del Estado de Puebla)

Tabla 1. Datos calcular y obtener las anomalías del aire, aire libre y Bouguer

Elaboración propia.

Resultados y análisis

Cartografía Satelital

Para los datos e imágenes multitemporales se emplearon las plataformas SkySat, Landsat y SkySat, para los cálculos geométricos se utilizó el software Math. En la Figura 2 y Figura 3 se muestra la combinación de bandas 4, 3, 2 (RGB) con resolución de 50 cm., registrándose un área promedio inicial de 11,281.46 m2 y un perímetro de 328.09 m. en la zona del colapso (Ramos et al, 2022).

Las coordenadas UTM utilizadas para el procesamiento de las imágenes son las siguientes:

Punto de control

Coordenada X

Coordenada Y

1

565780.42 E

2114947.23 N

Esquina superior izquierda

2

565940.12 E

2114948.50 N

Esquina superior derecha

3

565943.70 E

2114812.18 N

Esquina inferior derecha

4

565780.19 E

2114809.72 N

Esquina inferior izquierda

Figura 2. Área procesada con MultiSpec de imágenes de Landsat, RapidEye y SkySat de zona de colapso

Elaboración propia.

Figura 3. Área procesada con MultiSpec de imágenes de Landsat, RapidEye y SkySat de combinación de bandas en la zona de colapso

Elaboración propia.

Durante el desarrollo del estudio, se encontró que la geomorfología estructural colapsada se fue modificando en un lapso aproximado de un mes, observándose con los registros de las imágenes satelitales procesadas durante ese tiempo. La geodinámica terrestre y atmosférica de la zona, contribuyó para un mecanismo de modificación del lugar (Ramos et al, 2005).

Se analizaron y procesaron las variables meteorológicas de precipitación como factor principal de alteración geohidrológica, las salidas a campo sirvieron para corroborar con los sensores meteorológicos instalados las lecturas obtenidas con el desarrollo matemático analizado.

En las gráficas 1, 2 y 3 se muestran los datos del desarrollo estructural (perímetro, área y volumen) obtenidos por análisis espectral satelital, del 17 de mayo al 26 de junio.

Gráfica 1. Crecimiento del área en m2

Elaboración propia.

Gráfica 2. Crecimiento del perímetro en m

Elaboración propia.

Gráfica 3. Crecimiento del volumen en m3

Elaboración propia.

Estudio geológico y geomorfológico

Se identificaron las características de la zona de estudio en cuanto a sus geoformas encontradas (relieves y topografía), que indican el proceso o grado de evolución generado por procesos endógenos y exógenos, además del proceso antropogénico.

Se identificó que la zona desplomada ha sido utilizada en exceso para el desarrollo agrícola, extracción de materiales y agua, durante los últimos 100 años; el relieve identificado, muestra un paisaje plano producto de meteorización, generado por lluvia y viento principalmente; la lluvia y la radiación solar han contribuido al cambio de la geología dentro del proceso evolutivo.

Las geoformas identificadas son pliegues anticlinales y sinclinales producto de fallamiento originados por procesos tectónicos y volcánicos, originados por contacto dinámico, describiendo sus características en la Figura 4.

Figura 4. Características geológicas y geomorfológicas de la zona colapsada con procesamiento ortofotogramétrico: 1.-Efecto de compresión de la roca, 2.-Existe una importante cantidad de material metamórfico intemperizado, 3.-Posible deslizamiento por flujo, 4.-El almacenamiento de sedimentos suaves es producto de la agricultura, 5.- Se aglomeran las formaciones geológicas y planicies creadas por los procesos fluviales del río Actiopa

Elaboración propia.

Fotogrametría terrestre

En las prospecciones realizadas durante el análisis de la zona colapsada se realizó fotogrametría terrestre, algunas imágenes se presentan en la Figura 5.

Figura 5. Fotogrametría terrestre de las prospecciones realizadas durante el análisis de la zona colapsada

Elaboración propia.

Aceleración de la gravedad con respecto al agua, anomalías de aire-libre, Bouguer e incremento de la gravedad

En la Tabla 2 se muestran los resultados de la aceleración de la gravedad del agua (m/s2 y en mgal) y las anomalías de Bouguer y aire-libre, las dos últimas columnas son los resultados del incremento de la aceleración de la gravedad.

Tabla 2. Resultados de la aceleración de la gravedad del agua (m/s2 y en mgal) y las anomalías de Bouguer y aire-libre

Vértice

Latitud

Elevación

Gravedad (g) m/s2

Gravedad (g) (mgal)

Gravedad Reducida por Aire-Libre (mgal)

Módulo de la Gravedad Normal en el Elipsoide (m/s2)

Módulo de la Gravedad Normal en el Elipsoide (mgal)

Anomalía Aire-Libre (mgal)

Atracción Gravitatoria a la capa de Bouguer (mgal)

Gravedad Reducida por Bouguer (mgal)

Anomalía de Bouguer (mgal)

Anomalía de Bouguer (mgal)

1

19.125555

2224

9.77899995

977899.995

978586.3214

9.784544517

978454.4517

131.869651

248.8656

978337.4558

-116.995949

-116.995949

2

19.125447

2224

9.77899989

977899.989

978586.3154

9.784544397

978454.4397

131.8756848

248.8656

978337.4498

-116.9899152

-116.9899152

3

19.125363

2224

9.77899984

977899.984

978586.3104

9.7845443

978454.43

131.880378

248.8656

978337.4448

-116.985222

-116.985222

4

19.12523

2224

9.77899976

977899.976

978586.3024

9.784544146

978454.4146

131.8878088

248.8656

978337.4368

-116.9777912

-116.9777912

5

19.124913

2224

9.77899959

977899.959

978586.2854

9.784543799

978454.3799

131.9055185

248.8656

978337.4198

-116.9600815

-116.9600815

6

19.124813

2224

9.77899953

977899.953

978586.2794

9.784543683

978454.3683

131.9111055

248.8656

978337.4138

-116.9544945

-116.9544945

7

19.124924

2223

9.77900268

977900.268

978586.2858

9.784546897

978454.6897

131.5961291

248.7537

978337.5321

-117.1575709

-117.1575709

8

19.125297

2223

9.77900289

977900.289

978586.3068

9.784547315

978454.7315

131.5752903

248.7537

978337.5531

-117.1784097

-117.1784097

9

19.125563

2223

9.77900304

977900.304

978586.3218

9.784547614

978454.7614

131.5604291

248.7537

978337.5681

-117.1932709

-117.1932709

10

19.125647

2223

9.77900308

977900.308

978586.3258

9.784547701

978454.7701

131.5557365

248.7537

978337.5721

-117.1979635

-117.1979635

11

19.12553

2223

9.77900302

977900.302

978586.3198

9.784547575

978454.7575

131.5622729

248.7537

978337.5661

-117.1914271

-117.1914271

12

19.124866

2223

9.77900265

977900.265

978586.2828

9.784546834

978454.6834

131.5993693

248.7537

978337.5291

-117.1543307

-117.1543307

13

19.124952

2223

9.7790027

977900.27

978586.2878

9.784546932

978454.6932

131.5945646

248.7537

978337.5341

-117.1591354

-117.1591354

14

19.125097

2223

9.77900278

977900.278

978586.2958

9.784547093

978454.7093

131.5864638

248.7537

978337.5421

-117.1672362

-117.1672362

15

19.125277

2223

9.77900288

977900.288

978586.3058

9.784547294

978454.7294

131.5764076

248.7537

978337.5521

-117.1772924

-117.1772924

16

19.124835

2223

9.77900263

977900.263

978586.2808

9.784546797

978454.6797

131.6011013

248.7537

978337.5271

-117.1525987

-117.1525987

17

19.125266

2223

9.77900287

977900.287

978586.3048

9.784547278

978454.7278

131.5770223

248.7537

978337.5511

-117.1766777

-117.1766777

18

19.125091

2223

9.77900277

977900.277

978586.2948

9.78454708

978454.708

131.5867994

248.7537

978337.5411

-117.1669006

-117.1669006

Elaboración propia.

Cálculo y desarrollo del gasto hidrológico y cartográfico

El aforo o suma es el gasto hidrológico de un fluido que se traslada por un canal de forma geométrica regular o irregular, calculado en metros cúbicos sobre segundo (m3/s).

Para el análisis y estudio, se establecieron varios modelos matemáticos y geométricos terrestres, examinando diversos escenarios como límites de agua bajas y superficiales.

Para determinar el gasto del agua en la zona colapsada se aplicó el modelo general de gasto (13).

Q=AIC/360 (13)

Donde A = superficie de la cuenca en metros cuadrados, I = intensidad de la lluvia en milímetros por hora y C = es el coeficiente de escorrentía sin dimensiones.

Los datos se obtuvieron del Instituto Nacional de Estadística y Geografía (INEGI) con el Simulador de Flujos de Agua de Cuencas Hidrográficas (SIATL, 2019), las variables utilizadas dificultaron el desarrollo debido a la falta de datos meteorológicos en la zona, para el procesamiento de los datos se tomó con un modelo ideal geohidrológico de planicie; los datos se manejaron para la precipitación en mm/hr, se multiplicaron por 24 horas y después por 30 días.

Se tomó como base longitudinal del cauce el diámetro promedio del hundimiento geológico; cabe mencionar que para el estudio se tomaron en cuenta modelos ideales geométricos para calcular el gasto, tomando al colapso dentro del área de la subcuenca para fines prácticos.

Los datos de campo se compararon con los obtenidos por telemetría y estaciones terrestres meteorológicas.

Este proceso se realizó para los 18 puntos, tomado los datos promedios; en la tabla 3 se muestra el cálculo superficial para una escorrentía con un modelo de textura geomorfológica plana, la intensidad de la lluvia se toma como punto de referencia media para realizar el proceso, se realizaron una serie de cálculos con diversos valores para la intensidad de lluvia, el criterio tomado es utilizar escurrimientos parciales ya que se carece de datos previos meteorológicos, por lo que se tuvo que tomar lecturas de estaciones meteorológicas portátiles instaladas en 2 escuelas en el municipio de Calpan, Puebla, muy cercanas a la zona de estudio.

En la Tabla 3 se muestra el resultado del promedio superficial para el mes de junio de 2021, con intensidad de lluvia baja y escorrentía con factor fijo de 0.30, el gasto obtenido es de 0.00000273575 m3/s.

Tabla 3. Resultado del promedio superficial para el mes de junio de 2021

Datos para cálculo de promedio superficial de gasto de zona colapsada

Intensidad lluvia

1.05 mm/hr es igual a 0.0000002910000 m/s

Área subcuenca

11,281.46 m2

Coeficiente escorrentía

0.30

Longitud cauce

126 m

Gasto promedio superficial

0.00000273575 m3/s

Elaboración propia.

En la Gráfica 4 se aprecia el comportamiento del gasto en los primeros 4 puntos calculados, el punto 1 representa el lugar más bajo por lo que el gasto es mínimo, el punto 4 representa una ligera elevación por lo que el gasto respecto a la pendiente del terreno aumenta.

Gráfica 4. Comportamiento del gasto respecto a los primero 4 puntos procesados en la zona de estudio

Elaboración propia.

Es importante señalar que el proceso de filtración del agua superficial ocasionado por las lluvias que se presentaron en el lapso de un mes se procesó para analizar la escorrentía con un modelo cartográfico en 2d y 3d, además de analizar el perfil para determinar los puntos de mayor riesgo, donde se aprecia una dirección de norte a sur como se muestra en la figura 6.

Figura 6. El modelo cartográfico representa la dirección del agua (pendiente y aguas abajo) con el software libre Surfer, la trayectoria del agua tiene una longitud de 126 metros desde el punto más alto al más bajo y una dirección de norte a sur

Elaboración propia.

En la Figura 7, se apreciar el comportamiento del gasto en forma laminar de dos perfiles (3 y 4), la razón por la que se realizó este análisis fue para tener más información del comportamiento hidrológico dentro del colapso de la estructura.

Figura 7. El modelo estructural de la zona de estudio muestra una profundidad de 55 metros, en una longitud de entre 90 y 40 metros. Para analizar mejor el comportamiento laminar en la zona colapsada

Elaboración propia.

El proceso exógeno generado en el lugar es muy importante, ya que su dinámica atmosférica ha influido en la geomorfología del lugar, es decir el proceso evolutivo es producto no solo de la dinámica interna terrestre sino de la climática y meteorológica también Sinem et al, (2019).

Cálculo de la pendiente y su relación con la aceleración de la gravedad

Para la obtención de todos los datos, fue muy importante obtener las coordenadas de los 18 puntos analizados, así como calcular la pendiente o inclinación del punto superficial al punto más bajo del hundimiento; el procedimiento es dividirlo en cuatro cuadrantes colocados al norte, sur, este y oeste, y luego calcular la pendiente promedio de la estructura analizada.

El modelo o relación matemática para calcular la pendiente (14) es:

(14)

En la Gráfica 5 se presenta un modelo general de la aceleración de la gravedad de la zona analizada con relación a las elevaciones con la que se corroboró el adecuado procesamiento (a mayor altura menor la aceleración gravimétrica y a menor altura mayor la aceleración de la gravedad).

Gráfica 5. Los datos obtenidos corroboran que a mayor altura menor la aceleración gravimétrica y a menor altura mayor la aceleración de la gravedad

Elaboración propia.

Resultados y análisis

El análisis geocientífico del desplome estructural en el Municipio de Juan C. Bonilla en el estado de Puebla, México, pudo establecer la estratigrafía de la zona colapsada, con una planicie de 1000 metros a la redonda como zona de riesgo, referente a la zona derrumbada se determinó un mecanismo dinámico del basamento, con escurrimiento del acuífero promedio a 50 metros de profundidad; la zona de estudio se encuentra dentro de la cuenca del Río Alto Atoyac, con una geodinámica de varias décadas de desarrollo activo del manto acuífero, que con el paso del tiempo y el peso extra superficial y sobre explotación de agua, pudieron contribuir a generar un colapso.

Los datos obtenidos del estudio con modelos geométricos, matemáticos y satelitales son los siguientes: 44.5, 130.1 y 125.3 metros de profundidad, diámetro mayor y diámetro menor, respectivamente.

Los datos del perímetro, área y volumen se obtuvieron tomando un modelo geométrico elipsoidal, cuyos datos son: área superficial de 11,281.46 m2, perímetro de 401.036 m y volumen de 759,270.929 m3.

Se obtuvieron los promedios de: aceleración de la gravedad (9.77900181 m/s2), anomalía de Bouguer (-117.10757041 mgal o 0.00117108 m/s2), anomalía al aire libre (131.68342959 mgal o 0.00131683 m/s2), fuerza del agua (97.78857950 N/m3), gasto (0.00000273575 m3/s), presión atmosférica (773.138 hPa o mbar), pendiente en el terreno de la superficie (-1 m) y pendiente promedio del punto alto al bajo (-0.052631579 m).

Conclusiones

El análisis desde las geociencias del colapso estructural en el Municipio de Juan C. Bonilla en el estado de Puebla, México, pudo identificar de forma precisa la geología, geomorfología, hidrología y desarrollar la cartografía y geoestadística de la zona; se estableció que la estructura derrumbada tenía una morfología plana, y cuyo proceso geodinámico endógeno interrelacionó con procesos de intemperismo y meteorización como el agua, precipitación, viento, temperatura y radiación solar. Se concluye que la región analizada es un receptor de agua sustentada por la fractura de un ducto profundo natural que favoreció al hundimiento, generando oquedades (pequeñas cavernas y ductos hidrológicos de entrada y salida de agua).

La zona podría generar un riesgo importante para la población, ya que el proceso dinámico de la zona puede continuar durante varios meses o años, lo que podría ocasionar un nuevo colapso, siendo importante continuar con observaciones del lugar e instrumentar la zona con sensores de deformación para tener una idea de su proceso.

Este estudio puede complementarse con otros métodos geofísicos como los eléctricos y magnéticos, además de sismicidad inducida para determinar las características de vibración del suelo; hubo limitaciones al no permitir las autoridades volar drones por la cercanía del aeropuerto de Huejotzingo.

El presente estudio pretende contribuir a la prevención y gestión de riegos naturales en el estado, gracias al apoyo del Consejo de Ciencia y Tecnología del Estado de Puebla.

Referencias

Instituto Nacional de Estadística y Geografía, INEGI. (2019). Cálculo de gasto o caudal. http://antares.inegi.org.mx/analisis/red_hidro/PDF/CalculoGasto.pdf

Instituto Nacional de Estadística y Geografía, INEGI. (2020). Mapa digital. http://gaia.inegi.org.mx/mdm6/

Máximo-Romero, P., Ramos-Aguilar, R., Martínez-Ortiz, A. y Hernández-Sánchez, A. (2017). Daños estructurales en Ciudad de México, producidos por los sismos del 28 de julio de 1957 y el 19 de septiembre de 1985. Revista UNIMAR, 35(2), 161-186. https://revistas.umariana.edu.co/index.php/unimar/article/view/1533/pdf

Ramos-Aguilar, R., Rivera-Domínguez, A. y Mayorga-Raposo, R. (2005). Aplicación de imágenes satelitales para determinar el clima y la radiación solar en el estado de Puebla y la geohidrología de su zona norte. Ciencia Ergo Sum, 10(003), 283-294. https://www.redalyc.org/pdf/104/10410306.pdf

Ramos-Aguilar, R., Máximo-Romero, P., Narciso-Hernández, J., Mirón-Morales, M. y Beltrán-Cruz, M. (2012). Estudio geoestadístico para obtener la gravedad local, pendiente y cálculo hidrológico de las barrancas Xaltelulco, Tepeloncocone, Tenepanco, Colorada y Quimichule del volcán Popocatépetl. Boletín de Ciencias de La Tierra, (31), 65–84. https://revistas.unal.edu.co/index.php/rbct/article/view/31254

Ramos-Aguilar, R., Máximo-Romero, P., Soto-Cruz, B., Alcántara-Iniesta, S. y Pérez-Luna, J. (2013). Analysis of the Meteorological Variables for Puebla City 2011-2012 Applying the Modeling Ion-Wavelets in a Hypothetical Manner. Atmospheric and Climate Sciences, 3(4), 450–458. https://doi.org/10.4236/acs.2013.34046

Ramos-Aguilar, R., Máximo-Romero, P., Galindo-López, V., Ávila-Cruz, M., Yáñez-Pérez, G., Márquez-Máximo, R., y Covarrubias-Salvatori, V. (2021). Análisis técnico geoestadístico de los ríos Atoyac y Alseseca que atraviesan la ciudad de Puebla, México. Contactos, Revista de Educación en Ciencias e Ingeniería, (119), 5-22. https://contactos.izt.uam.mx/index.php/contactos/article/view/99

Ramos-Aguilar, R., Máximo-Romero, P., Galindo-López, V., Ávila-Cruz, M., Yáñez-Pérez, G., Márquez-Máximo, R., Hernández-Andrade, D. y Covarrubias-Salvatori, V. (2022). Análisis gravimétrico y cálculo de la fuerza de un líquido en el volcán Malinche, Puebla-Tlaxcala, México. Boletín de Ciencias de la Tierra, (49), 9-14. https://doi.org/10.15446/rbct.n49.91906

SIATL. (2019). Institutito Nacional de Estadística y Geografía. http://antares.inegi.org.mx/analisis/red_hidro/siatl

Sinem, E., Barthelmes, F., Reißland, S., Elger, K., Förste, C., Flechtner, F. y Schuh, H. (2019). ICGEM – 15 years of successful collection and distribution of global gravitational models, associated services, and future plans. Earth System Science Data, 11(2), 647–674. https://doi.org/10.5194/essd-11-647-2019


1 Profesor investigador, Facultad de Ingeniería, Benemérita Universidad Autónoma de Puebla, Puebla, México, rogelio.ramos@correo.buap.mx, https://orcid.org/0000-0001-5660-6608 (autor corresponsal)

2 Profesora investigadora, Facultad de Ingeniería, Benemérita Universidad Autónoma de Puebla, Puebla, México, patricia.maximo@correo.buap.mx, https://orcid.org/0000-0002-2698-842X

3 Profesor investigador, Facultad de Ingeniería, Benemérita Universidad Autónoma de Puebla, Puebla, México, máximo.avila@correo.buap.mx, https://orcid.org/0000-0002-8024-2119

4 Profesor investigador, Facultad de Ingeniería, Benemérita Universidad Autónoma de Puebla, Puebla, México, victor.galindo@correo.buap.mx, https://orcid.org/0000-0003-0112-2432

5 Profesora, Facultad de Ingeniería, Benemérita Universidad Autónoma de Puebla, Puebla, México, dolores.ramireza@correo.buap.mx, https://orcid.org/0009-0004-9875-6443

6 Becaria del Colegio de Ingeniería Topográfica y Geodesia, Benemérita Universidad Autónoma de Puebla, Puebla, México, lorena.cardenas@alumno.buap.mx, https://orcid.org/0009-0006-8762-7683

7 Becario del Colegio de Ingeniería Topográfica y Geodesia, Benemérita Universidad Autónoma de Puebla, Puebla, México, argenis.morales@alumno.buap.mx, https://orcid.org/0009-0006-8302-5108

8 Imágenes Geográficas S.A. de C.V., Guadalajara, México, mayanegui@icloud.com https://orcid.org/0009-0002-0601-2822

9 Director del Consejo de Ciencia y Tecnología del Estado de Puebla, Puebla, México, victoriano.covarrubias@puebla.gob.mx, https://orcid.org/0000-0001-5723-2601

EUNA UNA

Escuela de Ciencias Geográficas
Universidad Nacional, Campus Omar Dengo
Apartado postal: 86-3000. Heredia, Costa Rica
Teléfono: (506) 2562-3283
Correo electrónico revgeo@una.ac.cr