Tesis-2020-Ing. Civil-Huaman Martinez y Tito Ramos
Tesis-2020-Ing. Civil-Huaman Martinez y Tito Ramos
HUANCAVELICA
(Creada por Ley N° 25265)
Portada
TESIS
LÍNEA DE INVESTIGACIÓN
INFRAESTRUCTURA HIDRÁULICA SANEAMIENTO Y MEDIO AMBIENTE
PRESENTADO POR:
Bach. Huaman Martinez, Pool Jonhni
Bach. Tito Ramos, Nils Meyer
HUANCAVELICA, PERÚ
2020
Título
iii
Autores
iv
Asesor
v
Huaman Martinez, Pool Jonhni
vi
Agradecimiento
Al Msc. Ing. Leonardo Franco Castillo Navarro; por las excelentes clases
magistrales dictadas y sobre todo por el apoyo desinteresado en resolver ciertas dudas
que existían para la presente investigación sobre el uso de software Flo - 2D.
Al Msc. Ing. Iván Arturo Ayala Bizarro, asesor de la presente tesis; por el
apoyo y orientación en la elaboración de la misma, por habernos brindado la
oportunidad de compartir sus conocimientos de investigación; por tener una visión del
potencial humano de la tesis de pre grado que va más allá de la formación básica
recibida; por habernos demostrado a lo largo de estos años paciencia, comprensión y
ayuda incondicional.
Los autores
vii
Tabla de contenido
Portada .....................................................................................................................i
Autores ................................................................................................................... iv
Asesor ...................................................................................................................... v
Agradecimiento .....................................................................................................vii
Resumen ................................................................................................................ xx
Introducción ........................................................................................................xxii
CAPÍTULO I ........................................................................................................ 24
viii
1.4.3. Justificación social .............................................................................28
CAPÍTULO II ....................................................................................................... 29
ix
2.6.2. Variable 02 ...................................................................................... 129
4.2.1. Simulación del Flo - 2D en la Subcuenca del río Sicra ..................... 184
x
REFERENCIAS BIBLIOGRÁFICAS ............................................................... 201
A.1 Curvas IDF de las microcuencas de la subcuenca del río Sicra................. 208
A.2 Hietogramas de las microcuencas de la subcuenca del río Sicra ............... 213
A.4 Sección transversal del cauce del río Sicra para calibración del modelo
hidrológico HEC-HMS en la subcuenca del río Sicra ..................................... 229
C.2 Simulación de flujos de lodo y escombros para un Tr = 100 años ............ 234
C.3 Planos generados en entorno GIS, con los datos obtenidos del modelo Flo -
2D y el post procesador MAPPER PRO. ........................................................ 237
xi
Tabla de contenido de cuadros
Tabla 1. Clasificación de pendiente en las cuencas .................................................45
Tabla 6. Tipos de flujo de acuerdo con las características de la mezcla y del canal –
cauce ......................................................................................................................85
Tabla 8. Propiedades de las matrices de flujos de lodo (Limo y Arcilla) ............... 106
Tabla 10. Parámetros de resistencia del flujo laminar (Manual Flo - 2D) .............. 110
Tabla 12. Definición de nivel de amenaza según García, Rodríguez y O’Brien ..... 121
Tabla 13. Definición de la intensidad de inundación de agua (Water flood) .......... 123
Tabla 18. Resumen de áreas y perímetros de las microcuencas de la subcuenca del río
Sicra ..................................................................................................................... 141
Tabla 19. Ubicación de las estaciones climáticas virtuales empleadas en este estudio
............................................................................................................................. 143
xii
Tabla 20. Prueba de ajuste de bondad para datos de precipitación de cada microcuenca
de la subcuenca del río Sicra ................................................................................. 143
Tabla 24. Tiempo de concentración para cada microcuenca de la subcuenca del río
Sicra ..................................................................................................................... 148
Tabla 25. Cálculo del tiempo de retardo para cada microcuenca de la subcuenca del
río Sicra ................................................................................................................ 148
Tabla 27. Uso actual de los suelos para cada microcuenca de la subcuenca del río Sicra
............................................................................................................................. 153
Tabla 28. Valores del hidrograma líquido en el punto de afluencia ....................... 155
Tabla 29. Parámetros del método onda cinemática del cauce del río Sicra ............ 155
Tabla 32. Grado de susceptibilidad del suelo según Instituto Geológico y Metalúrgico
............................................................................................................................. 172
Tabla 33. Resultados del análisis granulométrico por tamizado (ASTM D-422) ... 178
Tabla 38. Salidas del programa Flo - 2D (Q = 225.00 m3/s - Tr = 50 años) .......... 185
Tabla 39. Salidas del programa Flo - 2D (Q = 349.30 m3/s - Tr = 100 años) ........186
xiii
Tabla 40. Resultado de volumen de sedimento depositado en el cauce del río Sicra
estimados de campo .............................................................................................. 196
xiv
Tabla de contenido de figuras
Figura 1. Vista satelital de la subcuenca del río Sicra Lircay, Huancavelica ...........25
Figura 22. Precipitación media diaria por cuenca para modelado hidrológico .........66
xv
Figura 23. Variables en el método de abstracciones de precipitación del SCS ........68
Figura 32. Movimiento idealizado de un flujo viscos (Johnson y Rodine, 1984) .....88
Figura 34. Metodología propuesta por O’Brien (2000) para el cálculo del hidrograma
de los flujos de lodo y escombros ............................................................................93
Figura 38. Procesos físicos simulados por el modelo Flo - 2D .............................. 115
Figura 40. Descarga del flujo a través de los elementos de grilla en planta ...........117
Figura 42. Intensidades de inundación de flujos de lodo y escombros (mud and debris
flow) basadas en profundidades 𝒉(𝒎) y velocidades de flujo 𝒗𝒎/𝒔 ..................... 122
xvi
Figura 44. Vista satelital del cauce del río Sicra – Prog. 0+565 km al 1+744 km .. 132
Figura 46. Curva hipsométrica y frecuencia de altitudes - subcuenca del río Sicra,
(cuenca en fase de juventud) ................................................................................. 139
Figura 49. Validación por diferentes medidas bondad de ajuste del producto Pisco con
los datos registrados en la estación meteorológica Lircay – 000627....................... 145
Figura 50. Mapa temático número de curva (CN) en condiciones normales .......... 151
Figura 51. Mapa temático de número de curva para la subcuenca del río Sicra ..... 152
Figura 53. Hidrograma líquido para la subcuenca del río Sicra (Tr=50 años y Tr=100
años) ..................................................................................................................... 156
Figura 54. Hidrograma líquido calibrado para la subcuenca del río Sicra (Tr=50 años
y Tr=100 años) ..................................................................................................... 157
Figura 56. Vista de zonas expuesta a inundación por las aguas del río Sicra .........159
Figura 57. Vista de zonas expuesta a inundación por las aguas del río Sicra .........160
Figura 59. Área de inundación margen derecha e izquierda río Sicra, zona urbano
marginal ciudad de Lircay, altura puente carrozable Santa Rosa aguas arriba ........161
Figura 60. Vista del puente Santa Rosa, en momentos del incremento del caudal de
las aguas del río Sicra, en proceso de erosión de las bases y estribos ..................... 162
Figura 61. Erosión de talud y pérdida parcial de plataforma, altura puente Chahuarma
del río Sicra .......................................................................................................... 163
Figura 62. Erosión y pérdida de talud de terrenos forestales y pastizales lado izquierdo
del cauce y trayectoria del río Sicra ....................................................................... 163
xvii
Figura 63. Vista de punto activo a deslizamientos N° 01 (D1) sobre el tramo alterno
en la vía Lircay – Ayacucho .................................................................................. 165
Figura 64. Vista de punto activo a deslizamientos N° 02 (D2) sobre el tramo alterno
en la vía Lircay – Ayacucho .................................................................................. 166
Figura 65. Vista de punto activo a deslizamientos N° 03 (D3) sobre el tramo alterno
en la vía Lircay – Ayacucho .................................................................................. 166
Figura 66. Vista de punto activo a deslizamientos N° 04 (D4) sobre el tramo alterno
en la vía Lircay – Ayacucho .................................................................................. 167
Figura 67. Vista de punto activo a deslizamientos N° 05 (D5) sobre el tramo alterno
en la vía Lircay – Ayacucho .................................................................................. 168
Figura 68. Vista de punto activo a deslizamientos N° 06 (D6) sobre el tramo alterno
en la vía Lircay – Ayacucho .................................................................................. 168
Figura 71. Mapa de peligro por movimientos en masa para el verano 2020 .......... 170
Figura 72. Mapa de peligro por inundaciones para el verano 2020........................ 171
Figura 74. Modelo digital de elevaciones de la subcuenca del río Sicra, del distrito de
Lircay - Angaraes ................................................................................................. 175
Figura 75. Grillado del modelo digital de elevaciones del distrito de Lircay, de 5mx5m
............................................................................................................................. 176
Figura 76. Hidrograma de caudal líquido y concentración de sólidos para Tr=50 años
............................................................................................................................. 177
Figura 77. Hidrograma de caudal líquido y concentración de sólidos para Tr=100 años
............................................................................................................................. 177
Figura 78. Ventana de Flo - 2D, simulación del flujo y desarrolla del hidrograma en
el tiempo ............................................................................................................... 185
xviii
Figura 80. Medición en campo el ancho de deslizamiento ocurrido, depositado en el
cauce del río Sicra ................................................................................................. 192
Figura 81. Identificación de puntos de control en zonas críticos caída de rocas hacia
el cauce del río Sicra ............................................................................................. 192
Figura 82. Identificación de puntos de control, Caída de rocas mayores 2.00 metros
de longitud, al cauce del río Sicra .......................................................................... 193
xix
Resumen
La subcuenca simulada en la presente tesis fue la subcuenca del río Sicra del
distrito de Lircay. De fuerte pendiente en ambos lados del cauce principal, bastante
acumulación de material en su cauce y con gran incremento de caudal debido a la
escorrentía que permite que se originen deslizamientos y derrumbes. En general los
taludes del cauce son bastantes inestables. La problemática es que, al producirse un
huayco, el material que se arrastra cae directamente sobre el río produciendo su
represamiento y consecuente desborde de flujos de lodo y escombros sobre el área
urbana de la ciudad de Lircay.
La subcuenca del río Sicra ha sido clasificada por CENEPRED, como una
subcuenca en "Peligro". Los eventos ocasionados nos recuerdan los sucesos de 2017 y
2019 como el huayco de Challhuapuquio cercana a dicha subcuenca que arrasó
completamente Instituciones educativas y viviendas aledañas.
xx
Abstract
The purpose of this research is to provide criteria for modeling mud and debris
flows using the Flo - 2D numerical model. This model simulates flow of non
newtonian fluids, such as torrential avalanches, in complex topographies, such as
urbanized areas and floodplains, water flow and hyperconcentrated flows (mud flow
and debris) can be modeled.
The sub-basin simulated in this thesis was the sub-basin of the Sicra river of
the district of Lircay. With a steep slope on both sides of the main channel, a lot of
material accumulation in its channel and with a large increase in flow due to runoff
that allows landslides and landslides to originate. In general, the slopes of the channel
are quite unstable. The problem is that, when a huayco occurs, the material that creeps
falls directly on the river producing its damming and consequent overflow of mud and
debris flows over the urban area of the city of Lircay.
The Sicra River sub basin has been classified by CENEPRED, as a sub-basin
in "Danger". The events caused remind us of the events of 2017 and 2019 as the huayco
of Challhuapuquio near said sub-basin that completely swept educational institutions
and surrounding homes.
With this work, a contribution is made in the study of mud and debris flows in
our region, initiating a very important line of research in terms of disaster mitigation.
At the same time, promote the use of this tool to define hazard maps of
hyperconcentrated flows; which facilitates the task of creating an environment of
conscience in the face of natural phenomena and informing the population in an easy
and technical way.
Key words: Mud and debris flows, model Flo - 2D, hazard maps.
xxi
Introducción
El Fenómeno del Niño (FEN) es uno de los grandes fenómenos naturales que
azotan cada año a nuestro país, las temporadas de fuertes lluvias desencadenan
inundaciones y huaycos. Esto ocurre esporádicamente y el resultado de su presencia
son pérdidas de viviendas, obras civiles, zonas de cultivo y sobretodo pérdidas
humanas.
xxii
definición de términos y definición operativa de variables e indicadores. CAPITULO
III: MATERIALES Y MÉTODOS, se encuentra referido al ámbito de estudio, el tipo
de investigación, el nivel de investigación, el método de investigación, del diseño de
investigación, población y muestra, técnicas e instrumentos de recolección de datos y
el procedimiento de recolección de datos. CAPITULO IV: DISCUSIÓN DE
RESULTADOS, se expone la presentación de resultados y la discusión.
xxiii
CAPÍTULO I
Según Mayta (2006), debido a que los flujos de lodo y escombros son
concentrados casi en su totalidad por masas sólidas, representan un peligro
latente para la población que vive a orillas de las quebradas, siendo, estas
formas geográficas las más comunes en la sierra peruana, y el lugar más
frecuente por donde se deslizan los fluidos de rocas y lodo, además la zona baja
24
de toda quebrada es considerada el lugar donde el deslizamiento muestra su
fuerza más destructiva enterrando todo a su paso. Por ende, es importante
considerar la teoría de flujo de hiperconcentrados aplicando modelos
numéricos como el Flo - 2D para, de esta manera mitigar los efectos de estos
flujos.
25
los ingenieros civiles en el diseño de obras, como puentes, carreteras, represas,
entre otras.
1.3. Objetivos
26
1.3.2. Objetivos específicos
1.4. Justificación
27
1.4.3. Justificación social
1.5. Limitaciones
28
CAPÍTULO II
MARCO TEÓRICO
2.1. Antecedentes
29
Oliveira (2012), realizó una investigación titulada “Estudio del impacto
de los flujos de escombros en la cuenca del río San Antonio - Brasil”, donde el
objetivo fue investigar la caracterización de la cuenca del río San Antonio, para
ello se empleó el modelo Flo - 2D. Con este modelo se pudo simular el
movimiento de agua sedimentada en zonas distorsionadas, las ecuaciones
usadas son principios de conservación de masa y medición del movimiento.
30
Además, es importante destacar que procedimiento se realizó a posteriori y, por
ende, la cartografía fue empleada posevento.
31
Pinaud (2014), elaboró una tesis titulada “Estudio experimental de las
características dinámicas de un dispositivo disipador viscoso a escala reducida”
en la Universidad Nacional de Ingeniería, Facultad de Ingeniería Civil, para
optar el título profesional de ingeniero civil. Los resultados del ensayo de carga
monotónica, y ensayo dinámico, ambos confirman un comportamiento bilineal
respecto a la fuerza de amortiguamiento. Cabe resaltar que se obtuvo mayor
amortiguamiento en el ensayo monotónico, esto debido a dos causas, en primer
lugar, debido a que el flujo del líquido en el disipador en el ensayo de la mesa
vibradora no se desarrolla en toda la región anular, segundo es causa de las
características del flujo, ya que en la prueba sobre la mesa vibradora el
movimiento es cíclico desarrollándose turbulencias y discontinuidades, lo cual
hace que el fluido fluya en muchas direcciones, lo cual reduce la presión en la
dirección del movimiento del pistón dando como resultado la reducción del
amortiguamiento. En la otra prueba con movimiento sinusoidal en la mesa
vibradora no se encontró relación alguna entre las frecuencias de movimiento
y el desarrollo de la fuerza de amortiguamiento del disipador.
32
hiperconcentrados, es necesario definir el hidrograma líquido, la correcta
distribución de la concentración de sedimentos y los parámetros reológicos de
la mezcla, así mismo el modelo Flo - 2D emplea el modelo cuadrático de
O´Brien y Julien, concluyendo que las zonas más vulnerables ante un desastre
son los terrenos ubicados en la zona la progresiva la que generaría grandes
pérdidas económicas y humanas.
33
el Titán 2D, así mismo se emplearon instrumentos topográficos, pluviómetros
y la geotecnia.
34
La cuenca de drenaje de una corriente, es el área de terreno donde todas
las aguas caídas por precipitación, se unen para formar un solo curso de agua.
Cada curso de agua tiene una cuenca bien definida, para cada punto de su
recorrido (Villón Béjar, 2002).
[Link]. Clasificación
35
(pendiente, elevación, área, cauce). El efecto de almacenaje del cauce es
muy importante.
36
Fuente: (Aparicio Mijares, 2001)
37
Fuente: (Aparicio Mijares, 2001)
38
[Link]. Elementos
39
Fuente: Adaptado de (Aparicio Mijares, 2001)
40
Figura 8. Influencia de la forma de la cuenca en el hidrograma
41
En cuanto al perímetro de la cuenca que se muestra en la Figura 8,
informa sucintamente sobre la forma de la cuenca; para una misma
superficie, los perímetros de mayor valor corresponden a cuencas alargadas,
mientras que los de menor corresponden a cuencas redondeadas. La
distancia equivalente que recorre el río entre el punto de salida aguas abajo
y el punto situado a mayor distancia topográfica aguas arriba (Ibáñez
Asensio, Moreno Ramón, & Gisbert Blanquer, 2011).
Fuente: Adaptado de (Ibáñez Asensio, Moreno Ramón, & Gisbert Blanquer, 2011)
𝐴
𝑊= ( II.1 )
𝐿𝑐
Donde:
42
[Link].2. Pendiente media
Criterio de Alvord.
Criterio de Horton.
Criterio de Nash.
𝑛
1
𝐼𝑝 = ∑ √𝛽𝑖 (𝑎𝑖 − 𝑎𝑖−1 ) ∗ ( II.2 )
𝑖=2
√𝐿
𝐴𝑖
𝛽𝑖 = ( II.3 )
𝐴𝑇
Donde:
43
Ip : Índice de pendiente
44
Tabla 1. Clasificación de pendiente en las cuencas
Pendiente media
Terrenos
(%)
0a2 Llano
2a5 Suave
5 a 10 Accidentado medio
10 a 15 Accidentado
15 a 25 Fuertemente accidentado
25 a 50 Escarpado
>50 Muy escarpado
45
Cabe mencionar que, ingresando el 50% del área en el eje de las
abscisas se obtiene la altitud media de la cuenca que intercepta con la curva
hipsométrica Figura 10.
𝐴𝑠
𝑅𝐻 = ( II.4 )
𝐴𝑏𝑖
Donde:
Donde:
47
mecanismos. Si un paisaje muestra un carácter permanente, estos dos
procesos opuestos están en equilibrio dinámico.
48
Fuente: Adaptado de (Agustín, 2009)
ℎ
𝑅𝑟 = ( II.5 )
𝐿
49
Por tener el concepto de tiempo de concentración una cierta base
física, han sido numerosos los autores que han obtenido formulaciones del
mismo, a partir de características morfológicas y geométricas de la cuenca.
A continuación, se muestran algunas de esas fórmulas empíricas
(Montenegro, 2005):
Kirpich:
0.385
𝐿3
𝑇𝑐 = 0.0195 ∗ ( ) ( II.6 )
𝐻
Temez:
0.75
𝐿𝑝
𝑇𝑐 = 0.126 ∗ ( ) ( II.7 )
𝑆𝑃 0.35
Pasini:
0.5
𝐴 ∗ 𝐿𝑝
𝑇𝑐 = 0.023 ∗ ( ) ( II.8 )
𝑆𝑝
Pizarro:
0.77
𝐿2
𝑇𝑐 = 13.548 ( ) ( II.9 )
𝐻
Donde:
50
H : Diferencia de cotas (m)
𝐿𝑡
𝐷𝑑 = ( II.10 )
𝐴
51
el área y cuatro veces la longitud de todos los cauces de la cuenca, o bien,
la inversa de cuatro veces la densidad de drenaje.
𝐴 1
𝐿𝑜 = = ( II.11 )
4 ∑ 𝐿𝑖 4𝐷
𝐻𝑚á𝑥 − 𝐻𝑚í𝑛
𝑆𝑚 = ( II.12 )
𝐿𝑐
52
Fuente: Adaptado de (Cahuana A. & Yugar M., 2009)
2.2.2. Precipitación
53
[Link]. Precipitaciones (PISCO - SENAMHI)
54
terrenas, obteniendo resultados óptimos para algunas zonas del país y bajos
para otras. La meta propuesta que inspiró dichos trabajos fue mejorar la
representación espacial de lluvias en el Perú usando los datos del
sensoramiento remoto como covariables para su asimilación en modelos
hidrológicos y desarrollo de productos para monitoreo de sequías e
inundaciones.
55
(SONICS), monitoreo decadiario de precipitaciones del SENAMHI
(MIDAS), elaboración de los mapas de umbrales de precipitaciones
extremas a nivel nacional, mapas hidroclimáticos de cuencas, vigilancia
nacional de la sequía, para la elaboración de estudios e investigaciones
hidroclimáticas entre otros.
56
quiebres evidentes en las series de tiempo y eliminación de segmentos
que no eran homogéneos. Además, la completación de datos de las series
mensuales y diarias fue llevado a cabo basado en una modificación del
algoritmo CUTOFF (Feng et al., 2014) y mapeos de cuantiles ("quantile
mapping") respectivamente. Finalizando el control de calidad de datos,
solo 443 de las 681 estaciones pluviométricas fueron seleccionadas para
la generación de datos PISCOp. La cantidad de años de información en
% considerando el período 1981-2016 por cada estación pluviométrica
antes de la completación de datos puede apreciarse en la Figura 14.
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.;
Sabino, E. & Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de
precipitación en Estudios, Investigaciones y Sistemas Operacionales de Monitoreo y
Pronóstico Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
57
datos únicamente en las zonas terrestres, comprendida entre las latitudes
50° S y 50° N y en longitudes 0° a 360°; teniendo una resolución espacial
de 0.05°. Las resoluciones temporales son diarias, pentadiarias,
decadiarias, mensuales y anuales. Las fuentes de datos utilizadas en la
creación de CHIRPS fueron: (1) CHPClim (Climate Hazards
Precipitation Climatology); (2) Infrarrojo (IR) térmico geoestacionario;
(3) TRMM (Tropical Rainfall Measuring Mission); (4) Modelo
atmosférico de campos de precipitación del sistema de predicción
climático de la NOAA (CFSv2); y (5) Observaciones in situ de
precipitación obtenidas de varias fuentes incluidos servicios
meteorológicos estadounidenses y nacionales.
Donde:
58
precipitación en Perú, se optó por utilizar una climatología basada en
TRMM2A25, obteniéndose de este modo una modificación del CHIRP
(CHIRPM) que fue llevada a cabo a paso diario (CHIRPMd) y mensual
(CHIRPMm).
𝑇𝑅𝑀𝑀2𝐴25 + 𝐸
𝐶𝐻𝐼𝑅𝑃𝑀 = (𝑏0 + 𝑏1 . (𝐶𝐶𝐷)). 𝐶𝐻𝑃𝑐𝑙𝑖𝑚 . ( II.14 )
𝐶𝐻𝑃𝑐𝑙𝑖𝑚 + 𝐸
Donde:
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.;
Sabino, E. & Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de
59
precipitación en Estudios, Investigaciones y Sistemas Operacionales de Monitoreo y
Pronóstico Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
60
PISCOpm y PISCOpd evidencia una mejora significativa del
coeficiente de correlación (CC) frente a CHIRPd y CHIRPm sobretodo en
la costa y sierra, esto está relacionado a la distancia entre las estaciones
utilizadas en la validación y las estaciones utilizadas en PISCOp, el bajo
desempeño de PISCOp en la selva es debido a la baja densidad de estaciones
que imposibilita la caracterización de la estructura espacial (mediante el
variograma) de forma adecuada. A diferencia del CC que nos permite definir
la variabilidad de las precipitaciones, tanto el BIAS % como el RMSE nos
permite obtener información sobre las intensidades. Considerando ello,
dichos estadísticos presentan un mejor rendimiento principalmente en la
sierra y selva, presentando fuertes sobreestimaciones en toda la costa, sin
embargo, esta región está caracterizada por presentar precipitaciones por
debajo de 100 mm anuales (Manz et. al, 2016) por lo que el bajo rendimiento
puede ser despreciado.
61
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.;
Sabino, E. & Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de
precipitación en Estudios, Investigaciones y Sistemas Operacionales de monitoreo y
Pronóstico Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
CHIRP
Mensual KED y Quiebres leves
2 TRMM2A2 SI 1981-2016
y diario RIDWO artificiales
5
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.; Sabino, E.
& Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de precipitación en
Estudios, Investigaciones y Sistemas Operacionales de monitoreo y Pronóstico
Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.; Sabino, E.
& Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de precipitación en
Estudios, Investigaciones y Sistemas Operacionales de monitoreo y Pronóstico
Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
63
inundaciones del SENAMHI – SONICS, monitoreo decadiario de la
precipitación – MIDAS, modelamiento hidrológico de cuencas, etc).
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.; Sabino, E.
& Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de precipitación en
Estudios, Investigaciones y Sistemas Operacionales de monitoreo y Pronóstico
Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
64
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.; Sabino, E.
& Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de precipitación en
Estudios, Investigaciones y Sistemas Operacionales de monitoreo y Pronóstico
Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.; Sabino, E.
& Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de precipitación en
Estudios, Investigaciones y Sistemas Operacionales de monitoreo y Pronóstico
Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
65
Fuente: Aybar, C.; Lavado-Casimiro, W.; Huerta, A.; Fernández, C.; Vega, F.; Sabino, E.
& Felipe-Obando, O. (2017). Uso del Producto Grillado “PISCO” de precipitación en
Estudios, Investigaciones y Sistemas Operacionales de monitoreo y Pronóstico
Hidrometeorológico. Nota Técnica 001 SENAMHI-DHI-2017, Lima-Perú
Figura 22. Precipitación media diaria por cuenca para modelado hidrológico
Fuente: Felipe, O
66
una cuenca mediante las características de la misma y la precipitación
(Aparicio Mijares, 2001).
67
Figura 23. Variables en el método de abstracciones de precipitación del SCS
Ia : Abstracción inicial
Pe : Exceso de precipitación
Fa : Abstracción continuada
P : Precipitación total
𝐹𝑎 𝑃𝑒
= ( II.15 )
𝑆 𝑃 − 𝐼𝑎
𝑃 = 𝑃𝑒 + 𝐼𝑎 + 𝐹𝑎 ( II.16 )
(𝑃 − 𝐼𝑎 )2
𝑃𝑒 = ( II.17 )
𝑃 − 𝐼𝑎 + 𝑆
68
NRCS. La abstracción inicial consiste principalmente en intercepción,
infiltración durante los comienzos de la tormenta, y las depresiones o
almacenamientos en la superficie, establecen una relación empírica para
estimar esta abstracción inicial Ia, que está en función al potencial máximo
de retención S.
𝐼𝑎 = 0.2𝑆 ( II.18 )
(𝑃 − 0.2 ∗ 𝑆)2
𝑃𝑒 = ( II.19 )
𝑃 − 𝐼𝑎 + 𝑆
69
Fuente: Adaptado de (Ven te Chow , 1994)
1000
𝑆= − 10 ( II.20 )
𝐶𝑁
25400
𝐶𝑁 = ( II.21 )
25 + 𝑆
Donde:
4.2𝐶𝑁(𝐼𝐼 )
𝐶𝑁(𝐼 ) = ( II.22 )
10 − 0.058 ∗ 𝐶𝑁(𝐼𝐼 )
23𝐶𝑁(𝐼𝐼)
𝐶𝑁(𝐼𝐼𝐼) = ( II.23 )
10 − 0.13 ∗ 𝐶𝑁(𝐼𝐼)
70
Estación de
Estación inactiva
crecimiento
GRUPO DESCRIPCIÓN
71
Alto potencial de escurrimiento. Suelos que tienen muy bajas
velocidades de infiltración cuando están mojados y consisten
principalmente de suelos arcillosos con alto potencial de
GRUPO D
hinchamiento, suelos con estratos arcillosos cerca de su
superficie, o bien, suelos someros sobre horizontes
impermeables.
Principios básicos:
73
Fuente: Adaptado de (Ven te Chow , 1994)
Donde:
74
Pe : Es el volumen de lluvia efectiva.
Ie : Intensidad efectiva.
75
2.2.3. Hidráulica de flujos de escombros
76
En general, múltiples autores concuerdan en que la estructura de un
flujo de escombros es fácilmente reconocible (Figura 27), en la cual se pueden
identificar de forma sencilla tres regiones claramente definidas (Prasad &
Rajaratnam, 1991). Asimismo, la concentración de los sólidos varía
notablemente en función a la región que se evalúa, así como también su
comportamiento.
Cola: El desarrollo del flujo en esta zona suele ser un flujo inmaduro,
que tiende a convertirse en un flujo de lodos o un cauce en estado de
hiperconcentración. Esta suele depositar los sedimentos más finos
acarreados por el paso de las regiones previamente descritas.
77
[Link]. Tipos de fluidos
78
newtonianos “hidráulica torrencial”, estos se dividen en tres tipos de
fluidos dilatantes, pseudoplasticos y de tipo Bingham.
Donde:
79
B : El fluido escurre más fácilmente cuanto mayor es el esfuerzo.
(ejm: barro)
Existe una saturación del agua por parte del material solido que
presenta una alta gradación de diámetros. Las partículas más finas de arcilla
generan una matriz de soporte que ofrece continuidad al elemento sólido, el
cual se lubrica con el agua al moverse.
Flujos de lodo:
Flujos de escombros:
80
Los flujos de escombros es el resultado de alguna forma de colapso
en el talud. Los escombros que caen como un deslizamiento colectan
humedad y se mueven a lo largo de la pendiente, éste se licua o se dilata
conforme avanza, aumentando la movilidad del fluido (Hamptom, 1972).
81
Fuente: Adaptado de (O’Brien, 2000)
82
Densidad τy
Concentración
Flujo los sólidos 2 Tipo de Fluido
de sedimentos 3 (dn/cm )
(g/cm )
1-40% en peso
Avenida de Agua 0.4-20% en 1.01-1.33 0-100 Newtoniano
volumen
40-70% en
Flujo peso No
1.33-1.80 100-400
Hiperconcentrado 20-47% en Newtoniano
volumen
70-90% en
Flujo de peso
1.80-2.30 >400 Viscoplástico
escombros 47-775 en
volumen
Mayor Perfil de
mecanismo de Viscosidad concentración Tipo de flujo
Flujo
soporte de los de predominante
sedimentos (poise) sedimentos
Fuerzas
Avenida de agua electrostática, 0.01-20 No uniforme Turbulento
turbulencia
Empuje,
Flujo esfuerzo No uniforme Turbulento a
20-200
Hiperconcentrado dispersivo, a uniforme Laminar
turbulencia
Cohesión,
empuje,
Flujo de esfuerzo
>>200 Uniforme Laminar
escombros dispersivo,
soporte
estructural
83
Fuente: Adaptado de (Castillo, 2006)
n
lia
Ju
an ela)
Altura con respecto
Shiramizudani o S u
Rí nez
(China) nt (V
e
rre
To 2
ra 7 s ta
al pie del flujo
Xikon a lia km ua
ag 99 me .6 aig la)
(China) m )1 ca ia) 23
G a pó n as omb o N ue
a L Rí nez
(J ol m2
562 k
e
(C 2 (V
m2 2 km
4k km .4 olombia)
.5 33 (C
2 3 Paez 2
84
cauce principal. Para resumir los criterios generales y para determinar
el tipo de flujo, se presenta la Tabla 6.
Tabla 6. Tipos de flujo de acuerdo con las características de la mezcla y del canal
– cauce
85
gruesas, así como envolver a todo el material que se presente en su
camino.
Arena Arena
Arcilla Limo Grava Cantos Bolones
Fina Gruesa
2um 20um 200um 2mm 2cm 0.2 m
86
[Link]. Propiedades de los flujos de escombros
[Link].1. Velocidad
Donde:
𝜔 : Velocidad promedio.
ѱ : Radio de la curva.
87
Figura 32. Movimiento idealizado de un flujo viscos (Johnson y Rodine, 1984)
[Link].2. Pendiente
88
Por medio de un análisis de estabilidad simple para un ancho y
longitud unitario completamente saturado, de altura h, con pendiente
𝑠𝑒𝑛(𝜃) y con ángulo de fricción ∅, se obtiene el criterio de iniciación del
movimiento de flujos de escombros.
De la mecánica de suelos:
𝜎̅ = 𝜎 − 𝑢𝑝 ( II.26 )
𝑢𝑝 = 𝛾𝑤 . ℎ. cos(𝜃) ( II.27 )
W sin
W cos
Del gráfico:
𝑤 = 𝛾𝑠𝑎𝑡 . ℎ ( II.28 )
89
𝜎 = 𝑊. cos(𝜃) ( II.29 )
𝜏 = 𝑊. 𝑠𝑒𝑛(𝜃) ( II.30 )
𝛾𝑠𝑎𝑡 − 𝛾𝑤
tan(𝜃) = . tan(𝜙) ( II.31 )
𝛾𝑠𝑎𝑡
𝐶𝑣 . (𝛾𝑠𝑎𝑡 − 𝛾𝑤 )
tan(𝜃) = . tan(𝜙) ( II.32 )
𝐶𝑣 . (𝛾𝑠𝑎𝑡 − 𝛾𝑤 ) + 𝛾𝑤
90
Finalmente se sugiere asumir un rango de ángulos de deposición
entre 10° y 14° para flujos abiertos o no confinados y de 8° a 12° para flujos
canalizados o confinados.
𝐶𝑑
𝑉𝑠 = 𝑉𝑎 . ( ) ( II.33 )
𝐶∗ − 𝐶𝑑
Donde:
𝑉𝑠 : Volumen de escombros en 𝑚3 .
𝑉𝑎 : Volumen de agua en 𝑚3 .
𝜌. tan(𝜃)
𝐶𝑑 = ( II.34 )
(𝜎 − 𝜌)(tan(𝜙) − tan(𝜃))
Donde:
𝜙 : Angulo de fricción.
𝐶∗
𝑄𝑠 = ( ) . 𝑄𝑝 ( II.35 )
𝐶 ∗ − 𝐶𝑑
𝜌−1
𝑄𝑓 = ( ) . 𝑄𝑝 ( II.36 )
𝜌−𝜎
𝑄𝑡 = 𝑄𝑠 + 𝑄𝑓 ( II.37 )
𝑉𝑠 𝜌−1
𝑄𝑡 = ((1 + )+( )) . 𝑄𝑝 ( II.38 )
𝑉𝑎 𝜎−𝜌
Donde:
𝑉𝑠 : Volumen de escombros en 𝑚3 .
𝑉𝑎 : Volumen de agua en 𝑚3 .
92
hidrograma de flujo de agua multiplicado por un factor que es una
función de la concentración volumétrica de sedimentos (Cv):
𝑄𝑡 = 𝐵𝐹. 𝑄𝐿 ( II.39 )
1
𝐵𝐹 = ( II.40 )
1 − 𝐶𝑣
Donde:
Figura 34. Metodología propuesta por O’Brien (2000) para el cálculo del hidrograma de
los flujos de lodo y escombros
93
1400
1200
1000
600
400
200
Time (hours)
94
La otra posibilidad considera el agua (en algunos casos con finos) y
material grueso por separado. A menudo se refiere a modelos de
flujo granular, estos modelos se basan íntimamente en el
conocimiento de las partículas y de su distribución de tamaño.
𝜕𝑉
𝜏 = 𝜂. ( ) ( II.41 )
𝜕𝑦
Donde:
𝜏 : Esfuerzo cortante.
𝜂 : Es el parámetro de viscosidad.
𝜕𝑉
: Tasa de deformación al corte o gradiente de velocidades.
𝜕𝑦
95
𝜏 = 𝑐 + 𝜎𝑛 . tan(𝜙) ( II.42 )
Donde:
𝑐 : Cohesión.
𝜎𝑛 : Esfuerzo normal.
𝜕𝑉
𝜏 = 𝑐 + 𝜎𝑛 . tan(𝜙) + 𝜂. ( ) ( II.43 )
𝜕𝑦
Donde:
𝜏 : Esfuerzo cortante.
𝜂 : Es el parámetro de viscosidad.
𝜕𝑉
: Tasa de deformación al corte o gradiente de velocidades.
𝜕𝑦
𝜕𝑉
𝜏(𝑦) = 𝜏𝑦 + 𝜂. ( ) ( II.44 )
𝜕𝑦
Donde:
96
𝜏 : Esfuerzo cortante.
𝜂 : Es el parámetro de viscosidad.
𝜕𝑉
: Tasa de deformación al corte o gradiente de velocidades.
𝜕𝑦
2
𝜕𝑉
𝜏 = 𝜂. ( ) ( II.45 )
𝜕𝑦
Donde:
𝜏 : Esfuerzo cortante.
𝜂 : Es el parámetro de viscosidad.
𝜕𝑉
: Tasa de deformación al corte o gradiente de velocidades.
𝜕𝑦
97
2
𝜕𝑉 𝜕𝑉
𝜏 = 𝜏𝑦 + 𝜂. ( ) + 𝐶1. ( ) ( II.46 )
𝜕𝑦 𝜕𝑦
Donde:
𝜏 : Esfuerzo cortante.
𝜂 : Es el parámetro de viscosidad.
𝐶1 : Es el parámetro turbulento-dispersivo.
𝜕𝑉
: Tasa de deformación al corte o gradiente de velocidades.
𝜕𝑦
98
físicos que gobiernan el intercambio de flujo entre el canal y la planicie de
inundación. Así mismo es un modelo de diferencias finitas que integra la
ecuación de momentum en forma completa, utilizando un esquema geométrico
con celdas de dimensión constante y finalizado en el cálculo celda por celda y
en 8 direcciones de las variables hidráulicas que como resultado brinda
parámetros aproximados como:
1. Profundidad.
2. Velocidad.
3. Caudal.
Ecuación de continuidad:
𝜕ℎ 𝜕ℎ𝑉𝑥 𝜕ℎ𝑉𝑦
+ + =𝑖 ( II.47 )
𝜕𝑡 𝜕𝑥 𝜕𝑦
Donde:
ℎ : Es la profundidad hidráulica.
99
𝑡 : Es el tiempo.
Ecuación de Movimiento:
𝜕ℎ 𝑉𝑥 𝜕𝑉𝑥 𝑉𝑦 𝜕𝑥 1 𝜕𝑥
𝑆𝑓𝑥 = 𝑆𝑜𝑥 − − . − . − . ( II.48 )
𝜕𝑥 𝑔 𝜕𝑥 𝑔 𝜕𝑦 𝑔 𝜕𝑡
𝜕ℎ 𝑉𝑦 𝜕𝑉𝑦 𝑉𝑥 𝜕𝑦 1 𝜕𝑦
𝑆𝑓𝑦 = 𝑆𝑜𝑦 − − . − . − . ( II.49 )
𝜕𝑦 𝑔 𝜕𝑦 𝑔 𝜕𝑥 𝑔 𝜕𝑡
Donde:
𝑔 : Aceleración de la gravedad.
(O’Brien & Julien, 1988), Majo & Pierson, (1990) y (Lan, 1991),
investigaron flujos de lodo con altas concentraciones de sedimento fino en la
matriz fluida. Estos estudios mostraron que los flujos de lodo se comportan
100
como los fluidos de Bingham con bajas tasas de corte. En matrices fluidas
con bajas concentraciones de sedimento, los esfuerzos turbulentos dominan
en el flujo. Altas concentraciones de partículas gruesas combinadas con bajas
concentraciones de partículas finas son requeridas para generar esfuerzos
dispersivos. El modelo cuadrático de esfuerzo al corte propuesto por (O’Brien
& Julien, 1985), describe el régimen continuo de flujo viscoso a flujo
turbulento/dispersivo.
101
𝜏 = 𝜏𝑐 + 𝜏𝑚𝑐 + 𝜏𝑉 + 𝜏𝑡 + 𝜏𝑑 ( II.50 )
Donde:
𝜕𝑉 𝜕𝑉 2
𝜏 = 𝜏𝑦 + 𝜂 ( ) + 𝐶 ( ) ( II.51 )
𝜕𝑌 𝜕𝑌
Donde:
𝜏𝑦 = 𝜏𝑐 + 𝜏𝑚𝑐 ( II.52 )
102
𝐶𝑣 : Es la concentración volumétrica de sedimentos.
𝐶∗ 1/3
𝑓(𝜌𝑚 , 𝐶𝑣 ) = 𝛼𝑖 ∗ 𝜌𝑚 [( ) − 1] ( II.54 )
𝐶𝑣
103
Fuente: Adaptado de (Ccama Lipa, 2018)
𝑆𝑓 = 𝑆𝑦 + 𝑆𝑉 + 𝑆𝑡𝑑 ( II.55 )
Donde:
𝑆𝑦 : Gradiente de cedencia.
𝑆𝑉 : Gradiente viscoso.
104
Los términos de gradiente viscoso y turbulento pueden escribirse en
términos de profundidad y velocidad media, donde la pendiente viscosa se
puede escribir como.
𝐾. 𝜂 𝑉
𝑆𝑣 = . ( II.56 )
8. 𝑦𝑚 ℎ2
Donde:
𝑛𝑡𝑑 2 . 𝑉 2
𝑆𝑡𝑑 = ( II.57 )
ℎ4/3
𝜏𝑦 𝐾𝜂 𝑉 𝑛𝑡𝑑 2 . 𝑉 2
𝑆𝑓 = + . + ( II.58 )
𝑦𝑚 . ℎ 8. 𝑦𝑚 ℎ2 ℎ4/3
105
Tabla 8. Propiedades de las matrices de flujos de lodo (Limo y Arcilla)
𝜏𝑦 = 𝛼1 𝑒 𝛽1 𝐶𝑣 (II.59)
𝜂 = 𝛼2 𝑒 𝛽2 𝐶𝑣 (II.60)
106
Donde:
𝜂𝑡 : Es la viscosidad turbulenta.
𝜂𝑑 : Es la viscosidad dispersiva.
107
Fuente: Adaptado de (Castillo, 2006)
𝜏𝑦 = 𝛼1 . 𝑒 𝛽1.𝐶𝑣 𝜂 = 𝛼2 . 𝑒 𝛽2.𝐶𝑣
Source
𝛼1 𝛽1 𝛼2 𝛽2
Field Data
Aspen pit 1 0.181 25.7 0.036 22.1
Aspen pit 2 2.72 10.4 0.0538 14.5
Aspen natural soil 0.152 18.7 0.00136 28.4
108
Aspen mine fill 0.0473 21.1 0.128 12
Aspen watershed 0.0383 19.6 0.000495 27.1
Aspen mine source area 0.291 14.3 0.000201 33.1
Glenwood 1 0.0345 20.1 0.00283 23
Glenwood 2 0.0765 16.9 0.0648 6.2
Glenwood 3 0.000707 29.8 0.00632 19.9
Glenwood 4 0.00172 29.5 0.000602 33.1
Relationships Available from the Literature
Iida (1938) - - 0.0000373 36.6
Dai et al. (1980) 2.6 17.48 0.0075 14.39
Kang and Zhang (1980) 1.75 7.82 0.0405 8.29
0.00136 21.2 - -
Qian et al. (1980)
0.05 15.48 - -
Chien and Ma (1958) 0.0588 19.1-32.7 - -
0.166 25.6 - -
Fei (1981)
0.0047 22.2 - -
109
Si el valor de K es ingresado, éste será utilizado por el modelo. Sin
embargo si K = 0, el modelo computará K en función de la rugosidad de
Manning:
Tabla 10. Parámetros de resistencia del flujo laminar (Manual Flo - 2D)
Superficie Rango de K
Concreto/asfalto 24 – 108
Arena escasa 30 – 120
Superficie gradual 90 – 400
Arcilla escasa – suelo de marga erocionada 100 – 500
Poca vegetación 1000 – 4,000
Hierba corta de pradera 3,000 – 10,000
Cesped de hierba añil 7,000 – 50,000
110
pico del hidrograma líquido. Una clasificación dada por O'Brien en su manual
de Flo - 2D se observa en la Tabla 11. Eventos de avenidas grandes tal como
avenidas de 100 años puede contener demasiada agua para producir eventos
de flujos de lodo viscosos. Eventos de lluvias más pequeñas tal como
tormentas con período de retorno de 10 a 25 años puede tener mayor tendencia
a crear flujos de lodo viscosos. La mayoría de cuencas con una historia de
eventos de flujos de lodo y escombros tendrán un suministro de sedimento
sustancial. Aún, tormentas muy pequeñas pueden generar oleadas de flujo de
lodo.
Concentración de sedimentos
Descripción del
En volumen En peso Características del flujo
tipo de flujo
No hay flujo; falla por deslizamiento de
0,65 – 0,80 0,83 – 0,91
bloques.
Derrumbe de bloques con deformación
Deslizamientos
interna durante el deslizamiento,
0,55 – 0,65 0,76 – 0,83
movimiento paulatino del terreno antes
de fallar.
Flujo evidente; deslizamiento lento
sostenido por el flujo de lodo;
0,48 – 0,55 0,72 – 0,76 deformación plástica bajo su propio
Flujos de lodo peso; cohesivo; no se expande en la
(Mud flow) superficie.
Flujo se extiende sobre la superficie;
0,45 – 0,48 0,69 – 0,72
flujo cohesivo; algo de mezcla.
111
El flujo se mezcla fácilmente; muestra
las propiedades fluidas en la
deformación; distribuido en la
superficie horizontal, pero mantiene una
0,40 – 0,45 0,65 – 0,69
superficie fluida inclinada; partículas
grandes se depositan (peñones);
aparecen ondas, pero se disipan
rápidamente.
Deposición marcada de gravas y cantos
rodados; se expande casi
Avenida de lodo completamente sobre la superficie
0,35 – 0,40 0,59 – 0,65
(Mud flood) horizontal; aparece la superficie líquida
con dos fases del fluido; las olas viajan
en la superficie.
Separación del agua en la superficie; las
olas viajan fácilmente; la mayoría de las
0,30 – 0,35 0,54 – 0,59
arenas y gravas se han sedimentado y se
mueven como arrastre de fondo
Se distinguen claramente las olas y
0,20 – 0,30 0,41 – 0,54 ondas; todas las partículas descansando
en el fondo en condición inmóvil.
Inundación de agua con carga
Inundación de
< 0,20 < 0,41 suspendida convencional y arrastre de
agua (Water flood)
fondo.
𝐶𝑤 . 𝛾
𝐶𝑣 = ( II.66)
𝛾𝑠 − 𝐶𝑤 . (𝛾𝑠 − 𝛾)
Donde:
112
𝛾𝑠 : Peso específico del sedimento.
𝛾𝑚 = 𝛾 + 𝐶𝑣 . (𝛾𝑠 − 𝛾) ( II.67)
∆𝑥
∆𝑡 = 𝐶. ( II.68)
(𝑣 + 𝑐)
Donde:
113
∆𝑥 : Es el ancho del elemento de la malla cuadrada.
𝜁𝑆𝑜 ∆𝑥 2
∆𝑡 < ( II.69)
𝑞𝑜
114
El modelo hidráulico Flo - 2D ha sido aprobado por FEMA para
estudios ribereños y análisis de avenidas. Entre los principales atributos de
este modelo se encuentran:
Transporte de sedimentos.
Estructuras hidráulicas.
115
Fuente: User’s manual Flo - 2D
116
La topografía de la superficie se discretiza en una malla formada por
celdas o elementos cuadrados de tamaño uniforme para toda el área de
estudio y a cada elemento se le asigna una posición en la malla, una
elevación o cota de terreno, un coeficiente de rugosidad (𝑛 de Manning) y
factores de reducción del flujo que atraviesa la celda.
Figura 40. Descarga del flujo a través de los elementos de grilla en planta
117
Fuente: Adaptado de ( (Escobar Soto, 2018)
118
4. Si la ecuación de onda dinámica completa es seleccionada por el
usuario, la velocidad de onda difusiva prevista para ese paso de
tiempo es usado como un dato inicial en la solución de Newton -
Raphson. Se hace nota que para flujos de sedimentos
hiperconcentrados, como flujos de lodo y escombros, el cálculo de
la velocidad incluye los términos adicionales de viscosidad y
esfuerzo de cedencia.
119
10. La simulación prosigue incrementándose el paso de tiempo hasta
que el criterio de estabilidad sea excedido.
USA : 𝑓(ℎ)
𝑣2
Austria : 𝑓(ℎ, 2.𝑔)
Switzerland : 𝑓(ℎ, 𝑣. ℎ)
120
Figura 41. Niveles de peligro de inundación basados en l frecuencia e intensidad de
inundación
INTENSIDAD
Alta
Media
Baja
Muy
FRECUENCIA Alta Media Baja
baja
Periodo de retorno 10 100 500 >>500
Probabilidad de excedencia 10% 1% 0.20% <<0.20%
Niveles de Color en
Descripción
amenaza el mapa
Las personas están en peligro, tanto dentro como
Alto Rojo fuera de las casas. Las edificaciones pueden ser
destruidas
121
Fuente: Adaptado de Office Federal de i OEconomie Des Aux (OFEE), 1997
Figura 42. Intensidades de inundación de flujos de lodo y escombros (mud and debris
flow) basadas en profundidades 𝒉(𝒎) y velocidades de flujo 𝒗(𝒎/𝒔)
H
2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.9 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.8 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.7 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.6 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.5 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.4 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1.1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
0.9 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
0.8 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
0.7 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
0.6 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3
0.5 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3
0.4 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
0.3 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
0.2 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
0.1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
VxH 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 2.05 2.1 2.15
122
En MAPPER podemos introducir valores de profundidad de flujo y
velocidades que definen los umbrales de intensidad. Para el caso de las
inundaciones de agua, las intensidades de inundación podrían definirse por
los valores de la Tabla 13, mientras que los flujos de lodos y escombros que
son más destructivos que las inundaciones de agua se definen por los valores
de la Tabla 14, los criterios de intensidad de los flujos de lodo y escombros
son más complejos.
123
elementos de la cuadrícula representa el riesgo de inundación bajo, medio y
alto.
2.3.1. Subcuenca
[Link]. Geomorfológicas
125
[Link]. Geodinámicos
126
y arcillas), las cuales son descritos como flujos viscosos además la
deposición final del flujo mantendrá una estructura dinámica en su
movimiento originando filas de dunas en dirección del flujo, durante la
sedimentación las partículas se mueven de manera integral.
127
y útil de los modelos matemáticos de muchos sistemas naturales de
ciencias como la física, la astrofísica, la química y la biología, etc. Las
simulaciones por computadora abarcan desde programas informáticos
cuya ejecución dura unos minutos hasta conjuntos de ordenadores
conectados en red cuya ejecución dura horas, e incluso hay
simulaciones que se extienden varios días.
128
2.5. Hipótesis
2.6. Variables
2.6.1. Variable 01
- Geomorfológicas.
- Hidrológicas.
- Geodinámicos.
2.6.2. Variable 02
- Profundidad de sedimento.
129
2.7. Operacionalización de variables
Definición Definición
Dimensiones Indicadores Escala
conceptual operacional
Variable 01
Parámetros de la
subcuenca
Concentración
Geomorfológicos volumétrica
(Cv)
Superficie
(Villon, 2002) Numero de
terrestre en
Afirma que es curva (CN)
donde las
toda área que Precipitaciones
precipitaciones
Subcuenca desarrolla su
que caen sobre Máximas P24 Ordinal
del río Sicra drenaje (mm)
ella tienden a Hidrológicos
directamente al
ser drenadas Caudal Máximo
curso principal
por el sistema (m3/s)
de la cuenca
de corrientes.
Volumen de
deslizamientos
Geodinámicos (m3)
Áreas con
erosión (m2)
Variable 02
Cantidad de
sedimento Metros cúbicos
transportado
Profundidad de
Metros
sedimento
(Rodríguez,
2017) Define Magnitud de la
Metros por
que los flujos Son mezclas velocidad del
segundo
de lodo y de partículas sedimento
escombros son gruesas y finas Fuerzas de
Flujos de eventos de con agua Newton por
impacto del
lodo y crecidas con donde metro Ordinal
sedimento
escombros concentraciones predominan
de sedimentos arenas, gravas, Presiones
Newton por
estáticas del
que superan el cantos y metro
sedimento
20% en bloques.
concentración Energías
volumétrica. especificas del Metros
sedimento
Kilómetros
Zonas de peligro
cuadrados,
(mapas de
metros
peligro)
cuadrados
130
CAPÍTULO III
MATERIALES Y MÉTODOS
131
Fuente: Adaptado de (Imagen Satelital Landsat, 2015)
3.4.1. Población
3.4.2. Muestra
Figura 44. Vista satelital del cauce del río Sicra – Prog. 0+565 km al 1+744 km
132
Fuente: Adaptado de (Google earth, 2019)
133
búsqueda y la recopilación de información bibliográfica, luego el trabajo en
campo, finalmente el procesamiento de datos en gabinete.
134
CAPÍTULO IV
DISCUSIÓN DE RESULTADOS
[Link]. Generalidades
136
Fuente: Elaboración propia
137
Tabla 16. Cuadro resumen de las características morfológicas de la subcuenca Sicra
138
4 3657 3774 3715.5 9.70 621.35 98.55 1.54
5 3775 3893 3834 15.69 611.65 97.01 2.49
6 3894 4012 3953 24.32 595.96 94.52 3.86
7 4013 4130 4071.5 29.96 571.63 90.66 4.75
8 4131 4249 4190 53.20 541.67 85.91 8.44
9 4250 4368 4309 80.16 488.47 77.47 12.71
10 4369 4486 4427.5 121.40 408.31 64.76 19.25
11 4487 4605 4546 128.57 286.91 45.50 20.39
12 4606 4724 4665 97.88 158.35 25.11 15.52
13 4725 4842 4783.5 42.07 60.47 9.59 6.67
14 4843 4961 4902 15.16 18.40 2.92 2.41
15 4962 5080 5021 2.95 3.23 0.51 0.47
16 5081 5199 5140 0.28 0.28 0.04 0.04
Total = 630.52
Figura 46. Curva hipsométrica y frecuencia de altitudes - subcuenca del río Sicra,
(cuenca en fase de juventud)
15.52%
4550 20.39%
19.25%
4350 12.71%
4150 8.44%
4.75%
3950 3.86%
2.49%
3750 1.54%
0.78%
3550 0.44%
0.23%
3350
0 100 200 300 400 500 600
139
Hydrological Stations) debido a que solamente se cuenta con una sola estación
meteorológica en el distrito de Lircay y provincia de Angaraes, las estaciones
cercanas se encuentran alejadas para poder hacer uso de métodos tradicionales
y poder determinar una precipitación real que represente a todo la subcuenca
del río Sicra.
140
Fuente: Elaboración propia
Tabla 18. Resumen de áreas y perímetros de las microcuencas de la subcuenca del río Sicra
141
Fuente: Elaboración propia
Figura 48. Estaciones climáticas virtuales ubicadas en los centroides de las microcuencas de
la subcuenca del río Sicra
142
Fuente: Elaboración propia
Tabla 19. Ubicación de las estaciones climáticas virtuales empleadas en este estudio
Centroide Centroide
Estación
(Coordenadas (Coordenadas
climática
UTM) geográficas)
virtual
X_Centroide Y_Centroide Latitud Longitud
Microcuenca 01 532686.52 8558177.75 -13.04 -74.7
Microcuenca 02 534106.69 8553527.77 -13.09 -74.68
Microcuenca 03 542575.41 8550041.13 -13.12 -74.61
Microcuenca 04 531447.82 8547179.73 -13.14 -74.71
Microcuenca 05 539032.27 8543454.98 -13.18 -74.64
Microcuenca 06 536379.59 8536794.41 -13.24 -74.67
Microcuenca 07 532278.41 8537443.08 -13.23 -74.7
Microcuenca 08 526288.13 8535715.69 -13.25 -74.76
Microcuenca 09 520689.16 8541850.92 -13.19 -74.81
Tabla 20. Prueba de ajuste de bondad para datos de precipitación de cada microcuenca de
la subcuenca del río Sicra
Delta Delta
Microcuencas Distribución
teórico tabular
Microcuenca 01 Gamma 3P 0.0407 0.2267
143
Microcuenca 02 Log normal 2P 0.0629 0.2267
Microcuenca 03 Log gumbel 0.0877 0.2267
Microcuenca 04 Log normal 2P 0.0759 0.2267
Microcuenca 05 Log gumbel 0.0877 0.2267
Microcuenca 06 Log gumbel 0.0591 0.2267
Microcuenca 07 Gamma 3P 0.0820 0.2267
Microcuenca 08 Gamma 3P 0.0820 0.2267
Microcuenca 09 Log normal 3P 0.0549 0.2267
Tabla 21. Precipitaciones máximas en 24 horas para diferentes períodos de retorno de cada
microcuenca de la subcuenca del río Sicra
144
[Link]. Validación del producto PISCO
ESTACIÓN LIRCAY
Código Latitud Longitud
Lircay – 000657 12°58’55’’ 74°43’5.1’’
Fuente: [Link]
Figura 49. Validación por diferentes medidas bondad de ajuste del producto Pisco con los
datos registrados en la estación meteorológica Lircay – 000627
145
Fuente: Elaboración propia
146
[Link]. Celdas de precipitación concentrada
147
concentración, se ha optado por utilizar la ecuación de Kirpich con el método
de California Culverts Pratice, 1942.
Tabla 24. Tiempo de concentración para cada microcuenca de la subcuenca del río Sicra
Tabla 25. Cálculo del tiempo de retardo para cada microcuenca de la subcuenca del río
Sicra
Tiempo de
Microcuencas
retardo
148
Microcuenca 08 67.54
Microcuenca 09 70.41
149
GRUPO
HIDROLÓGICO
DESCRIPCIÓN DEL USO DE LA TIERRA DEL SUELO
A B C D
Sin tratamiento de
72 81 88 91
conservación.
Tierra cultivada (1):
Con tratamiento de
62 71 78 81
conservación.
Condiciones pobres. 68 79 86 89
Pastizales
Condiciones óptimas. 39 61 74 80
Vegas de río Condiciones óptimas. 30 58 71 78
Troncos delgados,
45 66 77 83
Bosques cubierta pobre, sin hierba.
Parqueaderos
pavimentados, techos, 98 98 98 98
accesos, etc.
150
4. Las áreas permeables restantes (césped) se consideran como pastizales en
buena condición para estos números de curva.
5. En algunos países con climas más cálidos se puede utilizar 95 como número.
de curva.
Fuente: Adaptado de (Ven te Chow , 1994)
151
Fuente: Adaptado de Geo GPS Perú; [Link]
[Link]
Figura 51. Mapa temático de número de curva para la subcuenca del río Sicra
152
Finalmente se obtuvo el CN, para la subcuenca del río Sicra Tabla 27,
que será utilizado para estimar la precipitación neta o en exceso, a partir de
datos de precipitación de cada microcuenca, una vez que se obtiene como
dato la precipitación neta, esta es la que escurre en la superficie formando un
caudal.
Tabla 27. Uso actual de los suelos para cada microcuenca de la subcuenca del río Sicra
153
Figura 52. Modelo de la subcuenca del río Sicra
154
[Link].1. Datos de ingreso al modelo hidrológico HEC-HMS
Tabla 29. Parámetros del método onda cinemática del cauce del río Sicra
155
[Link].2. Simulación del modelo hidrológico HEC-HMS
Figura 53. Hidrograma líquido para la subcuenca del río Sicra (Tr=50 años y Tr=100
años)
156
coeficiente de Manning calibrado para dicho caudal, el cual se observa en la
Tabla 30 y Figura 54.
Figura 54. Hidrograma líquido calibrado para la subcuenca del río Sicra (Tr=50 años y
Tr=100 años)
Los caudales máximos calibrados obtenidos son: 225.00 m3/s y 349.30 m3/s
para un Tr=50 años y Tr=100 años respectivamente, ocurrido el 25 de enero
del 2017 a horas 14:00 horas.
157
observar fuertes pendientes presentando un relieve escarpado y abrupto
propiciando un flujo torrentoso y altamente turbulento, principalmente durante
el período de lluvias torrenciales, por la disminución brusca de la pendiente se
ha formado un cono de deyección, producto de la deposición del material solido
transportando por las aguas del río Sicra principalmente durante el período de
grandes avenidas.
158
Su origen se debe a varios factores donde el principal factor, para la
ocurrencia de la inundación de viviendas y erosión fluvial de vías de
comunicación, son las lluvias torrenciales de temporada, principalmente en el
cauce del río Sicra, generándose ruptura y/o colapso de muros rústicos de
contención. La afectación de viviendas por efecto de las lluvias continuas,
con intervalos de duración de 12, 24 o 48 horas y el volumen de las
precipitaciones que, por lo general se concentraba en unas 5 a 8 horas diarias,
tal como se detalla los antecedentes en las siguientes figuras.
Figura 56. Vista de zonas expuesta a inundación por las aguas del río Sicra
159
2.85 m.
Figura 57. Vista de zonas expuesta a inundación por las aguas del río Sicra
160
Área inundable N° 04:
Figura 59. Área de inundación margen derecha e izquierda río Sicra, zona urbano marginal
ciudad de Lircay, altura puente carrozable Santa Rosa aguas arriba
161
Fuente: Elaboración propia
Figura 60. Vista del puente Santa Rosa, en momentos del incremento del caudal de las
aguas del río Sicra, en proceso de erosión de las bases y estribos
162
Fuente: Elaboración propia
Figura 61. Erosión de talud y pérdida parcial de plataforma, altura puente Chahuarma del
río Sicra
Figura 62. Erosión y pérdida de talud de terrenos forestales y pastizales lado izquierdo del
cauce y trayectoria del río Sicra
163
Fuente: Elaboración propia
[Link]. Deslizamiento:
164
identificado 6 puntos activos a deslizamientos y huaycos, los cuales se
resumen en la Tabla 32 siguiente.
Zona
Progresiva Coordenadas
vulnerable a
(km) UTM
deslizamiento
Este Norte
N°
(m) (m)
01 4+864 530,946 8,560,338
02 4+896 531,038 8,560,286
03 4+885 530,162 8,562,783
04 9+904 534,200 8,558,000
05 10+063 534,287 8,558,295
06 10+190 534,388 8,558,523
Figura 63. Vista de punto activo a deslizamientos N° 01 (D1) sobre el tramo alterno en la
vía Lircay – Ayacucho
165
Progresiva: 4+896 km (UTM: 531038 E, 8560286 S).
Figura 64. Vista de punto activo a deslizamientos N° 02 (D2) sobre el tramo alterno en la
vía Lircay – Ayacucho
Figura 65. Vista de punto activo a deslizamientos N° 03 (D3) sobre el tramo alterno en la
vía Lircay – Ayacucho
166
Fuente: Elaboración propia
Figura 66. Vista de punto activo a deslizamientos N° 04 (D4) sobre el tramo alterno en la
vía Lircay – Ayacucho
167
Fuente: Elaboración propia
Figura 67. Vista de punto activo a deslizamientos N° 05 (D5) sobre el tramo alterno en la
vía Lircay – Ayacucho
Figura 68. Vista de punto activo a deslizamientos N° 06 (D6) sobre el tramo alterno en la
vía Lircay – Ayacucho
168
Fuente: Elaboración propia
Provincia de
Angaraes,
Lircay
169
Fuente: Instituto Geológico, Minero y Metalúrgico “INGEMMET”.
Figura 70. Mapa de suelos con recurrencia de eventos de caídas, deslizamientos y flujos de
masas, en la provincia de Angaraes, Lircay
Provincia de
Angaraes, Lircay
Figura 71. Mapa de peligro por movimientos en masa para el verano 2020
170
Fuente: Adaptado de CENEPRED
171
Fuente: Adaptado de CENEPRED
Tabla 32. Grado de susceptibilidad del suelo según Instituto Geológico y Metalúrgico
172
Fuente: Instituto Geológico, Minero y Metalúrgico “INGEMMET”
GRADO/CRITERÍO
Alta: Laderas que tiene zonas de falla, masas de roca con meteorización alta
a moderada, fracturadas con discontinuidades desfavorables; deposito
superficiales inconsolidados, materiales parcialmente a muy saturados,
laderas con pendientes entre 25° a 45°, donde han ocurrido movimientos en
masa o existe la posibilidad de que ocurra.
173
4.1.4. Pre - procesamiento de flujos de lodo y escombros con el modelo numérico
Flo - 2D
[Link].1. Topografía
174
- 2D, consistió en la construcción de un modelo digital de elevaciones
(DEM) que fue generado a partir de vuelo de un drone (Phantom 4 RTK),
del cual se obteniene un DEM de alta resolución de 0.5mx0.5m, una imagen
ortofoto de alta resolución 4K, 0.1mx0.1m.
Figura 74. Modelo digital de elevaciones de la subcuenca del río Sicra, del distrito de
Lircay - Angaraes
175
Este modelo fue ingresado al software Flo - 2D para realizar el
grillado que consiste en dividir el modelo digital de elevaciones (DEM), en
grillas en el cual el modelo calculará terminada la simulación, los tirantes y
velocidades del flujo de lodo y escombros, tal como se puede observar en la
Figura 75.
Figura 75. Grillado del modelo digital de elevaciones del distrito de Lircay, de 5mx5m
176
estudios O’Brien y Julien (1983) recomienda que la concentración
volumétrica empiece con un valor de concentración cercano a 0.22 y luego
aumentar gradualmente hasta llegar a su máximo valor, además este valor
máximo debe ubicarse minutos antes del pico del hidrograma líquido.
Figura 76. Hidrograma de caudal líquido y concentración de sólidos para Tr=50 años
Q
(m3/s) 0.45
200
0.4
Caudal (m3/s)
150
Cv
0.35
100
0.3
50
0.25
0 0.2
[Link]
Tiempo (hr)
Figura 77. Hidrograma de caudal líquido y concentración de sólidos para Tr=100 años
177
Caudal líquido y Concentración de solido
Tr = 100 años
400 0.5
350 Q
(m3/s) 0.45
300
0.4
Caudal (m3/s)
250
Cv
200 0.35
150
0.3
100
0.25
50
0 0.2
[Link]
Tiempo (hr)
Tabla 33. Resultados del análisis granulométrico por tamizado (ASTM D-422)
Caracterización de estratos
Límite Índice
Muestra Humedad Clasificación liquido plástico
(%) (SUCS)
178
Arcilla (%)
D200
En el cauce principal
CL-01 0.21 - SW-SP NP NP
CL-02 0.18 - SW-SP NP NP
En el abanico aluvial
CL-03 0.23 8.19 GW-GP 0.26 0.062
CL-04 0.31 10.45 GW-GP 0.27 0.087
Tabla 34. Selección de matriz de flujos de lodo (limo y arcilla), de propiedades similares
Aspen
natural 27 0.001 0.012 0.028 0.25 0.06
soil
179
Tabla 35. Resultados de los parámetros reológicos de esfuerzo de cedencia y viscosidad
dinámica
𝝉𝒚 = 𝜶𝒆𝜷𝐂𝐯 𝜼 = 𝜶𝒆𝜷𝐂𝐯
𝜶 𝜷 𝜶 𝜷
0.152 18.7 0.00136 28.4
180
Tabla 36. Valores de cobertura de usos de suelos representativos en cauce y zonas
urbanas
181
las grillas para la simulación. En el ambiente del pre-procesador (GDS) una
vez obtenidos las grillas y delimitado el dominio computacional se procede
asignar los coeficientes de Manning y la ubicación del hidrograma de
entrada (Inflow) en la grilla del dominio computacional. Luego de definir,
asignar los valores de Manning, también se define las grillas de salidas de
flujo (Outflow).
[Link]
[Link]
[Link]
[Link]
[Link]
[Link]
[Link]
182
nuestro caso el proceso a ser simulado es Mud/Debris o flujo de lodo
y escombros.
183
Finalmente, los resultados son graficados en los post-procesadores
MAPPER ++ y MAPPER Pro pertenecientes al programa Flo - 2D.
Los resultados obtenidos por el programa son las elevaciones del
terreno, tirantes, velocidades, presión de impacto, volúmenes de agua
y sedimento del flujo.
Se creó una grilla con celdas de 5x5 metros, teniendo un total 1,446
celdas para el área urbana de la subcuenca del río Sicra.
184
pico, aumenta demasiado rápido con respecto al tiempo, por lo que la
simulación se hace más lenta.
Figura 78. Ventana de Flo - 2D, simulación del flujo y desarrolla del hidrograma en el
tiempo
Bulked w/Sediment
Flujo (m3) Agua (m3)
(𝒎𝟑 )
Hidrograma de entrada (Inflow) 4’235,896.07 5’737,539.54
Almacenamiento dentro del área de
5,679.73 7,473.18
análisis
Flujo fuera del área de análisis
4’230,216.34 5’730,066.35
(Outflow)
185
Fuente: Elaboración Propia
Tabla 39. Salidas del programa Flo - 2D (Q = 349.30 m3/s - Tr = 100 años)
Bulked w/Sediment
Flujo (m3) Agua (m3)
((𝒎𝟑 )
Hidrograma de entrada (Inflow) 5’927,259.09 8’113,853.35
Almacenamiento dentro del área de
6,223.72 8,278.04
análisis
Flujo fuera del área de análisis
5’921,035.37 8’105,575.31
(Outflow)
186
La herramienta esencial del MAPPER PRO es la generación de mapas
de peligro, en nuestro caso fue usado para generar mapas de peligros digitales
según los parámetros del escenario simulados,
187
188
189
De donde se observa zonas de peligro:
190
La validación es el proceso de establecer evidencia que proporciona
información razonable de que la salida del software es lo suficientemente
buena como para representar los hechos en el terreno al nivel aceptable. La
validación del modelo de este estudio empleó la recolección de puntos de
control de tierra (GCP) en zonas de deslizamientos activos, evidencias por
erosión fluvial en zona urbana de Lircay. Ambos se utilizaron para identificar
la inundación máxima que pudiera causar en eventos de precipitaciones
máximas. Debido a evidencias ocasionadas en los últimos años por inundación
se ha optado realizar el modelamiento en tramo crítico del cauce del río Sicra,
de la zona urbana del distrito de Lircay, progresiva 0+565 km al 1+744 km.
Figura 79. Identificación de puntos de control en zonas críticos de suelos inestables con alta
presencia de grava suelto
191
Fuente: Elaboración Propia
Figura 81. Identificación de puntos de control en zonas críticos caída de rocas hacia el cauce
del río Sicra
192
Fuente: Elaboración Propia
Figura 82. Identificación de puntos de control, Caída de rocas mayores 2.00 metros de
longitud, al cauce del río Sicra
193
Fuente: Elaboración Propia
Figura 83. Identificación de puntos de control en zonas críticos y con mayor vulnerabilidad a
deslizamientos, caída de rocas y desprendimiento de escombros hacia el cauce principal del
río Sicra
194
Fuente: Elaboración Propia
195
Fuente: Elaboración Propia
Tabla 40. Resultado de volumen de sedimento depositado en el cauce del río Sicra estimados
de campo
Zonas Volumen
Fid* Área (m2) Code
Críticos (m3)
17 Zona 01 341,157.42 41,736.76 -1
21 Zona 02 549,154.46 80,935.23 -1
56 Zona 03 1’549,154.46 228,316.76 -1
48 Zona 04 854,689.23 81,750.95 -1
78 Zona 05 1’879,254.25 276,967.37 -1
23 Zona 06 978,157.42 138,709.22 -1
Volumen Total = 6’151,567.24 848,416.28
196
HUANCAVELICA” y dar cumplimiento a mi objetivo general que es el de
“Definir zonas de peligro de flujos de lodo y escombros generados por
parámetros geomorfológicos, hidrológicos y geodinámicas en la subcuenca del
rio Sicra, Lircay – Huancavelica”, se alcanza los siguientes resultados:
197
caudales, que después de ocurrido los deslizamientos y depositados en el cauce
del rio, estas tienen alto grado de transporte de estos materiales (flujos de lodo
y escombros).
198
CONCLUSIONES
199
RECOMENDACIONES
200
REFERENCIAS BIBLIOGRÁFICAS
Aguirre Pe, J., & T. Moncada , A. (1999). Regimenes en el flujo de barros y escombros.
Revista cientifica, 8-8.
Aguirre, J., Moncada, A., & Olivero, M. (2000). Flujo y deposición de barros y
escombros. Venezuela: Universidad de los Andes.
Bartali, R., Sarocchi, D., Nahmad, Y., & Rodríguez, A. (2012). Estudio de flujos
granulares de tipo geológico por medio de simulador multisensor
GRANFLOW-SIM. México: Boletín de la sociedad Geológica Mexicana.
Bateman, A., Medina, V., Hurlimann, M., & Velasco, D. (2017). Modelo
bidimensional para simulación de flujos detríticos FLATmodel. Aplicación a
una cuenca del Pirineo Catalán. Barcelona: Universidad Politécnica de
Cataluña.
Breña Puyol, A., & Jacobo Villa, M. (2006). Principios y fundamentos de hidrología
superficial. Tlalpan, D.F. México: Rectoría General de la Universidad -
México.
Breña Puyol, A., & Jacovo Villa, M. (10 de Enero de 2006). [Link].
Obtenido de Principios y Fundaamnetos de Hidrlogía Superficial:
201
[Link]
_fundamentos/[Link]
Build a free website of your own on. (14 de Noviembre de 2018). TRIPOD. Obtenido
de [Link]
Cahuana A., A., & Yugar M., W. (01 de septiembre de 2009). Scribd. Obtenido de
Material de apoyo didáctico para la enseñanza y aprendizaje de la asignatura
de hidrología civ-233: [Link]
30-09-09-Adscripicion-
MaterialDeApoyoDidacticoParaLaEnsenanzaYAprendizajeDeLaAsignatura-
pdf
Cahuana Andia, A., & Yugar Morales, W. (2009). Materialde apoyo didactico para la
enseñanza y aprendizaje de la asignatura de la hidrología civ-233. Bolivia .
Campos Aranda, D. (1998). Procesos del ciclo hidrológico. San luis Potosí S.L.P.
Mexico,: Universitaria Potosina.
Eder, M., Barbeito, O., Ambrosino, S., Hillman, G., Moya, G., & Rodriguez, A.
(2016). Análisis y modelo hidráulico (2D) de flujos hiperconcentrados en
piedemonde occidental de la sierra de ambato, catamarca, argentina.
SCIPEDIA, 65-71.
Escobar Soto, P., & Poma Yaranga, V. (2018). “Analisis de amenaza por el transporte
de detritos aplicado a la quebrada de Tantará del distrito de Tantará,
provincia de Castrovirreyna de la región de Huancavelica”. Huancavelica:
Universidad Nacional De Huancavelica.
Gamion, J. (2014). Modelo de flujo de escombros y lodo aplicando FLO-2D, caso sub
cuenca del rio Shullcas. Huancayo: Universidad Nacional del Centro del Perú.
Hungr, O. (2000). Analysis of Debris Flow Surges Using the Theory of Uniformly
Progressive Flow. Earth Surface Processes and Landforms.
Ibáñez Asensio, S., Moreno Ramón, H., & Gisbert Blanquer, J. M. (2011). Morfología
de las cuencas hidrológicas. riunet, 7-9. Obtenido de Morfología de las
cuencas.
203
Iparraguirre Ortiz, C. (01 de 10 de 2019). SCRIBD. Obtenido de Flujo De Escombros:
Aspectos Hidrologicos E Hidraulicos:
[Link]
Jakob, M., & Hungr, O. (2005). Debris-Flow Hazards and Related Phenomena
(Peligros de flujo de escombros y fenómenos relacionados). 2329 West Mall,
Vancouver, BC V6T 1Z4, Canadá.
Luque Pinto, R., & Alva Hurtado, J. (1999). Estimación y Evalución de Flujo de
Escombros. [Link], 2-2.
O´Brien, J., & Julien, P. (1997). On the importance of mudflow routing. Proceedings
of the first international Conference on Debris Flow Hazards
Mitigation,ASCE.
O'Brien, J., & Julien, P. (1988). Laboratory Analysis of Mudflow Properties. ournal
of Hydraulic Engineering,, 877-887.
204
Oliveira, J. (2012). Estudio del impacto de los flujos de escombros en la cuenca del
río San Antonio-Brasil. Brasil.
S Ibañez, A., & H Moreno, R. (2011). “El uso de las fotografías aéreas en los estudios
de suelos”. Riunet de la UPV, 3-5.
Sepúlveda, A., Patiño, J., & Rodriguez, C. (2016). Metodología para evaluación de
riesgo por flujo de detritos detonados por lluvias caso Útica-Colombia.
Bogotá: Pontificia Universidad Javeriana, Departamento de Ingeniería Civil.
205
Takahashi, T. (2009). A Review of Japanese Debris Flow Research. International
Journal of Erosion Control Engineering,, 1-14.
206
APÉNDICE
207
ANEXO A – HIDROLOGÍA
208
Figura A.1.3. Curva IDF de la microcuenca 03
209
Figura A.1.5. Curva IDF de la microcuenca 05
210
Figura A.1.7. Curva IDF de la microcuenca 07
211
Figura A.1.9. Curva IDF de la microcuenca 09
212
A.2 Hietogramas de las microcuencas de la subcuenca del río Sicra
213
Figura A.2.2. Hietogramas de la microcuenca 02 para un Tr = 50 años y Tr = 100 años
214
Figura A.2.3. Hietogramas de la microcuenca 03 para un Tr = 50 años y Tr = 100 años.
215
Figura A.2.4. Hietogramas de la microcuenca 04 para un Tr = 50 años y Tr = 100 años
216
Figura A.2.5. Hietogramas de la microcuenca 05 para un Tr = 50 años y Tr = 100 años
217
Figura A.2.6. Hietogramas de la microcuenca 06 para un Tr = 50 años y Tr = 100 años
218
Figura A.2.7. Hietogramas de la microcuenca 07 para un Tr = 50 años y Tr = 100 años
219
Figura A.2.8. Hietogramas de la microcuenca 08 para un Tr = 50 años y Tr = 100 años
220
Figura A.2.9. Hietogramas de la microcuenca 09 para un Tr = 50 años y Tr = 100 años
221
A.3 Hidrograma líquido con el modelo hidrológico HEC-HMS en la subcuenca
del río Sicra
222
25-Ene-17 12:50 10.9 25-Ene-17 20:20 56.7
25-Ene-17 13:00 12 25-Ene-17 20:30 55.6
25-Ene-17 13:10 15.4 25-Ene-17 20:40 54.5
25-Ene-17 13:20 25.6 25-Ene-17 20:50 53.5
25-Ene-17 13:30 47.1 25-Ene-17 21:00 52.6
25-Ene-17 13:40 88.8 25-Ene-17 21:10 51.7
25-Ene-17 13:50 151.7 25-Ene-17 21:20 50.8
25-Ene-17 14:00 187.9 25-Ene-17 21:30 50
25-Ene-17 14:10 208 25-Ene-17 21:40 49.2
25-Ene-17 14:20 221 25-Ene-17 21:50 48.5
25-Ene-17 14:30 225 25-Ene-17 22:00 47.7
25-Ene-17 14:40 220.3 25-Ene-17 22:10 47
25-Ene-17 14:50 210.8 25-Ene-17 22:20 46.4
25-Ene-17 15:00 199.6 25-Ene-17 22:30 45.7
25-Ene-17 15:10 188.3 25-Ene-17 22:40 45.1
25-Ene-17 15:20 177.6 25-Ene-17 22:50 44.5
25-Ene-17 15:30 167.4 25-Ene-17 23:00 43.9
25-Ene-17 15:40 158 25-Ene-17 23:10 43.3
25-Ene-17 15:50 149.1 25-Ene-17 23:20 42.8
25-Ene-17 16:00 140.9 25-Ene-17 23:30 42.3
25-Ene-17 16:10 133.4 25-Ene-17 23:40 41.8
25-Ene-17 16:20 126.5 25-Ene-17 23:50 41.3
25-Ene-17 16:30 120.2 26-Ene-17 00:00 40.8
25-Ene-17 16:40 114.4 26-Ene-17 00:10 40.3
25-Ene-17 16:50 109.1 26-Ene-17 00:20 39.8
25-Ene-17 17:00 104.2 26-Ene-17 00:30 39.2
25-Ene-17 17:10 99.8 26-Ene-17 00:40 38.6
25-Ene-17 17:20 95.7 26-Ene-17 00:50 37.9
25-Ene-17 17:30 91.9 26-Ene-17 01:00 37.1
25-Ene-17 17:40 88.5 26-Ene-17 01:10 36.4
25-Ene-17 17:50 85.3 26-Ene-17 01:20 35.5
25-Ene-17 18:00 82.3 26-Ene-17 01:30 34.6
25-Ene-17 18:10 79.6 26-Ene-17 01:40 33.6
25-Ene-17 18:20 77 26-Ene-17 01:50 32.5
25-Ene-17 18:30 74.6 26-Ene-17 02:00 31.3
25-Ene-17 18:40 72.4 26-Ene-17 02:10 30.1
25-Ene-17 18:50 70.4 26-Ene-17 02:20 28.7
25-Ene-17 19:00 68.5 26-Ene-17 02:30 27.3
25-Ene-17 19:10 66.7 26-Ene-17 02:40 25.8
25-Ene-17 19:20 65 26-Ene-17 02:50 24.3
25-Ene-17 19:30 63.4 26-Ene-17 03:00 22.9
25-Ene-17 19:40 61.9 26-Ene-17 03:10 21.4
25-Ene-17 19:50 60.5 26-Ene-17 03:20 20.1
25-Ene-17 20:00 59.2 26-Ene-17 03:30 18.8
25-Ene-17 20:10 57.9 26-Ene-17 03:40 17.5
223
26-Ene-17 03:50 16.4 26-Ene-17 11:20 1.6
26-Ene-17 04:00 15.3 26-Ene-17 11:30 1.5
26-Ene-17 04:10 14.3 26-Ene-17 11:40 1.5
26-Ene-17 04:20 13.3 26-Ene-17 11:50 1.4
26-Ene-17 04:30 12.5 26-Ene-17 12:00 1.4
26-Ene-17 04:40 11.6 26-Ene-17 12:10 1.3
26-Ene-17 04:50 10.9 26-Ene-17 12:20 1.3
26-Ene-17 05:00 10.2 26-Ene-17 12:30 1.2
26-Ene-17 05:10 9.6 26-Ene-17 12:40 1.2
26-Ene-17 05:20 9 26-Ene-17 12:50 1.1
26-Ene-17 05:30 8.4 26-Ene-17 13:00 1.1
26-Ene-17 05:40 7.9 26-Ene-17 13:10 1.1
26-Ene-17 05:50 7.4 26-Ene-17 13:20 1
26-Ene-17 06:00 7 26-Ene-17 13:30 1
26-Ene-17 06:10 6.6 26-Ene-17 13:40 1
26-Ene-17 06:20 6.2 26-Ene-17 13:50 1
26-Ene-17 06:30 5.9 26-Ene-17 14:00 0.9
26-Ene-17 06:40 5.5 26-Ene-17 14:10 0.9
26-Ene-17 06:50 5.2 26-Ene-17 14:20 0.9
26-Ene-17 07:00 5 26-Ene-17 14:30 0.8
26-Ene-17 07:10 4.7 26-Ene-17 14:40 0.8
26-Ene-17 07:20 4.4 26-Ene-17 14:50 0.8
26-Ene-17 07:30 4.2 26-Ene-17 15:00 0.8
26-Ene-17 07:40 4 26-Ene-17 15:10 0.8
26-Ene-17 07:50 3.8 26-Ene-17 15:20 0.7
26-Ene-17 08:00 3.6 26-Ene-17 15:30 0.7
26-Ene-17 08:10 3.5 26-Ene-17 15:40 0.7
26-Ene-17 08:20 3.3 26-Ene-17 15:50 0.7
26-Ene-17 08:30 3.1 26-Ene-17 16:00 0.7
26-Ene-17 08:40 3 26-Ene-17 16:10 0.6
26-Ene-17 08:50 2.9 26-Ene-17 16:20 0.6
26-Ene-17 09:00 2.7 26-Ene-17 16:30 0.6
26-Ene-17 09:10 2.6 26-Ene-17 16:40 0.6
26-Ene-17 09:20 2.5 26-Ene-17 16:50 0.6
26-Ene-17 09:30 2.4 26-Ene-17 17:00 0.6
26-Ene-17 09:40 2.3 26-Ene-17 17:10 0.6
26-Ene-17 09:50 2.2 26-Ene-17 17:20 0.5
26-Ene-17 10:00 2.1 26-Ene-17 17:30 0.5
26-Ene-17 10:10 2 26-Ene-17 17:40 0.5
26-Ene-17 10:20 2 26-Ene-17 17:50 0.5
26-Ene-17 10:30 1.9 26-Ene-17 18:00 0.5
26-Ene-17 10:40 1.8 26-Ene-17 18:10 0.5
26-Ene-17 10:50 1.7 26-Ene-17 18:20 0.5
26-Ene-17 11:00 1.7 26-Ene-17 18:30 0.5
26-Ene-17 11:10 1.6 26-Ene-17 18:40 0.5
224
26-Ene-17 18:50 0.4 26-Ene-17 22:40 0.3
26-Ene-17 19:00 0.4 26-Ene-17 22:50 0.3
26-Ene-17 19:10 0.4 26-Ene-17 23:00 0.3
26-Ene-17 19:20 0.4 26-Ene-17 23:10 0.3
26-Ene-17 19:30 0.4 26-Ene-17 23:20 0.3
26-Ene-17 19:40 0.4 26-Ene-17 23:30 0.3
26-Ene-17 19:50 0.4 26-Ene-17 23:40 0.3
26-Ene-17 20:00 0.4 26-Ene-17 23:50 0.2
26-Ene-17 20:10 0.4 27-Ene-17 00:00 0.2
26-Ene-17 20:20 0.4 27-Ene-17 00:10 0.2
26-Ene-17 20:30 0.4 27-Ene-17 00:20 0.2
26-Ene-17 20:40 0.4 27-Ene-17 00:30 0.2
26-Ene-17 20:50 0.3 27-Ene-17 00:40 0.2
26-Ene-17 21:00 0.3 27-Ene-17 00:50 0.2
26-Ene-17 21:10 0.3 27-Ene-17 01:00 0.2
26-Ene-17 21:20 0.3 27-Ene-17 01:10 0.2
26-Ene-17 21:30 0.3 27-Ene-17 01:20 0.2
26-Ene-17 21:40 0.3 27-Ene-17 01:30 0.2
26-Ene-17 21:50 0.3 27-Ene-17 01:40 0.2
26-Ene-17 22:00 0.3 27-Ene-17 01:50 0.2
26-Ene-17 22:10 0.3 27-Ene-17 02:00 0
26-Ene-17 22:20 0.3
26-Ene-17 22:30 0.3
TIME SERIES RESULT Tr = 100 YEAR
225
25-Ene-17 06:10 0 25-Ene-17 13:40 279.8
25-Ene-17 06:20 0 25-Ene-17 13:50 314.8
25-Ene-17 06:30 0 25-Ene-17 14:00 338.2
25-Ene-17 06:40 0 25-Ene-17 14:10 349.3
25-Ene-17 06:50 0 25-Ene-17 14:20 345.1
25-Ene-17 07:00 0 25-Ene-17 14:30 330.8
25-Ene-17 07:10 0 25-Ene-17 14:40 312
25-Ene-17 07:20 0 25-Ene-17 14:50 292.2
25-Ene-17 07:30 0 25-Ene-17 15:00 273
25-Ene-17 07:40 0 25-Ene-17 15:10 255.4
25-Ene-17 07:50 0 25-Ene-17 15:20 239.2
25-Ene-17 08:00 0 25-Ene-17 15:30 224.4
25-Ene-17 08:10 0 25-Ene-17 15:40 210.6
25-Ene-17 08:20 0 25-Ene-17 15:50 198
25-Ene-17 08:30 0 25-Ene-17 16:00 186.4
25-Ene-17 08:40 0 25-Ene-17 16:10 175.8
25-Ene-17 08:50 0 25-Ene-17 16:20 166.2
25-Ene-17 09:00 0 25-Ene-17 16:30 157.5
25-Ene-17 09:10 0 25-Ene-17 16:40 149.6
25-Ene-17 09:20 0 25-Ene-17 16:50 142.3
25-Ene-17 09:30 0 25-Ene-17 17:00 135.6
25-Ene-17 09:40 0 25-Ene-17 17:10 129.6
25-Ene-17 09:50 0 25-Ene-17 17:20 124
25-Ene-17 10:00 0 25-Ene-17 17:30 119
25-Ene-17 10:10 0 25-Ene-17 17:40 114.4
25-Ene-17 10:20 0 25-Ene-17 17:50 110.1
25-Ene-17 10:30 0 25-Ene-17 18:00 106.2
25-Ene-17 10:40 0 25-Ene-17 18:10 102.5
25-Ene-17 10:50 0 25-Ene-17 18:20 99.2
25-Ene-17 11:00 0 25-Ene-17 18:30 96
25-Ene-17 11:10 0 25-Ene-17 18:40 93.1
25-Ene-17 11:20 0 25-Ene-17 18:50 90.5
25-Ene-17 11:30 0 25-Ene-17 19:00 87.9
25-Ene-17 11:40 0.2 25-Ene-17 19:10 85.6
25-Ene-17 11:50 0.8 25-Ene-17 19:20 83.4
25-Ene-17 12:00 2.6 25-Ene-17 19:30 81.3
25-Ene-17 12:10 5.5 25-Ene-17 19:40 79.4
25-Ene-17 12:20 9.3 25-Ene-17 19:50 77.6
25-Ene-17 12:30 12.7 25-Ene-17 20:00 75.8
25-Ene-17 12:40 15 25-Ene-17 20:10 74.2
25-Ene-17 12:50 17.4 25-Ene-17 20:20 72.6
25-Ene-17 13:00 25.4 25-Ene-17 20:30 71.1
25-Ene-17 13:10 49 25-Ene-17 20:40 69.7
25-Ene-17 13:20 106.8 25-Ene-17 20:50 68.4
25-Ene-17 13:30 209.4 25-Ene-17 21:00 67.2
226
25-Ene-17 21:10 66 26-Ene-17 04:40 13.2
25-Ene-17 21:20 64.9 26-Ene-17 04:50 12.3
25-Ene-17 21:30 63.8 26-Ene-17 05:00 11.4
25-Ene-17 21:40 62.8 26-Ene-17 05:10 10.7
25-Ene-17 21:50 61.9 26-Ene-17 05:20 10
25-Ene-17 22:00 60.9 26-Ene-17 05:30 9.3
25-Ene-17 22:10 60 26-Ene-17 05:40 8.7
25-Ene-17 22:20 59.2 26-Ene-17 05:50 8.2
25-Ene-17 22:30 58.4 26-Ene-17 06:00 7.7
25-Ene-17 22:40 57.6 26-Ene-17 06:10 7.2
25-Ene-17 22:50 56.8 26-Ene-17 06:20 6.8
25-Ene-17 23:00 56 26-Ene-17 06:30 6.4
25-Ene-17 23:10 55.3 26-Ene-17 06:40 6
25-Ene-17 23:20 54.6 26-Ene-17 06:50 5.7
25-Ene-17 23:30 54 26-Ene-17 07:00 5.3
25-Ene-17 23:40 53.3 26-Ene-17 07:10 5.1
25-Ene-17 23:50 52.7 26-Ene-17 07:20 4.8
26-Ene-17 00:00 52.1 26-Ene-17 07:30 4.5
26-Ene-17 00:10 51.5 26-Ene-17 07:40 4.3
26-Ene-17 00:20 50.8 26-Ene-17 07:50 4.1
26-Ene-17 00:30 50.1 26-Ene-17 08:00 3.9
26-Ene-17 00:40 49.3 26-Ene-17 08:10 3.7
26-Ene-17 00:50 48.4 26-Ene-17 08:20 3.5
26-Ene-17 01:00 47.4 26-Ene-17 08:30 3.3
26-Ene-17 01:10 46.4 26-Ene-17 08:40 3.2
26-Ene-17 01:20 45.2 26-Ene-17 08:50 3
26-Ene-17 01:30 44 26-Ene-17 09:00 2.9
26-Ene-17 01:40 42.6 26-Ene-17 09:10 2.8
26-Ene-17 01:50 41 26-Ene-17 09:20 2.6
26-Ene-17 02:00 39.4 26-Ene-17 09:30 2.5
26-Ene-17 02:10 37.5 26-Ene-17 09:40 2.4
26-Ene-17 02:20 35.6 26-Ene-17 09:50 2.3
26-Ene-17 02:30 33.6 26-Ene-17 10:00 2.2
26-Ene-17 02:40 31.5 26-Ene-17 10:10 2.1
26-Ene-17 02:50 29.5 26-Ene-17 10:20 2
26-Ene-17 03:00 27.5 26-Ene-17 10:30 2
26-Ene-17 03:10 25.6 26-Ene-17 10:40 1.9
26-Ene-17 03:20 23.8 26-Ene-17 10:50 1.8
26-Ene-17 03:30 22.1 26-Ene-17 11:00 1.7
26-Ene-17 03:40 20.5 26-Ene-17 11:10 1.7
26-Ene-17 03:50 19 26-Ene-17 11:20 1.6
26-Ene-17 04:00 17.7 26-Ene-17 11:30 1.6
26-Ene-17 04:10 16.4 26-Ene-17 11:40 1.5
26-Ene-17 04:20 15.2 26-Ene-17 11:50 1.5
26-Ene-17 04:30 14.2 26-Ene-17 12:00 1.4
227
26-Ene-17 12:10 1.4 26-Ene-17 19:10 0.4
26-Ene-17 12:20 1.3 26-Ene-17 19:20 0.4
26-Ene-17 12:30 1.3 26-Ene-17 19:30 0.4
26-Ene-17 12:40 1.2 26-Ene-17 19:40 0.4
26-Ene-17 12:50 1.2 26-Ene-17 19:50 0.4
26-Ene-17 13:00 1.2 26-Ene-17 20:00 0.4
26-Ene-17 13:10 1.1 26-Ene-17 20:10 0.4
26-Ene-17 13:20 1.1 26-Ene-17 20:20 0.4
26-Ene-17 13:30 1 26-Ene-17 20:30 0.4
26-Ene-17 13:40 1 26-Ene-17 20:40 0.4
26-Ene-17 13:50 1 26-Ene-17 20:50 0.4
26-Ene-17 14:00 1 26-Ene-17 21:00 0.3
26-Ene-17 14:10 0.9 26-Ene-17 21:10 0.3
26-Ene-17 14:20 0.9 26-Ene-17 21:20 0.3
26-Ene-17 14:30 0.9 26-Ene-17 21:30 0.3
26-Ene-17 14:40 0.9 26-Ene-17 21:40 0.3
26-Ene-17 14:50 0.8 26-Ene-17 21:50 0.3
26-Ene-17 15:00 0.8 26-Ene-17 22:00 0.3
26-Ene-17 15:10 0.8 26-Ene-17 22:10 0.3
26-Ene-17 15:20 0.8 26-Ene-17 22:20 0.3
26-Ene-17 15:30 0.7 26-Ene-17 22:30 0.3
26-Ene-17 15:40 0.7 26-Ene-17 22:40 0.3
26-Ene-17 15:50 0.7 26-Ene-17 22:50 0.3
26-Ene-17 16:00 0.7 26-Ene-17 23:00 0.3
26-Ene-17 16:10 0.7 26-Ene-17 23:10 0.3
26-Ene-17 16:20 0.6 26-Ene-17 23:20 0.3
26-Ene-17 16:30 0.6 26-Ene-17 23:30 0.3
26-Ene-17 16:40 0.6 26-Ene-17 23:40 0.3
26-Ene-17 16:50 0.6 26-Ene-17 23:50 0.3
26-Ene-17 17:00 0.6 27-Ene-17 00:00 0.2
26-Ene-17 17:10 0.6 27-Ene-17 00:10 0.2
26-Ene-17 17:20 0.6 27-Ene-17 00:20 0.2
26-Ene-17 17:30 0.5 27-Ene-17 00:30 0.2
26-Ene-17 17:40 0.5 27-Ene-17 00:40 0.2
26-Ene-17 17:50 0.5 27-Ene-17 00:50 0.2
26-Ene-17 18:00 0.5 27-Ene-17 01:00 0.2
26-Ene-17 18:10 0.5 27-Ene-17 01:10 0.2
26-Ene-17 18:20 0.5 27-Ene-17 01:20 0.2
26-Ene-17 18:30 0.5 27-Ene-17 01:30 0.2
26-Ene-17 18:40 0.5 27-Ene-17 01:40 0.2
26-Ene-17 18:50 0.5 27-Ene-17 01:50 0.2
26-Ene-17 19:00 0.4 27-Ene-17 02:00 0
228
A.4 Sección transversal del cauce del río Sicra para calibración del modelo
hidrológico HEC-HMS en la subcuenca del río Sicra
229
ANEXO B – EMS
230
ANEXO C – Flo - 2D
[Link] FILE
Pro Model - Build No. 18.12.20
231
40.000 1.819 -0.000001 0.000000
41.000 1.852 -0.000001 0.000000
42.000 1.864 -0.000001 0.000000
43.000 1.888 -0.000001 0.000000
44.000 1.899 -0.000001 0.000000
45.001 1.947 -0.000001 0.000000
46.000 1.947 -0.000001 0.000000
47.001 1.955 -0.000001 0.000000
48.001 1.964 -0.000001 0.000000
49.001 1.970 -0.000001 0.000000
50.000 1.973 -0.000000 0.000000
51.000 1.965 -0.000000 0.000000
52.000 2.001 -0.000001 0.000000
53.000 2.052 -0.000000 0.000000
54.000 2.065 -0.000000 0.000000
55.000 2.091 -0.000000 0.000000
=====================================================================
=============================================================
MASS BALANCE INFLOW - OUTFLOW VOLUME
=====================================================================
=============================================================
*** INFLOW (CU METERS) ***
WATER BULKED W/SEDIMENT
=====================================================================
=============================================================
*** SURFACE OUTFLOW (CU METERS) ***
=====================================================================
=============================================================
*** TOTALS ***
232
(FOR FLOW DEPTHS GREATER THAN THE "TOL" VALUE TYPICALLY 0.1 FT OR 0.03
M)
THE MAXIMUM INUNDATED AREA (DEPTH > 0.15 M) IS: 74325.000 SQ.
METERS
=====================================================================
=============================================================
233
C.2 Simulación de flujos de lodo y escombros para un Tr = 100 años
[Link] FILE
Pro Model - Build No. 18.12.20
234
44.000 1.893 -0.000000 0.000000
45.000 1.942 -0.000000 0.000000
46.000 1.971 -0.000000 0.000000
47.001 1.993 -0.000000 0.000000
48.000 1.987 -0.000000 0.000000
49.000 1.988 -0.000000 0.000000
50.000 2.004 -0.000000 0.000000
51.000 2.005 -0.000000 0.000000
52.000 2.039 -0.000000 0.000000
53.000 2.090 -0.000000 0.000000
54.001 2.105 -0.000000 0.000000
55.000 2.123 -0.000000 0.000000
=====================================================================
=============================================================
MASS BALANCE INFLOW - OUTFLOW VOLUME
=====================================================================
=============================================================
*** INFLOW (CU METERS) ***
WATER BULKED W/SEDIMENT
=====================================================================
=============================================================
*** SURFACE OUTFLOW (CU METERS) ***
=====================================================================
=============================================================
*** TOTALS ***
235
THE MAXIMUM INUNDATED AREA (DEPTH > 0.15 M) IS: 82625.000 SQ.
METERS
=====================================================================
=============================================================
236
C.3 Planos generados en entorno GIS, con los datos obtenidos del modelo Flo -
2D y el post procesador MAPPER PRO.
WATER EVENT
(TR = 50 AÑOS)
237
238
239
240
241
242
243
WATER EVENT
(TR = 100 AÑOS)
244
245
246
247
248
249
250
MUD EVENT (TR
= 50 AÑOS)
251
252
253
254
255
256
MUD EVENT (TR
= 100 AÑOS)
257
258
259
260
261
262
ANEXO D – PANEL FOTOGRAFICO
263
Figura D.2. Identificación de zonas de deslizamientos en la subcuenca del río Sicra
264
Figura D.3. Analisis de muestras de suelo en laboratoria de E.P Ingenieía Civil – Hvca
265
Figura D.4. Vuelo del dron sobre el area de estudio (Tramo del cauce del río Sicra)
266
ANEXO E – MATRIZ DE CONSISTENCIA
Problema general Objetivo general Antecedentes Hipótesis general Variable 01 Tipo: Aplicada.
¿Cuáles son las zonas de Definir zonas de peligro Internacionales La presente investigación por Subcuenca del río Nivel: Descriptivo.
peligro por flujos de lodo y por flujos de lodo y el nivel en el que se desarrolla Sicra.
escombros generados por escombros generados Aguirre et al. (2000), “Flujos no presenta hipótesis. Según
Método: Deductivo.
y deposición de barros y Dimensiones:
parámetros por parámetros Hernández, et al. (2014,
Diseño: No
escombros”. p.104), “No en todas las
geomorfológicos, geomorfológicos, Geomorfológicas
experimental –
hidrológicos, y hidrológicos, y investigaciones cuantitativas Hidrológicas
Oliveira (2012), “Estudio del longitudinal.
geodinámicos en la geodinámicos en la se plantean hipótesis. El hecho Geodinámicos
impacto de los flujos de
subcuenca del rio Sicra, subcuenca del río Sicra, de que planteemos o no M O
escombros en la cuenca del Variables 02
Lircay – Huancavelica? Lircay – Huancavelica. hipótesis depende de un factor
río San Antonio - Brasil”.
esencial. El alcance inicial del Flujos de lodo y Donde:
Problemas específicos Objetivos específicos escombros.
Sepúlveda et al. (2016), estudio”.
M: Muestra.
¿Qué parámetros Determinar los “Metodología para Dimensiones:
evaluación de riesgo por O: Observación
geomorfológicos parámetros Cantidad de
influyen en el geomorfológicos flujo de detritos detonados
sedimento Población y muestra
transporte de flujos de que influyen en el por lluvias caso Útica -
transportado
lodo y escombros, los transporte de flujos Colombia”.
267
cuales generan zonas de de lodo y Profundidad de
Bateman et al. (2017), Población:
peligro, en la subcuenca escombros, los sedimento
desarrolló una investigación Subcuenca del rio
del río Sicra, Lircay – cuales generan Magnitud del
titulada “Modelo Sicra.
Huancavelica? zonas de peligro, en caudal de flujo
bidimensional para
la subcuenca del río hiperconcentrado Muestra: No
simulación de flujos
¿Qué parámetros –
Sicra, Lircay
detríticos FLATModel. Zonas de peligro probabilística, de
hidrológicos influyen Huancavelica. (mapas de peligro progresiva 0+565 km
Aplicación a una cuenca del
en el transporte de
al 1+744 km del
Pirineo Catalán”.
flujos de lodo y Identificar los
cauce principal del
escombros, los cuales parámetros
Nacionales río Sicra (área urbana
generan zonas de hidrológicos que
del distrito de
peligro, en la subcuenca influyen en el Silva (2001), “Estudio de Lircay).
del río Sicra, Lircay – transporte de flujos control de protección de
Huancavelica? de lodo y flujos de lodo en las Técnica e
escombros, los microcuencas de Quirío y instrumentos:
¿Qué parámetros cuales generan Pedregal-Chosica”.
geodinámicos influyen zonas de peligro, en Técnicas:
en el transporte de la subcuenca del río Castillo (2006), “Aplicación Observación
flujos de lodo y Sicra, Lircay – de un modelo numérico de sistemática.
escombros, los cuales flujos de escombros y lodo en
Huancavelica. Acopio de datos
generan zonas de una quebrada en el Perú”.
peligro, en la subcuenca Analizar los Instrumentos:
Pinaud (2014), “Estudio
parámetros Fichas
geodinámicos que experimental de las
bibliográficas
268
del río Sicra, Lircay – influyen en el características dinámicas de Datos de campo
Huancavelica? transporte de flujos un dispositivo disipador
Técnicas de
de lodo y viscoso a escala reducida”.
procesamiento y
escombros, los
análisis de datos:
cuales generan Gamión (2014), “Modelo de
269