Guía de Trabajo de Investigación #01
Guía de Trabajo de Investigación #01
Guía de Trabajo de Investigación #01
Facultad de Ingeniería
Informe
“Elaboración de un Estudio Hidrológico sobre la determinación de los Parámetros
Geomorfológicos y precipitación media de la microcuenca Muylo Mullucro”
Integrantes
Bustamante Delgado, Alex Abel (100%)
Delgado Irene, Kevin Omar (100%)
Figueroa Díaz, Fernando Raúl (100%)
Flores Campoverde, Jhon Anderson (100%)
Fuentes Segura, Deybi Yhoseht (100%)
Curso
Hidrología Ambiental - A
Docente
Reyes Aspiros, Roberto Bruno
Ciclo Académico
VI
ÍNDICE
I. Introducción .................................................................................................................................................. 3
II. Objetivos ................................................................................................................................................... 4
2.1. Objetivo General .................................................................................................................................. 4
2.2. Objetivos Específicos ............................................................................................................................ 4
III. Contenido Temático ................................................................................................................................. 5
3.1. Parámetros Geomorfológicos de cuencas hidrográficas .................................................................... 5
3.2. Precipitación Media de cuencas hidrográficas ................................................................................... 6
IV. Descripción general de la microcuenca Muylo ....................................................................................... 7
4.1. Ubicación ............................................................................................................................................... 7
4.1.1. Ubicación Geográfica ................................................................................................................... 7
4.1.2. Ubicación Hidrográfica ................................................................................................................ 7
4.1.3. Vías de acceso ................................................................................................................................ 7
V. Metodología de cálculo ................................................................................................................................. 8
5.1. Cálculo de los Parámetros Geomorfológicos ...................................................................................... 8
5.1.1. Delimitación de la divisoria de aguas .......................................................................................... 8
5.1.2. Cálculo del área de la divisoria de aguas .................................................................................. 28
5.1.3. Cálculo del perímetro de la divisoria de aguas ......................................................................... 29
5.2. Cálculo de la Pendiente Media........................................................................................................... 31
5.2.1. Método de Albord ....................................................................................................................... 31
5.2.2. Método de Horton ....................................................................................................................... 33
5.2.3. Método del Rectángulo Equivalente .......................................................................................... 43
5.2.4. Tiempo de concentración ........................................................................................................... 47
5.2.4.1. Tiempo de concentración por el método de Albord ............................................................. 47
5.2.4.2. Tiempo de concentración por el método de Horton ............................................................. 49
5.2.4.3. Tiempo de concentración por el método de Rectángulo Equivalente ................................. 51
5.3. Cálculo de la Precipitación Media ..................................................................................................... 54
5.3.1. Método Aritmético ...................................................................................................................... 54
5.3.2. Método Poligonal de Thiessen .................................................................................................... 55
5.3.3. Método de Isoyetas ..................................................................................................................... 82
VI. Conclusiones ......................................................................................................................................... 105
VII. Recomendaciones.................................................................................................................................. 106
VIII. Referencias Bibliográficas ................................................................................................................ 107
IX. Anexos ................................................................................................................................................... 108
“Elaboración de un Estudio Hidrológico sobre la determinación de los
Parámetros Geomorfológicos y Precipitación Media de la microcuenca
Muylo Mullucro”
I. Introducción
Un Ingeniero Civil Ambiental debe ser capaz de evaluar, analizar, interpretar y resolver todo aquello
que compete a la Hidrología, ciencia natural que estudia el agua, para que se tenga un correcto
abastecimiento de agua potable a una población o a una industria y pueda ser aprovechadas por todas
las personas. Es por eso que, en el presente estudio hidrológico, se abordará la determinación de los
Parámetros Geomorfológicos de la microcuenca Muylo Mullucro y también se determinará la
Precipitación Media de la microcuenca ya mencionada. Dicho estudio nos permitirá saber cómo drenan
o fluyen las aguas de una cuenca y cuál es el final de estas, debido a que los parámetros
geomorfológicos: área de la cuenca, perímetro de la cuenca, forma de la cuenca, relieve de la cuenca y
la red hidrográfica de la cuenca, permiten el estudio cuantitativo de las características físicas de la
microcuenca, analizando las redes de drenaje y sus pendientes. Asimismo, el aprovechamiento de agua
que se puede obtener del análisis hidrológico puede y se pretende beneficiar a las personas de la
comunidad de donde se está realizando dicho estudio, que se encuentra ubicado en el distrito de Tarma,
Provincia de Tarma, en el departamento de Junín. Dichos beneficios se pueden expresar en el regadío
continuo de las plantas, en el abastecimiento poblacional, pecuario y agrícola, entre otros beneficios.
Por último, para llevar a cabo este estudio hidrológico, se tendrá el apoyo de diferentes softwares
disponibles como: Google Earth, Global Mapper, AutoCad, Hidroesta y ArcGis; además, haremos uso
de los métodos: Albord, Horton y Rectángulo Equivalente para calcular la pendiente media de la
microcuenca y el uso de los métodos: Aritmético, Poligonal de Thiessen e Isoyetas para calcular la
precipitación media de la misma, con el fin de lograr un excelente aprovechamiento y conocimiento
hídrico de la microcuenca hidrográfica Muylo Mullucro.
II. Objetivos
❖ Área de la cuenca
El área de la cuenca se define como la superficie encerrada por la divisoria de aguas. La divisoria de
aguas es una línea imaginaria formada por los puntos de nivel topográfico más alto, marcando así, el
perímetro de la cuenca.
❖ Perímetro de la cuenca
Es la medición de la línea envolvente de la cuenca hidrográfica (Gaspari, et at., 2013). Por tanto, es la
longitud de la línea divisoria de aguas.
❖ Forma de la cuenca
Este indicador pretende parametrizar la cuenca para obtener una aproximación a su forma geométrica
(Fierro & Jiménez, 2011), los índices principales son:
b) Factor de forma
Es un factor de proporción entre el área de una cuenca rectangular (ideal) y el área de una
cuenca real, que establece la relación entre el ancho promedio de la cuenca y su longitud axial,
calculada mediante la fórmula Ff = Área / Longitud axial. La principal ventaja de este índice
cuantitativo es la de permitir obtener una idea acerca de la posibilidad de la ocurrencia de
crecidas. (Páez, et at., 2007). Por tanto, se deduce la siguiente ecuación:
𝑨
𝑭𝒇 = 𝟐
𝑳
𝐹𝑓 = 𝐹𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑓𝑜𝑟𝑚𝑎; 𝐴 = Á𝑟𝑒𝑎 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑒𝑛𝑐𝑎; 𝐿 = 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑒𝑛𝑡𝑎
𝑰𝒄 ∗ √𝑨 𝟏, 𝟏𝟐𝟖 𝟐 𝑰𝒄 ∗ √𝑨 𝟏, 𝟏𝟐𝟖 𝟐
𝑳= √
[𝟏 + 𝟏 − ( ) ] 𝒍= √
[𝟏 − 𝟏 − ( ) ]
𝟏. 𝟏𝟐𝟖 𝑰𝒄 𝟏. 𝟏𝟐𝟖 𝑰𝒄
𝐴 = Á𝑟𝑒𝑎 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑒𝑛𝑐𝑎; 𝐼𝑐 = 𝐼𝑛𝑑𝑖𝑐𝑒 𝑑𝑒 𝑐𝑜𝑚𝑝𝑎𝑐𝑖𝑜𝑑𝑎𝑑 𝑑𝑒 𝐺𝑟𝑎𝑣𝑒𝑙𝑖𝑜𝑢𝑠
❖ Relieve de la cuenca
a) Pendiente media de la cuenca
La pendiente de una cuenca tiene una relación con la infiltración, la escorrentía superficial, la
humedad del suelo y la contribución del agua subterránea a la escorrentía (Quispe, 2016).
• Criterio de Albord
𝑫∗𝑳
𝑺= ; 𝒅𝒐𝒏𝒅𝒆 𝑫: 𝑫𝒆𝒔𝒏𝒊𝒗𝒆𝒍 𝒆𝒏𝒕𝒓𝒆 𝒄𝒖𝒓𝒗𝒂𝒔, 𝑳: 𝒍𝒐𝒏𝒈𝒊𝒕𝒖𝒅 𝒕𝒐𝒕𝒂𝒍 𝒅𝒆 𝒍𝒂𝒔 𝒄𝒖𝒓𝒗𝒂𝒔 𝒚 𝑨: á𝒓𝒆𝒂 𝒅𝒆 𝒍𝒂 𝒄𝒖𝒆𝒏𝒄𝒂
𝑨
• Criterio de Horton
• Rectángulo Equivalente
𝑯
𝑺= ; 𝒅𝒐𝒏𝒅𝒆 𝑯: 𝒅𝒆𝒔𝒏𝒊𝒗𝒆𝒍 𝒕𝒐𝒕𝒂𝒍 𝒚 𝑳: 𝒍𝒂𝒅𝒐 𝒎𝒂𝒚𝒐𝒓 𝒅𝒆𝒍 𝒓𝒆𝒄𝒕á𝒏𝒈𝒖𝒍𝒐
𝑳
b) Curva Hipsométrica
Es la representación gráfica del relieve de una cuenca. Es la curva que puesta en coordenadas
rectangulares representa la relación entre la altitud y la superficie parcial de la cuenca que
queda sobre esa altitud (Quispe, 2016).
4.1. Ubicación
La microcuenca de Muylo, se encuentra en la provincia y distrito de Tarma, departamento de Junín. Su
extensión o superficie que abarca es de 23 311 Ha.
Cuadro de rutas
Medio de
De a Tipo de vía Tiempo Frecuencia
transporte
Centro
Av.
Cívico de Pavimentada 6 min Vehicular Particular
Bermudez
Tarma
Av. Jr. Malecón
Pavimentada 12 min Vehicular Particular
Bermudez Progreso
Jr. Malecón Microcuenca Trocha
30 min Camino a pie Particular
Progreso Muylo carrozable
V. Metodología de cálculo
En primer lugar, lo que realizaremos es crear una carpeta en Google Earth con el nombre de “PRODUCTO
ACREDITABLE 1 – HIDROLOGÍA”
Como segundo paso abriremos el punto de control, en este caso se llama “Punto de Control Tarma”,
seleccionamos y damos en abrir
Como podemos observar, ya el Google Earth nos ubicó el Punto de Control el cual será puesto en estudio
Ahora abriremos las curvas del nivel del lugar, para ello nosotros ya deberíamos haber descargado la
información del lugar de las Cartas Nacionales IGN. Damos Ctrl + O, y primero en la parte inferior derecha
vamos a cambiar de “*.kml*.kmz…” a ERSI Shape (*.shp), con la finalidad que nos deje abrir las curvas
Como vemos ya apareció con el nombre de “23l-curvas.shp”. Entonces damos clic en “23l-curvas.shp” y abrir
Ahora cambiaremos el color de las cotas, para nuestro proyecto lo pondremos de color ROJO
En Altura tenemos que verificar que este seleccionado “Sujetar elementos al suelo”
Podemos observar que las Curvas del nivel ya están en nuestro Google Earth
Una vez puesto las curvas, también procederemos a colocar los Rios y Lagos del lugar, para ello iremos a la carpeta Hd y
clic en 23l-rios.shp, y abrir
Al dar abrir, como en las curvas de nivel crearemos una nueva plantilla, en primer lugar, desplegamos las opciones que nos
brinda “Establecer campo de nombre” y seleccionamos el que dice “TIPO”
En el color lo pondremos celeste para que se acerque a la realidad
En altura siempre debemos darnos cuenta de que debe estar seleccionado “Sujetar elementos al suelo”
Como vemos ya se nos insertó los ríos del lugar
Ahora crearemos nuestro polígono que sería nuestra divisoria de agua entre Cuencas, para ello damos clic en polígono, y
cambios de nombre a “Divisoria de Agua”
Cambiamos el estilo, color de las líneas o perímetro de la divisora del agua y el área de igual forma, movemos a un lado el
recuadro y comenzamos a hacer nuestro polígono
Como vemos este es el proceso de la creación de nuestro polígono, tenemos que ser detallistas e ir calmados moviéndose de
diferentes lugares para que nos ubiquemos en las partes más altas de nuestra divisoria de agua
Como podemos observar logramos hacer nuestra Divisora de Agua
Una vez hecho la divisoria de agua, procedemos a guardarlo para luego trabajar en el programa Globas Mapper
Ahora ingresamos al programa GLOBAL MAPPER y abrimos el archivo de Divisoria de agua que antes guardamos. Damos
clic en “Open Data Files”
Como podemos observar ya lo proceso el programa, ahora daremos clic en ajustes para cambiar la Georreferencia en UTM
Damos clic en Projection y cambiamos a UTM
Después de ello conectamos el programa dando clic en el mundillo que hay ahí
En este apartado tenemos que darnos cuenta que este seleccionado “ASTER GDEMN…”, y damos clic en “CONNECT”
Como vemos ya fue procesado el lugar, ahora procederemos a crear las Curvas de Nivel.
Para crear las Curvas de Nivel daremos clic en “Analysis”, luego clic en “Generate Contours …”
Haciendo el paso anterior colocaremos que las Curvas de Nivel vayan en un intervalo de 10m, siendo que entre curvas
menores vayan entre 5m y las curvas mayores vayan entre 10m, damos OK
Procesamos y verificamos que las Curvas han sido puestas
Ahora abriremos los ríos dando Clic en “File”, luego “Open Date Files”
Damos Clic en “23I.ríos.shp” y abrir
Los ríos han sido procesados, es hora de agregar detalles, como cambiando de color a las curvas y a los ríos
Damos anti-clic en donde editaremos, en este caso serán los ríos, luego de ello damos clic en “optiones…”
Cambiaremos a color celeste para que se vea más real, y así de igual forma lo haremos con las curvas y la Divisoria de agua
Después de ello vamos a exportarlo en DWG para realizar el plano en Autocad, para ello damos clic en File
Abrimos AutoCad y editamos hasta que quede solo la Divisora de Agua, cambiando de capas y borrando lo innecesario
Abrimos el archivo en Google Earth, hacemos clic derecho en divisoria de agua y luego clic en
propiedades, vamos a medidas y vemos cuando de área tiene, de igual manera en AutoCad, hacemos
clic en la divisoria de aguas, nos dirigimos a propiedades y vemos cuánto de área tiene.
Obteniendo, tanto en el Google Earth como AutoCad, que la divisoria de agua tiene un área de 11.3
km2.
Abrimos el archivo en Google Earth, hacemos clic derecho en divisoria de agua y luego clic en
propiedades, vamos a medidas y vemos cuando de perímetro tiene, de igual manera en AutoCad,
hacemos clic en la divisoria de aguas, nos dirigimos a propiedades y vemos cuánto de perímetro tiene.
Obteniendo, tanto en el Google Earth como AutoCad, la divisoria de agua tiene un perímetro de 14.5
km.
5.2. Cálculo de la Pendiente Media
Comenzamos marcando estas franjas de 40m en 40m y teniendo nuestra longitud media cada 20m.
Iniciando desde nuestra curva de nivel mas baja 3110 junto con la curva 3150 formamos una zona que
sería nuestra Area N°1 y en este caso nuestra longitud media sería en la curva 3130. Este procedimieto
se repetirá hasta tener todas las demás franjas.
Estos datos lo vamos anotando en la siguiente tabla, por ejemplo:
Para A1, por programa tenemos:
✓ 𝐶𝑜𝑡𝑎 𝑀𝑒𝑛𝑜𝑟 = 3110, 𝐶𝑜𝑡𝑎 𝑀𝑒𝑑𝑖𝑎 = 3130, 𝐶𝑜𝑡𝑎 𝑀𝑎𝑦𝑜𝑟 = 3150
693.59
✓ 𝐿𝑜𝑛𝑔. 𝐶𝑜𝑡𝑎 𝑀𝑒𝑑𝑖𝑎 = 693.59𝑚 → 1000 = 0.69𝐾𝑚
✓ 𝐴𝑟𝑒𝑎 𝑑𝑒 𝐼𝑛𝑓𝑙𝑢𝑒𝑛𝑐𝑖𝑎 = 0.11 𝐾𝑚2
✓ 𝐷 = 40𝑚 → 0.04𝐾𝑚
Entonces por fórmula para la pendiente de la faja analizada:
𝐷×𝐿 0.04×0.69
✓ 𝑆𝑖 = 𝐴 , por lo tanto, 𝑆1 = 0.11 = 0.25 × 100 = 25%
AREA DE
N° DE COTA COTA COTA LONG.COTA LONG.COTA
INFLUENCIA D(km) Si (%)
AREA MEDIA MENOR MAYOR MEDIA (m) MEDIA (km)
(km2)
A1 3130 3110 3150 693.59 0.69 0.11 0.04 25.61%
• Como observamos en la imagen, una vez realizado el procedimiento, tenemos que en la línea 2 de la
malla hay 104 intersecciones en total. Y este procedimiento se realiza para las líneas restantes de la
malla de cuadrados.
• En la imagen podemos observar que ya realizamos las intersecciones con las curvas de nivel de todas las
líneas de la malla en el sentido x.
• Ahora realizamos el mismo procedimiento para las intersecciones, pero ahora lo realizaremos en el
sentido Y.
• Iniciamos con la línea 1 de la malla en el sentido Y, con la ayuda de los pequeños círculos rojos
enumerados, lo ubicamos en cada intersección de la línea de la malla con las curvas de nivel.
• Como podemos observar en la imagen, en la línea 1 de la malla en el sentido Y, hay un total de 84
intersecciones. Realizamos el mismo procedimiento para las demás líneas de la malla en el sentido Y.
• Como podemos observar en la imagen, ya hemos realizado todas las intersecciones de las lineas de la
malla con las curvas de nivel, tanto en el eje X como en el eje Y.
• Ahora procedemos a medir la longitud de cada línea de la malla, el cual se medirá desde el inicio de las
intersecciones hasta el fin. En la imagen observamos la medición de la longitud de la línea 1 de la malla
en el eje X.
• Podemos observar en la imagen que la longitud de la línea 1 en el eje x es de 2922.6642m. Realizamos
este proceso para las líneas restantes en el eje x.
• En esta imagen podemos observar la longitud de 2328.5304m, que pertence a la línea 1 de la malla en el
eje Y. Realizamos este mismo proceso para las líneas restantes en el eje Y.
𝐼𝑐√𝐴 1.12 2
𝐿= √
[1 + 1 − ( ) ]
1.12 𝐼𝑐
1.2164√11.3 1.12 2
𝐿= √
[1 + 1 − ( ) ]
1.12 1.2164
𝐿 = 5.0753
Lado l:
𝐼𝑐√𝐴 1.12 2
𝑙= √
[1 − 1 − ( ) ]
1.12 𝐼𝑐
1.2164√11.3 1.12 2
𝑙= √
[1 − 1 − ( ) ]
1.12 1.2164
𝑙 = 2.2265
Para las características topográficas de la microcuenca. Tomamos las curvas de nivel cada 40 metros y
el área en km2.
Área de la curva de nivel 3110 – 3150:
Después de tener todas las áreas, calculamos los segmentos del lado mayor Li (dividiendo cada área
parcial, entre el lado menor l del rectángulo equivalente):
Para el 𝑨𝟏 :
𝐴1 0.108
𝐿𝑖1 = => 𝐿𝑖1 = => 𝑳𝒊𝟏 = 𝟎. 𝟎𝟒𝟗
𝑙 2.2265
Para el 𝑨𝟐 :
𝐴2 0.129
𝐿𝑖2 = => 𝐿𝑖2 = => 𝑳𝒊𝟐 = 𝟎. 𝟎𝟓𝟖
𝑙 2.2265
Así sucesivamente calculamos todos los segmentos del lado mayor Li y lo ubicamos en la tabla:
A (km2) Li (km)
𝐴1 = 0.108 0.049
𝐴2 = 0.129 0.058
𝐴3 = 0.138 0.062
𝐴4 = 0.166 0.074
𝐴5 = 0.237 0.106
𝐴6 = 0.292 0.131
𝐴7 = 0.333 0.150
𝐴8 = 0.340 0.153
𝐴9 = 0.319 0.143
𝐴10 = 0.323 0.145
𝐴11 = 0.331 0.149
𝐴12 = 0.351 0.158
𝐴13 = 0.404 0.181
𝐴14 = 0.437 0.196
𝐴15 = 0.561 0.252
𝐴16 = 0.624 0.280
𝐴17 = 0.638 0.286
𝐴18 = 0.697 0.313
𝐴19 = 0.729 0.327
𝐴20 = 0.727 0.326
𝐴21 = 0.726 0.326
𝐴22 = 0.705 0.317
𝐴23 = 0.654 0.294
𝐴24 = 0.616 0.277
𝐴25 = 0.383 0.172
𝐴26 = 0.208 0.093
𝐴27 = 0.109 0.049
𝐴28 = 0.027 0.012
Donde: H = desnivel total de la microcuenca (m) y L = lado mayor del rectángulo equivalente (m)
𝐾𝑐 𝑥 √𝐴 1.128 2
𝐿= √
[1 + 1 − ( ) ]
1.128 𝐾𝑐
𝐿 = 5.08 𝐾𝑚
Longitud del lado menor del rectángulo
𝐾𝑐 𝑥 √𝐴 1.128 2
𝑙= √
[1 − 1 − ( ) ]
1.128 𝐾𝑐
𝐿 = 2.23 𝐾𝑚
A.08 ALTITUD MAS BAJA DE LA CUENCA
𝐶. 𝑚 = 3 110.00 𝑚. 𝑠. 𝑛. 𝑚
A.09 ALTITUD MAS ALTA DE LA CUENCA
𝐶. 𝑀 = 4 390.00 𝑚. 𝑠. 𝑛. 𝑚
A.10 DESNIVEL DE LA CUENCA
𝐻. 𝑐 = 4 390.00 𝑚. 𝑠. 𝑛. 𝑚 − 3 110.00 𝑚. 𝑠. 𝑛. 𝑚
𝐻. 𝐶 = 1 280.00 𝑚
A.11 PENDIENTE DE LA CUENCA. APLICANDO METODO ALBORD (m/m)
𝐷
𝑆= (𝑙 + 𝑙2 + 𝑙3 + ⋯ + 𝑙𝑛 )
𝐴 1
Donde:
✓ S = Pendiente de la cuenca
✓ D = Desnivel entre las líneas medias, aceptado como desnivel entre curvas (equidistancia)
✓ l = Longitud de la curva de nivel correspondiente a la faja analizada i.
✓ A = Área total de la cuenca
0.04
𝑆= (141.156)
11.312
𝑆 = 0.4991 𝑚/𝑚
A.12 PENDIENTE DE LA CUENCA (%)
𝑆 = 0.4991 𝑥 100% = 49.91 %
A.13 TIEMPO DE CONCENTRACION (METODO DE KIRPICH)
2 0.385
(𝐿𝑝 )
KIRPICH 𝑇𝑐1 = 0.06626 ( )
𝑆
Donde:
✓ Tc= Tiempo de concentración (hr)
✓ (Lp o L) = Longitud del curso o cauce principal (Km)
✓ (Sp o S)= Pendiente de media de la cuenca y/o microcuenca (m/m)
0.385
(3.81)2
𝑇𝑐1 = 0.06626 ( )
0.4991
𝑇𝑐1 = 0.24 𝐻𝑟
A.14 TIEMPO DE CONCENTRACION (METODO DE TEMEZ)
0.75
𝐿𝑝
TEMEZ 𝑇𝑐2 = 0.126 ( 0.35 )
(𝑆𝑝 )
0.75
3.81
𝑇𝑐2 = 0.126 ( )
(0.4991)0.35
𝑇𝑐2 = 0.41 𝐻𝑟
A.15 TIEMPO DE CONCENTRACION (U.S.CORP.ING)
𝐿 0.76
U.S.CORP.ING 𝑇𝑐3 = 0.3 ((𝑆)0.25)
3.81 0.76
𝑇𝑐3 = 0.3 ((0.4991)0.25)
𝑇𝑐3 = 0.95 𝐻𝑟
A.16 TIEMPO DE CONCENTRACION (PROMEDIO)
𝑇𝑐1 + 𝑇𝑐2 + 𝑇𝑐3
𝑇𝑐 =
3
0.24 + 0.41 + 0.95
𝑇𝑐 = = 0.53 𝐻𝑟
3
𝐿 = 5.08 𝐾𝑚
Longitud del lado menor del rectángulo
𝐾𝑐 𝑥 √𝐴 1.128 2
𝑙= [1 − √1 − ( ) ]
1.128 𝐾𝑐
𝐿 = 2.23 𝐾𝑚
A.08 ALTITUD MAS BAJA DE LA CUENCA
𝐶. 𝑚 = 3 110.00 𝑚. 𝑠. 𝑛. 𝑚
A.09 ALTITUD MAS ALTA DE LA CUENCA
𝐶. 𝑀 = 4 390.00 𝑚. 𝑠. 𝑛. 𝑚
A.10 DESNIVEL DE LA CUENCA
𝐻. 𝑐 = 4 230.00 𝑚. 𝑠. 𝑛. 𝑚 − 3 110.00 𝑚. 𝑠. 𝑛. 𝑚
𝐻. 𝐶 = 1 280.00 𝑚
A.11 PENDIENTE DE LA CUENCA. APLICANDO METODO HORTON (m/m)
(𝑆𝑥 + 𝑆𝑦 )
𝑆=
2
Donde:
✓ S = Pendiente de la cuenca
✓ Sx = Pendiente en el sentido x
✓ Sy = Pendiente en el sentido y
0.3043 + 0.3326
𝑆=
2
𝑆 = 0.3184 𝑚/𝑚
A.12 PENDIENTE DE LA CUENCA (%)
𝑆 = 0.3184 𝑥 100% = 31.84 %
A.13 TIEMPO DE CONCENTRACION (METODO DE KIRPICH)
2 0.385
(𝐿𝑝 )
KIRPICH 𝑇𝑐1 = 0.06626 ( )
𝑆
Donde:
✓ Tc= Tiempo de concentración (hr)
✓ (Lp o L) = Longitud del curso o cauce principal (Km)
✓ (Sp o S) = Pendiente de media de la cuenca y/o microcuenca (m/m)
0.385
(3.81)2
𝑇𝑐1 = 0.06626 ( )
0.3184
𝑇𝑐1 = 0.29 𝐻𝑟
A.14 TIEMPO DE CONCENTRACION (METODO DE TEMEZ)
0.75
𝐿𝑝
TEMEZ 𝑇𝑐2 = 0.126 ( 0.35 )
(𝑆𝑝 )
0.75
3.81
𝑇𝑐2 = 0.126 ( )
(0.3184)0.35
𝑇𝑐2 = 0.46 𝐻𝑟
A.15 TIEMPO DE CONCENTRACION (U.S.CORP.ING)
𝐿 0.76
U.S.CORP.ING 𝑇𝑐3 = 0.3 ((𝑆)0.25)
3.81 0.76
𝑇𝑐3 = 0.3 ((0.3184)0.25)
𝑇𝑐3 = 1.03 𝐻𝑟
A.16 TIEMPO DE CONCENTRACION (PROMEDIO)
𝑇𝑐1 + 𝑇𝑐2 + 𝑇𝑐3
𝑇𝑐 =
3
0.29 + 0.46 + 1.03
𝑇𝑐 = = 0.59 𝐻𝑟
3
𝐾𝑐 𝑥 √𝐴 1.128 2
𝐿= [1 + √1 − ( ) ]
1.128 𝐾𝑐
𝐿 = 5.08 𝐾𝑚
Longitud del lado menor del rectángulo
𝐾𝑐 𝑥 √𝐴 1.128 2
𝑙= [1 − √1 − ( ) ]
1.128 𝐾𝑐
𝐿 = 2.23 𝐾𝑚
A.08 ALTITUD MAS BAJA DE LA CUENCA
𝐶. 𝑚 = 3 110.00 𝑚. 𝑠. 𝑛. 𝑚
A.09 ALTITUD MAS ALTA DE LA CUENCA
𝐶. 𝑀 = 4 230.00 𝑚. 𝑠. 𝑛. 𝑚
A.10 DESNIVEL DE LA CUENCA
𝐻. 𝑐 = 4 230.00 𝑚. 𝑠. 𝑛. 𝑚 − 3 110.00 𝑚. 𝑠. 𝑛. 𝑚
𝐻. 𝐶 = 1 120.00 𝑚
A.11 PENDIENTE DE LA CUENCA. APLICANDO METODO RECTANGULO EQUIVALENTE (m/m)
𝐻
𝑆=
𝐿
Donde:
✓ H = Desnivel total de la microcuenca (m)
✓ L = Lado mayor del rectángulo equivalente (m)
4230 − 3110
𝑆=
5.08 𝑥 1000
𝑆 = 0.2204 𝑚/𝑚
A.12 PENDIENTE DE LA CUENCA (%)
𝑆 = 0.2204 𝑥 100% = 22.04 %
A.13 TIEMPO DE CONCENTRACION (METODO DE KIRPICH)
2 0.385
(𝐿𝑝 )
KIRPICH 𝑇𝑐1 = 0.06626 ( )
𝑆
Donde:
✓ Tc= Tiempo de concentración (hr)
✓ (Lp o L) = Longitud del curso o cauce principal (Km)
✓ (Sp o S)= Pendiente de media de la cuenca y/o microcuenca (m/m)
0.385
(3.81)2
𝑇𝑐1 = 0.06626 ( )
0.2204
𝑇𝑐1 = 0.33 𝐻𝑟
A.14 TIEMPO DE CONCENTRACION (METODO DE TEMEZ)
0.75
𝐿𝑝
TEMEZ 𝑇𝑐2 = 0.126 ( 0.35 )
(𝑆𝑝 )
0.75
3.81
𝑇𝑐2 = 0.126 ( )
(0.2204)0.35
𝑇𝑐2 = 0.51 𝐻𝑟
3.81 0.76
𝑇𝑐3 = 0.3 ((0.2204)0.25)
𝑇𝑐3 = 1.10 𝐻𝑟
A.16 TIEMPO DE CONCENTRACION (PROMEDIO)
𝑇𝑐1 + 𝑇𝑐2 + 𝑇𝑐3
𝑇𝑐 =
3
0.33 + 0.51 + 1.10
𝑇𝑐 = = 0.65 𝐻𝑟
3
Se cálcula la precipitación media como la media aritmética de las precipitaciones. Para ello
utilizaremos las precipitaciones encontradas en ArcGis.
Para ello haremos un promedio o media aritmética de nuestras precipitaciones la respuesta será dada en
Milímetros (mm).
300 + 400 + 200 + 850 + 750 + 620 + 980 + 610 + 290 + 750 + 120
𝑃𝑟𝑒𝑐𝑖𝑝𝑖𝑡𝑎𝑐𝑖ó𝑛 =
11
𝑃𝑟𝑒𝑐𝑖𝑝𝑖𝑡𝑎𝑐𝑖ó𝑛 = 533.64 𝑚𝑚
Como podemos observar por el programa Hidroesta nos arrojó el mismo resultado de precipitación por el
método aritmético
En primer lugar, exportamos las estaciones de nuestro trabajo, estas se exportarán en el programa Google Earth
En segundo lugar, guardamos nuestra Divisoria de agua o el área de nuestra Microcuenca en la Carpeta Gis, ahí
es donde estaremos trabajando todo el proceso
Abrimos el programa Global Mapper y abrimos el archivo guardado del área de la Microcuenca, esta lleva por
nombre “Area_quebrada_Tarma”
Una vez abierto damos clic en la llave de configuración, para así configurar nuestro proyecto en coordenadas
UTM
Guardamos el archivo presionando Ctrl+S, este archivo se guardará en la misma carpeta donde estamos
trabajando “GIS”
Exportamos en versión “Shapefile”, damos clic en OK
Nos saldrá esta pestaña, damos clic en “Select File”, para ver la carpeta donde lo guardaremos.
Lo guardaremos en la carpeta Gis, pero en el apartado de la carpeta Archivos SHP, tenemos que asegurarnos que
este en archivos ERSI SHAPE FILES
Ahora acomodamos las coordenadas y la precipitación donde corresponde y damos clic en “Edit”
Configuramos la ubicación del proyecto, dando Clic en “WGS 1984 UTM ZONE 18S”
Ahora vamos a exportar esta Data, dando Anticlic en “Estacione$ Events”, luego “Data”, y luego “Export Data”
Nuevamente lo guardaremos en la carpeta “Archivos Shp”
Clic en el símbolo de carpeta para que nos salga la ventana, seguido un clic en “Estaciones.shp” y por último en “Add”
Clic en la segunda carpeta y guardamos nuestro archivo “THIESSEN” en la carpeta “Shp”
Clic en “Environments…”
Clic en “Processing Extent” y buscamos “Same as layer Area_de_estudio” luego procesamos
Clic en “Imput_FID” y OK
Logramos procesar nuestras Estaciones con nuestro Plano
Para poder aplicar este método es necesario conocer la localización de las estaciones dentro de la zona
en estudio, ya que para su aplicación se requiere determinar la zona de influencia en cada una de ellas.
Es síntesis el método se basa en asignar cada punto de la cuenca a la estación más próxima, se deben
unir las estaciones de dos en dos y dibujar las mediatrices de estos segmentos, asignando a cada
estación el área limitada por las poligonal es que forman las mediatrices. Este proceso como vemos ya
lo hemos desarrollado utilizando ArcGIS donde nos arrojó los siguientes resultados:
Met. Thiessen
NOMBRE ESTE__X NORTE__Y ALTURA__Z PP_anual (AREA
Ai)km2
Estación 01 425793 8737790 3080 300 0.807
Estación 02 427168 8735294 3624 400 0.685
Estación 03 428120 8736563 3611 200 3.044
Estación 04 427948 8737907 3581 850 2.436
Estación 05 428968 8738743 4033 750 1.439
Estación 06 429567 8735697 4063 620 1.671
Estación 07 426593 8738616 3440 980 0.442
Estación 08 428845 8734473 3840 610 0.236
Estación 09 427559 8739809 3742 290 0.156
Estación 10 429178 8740111 4080 750 0.010
Estación 11 429846 8738534 4147 120 0.388
Después de ello aplicaremos la siguiente fórmula:
∑𝑛𝑖=1 ℎ𝑝𝑖 𝐴𝑖
ℎ𝑝𝑚 =
𝐴
𝐴 = Á𝑟𝑒𝑎 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑧𝑜𝑛𝑎 (𝑘𝑚2 )
𝐴𝑖 = Á𝑟𝑒𝑎 𝑡𝑜𝑡𝑎𝑙 𝑡𝑟𝑖𝑏𝑢𝑡𝑎𝑟𝑖𝑎 𝑑𝑒 𝑙𝑎 𝑒𝑠𝑡𝑎𝑐𝑖ó𝑛 𝑖 (𝑘𝑚2 )
ℎ𝑝𝑖 = 𝐴𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑝𝑟𝑒𝑐𝑖𝑝𝑖𝑡𝑎𝑐𝑖ó𝑛 𝑟𝑒𝑔𝑖𝑠𝑡𝑟𝑎𝑑𝑎 𝑒𝑛 𝑙𝑎 𝑒𝑠𝑡𝑎𝑐𝑖ó𝑛 𝑖 (𝑚𝑚)
ℎ𝑝𝑚 = 𝐴𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑝𝑟𝑒𝑐𝑖𝑝𝑖𝑡𝑎𝑐𝑖ó𝑛 𝑚𝑒𝑑𝑖𝑎 𝑒𝑛 𝑙𝑎 𝑧𝑜𝑛𝑎 𝑒𝑛 𝑒𝑠𝑡𝑢𝑑𝑖𝑜 (𝑚𝑚)
𝑁 = 𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑒𝑠𝑡𝑎𝑐𝑖𝑜𝑛𝑒𝑠 𝑙𝑜𝑐𝑎𝑙𝑖𝑧𝑎𝑑𝑎𝑠 𝑑𝑒𝑛𝑡𝑟𝑜 𝑑𝑒 𝑙𝑎 𝑧𝑜𝑛𝑎 𝑒𝑛 𝑒𝑠𝑡𝑢𝑑𝑖𝑜.
Entonces procederemos a Multiplicar cada área por cada precipitación, hacemos una sumatoria y
dividimos entre el área total de la microcuenca.
Met. Thiessen
PP_anual (AREA Ai * Ppi
Ai)km2
300 0.807 242.236
400 0.685 274.029
200 3.044 608.813
850 2.436 2070.290
750 1.439 1079.001
620 1.671 1036.172
980 0.442 432.863
610 0.236 144.151
290 0.156 45.234
750 0.010 7.670
120 0.388 46.607
En primer lugar, lo que haremos es abrir el ultimo archivo realizado por Thiessem. Y asimismo ocultamos el
procesamiento Inicial de Thiessen.
Como segundo paso procedemos a arrastrar THISSEN QUEBRADA02 por encima de estaciones, y
desactivamos ese mismo Layer. Procedemos a guardar un primer respaldo. Nos vamos a File-Save as…
Empezamos a interpolar los puntos para Isoyetas. Para lo cual activamos el ArcToolbox
Nos dirigimos a Spatial Analyst Tools - Interpolación – Spline (Hacemos doble click)
Desplegamos el apartado de “Input point features”, nos aparecerá una carpeta llamada “Estaciones” como lo
hayamos guardado anteriormente. Abrimos la carpeta y esperamos un momento.
Desplegamos el “Z value field” y colocamos precipitación Anual “PP_Anual”
Luego nos dirigimos al “Output raster”, nos vamos a la carpeta “Archivos Raster_isoyetas” y le colocamos de
nombre “INTERP” y le damos a Save (Guardar).
Modificamos el nombre de “INTERPOLACIÓN02” a “INTERP”, luego en “Output cell size” colocamos el
valor de 30, y en “Spline type” ponemos a “TENSION”
Nos dirigimos a “Evironment Settings” luego a “Processing Extent”, le damos a “Area de Estudio” y damos OK.
En la tabla de contenidos podemos ver los Rangos de Interpolacion
Luego de ello guardamos nuestro respaldo con el nombre de “02. Interpolacion de puntos para Isoyetas.mxd”
Procedemos con la Reclasificación de las Isoyetas. Nos dirigimos al “ArcToolbox” luego a “Spatial Analystic
Tools”, a continuación a “Reclass” y finalmente a “Reclassify” (Doble click), en “Input Raster” buscamos el
“INTERP” , para luego ir a “Classify”
Otra vez nos dirigimos al “ArcToolbox” luego a “Conversion Tools”, después a “From Raster” y finalmente a
“Raster to Polygon” (Damos doble click). Luego en “Input Raster” desplegamos y seleccionamos el “reclasif”,
en “Field” lo dejamos en “VALUE” , en el apartado de “Output polygon features” buscamos la carpeta de
nuestros “archivos Shp”, Guardamos como “RECLASS_POLYGON02” y save (Guardar) y damos OK.
Obtuvimos esto
Nos dirigimos al “ArcToolsBox”, luego a “Spatial Anslyst Tools”, siguiente a “Extraction”, a continuación a
“Extract by Mask”, en “Input raster” nos dirigimos a “reclasif”, en el apartado de “Input raster or feature mask
data” nos dirigimos a la carpeta de “Area_quebrada_tarma” y damos en “Output Raster” colocamos en nuestros
Archivos Raster “Reclass_quebradatar”, damos en Save y Ok.
Colocamos “.tif” en “Output Raster”
Obtuvimos esto.
Desactivo “reclasif” , nos dirigimos a “Geoprocessing” luego a “Clip”, en el “Input features” vamos a
“RECLASS_POLYGON02”, en “Clip features” vamos a “Area_quebrada_tarma”, luego vamos a “Output
Feature Class” y nos dirigimos a nuestros archivos Shp y luego creamos un “quebrada_clip02” damos en Save y
Ok.
En “quebrada:clip02” le damos click derecho “propierties” nos dirigimos a “Symbology” luego a “Categories”
enseguida a “Unique values”, vamos a “Value field” y colocamos “gridcode” luego le damos a “Add All
Values”, luego le damos a aceptar.
Distinguimos las zonas de precipitaciones.
Podemos cambiar el color dando click sobre el color, y seleccionando el que mas convenga para su distinción.
Continuamos nos dirigimos al “ArcToolbox” luego a “Spatial Analyst Tools” enseguida a “Zonal”, colocamos
“Zonal Statistics as Table” (Damos doble click). En el apartado de “Input raster or feature sone data” le damos a
“Reclass_quebradatar.tif”, en “Zone Field” ponemos “VALUE” , en el “Input Value raster” colocamos
“INTERP”, enseguida en “Output Table” nos vamos a la carpeta de Archivos Tabla y creamos “Zonal03.dbf”
y damos en Ok.
Por ultimo nos vamos a “quebrada_clip02” damos click derecho y le damos a “data” luego a “Export to CAD”,
luego en “Input Features” buscamos Gis base tarma “Archivos shp” seleccionamos “área_quebrada_tarma02”,
“estaciones” y también “quebrada_clip02”, en “Output File” guardamos en “archivos CAD” Le colocamos un
nombre “Procesamiento inicial en clase_Isoyetas” y damos en OK.
Hacemos el ultimo respaldo
Lo que nos trata de brindar el método de Isoyetas es básicamente las curvas que unen los puntos de
iguales precipitaciones. Lo importante de este método es que es más exacto, esto va ligado a su vez del
número de estaciones que se encuentren presentes sobre el área que se desea estudiar.
Cuando procesamos todos lo datos en el ArcGis nos va a dar los siguientes datos, que completamos con
las áreas que nos brinda el CAD, en el dado caso que un color se repita, lo que procedemos a hacer es a
sumar dichas áreas. Obtuvimos el siguiente recuadro completado.
AREA
VALUE MIN MAX PROMEDIO AREA (m2) (km2) AREA*PROMEDIO
3 190.028183 295.094574 242.561379 1163458.06 1.163 282.210
4 295.209625 411.081818 353.145722 2756139.28 2.756 973.319
5 411.155945 527.043152 469.099549 2184956.45 2.185 1024.962
6 527.145203 642.978333 585.061768 1494260.47 1.494 874.235
7 642.983154 758.842224 700.912689 1361763.77 1.362 954.478
8 758.920349 874.818542 816.869446 1521013.84 1.521 1242.470
9 874.893311 970.974548 922.93393 833122.85 0.833 768.917
TOTAL 11314714.72 11.315 6120.590
PRECIPITACION PROMEDIO 540.94
Como observamos en la Tabla nos arroja una precipitación promedio de 540.94 y para estar seguros de
que esta precipitación es correcta la comprobamos con el HIDROESTA.
Para ello solo vamos a necesitar las precipitaciones mínimas y de las máximas la mayor, que luego
ubicaremos en una sola columna. Y otro item necesario es el Area, ya sea en m2 o km2, como mejor lo
podamos trabajar. De ello generamos un Excel y lo guardamos.
Abrimos el Hidroesta 2 y nos dirigimos a “Precipitación” luego a “Isoyetas”
Una vez aquí cargamos el Excel que previamente guardamos y nos aparecerá lo siguiente. Y
verificamos cuanto nos sale la Precipitación promedio. Si es la misma que nos arroja el Excel entonces
el método está correcto.
Instituto Nacional de Ecología. 2004. Análisis morfométricos de Cuencas: Caso de estudio del Parque
Nacional Pico de Tancítaro. Dirección General de Investigación de Ordenamiento Ecológico y
conservación de ecosistemas. México. http://agua.org.mx/wp-
content/uploads/2009/02/morfometria_pico_tancitaro.pdf
Fierro, D. & Jiménez, L. (2011). Caracterización de la Microcuenca del río Manzano, Cantón Alausí,
Provincia de Chimborazo y Propuesta de Plan de Manejo, Utilizando Herramientas SIG.
Escuela Politécnica del Ejército. Sangolqui, Ecuador.
http://repositorio.espe.edu.ec/handle/21000/4288
Gaspari, F., Rodriguez, A., Senisterra, G., Denegri, G., Delgado, I., & Besteiro, S. (2012).
"Caracterización morfométrica de la cuenca alta del río Sauce Grande, Buenos Aires,
Argentina". 7mo Congreso de Medio Ambiente, 1-25.
http://sedici.unlp.edu.ar/handle/10915/25777
Páez, Gustavo & Ayala, Rubén & Araque, Fany (2007). Análisis geomorfológico de la microcuenca El
Guayabal, a propósito de las ocurrencias de las lluvias excepcionales de febrero de 2005.
Cuenca del río Mocotíes, estado Mérida-Venezuela. Revista Geográfica Venezolana, 48
(1),59-82. ISSN: 1012-1617. Disponible en:
https://www.redalyc.org/articulo.oa?id=347730365004
Villón, M. 2002. Hidrología. Segunda edición. Editorial Villón. Lima, PE. 431 p.
https://books.google.com.pe/books?id=-
JjGDwAAQBAJ&printsec=frontcover&dq=Libro+de+hidrolog%C3%ADa+M%C3%A1ximo
+Vill%C3%B3n&hl=es&sa=X&ved=2ahUKEwjDzIOa9YX3AhVlppUCHdAhD3MQ6wF6B
AgGEAE#v=onepage&q=Libro%20de%20hidrolog%C3%ADa%20M%C3%A1ximo%20Vill
%C3%B3n&f=false
IX. Anexos