SISTEMA NACIONAL DE PROTECCIÓN CIVIL

COORDINACIÓN NACIONAL DE PROTECCIÓN CIVIL

CENTRO NACIONAL DE PREVENCIÓN DE DESASTRES

Proyecto

ELABORACIÓN DE MAPAS DE PROPENSIÓN A LA INESTABILIDAD DE LADERAS UTILIZANDO EL

MÉTODO DEL TALUD INFINITO

Dirección de Investigación

Subdirección de Dinámica de Suelos y Procesos Gravitacionales

Ing. Leobardo Domínguez Morales

Geog. Antulio Zaragoza Alvares

Pas. en Geog. Humberto Ortiz Pérez

Pas. en Ing. Juan Antonio Hernández Lozano

Marzo 2020
 RESUMEN

La estabilidad o inestabilidad de una ladera o talud se revisa a partir del

conocimiento de las propiedades mecánicas de los suelos y/o de las rocas que la
conforman, de su distribución estratigráfica y de la topografía del terreno, las
cuales, con base en modelos y operaciones matemáticas, proporcionan lo que en
ingeniería se conoce como Factor de Seguridad (FS); el cual representa la alta, baja
o nula propensión del terreno a presentar deslizamientos, por lo cual es utilizado
para representar el grado de estabilidad o inestabilidad de una ladera en la que se
conocen tales propiedades mecánicas y topográficas.

En este proyecto se propone una metodología para elaborar mapas de propensión

a la inestabilidad de laderas, basados en la aplicación del método de talud infinito
en un Sistema de Información Geográfica (SIG), tomando como base el Continuo
Nacional de Geología de la República Mexicana, con datos estadísticos de las
propiedades mecánicas, y el modelo de inclinaciones del terreno (pendientes)
generado a partir del Continuo Nacional de Elevaciones (CEM 3.0), ambos del
Instituto Nacional de Estadística y Geografía (INEGI).
Contenido
1. Introducción ................................................................................................................................. 1
2. Justificación...................................................................................................................................... 3
3. Marco Teórico (fundamentos teóricos) .................................................................................. 4
3.1 Susceptibilidad por inestabilidad de laderas ................................................. 4

3.2 Definición de ladera y/o talud.................................................................................. 4

3.3 Procesos de movimiento .............................................................................................6

3.3.1 Tipos de Movimientos ................................................................................................. 7

4. Métodos de análisis ...................................................................................................................... 8
4.1 Factores condicionantes y desencadenantes................................................9

4.2 El método del Talud Infinito .................................................................................... 15

4.2.1 Equilibrio límite ............................................................................................................. 16

4.3 Factor de seguridad (FS) ........................................................................................... 17

4.4 Propiedades mecánicas de los suelos ............................................................. 18

4.4.1 Cohesión y ángulo de fricción............................................................................. 19

4.4.2 Peso volumétrico ........................................................................................................ 21

4.5 Circulo de Mohor ............................................................................................................ 21

4.6 Envolvente de Falla ...................................................................................................... 22

5. Metodología. .................................................................................................................................. 24
5.1 Trabajo de gabinete ..................................................................................................... 25

5.1.1 Recopilación bibliográficas .................................................................................... 25

5.2 Software ............................................................................................................................... 25

 5.3 Procesamiento de la información ...................................................................... 26

5.4 Elaboración de mapas en el software ArcGis ............................................. 30

5.5 Conversión a raster ....................................................................................................... 34

5.6 Conversión de grados a radianes ........................................................................ 37

5.7 Programación de la fórmula del (FS) en el SIG ArcGIS. ....................... 39

6. Resultados ...................................................................................................................................... 42
7. Conclusiones .................................................................................................................................. 46
8. Trabajo futuro ................................................................................................................................ 47
9. Referencias ..................................................................................................................................... 48

INDICE DE FIGURAS

Figura 1. Representación de un talud Artificial (Corte o relleno), B)

Ladera natural (Suarez 1998). ............................................................................................6

Figura 2. Se muestra un deslizamiento o movimiento, con sus

diversas partes (Suarez 1998). ...........................................................................................6

Figura 3. Deslizamiento en los suelos blandos (Sánchez 2009). ............ 8

Figura 4. Esquema simplificado de las variables que intervienen en

un análisis de estabilidad, donde la superficie de falla corresponde a
un arco de circunferencia, para el suelo puramente cohesivo. .................9

Figura 5. Erosión por afloramiento de agua subterránea (Copy of

Lluvias, Presión de poros y sus efectos, s/f). ........................................................... 11

Figura 6. Flujo de lodo provocado por las lluvias intensas durante el

paso del huracán Mitch en Nicaragua, en el 1998 (Conde, 2002). .......... 12
Figura 7. Terremoto provoco un derrumbe en Indonesia, (News
Mundo2008). ............................................................................................................................. 13

Figura 8. Esta imagen representa un factor condicionante en la

intervención de laderas ocasionadas por actividad antrópica
(Sostenible 2019)...................................................................................................................... 14

Figura 9. El Deslizamiento del Terreno corresponde a un movimiento

del suelo, generalmente por acción de sobrecargas de edificios
(Construpedia S.F). ................................................................................................................. 15

Figura 10. Esquema simplificado en el Método de Talud Infinito

(Jaimes, et. al 2008). .............................................................................................................. 16

Figura 11. Envolvente de falla y circulo de Mohr (Suarez, 1998) ............... 23

Figura 12. Datos de la capa de información del INEGI tomando en

cuenta el espesor de 1m, anexando la cohesión, peso volumétrico y
fricción convertida en radianes con diferentes ángulos de inclinación
aplicando la fórmula del FS. ........................................................................................... 29

Figura 13. Graficas de del factor de seguridad y el ángulo de

inclinación tomando el espesor de 1 metro. ..................................................... 29

Figura 14. Tipos de rocas y suelos con la fórmula del factor de

seguridad con diferentes ángulos de inclinación tomando espesor de
1 m. ................................................................................................................................................... 30

Figura 15. Mapa de cohesión de las rocas del territorio mexicano

integrado dentro de un SIG, tomando los datos de la tabla del
Instituto de Ingeniera. Fuente: Elaboración propia......................................... 31
Figura 16. Mapa de del ángulo de fricción obtenido de la capa de
información del INEGI tomando en cuenta lo valores del Instituto de
Ingeniera de la tabla 1. Fuente: Elaboración propia ....................................... 32

Figura 17. Mapa de peso volumétrico de la geología de las distintas

formaciones rocosas del país. Fuente: Elaboración propia ....................... 33

Figura 18. Mapa de pendientes de la República Mexicana a partir del

Modelo Digital de Elevación. Fuente: Elaboración propia ........................ 34

Figura 19. En la siguiente figura se muestra la capa de litología del

talud infinito cargado el software QGis. Fuente: Elaboración propia 35

Figura 20. Representación para la convención a raster en el software

Qgis. ................................................................................................................................................ 36

Figura 21. Mapa de fricción convertida a raster, donde la cruz que se

muestra dentro del mapa representa a un punto en específico y la
parte circulada de color rojo es el valor del pixel que arroja el punto
seleccionado. ............................................................................................................................ 37

Figura 22. Herramienta de SIG Raster calculator según Osorio 2010

permite crear y ejecutar expresiones de álgebra de mapas en una
herramienta .............................................................................................................................. 38

Figura 23. Cuadro de dialogo de la herramienta calculadora raster del

SIG ArcGIS. .................................................................................................................................. 39

Figura 24. Representación de la ventana con la herramienta “Add

Data “donde se seleccionó la carpeta de pendiente, dentro de la
carpeta se seleccionó la capa de “PenRan_NALCCL.img” ........................ 40
Figura 25. Mapa de pendientes convertido de grados a radianes ........ 41

Figura 26. Ventana de raster calculator donde se ingresó la ecuación

del factor de seguridad (FS). ........................................................................................... 41

Figura 27. Imagen del factor de seguridad (FS) ................................................ 42

Figura 28. Histograma que muestra la variación del FS, donde se

aprecia que los valores que más se repiten son los que están más
cercanos a la unidad............................................................................................................ 43

Figura 29. Comparación del Mapa del (FS) y del Mapa de

Susceptibilidad ..................................................................................................................... 45

ÍNDICE DE ECUACIONES

Ecuación 1 .................................................................................................................................................................. 17

Ecuación 2 ................................................................................................................................... 17

Ecuación 3 ................................................................................................................................... 18

Ecuación 4 ................................................................................................................................... 19

Ecuación 5 ................................................................................................................................... 19

Ecuación 6 ................................................................................................................................... 21

Ecuación 7 .................................................................................................................................. 22

Ecuación 8 .................................................................................................................................. 22

Ecuación 9 .................................................................................................................................. 24

Ecuación 10 ................................................................................................................................ 28
ÍNDICE DE TABLAS

Tabla 1. Valores promedio de propiedades geotécnicas del grupo

litológico del país propuestas (Suarez, 1998). ...................................................... 27

Tabla 2. Ejemplo de cálculos de factor de estabilidad (Burns 2012). 44

 1. Introducción
El problema de inestabilidad de laderas es complejo debido al número de
variables que intervienen en la determinación del grado de estabilidad, o bien, en
su propensión a la inestabilidad. Para el caso de México, el problema cobra gran
relevancia debido a su ubicación geográfica, topografía accidentada, sus
características geotécnicas y geológicas y la existencia de centros de población en
laderas de cerros, barrancas y montañas.
De acuerdo con casos documentados en México, la mayoría de fenómenos de
inestabilidad suceden durante lluvias intensas o de larga duración y en menor
proporción por sismos y/o erupciones volcánicas. Sin embargo, esta situación se ve
aumentada o exacerbada por las condiciones de marginación, pobreza y
crecimiento acelerado, y en ocasiones desordenado, de las ciudades y de las zonas
rurales que se desarrollan en o cerca de zonas montañosas, mismas que
contribuyen enormemente a detonar deslizamientos, caídos, derrumbes o flujos
de lodo y escombro.
Como parte de los objetivos, misión y visión del CENAPRED 2018, para identificar
las zonas de mayor peligro a este tipo de fenómenos, se propuso desarrollar una
metodología para elaborar mapas de susceptibilidad a inestabilidad de laderas
mediante la aplicación del método del talud infinito, el cual es una simplificación
de los métodos cualitativos o determinísticos, mismos que requieren valores de las
propiedades mecánicas de los suelos y/o de las rocas, así como de la topografía del
terreno.
Para el desarrollo de este proyecto se asignaron valores a las propiedades
mecánicas de las distintas litologías que componen en el territorio nacional, según
el Continuo Nacional de Geología de la República Mexicana de 1995, del Instituto
Nacional de Estadística y Geografía (INEGI) y los valores de resistencia sugeridos

1
por el Instituto de Ingeniería de la UNAM, mismos que se integraron a un Sistema
de Información Geográfica (SIG). Con estos datos y el mapa de pendientes
(inclinaciones) generado con el Continuo de Elevaciones Mexicano (CEM 3.0) de 15
m de resolución, también del INEGI, se programó en ArcGis la ecuación del
Método del Talud Infinito, obteniendo así el mapa resultante, a nivel nacional, de
valores del FS, los cuales representan el grado de propensión del terreno a
presentar inestabilidad.
Este trabajo tiene como finalidad identificar las zonas más propensas a
inestabilidad de laderas que sirviera como punto de comparación con el Mapa
Nacional de Susceptibilidad a la Inestabilidad de Laderas (cualitativo), pero en este
caso con un método cuantitativo o determinista, como el método de talud infinito,
que tiene como base la teoría de equilibrio limite. Es decir, que ningún valor del FS
sea menor a la unidad. Para este caso, se requieren valores de las propiedades
mecánicas de suelos, tales como cohesión, ángulo de fricción, peso volumétrico y
la topografía del terreno. Dichas variables son susceptibles de asignar a los
principales tipos de suelos y rocas que componen el territorio nacional y generar
un mapa de pendientes, a través del CEM 3.0. Dichas variables pueden interactuar
programando la fórmula del Talud Infinito en ArcGis y/o Qgis.

2
2. Justificación
A pesar que se han realizado varios estudios de inestabilidad de laderas en nuestro
país es uno de los temas más frecuentes de investigación, en particular en época
de lluvias y cuando ocurren sismos. Sin embargo, las estadísticas de su ocurrencia,
su estudio y la difusión de medidas preventivas son, quizá, las menos desarrolladas
entre los fenómenos perturbadores.
México es un país donde ocurren muchos fenómenos naturales que generan la
inestabilidad de laderas en muchas partes del territorio mexicano, debido a dichos
fenómenos existen remociones de masas en tramos carreteros que ocasionan
bloqueo en las vías de comunicación; otro factor que interviene en los
deslizamientos de masas es la actividad humana tanto en la deforestación de
árboles así como en la minería y obras civiles; generando todo esto un impacto
negativo en el ámbito socio-económico y ambiental ya que dan apertura a
pérdidas de vida (tanto humana, de flora y fauna), así como afectaciones a las
poblaciones dañando sus viviendas e infraestructura.
Debido a esta problemática se realiza una estimación de la susceptibilidad a la
inestabilidad de laderas mediante el método de talud infinito calculando el factor
de seguridad (FS) para generar mapas en el software ArcGis y su integración
futura al Atlas Nacional de Riesgos, que proporcionarán una base de
conocimientos del territorio y de los peligros que pueden afectar a la población, ya
que sirve como herramienta para tener una mejor planeación del desarrollo de
infraestructura más segura y de esta forma contribuir a la toma de decisiones para
la reducción de riesgos de desastres.

3
3. Marco Teórico (fundamentos teóricos)

3.1 Susceptibilidad por inestabilidad de laderas

La susceptibilidad es una propiedad del terreno que indica que tan favorables o
desfavorables son las condiciones de éste, para que pueda ocurrir inestabilidad, y
se refiere solamente a factores intrínsecos (condicionantes) a los materiales
naturales de la ladera, sin considerar factores desencadenantes, como sería el caso
de la precipitación o la sismicidad.
El análisis de susceptibilidad se refiere al proceso de cuantificación relativa de la
importancia que puede tener cada uno de los factores intrínsecos o
condicionantes en la posibilidad que se produzcan fenómenos de inestabilidad de
laderas. Con la suma de los pesos (contribución) de cada uno de los factores
considerados se obtienen diferentes valores de susceptibilidad en una
determinada área del territorio. Estos valores de susceptibilidad se clasifican
generalmente en categorías genéricas: susceptibilidad muy alta, alta, moderada,
baja y muy baja.
La forma más habitual de exponer este análisis es a través de mapas de
susceptibilidad por inestabilidad de laderas. Estos mapas presentan la clasificación
de la estabilidad relativa o la propensión a la inestabilidad de un área en diferentes
categorías, a las que se asigna una escala de color, siguiendo los criterios
habituales para los análisis de peligros y riesgos (rojos y naranjas para las
categorías de susceptibilidad alta y amarillos y verdes para los rangos de
susceptibilidad media y baja), y muestra donde hay o no, condiciones para que
pueda ocurrir inestabilidad de ladera (Gonzales et al., 2005).

3.2 Definición de ladera y/o talud

Un talud o ladera es una masa de tierra que no es plana sino que posee pendiente
o cambios de altura significativos. En la literatura técnica se define como ladera
4
cuando su conformación actual tuvo como origen un proceso natural y talud
cuando se conformó artificialmente (Figura 1).
Las laderas que han permanecido por muchos años pueden fallar en forma
imprevista debido a cambios topográficos, sismicidad, flujos de agua subterránea,
cambios en la resistencia del suelo, meteorización o factores de tipo antrópico o
natural que modifiquen su estado natural de estabilidad (Suarez 1998).

En el talud o ladera se delimitan los siguientes elementos esenciales:

 Altura: Es la distancia vertical entre el pie y la cabeza, la cual se presenta
claramente definida en taludes artificiales pero es complicada de cuantificar
en las laderas debido a que el pie y la cabeza no son accidentes topográficos
bien marcados.
 Pie: Corresponde al sitio de cambio brusco de pendiente en la parte inferior.
 Cabeza o escarpe: Se refiere al sitio de cambio brusco de pendiente en la
parte superior.
 Altura de nivel freático: Distancia vertical desde el pie del talud o ladera
hasta el nivel de agua medida debajo de la cabeza.
 Pendiente: Es la medida de la inclinación del talud o ladera. Puede medirse
en grados, en porcentaje o en relación m/1, en la cual m es la distancia
horizontal que corresponde a una unidad de distancia vertical (Sánchez
2009).

Todos éstos son datos que se pueden determinar en estudios de mecánica de

suelos y levantamientos topográficos. Así mismo, se pueden estimar a partir de
información estadística y/o casos documentados.

5
 Figura 1. Representación de un talud Artificial (Corte o relleno), B) Ladera natural (Suarez 1998).

3.3 Procesos de movimiento

Los movimientos se presentan comúnmente, a lo largo de superficies de falla, por

efecto de la gravedad, en forma de caída libre, movimientos de masa, erosión o
flujos. Al ocurrir la falla de una ladera y dependiendo del tipo de suelos o de rocas
que los conforman, algunos segmentos del material fallado pueden moverse hacia
arriba, mientras otros se mueven hacia abajo (Figura 2). En el siguiente esquema
se indican las principales características que conforman un deslizamiento.

Figura 2. Se muestra un deslizamiento o movimiento, con sus diversas partes (Suarez 1998).
6
3.3.1 Tipos de Movimientos
La inestabilidad de taludes se puede clasificar en base a distintos criterios.
Los tipos de movimientos principales son:
 Desprendimientos
Vuelcos
Reptación
Deslizamientos: Rotacionales y Trasnacionales
Flujos
Los movimientos tienden a comportarse como de caída libre cuando la pendiente
superficial es de más de 75 grados. En taludes de ángulo menor generalmente, los
materiales rebotan y en los taludes de menos de 45 grados los materiales tienden
a rodar.
 Deslizamientos
Son movimientos que se generan al superarse la resistencia al corte del material y
tienen lugar a lo largo de una o varias superficies o a través de una franja
relativamente estrecha del material. Según Suarez (1998) las superficies de
deslizamiento son notorias (figura 3).
Los deslizamientos pueden obedecer a procesos naturales o a desestabilización de
masas de tierra por el efecto de cortes, rellenos, deforestación, etc.

7
 Figura 3. Deslizamiento en los suelos blandos (Sánchez 2009).

 Deslizamientos Rotacionales.
Se presentan a lo largo de una superficie de deslizamiento interna, de modo
circular o cóncava. El movimiento tiene una índole más o menos rotacional
alrededor de un eje dispuesto en forma paralela al talud.

 Deslizamientos Traslacionales.
En este tipo de deslizamiento de traslación el movimiento de la masa tiene un
movimiento hacia fuera o hacia abajo, a lo largo de una superficie más o menos
plana o ligeramente ondulada y tiene muy poco o nada de movimiento de
rotación o volteo.

4. Métodos de análisis
El análisis de la estabilidad de una ladera o de un talud no es una tarea fácil y
requiere el empleo de conocimientos teóricos y experiencia para definir el o los
métodos a utilizar, tanto para el análisis numérico de la estabilidad como para la
selección de las técnicas de exploración y determinación de las propiedades de los

8
suelos. Los métodos de análisis y los cálculos de estabilidad se basan en un
planteamiento físico-matemático, en el que intervienen las fuerzas que se oponen
a la falla (fuerzas resistentes) y las que provocan la falla (fuerzas actuantes) sobre la
ladera o talud, mismas que determina su comportamiento y condiciones de
estabilidad (Figura 4).

Figura 4. Esquema simplificado de las variables que intervienen en un análisis de estabilidad, donde la
superficie de falla corresponde a un arco de circunferencia, para el suelo puramente cohesivo.

Los métodos de análisis se pueden agrupar en: métodos basados en la teoría de

equilibrio límite, que a su vez se basan en la teoría de la plasticidad y los que están
basados en la investigación del estado de esfuerzo.

4.1 Factores condicionantes y desencadenantes

Para este estudio de los movimientos es fundamental recalcar que factores
condicionan la estabilidad de los taludes y aquellos que se presentan como
desencadenantes.
Ejemplos de factores que influyen en la inestabilidad de taludes.
 Factores Naturales.

9
Se trata de agentes que integran la meteorización, siendo los agentes erosivos los
que tienen una mayor incidencia como factores desencadenantes, así como los
fenómenos de carácter tectónico que influyen en la corteza terrestre.
 El Agua.
El agua forma parte de uno de los factores naturales más importantes de mayor
incidencia como factor condicionante y desencadenante en la aparición de
inestabilidades.

Aguas Subterráneas:
Las aguas subterráneas influyen en el interior de la red de fracturación de un
macizo rocoso o de manera intersticial en los suelos. La absorción del agua de
Infiltración produce una continua expansión y contracción de los suelos, en
períodos alternantes de lluvia y sequía, que varía las características resistentes de
los mismos.
Cuando el agua subterránea aflora a la superficie del terreno puede producir el
desprendimiento de las partículas de suelo generando cárcavas. Cuando en el
avance de una dicha cárcava de erosión su superficial ésta captura un contacto
con afloramiento de agua subterránea éste contacto trata de ampliarse (figura 5)
en un proceso de deslizamientos laterales progresivos. Se presentan dos casos
diferentes:
a) En formaciones sedimentarias y aluviales sobre los planos de estratificación.
b) En suelos residuales sobre las superficies de cambio de meteorización.
Los procesos de erosión son similares en suelos de origen sementario o residual.
Estos procesos de erosión pueden ser activados en épocas de lluvia por el agua
infiltrada, la cual produce corrientes de agua subterránea no permanente en los
mantos de suelo de mayor permeabilidad (Suárez, 1998).

10
 Figura 5. Erosión por afloramiento de agua subterránea (Suarez 1998).

 Lluvia
Se presenta como un factor desencadenante de inestabilidades, contribuyendo a
aumentar la acción de diversos factores condicionantes, debido que comienza a
meteorizar el material que presenta la ladera, según Suárez (1998), Cuando el
régimen de lluvias es torrencial, el agua caída se canaliza en zonas deprimidas
causando acciones similares a las corrientes de agua (figura 6).

11
 Figura 6. Flujo de lodo provocado por las lluvias intensas durante el paso del huracán Mitch en
Nicaragua, en el 1998 (Conde, 2002).

 Sismicidad y Volcanismo.

Estos se cuentan como uno de los factores desencadenantes más importantes en

presencia de deslizamientos que pueden ocasionar daños graves. Cuando ocurre
un sismo se generan una serie de vibraciones, que se propagan como ondas de
diferente frecuencia, la aceleración, vertical y horizontal asociada a esas ondas,
origina la intensidad del estado de esfuerzos en el interior del terreno afectando al
equilibrio de los taludes, un claro ejemplo de un sismo ocurrió en indonesia con
una magnitud de 6.4, afecto la isla central de Lombok (figura 7).

12
 Figura 7. Terremoto provoco un derrumbe en Indonesia, (News Mundo 2008). Fuente recuperado de:
https://www.bbc.com/mundo/noticias-internacional-44997228

 Actividad Biológica.
La cubierta vegetal, por medio de sus raíces ayuda a mantener la estabilidad de los
taludes debido a que sirve de unión entre los componentes del suelo. También la
vegetación contribuye al drenaje, absorbiendo parte del agua contenida en el
terreno y atenúa la degradación superficial del mismo. En su aspecto negativo,
existen raíces capaces de disgregar el terreno u desplazar bloques de magnitud
considerable, debido a sus efectos mecánicos en grietas y fisuras, actuando como
cuñas (Suárez, 1998).

 Actividad Humana.
La actividad humana influye mucho en la inestabilidad de laderas debido que
modifica el equilibrio que existe en la naturaleza (deforestaciones) (figura 8). Dos
de las actividades más importantes son la Minería y las Obras Civiles.

13
 Figura 8. Esta imagen representa un factor condicionante en la intervención de laderas ocasionadas por
actividad antrópica (Sostenible 2019). Fuente recuperado de: https://sostenibilidad.semana.com/medio-
ambiente/articulo/15-millones-de-arboles-fueron-talados-en-remedios-antioquia/45377

Excavaciones:
Las excavaciones forman parte de uno de los factores desencadenantes más
extendidos debido a la necesidad de las mismas en las obras civiles como vialidad,
túneles. Producen una variación del estado de equilibrio del terreno, que se
traducen en planos de potenciales deslizamientos, desequilibrio de masas, etc.
Sobrecargas:
Es el resultado del incremento de peso debido a diversos tipos de construcciones
sobre el terreno natural, como se muestra en la figura 9. También puede
producirse por el peso del agua infiltrada en el terreno, como consecuencia de
fugas en conducciones de alcantarillado, canales, etc (Sánchez, 2009).

14
 Figura 9. El Deslizamiento del Terreno corresponde a un movimiento del suelo, generalmente por acción de
sobrecargas de edificios (Construpedia S.F).Fuente recuperado de:
https://www.construmatica.com/construpedia/Deslizamiento_del_Terreno

4.2 El método del Talud Infinito

Con frecuencia, en los deslizamientos de gran magnitud, la mayor parte de la

masa deslizada se mueve aproximadamente en forma paralela a la superficie del
terreno. La naturaleza del movimiento está controlada por algún elemento
geológico como una capa de roca o una capa de materiales poco resistentes. Si la
longitud relativa del deslizamiento es muy grande en relación con su espesor, la
contribución de la resistencia en la cabeza y el pie del deslizamiento, es menor
comparada con la resistencia del resto de la superficie de falla.
Se muestra una falla paralela a la superficie del talud, a una profundidad somera y
la longitud de la falla es mayor comparada con su espesor. Este tipo de
deslizamiento se puede analizar suponiendo un talud infinito (figura 10).
El principal uso del método del talud infinito es la elaboración de planos de
amenaza a los deslizamientos mediante el uso de SIGs.
Para un talud uniforme y relativamente largo, en el cual el mecanismo de falla
esperado no es muy profundo, los efectos de borde son despreciables y el factor de
15
seguridad puede calcularse (para un talud infinito) a partir de una unidad de área
con base en el criterio Mohr Coulomb (Suárez, 1998).

Figura 10. Esquema simplificado en el Método de Talud Infinito (Jaimes, et. al

2008).

4.2.1 Equilibrio límite

En el análisis de estabilidad de taludes, el método de equilibrio límite se usa

generalmente para calcular el factor de seguridad de taludes basado en el criterio
de Mohr-Coulomb. Sin embargo, el criterio de Mohr-Coulomb se limita a la
descripción de la masa rocosa. Para superar sus deficiencias, este documento
combinó el criterio de Hoek-Brown y el método de equilibrio límite y propuso una
ecuación para calcular el factor de seguridad de la pendiente con el método de
equilibrio límite en el criterio de Hoek-Brown a través de una resistencia cohesiva
equivalente y el ángulo de fricción
Este análisis de estabilidad consiste en determinar si existe suficiente resistencia
en los suelos del talud para soportar los esfuerzos de cortante que tienden a
causar la falla o deslizamiento.

16
La mayoría de los métodos de límite de equilibrio tienen en común, la
comparación de las fuerzas o momentos resistentes y actuantes sobre una
determinada superficie de falla. Las variaciones principales de los diversos
métodos son, el tipo de superficie de falla y la forma cómo actúan internamente
las fuerzas sobre la superficie de falla.

4.3 Factor de seguridad (FS)

El factor de seguridad es empleado para conocer cuál es el factor de amenaza

para que el talud falle en las peores condiciones de comportamiento para el cual
se diseña. Fellenius (1922) presentó el factor de seguridad como la relación entre la
resistencia al corte real, calculada del material en el talud y los esfuerzos de corte
críticos que tratan de producir la falla, a lo largo de una superficie supuesta de
posible falla (ecuación 1).

Ecuación 1

En las superficies circulares donde existe un centro de giro y momentos

resistentes y actuantes (ecuación 2).

Ecuación 2

La mayoría de los sistemas de análisis asumen un criterio de límite de equilibrio

donde el criterio de falla de Coulomb es satisfecho a lo largo de una determinada
superficie. Se estudia un cuerpo libre en equilibrio, partiendo de las fuerzas
actuantes y de las fuerzas resistentes que se requieren para producir el equilibrio.

17
Calculada esta fuerza resistente, se compara con la disponible del suelo o roca y se
obtiene una indicación del factor de seguridad (Suárez, 1998).

Para calcular el factor de seguridad de un talud infinito se utiliza la siguiente

ecuación 3:

Ecuación 3

Dónde:
c’= Cohesión (kg/m2)
=Peso Volumétrico (kg/m3)
’ =Ángulo de Fricción
=Angulo idealizado variado
= Espesor
Esta ecuación se divide en tres términos, el primero es la componente cohesiva de
la resistencia, el segundo es la componente friccionante y el tercer término toma
en cuenta la disminución de la resistencia debido a la presión de poro (Luzi et al.,
2000).

4.4 Propiedades mecánicas de los suelos

La representación matemática del fenómeno de falla al cortante en un

deslizamiento se realiza utilizando las teorías de la resistencia de materiales. Las
rocas y los suelos al fallar al corte se comportan de acuerdo a las teorías
tradicionales de fricción y cohesión, según la ecuación generalizada de Coulomb
(ecuación 4).
18
  c´+(Tan ’ (Para suelos saturados)
Ecuación 4

 c´+(Tan ’ + (a) Tan ’ (Para suelos parcialmente saturados)

Ecuación 5

Dónde:
 = Esfuerzo de resistencia al corte
c´ = Cohesión o cementación efectiva
 = Esfuerzo normal total
= Presión del agua intersticial o de poros 
a = Presión del aire intersticial
’= Ángulo de fricción interna del material
”= Ángulo de fricción del material no saturado

El análisis de la ecuación de Coulomb requiere predefinir los parámetros, ángulo

de fricción y cohesión, los cuales se consideran como propiedades intrínsecas del
suelo. La presencia del agua reduce el valor de la resistencia del suelo
dependiendo de las presiones internas o de poros de acuerdo a la ecuación de
Coulomb, en la cual el factor u está restando al valor de la presión normal. La
presión resultante se le conoce con el nombre de presión efectiva (Suárez 1998).

4.4.1 Cohesión y ángulo de fricción

La cohesión es una medida de la cementación o adherencia entre las partículas de
suelo. La cohesión en mecánica de suelos es utilizada para representar la
resistencia al cortante producida por la cementación, mientras que en la física este

19
término se utiliza para representar la tensión. En suelos eminentemente
granulares en los cuales no existe ningún tipo de cementante o material que
pueda producir adherencia, la cohesión se supone igual a 0 y a estos suelos se les
denomina Suelos no Cohesivos (Suárez, 1998).
Es el mayor o menor grado de resistencia que ofrece un mineral a la rotura,
deformación, aplastamiento, curvatura o pulverización. Se distinguen las
siguientes clases de tenacidad:
Frágil: es el mineral que se rompe o pulveriza con facilidad. Ejemplos: cuarzo y el
azufre.
Maleable: el que puede ser batido y extendido en láminas o planchas. Ejemplos:
oro, plata, platino, cobre, estaño.
Dúctil: el que puede ser reducido a hilos o alambres delgados. Ejemplos: oro, plata
y cobre.
Flexible: si se dobla fácilmente pero, una vez deja de recibir presión, no es capaz de
recobrar su forma original. Ejemplos: yeso y talco.
Elástico: el que puede ser doblado y, una vez deja de recibir presión, recupera su
forma original. Ejemplo: la mica.

El ángulo de fricción es la representación matemática del coeficiente de

rozamiento, el cual es un concepto básico de la física:
Coeficiente de rozamiento = Tan 
El ángulo de fricción depende de varios factores, entre ellos algunos de los más
importantes son:
 Tamaño de los granos
 Forma de los granos
 Distribución de los tamaños de granos
 Densidad

20
La fricción interna de un suelo, está definida por el ángulo cuya tangente es la
relación entre la fuerza que resiste el deslizamiento, a lo largo de un plano, y la
fuerza normal “p” aplicada a dicho plano. Los valores de este ángulo llamada
ángulo de fricción interna f, varían de prácticamente 0º para arcillas plásticas, cuya
consistencia este próxima a su límite líquido, hasta 45º o más, para gravas y arenas
secas, compactas y de partículas angulares. Generalmente, el ángulo f para arenas
es alrededor de 30º.

4.4.2 Peso volumétrico

O peso específico de la masa del suelo a relación entre el peso del mismo y el
volumen que ocupa, está dado por la ecuación 6.

Ecuación 6

Dónde:
Wm = La masa de suelo.
Vm= Volumen total de la muestra de suelos.

4.5 Circulo de Mohor

En un análisis en dos dimensiones, los esfuerzos en un punto pueden ser
representados por un elemento infinitamente pequeño sometido a los esfuerzos
, , y . Si estos esfuerzos se dibujan en unas coordenadas -, se puede trazar
el círculo de Esfuerzos de Mohr. En este círculo se definen los valores de σ máximo
(1) y σ mínimo (), conocidos como Esfuerzos principales. Para interpretar
correctamente el fenómeno de falla al cortante en un talud debe tenerse en
cuenta cuál es la dirección de los esfuerzos principales en cada sitio de la
superficie de falla. El esfuerzo σ1 es vertical en la parte superior de la falla y
horizontal en la parte inferior (Suárez, 1998).

21
4.6 Envolvente de Falla

El círculo de Mohr se utiliza para representar o describir la resistencia al cortante

de los suelos, utilizando la envolvente de falla Mohr – Coulomb, lo cual equivale a
que una combinación crítica de esfuerzos se ha alcanzado. Los esfuerzos por
encima de la envolvente de falla no pueden existir. La envolvente de falla Mohr -
Coulomb es generalmente una línea curva que puede representarse en la
ecuación 7.

S= A ( ‘) b
Ecuación 7

Dónde:
s = Resistencia al cortante
’= Esfuerzo normal efectivo
A y b = Constantes
Esta curva se define como una recta aproximada dentro de un rango seleccionado
de esfuerzos (figura 11).
s=c’+tan ’
Ecuación 8

22
 Figura 11. Envolvente de falla y circulo de Mohr (Suarez, 1998)

Dónde:
c’= Intercepto del eje de resistencia (cohesión) y
’= Pendiente de la envolvente (ángulo de fricción).

En la mayoría de los suelos, la envolvente de falla para niveles de esfuerzos

pequeños no es recta sino curva y el error de asumirla como recta puede modificar
sustancialmente los resultados de un análisis. En la realidad, no existe un ángulo
de fricción para esfuerzos normales bajos y es preferible utilizar todos los valores
de la envolvente. Sin embargo, los ensayos normales de resistencia al cortante no
se realizan con suficientes puntos para determinar las características de la curva
en el rango de esfuerzos bajos. Hawkins (1996) indica que es recomendable
presentar los ángulos de fricción como una función de las presiones normales
(ecuación 9).

23
 ’= f(’)
Ecuación 9

Y ’ (ultimo) = pendiente de la parte recta de la envolvente el circulo de Mohr

puede extenderse también al análisis de suelos parcialmente saturados, teniendo
en cuenta las presiones en el agua y el aire en los vacíos (Fredlund 1978).

5. Metodología.

24
5.1 Trabajo de gabinete
Esta etapa se realizó mediante una serie de actividades las cuales estuvieron
dirigidas a la obtención y recopilación de toda la información referente al área de
estudio, como los diferentes programas que se implementaron.

5.1.1 Recopilación bibliográficas

Se recabaron artículos técnicos, capítulos de libros y presentaciones relacionados
con los métodos para calcular el Factor de Seguridad (FS). Particularmente,
aquellos métodos que hacen referencia al método del Talud Infinito. Se trata de
uno de los métodos más simples que existen para calcular el FS en laderas
potencialmente inestables.
Se recabó el archivo Shape del Mapa Nacional de Geología, al cual se le asignaron
los campos con los atributos del ángulo de fricción de los suelos, los valores de la
cohesión y el peso volumétrico de los tipos de roca, tomados de las propuestas
realizadas por los investigadores del Instituto de Ingeniería de la UNAM. Dichas
modificaciones, permitirán implementar las operaciones que se programarán en
el software del sistema de información geográfica del CENAPRED, aplicando el
método del Talud Infinito.

5.2 Software
Se utilizó el software como ArcGis y QGIS para la elaboración de capas o polígonos
como también para realizar la conversión de vectorial a raster utilizando
información recopilada durante la investigación.
ArcGis: Es un sistema de información geográfica, (SIG) posee una interface gráfica
de usuario, de manejo sencillo, posibilita cargar con facilidad datos espaciales y
tabúlales para ser mostrados en forma de mapas, tablas y gráficos, también provee
las herramientas necesarias para consultar y analizar los datos y presentar los
25
resultados como mapas con calidad de presentación y se puede extender con
software adicional para dispositivos Windows (Von, 2006).
Q-GIS: Es un Sistema de Información Geográfica (SIG) de Código Abierto licenciado
bajo GNU - General Public License. QGIS es un proyecto oficial de Open Source
Geospatial Foundation (OSGeo). Corre sobre Linux, Unix, Mac OSX, Windows y
Android y soporta numerosos formatos y funcionalidades de datos vector, datos
ráster y bases de datos (Felpeto et al., 2001).

5.3 Procesamiento de la información

La ecuación del (FS) depende de varios parámetros del suelo, como la cohesión, el
ángulo de fricción y el Peso volumétrico que están relacionados con cada tipo
específico de roca / suelo.
Este modelo hace grandes suposiciones, pero es un enfoque más analítico y
debería mejorar el análisis de riesgos en comparación con la asignación arbitraria
de pesos a la pendiente y al tipo de suelo
Con la capa de información que proporcionó el Instituto Nacional de Estadística y
Geografía (INEGI), a la cual se le asignaron valores promedio de propiedades
mecánicas de suelo representativos de cada grupo litológico usados para los fines
de este estudio (peso específico cohesión c’ , y ángulo de fricción, ´) estos
valores fueron creados por diferentes autores ( Barton, 1974, Hoek y Bray, 1981,
Suárez, 1998 Jibson et al., 2000, González et al., 2002) con la finalidad de asignarle
propiedades a los materiales del territorio mexicano (tabla 1) (Suarez, 1998).

26
Tabla 1. Valores promedio de propiedades geotécnicas del grupo litológico del país propuestas (Suarez, 1998).

Se elaboraron graficas en las que se observa la variación de los valores del FS en

función del ángulo de fricción y de la cohesión para diferentes profundidades (1m,
5m, 10m, 20m, 45m, 85m); así como la influencia del ángulo e inclinación del
terreno que fueron de 10°,15°,25°,35°,45°,55°,65°,75°,85°,88°,89°.
Se convirtieron el ángulo de fricción y el ángulo variado ideal de la ladera en
Radianes utilizando la ecuación 10.

27
 
Ecuación 10

Se emplearon la ecuación del Factor de Seguridad (FS) para cada uno de los tipos
de rocas, posteriormente se realizan las gráficas de acuerdo a dichos tipo de
litología (Sedimentarias, Ígneas, Metamórficas), únicamente tomamos en cuenta
las dos primeras partes de la ecuación 3, la primera parte fue la componente
cohesiva de la Resistencia y la segunda la componente friccionante, omitiendo la
tercer parte debido que no se tuvieron datos de suelos saturados.
Con estos datos se pudo obtener un mejor panorama de las variaciones del Factor
de Seguridad (FS), además de analizar el efecto de las variaciones del ángulo de
inclinación del terreno, el cual está determinado por los diferentes tipos de roca
que componen el territorio nacional, a fin de conocer su influencia en el grado de
estabilidad global.
Con base en el análisis anterior, se determinó que el valor unitario hace que la
fórmula del Talud Infinito converja mejor con la profundidad de 1 metro tomando
en cuenta la roca más susceptible, ya que con éste, los valores del FS se acercaban
más a la unidad (figura 12).
Antes de exportar los valores de las variables geotécnicas y geológicas del
subsuelo, se asignaron valores a los campos que no tenían valor por tratarse de
cuerpos de agua. Con el propósito de identificarlos claramente durante los análisis,
se les asignaron valores para que la formula convergiera.

28
 Figura 12. Datos de la capa de información del INEGI tomando en cuenta el espesor de 1m, anexando la
cohesión, peso volumétrico y fricción convertida en radianes con diferentes ángulos de inclinación aplicando
la fórmula del FS.

Con los datos de la figura 11 se realizaron graficas tomando en cuenta el factor de

seguridad y el ángulo de inclinación para verificar los valores que se acercaban
más a la unidad como se muestra en la figura 13 y 14.

Figura 13. Graficas de del factor de seguridad y el ángulo de inclinación tomando el espesor de 1 metro.

29
Figura 14. Tipos de rocas y suelos con la fórmula del factor de seguridad con diferentes ángulos de inclinación
tomando espesor de 1 m.

Respecto al ángulo de inclinación se obtuvieron valores muy altos en el factor de

seguridad, motivo por el cual se examinaron los valores máximo y mínimo, los
cuales fueron 25.67348253 y 1.0001523, respectivamente.

5.4 Elaboración de mapas en el software ArcGis

La información de referencia fue retomada de los valores propuestos por el
Instituto de Ingeniería de la UNAM, con los cuales se tuvo una referencia general
para asignar propiedades a los materiales; así como de estudios de mecánica de
suelos disponibles en México, también retomados de los trabajos realizados por
dicho instituto.

30
De acuerdo con la teoría de equilibrio límite y las hipótesis del Método de Talud
Infinito, se estableció que ningún valor resultante, una vez aplicado el método,
debería arrojar valores del FS menores a la unidad. Con este propósito, en tabla de
litología del INEGI, a los datos que daban menor a la unidad, se les asigno valores
mínimos de cohesión (basándose en la tabla del Instituto de Ingeniería) para
obligar que el FS no fuera menor a la unidad. Se tomaron valores mínimos y
valores máximos. Dado que el agua no tiene valores de resistencia al esfuerzo
cortante, se les asignó un valor máximo a fin de que en el resultado pudiera
identificarse su presencia y posteriormente colocar una máscara incluir a todos los
cuerpos de agua como de muy baja susceptibilidad.
Teniendo todos los campos con valores, de la tabla de Excel, se convirtió el archivo
a formato DBF, y se ejecutó el comando Join (unión), a fin de agrupar los datos de
geología con propiedades mecánicas de los suelos. Con el propósito de poder
visualizar los resultados en ArcGis. En la figura 15 se muestra el mapa resultante de
Cohesión.

Figura 15. Mapa de cohesión de las rocas del territorio mexicano integrado dentro de un SIG, tomando los
datos de la tabla del Instituto de Ingeniera. Fuente: Elaboración propia
31
Par realizar el mapa de fricción se llevó a cabo un proceso similar al anterior se
realizó para obtener el mapa nacional de ángulo de fricción de los suelos que
componen el territorio nacional cargando la capa en el software ArcGis, (figura 16).

Figura 16. Mapa de del ángulo de fricción obtenido de la capa de información del INEGI tomando en cuenta lo
valores del Instituto de Ingeniera de la tabla 1. Fuente: Elaboración propia

Teniendo estos mapas se procedió a realizar el mapa de peso volumétrico en el

software en ArcGis, basándose en la tabla de Valores promedio de propiedades
geotécnicas propuestas para la estimación de deslizamiento de laderas en función
del grupo litológico del territorio mexicano (Figura 17).

32
 Figura 17. Mapa de peso volumétrico de la geología de las distintas formaciones rocosas del país. Fuente:
Elaboración propia

Para la elaboración del mapa de pendientes es de suma importancia considerar la

resolución del MDE. A partir de éste, con resolución de 15 metros, es decir que cada
píxel representa una distancia de 15 metros de lado, se puede definir como la
relación que existe entre el desnivel y la distancia horizontal entre un punto y otro,
se expresa en porcentaje o en grados, el cual la unidad de medida que se utilizo es
en grados. El mapa de pendientes representa por medio de una gama de colores,
diferentes rangos de pendiente en una zona (figura 18).

33
 Figura 18. Mapa de pendientes de la República Mexicana a partir del Modelo Digital de Elevación. Fuente:
Elaboración propia

5.5 Conversión a raster

Teniendo los mapas anteriores se procedió a realizar la conversión raster dentro
del programa QGis se utilizó la herramienta raster seleccionando conversión, en la
opción rasterizar (vectorial a raster) (figura 19).

34
 Figura 19. En la siguiente figura se muestra la capa de litología del talud infinito cargado el software QGis.
Fuente: Elaboración propia

Dentro de la siguiente ventana en el primer campo se carga la capa de

“Litología_base_Talud_Infinito”, en el siguiente se tomara el campo de “fricción” ,
no se agrega ningún valor fijo se queda en 0.000000, posteriormente se asigna la
unidad de tamaño de raster de salida como “pixel” se procede agregar una
resolución del ancho, alto, horizontal y vertical dándole un valor de 250 000 000
000,de igual forma se le dio una extensión de salida xmin, xmax, ymin, ymax, estos
valores lora roja por de faul al seleccionar la capa de fricción, con la finalidad de
que SIG realice un recorte en todo el contorno del mapa, no se asigna ningún
valor especificado a las bandas de salida y después se seleccionó la opción ejecutar
para que comience a realizare el proceso (figura 20).

35
 Figura 20. Representación para la convención a raster en el software Qgis.

Con la finalidad de que cada celda (también conocida como píxel) tenga un valor
distinto, si se selecciona un punto en específico dentro del mapa arrojara un valor
(Figura 21) y esto sirva para obtener un dato espacial mucho más adecuado para
el análisis, de la misma manera se realizó la conversión a raster del mapa de
cohesión y peso volumétrico.

36
Figura 21. Mapa de fricción convertida a raster, donde la cruz que se muestra dentro del mapa representa a un
punto en específico y la parte circulada de color rojo es el valor del pixel que arroja el punto seleccionado.

5.6 Conversión de grados a radianes

Para seguir con el procedimiento es necesario realizar la conversión del mapa de

pendientes de grados a radianes , utilizando la ecuación 10, para esta conversión
se utilizo el SIG Arcgis, Dentro de la Herramienta ArcToolbox se selecciono la
opción “Raster calulator” (figura 22).

37
 Figura 22. Herramienta de SIG Raster calculator según Osorio 2010 permite crear y ejecutar expresiones de
álgebra de mapas en una herramienta

Después de haber seleccionado dicha opción de Raster calculator, se desplaza un

cuadro donde trae varias funciones y pide agregar capas y variables, se asignó en
el cuadro de dialogo la capa de “mapa de pendientes nacional” y en la parte de
botones de operación se escribió la operación (“Mapa de pendiente NACIONAL.tif”
*3.141592653)/180, tomando en cuenta la ecuación 10, se procede a darle un raster
de salida (output raster) (figura 23).

38
 Figura 23. Cuadro de dialogo de la herramienta calculadora raster del SIG ArcGIS.

5.7 Programación de la fórmula del (FS) en el SIG ArcGIS.

Teniendo todos los mapas convertidos a raster y el mapa de pendientes
convertido de grados a radianes se procedió a programar la fórmula del (FS), para
eso se utilizó el SIG ArcGIS, se cargó la capa en la herramienta “Add Data” al
seleccionarla apareció una ventana donde se cargaron las capas se seleccionó la
capa de pendientes (figura 24).

39
 Figura 24. Representación de la ventana con la herramienta “Add Data “donde se seleccionó la carpeta de
pendiente, dentro de la carpeta se seleccionó la capa de “PenRan_NALCCL.img”

Al seleccionar dicha capa de “PenRan_NACCL.img” se cargó dentro del SIG el

mapa de pendientes en una forma cónica debido que es una de las más indicadas
para las latitudes medias, basada normalmente en dos paralelos estándar, es decir,
se trata de una proyección secante, El espaciado entre latitudes aumenta fuera de
los paralelos estándar (figura 25).

40
 Figura 25. Mapa de pendientes convertido de grados a radianes

Teniendo la capa cargada se procedió a programar la ecuación 3 para eso se

seleccionó la herramienta “calculadora raster” en “ArcToolbox” > “Spatial Analyst
Tools” > “Map Algebra” > “Raster Calculator” abrirá una ventana donde
ingresaremos dicha ecuación (figura 26), agregando una ruta de salida y un
nombre.

Figura 26. Ventana de raster calculator donde se ingresó la ecuación del factor de seguridad (FS).

41
6. Resultados
Una vez programada y aplicada la fórmula del Talud Infinito, con las capas de
información de topografía (pendientes) y geología (cohesión, fricción y peso
volumétrico), se obtuvo el resultado que se muestra en la figura 27. En éste se
muestra la variación a nivel nacional de los valores del FS calculados con el Método
del Talud Infinito. Destaca el gran número de polígonos que se apegan a los
colores rojo y anaranjado, que corresponden a las zonas donde el valor del FS es
muy cercano a uno. Es decir que, de acuerdo con la definición del FS, se trata de
laderas en las cuales el grado de estabilidad es bajo o muy bajo, por lo que pueden
considerarse como de Alta y Muy Alta propensión a deslizamientos.

Figura 27. Imagen del factor de seguridad (FS). Fuente. Elaboración propia.

42
Para conocer la variación de los valores del FS, se generó el histograma que se
muestra en la figura 28, en el cual se puede observar que el valor máximo del FS es
de 397.57, el cual corresponde a los valores verdes y claros del mapa nacional de
factores de seguridad. Así mismo, se aprecia que los valores que más se repiten
son los que se acercan a la unidad y dicha repetición disminuyen gradualmente
hasta los valores más altos del FS, que son los que menos abundan.

Figura 28. Histograma que muestra la variación del FS, donde se aprecia que los valores que más se repiten
son los que están más cercanos a la unidad. Fuente. Elaboración propia.

En un análisis previo, realizado en Excel con los valores de las propiedades

mecánicas (cohesión, fricción y peso volumétrico utilizados en el mapa) para
diferentes grados de inclinación, se pudo observar que esta tendencia era similar,
por lo que los resultados obtenidos con este método son confiables. No obstante,
por las limitaciones de tiempo y de las características de los equipos de cómputo
con los que se cuenta en CENAPRED, no fue posible, al menos para esta etapa del
proyecto, hacer más iteraciones, cambiando otras variables como la profundidad,
que también juega un papel relevante en los cálculos del FS.

43
Finalmente, para observar las diferencias entre esta versión y el Mapa Nacional de
Susceptibilidad a la Inestabilidad de Laderas, disponible en el portal del Atlas
Nacional de Riesgos, se muestra la siguiente figura en la que se comparan ambos
mapas en una escala a nivel nacional (Fig. 29). A primera vista, ambos mapas son
muy similares y tienen prácticamente las mismas tendencias, sin embargo, en el
del Talud Infinito los colores naranja y rojo son más abundantes y extensos que los
de Susceptibilidad.
Al revisar con mayor detalle los valores del mapa de Talud Infinito, se pudo
constatar que si se utilizan los intervalos del FS de la Tabla 2 que proponen Burns y
coautores (2012), el mapa puede cambiar drásticamente, ya que, si bien los valores
más frecuentes del mapa están cercanos a la unidad, también es cierto que
rebasan el valor de 1.5, motivo por el cual, si se aplica dicho criterio, el mapa
tendería hacia los colores más claros. Esta posibilidad se revisará en la segunda
etapa de este proyecto.

Tabla 2. Ejemplo de cálculos de factor de estabilidad (Burns 2012).

44
 Mapa Nacional de Susceptibilidad a la Inestabilidad de Laderas
Método Cualitativo a partir de datos documentados

Figura 29. Comparación del Mapa del (FS) y del Mapa de Susceptibilidad

45
7. Conclusiones
De acuerdo con los resultados obtenidos en esta etapa del proyecto, se pudo
comprobar que es factible elaborar mapas, a nivel nacional, utilizando el la fórmula
del Talud Infinito que es considerado como uno de los métodos más simples para
representar el grado de estabilidad de una ladera o talud, ya que representa
adecuadamente la posibilidad de falla superficial en laderas. Sin embargo, la
certidumbre del método y, consecuentemente, de los mapas, depende en gran
medida de los parámetros de resistencia de los suelos, mismos que, para futuras
versiones, deben verificarse con más datos obtenidos de pruebas de laboratorio
y/o de campo.
Por otra parte, un aspecto muy importante a considerar en la aplicación espacial
de éste método con SIG´s, es la capacidad de procesamiento de los equipos de
cómputo que se utilicen, ya que para análisis a nivel nacional, se tienen serias
dificultades en el manejo de la información por el tamaño de los archivos que se
generan, especialmente de los modelos digitales de elevación, los cuales forman
parte de los insumos para obtener la inclinación del terreno.
Desde el punto de vista ingenieril, el método del talud infinito es viable y
representa adecuadamente muchos de los casos de inestabilidad de laderas que
suceden frecuentemente en el territorio nacional, por lo tanto, su utilización es
recomendada, especialmente en regiones donde se tengan datos más precisos de
las propiedades mecánicas de los suelos y/o de las rocas que componen las laderas
o taludes.
Cabe resaltar que el caso estudiado en este proyecto, sólo consideró dos de los tres
componentes que componen la fórmula. Por una parte ello se debió a que en esta
etapa sólo se consideraron las componentes de cohesión y fricción, despreciando
el efecto de la humedad del suelo, ya que, para fines de susceptibilidad, son las
variables propias del terreno las que definen la propensión a los deslizamientos.

46
Para fines de la estimación del peligro, la humedad debida a fugas o lluvia, es el
factor detonante y puede ser incluida en futuros análisis, siempre y cuando se
tome en cuenta que esto puede requerir una mayor capacidad de cómputo, para
lo cual, los equipos que se tenían disponibles dificultaban este proceso.
Finalmente, la comparación del mapa nacional de susceptibilidad obtenido con el
método del Talud Infinito, sí guarda relación con las zonas de mayor propensión a
la inestabilidad de laderas del Mapa Nacional de Susceptibilidad, elaborado con
métodos cualitativos. Por lo tanto, puede ser utilizado para dar certidumbre a los
métodos menos sofisticados, como los no deterministas o heurísticos.

8. Trabajo futuro
Para las siguientes etapas del proyecto, se deberán hacer iteraciones con valores
diferentes de la profundidad de la superficie de falla, a fin determinar y representar
los casos más frecuentes de deslizamientos que se suceden en México. Por
cuestiones de tiempo y limitaciones técnicas de los equipos de cómputo, no fue
posible elaborar, en esta etapa, mapas de FS para valores diferentes a la unidad, no
obstante, se hicieron cálculos en Excel, donde se pudo verificar que a mayor
profundidad de la superficie de falla, la susceptibilidad o propensión del terreno a
presentar deslizamientos es mayor. Esta es una variable que se deberá incluir para
la obtención de mapas de FS en futuras etapas.
El tema de la humedad del suelo es otra variable que se deberá considerar en
futuros análisis, a fin de tomar en consideración este importante factor. Si bien la
fórmula del talud infinito toma en cuenta esta variable, habrá que hacer
suposiciones de su variación con el tiempo. Es decir, en caso de lluvias o fugas el
valor del FS adquirirá valores cambiantes con el tiempo, por lo tanto, se podrían
elaborar mapas dinámicos, en los que se observe la modificación del FS en función
de la distribución de las lluvias. Esta observación ya había sido considerada por el
Dr. Terzaghi en su célebre artículo de 1950: Mechanism of Landslides.
47
9. Referencias
Alcántara, Ayala, I. (2008). On the historical account of disastrous landslides in
Mexico: the challenge of risk management and disaster prevention.

Alcantara, I. et, al,. (2014). Fasiculo de Inestabilidad de Laderas. 2ª ed. CENAPRED,

México pp. 4 - 8.
Barton, R. (1974), “A review of the shear strength of filled discontinuities in rock”,
Norwegian Geotech. Inst. Publ. No. 105. Oslo: Norwegian Geotech. Inst.

CENAPRED. (2016). Mapa Nacional de Suceptibilidad por Inestabilidad de Laderas.

Obtenido de www.atlasnacionalderiesgos.gob.mx.

CENAPRED. (08 de Mayo de 2018). MANUAL DE ORGANIZACIÓN ESPECÍFICO,

Versión observada por la UGAJ-SEGOB .p.25.

Fellenius W. (1922). “Statens Jarnjvagars Geoteknniska Commission, Stockholm,

Sweden”

Felpeto, A., Araña, V., Ortiz, R., Astiz, M., García, A., 2001. Assessment and modelling
of lava flow hazard on Lanzarote.

Fredlund. D.G., Morgenstern N., , Widger A, (1978), “Shear strength of unsaturated

soils”, Canadian geotechnical Journal, , pp. 311.

González, L., Ferrer, M., Ortuño, L. y Oteo, C. (2002), “Ingeniería geológica”, Prentice
Hall, Pearson Education, Madrid, 715.

48
Gonzáles, Huesca, A., Dominguez, Morales, L., & Castañeda, Martinez, A. (2015).
Generacion de Insumos para el Atlas Ncional de Riesgos. Proyecto,
CENAPRED, MÉXICO, p. 4.

Hoek, E. y Bray, J. W. (1981), “Rock slope engineering”, Institution of Mining and

Metallurgy, 3a edición, 358 p.

Jibson, R.W., Harp, E.L. y Michael, J.M. (2000), “A method for producing digital
probabilistic seismic landslide hazard maps”, Engineering Geology, Vol. 58,
No. 3-4, pp. 271-289.

Jaimes, M. A., Niño, M., & Reinoso, E. (27 de 08 de 2013). Una aproximación para la
obtención de mapas de desplazamiento traslacional de laderas a nivel
regional inducido por sismos. Revista de Ingeniería Sísmica(89), 24.
Recuperado el 15 de 02 de 2011, pp 6-13

Sánchez, Villafuerte, M. A. (Enero de 2009). Modelo de aplicaciones tecnológicas

para la protección de taludes desestabilizados por exceso de humedad en el
distrito metropolitano de quito a fin de disminuir la accidentalidad y
congestión vehicular. Obtenido de vias terrestres. Ambato Ecuador .pp 17-21.

Suárez, J., (1998), “Deslizamientos y estabilidad de taludes en zonas tropicales”,

Instituto de Investigaciones sobre erosión y deslizamientos, Ingeniería de
Suelos Ltda., Bucaramanga Colombia, pp 6-12

Wiliam, J. Burns, P. (2012). “ Protocol for Shallow – Landslide Susceptibility

Mapping”, Special Paper. P. 12

49