UNIVERSIDAD MAYOR DE SAN ANDRES

FACULTAD DE TECNOLOGIA
CARRERA DE TOPOGRAFIA Y GEODESIA
CENTRO DE INVESTIGACIONES Y APLICACIONES GEOMÁTICAS (CIAG)

"EVALUACION HIDRAULICA DEL RIO LA PAZ; TRAMO PUENTE

ENCUENTRO DE LOS RÍOS CHOQUEYAPU E IRPAVI HASTA LA
CONFLUENCIA CON EL RIO HUAYÑAJAHUIRA, PARA UNA AVENIDA
MAXIMA DE 50 AÑOS DE RETORNO"

TESIS PRESENTADA PARA LA OBTENCIÓN DEL GRADO DE MAESTRÍA EN:

“CIENCIAS GEOMÁTICAS APLICADA A LA GESTIÓN TERRITORIAL,
RECURSOS NATURALES Y MEDIO AMBIENTE”

POSTULANTE: Lic. PEDRO MEDINA QUISPE

TUTOR : M.Sc. Ing. JOSE LUIS DELGADO ALVAREZ

LA PAZ – BOLIVIA
2018
AGRADECIMIENTO

Agradezco en primer lugar a la Dirección de la Carrera de Topografía y Geodesia de la Facultad

de Tecnología de la Universidad Mayor de San Andrés por haberme dado la oportunidad de ser
parte del Programa de Maestría en Ciencias Geománticas.

Al Ingeniero José Luis Delgado, por su tutoría en el desarrollo del trabajo de investigación,
brindándome su valioso tiempo y apoyo en base a sus conocimientos sobre temas de
Modelamiento Hidrológico e Hidráulica de sistemas de drenaje de la cuenca objeto de estudio.

Un agradecimiento profundo, a mi esposa Elena Torrez y mis hijas Mariela, Helen y Madeline.
Gracias por su incondicional apoyo y soporte cuando yo más lo necesitaba. También agradezco a
mi papa Benedicto Medina Torrez por su apoyo incondicional. Me enseñaron que no hay obstáculo
más grande en la vida que nuestras ganas y deseos de superación puedan superar.

Finalmente agradezco a Dios y la Virgencita, por acompañarme día a día en cada paso que doy;
sin su fortaleza y guía nada de esto sería posible, teniendo presente que así como con este logro
podré siempre de su mano alcanzar otros que me proponga durante mi carrera profesional.
Contenido

CAPITULO I......................................................................................................................................... 1
1.1 INTRODUCCIÓN .......................................................................................................................... 1
1.2 PLANTEAMIENTO DEL PROBLEMA ........................................................................................... 3
1.3 JUSTIFICACIÓN ..................................................................................................................... 3
1.4. OBJETIVOS ................................................................................................................................. 4
1.4.1 OBJETIVO GENERAL ........................................................................................................ 4
1.4.2 OBJETIVOS ESPECÍFICOS ............................................................................................... 4
1.5 ALCANCE...................................................................................................................................... 5
CAPITULO II........................................................................................................................................ 7
MARCO TEORICO .............................................................................................................................. 7
2.1 CUENCAS ..................................................................................................................................... 7
2.2 HIDROGRAFÍA DE SUPERFICIE ................................................................................................. 8
2.2.1 HIDROGRAFÍA .......................................................................................................................... 8
2.2.2 HIDROLOGÍA ............................................................................................................................. 8
2.2.3 CICLO HIDROLÓGICO DEL AGUA .......................................................................................... 9
2.2.4 PRECIPITACIÓN ..................................................................................................................... 12
2.2.5 ESTACIONES PLUVIOMETRICAS O PLUVIOGRAFOS ........................................................ 13
2.2.6 BALANCE HÍDRICO ................................................................................................................ 14
2.2.7 EL ESCURRIMIENTO SUPERFICIAL ..................................................................................... 15
2.2.8 NUMERO DE CURVA CN ........................................................................................................ 17
2.2.9 CAUDAL ................................................................................................................................... 17
2.2.10 AVENIDA ................................................................................................................................ 18
2.2.11 PERIODO DE RETORNO ...................................................................................................... 18
2.2.12 ESTACIONES HIDROMETRICAS ......................................................................................... 18
2.3 MODELOS HIDROLÓGICOS ..................................................................................................... 19
2.3.1 MODELADO HIDROLÓGICO HEC-HMS ................................................................................ 19
2.3.2 MODELO SWAT ...................................................................................................................... 21
2.3.3 MODELO HIDRÁULICO........................................................................................................... 22
2.3.4 MODELO IBER ........................................................................................................................ 22
2.4 SISTEMAS DE INFORMACIÓN GEOGRÁFICA (SIG) ............................................................... 23
2.4.1 FUNCIONES DE UN SIG ......................................................................................................... 24
2.4.2 SISTEMAS DE INFORMACIÓN GEOGRÁFICA APLICADO A PLANOS DE INUNDACIÓN . 25
2.5 LEVANTAMIENTO TOPOGRÁFICO POR EL METODO INDIRECTO ...................................... 28
2.5.1 UNMANNED AERKAL VEHICLE UAV O DRON ..................................................................... 29
2.5.2 PUNTOS DE CONTROL O DE APOYO FOTOGRAMETRICO .............................................. 30
2.5.3 PIX4D MAPPER ....................................................................................................................... 31
2.6 CARTOGRAFIA DE RIESGOS ................................................................................................... 31
2.6.1 RIESGO NATURAL .................................................................................................................. 31
2.6.2 RIESGO DE INUNDACION...................................................................................................... 32
2.6.3 CARTOGRAFIA DE RIESGOS DE INUNDACION .................................................................. 32
CAPITULO III..................................................................................................................................... 34
MATERIALES Y METODOLOGIA .................................................................................................... 34
3.1 MATERIALES .............................................................................................................................. 34
3.2 METODOLOGÍA .......................................................................................................................... 35
3.2.1 DELIMITACIÓN DEL ÁREA DE LA CUENCA Y SUB-CUENCAS DEL RIO LA PAZ NORTE
CON EL SOFTWARE SWAT. ........................................................................................................... 36
3.2.2 CARACTERIZACIÓN BIOFÍSICA DE LA CUENCA Y SUB-CUENCAS. ................................ 36
3.2.3 CARACTERIZACIÓN CLIMÁTICA DE LA CUENCA DEL RIO LA PAZ NORTE. ................... 36
3.2.4 DETERMINACIÓN DE LA LLUVIA DEL PROYECTO CON TIEMPO DE RETORNO 50
AÑOS. ............................................................................................................................................... 37
3.2.5 DETERMINACIÓN DEL NUMERO DE CURVA. ..................................................................... 37
3.2.6 GENERACIÓN DEL MODELO HEC-HMS. .............................................................................. 37
3.2.7 LEVANTAMIENTO TOPOGRÁFICO EN DETALLE DEL CURSO DEL RIO .......................... 37
3.2.8 GENERACIÓN DE MDT PARA EL SECTOR DEL RIO OBJETO DE ESTUDIO. ................... 38
3.2.9 DETERMINACIÓN DE LOS PARÁMETROS MORFOMÉTRICOS DEL RIO. ........................ 38
3.2.10 GENERACIÓN DEL MODELO DE DINÁMICA FLUVIAL CON IBER.................................... 38
3.2.11 DETERMINACIÓN DE LA CAPACIDAD HIDRÁULICA DEL CURSO DEL RIO ................... 38
3.2.12 CARTOGRAFÍA DE RIESGOS DE DESBORDES, COTA Y VELOCIDAD DE AGUA. ........ 38
3.2.13 PROPUESTA DE MEDIDAS DE MITIGACIÓN. .................................................................... 39
CAPITULO IV .................................................................................................................................... 40
RESULTADOS Y DISCUCIONES..................................................................................................... 40
4.1 CUENCA DE ESTUDIO .............................................................................................................. 40
4.1.1 MODELO DIGITAL DEL TERRENO ........................................................................................ 40
4.1.2 CARTOGRAFÍA ....................................................................................................................... 41
4.1.3 DELIMITACION DEL AREA DE LA CUENCA DEL RIO LA PAZ (NORTE) ............................ 44
4.1.4 DETERMINACIÓN DE PARÁMETROS MORFOMÉTRICOS ................................................. 44
4.2 CARACTERIZACION BIOFISICA DE LA CUENCA DEL RIO LA PAZ (NORTE) ...................... 45
4.2.1 CARACTERIZACION HIDROLOGICA ..................................................................................... 47
4.2.2 CARACTERIZACION GEOLOGICA ........................................................................................ 53
4.2.3 GEOMORFOLOGÍA ................................................................................................................. 55
4.2.4 PENDIENTE ............................................................................................................................. 58
4.2.5 CARACTERIZACION DE LA COBERTURA VEGETAL .......................................................... 61
4.2.6 CARACTERIZACION DE LA COBERTURA DE USO DE LA TIERRA ................................... 63
4.3 CARACTERIZACION CLIMATICA Y CURVAS IDF PARA 50 AÑOS DE RETORNO ............... 65
4.3.1 TEMPERATURA ...................................................................................................................... 66
4.3.2 EVAPORACIÓN ....................................................................................................................... 67
4.3.3 PRESIÓN BAROMÉTRICA (ATMOSFÉRICA) ........................................................................ 67
4.3.4 HUMEDAD RELATIVA ............................................................................................................. 67
4.3.5 VIENTOS .................................................................................................................................. 68
4.3.6 PRECIPITACIONES ................................................................................................................. 68
4.3.7 CURVAS DE INTENSIDAD DURACION Y FRECUENCIA ..................................................... 70
4.4 DETERMINACION DE LA LLUVIA DEL PROYECTO ................................................................ 71
4.5 DETERMINACIÓN DE NÚMEROS DE CURVA (NC) ................................................................ 72
4.6 MODELO DE TRANSFORMACION DE LA LLUVIA EN ESCORRENTIA HEC-HMS ............... 75
4.6.1 ANÁLISIS DE HOMOGENEIDAD DE PRECIPITACIONES .................................................... 77
4.6.2 CÁLCULO DE CAUDALES CON MODELO HEC-HMS .......................................................... 79
4.6.3 PROCESAMIENTO DE LOS DATOS DE LA CUENCA HIDROGRAFICA ............................. 82
4.7 LEVANTAMIENTO TOPOGRÁFICO EN DETALLE ................................................................... 86
4.7.1 TOMA DE IMÁGENES FOTOGRAMÉTRICAS CON UAV o DRON. ...................................... 86
4.7.2 DETERMINACIÓN DE PUNTOS DE CONTROL .................................................................... 88
4.8 GENARACION DEL MODELO DIGITAL DE TERRENO y ORTOFOTOS ................................. 89
4.8.1 GENERACIÓN DE MDT/MDS ................................................................................................. 90
4.8.2 GENERACIÓN DE ORTOFOTOS ........................................................................................... 91
4.9 MODELO DE DINAMICO FLUVIAL IBER ................................................................................... 92
4.9.1 CRITERIOS HIDRÁULICOS .................................................................................................... 93
4.9.2 MODELAMIENTO DINAMICO ................................................................................................. 95
4.9.3 ANALISIS DE RESULTADOS ................................................................................................ 100
4.10 CARTOGRAFIA DE RIESGOS DE DESBORDE.................................................................... 107
4.10.1 FACTORES DE RIESGO ..................................................................................................... 108
3
4.10.2 PRESENTACION GRAFICA DE LOS RESULTADOS CAUDAL PICO 57m /s .................. 108
4.11 PROPUESTA DE MEDIDAS DE MITIGACION ...................................................................... 110
4.11.1 OBJETIVO DE LA PROPUESTA DE PREVENCION Y MITIGACION DE RIESGO ........... 111
4.11.2 ACTIVIDADES DE PREVENCION Y MITIGACION DE RIESGOS ..................................... 111
4.11.3 PRESUPUESTO .................................................................................................................. 112
4.11.4 CRONOGRAMA ................................................................................................................... 112
CAPITULO V ................................................................................................................................... 113
CONCLUSIONES Y RECOMENDACIONES .................................................................................. 113
CAPITULO VI .................................................................................................................................. 116
BIBLIOGRAFIA................................................................................................................................ 116
 RESUMEN

La evaluación hidráulica del rio La Paz (Norte); tramo puente encuentro de los ríos Choqueyapu e
Irpavi hasta la confluencia con el rio Huayñajahuira, para una avenida máxima de 50 años de
retorno, se realiza para tener conocimiento pleno del riesgo de desborde del rio ante una probable
falta de capacidad hidráulica del cauce del rio.

Para el cumplimiento del objetivo de la evaluación, la investigación parte de la delimitación del área
de la cuenca y sub-cuencas del rio La Paz (Norte) a partir del Modelo Digital de Terreno (SRTM
30m); sub-cuencas correspondiente a los ríos: Choqueyapu, Orkojahuira, Irpavi, Achumani y
Huayñajahuira.

Las sub-cuencas son caracterizadas biofísicamente a partir de los siguientes mapas: Cobertura y
Uso Actual del Suelo, Suelo, Hidrológico, Morfológicos y Pendientes, Vegetación y la información
Climática para el cálculo del Número de Curva ponderado de las unidades hidrológicas
relaconadas con los parámetros del Soil Conservancy Service. Asimismo, a partir de los datos de la
Máxima Precipitación en 24 horas de la estación San Calixto, se deriva la curva (IDF) para 50 años
de retorno y obtener la lluvia del proyecto.

Con la información de: La caracterización biofísica y climática reflejada en el Numero de Curva de

las sub-cuencas. Así como la curva IDF de la estación San Calixto de la que deriva la lluvia del
proyecto (TR 50 años) y Los parámetros morfométricos de las sub-cuencas, se realiza el
modelamiento hidrográfico con HEC-HMS para la transformación de la precipitación a
escurrimiento superficial a efecto de obtener el caudal pico.

Para establecer la morfología del cauce del rio, se ha realizado el levantamiento topográfico por el
método indirecto. Esto comprende: Toma de imágenes fotogramétricas con sistemas UAV (DRON),
Determinación de Puntos de Control, Generación de la Ortofoto y Modelo Digital del Terreno. Este
último refleja la topografía del área de estudio.

Con la información de: Caudal pico y el MDT del tramo reflejada en la topografía, se realiza el
modelamiento hidráulico en el software IBER. El análisis y evaluación de los resultados del
modelamiento de los factores de riesgo identificados como: Cota de agua y Velocidad de agua, se
realiza respecto a la topografía, esta información temática es la base para la conformación de la
cartografía de riesgos. Finalmente se propone las acciones de prevención y mitigación.
LISTA DE FIGURAS

Figura 1.1 Ubicación del área de estudio

Figura 2.1 Tipos de cuencas hidrológicas.
Figura 2.2 Ciclo hidrológico o del agua.
Figura 2.3 Pluviómetro
Figura 2.4 Pluviógrafo
Figura 2.5 Modelo de datos raster y vectorial
Figura 2.6 EBEE Sensefly
Figura 3.1 Metodología de Investigación (elaboración propia)
Figura 4.1 Modelo Digital del Terreno de la cuenca del rio La Paz (Norte).
Figura 4.2 Mapa de delimitación de la cuenca del rio La Paz (Norte).
Figura 4.3 Mapa hidrográfico de la cuenca del rio La Paz (Norte).
Figura 4.4 Mapa de ubicación de la cuenca y sub-cuencas del rio La Paz (Norte)
Figura 4.5 Mapa Geológico.
Figura 4.6 Mapa de Geomorfológico (información ZONISIG)
Figura 4.7 Mapa de pendientes (Elaboración propia a partir del MDT SRTM 30m)
Figura 4.8 Mapa de vegetación de la cuenca del rio La Paz (Norte)
Figura 4.9 Mapa Cobertura y Uso actual del Suelo de la cuenca del rio La Paz (Norte).
Figura 4.10 Ingreso de masas de aire en eventos de baja intensidad (Ramirez, 1995).
Figura 4.11 Régimen de precipitaciones mensual; cuenca del río La Paz (Norte)
Figura 4.12 Mapa de ubicación de estaciones meteorológicas de la cuenca
Figura 4.13 Presentación gráfica del modelo
Figura 4.14 Datos de la serie de tiempos
Figura 4.15 Sumario global de los resultados del procesamiento de los datos.
Figura 4.16 Grafico del hidrograma y Tabla Sumario resultado
Figura 4.17 Grafico y tabla sumario de la descarga y flujo hidrográfico
Figura 4.18 Grafico y tabla sumaria de la máxima descarga de la unión de los ríos Choqueyapu y
Orkojahuira
Figura 4.19 Grafico y tabla sumaria de la máxima descarga de la unión de los ríos Choqueyapu e
Irpavi
Figura 4.20 Modelo Digital Terreno formato raster
Figura 4.21 Ortofoto del área de proyecto
Figura 4.22 Polígonos para espaciamiento diferenciado de la malla
Figura 4.23 Malla con espaciamiento diferenciado del cauce y las franjas de seguridad del rio
Figura 4.24 Malla del área de estudio fondo ortofoto
Figura 4.25 Cota de Agua (m)
Figura 4.26 Velocidad de agua (m/s)
Figura 4.27 Remoción de material de arrastre
3 3 3
Figura 4.28 Imágenes de la cota de agua para los caudales 57m /s, 75m / y 100m /s
3 3 3
Figura 4.29 Imágenes de velocidad de agua para los caudales 57m /s, 75m /s y 100m /s
3
Figura 4.30 Perfil de Cota de Agua para el caudal de 15m /s
3
Figura 4.31 Perfil de Cota de Agua para el caudal de 40m /s
3
Figura 4.32 Perfil de Cota de Agua para el caudal de 60m /s
3
Figura 4.33 Perfil de Cota de Agua para el caudal de 80m /s
Figura 4.34 Velocidades de agua para los caudales de 15m3/s, 40m3/s, 60m3/s y 80m3/s.
3
Figura 4.35 Mapa de Velocidad de Agua para un caudal pico de 57 m /s
3
Figura 4.36 Mapa de cota de Agua para un caudal pico de 57 m /s

LISTA DE IMÁGENES

Imagen 4.1; Río Choqueyapu, Río Irpavi

Imagen 4.2; Nacientes del Rio Choqueyapu, Río Kaluyo
Imagen 4.3; Río Orkojahuira, Confluencia con Choqueyapu.
Imagen 4.4; Río Palcoma, Río Irpavi
Imagen 4.5; Río Huayllani
Imagen 4.6; Río Achumani
Imagen 4.7; Rio Huayñajahuira
Imagen 4.8 PCT1
Imagen 4.9 PCT2
Imagen 4.10 PCT3
Imagen 4.11 PCT4
LISTA DE TABLAS

Tabla 4.1 Parámetros Morfométricos de la cuenca y sub-cuencas del rio La Paz.

Tabla 4.2 Temperatura Estación San Calixto La Paz (con datos de SENAMHI)
Tabla 4.3 Evaporación media mensual: Estación San Calixto La Paz (datos SENAMHI).
Tabla 4.4 Humedad relativa media mensual: Estación San Calixto La Paz (datos SENAMHI).
Tabla 4.5 Intensidad de la lluvia con 50 años de retorno; Estación San Calixto.
Tabla 4.6 Resumen del Número de Curva Ponderado de la sub-cuenca
Tabla 4.7 Lista de coordenadas de las estaciones meteorológicas
Tabla 4.8 Primera corrida del Vector Regional
Tabla 4.9 Plan de vuelo para la realización del vuelo y la toma de imágenes a color
Tabla 4.10 Coordenadas de Puntos de Control fotoidentificables
 CAPITULO I

1.1 INTRODUCCIÓN

El informe de la cumbre del cambio climático que se realizó en Cochabamba (mayo 2010)
señala que: debido al cambio climático, se incrementaron las lluvias intensas causantes
de inundaciones que podrían llevar a daños significativos en las infraestructuras civiles y
sembradíos, de acuerdo a escenarios regionales, se espera una reducción de los
periodos de lluvia, con una probabilidad creciente de fuertes lluvias en periodos de tiempo
cortos. En este sentido, Bolivia aún carece de métodos para estudiar los patrones de las
lluvias y las respuestas hidráulicas de los ríos, lo que incrementa la vulnerabilidad de la
población, en especial de los asentamientos humanos no planificados. Asimismo, es casi
generalizado que las cuencas hidrográficas en nuestro país no cuentan con sistemas de
manejo integrado de recursos hídricos en marcha, ni contemplan actividades tendientes a
la protección de las cuencas, ni los cauces. Este descuido se traduce en la degradación
de suelos, la declinación de la recarga de acuíferos y la reducción de la disponibilidad de
agua, que en su conjunto reduce la productividad y los ingresos de la población.

La ciudad de La Paz, está rodeada por la cadena montañosa de los Andes, donde
sobresalen imponentes y altivos el nevado lllimaní, Mururata y Milluni, además de ser una
ciudad cosmopolita, multiétnica y pluricultural. El 19 de febrero del 2002, la ciudad de La
Paz se convirtió en el escenario de una pesadilla, sus calles se transformaron en ríos
caudalosos, varios inmuebles se desplomaron, hubo vendedores callejeros que perdieron
todo, vehículos que avanzaban descontrolados sobre el granizo e impulsados por la
fuerza de los elementos que forman parte del caudal. El servicio de energía eléctrica se
interrumpió, 59 personas murieron arrastradas por el caudal con material de arrastre o
atrapadas entre los escombros de los edificios que sufrieron daños estructurales. (El
Diario, febrero del 2002).

El presente estudio, permitirá la realización de la caracterización biofísica de la cuenca y

sub-cuencas hidrográficas que cursan la ciudad de La Paz. Así como el modelamiento del
caudal y evaluación de la capacidad hidráulica del tramo de rio La Paz (Norte), aplicando
para este propósito sistemas digitales y software aplicativo en la plataforma SIG y
fundamentará la propuesta de ”Evaluación de la capacidad hidráulica del rio La Paz
(Norte) para una avenida máxima de 50 años de posibilidad de ocurrencia; tramo puente

Página 1
encuentro de los ríos Choqueyapu e Irpavi hasta la confluencia con el rio Huayñajahuira”,
estimando el caudal hidrográfico de superficie máximo y evaluando la capacidad
hidráulica del curso del rio en el indicado tramo. Ver figura -1.1; Ubicación del área de
estudio.

Figura 1.1. Ubicación del área de estudio (Elaboración propia).

El estudio se realizará utilizando la tecnología de sensor remoto Unmanned Aerial Vehicle

UAV (DRON) para el levantamiento topográfico en detalle (método indirecto) y sistema de
información geográfica (SIG); tecnología disponible para; Generación de información de

Página 2
levantamiento topográfico, Manejo de información geográfica, Modelamiento del flujo
hidrográfico y Evaluación de la capacidad hidráulica del mencionado tramo.

La información geográfica en formato digital, así como los resultados del procesamiento y
análisis de la información geográfica de la evaluación es de vital importancia ya que
permitirá tomar acciones preventivas a riesgos naturales de desborde e inundaciones en
las áreas circundantes al indicado sector del rio La Paz.

1.2 PLANTEAMIENTO DEL PROBLEMA

Las intensas precipitaciones anuales en la cabecera de la cuenca del rio La Paz (Norte),
generan una escorrentía que arrastra gran cantidad de áridos hacia la parte baja de la
cuenca. Esto se puede evidenciar en la zona del puente encuentro entre los ríos
Choqueyapu e Irpavi; espacio donde se depositan y se dispersan gran cantidad de
material de arrastre formando un abanico aluvial, esto genera un descontrol del flujo
hidrográfico y una eventual insuficiencia de la capacidad hidráulica del cauce del rio
compuesto por vertedero natural que cursa aguas abajo. Anualmente el GAMLP tiene que
realizar la remoción de sedimentos compuesto por bloques de piedra y áridos para
mantener la capacidad hidráulica del cauce del rio en el sector del puente encuentro de
los ríos Choqueyapu e Irpavi hasta la confluencia con el rio Huayñajahuira ubicado en la
Zona de la Florida de la ciudad de La Paz.

1.3 JUSTIFICACIÓN

La cuenca hidrográfica del río La Paz (Norte), está conformada por las sub-cuencas
hidrográficas correspondientes a los ríos Choqueyapu, Orkojahuira, Irpavi, Achumani y
Huayñajahuira. La cuenca hidrográfica del rio La Paz (Norte) tiene un área de 498.5 Km2
y de acuerdo al modelo digital de terreno SRTM 30m. (Shuttle Radar Topography
Mission), presenta un desnivel aproximado de 2170 metros entre el punto más alto en la
naciente de los ríos y el punto más bajo ubicado en el puente cruce a Mallasa zona
Aranjuez.

El presente estudio se justifica técnicamente, ya que éste generará información

especializada, así como el modelo de flujo hidrográfico de superficie basado en la
caracterización hidrológica y biofísica de la cuenca hidrográfica del rio La Paz (Norte),

Página 3
identificando las causas para la acumulación y dispersión de material de arrastre para la
evaluación de la capacidad hidráulica del tramo sujeto de estudio. Además el estudio
servirá también como herramienta para la toma de decisiones y acciones necesarias,
evitar asentamientos humanos en las áreas circundantes al indicado tramo del rio con
riesgos de desbordes e inundaciones. Para cumplir dicho cometido, se realizará el
procesamiento de datos hidrográficos de superficie, de lo que se obtendrá el conocimiento
necesario para la conformación del sistema de información geográfica en 3D y el modelo
de simulación para la evaluación del flujo de caudales de superficie y capacidad hidráulica
del vertedero natural para una avenida con tiempo de retorno de 50 años.

Los desbordes e inundaciones son un fenómeno natural que en su mayoría son causados
por las intensas precipitaciones en la cuenca hidrográfica del rio La Paz (Norte), éstas
afectan de manera significativa a la estabilidad de las obras civiles como los muros de
contención precarias existentes en el sector. Los riesgos de desborde en periodo de lluvia
están latentes y causan profunda preocupación en la población circundante. Por lo que es
necesario conformar el sistema de información geográfica en 3D, el modelo de simulación
y control del flujo hidrográfico de superficie y asegurar la capacidad hidráulica del
vertedero natural entre el puente encuentro de los ríos Choqueyapu e Irpavi hasta la
confluencia con el rio Huayñajahuira (zona la Florida), de manera que ayude a las
autoridades en la toma de decisiones oportunas para minimizar los efectos y daños
colaterales causados por los desbordes e inundación del rio La Paz (Norte).

1.4. OBJETIVOS

1.4.1 OBJETIVO GENERAL

Realizar la evaluación hidráulica del curso del rio La Paz (Norte); sector puente encuentro
de los ríos Choqueyapu e Irpavi hasta la confluencia con el rio Huayñajahuira, para una
avenida máxima de 50 años de retorno y propuesta de medidas de prevención de
desbordes e inundaciones.

1.4.2 OBJETIVOS ESPECÍFICOS

- Realizar la caracterización biofísica e hidrológica de la cuenca y sub-cuencas

hidrográficas del rio La Paz (Norte) que afectan a la ciudad de La Paz.

Página 4
 - Realizar el modelamiento hidrológico para estimar los caudales de los ríos que
afectan a la ciudad de La Paz; ocasionadas por precipitaciones máximas para un
periodo de retorno de 50 años.
- Realizar el levantamiento topográfico del tramo puente encuentro de los ríos
Choqueyapu e Irpavi hasta la confluencia con el rio Huayñajahuira (Zona La
Florida).
- Realizar el modelamiento de dinámica fluvial para evaluar la capacidad hidráulica
del vertedero natural del rio La Paz (Norte); entre el puente encuentro de los ríos
Choqueyapu e Irpavi hasta la confluencia con el rio Huayñajahuira (Zona la
Florida).
- Cartografiar los riesgos de desborde e inundaciones y elaborar la propuesta de
medidas preventivas.

1.5 ALCANCE

La presente investigación comprende la realización de las siguientes actividades en el

área de investigación de la cuenca y sub-cuencas del Rio La Paz (Norte) entre el
puente encuentro entre los ríos Choqueyapu e Irpavi hasta la confluencia con el rio
Huayñajahuira a efectos de proponer las medidas preventivas:
- Delimitación de la cuenca y sub-cuencas del rio La Paz (Norte) utilizando el Modelo
Digital de Terreno.
- Caracterización biofísica e hidrología de las sub-cuencas del rio La Paz (Norte) con
la información temática aplicable.
- Caracterización climática y cálculo de curvas IDF para determinar la lluvia del
proyecto para un tiempo de retorno de 50 años.
- Determinación de los parámetros morfométricos de la cuenca y sub-cuencas del rio
La Paz (Norte).
- Determinación del Numero de Curva (NC) de las sub-cuencas.
- Modelamiento hidrográfico para la transformación de la lluvia en escorrentía con el
software HEC-HMS.
- Levantamiento Topográfico a detalle del curso del rio La Paz (Norte) del sector
puente encuentro entre los ríos Choqueyapu e Irpavi hasta la confluencia con el rio
Huayñajahuira.

Página 5
- Determinación de los parámetros morfométricos del rio La Paz (Norte) del indicado
tramo.
- Modelamiento de dinámica fluvial para la evaluación de la capacidad hidráulica del
curso del rio La Paz (Norte) del indicado tramo utilizando el modelo dinámico IBER.
- Generación de la cartografía de riesgos de desbordes que sustentan la toma de
medidas de mitigación.

Página 6
 CAPITULO II
MARCO TEORICO

En el presente capítulo se enuncia las definiciones y conceptos que sustentan a la

fundamentación teórica aplicada en el presente estudio de evaluación de la capacidad
hidráulica del vertedero natural del rio La Paz (Norte); tramo puente encuentro entre los
ríos Choqueyapu e Irpavi hasta la confluencia con el rio huayñajahuira, este apartado
brinda al lector las definiciones, conceptos fundamentales y conocimiento previos
descritas en la literatura que ayudan al desarrollo de la investigación. En la construcción
del marco teórico de la presente investigación se ha considerado los siguientes aspectos:
- Cuencas hidrográficas y caracterización biofísica.
- Hidrografía de superficie y de dinámica fluvial (HEC-HMS e IBER).
- Modelos Hidrológicos.
- Levantamiento Topográfico por el método indirecto utilizando toma de imágenes con
DRON y georeferenciación con sistema GPS.
- Sistemas de Información Geográfica.
- Cartografía de riesgos de inundación.

2.1 CUENCAS

"La cuenca o unidad hidrogeológica es aquella superficie en la cual el agua precipitada se

transfiere a las partes topográficas bajas por medio del sistema de drenaje,
concentrándose generalmente en un colector que descarga a otras cuencas aledañas, o
finalmente al océano"; Agustin F. Breña y Marco A. Jacobo (2006). La cuenca hidrológica,
junto con los acuíferos, son las unidades fundamentales de la hidrología.

Desde el punto de vista de su salida o desagüe existen dos tipos de cuencas: endorreicas
(cerradas) y exorreicas (abiertas).
a) En el primer tipo, el punto de salida o desagüe se ubica dentro de los límites de la
cuenca y generalmente es un lago. Es decir, el agua concentrado en un lago se
sale a través del proceso de infiltración por debajo de la superficie terrestre.
b) En el segundo tipo, el punto de salida se localiza en parte baja de los límites de la
cuenca y a su vez la descarga se vierte en una corriente o en el mar; Ver figura 2.1
Tipos de cuencas hidrológicas.

Página 7
Figura 2.1. Tipos de Cuencas hidrológicas (fuente https://www.google.com/search).

2.2 HIDROGRAFÍA DE SUPERFICIE

2.2.1 HIDROGRAFÍA

Llamamos Hidrografía al estudio de todas las masas de agua de la Tierra, más

propiamente a la medida, recopilación y representación gráfica de los datos relativos a;
fondo del océano, las costas, las mareas y los ríos; de manera que se puedan representar
gráficamente sobre un mapa o sobre una carta hidrográfica, emparentada con ella está la
hidrología que se dedica al estudio de la distribución y las propiedades del agua de la
atmósfera y la superficie terrestre.

2.2.2 HIDROLOGÍA

"La hidrología estudia los procesos hidrológicos influenciados por el clima, superficie y
sub-superficie del suelo. La comprensión de las interacciones entre la hidrología y el
ecosistema, permite transferir predicciones hidrológicas de una zona a otra, aunque la
escala espacial juega un rol muy importante en la comparación de los procesos
hidrológicos"; Juan C Mecca (2008).

"La Hidrología es la ciencia que se enfoca al ciclo hidrológico global y a los procesos
involucrados en la parte continental de dicho ciclo"; Dingman (1994). Agustin F. Breña y
Marco A. Jacobo (2006) amplia indicando que "la hidrología es la geociencia que describe
y predice:

Página 8
- Las variaciones espaciales y temporales del agua en las etapas terrestre, oceánica y
atmosférica del sistema hídrico global, y

- El movimiento del agua sobre y debajo de la superficie terrestre, incluyendo los

procesos químicos, físicos y biológicos que tienen lugar a lo largo de su trayectoria".

"La hidrología clásica ha sido separada en el campo superficial y el subterráneo, donde el

concepto de infiltración se reduce al nivel de la superficie topográfica, entendiendo como
variables externas al sistema en estudio los intercambios de flujo en la zona no saturada.
Pero, la infiltración depende fuertemente del estado de redistribución de humedades en el
perfil de suelo y del almacenamiento de humedad en la zona no saturada". Juan C Mecca
(2008).

Entonces, la hidrología incluye las precipitaciones, la escorrentía, la humedad del suelo, la

evapotranspiración, el agua subterránea y el equilibrio de las masas glaciares.

Un rio es una corriente de agua que fluye por un cauce desde las tierras altas a las tierras
bajas y vierte las aguas en un rio colector, en el mar o en lago en una región endorreica.

2.2.3 CICLO HIDROLÓGICO DEL AGUA

” El ciclo hidrológico está compuesto por diferentes variables, las cuales se relacionan
entre sí por medio de los procesos hidrológicos. En general, se entiende por proceso a
una serie de acciones que producen un cambio o desarrollo en un sistema y para el caso
particular de la Hidrología, los procesos están asociados con aquellos fenómenos que
intervienen tanto en el movimiento del agua como en los cambios que sufre ésta en sus
características físicas, químicas y biológicas al desplazarse por diversos medios”; Agustin
F. Breña y Marco A. Jacobo (2006).

” El ciclo hidrológico se puede representar como un sistema, es decir, como una

estructura o volumen en el espacio, delimitada por una frontera, cuyos componentes
internos interactúan entre sí o con otros sistemas adyacentes” (Chow et al, 1988).

El ciclo hidrológico puede descomponerse en los siguientes términos:

El fenómeno de precipitación (P) es la variable que actúa sobre el sistema y ante la cual

Página 9
reacciona el resto. Una parte de ella se evapora antes de alcanzar la superficie terrestre.
El primer nivel de almacenamiento lo constituye la cubierta vegetal y se denomina
intercepción por el follaje (F). Parte del almacenamiento gotea desde las hojas llegando
luego al suelo y el resto se evapora.

”El agua que desde la precipitación directa y desde la intercepción por el follaje llega al
suelo, se llama precipitación efectiva (Pe)”; Juan C. Mecca (2008). La Pe que llega al
suelo constituye el segundo nivel de almacenamiento llamado superficial (S).

Las aguas subterráneas constituyen un tercer nivel de almacenamiento, que se toma en

cuenta en el modelamiento hidrológico de superficie.

”Los sistemas del ciclo hidrológico se pueden subdividir en tres: en un subsistema

atmosférico, en un subsistema de agua superficial y en un subsistema de agua
subterránea”; Agustin F. Breña y Marco A. Jacobo (2006). En cada subsistema se presenta
una capacidad de retención de volúmenes de agua, en cualquiera de sus fases, durante
un determinado intervalo de tiempo.

La capacidad de retención en un medio también recibe el nombre de almacenamiento y el

intervalo de tiempo que permanece un volumen recibe el nombre de tiempo de residencia.

- Subsistema atmosférico

”El subsistema atmosférico se abastece de la evaporación, fenómeno que se origina por

la incidencia de la energía proveniente del Sol y de la atmósfera alta”; Agustin F. Breña y
Marco A. Jacobo (2006), es decir se abastece de los volúmenes de vapor de agua que
llegan a la atmósfera desde la superficie del océano y/o desde la superficie de la tierra; en
este último caso, si existe una cubierta vegetal, se presenta el efecto combinado que
recibe el nombre de evapotranspiración. El vapor de agua podrá ser desplazado por los
procesos de circulación atmosférica a otras regiones geográficas donde, si se presentan
las condiciones adecuadas para la precipitación, es decir, se transformará en lluvia, nieve,
hielo (granizada) o rocío.

Página 10
- Subsistema del agua superficial

”Para el subsistema agua superficial la precipitación es una entrada, una vez eliminada el
agua interceptada, el volumen restante se transformará en escurrimiento sobre la
superficie del terreno”; Agustin F. Breña y Marco A. Jacobo (2006). La precipitación llega
al sistema de drenaje del área de captación para formar el escurrimiento superficial; parte
de estos escurrimientos o caudal de agua saldrán del subsistema como infiltración al
subsistema subterráneo o bien como escurrimiento por medio de los sistemas de drenaje
regionales que finalmente llegan al océano.

- Subsistema del agua subterránea

”La entrada a este subsistema de agua subterránea proviene básicamente de la

infiltración. Parte del volumen infiltrado al subsuelo percola a mayor profundidad, llegando
al almacenamiento del agua subterránea; la otra parte se almacena en la región
intermedia entre la superficie del terreno y el almacenamiento de agua subterránea de
mayor profundidad”; Agustin F. Breña y Marco A. Jacobo (2006).

El ciclo hidrológico también se puede conceptualizar como la sucesión de etapas que

atraviesa el agua al pasar de la tierra a la atmósfera y volver a la tierra: Evaporación;
desde los cuerpos de agua, mar o aguas continentales, Evapotranspiración; desde las
áreas cubiertas por vegetación, Condensación de nubes, Precipitación, Acumulación en el
suelo o masas de agua en mares o aguas continentales y luego la Reevaporación; Ver
figura 2.2; Ciclo hidrológico del agua.

El ciclo hidrológico involucra un proceso de transporte recirculatorio del agua

permanente, éste movimiento permanente del ciclo del agua se debe fundamentalmente a
tres causas; primero; el sol que proporciona la energía solar que eleva la temperatura
ambiente de la tierra y particularmente de los lagos y mares para elevar el agua en forma
de vapor de agua mediante los procesos de evaporación y evapotranspiración, segundo;
el viento, que transporta gran cantidad de vapor de agua desde las áreas húmedas a la
atmosfera a la altura de los espacios geográficos de las cuencas que conforman la tierra y
tercero la gravedad de la tierra, que hace que el agua condensada por efectos de la
acción del aire frio descienda al suelo mediante el proceso de precipitación y fluya por

Página 11
los cauces correspondientes dando lugar al escurrimiento.

Figura 2.2 Ciclo hidrológico o del agua. (Fuente https://www.significados.com/ciclo-del-agua/).

2.2.4 PRECIPITACIÓN

La precipitación constituye la principal entrada de agua dentro del Ciclo Hidrológico, este
fenómeno varía tanto espacial como temporalmente en una cuenca y sub-cuenca. ” Es la
cantidad de agua que cae en la superficie terrestre y proviene de la humedad atmosférica,
ya sea en estado líquido (llovizna y lluvias) o en estado sólido (escarcha, nieve, granizo)”
(Aparicio 1990). La precipitación es uno de los procesos meteorológicos más importantes
para la Hidrología, y junto a la evaporación constituyen la forma mediante la cual la
atmósfera interactúa con el agua superficial en el ciclo hidrológico del agua. La medición
de la lluvia se realiza en las estaciones meteorológicas a través de instrumentos llamados
pluviómetros o pluviografos y es uno de los datos necesarios para el cálculo del balance
hídrico; datos de lluvias que con mayor frecuencia se encuentran disponibles, si bien
puede variar la periodicidad y confiabilidad de éstos dependiendo del método de medición
y del instrumento de medición, así como la permanencia de las estaciones meteorológicas
a través del tiempo.

Página 12
2.2.5 ESTACIONES PLUVIOMETRICAS O PLUVIOGRAFOS

Sitios donde se encuentran ubicados e instalados los instrumentos de medición de la

intensidad de precipitaciones caídas, ya sea en forma de lluvia, nieve o granizo.
Expresada a través de la cantidad de litros o milímetros caídos por metro cuadrado.
Cualquier recipiente de boca ancha, cuya superficie sea conocida puede servir como
pluviómetro. Estos instrumentos con registradores automáticos instantáneos se
denominan Pluviógrafos.

Pluviómetro

Es un instrumento de precisión para poder leer el volumen de agua recogida en una

bureta de sección menor a la de la boca del pluviómetro. La lectura de agua recogida se
efectúa una vez al día. Ver figura 2.3; Pluviómetro.

Figura 2.3; Pluviómetro (fuente https://www.tiempo.com/ram/numero20/imagenes/instrumentos)

Pluviógrafo

También es un instrumento que tiene la capacidad de registrar la evolución de la

Página 13
precipitación en periodos de tiempo cada 5 minutos o cualquier otro intervalo de tiempo.
En algunos modelos el pluviógrafo está dotado de un flotador que hace subir a una
plumilla que registra gráficamente el llenado del recipiente a lo largo del tiempo. Ver figura
2.4 Pluviógrafo.

Figura 2.4; Pluviógrafo (fuente https://www.tiempo.com/ram/numero20/imagenes/instrumentos).

2.2.6 BALANCE HÍDRICO

”La cantidad de agua que circulan por el ciclo hidrológico se relacionan entre sí, conforme
al principio de conservación de la masa, debiendo existir un equilibrio entre entradas y
salidas. Las relaciones cuantitativas que representan este equilibrio constituyen las
expresiones del Balance Hidrológico. Para poder establecer el Balance Hidrológico es
imprescindible definir la unidad hidrogeológica o cuenca hidrográfica donde se aplica y el
intervalo de tiempo considerado”; Juan C. Mecca (2008). La disponibilidad de agua
producto de la precipitación anual en una cuenca está afectada fundamentalmente por la
evapotranspiración.

“La evapotranspiración es un proceso combinado que depende de la influencia de las

condiciones atmosféricas sobre el sistema agua-suelo-planta. Uno de los métodos más
sencillos para estimar la evapotranspiración, entre muchos existentes, es el de
Thomthwaite, el cual requiere, básicamente, de información climatológica y geográfica”
(Remenieras, 1974).

Página 14
2.2.7 EL ESCURRIMIENTO SUPERFICIAL

” El escurrimiento es la porción de la precipitación pluvial que ocurre en una zona o

cuenca hidrológica y que circula sobre o debajo de la superficie terrestre y que llega a una
corriente para ser drenada hasta la salida de una cuenca o bien alimentar un lago, si se
trata de cuencas abiertas o cerradas, respectivamente”; Agustin F. Breña y Marco A.
Jacobo (2006).

El cálculo de los escurrimientos superficiales se considera para dos objetivos: El

escurrimiento medio, para estimar el volumen de agua por almacenar o retener y los
escurrimientos máximos instantáneos para el diseño de obras civiles de conservación.

Escurrimiento Medio

En cuencas pequeñas o áreas de drenaje reducidas se aplica el Escurrimiento Medio,

para el cálculo del escurrimiento es necesario conocer el valor de la precipitación media a
partir del mapa de isoyetas, el área de drenaje y su coeficiente de escurrimiento en
función a la pendiente, uso y textura del suelo. La fórmula a utilizar sería la siguiente:

Vm = C Pm A

Donde: Vm = Volumen medio que puede escurrir (m3 )

A = Área de la cuenca (ha)
C = Coeficiente de escurrimiento (adimensional)
Pm = Precipitación media (mm)

Escurrimiento Máximo Instantáneo

El Escurrimiento Máximo Instantáneo se calcula para diseño de obras de excedencia

(vertederos, puentes y otros), pudiendo estimarse para diferentes periodos de retorno por
el Método Racional, Método Racional Modificado o por el método de Numero de Curvas
SCS (USA).

El Método Racional consiste en utilizar los valores de intensidad de la lluvia para

Página 15
diferentes periodos de retorno y el área de drenaje para estimar los escurrimientos
máximos instantáneos. La fórmula básica del Método Racional es:

Q=C.I.A
Donde:
Q = Caudal máximo expresado en m3/s
C = Coeficiente de escurrimiento.
I = Intensidad de la precipitación concentrada en m2/s en un período igual al
tiempo de concentración tc.
A = Área de la cuenca hidrográfica en m2.
Este método solo es aplicable en cuencas muy pequeñas, menores de 5 km2.

La modificación propuesta al método racional consiste en utilizar la lluvia máxima

precipitación en 24 horas para diferentes períodos de retorno, en lugar de la intensidad de
la lluvia. La siguiente ecuación corresponde al método racional modificado:

Qp = CLA / 360

Donde: Qp = Escurrimiento máximo instantáneo (m3 /s)

C = Coeficiente de escurrimiento
L = Lluvia máxima en 24 horas para un período de retorno dado (mm)
A = Área de drenaje (ha)
360 = Factor de ajuste de unidades.

Método de Numero de Curva del SCS (USA), se usa para estimar el escurrimiento medio
por evento y el máximo instantáneo utiliza el método de Numero de Curva, el cual utiliza
los datos de precipitación por evento o la precipitación máxima para un periodo de retorno
deseado y el máximo potencial de retención del agua del suelo. La siguiente ecuación
corresponde al método de Número de Curva (NC):

Q = (P- 0.2S)2 / (P+(0.8S))

Q > 0 Sí 2.0 S ≤ P
si no Q = 0

Página 16
Donde: Q = Escurrimiento medio (mm).
P = Precipitación por evento (mm).
S = Retención máxima potencial (mm).
El potencial máximo de retención de agua del suelo (S) depende de las condiciones del
suelo, vegetación y manejo del cultivo. Entonces es factible relacionarlo con el Número de
Curva, las cuales son función de los factores antes mencionados. Entonces el potencial
máximo de retención (S), se obtiene de la siguiente ecuación:

S = 25400/CN – 254

Donde: S = Potencial máximo de retención (mm).

CN = Curvas numéricas (adimensional).

2.2.8 NUMERO DE CURVA CN

”El CN es el número de curva de escurrimiento determinado para las diferentes

combinaciones hidrológicas suelo-vegetación. Es función del tipo de suelo, pendiente,
vegetación, porcentaje de cobertura, grado de humedad, temperaturas precedentes, etc.
Mediante este número, se determina el volumen de escorrentía como resultado de una
precipitación caída en un área dada”; Juan C Mecca (2008).

El número de curva, se estima inicialmente de acuerdo a la metodología elaborada por el

Servicio de Conservación de Suelos de los Estados Unidos de América, y luego se ajusta
durante el proceso de calibración de los valores utilizados para la trasformación de la
precipitación en escorrentía.

2.2.9 CAUDAL

Es la cantidad de volumen en metros cúbicos de agua que circula por un cauce, ducto o
canal por unidad de tiempo. Esto puede ser sobre o bajo la superficie terrestre, y que llega
a una corriente para finalmente ser drenada hasta la salida de la cuenca. El caudal se
mide en una estación de aforo por un aforador apropiado de acuerdo a la magnitud del
caudal.

Página 17
2.2.10 AVENIDA

Una avenida es una crecida del rio o elevación del nivel de un curso de agua,
significativamente mayor que el flujo medio del rio. Durante la crecida, el caudal de un
curso de agua aumenta en tales proporciones que la capacidad hidráulica del río puede
resultar insuficiente para contenerlo en particular cuando estas son torrenciales.

” Las avenidas torrenciales son fenómenos naturales que, aunque tienen una ocurrencia
relativamente baja, por su gran potencial destructivo, son de particular interés para
evaluar los riesgos naturales a los que se encuentra sometida una determinada cuenca.
Entre las razones que existen para que este tipo de evento sea tan dañino están su
naturaleza casi impredecible, la rapidez a la cual ocurre, su corta duración y su largo
período de retorno, así como su distribución poco uniforme en el espacio y el tiempo” J.
González, O. Chávez, and M. Hermelín (2005).

2.2.11 PERIODO DE RETORNO

"El período de retorno o de recurrencia (T) es el intervalo medio expresado en años en el

que un valor extremo alcanza o supera al valor (x)", al menos una sola vez" (Elías y Ruiz,
1979). Así, si la precipitación máxima en 24 horas para un período de retorno de 50 años
en el observatorio de San Calixto es de 60,3 mm/24 horas, significa que es posible que un
fenómeno de 60,3 mm de precipitación en 24 horas se repita o sea superado por lo
menos una vez en 50 años.

2.2.12 ESTACIONES HIDROMÉTRICAS

” La estación hidrométrica es un lugar fijo en una sección del rio donde se realizan un
conjunto de operaciones que permiten determinar el caudal circulante en momento y
tiempo determinado. Las estaciones hidrométricas se encuentran ubicadas en zonas de
fácil acceso donde el lecho del rio es estable. Es decir que no se erosiona, sedimente o
inunda fácilmente y donde el tramo del cauce sea lo más recto posible (50 metros aguas
arriba y aguas abajo de la estación para evitar distorsión en la información obtenidas”;
Vera Hernandez y Luis Edilberto (2000-2001).

La clasificación de una estación hidrométrica está basada en función al instrumental e

Página 18
implementación con que cuenta cada una de las estaciones hidrométricas. Si la estación
cuenta con un Limnimetro, se denomina Estación Limnimétrica. Si la estación cuenta
además del limnimetro con un limnigrafo se denomina Estación Limnigráfica.

2.3 MODELOS HIDROLÓGICOS

”Un modelo hidrológico es una representación física o matemática confiable de los

procesos hidrológicos reales que ocurren en una cuenca, cuyo objeto es cuantificar los
caudales entrantes a un área específica, provenientes desde sectores localizados
topográficamente en áreas más altas en tiempo y en forma que ellos ocurren”. Manuales
de HEC-HMS (US Army Corps of Engineers). Los modelos hidrológicos son ampliamente
utilizados para la transformación de la precipitación en escorrentía y determinación de la
capacidad hidráulica de canales y cauces de ríos, así como la simulación de afectaciones
e inundaciones por efectos de fuertes precipitaciones.

2.3.1 MODELADO HIDROLÓGICO HEC-HMS

El programa Hydrologic Modeling System (HEC-HMS), es producto de más de treinta y

cinco años de investigaciones del Hydrologic Engineering Center (HEC) del U.S. Army
Corps of Engineers. Es un programa de simulación hidrológica tipo evento, lineal,
semidistribuido, desarrollado para estimar los escurrimientos de salida en una cuenca ó
varias subcuencas, a partir de la transformación de una lluvia histórica ó hipotética,
aplicando para ello algunos de los métodos de cálculo de hietogramas de diseño,
pérdidas por infiltración, flujo base y conversión en escorrentía directa, mediante un
sistema que integra diferentes métodos hidrológicos, ya sea concentrados o distribuidos
para el tránsito de caudales.

”El modelo se compone de tres módulos básicos que definen en su conjunto el proyecto
de simulación de la cuenca:
- Módulo de precipitación: en el que se define una precipitación hipotética o histórica
para ser utilizada junto con el modelo de precipitación.
- Módulo de la cuenca: en el cual se define gráficamente la distribución de los
elementos hidrológicos (subcuencas, tramos de tránsito, uniones y otros). En cada
uno de ellos, se asignan atributos y parámetros. Los parámetros que utiliza este
modelo son: las pérdidas iniciales, la transformación de lluvia en escurrimiento y el

Página 19
 método del flujo base.
- Módulo de control de especificaciones: en el que se incluyen las fechas y hora de
inicio y fin de la simulación e intervalos de tiempo para realizar el cálculo del
hidrograma de escurrimiento”. (U.S. Army Corps of Engineers, 2007).

El programa incluye modelos por separado para representar los distintos componentes de
la escorrentía: modelos que computan el volumen escurrido, modelos de escorrentía
directa, modelos de flujo base y modelos de flujo de canales. Existen cinco grupos
básicos de información que deben suministrarse a HEC-HMS para efectuar las
simulaciones:
a - Información acerca de la precipitación histórica o de diseño.
b - Información acerca de las características del suelo.
c - Información morfométrica de las subcuencas.
d - Información hidrológica del proceso de transformación de lluvia en escorrentía.
e - Información hidráulica de los tramos de canal y de las capacidades de los
embalses (método de tránsito).

”El modelo Hydrologic Engineering Centers - Hydrology Modeling System (HEC-HMS)

está diseñado para simular los procesos de escorrentía, la precipitación de los sistemas
de cuencas dendríticas, también está diseñado para ser aplicable en una amplia gama de
áreas geográficas con el fin de resolver el mayor número posible de problemas. Los
hidrogramas producidos por el programa se utilizan directamente o en colaboración con
otros programas para los estudios de disponibilidad de agua, drenaje urbano, la previsión
de flujo, impacto sobre una futura urbanización, depósito de vertedero de diseño,
reducción de daños por inundaciones, la regulación de inundación y la operación de los
sistemas”. United State Army Corps of Engineers: HEC-HMS (2007).

El software HEC-HMS es un sistema de modelado generalizado capaz de representar

muchas cuencas diferentes, un modelo de la cuenca se construye mediante la separación
del ciclo hidrológico en segmentos manejables y la construcción de límites en tomo a la
cuenca de interés. Cualquier masa o flujo de energía en el ciclo puede ser representado
por un modelo matemático, en la mayoría de los casos, el modelo de varias opciones está
disponible para la representación de cada flujo, cada modelo matemático incluido en el
programa es adecuado en diferentes ambientes y bajo condiciones diferentes, hacer la

Página 20
elección correcta requiere un conocimiento de la cuenca, los objetivos del estudio
hidrológico y el juicio de ingeniería.

2.3.2 MODELO SWAT

“El Soil and Water Assessment Tool (SWAT), es un programa de modelamiento

hidrológico diseñado por el Departamento de Agricultura de los Estados Unidos en
conjunto con la Universidad de Texas”; Arnolt et al, (1990). Este modelo permite simular la
producción de agua y sedimentos en cuencas hidrográficas. “El SWAT se basa en un
balance hídrico para determinar la entrada, salida y almacenamiento de agua en la
cuenca. Para el modelamiento, la cuenca hidrográfica se divide en pequeñas subcuencas
con el fin de mejorar la exactitud de los cálculos, adicionalmente el SWAT trabaja por
unidades de respuesta hidrológica (HRU) obtenidas del cruce de los diferentes tipos de
suelo (textura), con el uso y cobertura del suelo”; Guzmán Edwin (2003).

La simulación de la Hidrología de la cuenca puede separarse en dos grandes divisiones:

la primera división es la fase terrestre del ciclo hidrológico que controla la cantidad de
agua y sedimento transportados al canal principal por cada sub-cuenca. La segunda
división es la del agua o la fase de rutina, la que puede definirse como el movimiento del
agua, sedimentos, etc., a través de la red del canal hasta el sitio de descarga de la
cuenca. El ciclo hidrológico simulado por el SWAT está basado en la ecuación del balance
hídrico:

t=i
SW t = SW + ∑(Ri - Q¡- ETi - Pi – QRi)
t=0

Donde:
SW t es la cantidad final de contenido de agua en el suelo (mm).
SW es la cantidad inicial de contenido de agua en el suelo (mm).
Ri es la cantidad de lluvia caída en el periodo de análisis (mm).
Qi es la cantidad de escorrentía (mm).
ETi Evapotranspiración (mm)
Pi Percolación (mm).
QRi Flujo de retorno (mm).

Página 21
2.3.3 MODELO HIDRÁULICO

“En hidráulica, el término modelo corresponde a un sistema que simula un objeto real
llamado prototipo, mediante la entrada de cierta información que se procesa y se presenta
adecuada para emplearse en el diseño y operación de obras de ingeniería. Un modelo
físico a escala reducida es una representación a escala del objeto real o prototipo, y
cumple ciertas condiciones matemáticas definidas”; Escuela de Ingeniería de Antoquia
(2000).

Los modelos hidráulicos unidimensionales que incorporen el efecto de remanso de agua

(típico de áreas de poca pendiente) son en principio adecuados para modelar las
estructuras de un sistema de control modular. Un modelo hidráulico unidimensional toma
como información de entrada los caudales evaluados por el modelo hidrológico. Con ellos
es posible calcular valores muy aproximados y físicamente consistentes de altura de
agua, velocidad, área inundada de cada sección de la estructura en función del tiempo.
(US Army Corps of Engineers, HEC-RAS).

2.3.4 MODELO IBER

"IBER es un modelo numérico de simulación de flujo turbulento en lámina libre en régimen

nopermanente, y de procesos medioambientales en hidráulica fluvial. El rango de
aplicación de IBER abarca la hidrodinámica fluvial, la simulación de rotura de presas, la
evaluación de zonas inundables, el cálculo de transporte de sedimentos y el flujo de
marea en estuarios"; IBER Modelización Bidimensional del Flujo en Lamina Libre en
Aguas Poco Profundos Manual de Referencia Hidráulico 07-06-2010. El modelo IBER
consta actualmente de 3 módulos fundamentalmente: un módulo hidrodinámico, un
módulo de turbulencia y un módulo de transporte de sedimentos.

"Todos los módulos trabajan sobre una malla no estructurada de volúmenes finitos
formada por elementos triangulares o cuadriláteros:
- El módulo hidrodinámico, que constituye la base de IBER resuelve las ecuaciones de
aguas someras bidimensionales promediadas en profundidad también conocidas
como 2D Shallow Water Equations (2D-SWE) o ecuaciones de St. Venant
bidimensionales. Dichas ecuaciones asumen una distribución de presión hidrostática y
una distribución relativamente uniforme de la velocidad en profundidad. La hipótesis

Página 22
 de presión hidrostática se cumple razonablemente en el flujo en ríos, así como en las
corrientes generadas por la marea en estuarios. Asimismo, la hipótesis de distribución
uniforme de velocidad en profundidad se cumple habitualmente en ríos y estuarios,
aunque pueden existir zonas en las que dicha hipótesis no se cumpla debido a flujos
locales tridimensionales o a cuñas salinas.

- El módulo de turbulencia permite incluir las tensiones turbulentas en el cálculo

hidrodinámico, pudiéndose utilizar para ello diferentes modelos de turbulencia para
aguas someras con diferente grado de complejidad. En la versión actual se incluyen
un modelo parabólico, un modelo de longitud de mezcla y un modelo k-ε (de
turbulencia).

- El módulo de transporte de sedimentos resuelve las ecuaciones de transporte de

fondo y transporte turbulento en suspensión, calculando a partir del balance de masa
de sedimento la evolución de la cota de fondo"; IBER Modelización Bidimensional del
Flujo en Lamina Libre en Aguas Poco Profundas Manual de Referencia Hidráulico 07-
06-2010.

2.4 SISTEMAS DE INFORMACIÓN GEOGRÁFICA (SIG)

El Centro Nacional de Información Geográfica (NCGI) de Estados Unidos define al SIG

como "Un sistema de hardware software y procedimientos elaborados para facilitar la
obtención, gestión, manipulación, análisis, modelado, representación y salida de datos
espacialmente referenciados, para resolver problemas complejos de planificación y
gestión".

Los sistemas de información geográfica (SIG), son manipulados por software o programas
informáticos y constan de las herramientas y funcionalidades necesarias para almacenar,
analizar y visualizar información geográfica.
Los componentes principales de un programa informático que permita construir un SIG
son los siguientes:
- Sistema de manejo de base de datos.
- Interfaz gráfica con los usuarios para el fácil acceso a las herramientas.
- Herramientas para captura y manejo de información georeferenciada.

Página 23
 - Herramientas para soporte de consultas, análisis y visualización de datos
geográficos.

2.4.1 FUNCIONES DE UN SIG

Ingreso de datos

El ingreso de datos se refiere a todas las operaciones por medio de las cuales los datos
espaciales de mapas, sensores remotos y otras fuentes son convertidos a un formato
digital. Entre los diferentes dispositivos comúnmente utilizados para esta operación están
los teclados, digitalizadores, barredores electrónicos y terminales interactivos o unidades
de despliegue visual Visual Display Units (VDU). Dado su costo relativamente bajo,
eficiencia y facilidad de operación; la digitalización es la mejor opción de ingreso de datos.

Se deben ingresar dos tipos diferentes de datos al SIG: referencias geográficas y

atributos. Los datos de referencias geográficas son las coordenadas (sea en términos de
latitud y longitud o columnas y líneas, en una determinada proyección cartográfica) que
fijan la ubicación de la información que se está ingresando.

Los datos de atributos asignan un código numérico a cada casilla o conjunto de

coordenadas y a cada variable, sea para representar los valores actuales (p.e., 200 mm
de precipitación, 1.250 metros de elevación) o para connotar tipos de datos categóricos
(usos del terreno, tipo de vegetación, etc.). La rutina de ingreso de datos requiere una
cantidad considerable de tiempo, ya sea el ingreso manual con teclado, digitalización, o
por barrido electrónico.

Almacenamiento de datos.
Almacenamiento de datos se refiere al modo como los datos espaciales son estructurados
y organizados dentro del SIG, de acuerdo a la ubicación, interrelación y diseño de
atributos.

Manipulación y procesamiento de datos.

La manipulación y procesamiento de datos se hace para obtener información útil de los
datos previamente ingresados al sistema. La manipulación de datos abarca dos tipos de
operaciones: (1) operaciones para eliminar errores y actualizar conjuntos de datos

Página 24
actuales y (2) operaciones que hacen uso de técnicas analíticas para dar respuesta a
preguntas específicas formuladas por el usuario. El proceso de manipulación puede ser
desde una simple sobreposíción de dos o más mapas, hasta una extracción compleja de
elementos de información dispares, de una gran variedad de fuentes.

La producción de datos se refiere a la exhibición o presentación de datos empleando

formatos comúnmente utilizados incluyendo mapas, gráficos, informes, tablas y cartas,
sea en forma impresa o como imagen en pantalla, o como un archivo de textos
trasladables a otros programas de cómputo para mayor análisis.

2.4.2 SISTEMAS DE INFORMACIÓN GEOGRÁFICA APLICADO A PLANOS DE

INUNDACIÓN

En la actualidad se han desarrollado herramientas que ayudan a estudiar la topología de

una zona con problemas de inundación, como son los SIG; Los SIG son una nueva
tecnología que permiten gestionar y analizar la información espacial y que surgió como
resultado de la necesidad de disponer rápidamente de información para resolver
problemas y contestar a preguntas de modo inmediato.

“Para la aplicación de resultados a la solución de problemas específicos, se combinan los

SIG con modelos hidrológicos e hidráulicos; los SIG tienen la ventaja de representar
aplicaciones como el de permitir utilizar los modelos hidrológicos desde una arquitectura
SIG, en la que gran parte de la información que precisan (datos de topología y cartografía)
le es suministrada por ésta. Igualmente, los datos resultantes del cálculo son
representados en este mismo entorno SIG, permitiendo análisis tanto gráficos como
alfanuméricos potentes e intuitivos” (Ayesa, 1993).

”La evaluación del riesgo por inundaciones mediante los SIG permiten: contar con una
visión de conjunto y multivariedad del riesgo, establecer relaciones espaciales y vincular
distintos tipos de información, contar con información digital de consulta directa, realizar
actualizaciones que respondan al dinamismo del problema”; Salgado (2005).

Base de Datos Espacial

Una base de datos espacial es una colección de datos referenciados espacialmente, que

Página 25
actúan como un modelo de la realidad; en el sentido de que ella representa una serie o
aproximación de fenómenos. Esos fenómenos seleccionados son considerados
suficientemente importantes para ser representados en forma digital. “La representación
digital puede ser para el pasado, presente o futuro” (IGAC, 2005).

Existen tres tipos básicos de estratos y diferentes combinaciones posibles entre ellos:
polígonos (llanuras de inundación, áreas de peligro de deslizamientos), líneas (líneas de
fallas, ríos, redes eléctricas) y puntos (epicentros, ubicación de pozos, instalaciones
hidroeléctricas, estaciones pluviométricas). La selección del tipo de estrato correcto para
una base de datos, depende de los usos anticipados y de la escala y resolución de los
datos de la fuente. Para distintos casos las áreas de inundación podrían estar
representadas por líneas a lo largo de los ríos; por polígonos en mapas a escala de
1:10.000.

Estructura de almacenamiento de datos geográficos.

“El componente espacial de un dato geográfico se refiere a una entidad geométrica sobre
la cual se observa un fenómeno. Se puede codificar de dos maneras distintas: el modelo
vectorial y el modelo raster“ Giovanni P. Guartatanga (2013); Ver figura 2.5 Modelo de
datos raster y vectorial.

Figura 2.5 Modelo de datos raster y vectorial (fuente web gabrielortiz.com)

Estos dos tipos de estructuras especiales son soportadas por los Estándares de

Página 26
Estructuras de Intercambio de Información Geográfica Digital cuya sigla en ingles es
DIGEST, estos datos son los que se listan a continuación: Dato Vector (en varios niveles
de topología) y Dato Ráster. Cada uno de estos tipos de datos es un diferente método de
representar un mapa o carta.

a) Dato vectorial.
Las propiedades constituyen el criterio de diferenciación de las entidades o elementos
geográficos que existen en el mundo real. Esas entidades son representadas por medio
de objetos de la base de datos, así por ejemplo en nuestro caso, un rio (la entidad) puede
representarse en un Sistema de Información Geográfica (SIG) mediante una linea (el
objeto).

b) Dato ráster.
Proporciona una aproximación basada en objetos elementales (celdas), que pueden
agruparse para constituir objetos complejos que representan el mundo real, como ser:
- Un punto: Se representa mediante una celda
- Una Linea: Se representa mediante una sucesión de celdas alineadas.
- Un polígono: Se representa mediante una agrupación de celdas contiguas.

“Existen dos grandes fórmulas para estructurar la información real en un sistema

informático, estas son el modelo vectorial y el modelo raster. A su vez ambas pueden
cubrir otros tipos de modelos que tienen más relación con la forma de concebir el espacio
que con el modelo lógico utilizado. Los modelos a los que nos referimos son el orientado a
capas y el orientado a objeto. Por último, hasta ahora se ha estado hablando de modelos
bidimensionales, es decir aquellos que no contemplan las tres dimensiones. La
incorporación de una tercera dimensión, es posible a partir de los Modelos Digitales de
Terreno (MDT) que pueden ser implementados tanto desde modelos raster como
vectoriales“; Giovanni P. Guartatanga (2013).

Modelamiento de datos.

Es aparente que una interpretación del mundo es necesaria, la que debe ser
suficientemente abstracta para que no sea afectada por la dinámica del mundo (los
pequeños cambios), y debe ser suficientemente robusta para poder representar cómo los

Página 27
datos y el mundo se relacionan. Una herramienta como esta es llamada modelo de datos,
el cual permite representar en forma más o menos razonable alguna realidad. “El modelo
de datos permite realizar abstracciones del mundo, permitiendo centrarse en los aspectos
macros, sin preocuparse de las particularidades; así nuestra preocupación se centra en
generar un esquema de representación, y no en los valores de los datos” (Obermeyer,
Nancy y Pinto, Jeffry, 1996).

Modelamiento Cartográfico.

El modelamiento cartográfico es el proceso por el cual a partir de un conjunto de datos

geográficos y empleando diferentes operaciones de análisis que siguen una secuencia
lógica, generamos nuevos datos geográficos, cuyas características se derivarán de las
que poseían los datos iniciales. Este es el proceso que empleamos para resolver
problemas del mundo real utilizando herramientas SIG, de tal forma que a partir de unos
datos de entrada obtengamos un resultado-solución a un problema planteado.
Generalmente los resultados de los procesos de modelado cartográfico suelen
representarse en forma de modelo cartográfico, que no es más que una representación a
modo de árbol en la que mediante la combinación, siguiendo una secuencia lógica, de
una serie de variables (mapas) de entrada se generan cartografías derivadas. “La forma
en que se combinan las variables, denominada como álgebra de mapas, nos ofrece
múltiples posibilidades mediante operaciones que permiten, a partir de uno, dos o más
mapas, obtener un nuevo mapa derivado. Estas operaciones pueden ser cualquier
herramienta de análisis disponible en un SIG. De esta forma disponemos de multitud de
operaciones, de las que podemos señalar entre otras operaciones aritméticas de
cualquier tipo, consultas, buffers, reclasificaciones, superposiciones, cálculos de
pendientes, etc.” (IGAC, 2005).

2.5 LEVANTAMIENTO TOPOGRÁFICO POR EL METODO INDIRECTO

El levantamiento topográfico por el método indirecto consiste en realizar las mediciones

de los objetos topográficos a cartografiar del espacio geográfico a través de sensor
remoto como; el vuelo aerofotogramétrico convencional o vuelo con DRON, equipado con
cámara digital para la captura de información geoespacial y post procesamiento de datos
para obtener el Modelo Digital de Elevación y el Mosaico de Ortofotos digitales

Página 28
fundamentalmente.

2.5.1 UNMANNED AERIAL VEHICLE UAV o DRON

Un Unmanned Aerial Vehcle se define con sus siglas en inglés UAV, se trata de un
vehículo volador no tripulado, también llamado “DRON”. Sus tipos, modelos y usos son
muy variados, para:
- Publicidad, Cine Documental, Eventos.
- Seguridad, Vigilancia de Fronteras, Monitoreo.
- Mantenimiento de Redes Eléctricas, Refinerías.
- Usos Militares, Cartografía, Topografía y Geodesia.
- Agronomía, Minería, Catastros, Control de Obras, GIS y otros.

Las UAVs o DRONs, deben ser pilotados AUTOMATICAMENTE, pero en la práctica, se

comportan en forma dual, es decir, pueden ser intervenidos y controlados por un operador
desde tierra. Los DRONs se clasifican básicamente en DRON de: Ala fija y Rotores.

DRON EBBE Sensefly

El eBee es el DRON profesional más fácil de usar para la captura de información puede
adquirir imágenes con una distancia de muestreo terrestre (GSD) de hasta 1,5 cm (0,6 in)
por píxel. Para hacerlo despegar, solo tiene que lanzarlo al aire. El DRON se encarga de
volar, capturar imágenes y aterrizar por sí solo. No obstante, puede alterar su plan de
vuelo siempre que quiera o hacerlo aterrizar de forma manual si es necesario. Gracias a
su estructura ultraligera el eBee tiene un peso de solo 700 g (1,5 lb), lo que reduce
enormemente su energía de impacto. También incluye una hélice de propulsión trasera
con sistema de seguridad y el innovador piloto automático de Sensefly, capaz de adoptar
una amplia variedad de comportamientos inteligentes a prueba de fallos. Ver figura 2.6;
EBEE Sensefly.

El sistema eBee contiene todo lo que necesita para comenzar a cartografiar: una cámara
RGB de alta resolución, baterías, un modem de radio y el software de planificación y
control de vuelo eMotion. Además, viene presentado en una sólida caja adecuada para su
transporte. Después solo se tiene que procesar y analizar las imágenes del eBee con un

Página 29
software profesional como el Pix4D Mapper.

Figura 2.6; EBBE Sensefly (www.google.com)

2.5.2 PUNTOS DE CONTROL O DE APOYO FOTOGRAMETRICO

Los puntos de control o de apoyo fotogramétrico son determinados mediante

procedimientos geodésicos de posicionamiento GPS y nivelación geométrica para
establecer las coordenadas (geográficas y planimétricas en UTM) y altura sobre el nivel
del mar.

Las coordenadas geodésicas de los puntos de control o puntos de apoyo fotogramétrico

para la georeferenciacion de las imágenes obtenidas con DRON son determinados
mediante observaciones GPS de precisión y enlazados con observaciones GPS
simultaneas a la red MARGEN-SIRGAS (estaciones permanentes del Instituto Geográfico
Militar IGM). Las alturas ortométricas de los puntos de control o puntos de apoyo
fotogramétrico son determinados mediante la nivelación geométrica a partir de los puntos
Bancos de Nivelación (Bench Mark BM) establecidos a nivel nación por el IGM.

Página 30
2.5.3 PIX4D MAPPER

”Pix4Dmapper es un software de procesamiento de imágenes que se basa en la

búsqueda automática de miles de puntos comunes entre las imágenes. Cada punto
característico encontrado en una imagen se llama un punto clave. Cuando se encuentra
que 2 puntos clave en 2 imágenes diferentes son iguales, se combinan los puntos clave.
Cada grupo de los puntos clave correctamente emparejados generarán un punto 3D.
Cuando hay una superposición alta entre 2 imágenes, el área común capturada es más
grande y más puntos clave se pueden combinar juntos. Cuantos más puntos clave haya,
se podrán calcular con mayor precisión los puntos 3D. Por lo tanto, la regla principal es
mantener una superposición alta entre las imágenes. Como el plan de adquisición de
imágenes tiene un gran impacto en la calidad de los resultados, es importante diseñarlo
cuidadosamente. Pix4Dmapper le permite al usuario procesar utilizando una Plantilla de
procesamiento”; Manual de Pix4D Mapper

Las plantillas con la etiqueta Rapid / LowRes producen resultados rápidos en baja
resolución que puede usarse en el campo como un indicador de qué tan bueno es el
conjunto de datos. Si el procesamiento Rapid / LowRes no da buenos resultados, muy
probablemente, el conjunto de datos no es adecuado y es necesario adquirir imágenes de
nuevo”. Manual de Pix4D Mapper.

El software Pix4D cuenta con módulo de planificación de vuelo Pix4D Capture y

procesamiento de datos Pix4D Mapper para obtener los productos; Modelo digital de
superficie (MDS), modelo digital de elevación (MDT) y Mosaico de Ortofotos.

2.6 CARTOGRAFIA DE RIESGOS

2.6.1 RIESGO NATURAL

”El riesgo natural es la mayor o menor probabilidad de exceder un valor específico de

consecuencias sociales o económicas en un sitio particular y en un tiempo determinado,
debido a la actividad de un proceso natural. Se obtiene de relacionar la amenaza con la
vulnerabilidad” Juan P. Celemin 2009.

Página 31
2.6.2 RIESGO DE INUNDACION

”Las inundaciones pueden ser consideradas como riesgo geoclimático, ya que su origen
se sitúa en la atmósfera, es decir, con las precipitaciones. Por otra parte, luego de las
lluvias, el movimiento del agua y los riesgos asociados se siguen por la geomorfología e
hidrología de la cuenca”; Ayala Carcedo 1988.

Los riesgos de inundación son cada vez más frecuentes, fundamentalmente por dos
factores: el primero la efectividad del cambio climático, con un posible aumento de la
frecuencia de las lluvias torrenciales y el potencial aumento del nivel del mar, como
consecuencia del calentamiento de la atmósfera; y el segundo factor, el impacto de las
actividades humanas, como las construcciones en los margines de los curso de los ríos y
obras hidráulicas de desviación y canalización de los ríos, o la construcción de puertos y/o
atracaderos sin medidas de evaluación y de corrección de su impacto ambiental. También
se incluye el factor de las actividades antrópicas, los procesos de desertización ocurridos
por talas masivas indiscriminadas de árboles, incendios y otras actividades contrarias a la
naturaleza.

”Durante un fenómeno de inundaciones, la predicción espacial es plenamente factible

para lo cual es necesario tener un registro de las precipitaciones y de los caudales de los
ríos y arroyos. A nivel temporal, la predicción es mucho más difícil y costosa. En algunas
cuencas, pueden emplearse valores umbrales de riesgo en función de registros
pluviométricos y proceder a la evacuación de la población afectada según las áreas
delimitadas en los mapas de riesgo”; Sanchez Miguel 2005.

2.6.3 CARTOGRAFIA DE RIESGOS DE INUNDACION

La cartografía de riesgos de inundación es un producto que representa gráficamente las

zonas de afectación por desbordes de los ríos. Generalmente se elabora la cartografía de
riesgos aplicando el proceso Multi Criterio Evaluación (MCE) en la que múltiples capas
son agregadas en un sistema de información geográfica (SIG) para obtener un solo mapa
de salida. Los criterios considerados pueden ser representados por medio de capas
(layers) de datos geográficos. El MCE es un método común para evaluar y combinar
diferentes criterios.

Página 32
”Para comenzar un MCE es necesario establecer el criterio. El mismo es la base para la
toma de una decisión y se caracteriza porque puede ser medido y evaluado. Cada criterio
puede ser representado por una imagen raster. Existen dos tipos de criterio: limitantes
(constraints) y factores (factors). El primero actúa como límite a las alternativas bajo
consideración. Se expresan en forma booleana (lógica), donde las áreas excluidas son
representadas con un 0, mientras que las consideradas tienen el valor de 1. Por otro lado,
un factor es un criterio que se utiliza para aumentar o disminuir la aptitud de una
alternativa que está siendo evaluada. Por lo tanto, está medida en una escala continua”;
Juan P Celemin 2009.

Página 33
 CAPITULO III
MATERIALES Y METODOLOGIA

3.1 MATERIALES

Los datos e información recopilada a distintos organismos son los siguientes:

- Modelo Digital de Elevación (MDT); Shuttle Radar Topography Mission SRTM 30m.
Este producto fue obtenido de la web UMSA (www.ciags.umsa.com).
- Datos pluviométricos de las estaciones meteorológicas ubicadas en la cuenca del Rio
La Paz (Norte), estos datos fueron obtenido de la página web del Servicio Nacional de
Meteorología e Hidrología (SENAMHI). Sistema SisMet : Precipitación máxima en 24
horas en mm. De la Estación San Calixto.
- Datos de Precipitación media anual de las estaciones meteorológicas de la cuenca del
Rio La Paz Norte de SENAMHI; Sistema SisMet.
- Datos de Precipitación Total anual de las estaciones del periodo 1917 – 2017 de
SENAMHI; Sistema SisMet.
- Datos mensuales de precipitación, temperaturas medias, evaporación, humedad
relativa y otros de SENAMHI; Sistema SisMet.
- Mapa Geológico en formato digital obtenido de la IDE Geobolivia.
http://geo.gob.bo/download/?w=GobDepartamento&l=lp_vi_Geomorfologia2000
- Mapa de uso de suelo urbano La Paz en formato digital obtenido de IDE Geobolivia
http://geo.gob.bo/download/?w=GobMunicipal&l=uso_suelo_urbanowgs84.
- Mapa de suelos de Bolivia en formato digital obtenido de la IDE Geobolivia:
http://geo.gob.bo/download/?w=onu&l=suelos.
- Mapa de Vegetación de Bolivia en formato digital obtenido de la IDE Geoboliva
http://geo.gob.bo/download/?w=mdpdd&l=vegetacion.
- Mapa de Cobertura y Uso Actual de la Tierra; en formato digital obtenido de la IDE
Geobolivia enlace http://www.vicetierras.gob.bo
- Software ArcGis y Software SWAT
- Software HEC-HMS
- Software Leica Geoffice
- DRON EBBE Sensefly
- Software PIX4D MAPPER
- Software IBER

Página 34
3.2 METODOLOGÍA

En este acápite, se describen las actividades a ser ejecutadas y los resultados que se
generaran en el desarrollo de la investigación utilizando los materiales como los datos e
información recopilada, la forma de llevarlas a cabo y los fines perseguidos en cada una
de ellas, para someter a prueba las hipótesis planteadas y ver si se ajustan a la realidad;
Ver figura 3.1 Metodología de Investigación.

Figura 3.1 Metodología de Investigación (elaboración propia)

Página 35
3.2.1 DELIMITACIÓN DEL ÁREA DE LA CUENCA Y SUB-CUENCAS DEL RIO LA PAZ
NORTE CON EL SOFTWARE SWAT.

Para este efecto se utiliza el modelo digital de terreno (MDT) obtenido en el marco del
proyecto Shuttle Radar Topography Mission (SRTM 30m), disponible en formato raster,
siendo la delimitación de la cuenca a partir del MDT automática. Procedimiento con la que
se ha delimitado la superficie del área de las sub-cuencas: Choqueyapu, Orkojahuira,
Irpavi, Achumani y Huyñajahuira. Esta actividad constituye la base de partida para la
descripción de las características biofísicas de la cuenca.

Delimitada espacialmente la cuenca y sub-cuencas, así como el corte del modelo digital
de elevación del área de investigación, se establece las características morfométricas de
la cuenca y sub-cuencas del sistema de drenaje como: Pendiente y Longitud del cauce,
Áreas, Tiempo de concentración y otros elementos requeridos para el modelamiento de la
hidrología de superficie.

3.2.2 CARACTERIZACIÓN BIOFÍSICA DE LA CUENCA Y SUB-CUENCAS.

Para este efecto se ha utilizado fundamentalmente la información geográfica temática en

formato digital como: Mapa de Cobertura y Uso Actual del Suelo y Mapa de Suelo;
complementariamente la información de: Mapa Hidrográfico, Mapa de Pendientes y Mapa
de Vegetación. La caracterización de las sub-cuencas permite establecer las
características cualitativas y cuantitativas en los aspectos biofísicos naturales y
construidas por el ser humano; predominantes existentes en el espacio de la cuenca y
sub-cuencas del rio la Paz (Norte).

3.2.3 CARACTERIZACIÓN CLIMÁTICA DE LA CUENCA DEL RIO LA PAZ NORTE.

Con los datos de la estación meteorología de San Calixto, se ha caracterizado las

condiciones climáticas de la cuenca del rio La Paz (Norte) como: Temperatura y
Evaporación media mensual, Presión barométrica (atmosférica), Humedad relativa,
Vientos. Asimismo, en base a los datos referidos a la Máxima Precipitación en 24 horas
de la estación San Calixto se deriva las curvas de Intensidad Duración y Frecuencias
(IDF) para 50 años de retorno.

Página 36
3.2.4 DETERMINACIÓN DE LA LLUVIA DEL PROYECTO CON TIEMPO DE RETORNO
50 AÑOS.

Utilizando la ecuación de la IDF derivada con los datos de la estación San Calixto
disponible en la Carrera de Topografía y Geodesia de la UMSA web
(www.ciags.umsa.com) y el tiempo de concentración calculado; se genera la lluvia del
proyecto a intervalos de 15 minutos. Estos datos constituyen la lluvia del proyecto para la
evaluación de la capacidad hidráulica del área de estudio.

3.2.5 DETERMINACIÓN DEL NUMERO DE CURVA.

Como método de evaluación de pérdidas o capacidad de retención de humedad de la

superficie de la tierra del área de la cuenca y sub-cuencas del rio La Paz (Norte).
Determina el NC con los resultados de la caracterización biofísica de las sub-cuencas
(cuantificación cuantitativa y cualitativa de cobertura del suelo) respecto a los parámetros
de Soil Conservance of USA, se estima el número de curva ponderado para cada una de
las sub-cuencas que aportan con cierto caudal en el puente encuentro entre loa ríos
Choqueyapu e Irpavi.

3.2.6 GENERACIÓN DEL MODELO HEC-HMS.

Con los valores de la lluvia del proyecto, parámetros morfométricos de las sub-cuencas y
el número de curva ponderado del área de las sub-cuencas se realiza el modelamiento
para la transformación de la lluvia en escurrimiento superficial (caudal) a efectos de
determinar el hidrograma de la cuenca para tiempo de retorno 50 años. Es decir, obtener
la máxima crecida del agua en el tramo puente encuentro entre los ríos Choqueyapu e
Irpavi aguas abajo hasta el encuentro con el rio Huayñajahuira a efectos de probar la
capacidad hidráulica del curso del rio.

3.2.7 LEVANTAMIENTO TOPOGRÁFICO EN DETALLE DEL CURSO DEL RIO

Del sector el puente encuentro de los ríos Choqueyapu e Irpavi hasta la confluencia con el
río Huayñajahuira; fundamentalmente por el método indirecto mediante de la ejecución de
las siguientes actividades:
- Toma de imágenes fotogramétricas con UAV (DRON)
- Determinación de Puntos de Control Fotoidentificables

Página 37
- Procesamiento de datos para la generación de Ortofotos y MDTs.

3.2.8 GENERACIÓN DE MDT PARA EL SECTOR DEL RIO OBJETO DE ESTUDIO.

A través del levantamiento topográfico con DRON y post-procesamiento de datos

fotogramétricos se genera el mosaico de ortofotos y el modelo digital de elevación que
refleja la geométrica del curso del rio objeto de estudio. Entonces con la información de
levantamiento topográfico reflejada en el modelo digital de terreno (MDT), se generan los
perfiles transversales para el modelamiento hidráulico del sector.

3.2.9 DETERMINACIÓN DE LOS PARÁMETROS MORFOMÉTRICOS DEL RIO.

Determinación de: Caudal base, Longitud del tramo, área cubierta por transversales,
coeficiente de rugosidad (coeficiente de Maning) y tipo de cobertura y otros elementos a
considerar en el modelamiento hidráulico.

3.2.10 GENERACIÓN DEL MODELO DE DINÁMICA FLUVIAL CON IBER.

Modelamiento Hidráulico del curso del rio; tramo entre el puente encuentro de los ríos
Choqueyapu e Irpavi hasta la confluencia con el río Huayñajahuira Utilizando como dato el
máximo caudal base obtenido con el modelo HEC-HMS y la geometría del curso de rio
reflejada en las secciones transversales obtenido a partir del modelo digital de terreno.

3.2.11 DETERMINACIÓN DE LA CAPACIDAD HIDRÁULICA DEL CURSO DEL RIO

En el sector puente encuentro de los ríos Choqueyapu e Irpavi hasta la confluencia con el
rio Huayñajahuira para una avenida con tiempo de retorno de 50 años. Se realiza el
análisis y evaluación de la avenida respecto a la capacidad hidráulica del indicado sector.

3.2.12 CARTOGRAFÍA DE RIESGOS DE DESBORDES, COTA Y VELOCIDAD DE

AGUA.

Conformación de un SIG con la información gráfica temática como uso de suelo y la

información geográfica de desbordes para diferentes situaciones Aplicando criterios,
restricciones y factores que inciden en el establecimiento de áreas de riesgo para la
elaboración de cartografía de riegos.

Página 38
3.2.13 PROPUESTA DE MEDIDAS DE MITIGACIÓN.

En base al SIG conformado con la información geográfica de los riesgos de desborde e

inundaciones identificados, proponer las medidas de mitigación inmediatas, de mediano y
de largo plazo a efectos de dar la sostenibilidad al uso actual del suelo y bienestar de los
habitantes de las áreas circundantes al tramo del rio objeto de estudio.

En base a los resultados obtenidos, se verificará el cumplimiento de la hipótesis.

Asimismo, se emitirán las conclusiones y recomendaciones correspondientes.

Página 39
 CAPITULO IV

RESULTADOS Y DISCUSIONES

4.1 CUENCA DE ESTUDIO

Para la delimitación del área de la cuenca y sub-cuencas del rio La Paz (Norte), se ha
realizado con el modelo digital de terreno (MDT) SRTM de 30 m. de resolución espacial y
el Software SWAT en la plataforma ArcGis (10.4) de forma automática, apoyado con la
información del mapa de jurisdicción municipal (COMLIT) y la información toponímica de
la cartografía nacional existente a escala 1:50000, que cubre el área de la cuenca y sub-
cuencas del rio La Paz (Norte).

4.1.1 MODELO DIGITAL DEL TERRENO

El modelo digital de terreno (MDT) constituye la información fundamental para la

delimitación del área de la cuenca y sub-cuenca hidrográfica de rio La Paz (Norte). Para
el propósito de la investigación se ha recopilado el MDT en formato raster SRTM de 30 m.
de resolución espacial. Ver figura 4.1; Modelo Digital del Terreno que cubre el área total
de la cuenca del rio La Paz (Norte); información base para la delimitación automática del
área de la cuenca y sub-cuencas del rio La Paz (Norte) y establecer los parámetros
morfométricos para el modelamiento y procesamiento de los datos meteorológicos en el
modelo de transformación de la precipitación en escurrimiento HEC-HMS, a efecto de
establecer flujo hidrográfico de la cuenca del rio la Paz (Norte).

El modelo digital de terreno (MDT) expresada en tres dimensiones (3D) ya sea en formato
vectorial o raster es la información básica fundamental para el manejo y gestión de
información hidrográfica de superficie de las cuencas, utilizada generalmente en la
realización de estudios hidrológicos de las cuencas, identificación y cuantificación del
sistema de drenaje hidrográfico, delimitación automática de la cuenca y sub-cuencas.
Asimismo, el MDT también es utilizado alternativamente en el procesamiento de datos
para la clasificación de suelos, zonificación climática y otros elementos temáticos que
hacen parte de la información requerida para múltiples aplicaciones, siempre y cuando no
exista suficiente información temática real recabada con trabajo de campo.

Página 40
Figura 4.1 Modelo digital del terreno de la cuenca del rio La Paz (Norte).

4.1.2 CARTOGRAFÍA

La cartografía nacional a escala 1:50000 (formato digital) con cobertura del área de la
cuenca del rio La Paz (Norte) y complementariamente la cartografía urbana escala
1:15000 (formato digital) existente del área urbana de la cuidad de La Paz, se han
utilizado como marco de referencia geográfica espacial y toponímica para la identificación

Página 41
de la cuenca y sub-cuencas del rio La Paz (Norte). La información altimétrica de la
cartografía digital existente solo se ha utilizado para relacionar la consistencia de la
información altimétrica generada a partir del modelo digital del terreno. Entonces, la
cartografía digital se ha constituido en una fuente de información geográfica y toponímica
para la delimitación e identificación de la cuenca y sub-cuencas (ríos Choqueyapu,
Orkojahuira, Irpavi, Achumani y Huayñajahuira).

La cartografía urbana digital a escala 1:15000 solo tiene cobertura del aérea urbana de la
ciudad de La Paz, razón por la que se ha utilizado de forma complementaria para
relacionar la información toponímica del área urbana y la caracterización biofísica de la
cuenca y sub-cuencas del rio La Paz (Norte).

4.1.3 DELIMITACIÓN DEL ÁREA DE LA CUENCA DEL RIO LA PAZ (NORTE)

La ubicación, forma y tamaño del área de la cuenca y sub-cuencas del rio La Paz (Norte),
se ha realizado mediante la delimitación automática del área a partir del modelo digital de
terreno (MDT SRTM de 30 m) con el Software Soil and Water Assessment Tool (SWAT)
en la plataforma ArcGis 10.4. Para la delimitación de la cuenca se ha definido como punto
origen (punto de salida de la escorrentía) del curso del río la Paz (Norte) el puente cruce a
Mallasa; Zona Aranjuez. Ver figura 4.2; Mapa de delimitación de la cuenca del rio La Paz
(Norte).

Figura 4.2 Mapa de delimitación de la cuenca y sub-cuencas del rio La Paz (Norte).

Página 42
Asimismo, se ha extractado del modelo digital del terreno con el software SWAT (Stream
definition, DEM based y stream área reducido 200 ha) la información gráfica de la
hidrografía de la cuenca y sub-cuencas del rio la Paz con suficiente densidad. Ver figura
4.3; Mapa hidrográfico de la cuenca del rio La Paz; posteriormente se ha realizado la
asignación de información toponímica de los ríos con la información de la cartografía
nacional a escala 1:50000.

Figura 4.3; Mapa hidrográfico de la cuenca del rio La Paz

Página 43
4.1.4 DETERMINACIÓN DE PARÁMETROS MORFOMÉTRICOS

Definida la red de drenaje a ser empleada en el modelo, se ha procedido a la

determinación de parámetros morfométricos tales como área de aporte, longitud del curso
principal, altura máxima y altura mínima. Información en base a la cual se determina para
cada sub-cuenca la pendiente, el tiempo de concentración y el tiempo de rezago
fundamentalmente.

El tiempo de concentración y de rezago son variables utilizadas extensivamente en el

diseño hidrológico para determinar la capacidad hidráulica máxima de diferente
estructura. El tiempo de concentración también es conocido como el tiempo de respuesta
o de equilibrio.

El tiempo de concentración es calculado a través de la siguiente fórmula de Témez, que

tiene la siguiente expresión:

Tc = 0.3 (L/(J0.25))0.76

Donde:
Tc = Tiempo de concentración en horas
L = Longitud del curso principal en Km.
J = Pendiente del curso principal en m/m

En cuanto al tiempo de rezago, se encuentran en la literatura la relación propuesta por el

“Soil Conservation Service” (SCS) USA, quien propone que es 0,6 veces el tiempo de
concentración.

Log T = 0.6 Tc

Aplicando las fórmulas mencionadas se procede al cálculo de los parámetros

morfométricos: Ver Tabla 4.1; Parámetros Morfométricos de la cuenca y sub-cuencas del
Rio La Paz (Norte).

Página 44
 Longitud TC Log time
Cuenca del Rio Superficie(Km2) (km) Pendiente (Hrs) (min)
La Paz 492.10 28.467 0.0714 6.02 216.72
Choqueyapu 140.53 24.267 0.0642 5.55 199.80
Choqueyaou 1 15.09 4.178 0.0340 1.69 60.84
Orkojahuira 91.68 23.189 0.0573 5.46 196.56
Irpavi 165.42 25.540 0.0630 5.64 203.04
Irpavi 1 0.75 1.319 0.0379 0.69 24.84
Achumani 62.28 15.200 0.0925 3.59 129.24
Huayñahuira 16.34 8.595 0.0904 2.30 82.81

Tabla 4.1 Parámetros Morfométricos de la cuenca y sub-cuencas del Rio La Paz (Norte).

Aplicando la anterior expresión a la cuenca del río La Paz, se determinó que el tiempo de
concentración de la cuenca es de 6.00 horas. Considerando que el curso principal de la
cuenca, atraviesa el área rural en una longitud de 18.13 km aproximadamente con un
desnivel de 1690 m, y fluye a través del área urbana a lo largo de 10. km
aproximadamente con una diferencia de altitudes de 680 m.

Es importante indicar, que el área de la cuenca del rio La Paz (Norte) está fuertemente
influenciado por la existencia de altas pendientes en la superficie de la tierra y la
existencia del área urbana de la ciudad de La Paz. Este último debido a la
impermeabilidad de la superficie de concreto y asfalto, el tiempo de concentración de las
sub-cuencas se reduce considerablemente en función a las superficies consideradas
como urbanas y no urbanas.

4.2 CARACTERIZACION BIOFISICA DE LA CUENCA DEL RIO LA PAZ (NORTE)

La cuenca del río La Paz (Norte) área de estudio se encuentra en los valles interandinos
de la cordillera de los andes; en la parte Sur de la jurisdicción del municipio de Nuestra
Señora de La Paz, el área de estudio tiene una superficie de 49847 hectáreas. El
Municipio de Nuestra Señora de La Paz limita al norte con el municipio de Guanay, Al Sur
con los municipios de Achocalla, Mecapaca y Palca, al Este con los municipios de
Caranavi, Coroico y Yanacachi, al Oeste con los municipios de Batallas, Pucarani y El
Alto; Ver figura 4.4; Mapa de Ubicación de la cuenca y sub-cuencas del rio La Paz (Norte).
La cuenca del rio La Paz (Norte) hidrológicamente está conformada por las sub-cuencas
que corresponden a los ríos Choqueyapu, Orkojahuira, Irpavi, Acumani y Huayñajahuira.

Página 45
Estos ríos son los principales afluentes que aportan al caudal hidrográfico del rio La Paz
(Norte).

Figura 4.4 Mapa de ubicación de la cuenca y sub-cuencas del rio La Paz (Norte).

Página 46
4.2.1 CARACTERIZACION HIDROLOGICA

El río La Paz nace en la confluencia de los rios Irpavi y Choqueyapu; pertenece a la

macro cuenca del Río Amazonas, Este río se constituye en las nacientes del río Boopi,
principal afluente del río Beni, que al confluir con los ríos Madre de Dios y Mamoré e
Itenez dan lugar al río Madera. La cuenca del río La Paz (Norte) hasta el sector de
Aranjuez, está conformada por las sub-cuencas de los ríos Choqueyapu, Orkojahuira,
Irpavi, Achumani y Huayñajahuira; la cuenca tiene una superficie de 498.4 km2, la
longitud del curso más largo es de 37 km y presenta un desnivel de 2175 m desde sus
nacientes a 5300 msnm aproximadamente hasta 3125 msnm en sector de Aranjuez
aproximadamente; arrastra gran cantidad de material en suspensión. Ver imagen 4.1 Río
Choqueyapo, Río Irpavi.

Imagen 4.1 Río Choqueyapo Río Irpavi

La topografía de la cuenca del río La Paz es sumamente accidentada, las laderas alta y
medias de la cuenca se caracterizan por presentar fuertes pendientes y quebradas
profundas, resultado de la acción de las lluvias en los suelos de baja cohesión; la parte
baja de la cuenca entre la unión de los ríos de La Paz y Huayñajahuira está compuesta
por relieves de pendiente baja y zonas de deposición o llanuras aluviales, sectores donde
es frecuente la acumulación del material de arrastre que transportan los ríos.
Principalmente los ríos de Irpavi, Achumani y Huañajahuira que anualmente acumulan
material de arrastre que disminuye la capacidad hidráulica del curso del rio La Paz (Norte)
requiriendo acciones de evacuación de áridos por el GAMLP para el control hidrológico.

Página 47
SUB-CUENCA DEL RÍO CHOQUEYAPU

El río Choqueyapu. Ver figura 4.4 Mapa de ubicación de la cuenca y sub-cuencas del rio
La Paz(Norte), nace en los cerros Wila Manquilizani, Charquerini y Chiar Kherini a alturas
mayores a los 5240 msnm aproximadamente con el nombre de Junitincut Jahuira. Ver
imagen 4.2 Nacientes del río Choqueyapu y Kaluyo; en su curso hacia el Sur recibe
aportes de afluentes que a su vez son el resultado de rebalses de varias lagunas, ya en
cercanías a la Mina Kaluyo, éste curso se denomina río Kaluyo hasta el sector de la
planta hidroeléctrica de Achachicala donde recibe el nombre de Choqueyapu. En este
tramo la cuenca está limitada por la Ceja, Achachicala y las serranías Condorjipiña; ya al
ingreso a la ciudad de La Paz el río Achachicala se denomina Choqueyapu.

Imagen 4.2; Nacientes del Río Choqueyapu, como Río Kaluyo

Junitincut Jahuira en ese sector.

El río Choqueyapu cursa un tramo aproximado de 10. km a través de embovedados hasta

el sector cruce de las avenidas Poeta y Hernando siles, a partir de este sitio el río discurre
a través de un vertedero (canal) hasta su confluencia primero con el río Orkojahuira en la
Curva de Holguín, luego confluye con el río Irpavi en la zona de Calacoto, a partir de esta
confluencia cambia de nombre a Río La Paz.

Los principales afluentes del río Choqueyapu en el área urbana por el margen derecho,
son los ríos Jacha, Apumalla, Carahuinchinca, San Pedro y Cotahuma, y por el margen
izquierdo los ríos Calle 8, Matadero, Zarzuela, 7 Enanos, Viscachani, Humahuaca,
Mejahuira, Zoqueri, Orkojahuira, Aguarani y Collpajahuira; entre los más importantes.

Página 48
SUB-CUENCA DEL RÍO ORKOJAHUIRA

El rio Orkojahuira. Ver figura 4.4 Mapa de ubicación de la cuenca y sub-cuencas del rio La
Paz(Norte), nace con el nombre de río Chuquiaguillo en el sector denominado La Cumbre;
cerro Jilaralpata a una altura de 5320 msnm aproximadamente; fluye en dirección Sur y
recibe aportes de la laguna Estrellani y otros afluentes, estos aportes son almacenados en
la represa Incachaca, aguas abajo el río Chuquiaguillo es alimentado por los rebalses de
la represa y por cursos de aguas provenientes de las serranías Almillanis, el cerro Khisa
Patilla y la serranía Condorjipiña por el Oeste y, del cerro Kinkillosa y la serranía Cuchilla
Chuquiaguillo por el Este. En este tramo el río discurre en dirección Suroeste; una vez
que atraviesa la zona de Chuquiaguillo cambia de nombre el rio a río Orkojahuira y en el
sector de la Curva de Holguín confluye con el río Choqueyapu. Ver imagen 4.3 Río
Orkojahuira, Confluencia con Choqueyapu.

Imagen – 4.3; Río Orkojahuira. confluencia con el río Choqueyapu.

La longitud del curso principal del rio Orkojahuira es de 23 km aproximadamente, el

desnivel desde el punto más alto hasta su confluencia con el río Choqueyapu es de 1715
m; desde el punto de vista morfológico, esta sub-cuenca tiene una pendiente de
accidentadas a accidentadas medias y presencia de área urbana. La respuesta a las
fuertes precipitaciones es rápida. En el área urbana, los afluentes más importantes son
los ríos Minasa, Fábrica de Fósforos, Chapuma, Guitarrani, y por la margen izquierda, los
ríos Venecia, Huallpajahuira.

Página 49
SUB-CUENCA DEL RÍO IRPAVI

El rio Irpavi; Ver figura 4.4 Mapa de ubicación de la cuenca y sub-cuencas del rio La
Paz(Norte), nace en los cerros Jamphaturi, Pumaya Pata, Mikhaya y Serkhe Kollu, a
alturas mayores a los 5440 msnm aproximadamente; el curso principal denominado Rio
Mikhaya fluye hacia el Oeste hasta desembocar en la laguna Ajuan Khota; a cuatro
kilómetros aguas abajo se encuentra la represa de Jampaturi. A partir de este punto el río
se denomina rio Anta, aguas abajo recibe por el Este el aporte del río Palcoma, luego de
esta confluencia el río toma el nombre de Rio Karpani; siguiendo la misma dirección de
flujo, cuando el río atraviesa el sector de Chicani cerca de la Zona Khallapa, adquiere el
nombre de rio Irpavi, a partir de este sector el cauce de este río se amplía y atraviesa una
extensa zona urbana. Aguas abajo, el río Irpavi confluye con el río Achumani manteniendo
su nombre, posteriormente se une al río Choqueyapu en el puente encuentro en la zona
de Calacoto de la ciudad de La Paz. En las siguientes fotografías: Ver imagen 4.4 Río
Palcoma y Río Irpavi; se puede ver las características del río Irpavi que arrastra gran
cantidad de bloques de piedra y áridos.

Imagen -4.4; Río Palcoma, aguas abajo de la presa Río Irpavi antes de su confluencia
Jampaturi. con el río Choqueyapu

La longitud del curso principal es de 27 km aproximadamente, su punto más alto se

encuentra a 5440 msnm y el punto de confluencia con el Achumani está a 3211 msnm.

Página 50
SUB-CUENCA DEL RÍO ACHUMANI

El rio Achumani; ver figura 4.4 Mapa de ubicación de la cuenca y sub-cuencas del rio La
Paz(Norte), nace a alturas mayores a los 4944 msnm aproximadamente en los cerros
Chachacomani, con el nombre de Kellumani, antes de su ingreso en la zona urbana éste
río se encuentra canalizado, aguas arriba de la urbanización Jardín confluye con el río
Huayllani también canalizado en este sector. Ver imágenes 4.5 Rio Huayllani y 4.6 Rio
Achumani. A partir de este punto este curso se denomina Achumani; aguas abajo recibe
el aporte del río Jakonta, que si bien tiene un área de aporte reducida su capacidad de
arrastre de material sólido es muy significativa debido a sus elevadas pendientes y
presencia de material suelto. Los afluentes más importantes del río Achumani son, por el
Oeste, las quebradas Jerenko, Karate y Koani, y por el Este, el río alemán y el Jillusaya.
El río Achumani, en la zona de Calacoto confluye con el río Irpavi.

Imagen 4.5.- Río Huayllani aguas arriba Río Huayllani en el sector canalizado, antes de
canalización. su confluencia con el río Kellumani.

Imagen 4.6.- Río Achumani, atravesando la zona Río Achumani antes de su confluencia con
Urbana. el río Irpavi.

Página 51
La longitud de su curso principal es de 15 km aproximadamente, las pendientes de este
curso son elevadas en la parte alta de la cuenca y las lluvias son de respuesta rápida.
Provocan fuerte erosión en las cabeceras, siendo el material erosionado transportado a la
parte baja de la sub-cuenca.

SUB-CUENCA DEL RÍO HUAYÑAJAHUIRA

El rio Huayñajahuira¸ ver figura 4.4 Mapa de ubicación de la cuenca y sub-cuencas del rio
La Paz(Norte), nace en los cerros Patarani y Pachajaya con el nombre de quebrada
Negra, fluye en dirección Oeste y recibe los aportes de pequeños ríos de las zonas de
Ovejuyo, Cota Cota, Los Pinos, San Miguel y La Florida, limitadas por los cerros Zapatuni,
Mamani y Chipaña. Al atravesar la zona de Cota Cota este río se encuentra canalizado
hasta el sector de La Florida donde confluye con el río La Paz.

La longitud de su curso principal del rio es de 8.59 km aproximadamente, tiene sus

pendientes que indican que se desarrolla en topografía accidentada a accidentada media,
siendo por tanto de respuesta también rápida a las precipitaciones. En las fotografías Ver
imagen 4.7 Rio Huayñajahuira, es posible observar las características de este cauce, que
tiene un alto potencial de transporte de material sólido, el mismo que es depositado al
ingresar al río La Paz.

Imagen 4.7; Río Huayñajahuira, Río Huayñajahuira, sector canalizado; área urbana.

En el área urbana de la ciudad de La Paz, el río Auquisamaña es el afluente más

importante del río Huayñajahuira.

Página 52
4.2.2 CARACTERIZACION GEOLOGICA

La cuenca del rio La Paz (Norte) está ubicada en tres unidades geomorfológicas
diferenciadas: La Cordillera, el Altiplano y el Valle interandino de La Paz. En el área de la
cuenca se encuentran las formaciones geológicas; relacionadas con: Formaciones de
Zócalo, Formaciones del Altiplano y las Formaciones Contemporáneas y depósitos
recientes, como grandes unidades genéricas:
- Formaciones del Zócalo (Aranjuez y Sica Sica).
- Formaciones del Altiplano (Gravas Milluni, Purapurani, Calvario y La Paz).
- Formaciones Contemporáneas y depósitos recientes (Terrazas Antiguas Miraflores e
Irpavi, Movimientos de remoción en masa antigua y reciente, agrupadas en 4 épocas
de erosión

Formaciones Rocosas del Zócalo

La cuenca del rio La Paz, es una depresión intermontano constituida por rocas plegadas y
falladas del Ordovícico, Silúrico, Devónico, Triásico, Jurásico, Cretácico (Formaciones del
zócalo), constituida en las estribaciones del pie de la Cordillera Real, que fue luego
rellenada por sedimentos del Terciario Tardío y del Cuaternario (Formaciones del
Altiplano).

La secuencia que compone las rocas del zócalo se inicia en el Ciclo Cordillerano
(Ordovicico Superior - Devónico Superior), que de base a tope comprende flujos
sedimentarios de la (Fm) Cancañiri, lutitas oscura de la (Fm) Uncía, areniscas y lutitas de
la (Fm) Catavi (Estas dos formaciones en la parte NE y en la parte sur de la cuenca) y
areniscas de la (Fm) Vila Vila (En la parte NE de la cuenca, continuando un proceso
regresivo representado por lutitas, limolitas y escasas areniscas de la (Fm) Sica Sica (en
la parte Sur de la cuenca de La Paz).

Posteriormente, durante la Fase Incaica del ciclo Andino (Triásico Inferior al presente),
continúa la depositación en ámbitos intermontanos con fases fluviales y de cono aluvial
correspondientes a la (Fm) Aranjuez, en la parte Sur de la cuenca. Durante el Triásico
Terminal se produce la intrusión del Batolito Taquesi, conformado por rocas graníticas (En
el borde NE de la cuenca). El análisis de esta fase sugiere que su desarrollo está

Página 53
vinculado al funcionamiento de un frente de cabalgamientos. Es en esta fase donde se
forma el relieve actual de las rocas del zócalo. Los movimientos de la fase Quechuana
reactivan la sedimentación, después de una discontinuidad litoestratigráfica, originando
los depósitos fluvio-lacustres de la (Fm) La Paz (Paleógeno Tardio y Neógeno),
compuestos por arcilitas, areniscas y conglomerados, que sobreyacen a las formaciones
del zócalo en forma discordante. También se producen manifestaciones magmáticas de la
Muela del Diablo.

Formaciones del Altiplano

En discordancia angular sobre las rocas del zócalo yace la Formación La Paz, que en su
parte inferior consiste de arenas, arcillas y limos de carácter lenticular y color grisáceo y la
parte superior compuesta por capas rojizas de gravas, arenas, limos y arcillas y también
un horizonte de cinerita (ceniza volcánica), muy continua. Esta Formación aflora en toda
la cuenca de La Paz y en el escarpe del Altiplano.

Seguidamente viene una serie de formaciones de origen glaciar formados por un till
(Formación Calvario) compuesto por una mezcla de pedrones, arenas y arcillas y
sedimentos fluvioglaciares, que constituye la Formación Purapurani y se caracteriza por
su buena selección y espesor, y constituye el mejor acuífero de toda la cuenca. Esta
unidad aflora en el escarpe del Altiplano, junto con las gravas Milluni, en el Cerro Calvario
y en la parte alta de Ovejuvo camino a Palca.

Formaciones Contemporáneas y Depósitos Recientes

Dentro del valle la cuenca de La Paz, de edad reciente y/o actual se hallan sedimentos de
menor extensión y espesor como son las terrazas fluvioglaciares como la grava
Miraflores, Irpavi, Obrajes, La Florida, Següencoma. Dentro de los depósitos
contemporáneos se contabilizan una serie de movimientos de remoción en masa que han
afectado y afectan a diferentes sectores de la ciudad de La Paz. Ver Figura 4.5; Mapa
Geológico.

Página 54
Figura 4.5; Mapa Geológico de la cuenca La Paz (Norte).

4.2.3 GEOMORFOLOGÍA

Debido a la configuración geológica de la cuenca de La Paz, los procesos

geomorfológicos cobran especial relevancia porque los terrenos que la conforman son

Página 55
esencialmente suelos en diversos grados de consolidación y diagénesis, abarcando
aproximadamente 60 % de su superficie, mientras que el restante 40 % corresponde a
rocas de diferentes edades geológicas.

Los factores erosivos y litológicos, como formadores del relieve, dominan las formas del
relieve. Las formas del relieve en los suelos de la cuenca, son de tipo glacis en las partes
altas y en el Altiplano, junto a las montañas y cordilleras. Farellones en las cimas, con
inclinaciones esporádicamente mayores que 45º; pendientes empinadas en las partes
intermedias, con inclinaciones de 20 a 30º y pendientes más suaves en las partes bajas e
inferiores. Valles glaciales y fluviales que cortan a los anteriores. Lechos y planicies
aluviales. Terrazas y abanicos aluviales. Planicies fluvioglaciales y planicies testigo en
diferentes alturas. Superficies de deslizamiento y flujos de barro, derrumbes o
deslizamientos secos.

La red fluvial, aunque localmente puede tener un diseño rectangular, como es la ladera
Oeste, en el conjunto general es de diseño dendrítico, pero con proceso de jerarquización,
como se evidencia por los cinco ríos importantes que constituyen las cinco sub-cuencas
que colectan el flujo de todos los afluentes que corren por los interfluvios, estas cinco sub-
cuencas son: Choqueyapu, Orkojahuira, Irpavi, Achumani y Huañajahuira.

Entre estos ríos existen los interfluvios o terrenos altos clasificados geomorfológicamente
en el marco del proyecto ZONISIG como; Área urbana, Pie de monte, Montañas
alta/media y Serranías alta/media que los separan y que constituyen la superficie de los
suelos utilizados para la edificación de la ciudad. Ver figura 4.6; Mapa Geomorfológico.

En estas áreas se desarrollan gran parte de las formas de relieve indicadas

anteriormente. Numerosos riachuelos, quebradas y torrentes abarrancan la superficie de
los interfluvios y de la ladera Oeste, debido a la naturaleza impermeable o casi
impermeable de su superficie, constituida principalmente por material limo-arcilloso donde
la escorrentía, en forma de riada, es concentrada.

El avenamiento dendrítico de la red fluvial, debido a las fuertes pendientes modela un

relieve quebrado e irregular, con algunas planicies aisladas elevadas en la parte central,
como Huaripampa y Wila Pampa por el Norte, Pampahasi al Este y planicies de fondo de

Página 56
valle como la de Miraflores hasta el Estado Mayor, la zona central desde la cervecería
hasta el cuartel de San Jorge, Obrajes y Alto Obrajes en parte, Alto Irpavi y Alto Achumani
y la planicie aluvial de Calacoto.

Choqueyapu, Orkojahuira e Irpavi, en cuyos nacientes montañosas y serranías se han

depositado morrenas laterales y frontales a lo largo de su curso. Los valles glaciales han
derivado en valles de río o fluviales que cortan a los primeros y se extienden hacia los
cursos bajos de los valles actuales, donde se desarrollan actualmente procesos de
profundización de cauces, debido a que se encuentran en un ciclo juvenil de erosión.

Los ríos Achumani y Huañajahuira se inician como arroyos de pendiente por las
precipitaciones pluviales y por alumbramiento de aguas subterráneas, encontrando aguas
abajo de su recorrido, terrenos más suaves y erosionables, lo mismo que el curso medio e
inferior del río Irpavi. Los ríos Achumani e Irpavi en su curso inferior, han podido
desarrollar un ciclo de erosión maduro, hasta situarse por debajo de su nivel de equilibrio,
con predominio de socavamiento lateral y agradación de su perfil longitudinal, propiciando
una tendencia a la elevación de su lecho por exceso de material grueso de arrastre, que
no puede ser evacuado por la escorrentía de agua.

El río Kellumani en particular, que juntamente con el río Huayllani constituyen las
nacientes del río Achumani, se inicia en la serranía de Wila Pampa, haciendo leves
incisiones en la planicie ligeramente inclinada, en la que además alumbran algunas
vertientes, cuando transmonta su escarpe entra en terrenos de fuerte pendiente pie de
monte con predominio de erosión concentrada que forma profundos barrancos. A medida
que desciende se unen otros afluentes que también erosionan, formando una amplia
cuenca con gran aporte de sedimentos gruesos provenientes de la grava. Cerca de su
confluencia con el río Huayllani entra en terrenos gravosos, que también aporta una
importante cantidad de gravas al lecho del río Achumani.

Página 57
Figura 4.6; Mapa de Geomorfológico (fuente ZONISIG 2000)

4.2.4 PENDIENTE

En el análisis del relieve el concepto de pendiente desempeña un papel esencial. Así, las
planicies descritas más arriba tienen ligera inclinación hacia el Sur o Suroeste (Planicie de
Calacoto). Estas planicies se han extendido hasta el Sudoeste de la cuenca, impactando

Página 58
las rocas de la serranía de Aranjuez; encima del barrio Amor de Dios se observan mantos
de estas terrazas actualmente cubiertas por torrentes de barro. La pendiente actual entre
Miraflores Alto (Plaza Villarroel) y la zona Amor de Dios es de 7.5 % aproximadamente, lo
mismo que la planicie de la Zona Central entre la cervecería y la zona Amor de Dios.

La planicie de Pampahasi es un caso especial, corresponde a un periodo interglaciar

anterior; constituyéndose en un cerro testigo a la altura de la culminación del frente, por
ello ha podido subsistir sin erosionarse durante largo tiempo. Tiene una pendiente de 11
%, que sumada a la acción de los acuíferos subterráneos que afloran en sus taludes, la
hace todavía inestable en sus bordes, como es el caso de Kupini en su talud Sur e
innumerables deslizamientos actuales y recientes en ambas vertientes – Chuquiaguillo e
Irpavi. Una pendiente está en equilibrio cuando, en las condiciones actuales, no tiende a
aumentar ni a disminuir sensiblemente su estabilidad.

Las pendientes de las vertientes de las sub-cuencas (vertientes en el sentido de laderas),

sean interfluvios o no (caso ceja de El Alto-río Choqueyapu), se da en terrenos no
rocosos, o sea compresibles en diverso grado. Se considera concentrada cuando tiende a
dividir la vertiente en múltiples barrancos, como producto de precipitaciones intensas para
que la altura de agua caída por unidad de tiempo, no pueda filtrarse. Tal el caso de la
mayoría de los riachuelos y torrentes de las laderas.

La escorrentía concentrada tiene lugar en las vertientes desnudas. Una cobertura vegetal
y sobre todo una cobertura forestal que se interponga entre la precipitación violenta y el
suelo, disminuye la posibilidad de abarrancamiento y la reduce casi a la nulidad.

Las pendientes de las laderas o vertientes alcanzan inclinaciones de 25 a 66 %

aproximadamente al pie de los farellones, disminuyendo a menos de 45 %
aproximadamente en la parte media y baja de su recorrido hacia los ríos o cauces
principales indicados. Ver figura 4.7; Mapa de Pendientes (expresado %) de la cuenca La
Paz (Norte).

Página 59
Figura 4.7 Mapa de pendientes (Elaboración propia a partir del MDT SRTM 30m)

Los interfluvios entre las sub-cuencas y además la ladera Oeste, tienen farellones que
alcanzan inclinaciones superiores a 45 %, llegando en algunos casos a superar el 70 %.

Página 60
Estas pendientes están sujetas al fenómeno de la solifluxión, en el sentido del descenso
de una formación en forma de desprendimientos y de barro, por la imbibición de los
terrenos gravosos por encima de capas más impermeables como son las formaciones La
Paz y la Cinerita Chijini. Esta imbibición de los terrenos es el resultado del alumbramiento
de aguas subterráneas de los acuíferos como el del Altiplano en la ladera Oeste,
acuíferos alrededor de Pampahasi, entre la grava Miraflores y formación La Paz en los
ríos Choqueyapu y Orkojahuira, fundamentalmente. Esta imbibición también es producida
por el agua de lluvias torrenciales que se precipitan sobre la cuenca y por las filtraciones
de aguas servidas de urbanizaciones asentadas en las partes altas, por rotura de tuberías
y alcantarillas. La imbibición de las unidades gravosas superiores provoca en las unidades
infrayacentes, como la formación La Paz, un cierto estado de plasticidad o aún de
liquidez, produciéndose un deslizamiento fangoso lento o flujos de barro respectivamente.

4.2.5 CARACTERIZACION DE LA COBERTURA VEGETAL

El mapa de vegetación muestra los diferentes tipos de vegetación existentes en la región,

para su elaboración recientemente se ha incorporado nueva información de campo, con
una base más homogénea para la interpretación, como los mosaicos GEOCOVER y el
modelo digital de elevación de la NASA; en este contexto, se ha vuelto a recorrer la mayor
parte del país para la toma de numerosos inventarios adicionales de la vegetación en el
campo, siendo una herramienta que ofrece información precisa sobre el potencial forestal
ecológico y las orientaciones para la protección de la vegetación. El Mapa de Vegetación
de Bolivia a escala 1:250 000 fue elaborado por W. Ferreira y G. Navarro en el año 2007

El mapa de vegetación utilizada para la caracterización de la cuenca del rio La Paz,

describe las características de las plantas predominantes en la cuenca. Fue elaborado en
varios años y en diferentes regiones y zonas del país: Pando, Yungas y sur del Beni
(Corredor Amboró-Madidi), Parque Nacional “Kaa-Iya” del Gran Chaco, TCO Monteverde
(Chiquitanía), Lomerío (Chiquitanía), Reserva de Vida Silvestre Ríos Blanco y Negro
(Chiquitanía), norte de los parques nacional Isiboro-Sécure y Amboró.

La leyenda del mapa de vegetación de Bolivia en general está basada en la agrupación

de las grandes regiones fisiográfico-biogeográficas: Cordillera, Amazonía, Beni,
Chiquitanía. Entonces el presente documento contiene una leyenda y una breve

Página 61
explicación de las características de la cobertura de vegetación del área de la cuenca del
Rio La Paz Norte. Ver figura 4.8; Mapa de vegetación de la cuenca del Rio La Paz (Norte).

Figura 4.8 Mapa de vegetación de la cuenca del rio La Paz

Descripción de las unidades de vegetación de la cuenca del Rio La Paz (Norte).

Página 62
- p2'+p3a+p3b+p26; p2'= Pajonales altoandinos de la Puna Húmeda sobre glacis y
piedemontes con suelos profundos + p3a= Pajonales higrofíticos. Herbazales
graminoides dominados por gramíneas y ciperáceas amacolladas, cespitosas y
rizomatosa.

- p26= Antrópico: cultivos extensos, barbechos y prados extensos.

- p26+p17+p3+p8+p5f; p26= Antrópico: cultivos extensos, barbechos y prados

extensos + p17= Matorrales y pajonales altimontanos de la Puna y Altiplano xerofíticos
sobre suelos bien drenados. Matorrales xeromórficos con gramíneas y cactácea.

- (p16d) +p17d+p23a+p24a+p25; (p16d) = Arbustal-cardonal puneño interandino del

sur de La Paz. Serie de Oreocereus fossulatus-Kageneckia lanceolata. 3200 – 3650
m. Cuenca interandina xérica de los ríos La Paz y Luribay.

- p1+p3+p8+p27r; p1= Vegetación subnival de la Puna subhúmedo-húmeda (pajonales,

prados, bofedales y vegetación geliturbada). Conjunto de tipos de vegetación del piso
bioclimático criorotropical que incluye: pajonales muy bajos o prados sobre los suelos.

- p3a= Pajonales higrofíticos. Herbazales graminoides dominados por gramíneas y

ciperáceas amacolladas, cespitosas y rizomatosas, acompañadas de hierbas
hemicriptofíticas rosuladas y postradas. Se desarrollan sobre suelos hidromorfos.

4.2.6 CARACTERIZACION DE LA COBERTURA DE USO DE LA TIERRA

La cobertura de la tierra es la cubierta fisca y/o biofísica que se observa sobre la

superficie terrestre, incluye la vegetación y elementos antrópicos (cultivos,
construcciones, infraestructura, así como afloramiento de roca, suelo desnudo, cuerpo de
agua, salares, aflorecencias salina, acumulación de arena, etc), se refiere a todo aquello
que ocupa un lugar sobre la superficie, cualquiera que fuese su origen físico o biológico,
natural o artificial tal cual se observa directamente o a través de sensores. “El uso actual
de la tierra se refiere a la actividad humana presente al momento de efectuar la
observación y/o la descripción de sus características en una época determinada sin tomar
en consideración su uso potencial o futuro”; Memoria Técnica del mapa de cobertura y

Página 63
uso actual de la tierra de Bolivia (2010).

La conceptualización del mapa indicada anteriormente, constituye la base de partida para

la realización de la caracterización de las unidades de cobertura y uso actual de la tierra
existentes en la cuenca y sub-cuencas del Rio la Paz (Norte).

De acuerdo al Mapa y Memoria Técnica del mapa de cobertura y uso actual de la tierra de
Bolivia (2010), el área de las sub-cuencas; Choqueyapu, Orkojahuira, Irpavi, Achumani y
Huayñajahuira están cubiertas por las siguientes unidades de cobertura de la tierra:
- Bosque denso Andino Montano Húmedo
- Campos de Nieve
- Cuerpos y Cursos de Agua
- Estructura urbana y rural, ciudades y pueblos
- Herbazal Graminoide Vivaz Altimontano Pluvionebular
- Herbazal Graminoide Vivaz Bofedal Puneño
- Herbazal Graminoide Vivaz Puna semiárida
- Plantaciones Forestales
- Vegetación Dispersa Vivaz en pedregales Puna Altoandina Semiarida

La cuenca del rio La Paz (Norte) está cubierta mayoritariamente por; Herbaceas
Graminoides amacollada vivaz sinusia arbustiva microfila, Dispersa de arbustos en
sustrato rocoso gramíneas y forvias efímeras, Matorral siempre verde y herbáceas
graminoide amacollada vivaz neblina y garua y una buena parte de área urbana. Ver
figura 4.9; Mapa Cobertura y Uso actual del Suelo de la cuenca del rio La Paz (Norte).

Las características de las sub-cuencas correspondientes a los ríos Irpavi, Achumani y

Huayñajahuira con superficie Herbaceas Graminoides amacollada viva sinusia arbustiva
microfila, hacen que estos ríos arrastren cantidad de bloques de piedra y material en
suspensión. Asimismo, indicar que la cobertura de área urbana de la cuenca del Rio La
Paz (Norte) como; piso pavimentado y techo de las viviendas, impermeabilizan el área
urbana lo cual produce escurrimiento inmediato de la precipitación.

Página 64
Figura 4.9 Mapa Cobertura y Uso actual del Suelo de la cuenca del rio La Paz (Fuente COBUSO
del Viceministerio de Tierras 2010).

4.3 CARACTERIZACION CLIMATICA Y CURVAS IDF PARA 50 AÑOS DE RETORNO

El clima en la cuenca del río La Paz se ha clasificado como mesotérmico con invierno

Página 65
seco frío (Cwb) según Koppen, que considera para su clasificación la temperatura y
humedad principalmente, y como clima de tierras altas (H) según la clasificación de
Trewartha – Robinson, que toma en cuenta la temperatura media ambiente y el déficit de
agua del terreno. Los parámetros meteorológicos que determinan el clima en la cuenca
del río La Paz son su latitud 16.5º al Sur del Ecuador, su altura que varía de 5400 msnm a
3125 msnm y su configuración topográfica.

4.3.1 TEMPERATURA

En función a la ubicación de la cuenca del Río La Paz en la zona tórrida ecuatorial, el

grado de radiación ultravioleta es elevado por la baja densidad atmosférica, situación que
se acentúa aún más por la altitud a la que se encuentra la cuenca. Las temperaturas,
siguen un régimen estacional, aunque con variaciones de una estación a otra de 4 a 5º C;
más importantes suelen ser las variaciones de temperatura entre el día y la noche. En la
siguiente tabla se presentan datos de temperatura de la estación SENAMHI – San Calixto
Paz (1995-2016), Ver Tabla 4.2: Temperatura Estación San Calixto La Paz.

TEMPERATURA MEDIA °C
JUL AGO SEP OCT NOV DIC ENE FEB MAR ABR MAY JUN

9,36 10,15 11,04 12,05 12,95 12,66 12,09 12,05 12,05 11,58 10,71 9,76
TEMPERATURA MAXIMA MEDIA °C
16,06 16,80 17,32 18,06 18,93 18,01 16,83 16,78 16,95 17,08 17,18 16,44
TEMPERATURA MINIMA MEDIA °C
2,71 3,52 4,75 6,08 7,01 7,34 7,34 7,33 7,17 6,08 4,25 3,10

Tabla 4.2 Temperatura Estación San Calixto La Paz (datos de SENAMHI)

También, es importante indicar que debido a las diferencias altitudinales que se presentan
en la cuenca, desde 5440 mnsm en la ladera Oeste y 3211 msnm en la zona de Calacoto,
existe también un gradiente de temperaturas en función a la altura. En el estudio
“Situations Météorologique et Variations Climatologiques en Bolivie-Rapport Provisoire”
(Ronchal, 1985), se indica que por la variación altitudinal en la cuenca del río La Paz, las
temperaturas medias tienen un gradiente térmico altitudinal medio de 0.64 °C/100 m, por
lo cual se establece que a medida que se baja en altitud las temperaturas medias son
más benignas, es así que generalmente existe una diferencia de temperaturas de 2.3°C
entre el centro de la ciudad y la zona Sur, y de 5.1°C entre ésta última y El Alto.

Página 66
4.3.2 EVAPORACIÓN

Los valores de evaporación en la cuenca son elevados, propios de las regiones de altura,
donde existen amplitudes térmicas tan marcadas en el transcurso diario y efectos de poca
nubosidad en gran parte del año; de acuerdo a datos de SENAMHI 2004 y 2012 la
evaporación media anual es de 2.7 mm por mes. Tabla 4.3; Evaporación media mensual:
Estación San Calixto La Paz.

EVAPORACION MEDIA MENSUAL (mm) CON DATOS DEL 2004 Y 2012

JUL AGO SEP OCT NOV DIC ENE FEB MAR ABR MAY JUN

3,05 2,80 2,80 2,95 2,70 2,75 1,45 1,20 1,85 2,50 3,95 3,60

Tabla 4.3; Evaporación media mensual: Estación San Calixto La Paz (Datos de SENAMHI).

4.3.3 PRESIÓN BAROMÉTRICA (ATMOSFÉRICA)

“En la ciudad de La Paz la presión barométrica tiene un valor promedio de 493.5 mm de

mercurio (657.8 milibares). Siendo estable estos valores, sin embargo, suele variar como
máximo 10.1 mm Hg; esta presión representa el 65% de la presión atmosférica al nivel de
mar que se considera normal”; Segaline Nieto – Cabré (1988). En las cabeceras de la
cuenca, la presión barométrica anual media alcanza 626.6 mb; datos registrados en la
estación de El Alto - AASANA, valor que es concordante con su inversa con la altitud de la
ciudad.

4.3.4 HUMEDAD RELATIVA

La disponibilidad de humedad es variable por factores estacionales y locales, siendo

mayor en verano por el aporte de masas de aire húmedo provenientes del Norte y Noreste
(Yungas), y también por la presencia de los glaciares cercanos de la cordillera Oriental,
que alimentan ríos y pequeños cauces que desembocan en lagunas o reducidas áreas de
inundación que se convierten en superficies potenciales de evaporación. “A nivel diario,
las diferencias de temperatura entre el día y la noche determinan notables variaciones en
la humedad relativa, es así que, en los meses de junio, julio y principios de agosto, por
ejemplo, la humedad relativa media durante el día es del 20%, mientras que por la noche
la humedad puede superar el 90%”; Segaline Nieto – Cabré (1988). A nivel mensual,
como se indicó anteriormente, las variaciones de este parámetro son estacionales,

Página 67
presentándose los valores anuales mayores en los meses de diciembre a abril. En la tabla
siguiente se muestran los valores medios de la estación San Calixto. Tabla 4.4; Humedad
relativa media mensual: estación San Calixto.

HUMEDAD RELATIVA MEDIA MENSUAL(%)

JUL AGO SEP OCT NOV DIC ENE FEB MAR ABR MAY JUN

47,23 51,51 57,82 56,39 55,44 62,26 67,56 68,35 66,53 62,04 49,68 46,17

Tabla 4.4 Humedad relativa media mensual: Estación San Calixto La Paz (datos 1995 – 2016 de SENAMHI).

4.3.5 VIENTOS

Con respecto a los vientos, éstos están condicionados por la configuración orográfica, la
baja regulación térmica y la poca densidad del aire; en la ciudad de La Paz. “El régimen
de calmas es notorio por las noches entre las 21:00 y 24:00 horas y es normal que se
extienda hasta la salida del sol, debido a la irradiación nocturna. Según datos del
Observatorio San Calixto”; Segaline Nieto-Cabré (1988), los vientos predominantes
durante el año son de dirección Sudeste, presentándose pequeñas variaciones entre los
meses de mayo a julio, donde los vientos predominantes son del Oeste.

4.3.6 PRECIPITACIONES

Las precipitaciones en la cuenca del río La Paz por lo general provienen de la región del
norte y son resultado de los rebalses de la cordillera oriental que ingresan a la cuenca y
van circulando desde la parte alta a la baja; las masas de aire ingresan a la cuenca a
través de dos aberturas que se forman entre los nevados Huayna Potosí-Cerro Rosario y
Cerro Rosario-Mururata, sectores próximos a los límites Norte y Este de la cuenca. “La
primera abertura corresponde a las cabeceras de los ríos Choqueyapu y Orkojahuira, y la
segunda a las partes altas de las cuencas de los ríos Irpavi y Achumani” (Ramirez, 1995).

Las masas de aire que ingresan por las aberturas mencionadas, sobretodo en eventos de
baja intensidad, producen precipitaciones independientes, que afectan a la cuenca de los
ríos Choqueyapu y Orkojahuira (A-B) en el primer caso y/o en el segundo a Irpavi,
Achumani y Huayñajahuira (C-D-E), tal como lo muestra en la figura 4.10. Ingreso de
masas de aire. Por tanto, con menor frecuencia se producen eventos simultáneos en toda

Página 68
la cuenca.

Figura 4.10; Ingreso de masas de aire en eventos de baja intensidad (Ramirez, 1995).

Con relación a eventos extremos, éstos al ser de origen convectivo, se caracterizan por su
corta duración y alta intensidad. Se puede observar también la variabilidad de las
precipitaciones en función a la topografía (origen orográfico), es así que las estaciones
ubicadas en puntos altos cercanos a la cordillera registran una precipitación anual más
alta que las estaciones ubicadas a cotas inferiores.

Las precipitaciones anuales en las diferentes estaciones se encuentran alrededor de los

500 a 600 mm/año aproximadamente. Las lluvias que se producen en la cuenca tienen un
carácter estacional, concentrándose las máximas precipitaciones en el período húmedo,
entre los meses de noviembre a marzo, periodo en el cual se produce aproximadamente
el 77% de la precipitación anual; el 17% de las precipitaciones se presenta generalmente
en las épocas de transición en los meses de abril, septiembre y octubre, mientras que en
la época seca suele precipitar solamente el 6.5% de la lluvia anual restante. En la figura
4.3, se muestra las precipitaciones mensuales de julio a junio, considerando registros de
estaciones representativas de las sub-cuencas de los ríos Choqueyapu, Orkojahuira,
Irpavi, Achumani y Huayñajahuira. Ver figura 4.11; Régimen de precipitación mensual en
la cuenca del rio La Paz Norte.

Página 69
Figura 4.11; Régimen de precipitaciones mensual; cuenca del río La Paz (datos SENAMHI).

En general, el mes de enero es el que presenta mayor número de días con lluvia.

4.3.7 CURVAS DE INTENSIDAD DURACION Y FRECUENCIA

La información pluviométrica y pluviografica de las estaciones meteorológicas de

SENAMHI ubicadas en el área de la cuenca del rio La Paz, fueron recopiladas del periodo
1961 a 2016, fundamentalmente en lo que se refiere a la máxima precipitación en 24
horas de la estación meteorológica referente San Calixto (Anexo -1) y la precipitación
media anual de las estaciones meteorológicas ubicadas en la cuenca del rio La Paz
(Anexo 2). Información necesaria para la formulación de las curvas de Intensidad
Duración y Frecuencia, modelamiento del escurrimiento hidrográfico y capacidad
hidráulica del vertedero natural del puente encuentro de los ríos Irpavi y Choqueyapu
aguas abajo.

Para obtener la intensidad de la lluvia (precipitación), tiempo de duraciones y lluvia

acumulada, se aplicó el Método Analítico (Aparicio 1997) los siguientes valores;
Coeficientes de discretización definidos en base a información pluviográfica máxima
precipitación en 24 horas registrada en la estación San Calixto. Con los datos
discretizados, se plantea la alternativa de obtener una ecuación que genere los
parámetros para la obtención de las curvas IDF a través de regresión lineal, de modo de
extrapolar la ecuación generada para toda la cuenca del rio La Paz (Norte) o a zonas que
carezcan de registros pluviográficos y que se encuentren relativamente cerca.

Página 70
Con los datos de la máxima precipitación en 24 horas registrada en la estación San
Calixto la Carrera de Topografía y Geodesia; Facultad Tecnológica de la UMSA ha
generado las curvas de Intensidad – Duración – Frecuencia; para un periodo de retorno
de 50 años disponible en la web (www.ciags.umsa.com).

Estación San Calixto (www.ciags.umsa.com);

I =122.0761 x T 0.149024/D 0.61629
Donde; I intensidad, T años de retorno y D tiempo de duración

4.4 DETERMINACION DE LA LLUVIA DEL PROYECTO

Para la generación de la lluvia del proyecto se ha utilizado la ecuación de la Intensidad-

Duración-Frecuencia para 50 años de retorno de la estación meteorológica de San
Calixto.
En el siguiente cuadro se muestra la intensidad de la lluvia calculada en base a datos de
la Estación San Calixto (ver tabla 4.5); Intensidad de la lluvia para 50 años de retorno.

No. Tiempo Intensidad Lluvia acumulada Incremento Lluvia Lluvia mm Ordenado

(mm)
1 15 41,195135 10,2987837 10,2987837 0,471
2 30 26,8716325 13,4358162 3,1370325 0,491
3 45 20,9292256 15,6969192 2,26110298 0,512
4 60 17,5283958 17,5283958 1,83147661 0,537
5 75 15,2759353 19,0949191 1,56652324 0,564
6 90 13,6521572 20,4782358 1,38331675 0,596
7 105 12,414692 21,725711 1,24747516 0,633
8 120 11,4337921 22,8675843 1,14187331 0,677
9 135 10,633112 23,9245021 1,05691778 0,730
10 150 9,9645096 24,911274 0,98677194 0,794
11 165 9,39597587 25,8389336 0,92765964 0,877
12 180 8,90531737 26,7159521 0,87701848 0,987
13 195 8,47661361 27,5489942 0,83304211 1,142
14 210 8,09811743 28,343411 0,79441677 1,383
15 225 7,76095274 29,1035728 0,76016177 1,831
16 240 7,45827537 29,8331015 0,72952869 3,137
17 255 7,18471457 30,5350369 0,70193545 10,299
18 270 6,93599083 31,2119588 0,67692184 2,261
19 285 6,70864789 31,8660775 0,65411871 1,567
20 300 6,49986073 32,4993036 0,63322616 1,247

Página 71
 21 315 6,30729551 33,1133014 0,61399779 1,057
22 330 6,12900554 33,7095305 0,59622904 0,928
23 345 5,96335288 34,289279 0,57974858 0,833
24 360 5,80894845 34,8536907 0,56441168 0,760
25 375 5,66460574 35,4037859 0,55009515 0,702
26 390 5,52930451 35,9404793 0,53669344 0,654
27 405 5,4021622 36,4645948 0,52411551 0,614
28 420 5,28241103 36,9768772 0,51228238 0,580
29 435 5,16937964 37,4780024 0,50112519 0,550
30 450 5,06247813 37,968586 0,49058355 0,524
31 465 4,96118584 38,4491903 0,48060431 0,501
32 480 4,86504134 38,9203307 0,47114045 0,481
33 495 4,77363405 39,3824809 0,4621502 0,462

Tabla 4.5; Intensidad de la lluvia con 50 años de retorno; Estación San Calixto.

4.5 DETERMINACIÓN DE NÚMEROS DE CURVA (NC)

El modelo de transformación de precipitación en escorrentía HEC-HMS y otros modelos

utilizados frecuentemente para el cálculo de caudal pico, requieren de la determinación
del Número de Curva (NC) ponderado o estimado para la cuenca. Valor que permite
estimar parte de las pérdidas o capacidad de retención de la humedad de la tierra para
transformar la precipitación total en escurrimiento.

La determinación de CN se puede realizar en base a la identificación de la combinación

de los grupos de suelos, pendiente del suelo y la cobertura vegetal y uso descrita en la
caracterización de la cuenca de estudio. Los factores se pueden estimar a través de la
superposición de los mapas de: Mapa hidrográfica, Mapas fisiográfico, Mapa de
pendientes, Mapa de cobertura y uso actual de la tierra y Mapa de suelos
fundamentalmente para conformación del mapa de complejos hidrológicos que permite
identificar las unidades y estimar el número de curva (NC) correspondiente.

El mapa geológico utilizado en la caracterización fisiográfica del presente estudio

identifica de forma aproximada las unidades de suelo en la cuenca del rio La Paz (Norte).
Sin embargo, el mapa no identifica la profundidad de las unidades de suelo que le hacen y
sustentan el tipo suelos. Por lo que, para la determinación del Número de Curva del
espacio geográfico de la cuenca del rio La Paz (Norte), se ha utilizado alternativamente el
mapa de pendientes clasificado y el mapa de cobertura y uso actual de la tierra.

Página 72
Considerando que el mapa de pendiente clasificada por rangos de pendientes tiene
relación con la clasificación de tipo suelo en terrenos ondulados y montañosos como el
área de estudio.

Para la clasificación del tipo de suelo en función de la pendiente del suelo y relacionada
con la clasificación de la Soil Conservance Service (A, B, C y D), se tiene los siguientes
parámetros de relación:
Suelos con pendiente <15%; suelo de tipo A
Suelos con pendiente de 15 a 25%; suelo de tipo B
Suelos con pendiente de 25 a 45%; suelo de tipo C
Suelos con pendiente >45% suelo de tipo D

La pendiente de la superficie topográfica de la cuenca hidrográfica juega un rol importante

en la formación de suelos y la retención de la humedad. Por lo que para fines del presente
estudio inicialmente se ha extractado el mapa de pendientes de la cuenca del rio La Paz
(Norte) en función al perímetro de cuenca y se ha clasificado el mapa de pendientes
relacionados con la clasificación de Soil Conservancey Service (A, B, C y D) indicadas
anteriormente.

Clasificación del mapa de pendiente RASTER con el software ArcGis Spatial Analisis-
Reclass
0 a 15% clase 1 (A)
15 a 25% clase 2 (B)
25 a 45% clase 3 (C)
45 a mayor clase 4 (D)

El mapa de pendiente en formato RASTER clasificado se convierte en la clasificación de

tipo de suelo y se guarda en formato GRID para que genere su propia tabla de atributos.
Entonces, en la tabla de atributos se genera la clasificación de suelo de; 1 (A), 2(B), 3(C)
y 4(D).

El mapa de cubertura y uso actual de la tierra disponible en formato RASTER lleva una
tabla de atributo numérico y descriptivo que codifica y describe la cobertura y uso de
suelo. La combinación de la información de Suelo (Pendiente) y Cobertura de uso

Página 73
(COBUSO) en formato RASTER con el software ArcGis Spatial Analisys – Local -
Combine, se genera el mapa combinado con los atributos de Tipo de Suelo y Cobertura
de Uso en una sola tabla de atributos. En esta tabla se agrega el atributo de Numero de
Curva correspondiente consignada por la Soil Conservance Service.

El número de curva (CN) ponderado está relacionado con la superficie de la unidad y el

número de curva de la matriz, la unidad con mayor superficie incide en el número de
curva ponderado correspondiente. Estos números de curva constituyen un elemento
importante para el cálculo del número de curva ponderado para cada cuenca o sub-
cuenca. El número de curva ponderado es un elemento fundamental para la
transformación de la lluvia en escurrimiento utilizando el software HEC-HMS.

CNponderado = (ΣCNi * Ai) / A

Donde:
CNi Número de curva para determinado complejo
Ai Area parcial que ocupa el complejo
A Area total de la cuenca
CNponderado Número de Curva ponderado para la cuenca

Utilizando la tabla de atributos del mapa combinado (Tipo de suelo y Cobertura de suelo)
y asignando los valores de numero de curva de Soil Conservancy Service a las unidades
de suelo a cada sub-cuenca del modelo correspondiente a los ríos; Choqueyapu,
Choqueyapu1, Orkojahuira, Irpavi, Irpavi1 y Achumani; se ha calculado el Numero de
Curva Ponderado aplicando la formula indicada anteriormente. Ver tabla 4.6 Resumen de
Numero de Curva Ponderado de las Sub-cuencas. El detalle del cálculo de número de
curva se encuentra en Anexo 4.
SUBCUENCA NUMERO DE CURVA
PONDERADA
Choqueyapu 75
Choqueyapu1 76
Orkojahuira 75
Irpavi 78
Irpavi1 75
Achumani; 75

Tabla 4.6 Resumen de Numero de Curva Ponderado de las Sub-cuencas

Página 74
4.6 MODELO DE TRANSFORMACION DE LA LLUVIA EN ESCORRENTIA HEC-HMS

El cálculo de caudal para la verificación de la capacidad hidráulica de la cuenca del rio La

Paz (Norte); considera el cálculo de caudal con el aporte de las cuatro sub-cuencas
(Chuqueyapu, Orkojahuira, Irpavi y Achumani) a través del modelo de transformación de
precipitación en escorrentía HEC-HMS.

El sistema de la red de medición meteorológica de SENAMHI, cuenta con varias

estaciones meteorológicas. La ubicación y las coordenadas de las estaciones
meteorológicas se detallan a continuación. Ver tabla 4.7; Lista de coordenadas de las
estaciones meteorológicas de SENAMHI ubicadas en el área de estudio.

Nro. Estación Departamento Provincia Municipio Latitud Longitud Altitud

1 Achachicala La Paz Murillo La Paz -16.449 -68.148 3960

2 Achumani La Paz Murillo La Paz -16.531 -68.071 3384
3 Alto Achachicala La Paz Murillo La Paz -16.347 -68.085 4383
4 Alto Irpavi - EMI La Paz Murillo La Paz -16.529 -68.080 3265
Alto Obrajes-B.
5 La Paz Murillo La Paz -16.523 -68.101 3480
Magisterio
6 Aranjuez - C SAT La Paz Murillo La Paz -16.559 -68.093 3000
7 Cahua - Zongo La Paz Murillo La Paz -16.041 -67.979 2500
Chacaltaya
8 La Paz Murillo La Paz -16.363 -68.134 4100
Vertedero
9 Chicani La Paz Murillo La Paz -16.482 -68.077 3420
10 Chuquiaguillo La Paz Murillo La Paz -16.450 -68.093 4024
11 Cota Cota - UMSA- La Paz Murillo La Paz -16.538 -68.064 3432
12 El Tejar - La Paz La Paz Murillo La Paz -16.496 -68.156 3820
13 Hampaturi La Paz Murillo La Paz -16.402 -68.015 4336
14 Holguin La Paz Murillo La Paz -16.521 -68.115 3200
15 Huaylipaya La Paz Murillo La Paz -16.043 -68.001 1314
Laycacota - C A
16 La Paz Murillo La Paz -16.504 -68.123 3635
SAT
17 Obrajes La Paz Murillo La Paz -16.531 -68.101 3805
18 Pasankeri La Paz Murillo La Paz -16.523 -68.144 3935
19 San Calixto La Paz Murillo La Paz -16.495 -68.133 3658
20 Villa Copacabana La Paz Murillo La Paz -16.483 -68.114 3680
21 Vino Tinto La Paz Murillo La Paz -16.481 -68.139 3865

Tabla 4.7; Tabla Lista de coordenadas de las estaciones meteorológicas en la cuenca del rio La
Paz (Norte)

Página 75
Las estaciones meteorológicas se encuentran concentradas en el área urbana de la
ciudad de La Paz. La mayoría están ubicadas en el área urbana intensiva de la ciudad de
La Paz, razón por la que las cabeceras de las cuencas tienen información muy limitada.
Ver figura 4.12 Mapa de ubicación de estaciones meteorológicas en la cuenca del rio La
Paz (Norte).

Figura 4.12; Mapa de ubicación de estaciones meteorológicas de la cuenca

Página 76
Asimismo; indicar que con excepción de la estación San Calixto, existen estaciones
meteorologías con registros discontinuos e incompletos, con periodos faltantes o
estaciones que ya no se encuentran en funcionamiento, como por ejemplo las estaciones
ubicadas en las partes altas de las cuencas de los ríos Irpavi, Achumani y Huayñajahuira.

El Anexo 2, se identifica claramente los años con información completa, periodos con
información incompleta y años sin información, en este caso los registros que se
consignan corresponden a precipitaciones máximas en 24 horas (P24) – denominación
oficial establecida por SENAMHI y complementariamente datos de lluvias diarias
registradas, información necesaria para realizar el cálculo de caudales.

De la revisión de la información disponible en el SIS MET de SENAMHI Bolivia, está

disponible la información pluviométrica de las estaciones San Calixto (periodo 1961-1991
de 1995 a 2016), Huaylipa (periodo 1961-2016), Laikacota (1961-2016), Pasankeri (1961-
2016), Villa Copacabana (1980 -2016) y Vino Tinto (1976 – 2003 y 2005 - 2016), siendo la
estación San Calixto la más importante considerando la calidad y longitud de sus
registros, estación que por las razones expuestas es empleada como base en el estudio.
Es importante indicar, que además de la información pluviográfica de San Calixto, se
dispone solamente de los registros de las estaciones meteorológicas de Vino Tinto y
Laykacota con registro de la fuerte tormenta con granizada de febrero de 2002.

4.6.1 ANÁLISIS DE HOMOGENEIDAD DE PRECIPITACIONES

El análisis de homogenidad de las precipitaciones a nivel anual tiene como base los
registros de precipitaciones anuales de las estaciones consignadas en (Tabla 4.13),
empleando el método estadístico denominado Vector Regional. Este método se basa en
los principios del método de las dobles acumuladas, que sustenta la hipótesis de que las
estaciones involucradas en el estudio pertenecen a una misma zona climática, en función
a precipitaciones totales anuales, cualquiera sea la fluctuación climática de la zona; los
totales anuales corresponden al periodo hidrológico anual.

El Vector Regional crea un vector de índices de precipitación anual, correspondientes a

una estación ficticia homogénea o promedia, en base a la información de las estaciones
consideradas en el estudio, esta serie ficticia comprende los efectos de persistencia,

Página 77
 tendencia y pseudo ciclos de la zona climática. La estación ficticia o promedia creada se
convierte en la estación pivote para el método de las dobles acumuladas, contrastando su
valor con cada una de las estaciones para verificar la homogeneidad temporal.

En el Anexo 3, se presentan los valores de Precipitación Total Anual, para el año

hidrológico de las diferentes estaciones que forma parte del estudio considerándose el
periodo 1961 – 2016; se hace notar que el método del Vector Regional requiere como
mínimo tres estaciones y tres años de información, razón por la cual el periodo de análisis
se inicia recién en 1961.

En la Tabla 4.8 se muestran los resultados de una primera corrida del Vector Regional,
donde se muestran, además de los años de información considerada por estación, la
precipitación total anual media observada, la precipitación anual calculada y el coeficiente
de correlación entre el Vector con cada una de las estaciones. El coeficiente de
correlación “r”, es el resultado más importante del proceso, debido a que permite
identificar a estaciones meteorológicas que tienen mala información o que pertenecen a
otra zona climática; un valor aceptable para este parámetro es que sea mayor o igual a
0.7. Las estaciones por debajo de este valor deben ser excluidas.

ESTACION PERIODO DE PRECIPITACION PRECIPITACION COEFICIENTE DE

REGISTROS MEDIA ANUAL MEDIA ANUAL CORRELACION "r"
(AÑOS) X OBSERVADA (mm) Y CALCULADA (mm)

SAN CALIXTO 53 578.5 548.299 1.055081449

ACHACHICALA 16 681.6 596.779 0.875555669

ACHUMANI 17 533.5 595.468 1.116154598

ALTO ACHACHICALA 25 552.6 584.986 1.05860722

ALTO IRPAVI 25 552.6 584.986 1.05860722

ALTO OBRAJES B M 13 562.7 600.710 1.067548502

CHUQUIAGUILLO 7 625.3 608.571 0.973246662

HUAYLIPAYA 56 521.8 544.368 1.043250487

LAYKACOTA 56 521.8 544.368 1.043250487

PASANKERI 56 521.8 544.368 1.043250487

VILLA COPACABANA 37 617.5 569.263 0.921883657

VINO TINTO 40 617.8 565.332 0.915073422

Tabla 4.8; Primera corrida del Vector regional

Página 78
4.6.2 CÁLCULO DE CAUDALES CON MODELO HEC-HMS

Para el cálculo del caudal de diseño de la cuenca y sub-cuencas del rio La Paz (Norte), se
aplicó el modelo de simulación de crecidas HEC-HMS. Este modelo está diseñado para
simular la escorrentía superficial que resulta de una precipitación, mediante la
representación de la cuenca hidrográfica como un sistema de componentes
interconectados. Cada componente modela un aspecto del proceso lluvia – escorrentía
dentro de una sub-cuenca, los componentes del modelo incluyen la escorrentía superficial
de la sub-cuenca, el tránsito de caudales por los canales y los almacenamientos, los
cuales se representan por un conjunto de parámetros que especifican las características
particulares del componente y las relaciones matemáticas que describen sus procesos
físicos. Los resultados finales del proceso de modelación son los hidrogramas de
escorrentía directa para varias sub-cuencas y los hidrogramas de caudal en lugares
predeterminados de la cuenca hidrográfica, como el hidrograma de la escorrentía en el
puente encuentro de los ríos Irpavi y Choqueyapu.

El HEC-HMS para realizar los procesos indicados, requiere información que se introduce
en cuatro bloques: modelo de cuenca, modelo meteorológico, especificaciones de control
y serie de datos expresados en tiempo (ver figura 4.13) Presentación gráfica del modelo.

Figura 4.13; Presentación gráfica del modelo

Página 79
Modelo de cuenca (Basin Model)
En el modelo de cuenca del HEC-HMS, se aplica a las 5 sub-cuencas, con la siguiente
información básica:
a) Superficie de la sub-cuenca expresada en Km2
b) Tasa de pérdidas: Se emplea el método del Soil Conservation Service (SCS) que
sugiere en base a las características geoelogicas y vegetales expresadas en base
a la zonificación de suelo y cobertura vegetal el Numero de Curva ponderado CN’s
y pérdidas iniciales.
c) Transformación de precipitación en escorrentía: Se utiliza el método del SCS,
introduciendo como dato el tiempo de rezago (Log Time) = 0.6*Tc.
d) Tránsito de crecidas: En tramos canalizados o embovedados de sección regular se
empleó el método Muskingum Estándar que requiere geometría de la estructura y
rugosidad, en tramos de canales naturales se aplica el método de Muskingum que
permite introducir una sección variable y rugosidades características.

Modelo meteorológico (Meteorologic model)

En el modelo meteorológico del HEC-HMS, se emplea hietograma de la lluvia del
proyecto, considerándose lluvia calculada con curva de Intensidad, Duración y Frecuencia
(IDF) derivada con los datos de la Estación San Calixto para 50 años de retorno.

A efectos de sustentar lo indicado y para establecer el comportamiento de la intensidad,

duración y frecuencia de las tormentas en la cuenca del río La Paz, se ha seleccionado la
precipitación máxima en 24 horas mayores a los 30 mm, medidos en el pluviómetro de
San Calixto y registrados por el pluviógrafo de la misma estación. Para la selección de la
precipitación instantánea más fuerte se ha considerado que la precipitación total
registrada en un día (24 horas) puede estar compuesta de varias tormentas. El criterio
empleado para la separación de las tormentas es que si no existe precipitación por más
de una hora se divide la precipitación de las 24 horas de precipitación consignada como
dato meteorológico, esto debido a las características de las tormentas de la región que
son de origen conectivo, en ningún caso las tormentas superan las 4 horas de duración y
son poco frecuente, con característica de gran intensidad y de corta duración.

Las precipitaciones empleadas son las que corresponden a recurrencias de 10, 25, 50 y

Página 80
100 años, determinadas en base a las ecuaciones de Intensidad – Duración – Frecuencia
calculadas para cada estación, introduciéndose las mismas con un intervalo de 15
minutos hasta completar las horas del tiempo de concentración. Sin aplicación de factores
de reducción de precipitación por área de las sub-cuencas hidrográficas por no
considerarse necesarios, ya que principalmente en la parte media y baja de la cuenca
hidrográfica, las áreas de aporte son reducidas, por lo que cada estación utilizada
representa suficientemente bien a dicha área. En la parte alta de la cuenca, la reducción
resultante, considerando la superficie de las sub-cuencas hidrográficas, es en algunos
casos del 1%, por tanto, se opta por no aplicar estos factores por considerarlos “no
significativos” para el cálculo. Entonces en el presente caso se ha aplicado la opción de
hietograma derivado de lluvias con 50 años de retorno (datos del pluviografo san Calixto),
introducido por el usuario a través de hietogramas acumulados.

Control de especificaciones (Control Specification)

En el control de especificaciones se ha considerado un tiempo de 8.25 horas, periodo en
el cual se asume que el efecto de la crecida del caudal a la salida de la cuenca
hidrográfica del rio La Paz en la zona de La Florida ya es poco significativo.

Se ingresaron los datos de la lluvia del proyecto generados con el IDF estación
pluviografica de San Calixto en el modelo con un intervalo de 15 minutos hasta completar
la duración del tiempo de concentración de la lluvia. No se aplicaron factores de reducción
de precipitación por área a las sub-cuencas hidrográficas por no considerarse necesarios,
ya que principalmente en la parte media y baja de la cuenca hidrográfica, las áreas de
aporte son reducidas, por lo que la estación utilizada representa apropiadamente bien a
dicho área. En la parte alta de la cuenca hidrográfica, la reducción resultante,
considerando la superficie de las sub-cuencas hidrográficas, es en algunos casos del 1%,
por tanto, se optó por no aplicar estos factores por considerarlos “no significativos” para el
cálculo. Se aplicó la opción de hietograma introducido por el usuario a través de
hietogramas acumulados.

Datos en Serie de Tiempos (Time-Series Data)

En los datos; Serie de tiempos se ha ingresado al modelo los datos de hietograma
introducido manualmente por el usuario a través de hietogramas acumulados, Ver figura
4.14; Datos de la serie de tiempos.

Página 81
Figura 4.14) Datos de la serie de tiempos

4.6.3 PROCESAMIENTO DE LOS DATOS DE LA CUENCA HIDROGRAFICA

Con los datos introducidos en HEC-HMS: Modelo de Cuenca, Modelo Meteorológico,

modelo de Control de especificaciones y Datos en serie de tiempo, se ha procesado los
datos y se ha obtenido los resultados de los caudales de diseño de crecida obtenidos
como resultado de numerosas calibraciones, tanto en lo referente a los números de curva,
tiempos de concentración y rugosidades (Ver figura 4.15); Sumario global de los
resultados del procesamiento de datos de precipitación.

Figura 4.15); Sumario global de los resultados del procesamiento de los datos de precipitación.

Página 82
En el siguiente gráfico y tabla sumario, se muestran los resultados del procesamiento de
datos de la Sub-Cuenca del rio Chuqueyapu. El gráfico representa la precipitación
excedente y pérdida de la precipitación para el cálculo de la escorrentía. Esto debido a
las características de absorción del suelo y la cobertura vegetal de la sub-cuenca. Es
posible apreciar en el hidrograma aspectos importantes del comportamiento de la
precipitación. Ver figura 4.16; Gráfico del Hidrograma y Tabla Sumario resultado del
procesamiento de los datos de precipitación.

Figura 4.16; Gráfico del hidrograma y Tabla Sumario resultado del procesamiento de datos.

Página 83
Rio Choqueyapu

El gráfico y la tabla sumaria muestran el flujo hidrográfico que se presenta en el curso de

agua del rio Choqueyapu durante la precipitación, fundamentalmente el pico de la
descarga y el flujo hidrográfico. Ver figura 4.17; Gráfico y Tabla Sumario de la descarga y
flujo hidrográfico rio Choqueyapu.

Figura 4.17; Gráfico y tabla sumaria de la descarga y flujo hidrográfico del rio Choqueyapu.

Unión entre los ríos Choqueyapu y Orkojahuira

El gráfico y la tabla sumaria de la unión entre los ríos de las subcuencas de Choqueyapu
y Orkojahuira muestra el caudal hidrográfico que se presenta en la unión, la máxima
descarga y el volumen en un tiempo determinado durante la precipitación. Ver figura 4.18;
Gráfico y Tabla Sumario de la máxima descarga y el volumen en un tiempo determinado.

Página 84
Figura 4.18; Gráfico y tabla sumaria de la máxima descarga y el volumen en un tiempo
determinado de la unión de los ríos Choqueyapu y Orkojahuia.

Unión de los ríos Choqueyapú e Irpavi

El gráfico y la Tabla Sumario de la unión entre los ríos Irpavi y Choqueyapu muestra el
caudal hidrográfico que se presenta en la unión, la máxima descarga de 56.7 m3/s y el
volumen en un tiempo determinado durante la precipitación (ver figura 4.19) Gráfico y
Tabla Sumario de la máxima descarga y el volumen en un tiempo determinado.

Figura 4.19; Gráfico y tabla sumaria de la máxima descarga y el volumen en un tiempo

determinado de la unión de los ríos Choqueyapu e Irpavi.

Página 85
4.7 LEVANTAMIENTO TOPOGRÁFICO EN DETALLE

El tramo del curso del rio La Paz (Norte) desde el puente encuentro entre los ríos de
Chuqueyapu e Irpavi hasta el puente cruce a Mallasa (Zona Aranjuez), tiene
características particulares a considerar para la realización del levantamiento topográfico
correspondiente:
- A partir del puente de encuentro se expande el ancho con deposición de una gran
cantidad de material de arrastre compuesta por bloque de piedras de tamaño variado,
más abajo depósitos de áridos de diferentes tamaños, clastros y arena gruesa para
construcción que es explotado por los concesionarios.
- El rio anualmente lleva consigo aguas pluviales y permanentemente aguas servidas
que se generan por la intensa actividad humana en el área urbana de la ciudad de La
Paz. Es altamente contaminada y lleva elementos nocivos para la salud humana.
- La topografía del indicado tramo es bastante accidentada.

Con estas consideraciones se ha seleccionado el método para la realización del

levantamiento topográfico por el método indirecto mediante la ejecución de las siguientes
actividades:
- Toma de imágenes fotogramétricas con UAV (DRON)
- Determinación de puntos de control fotoidentificables
- Procesamiento de datos
o Aerotriangulación
o Generación de MDT
o Generación de Ortofotos

4.7.1 TOMA DE IMÁGENES FOTOGRAMÉTRICAS CON UAV o DRON.

La tecnología hace disponible para la realización del levantamiento topográfico por

método indirecto el sistema UAV (DRON), basado en la recolección de datos desde el
espacio aéreo sin entrar en contacto con la superficie topográfica por donde drenan las
aguas servidas. Una nave no tripulada UAV equipada y controlada con sistema de
posicionamiento global (GPS) e incorporada una cámara fotográfica digital a color permite
la toma de imágenes aéreas con las características fotogramétricas establecidas en las
especificaciones técnicas.

Página 86
Las imágenes aéreas fueron tomadas con las siguientes características:
- eBee UAV (DRON) de ala fija.
- Camera Model Name(s) DSC-WX220_4.4_4896x3672 Color (RGB).
- Average Ground Sampling Distance (GSD) 5.07 cm / 1.99 in.
- Area Covered 0.7808 km2 / 78.0836 ha / 0.3016 sq. mi. / 193.049 acres.
- Time for Initial Processing (without report) 44m:31s.

Plan de vuelo.
El sistema SenseFly cuenta con el software eMotion 2 para la planificación del vuelo con
las UAVs o DRONs, herramienta tecnológica en la que se consigna todos los parámetros
y características de la fuente de información a ser recolectados (datos e imágenes) para
satisfacer su aplicación en la producción del modelo digital de terreno y el mosaico de
ortofotos. En la siguiente tabla 4.9, se muestra el plan de vuelo para la realización del
vuelo y la toma de imágenes a color.

No. características Descripción

Numero de Vuelo 1
Tiempo de vuelo total continuo 19:09 minutos
Distancia total volada 9.5 Km.
Área Total 19.51 hectáreas
Numero de líneas de vuelo 6
Sobreposición longitudinal 80%
Sobreposición Transversal 60%
Espacio entre líneas de vuelo 88.1 m.
Distancia entre fotos 32.5
Cobertura por imagen 220.3 m x 162.5 m.
# Waypoints 10
Dias de Trabajo 1 día
Horas de procesamiento automático de información 3 días
Waypoint de referencia 68.0912091W; long
16.5522260S; lat
Altura promedio de vuelo 176.6 m.
Resolución espacial 5 cm/px

Tabla 4.9; Plan de vuelo para la realización del vuelo y la toma de imágenes a color

Página 87
4.7.2 DETERMINACIÓN DE PUNTOS DE CONTROL

El establecimiento de las coordenadas geodésicas y la altura de los puntos de control

fotoidentificables, se ha realizado mediante el posicionamiento GPS Diferencial (Doble
frecuencia) en los siguientes términos:
- La determinación de las coordenadas planimetrías (Este, Norte) y altura de los 4
puntos de control, se ha realizado en base a las coordenadas y datos de las
estaciones permanentes MARGEN-SIRGAS, administradas por el Instituto Geográfico
Militar consignados en el sistema de control como BLPZ y EMI respectivamente.

- Tiempo de observación de datos GPS; de 2 horas aproximadamente con registro de

datos GPS cada 15 segundos y mascara de elevación 15 grados sobre el horizonte.

- La ubicación de los cuatro puntos de control fotoidentificable fue seleccionado

utilizando las imágenes aéreas agrupadas en mosaico preliminar de imágenes
tomadas con el DRON, estos puntos fotoidentificables que se muestra a continuación:
Imagen 4.8 PCT1, Imagen 4.9 PCT2, imagen 4.10 PCT3 e imagen 4.11 PCT4.

Imagen 4.8 PCT1 Imagen 4.9 PCT2

Página 88
Imagen 4.10 PCT3 Imagen 4.11 PCT4

El procesamiento de datos GPS y el ajuste de coordenadas fuer realizado con el software

Leica Geoffice. Las coordenadas y alturas ajustadas están en la siguiente Tabla 4.10:
Punto Coord-Este Coord-Norte Altura Elip Altura geometrica
INGA 588751.180 8172270.659 4090.652 4040.941
EMI 597989.396 8171967.669 3411.337
PC1 596670.972 8170719.438 3277.357 3227.632
PC2 596808.988 8170661.599 3277.583 3227.846
PC3 596588.924 8168900.641 3236.958 3187.216
PC4 596660.496 8168980.477 3238.955 3189.226

Tabla 4.10; Coordenadas de Puntos de Control Fotoidentificables

La exactitud de las coordenadas geodésicas de los puntos de Control respecto a las

coordenadas de las estaciones permanentes administradas por el IGM es de RMSE= ±10
cm. La exactitud del desnivel de altura obtenida es RMSE=± 25 cm.

4.8 GENARACION DEL MODELO DIGITAL DE TERRENO y ORTOFOTOS

El sistema eBee SenseFly de ala fija controlada por sistema de posicionamiento global
(GPS). Cuenta con software eMotion 2 para planificación del vuelo y toma de datos, y el
Pix4D Mapper para el post procesamiento profesional y generación de resultados.

El sistema utiliza el software Pix4D Mapper para la georeferenciación y ajuste del bloque
de fotografías aéreas y calcular la ubicación exacta del centro de la fotografía digital y los

Página 89
parámetros de orientación externa (Omega, Phi y Kapa) del bloque de imágenes aéreas a
efectos de restablecer espacialmente el momento de la toma de fotografías aéreas para la
generación de los productos MDT/MDS y Ortodotos.

4.8.1 GENERACIÓN DE MDT/MDS

El software de postprocesamiento Pix4D Mapper a partir de las fotografías aéreas a color

permite generar la nube de puntos con cobertura del aérea de interés en formato .LAS.
Estos datos en formato .LAS son procesados con el software aplicativo indicado para
eliminar todo tipo de puntos al margen de la clasificación de la nube puntos .LAS en las
siguientes clases:
- Punto .LAS a nivel de suelo.
- Punto .LAS correspondientes a vegetación alta
- Punto .LAS correspondientes a vegetación media
- Punto .LAS correspondientes a las construcciones

La nube de puntos en formato .LAS, editadas y clasificadas a nivel de suelo, vegetación

alta, vegetación media y construcciones; juntas constituyen el modelo digital de
Superficie en formato vectorial. Esta información conformada en Triangular Irregular
Network (TIN) y la rasterización correspondiente del modelo digital se obtiene el Modelo
Digital de Superficie (MDS) en formato raster con resolución espacial de 0.50 m. y una
exactitud estimada de RMSEz ± 0.50 m.

Solamente los puntos .LAS a nivel de suelo constituyen la fuente de información para la
generación del modelo digital de terreno en formato vectorial. Esta información a través
de la conformación del TIN y la rasterización el modelo digital de terreno se obtiene el
MDT en formato de presentación raster. Ver figura 4.20; Modelo Digital del Terreno del
área de proyecto en formato raster.

Página 90
Figura 4.20; Modelo Digital Terreno (formato raster).

4.8.2 GENERACIÓN DE ORTOFOTOS

Finalmente, las ortofotos derivado del vuelo fotogramétrico con DRON son generadas
mediante el proceso de Rectificación Diferencial, que consiste en establecer al Modelo
Digital de Superficie como base para segmentación milimétrica de las imágenes
fotográficas para la eliminación del desplazamiento por relieve. El mosaicado de las
imágenes segmentadas constituyen la ortofoto con proyección ortogonal al igual que la

Página 91
cartografía vectorial producida mediante restitución fotogramétrica. Ver figura 4.21;
Ortofoto del área de proyecto.

Figura 4.21; Ortofoto del área de proyecto.

4.9 MODELO DE DINAMICO FLUVIAL IBER

Para simular la respuesta del vertedero natural a la acción hidrográfica de la cuenca del
rio La Paz (Norte) y de su red de drenaje pluvial como consecuencia de un evento con
fuertes precipitaciones, se utiliza la herramienta de modelación dinámica y simulación
matemática en hidráulica. En la parte hidráulica para determinar el comportamiento del

Página 92
sistema de drenaje se ha utilizado el software libre IBER desarrollado por el Grupo de
Ingeniería del Agua y del Medio Ambiente, GEAMA (Universidad de A Coruña, UDC) y el Instituto
FLUMEN (Universitat Politècnica de Catalunya, UPC, y Centro Internacional de Métodos
Numéricos en Ingeniería, CIMNE), es un modelo muy reconocido y de amplia aplicación en
todo el mundo.

Para determinar los parámetros hidráulicos es necesario realizar el relevamiento del

vertedero natural y superficial, determinándose los parámetros hidráulicos como ser
coeficientes de rugosidad expansión y contracción. Luego introducir la geometría del
canal natural en forma de malla del modelo digital del terreno. Una vez establecidas las
características físicas (geometría y coeficientes), se realiza el armado de la red para la
simulación del modelo.

4.9.1 CRITERIOS HIDRÁULICOS

Para fines de este estudio se ha llevado adelante la identificación de los parámetros

morfométrico del rio y los criterios del componente hidráulico del proyecto. A continuación,
se desarrollan los criterios hidráulicos desde el punto de vista del modelo Bidimensional
IBER, programa elegido para la modelación del tramo de la red seleccionada de la cuenca
de la ciudad de La Paz (puente encuentro de los ríos de Irpavi y Choqueyapu hasta la
confluencia con el rio Huayñajahuira (Zona la Florida).

Los criterios hidráulicos considerados son: el periodo de retorno, el régimen de flujo, el

tipo de flujo y los coeficientes de rugosidad. Los mismos que se desarrollan a
continuación.

PERIODO DE RETORNO

Para la evaluación de la capacidad hidráulica de vertederos natural/artificial y/o diseño de

las estructuras de conservación de suelos e hidráulicas, es necesario el estudio de las
precipitaciones máximas. El período de retorno a utilizarse será mayor cuanto mayor sea
la importancia y la repercusión social, ecológica y económica de la obra civil. Así la
necesidad de disponer de amplios períodos de retorno contrasta con la disponibilidad de
series de datos climatológicos, por lo que se debe recurrir a estimaciones estadísticas.

Página 93
El periodo de retorno seleccionado para la evaluación hidráulica del rio La Paz; tramo
puente encuentro de los ríos Choqueyapu e Irpavi hasta la confluencia con el rio
Huayñajahuira para una avenida máxima es de 50 años de retorno

RÉGIMEN DE FLUJO

Para fines de este estudio se ha considerado una simulación bajo régimen permanente
aplicando el caudal de diseño que resulta del modelamiento hidrológico de superficie con
el software HEC-HMS que deriva de los datos de la estación San Calixto; precipitación
pluvial máxima en 24 horas con un periodo de retorno de 50 años, datos que permiten la
reproducción de la condición más desfavorable. Esta simplificación posibilita representar
adecuadamente la onda de avenida con el caudal pico. Se asume por caudal pico como el
máximo caudal producido por la onda de avenida.

La simulación hidráulica del área de estudio bajo régimen permanente no considera

posibles obstrucciones y reducciones de sección ocasionadas por sedimentos y/o
escombros. La presencia de obstrucciones, sedimentos o escombros, no es algo que
debiera preocupar, siempre y cuando éstos tengan dimensiones menores a la sección de
conducto. Bajo tales circunstancias, el obstáculo tendrá un efecto temporal sobre el
comportamiento hidráulico de la conducción, siendo removido por la fuerza del agua y no
así cuando pierde la fuerza.

COEFICIENTES DE RUGOSIDAD

A lo largo de su cauce, el flujo sufre pérdidas de energía debidas a la resistencia al flujo

que presentan las paredes y el lecho del cauce. Una manera de describir la rugosidad del
perimetro mojado es mediante el número de Manning. El número de Manning n es un
coeficiente que describe la rugosidad de contorno. Se utiliza para tomar en cuenta la
pérdida de energía debido a la fricción entre el lecho y el agua. En hidráulica fluvial
(hidráulica de contorno movible), el coeficiente de Manning incluye los efectos de otras
pérdidas, como es el caso de la rugosidad de partícula, la rugosidad de forma del lecho,
las irregularidades de los bancos, la vegetación, las pérdidas por serpenteo, y las
pérdidas debido al efecto de confluencia. (USACE, 1993). Para los fines de estudio se
adoptan los siguientes valores de ¨Número de Manning¨ (Ven Te Chow):

Página 94
 • Para embovedados con solera revestida: n = 0.020
• Para canalizaciones abiertas con lecho móvil: n = 0.030 (piso)
n = 0.020 (pared)
• Para canalizaciones abiertas con solera revestida: n = 0.020 (piso)
n = 0.020 (pared)
• Para cursos naturales con o sin obras de control: n = 0.035

4.9.2 MODELAMIENTO DINAMICO

El modelamiento dinámico con el software IBER comprende; el establecimiento de las

condiciones del modelo en base a los datos y resultados requeridos, procesamiento de los
datos y generación de resultados para el análisis correspondiente.

DATOS REQUERIDOS

Para la conformación del modelo se requiere fundamentalmente de: Caudal pico resultado
del modelo hidrológico HEC-HMS, Modelo Digital de Terreno (formato txt) y los valores
estimados del Coeficiente de Rugosidad.

Caudal pico resultado del HEC-HMS

En el modelamiento hidrológico de superficie de la cuenca del rio La Paz (Norte) con los
datos de precipitación máxima en 24 horas registrada en la estación meteorológica San
Calixto y el uso del software HEC-HMS, se ha transformado precipitación en escorrentía,
derivándose de este proceso el caudal pico.

Modelo Digital de Terreno (formato txt)

El MDT con cobertura del tramo de estudio; puente encuentro de los ríos Choqueyapu e
Irpavi hasta la confluencia con el rio Huayñajahuira. El MDT con cobertura del curso del
rio y la franja de seguridad a partir del margen del rio.

Coeficiente de Rugosidad
Los valores de rugosidad (coeficiente de Maning) estimados de acuerdo a las
características del lecho del rio y la franja de seguridad. Para el presente caso se tiene
estimado el coeficiente de rugosidad del lecho del rio (0.030) y la franja de seguridad
ocupada por infraestructura urbana (0.025).

Página 95
PROCEDIMIENTO

Creación de un proyecto: En el software IBER crear el proyecto con el nombre apropiado

para el desarrollo del proyecto de modelamiento dinámico (ej. Corrida_q57_1_4).

Poligonizar el área de estudio: En el software ArcGis (formato shape) crear un elemento

polígono del tramo del rio objeto de estudio y polígonos de la franja de seguridad a partir
de ambos margines del rio. Estos polígonos deben ser topológicamente consistentes. Es
decir, no debe haber vacíos ni sobreposiciones entre los polígonos, previo a la
importación de la información gráfica en el software IBER. El modelo digital del terreno
para uso en el modelamiento dinámico con IBER debe estar en formato .txt
(topografía.txt). En el software IBER desde Geometría – edición- Colapsar; hacer el
colapse de los polígonos en línea y superficie respectivamente. Ver figura 4.22 Polígonos
para espaciamiento diferenciado de la malla.

Figura 4.22 Polígonos para espaciamiento diferenciado de la malla.

Definir la malla: En el software IBER definir malla-no estructurada-asignar. Polígono del

cauce del rio con espaciamiento de 1 m y los polígonos de franja de seguridad con
espaciamiento de 4 m. Luego con Generar Malla se muestra la malla definida. Ver figura
4.23 Malla con espaciamiento diferenciado del cauce y las franjas de seguridad del rio.

Página 96
Figura 4.23 Malla con espaciamiento diferenciado del cauce y las franjas de seguridad del rio.

Simulación: en el software IBER Datos-datos del problema-parámetros de tiempo,

establecer parámetros de tiempo; tiempo máximo de simulación 1500 seg e intervalo de
resultados cada 30 seg. Luego en Datos-hidrodinamica-condiciones de contorno
establecer las condiciones de entrada 2D (critico/subcrítico y el caudal pico de HEC-HMS)
y salida 2D; asignar la entrada y salida en el polígono del cauce del tramo de rio objeto de
estudio. Finalmente, en Datos-rugosidad-uso de suelo, asignar al cauce natural del rio el
coeficiente de Manning (rio 0.030) y a la franja de seguridad del rio donde se encuentra la
infraestructura urbana (infraestructura 0.025) y asignar a los respectivos polígonos.

Subir la topografía con el modelo digital de terreno a la malla generada: Herramientas

Iber-Malla-editar-asignar elevación desde archivo; subir la información topográfia.txt
(MDT) para que la malla adquiera las alturas del modelo digital del terreno con las que el
software IBER simulará la capacidad hidráulica del tramo de rio objeto de estudio.
Asimismo, visualizar de fondo la ortofoto del área de estudio; Vista-Imagen de fondo-
Mostrar. Ver figura 4.24 Malla del área de estudio fondo ortofoto.

Página 97
Figura 4.24 Malla del área de estudio fondo ortofoto

PROCESAMIENTO DE DATOS
Calcular con el software IBER los factores hidráulicos del cauce del rio; ejecutando en el
menú principal Calcular, corre el software IBER con los cálculos correspondientes para la
generación de los resultados como la cota de agua y velocidad del agua. Con el
Postproceso se generan los resultados verificables en formato ASCCI y exportables a
formato Raster para visualización en sistemas de información geográfica y la
conformación de la base gráfica para la cartografía de riesgos.

RESULTADOS
Resultados del procesamiento corrida q57m3/s_1_4 con caudal pico 57m3/s con
espaciamiento de la malla en el área del cauce del rio 1 m, la franja de seguridad de 4 m y
el caudal resultado de la transformación de la lluvia del proyecto en escorrentía por HEC-
HMS. Los principales resultados son: Cota de Agua y Velocidad de agua. Ver figura 4.25
Cota de Agua (m) y figura 4.26 Velocidad de agua (m/s).

Página 98
Figura 4.25 Cota de Agua (m).

Figura 4.26 Velocidad de agua (m/s)

Página 99
Asimismo, para el análisis se tiene resultados del procesamiento con caudal 75m3/s;
corrida q75_1_4, con caudal 100m3/s; corrida q100_1_4, con caudal 15m3/s; corrida
q15_1_4, con caudal 40m3/s; corrida q40_1_4, con caudal 60m3/s; corrida q60_1_4 y con
caudal 80m3/s; corrida q80_1_4; con espaciamiento de la malla en el área del cauce del
rio 1 m, la franja de seguridad del rio con 4 m. Los principales resultados también son:
Cota de Agua (m) y Velocidad de agua (m/s).

4.9.3 ANALISIS DE RESULTADOS

Para el análisis de los resultados del procesamiento de datos hidrológicos con el software
IBER, partimos de la identificación de las condiciones técnicas del modelo digital del
terreno (MDT) elaborado del área de proyecto, cuya característica fundamental es que
tiene una resolución espacial tamaño de pixel de 0.50 m. y tiene una exactitud estimada
en altura RMSEz ± 0.40 cm. Entonces el espaciamiento de la malla sugerido es mayor a
la resolución espacial. Por lo que se ha definido el espaciamiento de la malla diferenciada;
en el área del cauce del río el espaciamiento de la malla es de 1 m y en el área de la
franja de seguridad del rio (infraestructura urbana) el espaciamiento de la malla es de 4 m.

Los ríos de Choqueyapo e Irpavi antes del puente encuentro son canalizados a través de
vertederos construidos para el control del cauce y velocidad del agua correspondiente. El
último aliviador regulador de la velocidad de agua del vertedero, se encuentra
aproximadamente a 42 m. del margen sur del puente encuentro. A partir de este último
aliviador regulador de velocidad de agua, el cauce es natural y casi sin control ya que se
produce una caída de pendiente de -14% hasta unos 25 m aguas abajo. Lugar donde se
deposita gran cantidad de material de arrastre formando un abanico aluvial,
fundamentalmente formado por el material de arrastre del rio Irpavi.

En la Ortofoto del lugar del cambio de pendiente de -14% a -2% a 25 m aguas abajo del
último aliviador regulador de la velocidad de agua del vertedero. Ver figura 4.27 Remoción
de material de arrastre, se puede observar que el GAM La Paz, ha realizado los trabajos
de remoción del material de arrastre a ambos márgenes del rio, evitando la sedimentación
y consolidación de abanico aluvial. Asimismo, se puede evidenciar en la imagen de las
ortofotos que hay acceso al valle del rio por vehículo, medio por el que se ha evacuado

Página 100
gran cantidad de áridos generalmente requeridos como material de construcción de
infraestructura urbana pública y privada de la ciudad de La Paz.

Figura 4.27 Muestra de remoción de material de arrastre.

A partir de la ortoimagen elaborado para el presente estudio, se ha realizado la

digitalización del cauce del rio La Paz, considerando la huella de la crecida del rio por
ambos márgenes del rio a efectos analizar el comportamiento de la capacidad hidráulica
del tramo del rio para una avenida con 50 años de retorno (caudal pico de 57 m3/s
modelado con HEC-HMS) y el comportamiento para los caudales asumidos de 75 m3/s,
100 m3/s, 15 m3/s, 40 m3/s y 60 m3/s.

Los resultados de la corrida del modelo dinámico IBER con caudal 57 m3/s; corrida
q57_1_4; en lo que corresponde a cota de agua (m) y velocidad (m/s), se puede
evidenciar lo siguiente:

Cota de agua: la imagen de la cota de agua muestra que la cota del último regulador del
vertedero del encuentro de los ríos de Choqueyapu e Irpavi (Naciente del rio La Paz) tiene

Página 101
como cota 3224 m. y la unión con el rio Huayñajahuira tiene la cota de 3219 m. Ver figura
4.28 Imágenes de la cota de agua para los caudales 57m3/s, 75m3/ y 100m3/s. En general
las cotas de agua del tramo guardan una relación apropiada con el modelo digital del
terreno, es decir que conforme va bajando el rio en su cauce natural la cota también va
bajando y no se evidencia que pudiera haber riesgo de desborde del rio. Lo que significa
que la capacidad hidráulica del tramo del rio puede soportar una avenida con 50 años de
retorno o más.

3 3 3
Figura 4.28 Imagines de la cota de agua para los caudales 57m /s, 75m / y 100m /s

Velocidad del agua: la imagen de la velocidad de agua muestra que el agua por el
vertedero viene con una velocidad de 1.76 m/s aproximadamente y va por el curso natural
aguas abajo del último aliviador regulador de velocidad del vertedero con una velocidad
de 4.41 m/s. con excepción del salto hidráulico que se presenta aguas abajo del último
aliviador regulador de velocidad del vertedero a unos 25 m. En el salto hidráulica o cambio
de velocidad del agua se tiene una velocidad de hasta 7.93 m/s. Ver figura 4.29 Imágenes
de velocidad de agua para los caudales q57, q75 y q100. Al no haber en el curso del rio
La Paz más fuertes cambios de pendiente, el rio para el caudal de 57 m3/s tiende a
mantenerse a la velocidad constante de 4.41 m/s aguas abajo. Lo cual evidencia que el
rio por el caudal considerado llevará gran cantidad de material de arrastre.

Asimismo, indicar que, en la imagen de velocidad de agua en los márgenes del cuerpo de
agua, la velocidad es baja de hasta 0.88 m/s. Esto fenómeno se debe a la incidencia de la
rugosidad del suelo y también se reproduce este fenómeno eventualmente en los lugares

Página 102
de acceso de ríos contribuyentes menores al rio principal.

3 3 3
Figura 4.29 Imágenes de velocidad de agua para los caudales 57m /s, 75m /s y 100m /s

Los resultados de las corridas del modelo dinámico IBER con caudales de referencia de
75 m3/s y 100 m3/s; corrida q57_1_4 y corrida q100_1_4; en lo que corresponde a cota
de agua (m) y velocidad (m/s), en imágenes tienen características similares a la corrida
con el caudal de 57 m3/s. con un ligero incremento de velocidad en el salto hidráulico por
la caída de pendiente aguas abajo del último aliviador regulador de velocidad del
vertedero de hasta 8.28 m/s y 8.67 m/s respectivamente.

Para el análisis de los resultados de cota de agua y velocidad de agua, se ha procesado

los datos con el software IBER para caudales de 15m3/s (corrida q15_1_4), 40m3/s
(corrida q40_1_4), 60m3/s (corrida q60_1_4) y 80m3/s (corrida q80_1_4). La cota de
referencia de la máxima crecida del rio (2017); identificada a través de la huella de ambos
márgenes del rio La Paz (Norte) en ortoimagen y el MDT es de 3220.7 msnm. Este
elemento gráfico es la referencia básica para comparar el perfil de cota de agua para los
caudales 15m3/s, 40m3/s, 60m3/s y 80m3/s.

Perfil de Cota de Agua (IBER):

Para generar la sección transversal de cota de agua para el caudal de 15m3/s, ejecutar el
proyecto corrida q15_1_4 en el IBER. Crear secciones transversales-dibujar lamina de

Página 103
agua máxima-opciones de conjunto de gráficas (mínimo 3219 y máxima 3224), se ha
obtenido la sección transversal a 40 m aguas abajo del salto hidráulico. Ver figura 4.30
Perfil de Cota de Agua para el caudal de 15m3/s. En el gráfico el contorno Corte 4(7) Paso
0; es el lecho del rio reflejado por la topografía del área de estudio y el gráfico del
contorno corte 4(8) Paso 1500; es la cota de agua para un caudal de 15m3/s.

3
Figura 4.30 Perfil de Cota de Agua para el caudal de 15m /s

De esta presentación gráfica se puede evidenciar que con el caudal 15m3/s la cota de
agua no sobrepasa la cota de referencia de la máxima crecida del rio (año 2017).

Asimismo, para generar la sección transversal de cota de agua para el caudal de 40m3/s,
ejecutar el proyecto corrida q40_1_4 en IBER. Crear secciones transversales-dibujar
lamina de agua máxima-opciones de conjunto de gráficas (mínimo 3219 y máxima 3224),
se tiene la sección transversal a 40 m aguas abajo del salto hidráulico. Ver figura 4.31
Perfil de Cota de Agua para el caudal de 40m3/s. En el gráfico el contorno Corte 1(1) Paso
0; es el lecho del rio reflejado por la topografía del área de estudio y el gráfico del
contorno corte 1(4) Paso 1500; es la cota de agua para un caudal de 40m3/s.

Página 104
 3
Figura 4.31 Perfil de Cota de Agua para el caudal de 40m /s

De esta presentación gráfica se puede evidenciar que con el caudal 40m3/s la cota de
agua ya sobrepasa parcialmente en el sector Oeste de la sección transversal la cota de
referencia de la máxima crecida del rio (2017).

De la misma manera; para generar la sección transversal de cota de agua para el caudal
de 60m3/s, ejecutar el proyecto corrida q60_1_4 en IBER. Crear secciones transversales-
dibujar lámina de agua máxima-opciones de conjunto de gráficas (mínimo 3219 y máxima
3224), se tiene la sección transversal a 40 m aguas abajo del salto hidráulico. Ver figura
4.32 Perfil de Cota de Agua para el caudal de 60m3/s. En el gráfico el contorno Corte 1(1)
Paso 0; es el lecho del rio reflejado por la topografía del área de estudio y el gráfico del
contorno corte 1(4) Paso 1500; es la cota de agua para un caudal de 60m3/s.

3
Figura 4.32 Perfil de Cota de Agua para el caudal de 60m /s

Página 105
De esta presentación gráfica se puede evidenciar que con el caudal 60m3/s la cota de
agua ya sobrepasa casi totalmente la cota de referencia de la máxima crecida del rio
(2017).

Finalmente, para generar la sección transversal de cota de agua para el caudal de 80m3/s,
ejecutar el proyecto corrida q80_1_4 en IBER. Crear secciones transversales-dibujar
lamina de agua máxima-opciones de conjunto de gráficas (mínimo 3219 y máxima 3224),
se tiene la sección transversal a 40 m aguas abajo del salto hidráulico. Ver figura 4.33
Perfil de Cota de Agua para el caudal 80m3/s. En el gráfico el contorno Corte 2(3) Paso 0;
es el lecho del rio reflejado por la topografía del área de estudio y el gráfico del contorno
corte 2(4) Paso 1500; es la cota de agua para un caudal.

3
Figura 4.33 Perfil de Cota de Agua para el caudal de 80m /s

De esta presentación grafica se puede evidenciar que con el caudal 80m3/s la cota de
agua ya sobrepasa totalmente la cota de referencia de la máxima crecida del rio (gestión
2017).

Perfil de Velocidad de Agua (ArcGis-Exel)

Los resultados de la velocidad de agua del BER; correspondiente a los caudales de:
15m3/s, 40m3/s, 60m3/s y 80m3/s, para su visualización en ArcGis son convertidos a
formato raster .tif. Sobre esta información espacial de las velocidades de agua para el
caudal correspondiente y el MDT apoyado con las ortoimágenes, se tiene seleccionado la
sección transversal a 40 m aguas abajo del salto hidráulico.

Página 106
Con ArcGis 3D analisys Tool-Funcional surface-Interpolate shape, se realiza la
interpolación entre el vector de la sección transversal y la información espacial de
velocidad de agua resultante del IBER para los caudales de 15m3/s, 40m3/s, 60m3/s y
80m3/s. Los datos de la interpolación son exportados a formato EXEL para la
organización y representación de la información en EXEL. Ver figura 4.34 Velocidades de
agua para los caudales de 15m3/s, 40m3/s, 60m3/s y 80m3/s.

Figura 4.34 Velocidades de agua para los caudales de 15m3/s, 40m3/s, 60m3/s y 80m3/s.

De la representación grafica de la velocidad de agua para los caudales 15m3/s, 40m3/s,

60m3/s y 80m3/s, se puede evidenciar que la velocidad de agua es directamente
proporcional al caudal de agua que cursa en un vertedero natural o artificial.

4.10 CARTOGRAFIA DE RIESGOS DE DESBORDE

El modelo digital de terreno y las ortoimágenes generadas del levantamiento topográfico a

detalle del área objeto de estudio, constituyen la base para la conformación de la base de
datos con información geográfica para la generación de la cartografía de riesgos de
desborde. Adicionando a esta base las coberturas de la información geoespacial
generadas en el modelamiento y procesamiento de datos con IBER para el caudal pico de
57m3/s y la información geoespacial de los caudales de referencia; 75 m3/s y 100 m3/s.

Página 107
4.10.1 FACTORES DE RIESGO

Velocidad del agua; Es un factor de riesgo importante a considerar, particularmente por el

tema de erosión, material de arrastre, sedimentación y formación de abanicos aluviales en
lugares donde la velocidad del agua es menor por efectos de pendiente del lecho del rio.
Entonces, en el presente caso de estudio puente encuentro de los ríos Choqueyapy e
Irpavi luego del último aliviadero existe un salto hidráulico de unos 25 m.
aproximadamente. La sedimentación del material de arrastre en el lecho del rio podría
activar el riesgo de la disminución de la capacidad hidráulica del cauce por acumulación
de material de arrastre del rio.

Desborde de agua; El desborde del rio por insuficiencia de la capacidad hidráulica del
cauce del rio podría generar inundaciones que afectarían en el presente caso a la
infraestructura urbana y la convivencia de la sociedad colindante en su conjunto.
Debiendo realizarse periódicamente la remoción y/o evacuación del material de arrastre y
la evaluación de la capacidad hidráulica del cauce del rio La Paz (Norte), particularmente
en el tramo de objeto del presente estudio.

4.10.2 PRESENTACION GRAFICA DE LOS RESULTADOS CAUDAL PICO 57m3/s

El software IBER en el procesamiento de datos genera los resultados y/o los mapas
espaciales en formato .txt (cota de agua, velocidad del agua y otros). Esta información
para su incorporación a la base grafica debe ser convertidos a formato RASTER. En el
presente caso de estudio se ha utilizado el ArcGis Tools; Conversion Tool – to raster –
ascci to raster. Luego exportar dato en formato .tff para incorporar la información en la
base de datos grafica del área de estudio cuya presentación cartográfica de los resultados
se encuentra en Anexo 5. El resultado del procesamiento de datos del factor Velocidad de
Agua para una avenida máxima de 50 años de retorno con un caudal pico de 57m3/s, se
puede observar en el mapa. Ver figura 4.35 Mapa de Velocidad de Agua para un caudal
pico de 57 m3/s.

Página 108
Figura 4.30 Mapa de Velocidad de Agua para un caudal pico de 57 m3/s.

La presentación grafica del resultado referido al factor Cota de Agua para una avenida
máxima de 50 años de retorno con un caudal pico de 57m3/s, se puede observar en el
mapa. Ver figura 4.36 Mapa de cota de Agua para un caudal pico de 57 m3/s.

Página 109
Figura 4.31 Mapa de cota de Agua para un caudal pico de 57 m3/s.

4.11 PROPUESTA DE MEDIDAS DE MITIGACION

Conceptualmente la gestión de riesgo tiene como componente; El conocimiento de riesgo

y las estrategias de prevención y mitigación. En el presente caso de la evaluación

Página 110
hidráulica del tramo del rio objeto de estudio se tiene conocimiento del riesgo, el cual
sugiere en este caso la elaboración de la propuesta de medidas preventivas y mitigación
de riesgos.

La evaluación hidráulica del rio La Paz; tramo puente encuentro de los ríos Choqueyapu e
Irpavi hasta la confluencia con el rio Huayñajahuira, para una avenida máxima de 50 años
de retorno, Así como los resultados del modelamiento hidrodinámico de la capacidad
hidráulica del tramo del rio objeto de estudio constituye un elemento importante para
orientar la propuesta de gestión de riesgos por desborde del rio. De la evaluación de la
capacidad hidráulica se tiene conocimiento que el riesgo no es eminente y que se debe
formular una propuesta de prevención y mitigación de riesgos.

4.11.1 OBJETIVO DE LA PROPUESTA DE PREVENCION Y MITIGACION DE RIESGO

Realizar el seguimiento y control a los efectos del desarrollo de los factores de riesgo
identificados como; cota de agua que podría variar por efectos de la deposición de
material de arrastre y ocasionar el desborde del rio, y la velocidad de agua que incide en
la suspensión y deposición de material de arrastre en el lecho del rio.

4.11.2 ACTIVIDADES DE PREVENCION Y MITIGACION DE RIESGOS

- Mantenimiento del camino de acceso para la maquinaria pesada al cauce del rio
tramo puente encuentro de los ríos Choqueyapu e Irpavi hasta la confluencia con
el rio Huayñajahuira.
- Remoción de material de arrastre del cauce del rio hacia los margines del rio con
maquinaria pesada. Preferentemente cada año después de la época de lluvia.
- Evacuación y/o explotación de los áridos seleccionados como material de
construcción para la infraestructura urbana pública o privada.
- Evaluación de la capacidad hidráulica del tramo objeto estudio de forma periódica,
fundamentalmente antes y después de la época de lluvia anual.
- Alternativamente el GAM deberá prever la canalización del rio mediante la
construcción de vertedero con vías de acceso al lecho del rio para la remoción y/o
explotación de material de arrastre o áridos. Previo estudio de factibilidad.

Página 111
4.11.3 PRESUPUESTO

El Gobierno Autónomo del Municipio de Nuestra Señora de La Paz, deberá planificar en el

POA anual, la realización de los trabajos de prevención y mitigación de riesgos indicadas
anteriormente e inscribir el presupuesto correspondiente.

4.11.4 CRONOGRAMA

Actividad ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
Mantenimiento de xxxx xxxx
vías de acceso al rio
Remoción del xxxx xxxx
material de arrastre
Evacuación y/o xxxx xxxx xxxx xxxx xxxx xxxx xxxx xxxx
explotación de áridos.
Evaluación de la xxxx xxxx
capacidad hidráulica
del tramo

Página 112
 CAPITULO V

CONCLUSIONES Y RECOMENDACIONES

CONCLUSIONES

Mas del 50% de los espacios geográficos de las cuencas de los ríos Irpavi, Achumani y
Huayñajahuira se encuentran en zonas geológicas que resultan de procesos
geomorfológicos recientes Siluricos y Triasicos en la que los suelos son de diversos
grados de consolidación y diagénesis, donde el material predominante es suelto y frágil a
la acción erosiva de la precipitación pluvial como son los suelos compuestos en su
generalidad de arenisca, arcilitas, lutitas y lavas poco resistentes a la erosión pluvial. Por
estas características las cuencas de los ríos Irpavi y Achumani anualmente en época de
lluvia arrastran gran cantidad de material como; piedras, clastro, arena gruesa y otros
áridos.

De acuerdo a la caracterización de pendientes de la cuenca del Rio La Paz (Norte) la

superficie topográfica de la cuenca es bastante accidentada con pendientes superiores a
5%, lo cual hace que las sub-cuencas; ríos Choqueyapu, Orkojahuira, Irpavi, Achumani y
Huayñajahuira reaccionen de forma rápida a las fuertes precipitaciones, generando
escurrimiento caudaloso que erosiona fácilmente el suelo; estos tres últimos ríos
incorporan en la escorrentía material de suspensión y de arrastre.

Los ríos de Irpavi y Achumani arrastran gran cantidad de material a través de los
vertederos abiertos hasta unos 67 m. aguas abajo del puente encuentro de los ríos Irpavi
y Choqueyapu, lugar donde se presenta un salto hidráulico de 25 m aproximadamente,
aquí la escorrentía pierde fuerza debido a la reducción de la pendiente y dispersa el
material de arrastre formando abanico aluvial que aumenta el ancho del curso del rio.

Mediante el uso del modelo de transformación de las precipitaciones en escurrimiento

pluvial HEC-HMS con datos de la lluvia del proyecto generada con la curva IDF derivada
con datos de la Estación San Calixto (SENAMHI), parámetros morfométricos de las sub-
cuencas, se tiene el caudal pico de diseño en el puente encuentro de los ríos de
Choqueyapu e Irpavi de 56.7 m3/s.

Página 113
El modelamiento de la capacidad hidráulica con IBER del curso del rio La Paz (Norte) del
puente encuentro entre los ríos Choqueyapu e Irpavi hasta la confluencia con el rio
Huayñajahuira, se ha llevado a cabo bajo régimen permanente, flujo supercrítico
introduciendo parámetros del rio como la rugosidad y el caudal pico resultante del
modelamiento hidrográfico de superficie con HEC-HMS. Considerar que la lluvia del
proyecto deriva de los datos de máxima precipitación en 24 horas registrada en la
estación meteorológica San Calixto reflejada en la curva IDF para 50 años de retorno.

Las ortoimagenes elaboradas en el marco del presente estudio permiten identificar y

representar gráficamente el espacio del cuerpo de agua para el caudal máximo del año
2017 a través de las huellas de la crecida del rio en ambos márgenes. El resultado del
modelamiento dinámico con IBER con el caudal pico de 57 m3/s (corrida_q57_1_4) en lo
corresponde a la Cota de Agua y Velocidad de Agua sobrepuesta a la representación
gráfica del cuerpo de agua para el caudal máximo del año 2017 supera en pequeña
proporción y no muestra tendencia a desborde del rio. Entonces la capacidad hidráulica
del cauce del rio La Paz; tramo puente encuentro de los ríos Choqueyapu e Irpavi hasta la
confluencia con el rio Huayñajahuira, para una avenida máxima de 50 años de retorno es
satisfactorio y no existe riesgo eminente de desborde del rio.

RECOMENDACIONES

Para la remoción del material de arrastre hacia ambos márgenes del rio y una eventual
explotación anual de los áridos por parte de los concesionarios autorizados. El GAMLP
debe realizar el mantenimiento de las vías de accesos principales y alternas para ingreso
y salida de la maquinaria pesada. La remoción de material de arrastre del vale del rio es
sumamente importante.

Los resultados del modelamiento hidrodinámico de la capacidad hidráulica del tramo del
rio objeto de estudio constituye un elemento importante para orientar la propuesta de
gestión de riesgos por desborde del rio. De la evaluación de la capacidad hidráulica se
tiene conocimiento que el riesgo no es eminente y que se deben realizar acciones de
prevención y mitigación de riesgos como: Mantenimiento del camino de acceso para la
maquinaria pesada al cauce del rio, Remoción de material de arrastre del cauce del rio
hacia los márgenes del rio a efectos de rencausar y mantener el cauce natural del rio.

Página 114
En el lugar donde se presenta el salto hidráulico la escorrentía pierde la fuerza debido a la
reducción de la pendiente y se dispersa gran cantidad de material de arrastre formando
abanicos aluviales y eventualmente llanuras aluviales. La sedimentación del material de
arrastre en el cauce natural del rio aumenta el ancho del curso del rio y disminuye la
capacidad hidráulica. Por lo que se recomienda realizar la evaluación de la capacidad
hidráulica del tramo objeto estudio de forma periódica antes y después del periodo de
lluvia.

La alternativa de diseño y construcción del vertedero abierto o mantener el actual

vertedero natural del tramo objeto de estudio, debe contar con vías de acceso
permanentes debidamente consolidadas para el ingreso y salida de maquinaria pesada,
requerida para la realización de la actividad de remoción del material de arrastre y
eventual explotación de áridos por los concesionarios. Tomando en cuenta que el rio
Irpavi deposita anualmente gran cantidad material de arrastre.

Página 115
 CAPITULO VI

BIBLIOGRAFIA

Ahenke O. José Luis, “Test de Diversos Modelos Matemáticos para Determinar

Crecidas (Aplicado a la Cuenca de La Paz)”, UMSA, 1992.

Aparicio, Principios y Fundamentos de la Hidrología de Superficial, 1989

Arnolt et. Estimacion del Balance Hidrico y Producción de sedimentos bajo tres
escenarios de cobertura en la subcuenca del Rio Jabonal, Costa Rica mediante el empleo
de SWAT. 1990

Ayala Carcedo, F. J. (1987). “Introducción a los riesgos geológicos” Riesgos Geológicos;

I.G.M.E. Madrid. Vol. 1, 3-21.

Juan P. Celemin, Elaboración de la Cartografía de Riesgo de Inundaciones y Propuesta

de mejora de conservación en la Cuenca y Reserva MAB de Mar Chiquita, provincia de
Buenos Aires Argentina. 2009.

Francisco Coronado; Diseño y construcción de canales, Universidad Nacional de

Ingeniería de Lima Perú, 1992

E.P.T.I.S.A., “Anteproyecto de la Autopista La Paz – Lipari, Parte 2º Canalización

del río Choqueyapu”, Tomo 3, 1970.

J. González, O. Chávez, and M. Hermelín, «Aspectos geomorfológicos de la avenida

torrencial del 31 de enero de 1994 en la cuenca del río Fraile y sus fenómenos
asociados,» Desastres de origen natural en Colombia 1979-2004, M. Hermelín, ed.,
Medellin: Universidad Eafit, 2005

Etsa Ayesa G; Diagnostico Hidráulico, 1993

Jose L. Gomez Navarro; Saltos de agua y presas de agua, 1982.

IGAC, Producción y Actualización de la información geográfica Basica Digital

2005

HAM – GTZ, , “Proyecto de Regulación de las Cuencas Irpavi y Achumani”,

(Hidrología), 1990.

HAM – GTZ, “Plan de prevención contra inundaciones para la ciudad de La Paz,

periodo 1994-1995”, Cuencas Choqueyapu

Lahmeyer International GMBH – Prudencio, Claros y Asociados, “Agua Potable La

Paz – Ampliación de Embalses de la Ciudad de La Paz – Sistema Pampahasi –

Página 116
Nueva Presa Incachaca”, 1987.

Memoria de Proyecto; Metodología para la Elaboración de Cartografía de Peligrosidad por

Avenida e Inundaciones. 2009.

Nipponkoeico & PCA Consultores; Revisión y actualización del plan maestro de drenaje
para el área urbana de La Paz, 2007
Norma Boliviana NB688 -07; Diseño de Sistemas de. Alcantarillado Sanitario y Pluvial,
2001
Obermeyer y Pinto; Los Sistemas de Información Geográfica (SIG) y la hidrología, 1996

Santiago A Ochoa. Evaluación de modelos hidrodinámicos para representar flujos en

cauces naturales. Aplicación en un tramo del Río Suquía que incluye la confluencia con el
arroyo la cañada, Provincia de Córdoba; Argentina 2014

PCA Ingenieros Consultores S.A., “Estudio a Nivel de Diseño Final Rehabilitación

Plataforma Vial Av. Costanera Sector La Florida”, 2003.

Prudencio Claros y Asociados – SAE Sociedad Argentina de Estudios, “Diseño

Final Autopista La Paz – Aranjuez”, 1979.

Remeneiras G. Tratado de Hidrología Aplicada 1974.

Salgado Rabadan et; Hidrología superficial y Modelaje hidrológico de escenarios para

eventos de inundaciones, 2005

Salzgitter Consult GMBH, “Agua Potable La Paz Proyecto Saneamiento Sistema

Pampahasi – Estudio de Prefactibilidad Ampliación de la Capacidad de Embalses”,
1986.

Segaline Nieto Hugo, Cabré S.J. Hugo, “El Clima de La Paz”, Datos del
Observatorio San Calixto, 1988.

Universidad de las Américas, Recursos hídricos y la modelación hidrológica y de cuerpos

de agua, 2000

US Army Corps of Engineers, “Hydrologic Modeling System HEC-HMS – Technical

Reference Manual”, 2000.

US Army Corps of Engineers, “Hydrologic Modeling System HEC-RAS – Technical

Reference Manual”, 2000.

Veissman W. Introduction to Hidrology, 1989

Ven Te Chow, “Hidrología Aplicada”, 1987.

Ven Te Chow. Hidráulica de Canales Abiertos, 1994.

Página 117
 ANEXO -1
MÁXIMA PRECIPITACIÓN EN 24 HORAS
ESTACION SAN CALIXTO
Precipitacion maxima en 24 horas en
mm.
Datos
SENAMHI
AŃO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
1917 20.0 15.4 12.4 5.8 0.0 0.1 15.5 3.2 10.5 4.0 21.5 33.0 33.0
1918 18.0 22.5 15.5 0.9 7.0 9.7 19.2 1.0 19.5 11.0 9.5 38.0 38.0
1919 9.6 31.0 17.7 8.5 1.9 2.0 5.1 2.3 3.9 8.5 40.4 27.2 40.4
1920 17.5 21.3 27.1 8.9 3.8 6.3 9.5 13.0 16.2 9.9 16.5 16.0 27.1
1921 14.3 18.9 14.4 5.6 0.0 0.0 2.2 3.0 1.5 8.0 8.8 15.8 18.9
1922 14.0 16.2 11.3 14.0 2.0 1.5 0.0 6.5 15.4 14.0 20.0 32.7 32.7
1923 20.5 9.0 4.4 8.8 10.0 4.7 3.5 6.6 8.3 17.1 16.5 9.6 20.5
1924 15.2 20.0 4.6 6.8 8.0 5.5 0.0 0.2 4.5 8.4 5.7 11.0 20.0
1925 16.5 7.0 13.2 15.0 3.4 1.4 8.5 4.8 1.4 9.2 7.3 11.0 16.5
1926 27.2 16.5 25.2 19.5 4.1 1.7 0.0 5.0 9.7 2.4 4.0 14.0 27.2
1927 17.3 22.4 29.0 6.2 7.8 7.7 1.0 7.6 8.2 12.2 8.0 36.0 36.0
1928 16.6 25.3 24.5 13.0 10.4 3.5 0.2 7.0 28.5 7.8 14.3 18.0 28.5
1929 42.0 29.7 25.0 8.8 7.5 7.0 20.0 7.5 4.3 17.7 35.0 15.5 42.0
1930 27.5 25.8 11.0 29.7 10.8 4.6 5.5 12.0 19.5 8.1 12.5 12.0 29.7
1931 39.5 16.2 20.5 11.5 0.0 3.3 0.0 1.4 13.0 6.4 5.2 22.6 39.5
1932 36.5 16.7 13.5 20.0 18.7 0.0 12.5 6.4 6.4 12.5 25.0 32.6 36.5
1933 29.0 33.5 15.0 7.5 4.5 0.0 3.6 2.0 4.0 8.3 2.4 19.0 33.5
1934 29.5 25.0 10.4 11.5 5.2 1.4 6.5 7.5 10.0 5.0 22.5 10.0 29.5
1935 28.0 12.5 16.7 2.5 0.6 3.0 2.6 0.4 32.0 4.3 14.5 13.0 32.0
1936 14.1 15.3 5.8 3.5 6.9 0.0 4.8 5.2 14.6 11.0 26.5 17.2 26.5
1937 28.2 13.0 21.0 9.5 5.8 0.3 4.4 3.4 3.4 9.8 5.2 25.0 28.2
1938 26.0 40.0 11.0 7.0 0.3 6.4 2.0 1.3 7.8 4.3 14.0 19.6 40.0
1939 13.0 12.7 14.0 17.0 7.0 0.0 10.0 20.0 5.2 5.5 13.0 22.0 22.0
1940 20.0 5.0 9.0 4.0 5.5 2.4 0.6 1.5 4.5 12.5 5.2 21.5 21.5
1941 27.5 26.6 7.5 8.5 5.0 3.0 0.0 0.0 11.6 12.0 22.6 11.0 27.5
1942 10.6 18.5 7.0 30.0 0.0 **** 0.0 2.0 2.0 20.0 **** 10.0 ****
1943 20.0 12.7 8.7 26.3 16.5 1.9 0.0 4.8 5.0 6.3 10.6 27.5 27.5
1944 15.0 21.5 11.0 15.0 4.3 0.0 0.6 0.0 25.0 25.0 17.0 16.3 25.0
1945 19.0 21.8 16.0 11.7 0.0 0.0 5.0 4.8 28.0 18.0 22.5 16.0 28.0
1946 8.8 28.0 11.0 4.1 11.0 0.2 21.0 6.7 23.7 8.7 14.2 17.0 28.0
1947 49.0 33.0 13.0 7.8 7.0 0.1 3.0 2.0 11.0 14.0 8.0 10.0 49.0
1948 27.7 15.0 14.0 14.8 15.0 4.2 6.0 2.7 2.5 10.0 30.0 24.0 30.0
1949 25.5 12.5 16.2 11.5 0.5 12.0 0.3 3.0 26.0 13.7 24.0 11.6 26.0
1950 24.0 9.3 11.0 8.0 5.8 0.0 6.6 7.0 12.0 10.0 18.6 11.6 24.0
1951 26.0 16.3 16.0 5.3 2.6 0.0 0.0 8.0 8.5 20.5 12.8 24.5 26.0
1952 27.0 42.6 4.2 8.6 7.0 5.0 9.5 2.3 18.5 8.4 12.6 40.0 42.6
1953 29.0 13.0 10.0 8.0 1.0 2.5 0.0 8.0 5.1 18.0 14.6 14.5 29.0
1954 26.0 31.0 20.0 8.0 7.0 3.0 0.0 2.0 15.0 15.0 29.4 16.3 31.0
1955 21.2 20.0 30.0 6.4 4.5 2.7 4.0 4.0 12.5 10.0 7.0 23.7 30.0
1956 29.5 13.3 7.8 4.0 0.0 0.0 4.8 4.0 7.5 3.5 6.5 12.0 29.5
1957 38.0 18.2 18.0 17.6 2.7 20.0 0.0 4.2 16.5 20.0 27.0 23.0 38.0
1958 23.0 21.5 20.0 4.0 10.5 0.0 5.5 10.5 7.0 10.0 5.7 16.0 23.0
1959 17.0 20.0 32.0 12.0 11.0 2.0 0.0 3.5 15.9 6.5 8.5 21.5 32.0
1960 27.8 21.4 13.8 17.3 5.5 0.6 0.0 11.6 16.2 10.5 22.0 24.8 27.8
1961 35.2 23.1 10.3 23.5 9.6 1.0 0.0 2.2 14.1 12.2 9.1 26.0 35.2
1962 10.1 15.4 16.0 17.7 2.2 1.6 1.0 4.4 8.1 10.7 11.5 16.5 17.7
1963 18.0 16.2 12.6 16.0 4.2 1.9 0.0 3.8 10.3 5.3 8.9 17.8 18.0
1964 14.3 18.3 22.7 10.9 3.0 0.0 2.7 5.1 13.8 12.4 13.4 17.8 22.7
1965 41.0 37.2 17.9 8.3 11.0 0.0 3.0 2.2 3.4 8.2 27.6 18.4 41.0
1966 10.4 18.2 6.6 8.4 20.7 0.0 1.9 0.0 4.4 20.8 15.5 13.9 20.8
1967 10.8 15.5 23.7 8.3 4.3 0.8 9.5 10.0 11.3 19.0 7.1 14.4 23.7
1968 30.2 24.7 14.0 7.2 10.6 4.4 5.5 5.4 12.0 6.0 21.4 30.5 30.5
1969 16.2 16.8 24.3 10.7 3.6 1.6 6.1 3.0 4.0 5.4 13.5 15.0 24.3
1970 29.7 23.5 10.0 8.6 11.5 0.0 1.0 0.8 7.3 8.5 9.5 26.7 29.7
1971 13.7 22.0 3.4 6.5 1.4 3.7 0.0 4.0 0.9 13.5 28.1 22.9 28.1
1972 24.0 7.4 14.5 17.2 0.8 1.2 0.0 11.9 9.0 8.6 10.5 33.8 33.8
1973 20.5 17.4 15.1 12.0 1.8 1.9 18.9 9.0 16.0 8.9 18.3 23.1 23.1
1974 27.3 16.4 11.2 5.4 4.0 6.3 0.4 12.9 4.4 13.7 8.6 17.6 27.3
1975 24.4 15.3 24.1 6.4 7.7 3.3 0.0 2.2 5.3 7.7 10.0 16.2 24.4
1976 32.4 19.7 11.0 1.4 14.2 7.7 5.2 3.8 26.6 3.7 5.2 50.1 50.1
1977 17.0 45.8 32.7 1.7 9.1 0.1 13.0 7.3 8.7 10.3 30.5 39.4 45.8
1978 16.2 24.5 8.4 19.0 0.9 0.0 0.3 1.7 4.7 2.3 18.9 12.6 24.5
1979 16.2 22.2 22.3 6.8 4.4 0.0 5.1 1.3 8.8 15.7 21.9 24.8 24.8
1980 7.7 16.8 26.5 5.7 5.4 3.4 6.8 27.7 12.5 22.4 12.0 16.0 27.7
1981 22.8 29.1 17.8 20.5 2.2 0.0 0.0 9.4 22.8 14.2 7.7 7.0 29.1
1982 35.2 23.0 24.2 11.9 0.0 4.1 1.0 2.0 10.9 11.1 26.8 27.0 35.2
1983 19.1 8.8 30.0 7.5 6.6 7.0 9.8 7.0 18.2 9.7 2.7 25.3 30.0
1984 33.7 17.6 39.1 5.1 3.0 3.2 0.2 5.8 1.3 5.5 27.3 15.1 39.1
1985 19.7 20.5 19.2 20.0 2.9 9.0 0.3 1.0 21.0 18.7 17.9 21.2 21.2
1986 39.6 20.7 16.9 7.5 9.1 0.0 4.0 25.0 17.7 7.8 15.2 24.6 39.6
1987 25.8 7.8 24.4 3.0 6.3 0.3 4.8 8.1 10.7 26.0 17.7 10.1 26.0
1988 35.8 26.8 22.0 22.3 16.6 0.0 1.2 0.0 8.2 10.6 13.0 15.6 35.8
1989 16.0 30.1 7.1 14.1 11.6 2.0 5.0 1.6 7.5 2.3 8.1 8.1 30.1
1990 30.9 26.7 3.0 12.3 12.8 22.5 0.0 2.0 12.3 5.9 20.6 29.9 30.9
1991 16.7 14.5 12.3 7.3 13.0 19.1 0.5 0.6 6.8 1.5 17.0 34.1 34.1
1992 **** **** **** **** **** **** **** **** **** **** **** **** ****
1993 **** **** **** **** **** **** **** **** **** **** **** **** ****
1994 **** 10.2 17.0 6.8 7.5 8.3 2.2 0.0 6.9 6.0 18.5 14.0 ****
1995 22.0 41.2 20.3 7.9 2.6 0.0 4.8 10.4 6.7 9.5 16.8 43.4 43.4
1996 30.6 9.0 14.3 10.3 0.0 1.2 5.0 3.0 9.8 7.0 17.8 31.6 31.6
1997 44.8 20.3 18.5 18.6 9.8 0.0 0.0 2.5 25.6 21.2 19.0 19.0 44.8
1998 17.2 15.8 19.5 34.0 0.0 13.5 0.0 5.3 0.4 14.5 32.8 12.7 34.0
1999 25.5 17.4 18.0 14.8 1.8 1.5 2.7 1.0 28.5 13.0 10.0 27.8 28.5
2000 17.0 9.5 20.5 18.0 6.0 7.5 0.5 9.5 2.2 27.5 3.5 19.3 27.5
2001 22.5 25.7 15.1 8.0 9.5 4.0 8.5 11.5 9.5 19.5 5.7 14.5 25.7
2002 32.8 59.9 22.5 7.4 11.6 8.8 13.0 3.5 7.8 26.2 18.0 14.0 59.9
2003 54.0 18.0 20.7 5.8 1.7 6.6 4.2 6.3 5.1 3.6 3.6 31.5 54.0
2004 20.0 10.0 14.0 5.3 5.3 3.6 6.3 14.0 5.5 6.1 16.5 13.8 20.0
2005 25.5 11.9 14.3 9.8 0.0 0.0 1.5 0.8 10.9 14.2 13.9 15.2 25.5
2006 25.2 30.0 13.0 2.8 1.9 0.0 0.0 15.4 5.2 13.0 39.6 11.2 39.6
2007 20.0 21.0 9.8 15.2 4.5 0.0 17.8 0.0 10.5 27.2 19.8 36.0 36.0
2008 21.8 20.8 16.3 2.5 3.2 3.5 2.2 3.2 6.2 10.4 7.0 23.8 23.8
2009 20.3 33.7 9.6 20.0 0.5 0.0 12.6 2.2 13.4 11.0 20.0 28.2 33.7
2010 20.2 17.8 7.5 1.0 10.3 0.0 2.2 3.2 8.8 0.0 4.6 20.0 20.2
2011 18.8 34.3 8.9 0.7 4.1 0.2 5.9 0.3 6.8 8.3 21.9 35.0 35.0
2012 30.4 18.0 15.2 7.5 2.3 7.4 2.3 0.3 1.2 14.2 13.8 12.8 30.4
2013 20.2 41.0 7.8 11.3 5.0 5.6 7.2 12.0 1.6 13.8 10.2 24.3 41.0
2014 18.6 17.8 9.8 9.3 11.0 0.0 2.6 5.6 6.2 20.2 3.6 20.8 20.8
2015 26.7 14.6 18.8 14.3 5.2 0.0 10.6 9.2 12.9 18.9 9.1 18.4 26.7
2016 19.4 17.7 8.0 8.9 0.0 2.2 1.0 9.2 9.0 14.2 6.8 36.3 36.3
 ANEXO -2
PRECIPITACIÓN MEDIA ANUAL
ESTACION ACHACHICALA
Precipitacion media anual
Datos SENAMHI
AŃO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
1975 **** **** **** **** **** 0.0 0.0 4.2 43.8 37.3 56.0 136.3 0
1976 225.1 153.6 65.2 3.1 48.9 17.0 10.3 19.9 100.0 7.8 15.8 101.7 768
1977 88.9 160.9 209.6 9.0 35.1 0.0 19.6 11.8 55.5 43.0 135.4 194.3 963
1978 194.2 143.8 41.0 36.7 2.6 0.0 2.8 9.9 15.8 17.3 139.4 152.4 756
1979 265.8 67.5 160.3 29.5 4.3 0.1 2.2 2.2 16.4 112.0 68.0 186.1 914
1980 61.0 72.8 98.8 15.8 8.9 4.6 10.3 32.0 57.0 81.8 31.9 52.4 527
1981 212.5 159.6 76.6 74.8 0.0 0.0 0.0 34.8 52.1 78.3 59.7 73.5 822
1982 93.1 42.5 11.3 13.0 0.8 10.2 0.0 6.8 40.0 57.8 82.3 66.6 424
1983 69.1 65.5 55.4 32.1 27.6 9.1 13.1 14.3 42.9 31.8 24.3 141.0 526
1984 200.7 151.1 100.0 5.3 14.0 6.8 2.1 14.5 6.1 35.3 82.4 83.1 701
1985 179.0 136.5 46.1 89.3 19.6 8.5 2.0 2.3 38.1 50.2 90.8 143.9 806
1986 103.8 146.8 138.7 46.9 17.1 0.0 3.8 31.4 82.9 41.0 63.3 157.2 833
1987 194.4 16.7 50.3 11.4 13.6 6.5 14.1 8.3 27.5 48.1 103.9 44.3 539
1988 180.2 74.8 133.9 57.5 32.6 0.0 3.9 0.0 27.5 42.4 30.1 99.0 682
1989 90.3 67.6 36.8 52.1 21.5 4.3 7.8 9.0 19.2 12.4 27.8 113.6 462
1990 228.1 44.7 21.6 46.0 38.6 73.5 0.5 8.0 23.6 31.5 80.4 114.3 711
1991 115.2 93.6 40.2 9.9 15.4 38.3 0.6 1.3 18.3 16.0 48.4 72.5 470
1992 140.4 88.1 29.9 27.5 0.0 2.5 **** **** **** **** **** **** 0
ESTACION ACHUMANI
Precipitacion media anual
Datos SENAMHI
AŃO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
1991 105 32.4 60 14 19 21.5 0 7.9 12.7 6.8 47.1 50.4 377.2
1994 62.8 53.9 37.4 64.5 9.5 15.5 1.3 2.3 25 38.3 44.6 66.7 421.8
1995 87.5 66.6 99 20.6 3.8 0 5.2 15.6 0.8 14.4 52 177 542.4
2001 167 93.5 66 29.7 18.5 5.2 9.6 31 14.8 61 13.1 87.6 596.5
2002 79.6 104 109 29.7 13.6 9.7 37.5 9.2 37 95.2 37.5 70.9 633.3
2003 137 124 85.6 16.6 5.8 0.6 2.3 7.6 27.6 28.9 23.2 100 558.9
2004 212 69.4 49.6 8.7 6 5.3 15.3 13.3 26.2 18.3 77.9 68.5 570.3
2005 110 59.9 29.6 22.3 0 0 1 0 46.5 76.4 76.2 52.8 474.2
2006 174 89.6 69.7 18 2 0 0 17 13.2 34.1 77.3 122 616.5
2007 124 98.4 93.7 61.8 12.2 0 28 0.1 43.4 21.1 38.7 140 661.4
2008 212 66.4 37.7 3.4 12.6 9.8 5.7 4.5 10.6 35 20.8 118 536.2
2009 88.2 115 37.6 20.5 3.6 0 9.9 3.6 27.2 36.8 60.8 110 513.9
2010 194 191 66.6 21.4 58 0 19.6 4.8 5.1 39.6 4.3 53.5 657.7
2011 51.1 129 55.7 0.3 10.5 0 12.8 0 24.3 35 13.4 81.7 413.7
2012 83.5 133 50.9 43.6 1.7 18.1 5.7 1.2 1 24.6 66.6 80.6 510
2015 128 93.3 94.6 47.4 4.9 0 12 14.4 26.3 45.3 45.7 50.4 562.6
2016 63.9 154 4.6 20.6 0 4.7 1.4 14.1 12.6 60.5 15.9 70.9 422.7
ESTACION ALTO IRPAVI
Precipitacion media anual
Datos SENAMHI
AŃO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
1991 135.1 91.6 97.8 27.9 11.0 22.8 1.8 2.7 24.3 25.9 66.8 64.8 573
1992 112.2 53.9 20.9 21.4 0.0 21.4 3.5 31.1 8.1 31.2 55.5 87.7 447
1993 133.4 38.6 43.6 45.9 6.1 4.2 2.4 29.2 23.8 62.9 39.7 103.1 533
1994 56.5 73.0 69.2 37.1 8.6 6.8 4.9 3.8 28.9 36.2 54.0 86.8 466
1995 122.5 36.4 88.3 10.4 14.0 0.0 2.2 3.0 27.6 14.3 71.0 68.9 459
1996 94.9 57.6 53.5 32.0 0.1 0.0 3.6 11.8 40.8 25.2 80.0 74.3 474
1997 113.5 92.4 122.3 25.2 1.5 0.0 0.0 3.1 42.6 41.4 34.8 23.4 500
1998 38.8 81.7 69.2 40.6 0.0 24.9 0.0 0.0 7.5 57.2 74.0 49.9 444
1999 91.5 75.9 118.9 47.1 0.0 6.5 0.0 3.4 30.9 50.0 28.4 93.9 547
2000 90.1 98.2 73.8 33.2 2.0 10.5 0.0 28.1 0.0 70.2 13.0 113.9 533
2001 209.7 96.0 130.3 37.6 32.1 4.3 2.0 12.6 8.3 30.1 21.1 76.3 660
2002 59.1 86.3 102.8 18.0 6.8 33.4 21.3 15.1 23.6 68.6 58.6 49.8 543
2003 195.8 99.2 102.8 15.6 1.9 1.0 18.2 8.1 20.6 8.9 48.1 119.1 639
2005 142.8 101.7 16.3 16.1 1.3 0.0 0.0 2.3 25.5 78.6 79.2 113.0 577
2006 174.2 75.2 45.5 23.9 4.8 1.7 0.0 21.6 22.7 38.3 78.3 119.2 605
2007 108.5 107.7 60.1 56.9 7.8 0.0 9.2 2.9 52.0 33.5 16.3 148.9 604
2008 192.8 79.6 84.9 10.7 10.1 7.5 0.0 7.9 4.9 39.8 26.9 87.9 553
2009 87.8 85.3 52.7 27.1 6.2 0.0 18.7 0.0 15.8 51.2 37.7 61.8 444
2010 132.5 100.3 48.7 10.5 19.8 0.5 3.5 3.4 15.7 24.4 5.6 40.6 406
2011 66.9 199.8 57.0 17.8 2.5 0.0 6.8 2.8 40.1 13.0 25.3 95.0 527
2012 93.8 103.7 49.8 57.6 0.0 13.8 5.2 0.8 4.5 26.9 100.5 158.1 615
2013 71.2 138.8 94.1 17.7 39.2 47.2 15.5 46.7 9.5 50.4 30.9 84.9 646
2014 138.8 78.4 85.1 51.4 12.4 2.3 19.3 32.8 85.4 14.7 27.3 97.1 645
2015 225.0 83.3 141.5 122.0 6.6 0.0 4.2 21.9 42.9 54.0 82.5 67.5 851
2016 87.5 121.1 10.8 98.0 1.2 10.0 9.8 6.7 8.1 22.9 35.4 114.5 526
ESTACION SAN CALIXTO
Precipitacion media
anual
Datos SENAMHI
AŃO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
1945 71.1 75 74.7 30.7 0 0 2.8 7.5 32.2 37.5 107 74.5 513
1946 79.5 142 48 8.5 24.5 0 30 28 37.5 20 33 119 569.5
1947 133 154 30.7 33 5.2 0 8 0 **** **** **** **** 0
1948 **** **** 68.2 38.5 25.2 8 **** **** **** **** **** **** 0
1949 108 55.1 40.1 23.8 **** **** **** **** **** **** **** **** 0
1950 57.3 19.8 **** **** **** **** **** **** **** **** **** **** 0
1951 **** **** **** **** **** 0 0 6.3 3.5 12.7 19.7 **** 0
1952 7.5 104 15.6 0 **** 6.4 23.4 2.5 52.3 0 36.9 39.6 0
1953 113 0 15 13 0 0 **** 5.9 7.5 35 81.8 53.5 0
1954 86.3 98.4 82.5 18 2 2 1 0 37.8 4.8 82.9 38.1 453.8
1955 124 69.4 90.4 3.6 6 0.6 4.9 4.8 44.3 31.5 16.2 69.1 464.7
1956 110 49.7 30 1 2 0 3 8.5 18.3 16 22 68.5 329.2
1957 103 94.2 23 69.4 11.8 13 0 9 24.1 38.5 39.4 53.7 479.4
1958 91.4 61.9 51 2.3 16.5 0 8 0 13.4 21 37.3 69.7 372.5
1959 66.1 83.9 86 **** **** **** 0 3 26 20.4 30.5 105 0
1960 139 106 38.3 45.2 9.5 0 0 16.5 78.7 41.3 114 70 657.8
1961 96.5 107 23.9 50.6 17.9 1.9 0 3.2 24.8 25.2 41.3 117 509.3
1962 75.7 55.8 81.1 34 3.5 1.2 0.1 4.5 29.8 22.8 10.7 107 426.6
1963 137 83.5 61.9 36.4 5 2.1 0 0 70.2 18.6 23.9 75.1 513.3
1964 72.4 91.1 122 14.3 2.8 0 7.9 5.8 34 39.4 46.9 92 528.4
1965 128 111 69.8 29.8 10.3 0 5.6 0.3 33.5 19 78.6 99 585.2
1966 60.1 72.6 33.3 11 37.9 0 1.2 0 10.3 47.8 42.3 65.5 382
1967 38.2 38.8 68 1.3 8.4 0.5 15 19.1 32.3 44 6.9 109 381.8
1968 63.4 131 69.1 13.8 19.5 4.8 1.5 17 25.5 15 81.3 88 529.4
1969 116 65.8 51 21.8 2.5 2.6 6.6 2.8 3.1 16.2 25.4 65.2 379.1
1970 150 103 75.7 20.1 17.8 0 1.5 3 35.9 38 39.8 123 607.6
1971 92.7 138 31.9 17.3 2.3 8.1 0 5.2 0.7 46.4 42.2 80.7 465
1972 137 52.8 47.6 42.2 0.2 1.1 1.6 11.5 45.7 21.7 42.6 89.9 493.9
1973 110 81.7 80.2 26.2 4.7 1.6 26.7 13.9 36.4 20.3 30.1 83.1 514.4
1974 175 131 44.7 27.5 1.9 20 0.8 61.6 6.9 50.6 10.2 97.4 627.4
1975 165 105 96.9 23.6 19 4.2 0 7 36.4 33.2 34.1 102 625.9
1976 124 90.3 56.8 2.2 30.7 9.9 5.7 9.9 86.6 8.8 11.9 93.1 530.1
1977 63.8 136 171 5.7 10.7 0.6 22.3 9.7 29.4 31.2 103 121 703.5
1978 110 139 28.9 56 1.4 0 0.1 5.1 12.3 5.2 90.6 98 546.5
1979 167 80.6 141 16 5.2 0 7.4 1.5 20.8 110 63.4 142 753.8
1980 45.6 60.4 80.6 16.5 5.4 3.4 12.3 33.6 48.4 76.6 19.2 48.1 450.1
1981 146 145 61.5 44.9 3 0 0 35.7 56.5 62.6 33.9 33.5 622.2
1982 130 39.3 60.4 29 0 6.5 2.2 3.2 29.6 50.6 57.9 59.9 469
1983 61.7 25.2 39.1 28.8 19.2 7.5 8.7 14.2 29.7 25.3 7.1 80.3 346.8
1984 222 170 174 13 3 8.5 0 15.5 2.6 18.8 79 95.3 802.6
1985 100 139 62.2 68.8 2.7 16.8 0.6 1.6 56.8 23.9 91.2 108 672
1986 85.4 126 105 37.1 14.7 0 3 27.4 63.8 29.1 42.7 139 672.4
1987 185 26.7 36.4 9.2 15.9 0.4 9.5 8.4 23.1 57.4 57.2 28.6 457.9
1988 170 70.2 161 30.7 30.2 0 1.1 0 18.4 25.2 24.2 60.7 591.3
1989 88 84.5 28 35.6 20.7 2.2 7.8 8.1 14.7 8 14.3 39 350.9
1990 178 80.7 13.4 35.1 20.1 65.3 0.1 5 25.4 30.1 90.3 115 658.3
1991 111 85.9 58.3 15.9 17.6 27.7 0.5 1.7 15.7 3.3 44.6 55.4 437.4
1992 132 48.1 34.8 18.6 0 14.1 10.8 41.8 9.6 45.3 63.1 44.6 462.6
1993 114 20.1 45.2 42.6 1.8 4.4 2.6 37.8 25.9 43.7 34.6 115 487.3
1994 56 71.1 54.8 35.7 6.3 13.2 1.4 1.2 15.1 12.7 40 57.1 364.6
1995 96.5 86.4 57.8 11 3.7 0 4.9 11.8 8.4 10.8 38.8 114 444.2
1996 136 25.8 54.1 14.7 0.3 0.8 5.4 5.8 22.1 19.7 50.1 107 442
1997 160 75.3 129 43.5 17.6 0 0 5.9 51.4 37 49.3 43.3 612.3
1998 84.4 61.5 45.1 48 0 19.2 0 7.4 2.1 57.4 83.2 53.3 461.6
1999 106 48.4 109 33.8 5.5 1.2 1.2 1 51.1 65.1 26.5 73.4 522.5
2000 110 34.1 72 12.6 5 21.1 0.2 21.9 1.1 62.3 10.6 129 479.3
2001 180 107 77.7 20.4 25.7 6.5 9.3 31.9 17.8 67.3 9.8 75.7 628.9
2002 90.3 118 88.1 35.1 12.7 10.6 25.9 9.5 26.2 86.3 50.6 63.8 617.4
2003 135 97.9 82.8 6.9 5.8 1.6 2.4 8.9 27.2 44.3 9.3 106 527.4
2004 182 61.9 51 10.6 3.7 6.4 17.5 24.6 12.5 15.8 48.8 36.3 470.6
2005 125 74.8 35.9 35.9 0 0 2.3 0.4 30.8 51.7 60.2 63.4 480.6
2006 144 100 76.3 15.9 1.4 0.1 0 19.9 12.2 30.6 67.4 122 590.1
2007 104 79.3 60 42.9 4.5 0.1 21.8 0.1 46.1 49.2 35.6 91.3 534.4
2008 154 69.7 67.2 1.1 4.7 3.4 2.7 3.1 4.8 23.6 12.3 112 458.3
2009 61.4 116 53.4 30.5 0.3 0 15.1 1.5 25.7 27 75.7 137 544.1
2010 120 142 31.9 3.9 20.9 0 2.8 6.2 8.3 50.2 9.8 110 505
2011 57.1 157 52.4 0.8 8.2 0 10.7 1.7 32.8 22.5 25.5 128 496.6
2012 159 92.2 64.8 26.9 2.9 7.7 2.6 0.2 6.1 42.9 79.3 81.2 565.3
2013 64.4 130 30.8 11.2 23.3 14.3 7.8 20.4 4.8 37.9 33 99.9 477.6
2014 94.2 73.4 43.1 21.5 25.1 0 2.5 10.4 34.1 37.6 8.7 72.4 423
2015 149 64.5 77.9 63 6.1 0 10.2 17.4 33.3 51.4 40.9 91.2 604.6
2016 56.8 98.7 4.8 30.6 0 2.4 0.6 9.2 9.8 46.6 27.4 101 388
ESTACION VILLA COPACABANA
Precipitacion media
anual
Datos SENAMHI
AŃO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
1977 **** **** **** **** **** **** **** **** **** **** 62.8 120 0
1978 402 102 27.4 30 0 0 0 **** **** **** **** **** 0
1979 91.5 **** **** **** **** **** **** **** **** **** **** 153 0
1980 59 62.3 98.8 18.1 4.8 3.2 9.1 34.5 56.8 93.7 14.3 52.1 506.7
1981 195 157 71.8 52.4 3.7 0 0 40.6 61.2 68.5 38.2 71.7 760.3
1982 125 43.4 101 27.4 0 9.8 1 3 35.6 52.1 67.7 65.3 531.5
1983 76.5 42.5 41.9 32.3 28.4 6.7 8.7 13.4 37.5 20.2 13.6 133 454.5
1984 213 186 144 16.3 3 1.7 2.4 19.4 1.9 30.2 118 90.7 826.7
1985 135 135 65.8 87.2 4.3 18.3 0.5 4.6 66.6 34.2 98.3 144 794.1
1986 95 141 132 45.5 15.9 0 3.5 29.2 76.4 38.4 54.1 156 786.6
1987 235 32.5 50.3 19.1 17.3 2.9 13.9 11.1 36.1 56.1 66.6 49.1 590.3
1988 152 77 174 63 21 0 0.5 0 27.4 34.7 16.7 89.1 655.6
1989 104 93.2 21.1 54.1 25.6 0.5 12.9 12.1 19.8 8.7 29.5 87.2 468.7
1990 196 79.3 20.7 33.6 30 75.5 0 5.7 19.2 30.5 94.8 102 686.6
1991 145 75.3 50.1 23.8 26.4 34.5 0 1 27.5 10.8 60.4 85.5 540.4
1992 129 69.5 56 32.2 0 11.7 14.3 56.7 9.2 65.6 89.3 86.3 619.5
1993 149 39.3 93 74.2 6 4.4 5 49 38.8 62.5 54.3 113 688.1
1994 81.3 96.2 90.1 41.6 10.8 15.6 1.6 2.9 28.3 30.9 47.1 64.6 511
1995 138 95.9 90.1 27.5 6.3 0 5.6 14.4 10.6 15 66.4 110 579.6
1996 154 34.8 37 26.1 0 0.9 6.9 7.8 33 27.3 72.2 130 530.3
1997 234 113 192 47.5 13.2 0 0 9.7 62.7 52.5 54.1 51.4 830.6
1998 128 101 55.6 38.5 0 30.7 0 13.9 3.3 74.1 110 48.8 602.9
1999 103 70.4 128 35.1 3.8 0.8 7.9 1.3 60.9 72.9 22.2 89.3 595.6
2000 112 50.7 116 8.6 9.4 27.4 0.9 36.9 2.6 88.6 12 150 614.5
2001 202 152 87.1 20.8 29.5 8.4 9.9 41.5 21.3 77.9 9.1 94.4 754.4
2002 94.7 111 113 42 18.1 11.2 29.5 9.9 43.7 128 62.2 86.5 749.7
2003 134 129 129 12.8 6.1 4.9 2.1 3.7 30.1 33.1 13.8 130 628.4
2004 228 86.6 57.4 18 16.7 4.9 22.1 32.1 22 30.9 40.2 48.9 607.9
2005 144 106 37.8 28.4 0 0 0.6 0 48.6 48.7 66.4 88 568.2
2006 168 95.5 85.6 19 2 0 0 0 13.5 23.5 69 137 612.6
2007 129 87.9 60.7 47.2 0 0 23.7 0 32.8 46.7 67 105 600.4
2008 198 81.9 84.5 0.4 2.2 6.1 5.6 3.8 6.1 28.5 22 131 570.3
2009 78.4 116 36 23.1 1.5 0 18.5 0.1 26.6 25.9 84.5 140 550.9
2010 137 149 28.8 0.1 14 0 1.8 5 9.7 47 14.9 122 528.4
2011 69.8 149 67.3 1.8 5.1 0 8.2 0.1 24.3 19.4 26.3 154 525.9
2012 168 116 55.2 48.2 1.4 8.7 1.8 0.2 1.2 30.6 68.5 115 614.4
2013 93.5 188 67.9 13.4 25 21.4 9.6 22.2 0.9 51.5 54 119 667.1
2014 157 80.5 65.5 22.1 28.9 0 5.7 11.6 50.4 37.5 15.1 83 556.9
2015 165 77.2 89.8 55.2 8.6 0 11.9 22.2 42.7 53 46.5 95.6 668
2016 86.2 108 17.7 35.3 0 3.4 0 10.6 14.6 59.8 20.6 115 470.9
 ANEXO -3
PRECIPITACIÓN TOTAL ANUAL
ESTACIONES METEOROLOGICAS DE LA CUENCA LA PAZ (NORTE)
Precipitaciones Totales Anuales
Datoa
SENAMHI

CHUQUIAGUILLO
ALTO IRPAVI

VINO TINTO
ACHUMANI
ACHACHICALA

HUAYLIPAYA
ALTO OBRAJES B M

PASANKERI
ALTO ACHACHICALA

LAYKACOTA

VILLA COPACABANA
SAN CALIXTO
AŃO

1917 477.1
1918 587.8
1919 596.3
1920 680.5
1921 526.0
1922 650.2
1923 543.8
1924 355.6
1925 462.6
1926 510.6
1927 626.5
1928 595.0
1929 661.9
1930 741.9
1931 636.4
1932 813.4
1933 536.7
1934 643.6
1935 490.8
1936 466.6
1937 478.8
1938 489.4
1939 556.5
1940 344.9
1941 606.3
1942
1943 541.2
1944 621.7
1945 539.5 601.5 513.0 513.0 513.0
1946 573.2 569.5 569.5 569.5
1947 554.3
1948 706.7
1949 580.6
1950 412.2
1951 516.9
1952 576.7
1953 603.5
1954 686.5 171.1 453.8 453.8 453.8
1955 612.6 464.7 464.7 464.7
1956 396.4 329.2 329.2 329.2
1957 621.5 133.3 479.4 479.4 479.4
1958 487.8 372.5 372.5 372.5
1959 594.9
1960 730.9 657.8 657.8 657.8
1961 587.6 697.1 509.3 509.3 509.3
1962 483.9 480.1 426.6 426.6 426.6
1963 597.1 568.3 513.3 513.3 513.3
1964 590.6 575.1 528.4 528.4 528.4
1965 674.3 570.3 585.2 585.2 585.2
1966 425.5 382.0 382.0 382.0
1967 494.3 381.8 381.8 381.8
1968 608.7 529.4 529.4 529.4
1969 415.8 379.1 379.1 379.1
1970 621.4 607.6 607.6 607.6
1971 451.4 465.0 465.0 465.0
1972 583.8 493.9 493.9 493.9
1973 633.8 514.4 514.4 514.4
1974 683.6 627.4 627.4 627.4
1975 599.7 625.9 625.9 625.9
1976 716.3 768.4 530.1 530.1 530.1 636.3
1977 743.5 963.1 703.5 703.5 703.5 660.6
1978 535.6 755.9 546.5 546.5 546.5 602.9
1979 749.2 914.4 753.8 753.8 753.8 802.3
1980 453.0 527.3 450.1 450.1 450.1 506.7 586.3
1981 644.8 821.9 622.2 622.2 622.2 760.3 778.1
1982 481.5 424.4 469.0 469.0 469.0 531.5 512.7
1983 409.0 526.2 346.8 346.8 346.8 454.5 470.7
1984 813.0 701.4 802.6 802.6 802.6 826.7 883.0
1985 707.3 806.3 672.0 672.0 672.0 794.1 820.3
1986 750.9 832.9 672.4 672.4 672.4 786.6 820.3
1987 533.2 539.1 457.9 457.9 457.9 590.3 648.4
1988 643.8 681.9 591.3 591.3 591.3 655.6 705.2
1989 385.1 462.4 350.9 350.9 350.9 468.7 445.8
1990 648.3 710.8 658.3 658.3 658.3 686.6 678.0
1991 483.5 469.7 377.2 572.5 572.5 437.4 437.4 437.4 540.4 511.1
1992 446.9 446.9 462.6 462.6 462.6 619.5 532.4
1993 532.9 532.9 487.3 487.3 487.3 688.1 631.3
1994 421.8 465.8 465.8 364.6 364.6 364.6 511.0 468.3
1995 600.0 542.4 458.6 458.6 444.2 444.2 444.2 579.6 635.0
1996 507.4 473.8 473.8 442.0 442.0 442.0 530.3 569.0
1997 729.5 500.2 500.2 612.3 612.3 612.3 830.6 707.7
1998 550.7 443.8 443.8 461.6 461.6 461.6 602.9 515.5
1999 614.4 546.5 546.5 522.5 522.5 522.5 595.6 594.2
2000 488.7 533.0 533.0 479.3 479.3 479.3 614.5 556.7
2001 656.7 596.5 660.4 660.4 541.3 628.9 628.9 628.9 754.4 762.1
2002 747.3 633.3 543.4 543.4 498.4 617.4 617.4 617.4 749.7 777.1
2003 553.1 558.9 639.3 639.3 527.4 527.4 527.4 628.4 608.5
2004 526.9 570.3 501.0 470.6 470.6 470.6 607.9
2005 472.8 474.2 576.8 576.8 517.8 480.6 480.6 480.6 568.2 501.2
2006 623.2 616.5 605.4 605.4 985.4 590.1 590.1 590.1 612.6 677.3
2007 620.0 661.4 603.8 603.8 480.9 534.4 534.4 534.4 600.4 669.6
2008 543.6 536.2 553.0 553.0 561.2 458.3 458.3 458.3 570.3 567.6
2009 579.6 513.9 444.3 444.3 517.2 544.1 544.1 544.1 550.9 584.5
2010 473.9 657.7 405.5 405.5 626.5 505.0 505.0 505.0 528.4 550.6
2011 558.6 413.7 527.0 527.0 581.7 496.6 496.6 496.6 525.9 547.0
2012 570.5 510.0 614.7 614.7 614.4 565.3 565.3 565.3 614.4 597.2
2013 571.5 646.1 646.1 659.5 477.6 477.6 477.6 667.1 549.1
2014 463.8 645.0 645.0 494.9 423.0 423.0 423.0 556.9 477.0
2015 594.4 562.6 851.4 851.4 674.6 604.6 604.6 604.6 668.0 636.4
2016 463.5 422.7 526.0 526.0 546.4 388.0 388.0 388.0 470.9 435.8
 ANEXO 4
NUMERO DE CURVA DE LAS SUBCUENCAS
DEL RIO LA PAZ NORTE
ANEXO 4 NUMERO DE CURVA DE LAS SUBCUENCAS DEL RIO LA PAZ NORTE

Número de Curva Ponderado de la sub-cuenca Río Choqueyapu.

VALUE COUNT SUELOS COBUSO NC PONDERA COUNT*CN

1 26821 1 45 68 12,108 1823828

2 26033 2 45 79 13,653 2056607

3 10989 3 45 86 6,274 945054

4 115 3 54 91 0,069 10465

5 134 2 54 86 0,077 11524

6 151 1 54 77 0,077 11627

7 317 4 45 89 0,187 28213

8 26798 1 38 68 12,097 1822264

9 10497 2 38 79 5,505 829263

10 2556 3 38 86 1,459 219816

11 11 4 53 94 0,007 1034

12 3 3 53 91 0,002 273

13 10744 1 39 71 5,064 762824

14 8024 2 39 81 4,315 649944

15 1648 3 39 88 0,963 145024

16 39 1 8 45 0,012 1755

17 130 1 53 77 0,066 10010

18 1205 1 41 77 0,616 92785

19 64 2 41 86 0,037 5504

26 4808 2 47 84 2,681 403872

27 16049 1 47 74 7,884 1187626

28 762 3 47 90 0,455 68580

29 1557 2 58 60 0,620 93420

 30 876 1 58 36 0,209 31536

31 303 3 58 73 0,147 22119

TOTAL 150634 CN PONDERDO 74,585 74,585

Número de Curva Ponderado de la sub-cuenca río Choqueyapu 1.

VALUE COUNT SUELOS COBUSO NC PONDERA COUNT*CN

1 422 1 45 68 1,775 28696

2 226 2 45 79 1,104 17854

3 35 3 45 86 0,186 3010

8 111 1 38 68 0,467 7548

9 29 2 38 79 0,142 2291

10 4 3 38 86 0,021 344

13 4698 1 39 71 20,627 333558

14 3253 2 39 81 16,294 263493

15 937 3 39 88 5,099 82456

26 1560 2 47 84 8,103 131040

27 4776 1 47 74 21,855 353424

28 120 3 47 90 0,668 10800

TOTAL 16171 CN PONDERADO 76,341 76,341

Número de Curva Ponderado de la sub-cuenca río Orkojahuira.

VALUE COUNT SUELOS COBUSO NC PONDERA COUNT*CN

1 22053 1 45 68 15,260 1499604

2 13154 2 45 79 10,574 1039166

3 9488 3 45 86 8,303 815968

7 89 4 45 89 0,081 7921

8 17357 1 38 68 12,010 1180276

9 11879 2 38 79 9,549 938441

10 5544 3 38 86 4,852 476784

13 3464 1 39 71 2,503 245944

14 2565 2 39 81 2,114 207765

15 713 3 39 88 0,638 62744

16 65 1 8 45 0,030 2925

17 82 1 53 77 0,064 6314

18 953 1 41 77 0,747 73381

19 97 2 41 86 0,085 8342

21 8 3 41 91 0,007 728

22 6 4 38 89 0,005 534

26 2885 2 47 84 2,466 242340

27 7503 1 47 74 5,650 555222

28 150 3 47 90 0,137 13500

29 166 2 58 60 0,101 9960

30 33 1 58 36 0,012 1188

31 18 3 58 73 0,013 1314
 TOTAL 98272 CN PONDERADO 75,203 75,203

Número de Curva Ponderado de la sub-cuenca río Irpavi.

VALUE COUNT SUELOS COBUSO NC PONDERA COUNT*CN

1 33472 1 45 68 12,837 2276096

2 36666 2 45 79 16,337 2896614

3 41525 3 45 86 20,141 3571150

4 239 3 54 91 0,123 21749

5 149 2 54 86 0,072 12814

6 188 1 54 77 0,082 14476

7 1709 4 45 89 0,858 152101

8 11158 1 38 68 4,279 758744

9 14013 2 38 79 6,243 1107027

10 9959 3 38 86 4,830 856474

11 35 4 53 94 0,019 3290

12 169 3 53 91 0,087 15379

13 8586 1 39 71 3,438 609606

14 6071 2 39 81 2,773 491751

15 3024 3 39 88 1,501 266112

16 125 1 8 45 0,032 5625

17 405 1 53 77 0,176 31185

18 2678 1 41 77 1,163 206206

19 554 2 41 86 0,269 47644

 20 118 2 53 86 0,057 10148

21 214 3 41 91 0,110 19474

22 353 4 38 89 0,177 31417

23 3 2 8 66 0,001 198

24 33 4 41 94 0,017 3102

25 119 4 39 89 0,060 10591

26 1411 2 47 84 0,668 118524

27 3968 1 47 74 1,656 293632

28 365 3 47 90 0,185 32850

TOTAL 177309 CN PONDERADO 78,191 78,191

Número de Curva Ponderado de la sub-cuenca río Irpavi 1.

VALUE COUNT SUELOS COBUSO NC PONDERA COUNT*CN

1 9 1 45 68 0,745 612

13 6 1 39 71 0,518 426

26 84 2 47 84 8,584 7056

27 712 1 47 74 64,097 52688

28 11 3 47 90 1,204 990

TOTAL 822 CN PONDERDO 75,148 75,148

Número de Curva Ponderado de la sub-cuenca río Achumani.

VALUE COUNT SUELOS COBUSO NC PONDERA COUNT*CN

1 7262 1 45 68 7,397 493816

2 12693 2 45 79 15,020 1002747

3 4137 3 45 86 5,329 355782

8 11080 1 38 68 11,286 753440

9 7678 2 38 79 9,086 606562

10 2586 3 38 86 3,331 222396

13 10796 1 39 71 11,482 766516

14 3344 2 39 81 4,057 270864

15 656 3 39 88 0,865 57728

22 5 4 38 89 0,007 445

26 478 2 47 84 0,601 40152

27 5996 1 47 74 6,646 443704

28 49 3 47 90 0,066 4410

TOTAL 66760 CN PONDERADO 75,173 75,173

 ANEXO 5
CARTOGRAFIA DE FACTORES RIESGO

VELOCIDAD DE AGUA

COTA DE AGUA