Modelamiento Hidrológico con HEC-HMS
Modelamiento Hidrológico con HEC-HMS
MODELAMIENTO HIDROLOGICO CON EL PROGRAMA HEC‐HMS 4.2
EN EL PUNTO DE CONFLUENCIA
El modelamiento hidrológico del HEC HMS 4.2, determina la escorrentía causada por las lluvias de
diseño. Este programa fue creado por el Cuerpo de Ingenieros del Ejército de los Estados Unidos,
para simular el tránsito de avenidas causadas por la precipitación.
Al programa se ingresan tres datos de entrada:
Modelo de cuenca o base: modelo de pérdida de agua, modelo de transformación,
modelo de flujo de base.
Modelo meteorológico
Especificaciones de control.
MODELO DE LA CUENCA
La siguiente figura muestra el esquema de cálculo que se utilizara en la cuenca de estudio, el cual
consiste en 10 subcuencas.
Fig. N° 01 – Esquema de cálculo en el programa HEC‐HMS.
_ _
1
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
El siguiente cuadro se muestra los parámetros geomorfológicos de la cuenca estudiada.
Cuadro N° 1 – Parámetros geomorfológicos de la cuenca.
Para calcular el tiempo de concentración se utilizaron las siguientes fórmulas:
‐ Formula de Kirpich
.
0.000325 .
Dónde:
tc: Tiempo de concentración (hr).
L: Longitud del cauce principal (m).
S: Pendiente del cauce principal (m/m).
‐ Formula de U.S. Corps of Engineers
.
0.3 .
Dónde:
tc: Tiempo de concentración (hr).
L: Longitud del cauce principal (Km).
S: Pendiente del cauce principal (m/m).
El siguiente cuadro muestra el tiempo de concentración calculado:
Cuadro N° 2 – Tiempo de concentración
_ _
2
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
Kirpich U. S. C. E. Tc Tc
SUBCUENCA
(horas) (horas) (horas) (min)
04 4.17 9.42 9.42 565
05 5.05 10.82 10.82 649
06 6.95 13.73 13.73 824
07 4.48 10.34 10.34 621
08 5.01 11.11 11.11 667
09 2.88 6.78 6.78 407
10 2.60 5.98 5.98 359
Se consideró usar el tiempo de concentración obtenida por la fórmula de la U.S. Corps of
Engineers, debido a que la fórmula de Kirpich se admite más en cuencas urbanas, mientras que en
las cuencas rurales o naturales supone una notable infravaloración del concepto de tiempo de
concentración y más parece adecuarse al tiempo de retardo que separa los centros de gravedad
del hietograma y del hidrograma superficial, cuya relación con el tiempo de concentración es del
orden del 45%. En nuestros resultados si se comparan los resultados de Kirpich como tiempo de
retardo y de U.S.C.E. como tiempo de concentración se obtiene la relación mostrada en el
siguiente cuadro.
Cuadro N° 3 –Comparación Tc Kirpich y U.S.C.E.
Kirpich ‐ Tr U. S. C. E. – Tc
SUBCUENCA Tr/Tc
(horas) (horas)
01 4.16 9.53 0.44
02 6.00 12.40 0.48
03 3.29 7.59 0.43
04 4.17 9.42 0.44
05 5.05 10.82 0.47
06 6.95 13.73 0.51
07 4.48 10.34 0.43
08 5.01 11.11 0.45
09 2.88 6.78 0.42
10 2.60 5.98 0.43
Del cuadro anterior se tiene una relación promedio de 0.45, valor similar al mencionado. Además,
la fórmula de U.S.C.E. es aplicable en cuencas rurales, con lo cual se tomará como tiempo de
concentración los obtenidos por la fórmula de la U.S.C.E.
Modelo de pérdida de agua
Para el Modelo de pérdida de agua se utilizara el modelo del Número de Curva SCS, en el cual se
debe conocer 3 tipos de datos:
El número de curva CN
La abstracción inicial (mm)
Porcentaje de suelo impermeable en la cuenca (%)
_ _
3
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
El número de curva CN, varía entre 0 y 100 y se aplica para condiciones antecedentes de humedad
normales (AMC II), para condiciones seca (AMC I) o condiciones húmedas (AMC III) los números
de curva equivalentes pueden calcularse por:
4.2 23
;
10 0.058 10 0.13
El grupo de antecedente de humedad (AMC) se determina según la siguiente tabla:
Tabla N° 1 – Valores de número
GRUPO LLUVIA ANTECEDENTE TOTAL DE 5 DIAS (PULG)
AMC ESTACION INACTIVA ESTACION ACTIVA
I Menor que 0.5 Menor que 1.4
II 0.5 a 1.1 1.4 a 2.1
III Mayor que 1.1 Mayor que 2.1
Los números de curva han sido tabulados por el Soil Conservation Service con base en el tipo de
suelo y el uso de la tierra. Se definen cuatro grupos de suelos:
Grupo A: Arena profunda, suelos profundos depositados por el viento, limos agregados.
Grupo B: Suelos pocos profundos depositados por el viento, marga arenosa.
Grupo C: Margas arcillosas, margas arenosas poco profundas, suelos con bajo contenido
orgánico y suelos con altos contenidos de arcilla.
Grupo D: Suelos que se expanden significativamente cuando se mojan, arcillas altamente
pláticas y ciertos suelos salinos.
Los valores de CN para varios tipos de uso de la tierra en estos tipos de suelos se muestran en la
siguiente tabla:
Tabla N° 2 – Valores de número de curva CN.
GRUPO HIDROLÓGICO DEL SUELO
DESCRIPCIÓN DEL USO DE LA TIERRA
A B C D
Tierra cultivada:
Sin tratamiento de conservación 72 81 88 91
Con tratamiento de conservación 62 71 78 81
Pastizales:
Condiciones pobres 68 79 86 89
Condiciones óptimas 39 61 74 80
Vegas de ríos: condiciones óptimas 30 58 71 78
Bosques: 45 66 77 83
Troncos delgados, cubierta pobre, sin hierbas, cubierta buena. 25 55 70 77
Área abiertas, césped, parques, campos de golf, cementerios, etc.
Óptimas condiciones: cubierta de pasto en el 75% o más 39 61 74 80
condiciones aceptables cubierta de pasto en el 50 al 75% 49 69 79 84
Áreas comerciales de negocios (85% impermeables) 89 92 94 95
Distritos Industriales (72% impermeables) 81 88 91 93
_ _
4
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
GRUPO HIDROLÓGICO DEL SUELO
DESCRIPCIÓN DEL USO DE LA TIERRA
A B C D
Residencial:
Tamaño promedio del lote Porcentaje promedio impermeable
1/8 acre o menos 65 77 85 90 92
1/4 acre 38 61 75 83 87
1/3 acre 30 57 72 81 86
1/2 acre 25 54 70 80 85
1 acre 20 51 68 79 84
Parqueadores pavimentados, techos, accesos, etc. 98 98 98 98
Calles y carreteras:
Pavimentados con cunetas y alcantarillados. 98 98 98 98
Grava 76 85 89 91
Tierra 72 82 87 89
La abstracción inicial es la cantidad de lluvia en mm que se pierde antes de ocasionar el flujo
superficial, el cual se puede calcular según la siguiente relación:
5080
50.8
Dónde:
Ia : Abstracción inicial (mm)
CN: Número de curva
El porcentaje impermeable está dado por el área ocupada por lagunas y zonas rocosas que para
propósitos prácticos puede considerarse impermeable. Dicho porcentaje se determina según la
siguiente relación:
Á
% ∗ 100
Á
Para nuestro caso se cuenta con las siguientes lagunas y el siguiente porcentaje impermeable:
Cuadro N° 4 – Lagunas presentes en la cuenca de estudio y % impermeable.
ÁREA IMPERMEABLE ÁREA CUENCA %
SUBCUENCA LAGUNA
(Km2) (Km2) IMPERMEABLE
02 Sacacanicocha, Isquina Sacacani 4.24 872.52 0.49
03 S/N 0.86 503.36 0.17
04 Cochachaca, Cocha Uma, Sibinacocha 28.12 309.56 9.08
05 S/N 1.92 342.86 0.56
06 Asnacocha, Pampamarca 11.04 641.87 1.72
07 Coñocota 0.90 485.23 0.19
08 Langui Layo 54.11 326.86 16.55
_ _
5
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
Modelo de transformación
El modelo de transformación será el modelo del Hidrograma sintético triangular SCS, el cual
requiere como dato de ingreso el tiempo de retardo. El tiempo de retardo se puede relacionar de
la siguiente manera:
0.60 ∗
Dónde:
tr : Tiempo de retardo (mm)
tc : Tiempo de concentración (mm)
Modelo de flujo base
No se considerara el modelo de flujo base, pues se considerara que al inicio del modelamiento no
hay flujo superficial en la cuenca.
Datos considerados:
Respecto al CN se eligió el tipo III del grupo de antecedente de humedad (AMC) ya que si se
produce una gran avenida los días previos habrá estado lloviendo, y además de que es la más
desfavorable. El CN (III) se obtuvo en función del CN (II), según la relación mencionada líneas
arriba, mientras que el CN (II) se obtuvo en base a la Tabla N° 2, a la observación de imágenes
satelitales y mapas de cobertura vegetal, tal como se muestra en la siguiente figura:
Fig. N° 02 – Mapa de cobertura vegetal.
_ _
6
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
El siguiente cuadro muestra los datos utilizados en el modelado de la cuenca del programa HEC‐
HMS:
Cuadro N° 5 – Datos de entrada del modelo de cuenca
Modelo de transito de avenida
K: Medida del tiempo de transporte de una onda de un punto a otro, expresándose en
unidades de tiempo (horas), esta puede ser calculada como el tiempo de retardo del
tramo evaluado. Por tanto, se definirá según la siguiente ecuación:
.
0.18 .
Dónde:
K: Tiempo de transporte (hr).
L: Longitud del cauce principal (Km).
S: Pendiente del cauce principal (m/m).
X: Variable adimensional que define la atenuación del hidrograma y que se encuentra
entre 0.0 y 0.5 con un valor medio cercano a 0.3; X= 0 significa una máxima atenuación y
X = 0.5 significa no atenuación. Asimismo, no se necesita una gran precisión en la
determinación de X debido a que los resultados del método son relativamente insensibles
al valor de este parámetro (ref.: Universidad de Talca, Facultad de Ciencias Forestales,
Departamento de Gestión Forestal y Ambiental).
Número de tramos en que se dividirá el cauce.
_ _
7
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
El siguiente cuadro muestra los datos de entrada del módulo de transito de avenida, donde el
parámetro K se determinó según las características geométricas del tramo y el parámetro X se
asignará con un valor de 0.20 (hay que recordar que este valor no necesita una gran precisión
debido a que los resultados son relativamente insensibles a este parámetro):
Cuadro N° 6 – Datos de entrada del módulo de transito de avenida
MODELO METEOROLÓGICO
El modelo meteorológico contiene los datos de precipitación o bien registros históricos o
hipotéticos, respecto al modelo metrológico se emplearan los hietogramas de diseño obtenidos
en la estación Sicuani.
Cuadro N° 7 – Características de la estación pluviométrica seleccionada
30
24.37
25
20
Precipitación (mm)
15
10
4.61
3.09
5
2.39
1.98
1.70
1.50
1.35
1.22
1.13
1.04
0.98
0.92
0.87
0.82
0.78
0.74
0.71
0.68
0.66
0.63
0.61
0.59
0.57
0
0 ‐ 60
960 ‐ 1020
1020 ‐ 1080
1080 ‐ 1140
1140 ‐ 1200
1200 ‐ 1260
1260 ‐ 1320
1320 ‐ 1380
1380 ‐ 1440
60 ‐ 120
120 ‐ 180
180 ‐ 240
240 ‐ 300
300 ‐ 360
360 ‐ 420
420 ‐ 480
480 ‐ 540
540 ‐ 600
600 ‐ 660
660 ‐ 720
720 ‐ 780
780 ‐ 840
840 ‐ 900
900 ‐ 960
Tiempo (min)
Grafico N° 1 – Hietograma de Diseño para un periodo de retorno de 50 años.
_ _
8
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
35
27.87
30
25
Precipitación (mm)
20
15
10
5.27
3.54
2.74
5
2.26
1.94
1.71
1.54
1.40
1.29
1.20
1.12
1.05
0.99
0.94
0.89
0.85
0.81
0.78
0.75
0.72
0.70
0.67
0.65
0
0 ‐ 60
60 ‐ 120
120 ‐ 180
180 ‐ 240
240 ‐ 300
300 ‐ 360
360 ‐ 420
420 ‐ 480
480 ‐ 540
540 ‐ 600
600 ‐ 660
660 ‐ 720
720 ‐ 780
780 ‐ 840
840 ‐ 900
900 ‐ 960
960 ‐ 1020
1020 ‐ 1080
1080 ‐ 1140
1140 ‐ 1200
1200 ‐ 1260
1260 ‐ 1320
1320 ‐ 1380
1380 ‐ 1440
Tiempo (min)
Grafico N° 2 – Hietograma de Diseño para un periodo de retorno de 140 años.
ESPECIFICACIONES DE CONTROL
Las especificaciones de control contienen todo el calendario de información para el modelo,
incluyendo el tiempo de comienzo y el tiempo de término para el cálculo en la simulación
computacional. En las especificaciones de control se tendrá en cuenta que el tiempo de
modelamiento será mayor a la duración del Hietograma de diseño.
_ _
9
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
RESULTADOS DE LA SIMULACIÓN HIDROLÓGICA (Tr= 50 años)
De acuerdo a la combinación del modelo de la cuenca (model basin), modelo meteorológico y
especificaciones de control, se definen las respectivas corridas (run) para un tiempo de retorno de
50 años, se tiene como resultado el siguiente hidrograma donde se observa que el caudal máximo
obtenido es de 550.90 m3/s.
Fig. N° 03 – Resultados HMS para 50 años.
_ _
10
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
Fig. N° 04 – Hidrograma para 50 años.
_ _
11
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
RESULTADOS DE LA SIMULACIÓN HIDROLÓGICA (Tr= 140 años)
De acuerdo a la combinación del modelo de la cuenca (model basin), modelo meteorológico y
especificaciones de control, se definen las respectivas corridas (run) para un tiempo de retorno de
140 años, se tiene como resultado el siguiente hidrograma donde se observa que el caudal
máximo obtenido es de 765.50 m3/s.
Fig. N° 05 – Resultados HMS para 140 años.
_ _
12
CONSORCIO CHAKAKUNA QOSCO
ESTUDIO DE HIDROLOGIA E HIDRAULICA
_
Fig. N° 06 – Hidrograma para 140 años.
_ _
13
CONSORCIO CHAKAKUNA QOSCO