Universidad Católica de Salta – Carrera: Ingeniería Civil

Curso: Hidrología de Superficie – Año 2008

Unidad 5: Escorrentía y Análisis de Hidrogramas
5.1

UNIDAD 5: Escorrentía y Análisis de hidrogramas

5.1 HIDROGRAMA UNITARIO

El Hidrograma Unitario (HU) de una cuenca es el hidrograma de escurrimiento directo a la

salida de la misma, originado por una lluvia efectiva de altura unitaria (1 mm, 1 cm ó 1 plg)
distribuida uniformemente en el espacio (área de la cuenca) y en el tiempo (intensidad
constante), durante una duración especificada.

Si la duración de la lluvia efectiva es finita, se denomina Hidrograma Unitario de Tiempo

Discreto (HUT) y se simboliza HU (1, t u ), donde 1 es la altura de lluvia efectiva total y t u es
la duración de lluvia especificada. Si la duración es infinitesimal, se denomina Hidrograma
Unitario Instantáneo (HUI) y se simboliza u (t) (Figura Nº 5.1).

Figura Nº 5.1

5.1.1 Hipótesis simplificativas

a) La lluvia efectiva es uniformemente distribuida en el tiempo (intensidad constante) para

la duración especificada. Por lo tanto, para su determinación se deben elegir tormentas
de corta duración y alta intensidad, debido a que es más probable que éstas produzcan
una intensidad relativamente constante.

b) La lluvia efectiva es uniformemente distribuida en el espacio. Para que se cumpla

razonablemente la hipótesis de uniformidad espacial de la lluvia, el tamaño de la cuenca
debe ser pequeño. El tamaño máximo de aplicación depende del tipo de tormenta
(frontal o convectiva) y de la precisión deseada en los resultados. A pesar de que el
modelo fue desarrollado originalmente para cuencas grandes (Sherman, 1932), Brater
(1940) demostró que el HU es aplicable a cuencas chicas, desde 0.5 km2 hasta 25 km2.
Aunque otros autores admiten tamaños mayores, como 2000 km2 (OMM, 1983) y
5000km2 (Linsley, 1977), el HU debe aplicarse a cuencas pequeñas1, a menos que sea
aceptable una precisión inferior.

1
Cuencas cuya sensibilidad a lluvias de corta duración y alta intensidad y al uso del suelo no es suprimida por
las características del canal (Chow, 1964).

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.2

Para aplicar el HU a una cuenca grande, se debe desagregar la misma en subcuencas

pequeñas, aplicar HU a cada una de éstas y luego combinar los hidrogramas parciales
mediante una propagación de crecidas y suma de hidrogramas hacia aguas abajo.

c) El tiempo de base de los hidrogramas de escurrimiento directo (HED) originados por

lluvias efectivas de igual duración es constante. En un hidrograma de flujo de agua
superficial se cumple:

tb = tr + tc (5.1)

t b : tiempo de base del HED (depende del método de separación del flujo base), t c : tiempo de
concentración, constante según la hipótesis de sistema lineal, t r : tiempo de lluvia efectiva.

Para una duración t r dada, para que t b permanezca constante entre lluvias, t c debe
mantenerse constante. Esto implica una hipótesis de sistema lineal. Sherman clasificó la
escorrentía en superficial y subterránea y definió el HU para ser usado únicamente con
la escorrentía superficial. No obstante, en la práctica se extiende su aplicación al cálculo
de hidrogramas resultantes de lluvias efectivas.

d) Se asume que la cuenca es un sistema lineal, es decir se cumplen las propiedades de

superposición y proporcionalidad.

e) Se asume que la cuenca es un sistema invariante en el tiempo, es decir, sus parámetros

no cambian con el tiempo.

El modelo no es aplicable a la escorrentía originada por fusión de nieve.

5.2 HIDROGRAMA UNITARIO DE TIEMPO DISCRETO

Aplicando las propiedades de proporcionalidad y superposición, se puede determinar el HED

a la salida de una cuenca por convolución entre el hietograma de lluvia efectiva (HPE) y la
función núcleo HUT, de acuerdo al siguiente procedimiento (Figura Nº 5.2):

i) Se desagrega el HPE en intervalos tu y se considera una intensidad de lluvia constante

en cada intervalo.

ii) Se determinan los hidrogramas elementales generados por cada bloque del HPE
multiplicando las ordenadas del HUT por la lluvia efectiva, Pi, correspondiente a cada
bloque i (proporcionalidad).

iii) Se calcula el HED total sumando, para cada tiempo, las ordenadas de los hidrogramas
elementales (superposición).

Este procedimiento implica aplicar el sistema de ecuaciones lineales:

Q0 = 0
Q1 = P1 U1
Q 2 = P2 U1 + P1 U2
(5.2)
Q 3 = P3 U1 + P2 U 2 + P1 U3
.........................

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.3

Generalizando:

t
Q t = ∑ U t - i + 1 Pi (5.3)
i =1

Q t : ordenada del HED al tiempo t, con dimensiones [L/T] (ej. mm/h) ó [L3/T] (ej. m3/s); U t − i + 1 :
ordenada del HUT al tiempo t – i + 1, con dimensiones [1/T] (ej. 1/h) ó [L3 / (T L)] (ej. m3/s/mm), Pi :
ordenada del HPE al tiempo i [L] (ej. mm).

Figura Nº 5.2: Determinación del HED por convolución discreta entre HPE y HUT.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.4

La ecuación (5.3) se denomina ecuación de convolución discreta. Comparando esta

ecuación con (1.13), se observa que: a) HED es la función de salida, Y(t); b) HPE es la
función de entrada, X(t); c) HUT es la función de respuesta de la cuenca, h(t), y d) la
sumatoria de productos es la operación Φ, denominada convolución para modelos lineales.

Las funciones HPE y HUT deben estar definidas con un mismo intervalo de tiempo, igual a
t u . Se cumple la condición:

NU = NQ – NP + 1 (5.4)

NU: número de ordenadas no nulas del HUT, NQ: número de ordenadas no nulas del HED, NP:
número de ordenadas no nulas del HPE.

Las hipótesis de linealidad e invariancia en el tiempo son fundamentales en la teoría del HU.
Las otras hipótesis son hechas por conveniencia práctica pero no son esenciales desde el
punto de vista teórico.

El flujo de agua superficial es un fenómeno no lineal. Además, todo sistema hidrológico es

dinámico y sufre modificaciones de origen natural y/o antrópico. Las naturales son
originadas por variaciones periódicas y estacionales de las variables hidrometeorológicas
que influyen sobre ciertos parámetros del sistema (ej. cambios de la cobertura vegetal),
mientras que las antrópicas son debidas al uso y tratamiento del suelo y a la construcción de
obras, como canalizaciones, presas de regulación, urbanización y otras.

Por lo tanto, las hipótesis de linealidad e invariancia en el tiempo no son absolutamente

ciertas. Sin embargo, en muchas ocasiones los modelos lineales satisfacen la precisión
deseada de los resultados por una compensación de las distintas no linealidades del
sistema (sublineales y supralineales). Debido a que introducen una gran simplificación en el
tratamiento matemático del flujo directo, los modelos lineales son de frecuente aplicación
práctica en hidrología de superficie.

5.2.1 Duración de referencia de la lluvia efectiva

La duración de referencia de la lluvia efectiva que se adopta para el HUT se denomina

tiempo unitario ( t u ). El tiempo unitario es la mayor duración de la lluvia efectiva que no
produce un incremento apreciable del tiempo de base del hidrograma resultante. La Figura
Nº 5.3 muestra cuatro HPE y los correspondientes HED sobre una cuenca dada. En todos
los casos, la precipitación efectiva total se mantiene constante y varían la duración de la
lluvia tr y, en consecuencia, la intensidad h.

Figura Nº 5.3: Concepto de tiempo unitario.

Q: caudal directo, i: intensidad de lluvia efectiva.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.5

A partir de tr1, se observa que a medida que disminuye la duración de la lluvia efectiva,
disminuye el tiempo de base y aumenta el pico del hidrograma. Sin embargo, para
duraciones de lluvia menores o iguales a tr3 no se producen modificaciones apreciables del
tiempo de base de los hidrogramas resultantes. Por lo tanto, tr3 es el tiempo unitario de la
cuenca en estudio.

El tiempo unitario varía en general entre 1 y 24 h, según las características de la cuenca

(tamaño, pendiente y otras). La Tabla Nº 5.1 muestra criterios sugeridos por distintos
autores para adoptar t u .

Tabla Nº 5.1: Criterios para adoptar t u .

Autor Area [Km²] t u [h]

> 2600 12 a 24
260 - 2600 6 a 12
Sherman 50 - 260 2
< 50 1/3 a 1/4 t c
Remenieras t u = 1/3 a 1/5 t c
Linsley t u = 1/4 t lag

t c : tiempo de concentración de la cuenca, t lag : tiempo de retardo (tiempo comprendido entre los
centros de gravedad del HPE y del HED).

Es conveniente adoptar un t u suficientemente pequeño respecto al tiempo de concentración

de la cuenca ( t u ≤ t c / 5), de manera que la curva de crecida del hidrograma esté definida
por una cantidad suficiente de puntos.

5.2.2 Determinación de un HUT en base a los datos HPE – HED de un evento

Si se dispone de un HPE medio areal sobre la cuenca y del correspondiente HED a la

salida, se puede determinar un HUT.

Se pueden presentar dos casos, de acuerdo a la relación entre t r y t u :

a) Si t r < t u : duraciones de lluvia menores que t u no originan modificaciones significativas

del hidrograma resultante. Por lo tanto, se asume que el HED ha sido generado por una
lluvia de duración t u y el HUT se calcula dividiendo cada ordenada del HED por la
precipitación efectiva total (se aplica la propiedad de proporcionalidad):

Qt
Ut = (5.5)
P

U t : ordenada del HUT al tiempo t, Q t : ordenada del HED al tiempo t, P: lluvia efectiva total
del evento.

b) Si t r > t u : se aplica la propiedad de proporcionalidad y luego el método del hidrograma

en S para llevar la duración a t u (Figura Nº 5.4).

El hidrograma en S (HS) es el hidrograma resultante de una lluvia efectiva continua de

intensidad constante y duración mayor a t c .

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.6

Para determinar el HS, se calcula el hidrograma H (1, t r ) dividiendo las ordenadas de H (P,
t r ) por P (proporcionalidad). Luego, se repite n veces el hidrograma H (1, t r ) desplazándolo
en cada ocasión en un tiempo t r y se suman las ordenadas de los hidrogramas parciales
(propiedad de superposición). La cantidad n se adopta de modo que el HS alcance un
caudal máximo de equilibrio. El tiempo de concentración del flujo directo es ( t ′c = t b − t r ) ,
por lo que n se calcula:

n=
(t b − tr )
(5.6)
tr

En general, debe hacerse un ajuste gráfico de la curva en S pues los puntos no se alinean
perfectamente sobre ella, sino que tienden a fluctuar alrededor del caudal máximo, debido a
que la hipótesis de intensidad constante no se cumple exactamente en la realidad.

Figura Nº 5.4: Aplicación del hidrograma en S.

I: intensidad de lluvia efectiva (I = 1 / tr), Q: caudal directo.

Para un tiempo igual a t c se alcanza el caudal máximo o caudal de equilibrio, ya que a

partir del mismo aporta toda el área a la salida:

Q máx = i A (5.7)

i: intensidad de lluvia, A: área de la cuenca.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.7

Para obtener el HUT (1, t u ), se desplaza el HS en un tiempo t u , se calcula la diferencia

entre los HS original y desplazado y finalmente se convierte la lámina a la unidad.

Para cada evento observado, se puede determinar un HUT.

5.2.3 Determinación del HU característico de una cuenca

El HU característico de una cuenca (HUC) se calcula a partir de los HUT obtenidos de

distintas tormentas, lo cual minimiza posibles errores de los HU individuales por
compensación de los mismos.

Para ello se calculan los parámetros tiempo al pico, caudal pico y tiempo de base medios
t p , Q p y t b . Se desplazan los HUT en el sentido de los tiempos hasta poner en fase los
picos individuales con t p y se traza el HUT característico, pasando por los 3 puntos medios
y promediando las formas individuales (Figura Nº 5.5).

Figura Nº 5.5: Hidrograma unitario característico.

El HUC y los HUT deben verificar una lámina de escurrimiento unitaria.

Ejemplo Nº 5.1: A partir de los siguientes datos morfológicos de una cuenca y de los datos
HPE-HED de 2 eventos (Tabla Nº 5.2), se pide determinar el HUC (1 mm, tu) de la misma.

A= 33.84 km2
L= 5 km
ΔH = 2 m
S= 0.00040 m/m

A: área de la cuenca, L: longitud del cauce principal, ΔH: desnivel entre la cabecera y la salida de la
cuenca, S: pendiente media del cauce principal.

Tabla Nº 5.2:
EVENTO Nº 1 EVENTO Nº 2
Tiempo HPE Tiempo HPE
[h] [mm] [h] [mm]
0 0.0 0 0.0
1 5.0 1 6.0
Total 5.0 2 8.0
Total 15.0

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.8

EVENTO Nº 1 EVENTO Nº 2
Tiempo HED Tiempo HED
[h] [m3/s] [h] [m3/s]
0 0.0 0 0.0
1 1.5 1 1.4
2 5.5 2 5.6
3 11.5 3 11.2
4 8.5 4 22.4
5 6.0 5 26.6
6 5.5 6 21.0
7 3.5 7 16.8
8 2.5 8 11.2
9 1.5 9 7.0
10 1.2 10 5.2
11 0.8 11 2.8
12 0.0 12 1.4
Ed [hm3] 0.169 13 0.0
Ed [mm] 5.0 Ed [hm3] 0.474
Ed [mm] 15.0

Ed: escurrimiento directo.

Se adopta tu. Para ello, se calcula el tiempo de concentración (flujo superficial) de la cuenca
con la fórmula de Kirpich (1.5):

tc (flujo sup) = 5.7 h

tc/5 = 0.9 h
Se adopta tu = 1.0 h

Se determina el HU (1 mm, 1h) para el evento Nº 1:

tr = tu = 1 h , por lo que sólo se aplica proporcionalidad (5.5) (Tabla Nº 5.3).

Tabla Nº 5.3:
Tiempo HED HU
(5 mm, 1 h) (1mm, 1h)
[h] [m3/s] [m3/s/mm]
0 0.0 0.00
1 1.5 0.30
2 5.5 1.10
3 11.5 2.30
4 8.5 1.70
5 6.0 1.20
6 5.5 0.90
7 3.5 0.70
8 2.5 0.50
9 1.5 0.30
10 1.2 0.24
11 0.8 0.16
12 0.0 0.00
Verif. 1.0

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.9

Se determina el HU para el evento Nº 2. tr = 2 h > tu = 1 h , por lo que se aplica

proporcionalidad e Hidrograma en S (Tabla Nº 5.4).

Se calcula H (1mm, 2h), dividiendo cada ordenada del HED (14mm, 2h) por la lluvia total del
evento (col. 3).

Se repite n veces H (1 mm, 2 h), donde n = 6 (5.6) (columnas 3 a 8).

Se determina el HS como suma de los H (1mm,2h) repetidos (col. 9). Al tiempo (13 h – 2 h =
11 h) el HS alcanza el caudal máximo de equilibrio: Qmáx = 1 mm / 2 h * A [km2] = 5.7 m3/s
(5.7).

Se desfasa HS en un tiempo tu = 1 h (col. 10).

Se calcula el hidrograma Hdif, como diferencia entre el HS original y el desplazado (col. 11)
y se lleva la altura de lluvia de Hdif a la unidad, multiplicando las ordenadas por tr/tu (2h/1h =
2) (col. 12).

Tabla Nº 5.4:

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12)
Tiempo HED H H H H H H HS HS H dif. HU
(14 mm, (1 mm, (1 mm, (1 mm, (1 mm, (1 mm, (1 mm, (0.5 mm, (1 mm,
2 h) 2 h) 2 h) 2 h) 2 h) 2 h) 2 h) des. 1h) 1h)
[h] [m3/s] [*] [*] [*] [*] [*] [*] [*] [*] [*] [*]
0 0.0 0.0 0.0 0.0 0.0
1 1.4 0.1 0.1 0.0 0.1 0.2
2 5.6 0.4 0.0 0.4 0.1 0.3 0.6
3 11.2 0.8 0.1 0.9 0.4 0.5 1.0
4 22.4 1.6 0.4 0.0 2.0 0.9 1.1 2.2
5 26.6 1.9 0.8 0.1 2.8 2.0 0.8 1.6
6 21.0 1.5 1.6 0.4 0.0 3.5 2.8 0.7 1.4
7 16.8 1.2 1.9 0.8 0.1 5.0 3.5 0.5 1.0
8 11.2 0.8 1.5 1.6 0.4 0.0 5.3 5.0 0.3 0.6
9 7.0 0.5 1.2 1.9 0.8 0.1 5.5 5.3 0.2 0.4
10 5.2 0.3 0.8 1.5 1.6 0.4 0.0 5.6 5.5 0.1 0.2
11 2.8 0.2 0.5 1.2 1.9 0.8 0.1 5.7 5.6 0.1 0.2
12 1.4 0.1 0.3 0.8 1.5 1.6 0.4 5.7 5.7 0.0 0.0
13 0.0 0.0 0.2 0.5 1.2 1.9 0.8 5.7 5.7
14 0.1 0.3 0.8 1.5 1.6 5.7 5.7
15 0.0 0.2 0.5 1.2 1.9 5.7 5.7
16 0.1 0.3 0.8 1.5 5.7 5.7
17 0.0 0.2 0.5 1.2 5.7 5.7
18 0.1 0.3 0.8 5.7 5.7
19 0.0 0.2 0.5 5.7 5.7
20 0.1 0.3 5.7 5.7
21 0.0 0.2 5.7 5.7
22 0.1 5.7 5.7
23 0.0 5.7 5.7
Verif. 1.0

[*] = [m3/s/mm]

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.10

Figura Nº 5.6:

5.0

4.5

4.0

3.5
CAUDAL [m3/s]

3.0

2.5

2.0

1.5

1.0

0.5

0.0
0 2 4 6 8 10 12 14
TIEMPO [h]
HS HS des H dif. HU

Se determina el Hidrograma Unitario Característico (Tabla Nº 5.5). Para ello, se calculan: tp

med = 3.5 h (se adopta 5.0 h), Qp med = 2.25 [m3/s] y tb med = 12 h.

Tabla Nº 5.5:
Tiempo Evento 1 Evento 2
HU HU HUC
(1 mm, (1 mm, (1 mm,
1h) 1h) 1h)
[h] [*] [*] [*]
0 0.0 0.00
1 0.0 0.2 0.10
2 0.3 0.6 0.45
3 1.1 1.0 1.05
4 2.3 2.2 2.25
5 1.7 1.6 1.65
6 1.2 1.4 1.30
7 0.9 1.0 0.95
8 0.7 0.6 0.65
9 0.5 0.4 0.45
10 0.3 0.2 0.25
11 0.2 0.2 0.22
12 0.2 0.0 0.00
Verif. 1.0 1.0 1.0

[*] = [m3/s/mm]

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.11

Figura Nº 5.7:

2.5

2.0
CAUDAL [m3/s]

1.5

1.0

0.5

0.0
0 2 4 6 8 10 12 14
TIEMPO [h]
HU evento 1 HU evento 2 HUC

Se calcula el HED para el evento 2 aplicando (5.3) (columnas 3 a 5 de la Tabla Nº 5.6). Para
ello, se consideran el HPE de dicho evento y el HUC obtenido previamente (col. 2). Se
comparan los hidrogramas calculado y observado (HEDo).

Tabla Nº 5.6

(1) (2) (3) (4) (5) (6)

Tiempo HUC H1 H2 HEDc HEDo
(1 mm, (6 mm, (8 mm, (14 mm, (14 mm,
1h) 1h) 1h) 1h) 1h)
[h] [*] [m3/s] [m3/s] [m3/s] [m3/s]
0 0.0 0.0 0.0 0.0
1 0.1 0.6 0.0 0.6 1.4
2 0.5 2.7 0.8 3.5 5.6
3 1.1 6.3 3.6 9.9 11.2
4 2.3 13.5 8.4 21.9 22.4
5 1.7 9.9 18.0 27.9 26.6
6 1.3 7.8 13.2 21.0 21.0
7 1.0 5.7 10.4 16.1 16.8
8 0.7 3.9 7.6 11.5 11.2
9 0.4 2.7 5.2 7.9 7.0
10 0.3 1.5 3.6 5.1 5.2
11 0.2 1.3 2.0 3.3 2.8
12 0.0 0.0 1.8 1.8 1.4
13 0.0 0.0 0.0

[*] = [m3/s/mm]

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.12

Figura Nº 5.8:

30.0

25.0

20.0
CAUDAL [m3/s]

15.0

10.0

5.0

0.0
0 2 4 6 8 10 12 14
TIEMPO [h]
HEDc HEDo

5.3.3 Modelos Conceptuales de HUI

Modelo de Embalse Lineal Simple

Este modelo simula el flujo directo a la salida de la cuenca como el flujo a la salida de un
embalse lineal (Zoch, 1934) (Figura Nº 5.15).

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.13

Figura Nº 5.15

Utiliza la ecuación de continuidad y una relación lineal entre almacenamiento y caudal de

salida:

dS
= I− Q (5.51)
dt

S =K Q (5.52)

S: almacenamiento, t: tiempo, I: caudal de entrada al embalse, Q: caudal de salida del embalse, K:

coeficiente de embalse, representa el tiempo de vaciado del embalse para un caudal constante Q.

Derivando (5.52) con respecto al tiempo e igualando miembro a miembro con (5.51):

dQ
K + Q = I (5.53)
dt

Como se ha visto, un hietograma de intensidad constante y unitaria genera como salida la

función impulso unitario q (t). Si se aplica una entrada de intensidad constante y unitaria
(I=1), el caudal de salida es q (t):

dq
K + q =1 (5.54)
dt

La solución de esta ecuación es:

t
− (5.55)
q=c e K +1

Como condición inicial se tiene: t = 0, q = 0, de donde: c = -1. Por lo tanto:

t
− (5.56)
q =1− e K

Derivando con respecto al tiempo se tiene la expresión del HUI del modelo (Figura Nº 5.16):

t
1 −K (5.57)
u( t ) = e
K

El modelo tiene un único parámetro: K. La ordenada al origen de u(t) es igual a 1/K y u(t)
tiende a cero a medida que t tiende a ∞.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.14

Figura Nº 5.16

Como la expresión analítica del modelo es sencilla, es posible determinar las ordenadas del
HUT. Reemplazando (5.57) en (5.50):

1
(1- τ ) −
1
1 - K K (5.58)
U1 = ∫ e dτ = 1 − e
0 K

1
(2 - τ ) −
1 ⎛
−
1⎞
−
1
1 - K ⎜ ⎟
U2 = ∫ e dτ = e K ⎜1 − e K ⎟ = e K U1 (5.59)
0 K ⎜ ⎟
⎝ ⎠

de donde:
1
− (5.60)
Ui = e K Ui − 1

Válido para i ≥ 2.

Aplicando (5.3):
⎛ −
1⎞
⎜ K ⎟
Q1 = P1 ⎜1 − e ⎟ (5.61)
⎜ ⎟
⎝ ⎠

⎛ −
1⎞
−
1 ⎛ −
1⎞ ⎛ −
1⎞
−
1
⎜ K ⎟ K ⎜ K ⎟ ⎜ K ⎟ K
Q 2 = P2 ⎜1 − e ⎟ + P1 e ⎜1 − e ⎟ = P2 ⎜1 − e ⎟ + Q1 e (5.62)
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ ⎠ ⎝ ⎠ ⎝ ⎠

Generalizando:
⎛ −
1⎞
−
1
⎜ K ⎟ K
Q t + 1 = Pt + 1 ⎜1 − e ⎟ + Qt e (5.63)
⎜ ⎟
⎝ ⎠

Q: ordenada del HED [L/Δt], P: precipitación efectiva en el intervalo [L], K [Δt].

De acuerdo a (5.63), el modelo calcula Q t + Δt como la suma de dos términos, lo cual puede
visualizarse en la Figura Nº 5.17 para t = 1 (t + Δt = 2). El primer término del segundo
miembro (*1) es el generado por la lluvia efectiva caída en el intervalo (1, 2) y el segundo
(*2), es originado por la recesión del caudal al tiempo t = 1.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.15

Figura Nº 5.17

Se puede observar en (5.63) que para tiempos posteriores al fin de la lluvia, se anula el
primer componente y sólo actúa el segundo. En consecuencia, con este modelo no se
pueden representar hidrogramas que tengan el tiempo al pico mayor a la duración de la
lluvia. Esto limita la aplicación de este modelo a cuencas pequeñas, con poco
amortiguamiento, como es el caso de pequeñas cuencas urbanas.

Reemplazando (5.57) en (5.17), se obtiene el tiempo de retardo o esperanza matemática o

momento al origen de orden 1 del HUI:

t
∞ 1 -K
t L = E( t ) = MU1 = ∫ e t dt = K (5.64)
0 K

En sistemas lineales e invariantes, se cumple:

MU1 = MQ1 − MP1 (5.65)

MU2 = MQ 2 − MP2 (5.66)

MU1 , MQ1 y MP1 : momentos de orden 1 (al origen) del HUI, HED y HPE, respectivamente, MU2 ,
MQ 2 y MP2 : momentos de orden 2 (centrados respecto al centroide) del HUI, HED y HPE,
respectivamente.

∑ Qi t i
MQ1 = (5.67)
∑ Qi

∑ Pi t i
MP1 = (5.68)
∑ Pi

2
∑ Q i t ′i
MQ 2 = (5.69)
∑ Qi

2
∑ Pi t ′i
MP2 = (5.70)
∑ Pi

t i : tiempo genérico i medido al origen de los tiempos, t ′i : tiempo genérico i medido respecto del
centroide.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.16

Con (5.65) se puede estimar un valor inicial del parámetro K (Figura Nº 5.18).

Figura Nº 5.18

El valor final del parámetro se obtiene por calibración. El procedimiento de calibración más
frecuentemente usado es por prueba y error, variando el valor del parámetro en cada corrida
se procura alcanzar la mejor correspondencia del hidrograma calculado con el observado.

Con el objeto de auxiliar la calibración, para cada hidrograma simulado se pueden calcular el
coeficiente de eficiencia, E f (Nash y Sutcliffe, 1971) y el error estándar, E e :

n
o
(
c 2
∑ Q t −Q t
t =1
)
Ef = 1 − (5.71)
∑
n

t =1
(
Q ot − Q o 2)
n
o
( c
∑ Q t −Q t
t =1
)2
(5.72)
Ee =
n

Q ot y Q ct : caudales observado y calculado al tiempo t, Q o : caudal observado medio, n: número de

intervalos de tiempo,.

El hidrograma simulado representa adecuadamente al observado cuando E f se aproxima a la

unidad y E e es suficientemente pequeño respecto a los caudales observados. Por lo tanto, se
adopta como valor final del parámetro a aquél que maximiza E f . E f representa la fracción de
la varianza de los valores observados que es explicada por el modelo. Puede variar desde
menos infinito a 1. Un valor igual a 1 indica que el ajuste del modelo es perfecto, mientras
que un valor 0 o negativo indica que la media de los caudales observados es igual o mejor
estimador que el modelo, respectivamente. Valores superiores a 0.7, generalmente indican
un ajuste aceptable.

Si se considera a la cuenca como un sistema lineal e invariante, K es constante. Esta

hipótesis simplificativa no es totalmente cierta y K varía entre distintos eventos de una
misma cuenca.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.17

Ejemplo Nº 5.3: Dados el HPE y el HED (columnas 1, 2, 4 y 5 de la Tabla Nº 5.8), medidos

en una cuenca urbana de 200.2 ha de superficie, se pide calcular el hidrograma generado
por dicho evento aplicando el modelo de Zoch. Calcular el coeficiente de eficiencia y hacer
un análisis de sensibilidad del modelo variando el parámetro K en ± 50%. El tc de la cuenca
es 60 min.

Se calculan MQ1 (5.67) y MP1 (5.68) (cálculos auxiliares en columnas (3), (6) y (7) de la
Tabla Nº 5.8). Para los caudales (valores instantáneos) se consideran los tiempos a cada
caudal. Para las lluvias (valores incrementales cada 5 min), se consideran los tiempos al
instante medio de cada bloque de lluvia (tpi).

Tabla Nº 5.8
(1) (2) (3) (4) (5) (6) (7)
Tiempo HED Qi*ti Tiempo HPE tpi Pi * tpi
[h] [m3/s] [m3/s*h] [h] [mm] [h] [mm*h]
0.000 0.000 0.083 0.2 0.042 0.008
0.167 0.052 0.009 0.167 0.6 0.125 0.075
0.333 0.125 0.042 0.250 0.8 0.208 0.167
0.500 0.224 0.112 0.333 0.5 0.292 0.146
0.667 0.346 0.231 0.417 0.6 0.375 0.225
0.833 0.650 0.542 0.500 0.5 0.458 0.229
1.000 0.915 0.915 0.583 0.7 0.542 0.379
1.167 1.110 1.295 0.667 0.3 0.625 0.188
1.333 1.230 1.640 0.750 0.1 0.708 0.071
1.500 1.246 1.869 0.833 0.2 0.792 0.158
1.667 1.230 2.050 0.917 0.4 0.875 0.350
1.833 1.150 2.108 1.000 0.1 0.958 0.096
2.000 1.050 2.100
2.167 0.953 2.065
2.333 0.843 1.968
2.500 0.755 1.888
2.667 0.675 1.800
2.833 0.590 1.672
3.000 0.516 1.548
3.167 0.443 1.403
3.333 0.382 1.273
3.500 0.335 1.173
3.667 0.285 1.045
3.833 0.242 0.928
5.000 0.211 0.844
5.167 0.181 0.754
5.333 0.155 0.672
5.500 0.131 0.590
5.667 0.115 0.537
5.833 0.097 0.469
5.000 0.080 0.400
5.167 0.070 0.362
5.333 0.060 0.320
5.500 0.050 0.275
5.667 0.040 0.227
5.833 0.030 0.175
6.000 0.025 0.150
6.167 0.020 0.123

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.18

(1) (2) (3) (4) (5) (6) (7)

Tiempo HED Qi*ti Tiempo HPE tpi Pi * tpi
[h] [m3/s] [m3/s*h] [h] [mm] [h] [mm*h]
6.333 0.015 0.095
6.500 0.010 0.065
6.667 0.008 0.053
6.833 0.004 0.027
7.000 0.000 0.000

SUM 16.649 35.810 5.0 2.092

MQ1 = 2.15 h y MP1 = 0.42 h. Se calcula K (5.65): K = 1.73 h.

tc = 60 min, tc / 5 = 12 min ⇒ se adopta Δt (tu) = 10 min y se calcula HED (5.63) (col. 3 de

Tabla Nº 5.9):

Tabla Nº 5.9

(1) (2) (3) (4)

Análisis de sensibilidad
Qc (K - 50%) Qc (K + 50%)
K[h] = 1.73 0.9 2.6
K[Δt] = 10.4 5.2 15.6
exp(-1/K) = 0.908 0.825 0.938
Tiempo HPE Qc Qc Qc Qc Qc Qc
[h] [mm/Δt] [mm/Δt] [m3/s] [mm/Δt] [m3/s] [mm/Δt] [m3/s]
0.000 0.0 0.000 0.000 0.000 0.000 0.000 0.000
0.167 0.8 0.073 0.245 0.140 0.467 0.050 0.166
0.333 1.3 0.186 0.620 0.343 1.145 0.127 0.425
0.500 1.1 0.270 0.900 0.476 1.587 0.188 0.627
0.667 1.0 0.337 1.123 0.567 1.893 0.238 0.795
0.833 0.3 0.333 1.112 0.521 1.737 0.242 0.808
1.000 0.5 0.349 1.163 0.517 1.725 0.258 0.861
1.167 0.317 1.057 0.426 1.423 0.242 0.808
1.333 0.288 0.960 0.352 1.174 0.227 0.758
1.500 0.261 0.872 0.290 0.968 0.213 0.710
1.667 0.237 0.792 0.239 0.799 0.200 0.666
1.833 0.216 0.719 0.198 0.659 0.187 0.625
2.000 0.196 0.653 0.163 0.544 0.176 0.586
2.167 0.178 0.593 0.134 0.449 0.165 0.550
2.333 0.161 0.539 0.111 0.370 0.155 0.516
2.500 0.147 0.489 0.092 0.305 0.145 0.484
2.667 0.133 0.445 0.075 0.252 0.136 0.453
2.833 0.121 0.404 0.062 0.208 0.127 0.425
3.000 0.110 0.367 0.051 0.171 0.120 0.399
3.167 0.100 0.333 0.042 0.141 0.112 0.374
3.333 0.091 0.303 0.035 0.117 0.105 0.351
3.500 0.082 0.275 0.029 0.096 0.099 0.329
3.667 0.075 0.250 0.024 0.079 0.093 0.309
3.833 0.068 0.227 0.020 0.066 0.087 0.289
5.000 0.062 0.206 0.016 0.054 0.081 0.271
5.167 0.056 0.187 0.013 0.045 0.076 0.255
5.333 0.051 0.170 0.011 0.037 0.072 0.239

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.19

(1) (2) (3) (4)

Análisis de sensibilidad
Qc (K - 50%) Qc (K + 50%)
K[h] = 1.73 0.9 2.6
K[Δt] = 10.4 5.2 15.6
exp(-1/K) = 0.908 0.825 0.938
Tiempo HPE Qc Qc Qc Qc Qc Qc
[h] [mm/Δt] [mm/Δt] [m3/s] [mm/Δt] [m3/s] [mm/Δt] [m3/s]
5.500 0.046 0.154 0.009 0.030 0.067 0.224
5.667 0.042 0.140 0.008 0.025 0.063 0.210
5.833 0.038 0.127 0.006 0.021 0.059 0.197
5.000 0.035 0.116 0.005 0.017 0.055 0.185
5.167 0.031 0.105 0.004 0.014 0.052 0.173
5.333 0.029 0.095 0.003 0.012 0.049 0.163
5.500 0.026 0.087 0.003 0.010 0.046 0.152
5.667 0.024 0.079 0.002 0.008 0.043 0.143
5.833 0.021 0.071 0.002 0.007 0.040 0.134
6.000 0.019 0.065 0.002 0.005 0.038 0.126
...... ...... ...... ...... ...... ...... ......

Se calculan los errores relativos del caudal pico y del tiempo al pico, así como el coeficiente
de eficiencia (5.71) y el error estándar (5.72): ErQp = -6.6 %, Ertp = -33.3 %, Ef = 0.65 y Ee =
0.242 m3/s.

Figura Nº 5.19

2.000

1.800

1.600

1.400
CAUDAL [m3/s]

1.200

1.000

0.800

0.600

0.400

0.200

0.000
0.00 1.00 2.00 3.00 4.00 5.00 6.00
TIEMPO [h]
Qo Qc Qc (K - 50%) Qc (K + 50%)

Qc: HED calculado con el modelo, Qo: HED observado.

Se observa que el modelo no representa adecuadamente el hidrograma observado (Ef =

0.65). Esto es debido a que el modelo no puede representar tiempos al pico posteriores a la
duración de la lluvia, como ocurre en la cuenca en estudio (tp = 1.5 h, tr = 1.0 h). Del análisis
de sensibilidad, se obtiene que si K aumenta, el hidrograma simulado es más amortiguado.
Inversamente, si K disminuye, el hidrograma simulado es más empuntado.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.20

Modelo de Embalses Lineales en Cascada

Este modelo representa el escurrimiento superficial como el flujo a la salida de un conjunto

de n embalses lineales en cascada, todos con el mismo coeficiente de embalse K (Nash,
1957) (Figura Nº 5.20).

Figura Nº 5.20

Si al primer embalse de aguas arriba ingresa una precipitación efectiva unitaria instantánea,
la salida del mismo es (5.57):

t
1 −K (5.73)
u1( t ) = e
K

La entrada al segundo embalse es salida del primero. La salida del segundo se obtiene
aplicando (5.10):

τ (t − τ ) t
t1 −K 1 − K 1 t −K (5.74)
u2 (t) = ∫ e e dτ = e
0 K K K K

Para n reservorios:
(n − 1) −
t
1 ⎛t⎞ K 1 (5.75)
un ( t ) = ⎜ ⎟ e
K ⎝K ⎠ (n − 1)!
Para contemplar valores de n no enteros, se reemplaza el factorial de orden (n-1) por la
función gamma de orden n (Tabla Nº 5.10):

(n − 1) −
t
1 ⎛t⎞ 1 (5.76)
un ( t ) = ⎜ ⎟ e K
K ⎝K ⎠ Γn

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.21

Tabla Nº 5.10: Función Gamma

x Γ (x) x Γ ( x) x Γ (x) x Γ ( x)
1.00 1.00000 1.25 0.90640 1.50 0.88623 1.75 0.91906
1.01 0.99433 1.26 0.90440 1.51 0.88659 1.76 0.92137
1.02 0.98884 1.27 0.90250 1.52 0.88704 1.77 0.92376
1.03 0.98355 1.28 0.90072 1.53 0.88757 1.78 0.92623
1.04 0.97844 1.29 0.89904 1.54 0.88818 1.79 0.92877
1.05 0.97350 1.30 0.89747 1.55 0.88887 1.80 0.93138
1.06 0.96874 1.31 0.89600 1.56 0.88964 1.81 0.93408
1.07 0.96415 1.32 0.89464 1.57 0.89049 1.82 0.93685
1.08 0.95973 1.33 0.89338 1.58 0.89142 1.83 0.93969
1.09 0.95546 1.34 0.89222 1.59 0.89243 1.84 0.94261
1.10 0.95135 1.35 0.89115 1.60 0.89352 1.85 0.94561
1.11 0.94740 1.36 0.89018 1.61 0.89468 1.86 0.94869
1.12 0.94359 1.37 0.88931 1.62 0.89592 1.87 0.95184
1.13 0.93993 1.38 0.88854 1.63 0.89724 1.88 0.95507
1.14 0.93642 1.39 0.88785 1.64 0.89864 1.89 0.95838
1.15 0.93304 1.40 0.88726 1.65 0.90012 1.90 0.96177
1.16 0.92980 1.41 0.88676 1.66 0.90167 1.91 0.96523
1.17 0.92670 1.42 0.88636 1.67 0.90330 1.92 0.96877
1.18 0.92373 1.43 0.88604 1.68 0.90500 1.93 0.97240
1.19 0.92089 1.44 0.88581 1.69 0.90678 1.94 0.97610
1.20 0.91817 1.45 0.88566 1.70 0.90864 1.95 0.97988
1.21 0.91558 1.46 0.88560 1.71 0.91057 1.96 0.98374
1.22 0.91311 1.47 0.88563 1.72 0.91258 1.97 0.98768
1.23 0.91075 1.48 0.88575 1.73 0.91467 1.98 0.99171
1.24 0.90852 1.49 0.88595 1.74 0.91683 1.99 0.99581
1.25 0.90640 1.50 0.88623 1.75 0.91906 2.00 1.00000

Nota: Los valores de la función Gamma para x < 1 y x > 2 se pueden calcular mediante las
Γ (x + 1))
fórmulas Γ ( x ) = , Γ ( x ) = ( x − 1) Γ (x − 1) , respectivamente.
x

Los parámetros del modelo son n y K. Para estimar éstos, se plantean 2 momentos
estadísticos del HUI: el momento al origen de orden 1 ( MU1 ) o tiempo de retardo y el
momento centrado de orden 2 ( MU2 ) o varianza.

El MU1 es:
∞ (n − 1) −
t
1 ⎛t⎞ K 1 (5.77)
MU1 = ∫ ⎜ ⎟ e t dt = n K
0 K ⎝K ⎠ (n − 1)!
Como se ha descrito, t L representa el tiempo medio de llegada del escurrimiento a la salida
de la cuenca. La dispersión de los tiempos de llegada con respecto al valor medio se
expresa por el momento centrado de segundo orden o varianza:

( )
MU2 = VAR ( t ) = E t 2 − [E (t )]2 (5.78)

VAR(t): varianza o momento del HUI de orden 2 centrado respecto a la media, E(t²): momento al
origen de orden 2, E(t): momento al origen de orden 1.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.22

(n − 1) t
E (t ) = ∫ u ( t ) t
∞ ∞ 1 ⎛t⎞ − 1
2 2
dt = ∫ ⎜ ⎟ e K t 2 dt = K 2 n (n + 1) (5.79)
0 0 K ⎝K ⎠ (n − 1)!
Reemplazando (5.77) y (5.79) en (5.78) se obtiene:

MU2 = K 2 n (5.80)

Reemplazando MU1 y MU2 en (5.65) y (5.66) y operando, se obtienen las ecuaciones para
calcular los valores iniciales de los parámetros:

MQ 2 − MP2
K= (5.81)
MQ1 − MP1

n=
(MQ1 − MP1 )2 (5.82)
MQ 2 − MP2

Derivando u(t) respecto al tiempo e igualando a cero, se obtiene el tiempo al pico del HUI:

t p = (n − 1) K (5.83)

Figura Nº 5.21

La expresión del u(t) de Nash es una función Gamma de dos parámetros:

1
f ( x) = ma x a − 1 e − m x (5.84)
Γa

a: parámetro de forma (igual a n), m: parámetro de escala (igual a 1/K).

Una vez calculado el HUI, se puede obtener el HUT con (5.40). El modelo de Nash tiene una
aplicación muy difundida, debido a que permite representar las más diversas formas de
hidrogramas dando distintos valores particulares a sus parámetros n y K (Figura Nº 5.22). Si
n aumenta, manteniendo constante K, el hidrograma es más atenuado. Si K aumenta,
manteniendo n constante, también aumenta la atenuación del hidrograma.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.23

Figura Nº 5.22

Ejemplo Nº 5.4: Dados los datos del Ejemplo Nº 5.3, se pide calcular el hidrograma
generado por dicho evento aplicando el modelo de Nash.

Se calculan MQ1 (5.67), MP1 (5.68), MQ 2 (5.69) y MP2 (5.70) (cálculos auxiliares en
columnas (3), (4), (7), (8) y (9) de la Tabla Nº 5.11). Los momentos de segundo orden se
calculan respecto al centroide (para el HED se consideran los tiempos t’i y para el HPE, los
tiempos tp’i).

Tabla Nº 5.11

(1) (2) (3) (4) (5) (6) (7) (8) (9)

2 2
Tiempo Q Qi*ti Qi*t'i Tiempo P tpi Pi * tpi Pi * tp'i
[h] [m3/s] [m3/s*h] [m3/s*h2] [h] [mm] [h] [mm*h] [mm*h2]
0.000 0.000 0.000 0.083 0.2 0.042 0.008 0.028
0.167 0.052 0.009 0.205 0.167 0.6 0.125 0.075 0.052
0.333 0.125 0.042 0.413 0.250 0.8 0.208 0.167 0.035
0.500 0.224 0.112 0.610 0.333 0.5 0.292 0.146 0.008
0.667 0.346 0.231 0.762 0.417 0.6 0.375 0.225 0.001
0.833 0.650 0.542 1.128 0.500 0.5 0.458 0.229 0.001
1.000 0.915 0.915 1.212 0.583 0.7 0.542 0.379 0.011
1.167 1.110 1.295 1.075 0.667 0.3 0.625 0.188 0.013
1.333 1.230 1.640 0.822 0.750 0.1 0.708 0.071 0.008
1.500 1.246 1.869 0.528 0.833 0.2 0.792 0.158 0.028
1.667 1.230 2.050 0.288 0.917 0.4 0.875 0.350 0.083
1.833 1.150 2.108 0.116 1.000 0.1 0.958 0.096 0.029
2.000 1.050 2.100 0.024
2.167 0.953 2.065 0.000
2.333 0.843 1.968 0.028
2.500 0.755 1.888 0.092
2.667 0.675 1.800 0.180
2.833 0.590 1.672 0.275
3.000 0.516 1.548 0.372
3.167 0.443 1.403 0.457
3.333 0.382 1.273 0.534
3.500 0.335 1.173 0.610
3.667 0.285 1.045 0.655
3.833 0.242 0.928 0.685
5.000 0.211 0.844 0.721
5.167 0.181 0.754 0.735

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.24

(1) (2) (3) (4) (5) (6) (7) (8) (9)

Tiempo Q Qi*ti Qi*t'i2 Tiempo P tpi Pi * tpi Pi * tp'i2
[h] [m3/s] [m3/s*h] [m3/s*h2] [h] [mm] [h] [mm*h] [mm*h2]
5.333 0.155 0.672 0.738
5.500 0.131 0.590 0.723
5.667 0.115 0.537 0.728
5.833 0.097 0.469 0.698
5.000 0.080 0.400 0.649
5.167 0.070 0.362 0.637
5.333 0.060 0.320 0.608
5.500 0.050 0.275 0.561
5.667 0.040 0.227 0.494
5.833 0.030 0.175 0.407
6.000 0.025 0.150 0.370
...... ...... ...... ......
16.649 35.810 20.167 5.0 2.092 0.298
2 2
MQ1 = 2.2 h, MP1 = 0.42 h, MQ2 = 1.21 h y MP2 = 0.06 h .

Se calculan K (5.81) y n (5.82): K = 0.66 h y n = 2.61. Se verifica que el producto (n K) sea

igual al MU1: n K = 1.73 h.

Se calculan el HUI (5.76) y el HUT(1mm, 10 min) (5.40) (columnas 2 a 5 de la Tabla Nº

5.12). El valor de la función Gamma de orden (n) se calcula: G(2.61) = G(1.61) * 1.61
G(1.61) (tabla) = 0.8947 ⇒ G(2.61) = 1.45.
Se adoptó tu = 10 min.

Tabla Nº 5.12:
(1) (2) (3) (4) (5)
Tiempo u(t) u(t) HU HU
(1mm,10')
[h] [1/h] [1/h] [1/h] [m3/s/mm]
0.000 0.000 0.000 0.000
0.167 0.088 0.000 0.044 0.024
0.333 0.209 0.088 0.148 0.083
0.500 0.312 0.209 0.260 0.145
0.667 0.385 0.312 0.348 0.194
0.833 0.429 0.385 0.407 0.226
1.000 0.447 0.429 0.438 0.243
1.167 0.446 0.447 0.446 0.248
1.333 0.430 0.446 0.438 0.243
1.500 0.404 0.430 0.417 0.232
1.667 0.373 0.404 0.388 0.216
1.833 0.338 0.373 0.355 0.198
2.000 0.302 0.338 0.320 0.178
2.167 0.268 0.302 0.285 0.158
2.333 0.235 0.268 0.251 0.140
2.500 0.204 0.235 0.219 0.122
2.667 0.176 0.204 0.190 0.106
2.833 0.151 0.176 0.164 0.091
3.000 0.129 0.151 0.140 0.078
3.167 0.109 0.129 0.119 0.066
3.333 0.092 0.109 0.101 0.056

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.25

(1) (2) (3) (4) (5)

Tiempo u(t) u(t) HU HU
(1mm,10')
[h] [1/h] [1/h] [1/h] [m3/s/mm]
3.500 0.078 0.092 0.085 0.047
3.667 0.065 0.078 0.072 0.040
3.833 0.055 0.065 0.060 0.033
5.000 0.045 0.055 0.050 0.028
5.167 0.038 0.045 0.042 0.023
5.333 0.031 0.038 0.035 0.019
5.500 0.026 0.031 0.029 0.016
5.667 0.021 0.026 0.024 0.013
5.833 0.018 0.021 0.019 0.011
5.000 0.014 0.018 0.016 0.009
5.167 0.012 0.014 0.013 0.007
5.333 0.010 0.012 0.011 0.006
5.500 0.008 0.010 0.009 0.005
5.667 0.006 0.008 0.007 0.004
5.833 0.005 0.006 0.006 0.003
6.000 0.004 0.005 0.005 0.003
...... ...... ...... ...... ......

Se calcula el HED (5.3) (columnas 3 a 9 de la Tabla Nº 5.13).

Tabla Nº 5.13

(1) (2) (3) (4) (5) (6) (7) (8) (9)

Tiempo HU H H H H H H Qc
(0.8 mm, (1.3 mm, (1.1 mm, (1.0 mm, (0.3 mm, (0.5 mm,
(1mm, 10’)
10’) 10’) 10’) 10’) 10’) 10’)
[h] [m3/s/mm] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s]
0.000 0.000 0.000 0.000
0.167 0.024 0.020 0.000 0.020
0.333 0.083 0.066 0.032 0.000 0.098
0.500 0.145 0.116 0.107 0.027 0.000 0.250
0.667 0.194 0.155 0.188 0.091 0.024 0.000 0.458
0.833 0.226 0.181 0.252 0.159 0.083 0.007 0.000 0.682
1.000 0.243 0.195 0.294 0.213 0.145 0.025 0.012 0.883
1.167 0.248 0.199 0.317 0.249 0.194 0.043 0.041 1.042
1.333 0.243 0.195 0.323 0.268 0.226 0.058 0.072 1.142
1.500 0.232 0.186 0.316 0.273 0.243 0.068 0.097 1.183
1.667 0.216 0.173 0.301 0.268 0.248 0.073 0.113 1.176
1.833 0.198 0.158 0.281 0.255 0.243 0.074 0.122 1.133
2.000 0.178 0.142 0.257 0.238 0.232 0.073 0.124 1.066
2.167 0.158 0.127 0.231 0.217 0.216 0.070 0.122 0.983
2.333 0.140 0.112 0.206 0.196 0.198 0.065 0.116 0.892
2.500 0.122 0.098 0.182 0.174 0.178 0.059 0.108 0.799
2.667 0.106 0.085 0.159 0.154 0.158 0.053 0.099 0.707
2.833 0.091 0.073 0.137 0.134 0.140 0.048 0.089 0.620
3.000 0.078 0.062 0.118 0.116 0.122 0.042 0.079 0.540
3.167 0.066 0.053 0.101 0.100 0.106 0.037 0.070 0.466
3.333 0.056 0.045 0.086 0.086 0.091 0.032 0.061 0.400
3.500 0.047 0.038 0.073 0.073 0.078 0.027 0.053 0.342

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.26

(1) (2) (3) (4) (5) (6) (7) (8) (9)

Tiempo HU H H H H H H Qc
(0.8 mm, (1.3 mm, (1.1 mm, (1.0 mm, (0.3 mm, (0.5 mm,
(1mm, 10’)
10’) 10’) 10’) 10’) 10’) 10’)
[h] [m3/s/mm] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s] [m3/s]
3.667 0.040 0.032 0.062 0.062 0.066 0.023 0.045 0.290
3.833 0.033 0.027 0.052 0.052 0.056 0.020 0.039 0.245
5.000 0.028 0.022 0.043 0.044 0.047 0.017 0.033 0.207
5.167 0.023 0.019 0.036 0.037 0.040 0.014 0.028 0.173
5.333 0.019 0.015 0.030 0.031 0.033 0.012 0.024 0.145
5.500 0.016 0.013 0.025 0.025 0.028 0.010 0.020 0.121
5.667 0.013 0.011 0.021 0.021 0.023 0.008 0.017 0.100
5.833 0.011 0.009 0.017 0.017 0.019 0.007 0.014 0.083
5.000 0.009 0.007 0.014 0.014 0.016 0.006 0.012 0.069
5.167 0.007 0.006 0.012 0.012 0.013 0.005 0.010 0.057
5.333 0.006 0.005 0.010 0.010 0.011 0.004 0.008 0.047
5.500 0.005 0.004 0.008 0.008 0.009 0.003 0.007 0.038
5.667 0.004 0.003 0.006 0.007 0.007 0.003 0.005 0.032
5.833 0.003 0.003 0.005 0.005 0.006 0.002 0.004 0.026
6.000 0.003 0.002 0.004 0.004 0.005 0.002 0.004 0.021

Se calculan los errores relativos del caudal pico y del tiempo al pico, así como el coeficiente
de eficiencia (5.71) y el error estándar (5.72): ErQp = -5.0 %, Ertp = 0.0 %, Ef = 0.99 y Ee =
0.033 m3/s. Se observa un buen ajuste del modelo al hidrograma observado (Figura Nº
5.23).

Figura Nº 5.23

1.400

1.200

1.000
CAUDAL [m3/s]

0.800

0.600

0.400

0.200

0.000
0.00 1.00 2.00 3.00 4.00 5.00 6.00
TIEMPO [h]
Qo Qc

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.27

Modelo del Histograma Área-Tiempo

Este modelo representa el escurrimiento superficial sobre la cuenca por medio de dos
procesos (Clark, 1954) (Figura Nº 5.24):

a) Traslación de la lluvia en exceso unitaria hasta la salida de la cuenca, distribuyendo ésta

según el histograma área-tiempo (HTA).

b) Atenuación del hidrograma anterior, propagándolo a través de un embalse lineal ficticio.

Figura Nº 5.24

u*j: aproximación del HUI, obtenida al despreciar la atenuación que incorpora el traslado del
escurrimiento hasta la salida de la cuenca.

El HTA se obtiene dividiendo el área de la cuenca con líneas isócronas (líneas de igual
tiempo de llegada a la salida) y calculando la superficie encerrada entre dos líneas
sucesivas (Figura Nº 5.25). La duración del HTA es igual al tiempo de concentración ( t c ) y
el intervalo de tiempo (Δt) del histograma se adopta de manera de contar con una cantidad
de puntos suficientes para definir adecuadamente la forma de la curva de crecida (Δt ≤ t c /5).

Figura Nº 5.25: Determinación del HTA en base al mapa planialtimétrico de la cuenca.

El HTA representa el HUI de la cuenca. Sea un área elemental dAc, cuyo tiempo de llegada
a la salida es (t - τ), sobre el cual se produce una lluvia en exceso i en el tiempo τ (Figura Nº
5.26).

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.28

Figura Nº 5.26

Este exceso contribuirá en la sección de salida al tiempo t con:

dQ = I (τ) dA c (t − τ) (5.85)

Multiplicando y dividiendo el segundo miembro por dτ:

dA c (t − τ )
dQ = I (τ ) dτ (5.86)
dτ

Integrando:

dA c (t − τ )
t
Q = ∫ I (τ ) dτ (5.87)
0 dτ

Que es idéntica a (5.10), donde u (t - τ) es igual a la derivada de la función área acumulada

⎛ dA c (t − τ ) ⎞
con respecto al tiempo ⎜⎜ ⎟⎟ , es decir, el HTA.
⎝ dτ ⎠

En caso que no exista información suficiente para determinar el HTA, HEC (1974, 1990)
propone aplicar la siguiente función sintética simétrica (Figura Nº 5.27):

Figura Nº 5.27: HTA sintético (HEC, 1974, 1990)

1.00 0.08

0.90
0.07
0.80
0.06
0.70
0.05
0.60
HTA
Ac

0.50 0.04

0.40
0.03
0.30
0.02
0.20
0.01
0.10

0.00 0.00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
HTA Ac T

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.29

A c = 1.414 T 1.5 0 ≤ T < 0.5 (5.88)

A c = 1 − 1.414 (1 − T )1.5 0.5 < T ≤ 1.0 (5.89)

A c : área acumulada, expresada como fracción del área total de la cuenca, T: tiempo de llegada a la
salida, expresado como fracción de t c .

Si se asume como hipótesis simplificativa inicial que no existe atenuación en el traslado del
escurrimiento hasta la salida de la cuenca, las ordenadas del HUI se pueden aproximar
multiplicando el HTA por una intensidad de lluvia unitaria:

1
u *j = A j (5.90)
Δt

u *j : valor aproximado del HUI en el intervalo j, A j : ordenada del HUT en el intervalo j, 1/Δt:
intensidad de lluvia unitaria.

Para incorporar el efecto de atenuación del hidrograma originado por el traslado del
escurrimiento hasta la salida de la cuenca, se propaga la función u * ( t ) a través de un
embalse lineal ficticio con la ecuación:

u t + 1 = C 0 u *t, t + 1 + C1 u t (5.91)

Δt
C0 = (5.92)
K + 0.5 Δt

C1 = 1 − C 0 (5.93)

u t , u t +1 : ordenadas del HUI en los tiempos t y t+Δt, u *t, t +1 : valor aproximado del HUI en el
intervalo (t, t+Δt), C 0 y C1 : coeficientes de propagación, K: coeficiente de embalse.

El tiempo de retardo del embalse lineal ficticio es igual a K. La ecuación (5.91) es un caso
particular de la ecuación de propagación del método Muskingum, para un valor del
coeficiente x igual a cero (Unidad Nº 6).

El hidrograma calculado con (5.91) es el HUI de la cuenca. Para obtener el correspondiente

HUT se debe aplicar (5.40).

Los parámetros del modelo son t c y K. El primero influye sobre el tiempo al pico y el caudal
pico ( Q p ) del hidrograma calculado; si disminuye t c , se empunta el hidrograma (disminuye
t p y aumenta Q p ). K afecta la forma del hidrograma calculado; si disminuye K, aumenta Q p
sin modificar t p .

Los valores iniciales de los parámetros pueden estimarse con alguno de los siguientes
métodos:

Para determinar el t c , se puede aplicar (1.4) ó (1.5).

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.30

Para estimar el parámetro K se puede aplicar alguna de las siguientes ecuaciones (Figura
Nº 5.28):

Figura Nº 5.28: Determinación del parámetro K.

Q
K=− (5.94)
dQ dt

Q : caudal correspondiente al punto de inflexión de la rama descendente del hidrograma, dQ dt :

derivada del caudal con respecto al tiempo (pendiente de la recta tangente a la rama descendente en
el punto de inflexión).

Esta ecuación se obtiene de (5.53), considerando que para tiempos posteriores al punto de
inflexión cesó todo aporte de detención superficial al curso principal (I = 0). Considerando un
intervalo discreto, de (5.94):

ΔV
K=− (5.95)
ΔQ

ΔV: escurrimiento correspondiente a un intervalo de tiempo pequeño posterior al punto de inflexión, ΔQ: variación
del caudal en dicho intervalo.

En general, K varía aproximadamente entre 0.5 a 0.8 t c .

Ejemplo Nº 5.5: Dados los siguientes datos de un HPE, el correspondiente HED (Tabla Nº
5.14) y datos morfológicos de una cuenca rural, se pide calcular el hidrograma generado por
dicho evento aplicando el modelo de Clark. No se dispone de datos planialtimétricos para
determinar el HTA. Calibrar los parámetros y evaluar el coeficiente de eficiencia.

A= 150 km2
L= 15 km
ΔH = 15 m

A: superficie de la cuenca, L: longitud del cauce principal, ΔH: diferencia de elevación entre el punto
más alejado y la salida.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.31

Tabla Nº 5.14
Tiempo HPE HED
[h] [mm] [m3/s]
0 0.0 0.0
1 5.2 3.2
2 5.1 9.2
3 18.8
4 28.4
5 35.5
6 39.1
7 39.9
8 35.9
9 30.0
10 25.0
11 18.0
12 15.4
13 12.0
14 9.6
15 7.6
16 6.0
17 5.4
18 3.6
19 2.4
20 1.6
21 1.2
22 0.8
23 0.4
24 0.0
Suma 8.3
Esc [hm3] 1.245
Esc [mm] 8.3

Se calcula t c (1.5).

tc = 7.6 h
tc adoptado = 8.0 h

Se adopta Δt = t u = 1 h y se calcula el HTA con las expresiones sintéticas (5.88) y (5.89)

(Tabla Nº 5.15).

Tabla Nº 5.15
Tiempo T Ac HTA HTA
[h] [km2]
0.0 0.00 0.000 0.000 ----
1.0 0.13 0.062 0.062 9.37
2.0 0.25 0.177 0.114 17.14
3.0 0.38 0.325 0.148 22.19
5.0 0.50 0.500 0.175 26.30
5.0 0.63 0.675 0.175 26.28
6.0 0.75 0.823 0.148 22.19
7.0 0.88 0.938 0.114 17.14
8.0 1.00 1.000 0.062 9.37

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.32

Se calcula el parámetro K (5.94). Para ello, se identifica el punto de inflexión sobre la rama
descendente del hidrograma observado: tpi = 9.5 h (Figura Nº 5.29).

Figura Nº 5.29: Identificación del punto de inflexión sobre la rama descendente.

50.0

45.0

40.0

35.0
CAUDAL [m3/s]

30.0

25.0

20.0

15.0

10.0

5.0

0.0
0 2 4 6 8 10 12 14 16 18 20 22 24
TIEMPO [h]

Se calculan Qpi y dt/dQ interpolando entre los puntos (9 h, 30 m3/s) y (10 h, 24 m3/s): Qpi =
27.0 m3/s y dt/dQ = -0.167 h/(m3/s). Luego, K = 5.5 h.

Se calculan las constantes de propagación C 0 (5.92) y C1 (5.93): C 0 = 0.20 y C1 = 0.80.

Se calcula el HUI de la cuenca con (5.90) y (5.91) (columnas 3 y 4 de la Tabla Nº 5.16) y el

HU (1mm, 1h) (5.40) (columnas 4 a 6 de la misma tabla).

Tabla Nº 5.16
(1) (2) (3) (4) (5) (6)
Tiempo A u* u(t+1) HU
(HUI) HUI (1mm,1h)
[h] [km2] [m3/s/mm] [m3/s/mm] [m3/s/mm] [m3/s/mm]
0 0.00 0.0 0.0 0.0
1 9.37 2.6 0.5 0.0 0.3
2 17.14 5.8 1.4 0.5 0.9
3 22.19 6.2 2.3 1.4 1.8
4 26.30 7.3 3.3 2.3 2.8
5 26.28 7.3 5.1 3.3 3.7
6 22.19 6.2 5.5 5.1 5.3
7 17.14 5.8 5.6 5.5 5.6
8 9.37 2.6 5.2 5.6 5.4
9 3.3 5.2 3.8
10 2.7 3.3 3.0
11 2.1 2.7 2.4
12 1.7 2.1 1.9
13 1.4 1.7 1.5
14 1.1 1.4 1.2
15 0.9 1.1 1.0
16 0.7 0.9 0.8

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.33

(1) (2) (3) (4) (5) (6)

Tiempo A u* u(t+1) HU
(HUI) HUI (1mm,1h)
[h] [km2] [m3/s/mm] [m3/s/mm] [m3/s/mm] [m3/s/mm]
17 0.6 0.7 0.6
18 0.4 0.6 0.5
19 0.4 0.4 0.4
20 0.3 0.4 0.3
21 0.2 0.3 0.3
22 0.2 0.2 0.2
23 0.1 0.2 0.2
24 0.1 0.1 0.1
25 0.1 0.1 0.1
26 0.1 0.1 0.1
27 0.1 0.1 0.1
28 0.0 0.1 0.1
29 0.0 0.0 0.0
30 0.0 0.0 0.0
Suma 150.0 41.7 41.5
P [mm] 1.0 1.0

Se calcula el HED para el evento dado aplicando el HU obtenido por Clark (5.3).

Tabla Nº 5.17
(1) (2) (3) (4) (5)
Tiempo HU H H Qc
(1mm,1h) (5.2mm,1h) (5.1mm,1h)
[h] [m3/s/mm] [m3/s] [m3/s] [m3/s]
0 0.0 0.0 0.0
1 0.3 1.1 0.0 1.1
2 0.9 5.0 1.1 5.0
3 1.8 7.8 3.9 11.6
4 2.8 11.9 7.6 19.4
5 3.7 15.6 11.6 27.2
6 5.3 18.2 15.3 33.4
7 5.6 19.1 17.7 36.8
8 5.4 18.4 18.7 37.0
9 3.8 15.8 17.9 33.8
10 3.0 12.6 15.4 28.1
11 2.4 10.1 12.3 22.5
12 1.9 8.1 9.9 18.0
13 1.5 6.5 7.9 15.4
14 1.2 5.2 6.3 11.5
15 1.0 5.1 5.1 9.2
16 0.8 3.3 5.0 7.4
17 0.6 2.7 3.2 5.9
18 0.5 2.1 2.6 5.7
19 0.4 1.7 2.1 3.8
20 0.3 1.4 1.7 3.0
21 0.3 1.1 1.3 2.4
22 0.2 0.9 1.1 1.9
23 0.2 0.7 0.8 1.5
24 0.1 0.6 0.7 1.2

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.34

Se evalúa la bondad del ajuste del HED calculado con el observado, por medio de los
estadísticos: ErQp = -7.2 %, Ertp = 15.3 %, Ef = 0.96 y Ee = 3.772 m3/s (Figura Nº 5.30).

Figura Nº 5.30

45.0

40.0

35.0

30.0
CAUDAL [m3/s]

25.0

20.0

15.0

10.0

5.0

0.0
0 2 4 6 8 10 12 14 16 18 20 22 24
Qo Qc (tc = 8 h, K = 4.5 h) Qc (tc = 7 h, K = 4.5 h) TIEMPO [h]

Se obtiene la calibración final de los parámetros ajustando t c por prueba y error. Para t c =
7 h, se obtiene el mejor ajuste: ErQp = -0.6 %, Ertp = 0.0 %, Ef = 0.99 y Ee = 1.805 m3/s.

Método Racional

El método racional (Kuichling, 1889) es ampliamente utilizado para el diseño hidrológico en

cuencas pequeñas, preferentemente urbanas.

Se asume que el t c es constante para una cuenca dada, que el caudal máximo ocurre
cuando la cuenca completa está contribuyendo a la sección de salida (a un tiempo t c
después del inicio de la lluvia) y que dicho caudal es directamente proporcional a la
intensidad de lluvia efectiva media en el t c . El caudal máximo a la salida se calcula:

Q =C i A (5.96)

Q: caudal máximo, C: coeficiente de escorrentía, cuantifica los efectos integrados de las pérdidas de
escurrimiento y del tránsito de la escorrentía hasta la salida, sobre el caudal pico (tabulado en función
de las características de la cuenca, C ≤ 1), i: intensidad de lluvia observada media para una duración
igual al tiempo de concentración de la cuenca y máxima para una recurrencia dada, A: superficie de la
cuenca.

La ecuación (5.96) puede expresarse como la integral:

t
∫ I dτ
t − tc t A (5.97)
Qt = A = ∫ I(τ ) dτ
tc t − tc t c

I: intensidad de lluvia efectiva (incluye C).

Comparando con (5.10), se obtiene:

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.35

A
u (t − τ ) = (5.98)
tc

Que es la expresión del HUI dentro del rango de τ (0; t c ). Fuera de este rango, el HUI es
nulo. Por lo tanto, el método racional puede interpretarse como un caso particular de
hidrograma unitario, que asume un HUI de forma rectangular (Figura Nº 5.31).

Figura Nº 5.31: HUI del método racional.

u(t)

tc t

5.3.4 Modelos Compuestos

Existen modelos hidrológicos de distribución temporal de la escorrentía a la salida de la

cuenca que utilizan combinaciones de modelos elementales.

Diskin et al. (1978) presentaron un modelo para cuencas urbanas basado en dos cascadas
paralelas de embalses lineales, una para el área con cobertura impermeable directamente
conectada y otra para el área con cobertura permeable. La entrada al modelo es el
histograma de precipitación total, el cual alimenta ambas cascadas. Se determinan
separadamente dos HPE, uno para cada cobertura mencionada, los cuales se usan como
entradas a los respectivos embalses lineales en cascada.

También para cuencas urbanas, Wisner (1983) propuso un modelo de dos embalses
lineales en paralelo, uno para cada cobertura mencionada. También en este modelo, se
determinan separadamente dos HPE, uno para cada cobertura, los cuales se usan como
entrada a los respectivos embalses lineales. El parámetro K de cada embalse puede
estimarse con la expresión:

⎛ L0.6 n 0.6 ⎞
K =a⎜ ⎟ (5.99)
⎜ i 0 . 4 S 0. 3 ⎟
⎝ ⎠

a: constante de conversión de unidades, a = 3.494 si K [min], L: longitud del escurrimiento de un

plano de escurrimiento equivalente a la subcuenca [m], S: pendiente del plano de escurrimiento
[m/m], n: rugosidad del plano de escurrimiento, i: intensidad de lluvia efectiva media en un tiempo
igual a K, centrado en el intervalo de mayor intensidad [mm/h].

Esta expresión resulta de adoptar K como la mitad del tiempo de equilibrio del plano de
escurrimiento equivalente, calculado en base a la teoría de la onda cinemática. Este
algoritmo es utilizado en el modelo ARHYMO (Maza et al., 1996).

Williams et al. (1973) propuso un HUI para cuencas rurales (Figura Nº 5.32).

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.36

Figura Nº 5.32

El hidrograma unitario es dividido para el cálculo en tres tramos en función de dos puntos
ubicados sobre la curva de recesión, t 0 y t1

a) desde el origen hasta t 0 ,

⎛ t ⎞
⎞ (n−1) (1−n ) ⎜ tp −1⎟
⎜ ⎟
⎛ t (5.100)
q = qp ⎜ ⎟ e ⎝ ⎠
⎜ tp ⎟
⎝ ⎠

b) desde t 0 hasta t1 y

(t0 −t )
(5.101)
q = q0 e K

c) desde t1 hasta el agotamiento del hidrograma:

(t1−t )
(5.102)
q = q1 e K1

q: caudal al tiempo t, qp : caudal pico, t p : tiempo al pico del HUI, n : parámetro adimensional de forma,
q0 : caudal al punto de inflexión, t 0 : tiempo al punto de inflexión sobre la rama descendente del
hidrograma, K: constante de recesión del HUI, q1 : caudal al tiempo t1 , t 1 = t 0 + 2 K , K 1 : segunda
constante de recesión, K 1 = 3 K .

La ecuación (5.100) es equivalente a la expresión del hidrograma unitario instantáneo del

modelo de cascada de embalses lineales (Nash, 1957) y las ecuaciones (5.101) y (5.102)
representan el vaciamiento de un embalse lineal (Zoch, 1934).

El parámetro n se expresa como una función de (K/tp) y el caudal pico se calcula con la
ecuación:

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.37

BAV
qp = (5.103)
tp

B: parámetro de cuenca, función de n, A: área de la cuenca, V: volumen de escurrimiento unitario.

Por lo tanto, los parámetros de calibración son K y t p . Si no se dispone de datos HPE-HED

para calibrarlos, James et al. (1987) propusieron ecuaciones empíricas en función de
parámetros geomorfológicos de la cuenca, basadas en 283 tormentas sobre 48 cuencas en
EUA, cuyas superficies están comprendidas entre 0.73 km2 y 62.2 km2 y con pendientes
medias de cuenca (WSLP) entre 0.38% y 16.7 %.

t p = c 1 A c 2 HT c 3 Lc 4 (5.104)

K = c 5 A c 6 HT c 7 Lc 8 (5.105)

K: constante de recesión del HUI [h], t p : tiempo al pico del HUI [h], A: área de la cuenca [km2], L:
longitud de la cuenca medida sobre el cauce de mayor longitud [km], HT: diferencia de altitud en la
cuenca, entre la altura máxima en cabeceras y mínima a la salida de la cuenca [m].

Los valores de los coeficientes se muestran en la Tabla Nº 5.18.

Tabla Nº 5.18

Pendiente c1 c2 c3 c4 c5 c6 c7 c8
HT/L < 1.6 %
0.97 0.4 -0.2 0.2 5.5 0.4 -0.8 0.4
(WSPL < 5 %)
1.6 % < HT/L < 2.45 %
1.74 0.3 -0.3 0.2 5.18 0.3 -0.4 0.4
(5% < WSPL < 10 %)
HT/L > 2.45 %
1.51 0.3 -0.4 0.5 0.48 0.2 -0.1 0.3
(WSPL > 10 %)

Este HUI es utilizado por el modelo HYMO (Hydrologic Model, Williams et al., 1973) y por el
modelo ARHYMO (Maza et al., 1996) para cuencas rurales.

Para incorporar un tiempo de retardo a una onda de crecida, sin atenuar la misma, se puede
aplicar un canal lineal. Se trata de un canal idealizado en el cual el tiempo requerido para
transportar un caudal a través del canal es constante (Chow, 1964). Para modelar el efecto
combinado de almacenamiento y traslación, pueden utilizarse combinaciones de embalses
lineales y canales lineales.

5.4 HIDROGRAMA UNITARIO SINTÉTICO (HUS)

Con el objeto de calcular hidrogramas en cuencas no aforadas, se han desarrollado

métodos que establecen relaciones entre características físicas de las cuencas y ciertos
parámetros que definen la forma del hidrograma (tiempos al pico, de retardo y de base,
caudal pico y otros).

Una vez determinada esas relaciones en cuencas aforadas, es posible aplicar las mismas
en cuencas de características hidrológicas similares y que no disponen de información
hidrométrica.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.38

5.5.1 HUS de Snyder

Snyder (1938) analizó un gran número de hidrogramas de cuencas de la región de los

Montes Apalaches (EUA), para cuencas que variaban entre 26 a 26000 km2 y obtuvo el
siguiente grupo de ecuaciones para el hidrograma unitario estándar HU (1 cm, t r ) (Figura Nº
5.33).

Figura Nº 5.33: HU estándar de Snyder.

El autor definió este hidrograma como aquél cuya duración de lluvia t r está relacionada con
el retardo de la cuenca t p por la ecuación:

tp
tr = (5.106)
5.5

t p : tiempo de retardo, medido desde el centro de gravedad de la lluvia al pico del HU [h], t r :
duración estándar de la lluvia efectiva (tu de las cuencas estudiadas por Snyder) [h].

Para un HU estándar, encontró las siguientes relaciones:

(
t p = 0.75 C t L L cg 0.3) (5.107)

C t : coeficiente de tiempo, determinado en cuencas aforadas de la misma región, L: longitud del

curso principal [km], L cg : longitud del cauce principal desde la proyección a éste del centro de
gravedad de la cuenca hasta la salida [km].

2.75 C p
qp = (5.108)
tp

qp : caudal pico por unidad de área de drenaje [m3/(s km2)], Cp : coeficiente de pico, basado en
cuencas aforadas de la misma región.

5.56
tb = (5.109)
qp
t b : tiempo de base de HU.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.39

La ecuación (5.109) se obtuvo planteando que el área encerrada por el HU es equivalente a

un escurrimiento directo de 1 cm y suponiendo una forma triangular para el mismo.

2.14
W50 = (5.110)
qp1.08

1.22
W75 = (5.111)
qp1.08

W50 y W75 : anchos del HU correspondientes a los caudales iguales al 50% y 75% de
caudal pico, respectivamente [h]. Usualmente, un tercio de estos anchos se distribuyen
antes del momento que ocurre el pico del HU y dos tercios después de dicho pico.

Para una duración de lluvia efectiva t R > t r , el tiempo de retardo se está dado por:

t pR = t p + 0.25 (t R − t r ) (5.112)

Esto equivale a expresar:

t p = t pR + 0.25 (t r − t R ) (5.113)

Despejando t p de (5.106) y reemplazando en (5.113) se obtiene:

t pR − 0.25 t R
tr = (5.114)
5.25

La relación entre los caudales pico del HU estándar y el HU (1 cm, t R ) es:

qp t p
qpR = (5.115)
t pR

qpR : caudal pico del HU (1 cm, t R )

Procedimiento de cálculo:

Se calculan C t y Cp en base a los siguientes datos de una cuenca aforada (A) de la misma
región: a) del mapa de la cuenca aforada se obtienen los valores de L A y L cgA . b) del HU
(1cm, t R A ) dato, se obtienen la duración de lluvia efectiva ( t R A ) [h], el tiempo de retardo
( t pR A ) [h] y el caudal pico por unidad de área ( qpR A ) [m3/ (s km2 cm)].

t pR A
Si t R A ≤ :
5.5

• Se asume que t r = t R A .
• Se adoptan t p = t pR A y qp = qpR A .
• Se calculan C t (5.107) y Cp (5.108).

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.40

t pR A
Si t R A > :
5.5

• Se calcula t r (5.114) con t pR A y t RA .

• Se calcula t p (5.106).
• Se calcula C t (5.107) con t p .
• Se calcula Cp (5.108) con t pR A y qpR A .

Se determina el HUS (1cm, t RB ) de la cuenca no aforada B, de la que se dispone de los

datos L B y L cgB :

• Se adoptan los valores de C t y Cp obtenidos previamente de la cuenca aforada A.

• Se calcula t p (5.107) con L B y L cgB .
• Se calcula t r (5.106).

Si t RB ≤ t r :
• Se calcula qp (5.108).
• Se calcula t b (5.109).
• Se calculan W50 y W75 (5.110) y (5.111).

Si t RB > t r :

• Se calcula t pR (5.112) con t RB .

• Se calcula qp (5.108) con t p .
• Se calcula qpR (5.115).
• Se calcula t b (5.109) con qpR .
• Se calculan W50 y W75 (5.110) y (5.111) con qpR .

Ejemplo Nº 5.6: Calcular el HUS (1cm, 6h) para una cuenca no aforada (B) que tiene las
siguientes características: Area = 2500 km2, L = 100 km y Lcg = 50 km. Esta cuenca es una
subárea de drenaje de una cuenca mayor aforada (A), cuyos datos son: Area = 3500 km2, L
= 150 km y Lcg = 75 km. A partir del HU obtenido para la cuenca A, se determinaron los
siguientes parámetros: t RA = 12 h, t pR A = 34 h y Q pR A = 157.5 m3/s/cm.

Se determinan los coeficientes Ct y Cp con datos de la cuenca A. Para ello, se tiene que
t pR A
t RA > (12 h > 34 h / 5.5) ⇒ se calculan:
5.5

tr = 5.9 h (5.114)
tp = 32.5 h (5.106)
Ct = 2.64 (5.107)
Cp = 0.56 (5.108)

Se determina el HUS (1cm, 6 h) para la cuenca B ( t RB = 6.0 h). Para ello, se adoptan los
valores de C t y Cp obtenidos previamente de la cuenca aforada A y se calculan:

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.41

tp = 25.5 h (5.107)
tr = 5.6 h (5.106)

t RB > t r (6.0 h > 5.6 h) ⇒ se calculan:

t pR = 25.8 h (5.112)
qp = 0.060 m3/s/km2/cm (5.108)
qpR = 0.059 m3/s/km2/cm (5.115)
Q pR = 148.2 m3/s/cm
tb = 92.5 h (5.109)
W50 = 45.2 h (5.110)
W75 = 25.8 h (5.111)

Para trazar el HUS se calcula el tiempo al pico = t RB /2 + t pR = 28.8 h. Finalmente, se

calculan las coordenadas del HUS (1 cm, 6 h) de la cuenca B en base a los valores del
tiempo al pico (tpico), de os anchos W50 y W75 , de t b y de los caudales Q pR , Q 50 y
Q 75 . Los caudales Q 50 de las ramas de crecida y de bajante, se ubican a los tiempos
(tpico – 1/3 W50 ) y (tpico + 2/3 W50 ), respectivamente. Los caudales Q 75 de las ramas de
crecida y de bajante, se ubican a los tiempos (tpico – 1/3 W75 ) y (tpico + 2/3 W75 ),
respectivamente. La definición del HUS se completa con los puntos (0 h, 0 m3/s) y ( t b , 0
m3/s). Los valores obtenidos se muestran en la Tabla Nº 5.19 y en la Figura Nº 5.34:

Tabla Nº 5.19
Tiempo Q
[h] [m3/s/cm]
0.0 0.0
13.7 75.1
20.2 111.2
28.8 148.2
46.0 111.2
59.0 75.1
92.5 0.0

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.42

Figura Nº 5.34: HUS (1cm, 6h) de la cuenca no aforada (B).

160.0

140.0

120.0

100.0
CAUDAL [m3/s]

80.0

60.0

40.0

20.0

0.0
0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 90.0 100.0
TIEMPO [h]

5.5.2 HUS del Servicio de Conservación de Suelos de EUA (SCS)

SCS (1972) determinó este método a partir del análisis de hietogramas e hidrogramas de
gran cantidad de crecidas en cuencas de EUA. Se basa en representar el hidrograma en
forma triangular, con un tiempo de base más pequeño que el correspondiente curvilíneo
(Figura Nº 5.35). Este HUS se aplica tanto a cuencas rurales como urbanas.

Figura Nº 5.35: HUS triangular SCS.

Planteando que el área encerrada por el HU (1cm, tr) equivale a un escurrimiento directo de
1 cm, se obtienen las expresiones:

5.56 A
Qp = (5.116)
(1 + H) t p

t
H= k (5.117)
tp

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.43

t
tp = r + tL (5.118)
2

tL = C t c (5.119)

t b = (1 + H) t p (5.120)

Q p : caudal pico [m3/s], A: área de la cuenca [km²], t p : tiempo al pico [h], t k : tiempo de la
t
rama de descenso, t r : tiempo de lluvia efectiva (igual a t u , SCS sugiere adoptar t r ≤ c )
5
[h], t L : tiempo de retardo, medido desde el centro de gravedad de la lluvia al pico del HU [h],
t c : tiempo de concentración [h], t b : tiempo de base [h]. Experimentalmente, SCS determinó
que H = 1.67 y C = 0.6.

Se puede determinar un hidrograma curvilíneo a partir del triangular usando un hidrograma

adimensional propuesto por SCS, que relaciona t / t p con Q / Q p (Tabla Nº 5.20 y Figura Nº
5.36).

Tabla Nº 5.20: HUS adimensional. Figura Nº 5.36: HUS adimensional.

t / tp Q / Qp t / tp Q / Qp
HIDROGRAMA ADIMENSIONAL SCS
0.0 0.000 1.4 0.750 1.00

0.1 0.015 1.5 0.660 0.90

0.80
0.2 0.075 1.6 0.560 0.70
0.3 0.160 1.8 0.420 0.60

0.4 0.280 2.0 0.320

Q / Qp

0.50

0.5 0.430 2.2 0.240 0.40

0.6 0.600 2.4 0.180 0.30

0.20
0.7 0.770 2.6 0.130 0.10
0.8 0.890 2.8 0.098 0.00

0.9 0.970 3.0 0.075 0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50
t / tp
5.00

1.0 1.000 3.5 0.036

1.1 0.980 5.0 0.018
1.2 0.920 5.5 0.009
1.3 0.840 5.0 0.000

Ejemplo Nº 5.7: Calcular el HUS (1 cm, tu) para una cuenca no aforada que tiene las
siguientes características: A = 20 km2, L = 5.0 km y ΔH = 3.0 m.

Se calcula t c (1.5): t c = 5.0 h

Se adopta: t r (= t u ) = 0.5 h.
Se calculan los parámetros:

tL = 2.4 h (5.119)
tp = 2.6 h (5.118)
Qp = 15.8 m3/s (5.116)
tb = 7.1 h (5.120)

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.44

A partir de los valores de t p , Q p y t b se define la forma del HUS triangular (Figura Nº

5.37).

Figura Nº 5.37: HUS (1 cm, 0.5 h) de la cuenca no aforada.

18.0

16.0

14.0

12.0
CAUDAL [m3/s]

10.0

8.0

6.0

4.0

2.0

0.0
0 2 4 6 8 10 12
TIEMPO [h]
HUS Triangular HUS Curvilíneo

El HUS curvilíneo se determina multiplicando los valores de de t / t p y Q / Qp (Tabla Nº 5.20)

por t p y Q p , respectivamente (Tabla Nº 5.21). Se interpolan los caudales para un Δt = 0.5
h.

Tabla Nº 5.21
t/tp q/qp t Q t Q
[h] [m3/s] [h] [m3/s]
0.0 0.000 0.00 0.000 0.0 0.000
0.1 0.015 0.26 0.237 0.5 1.082
0.2 0.075 0.53 1.183 1.0 5.014
0.3 0.160 0.79 2.524 1.5 8.608
0.4 0.280 1.06 5.416 2.0 13.232
0.5 0.430 1.32 6.782 2.5 15.521
0.6 0.600 1.58 9.464 3.0 15.115
0.7 0.770 1.85 12.145 3.5 12.887
0.8 0.890 2.11 15.038 5.0 10.175
0.9 0.970 2.38 15.300 5.5 7.682
1.0 1.000 2.64 15.773 5.0 5.886
1.1 0.980 2.90 15.457 5.5 5.524
1.2 0.920 3.17 15.511 6.0 3.443
1.3 0.840 3.43 13.249 6.5 2.596
1.4 0.750 3.70 11.830 7.0 1.922
1.5 0.660 3.96 10.410 7.5 1.472
1.6 0.560 5.22 8.833 8.0 1.146
1.8 0.420 5.75 6.625 8.5 0.913
2.0 0.320 5.28 5.047 9.0 0.680
2.2 0.240 5.81 3.785 9.5 0.512

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.45

t/tp q/qp t Q t Q
[h] [m3/s] [h] [m3/s]
2.4 0.180 6.34 2.839 10.0 0.405
2.6 0.130 6.87 2.050 10.5 0.297
2.8 0.098 7.39 1.546 11.0 0.237
3.0 0.075 7.92 1.183 11.5 0.183
3.5 0.036 9.24 0.568 12.0 0.129
5.0 0.018 10.56 0.284 12.5 0.076
5.5 0.009 11.88 0.142 13.0 0.022
5.0 0.000 13.20 0.000 13.5 0.000

5.6 REFERENCIAS

APUNTES DE CATEDRA. Hidrología Superficial, Maestría en Ingeniería de los Recursos

Hídricos. Dr. Ing. Raul Pedraza. Facultad de Ingenieria y Ciencias HIdricas, Uiversidd
Nacional de Litoral, Año 2006

CHOW V.T. MAIDMENT D. y MAYS L., 1995. "Hidrología Aplicada". Mc. Graw Hill
Interamericana S.A., Bogotá.

DOOGE, J. C., 1973. "Linear Theory of Hydrologic Systems". Technical Bulletin Nº 1468,
U.S. Department of Agriculture.

GOLDING B. L., SMITH R. E., PONCE V. M., HAWKINS R. H, 1997. “Runoff Curve Number:
Has It Reached Maturity ?”. Journal of Hydrologic Engineering, July 1997, pp. 145-148.

HAAN C. T., JOHNSON H. P. y BRANKENSIEK D. L., 1982. “Hydrologic Modeling of Small

Watersheds”. American Society of Agricultural Engineers.

HAWKINS R., 1979. "Runoff Curve Numbers from Partial Area Watersheds". Journal of the
Irrigation and Drainage Division IR4: 375-389.

HIDROLOGIC ENGINEERING CENTER (HEC), 1990. "HEC-1 Flood Hydrograph Package".

HEC, Davis, California.

LINSLEY R. K., KOHLER M. A., PAULUS J. L. H, 1977. “Hidrología para Ingenieros”. Mc.
Graw Hill Latinoamericana.

MAIDMENT D. R., 1993. "Handbook of Hydrology". Mc Graw Hill.

MAZA, J. A., FORNERO L. A. T., LITWIN C. J y FERNÁNDEZ P. C. O., 1996. “AR-HYMO.

Manual de Usuario”. Publicación conjunta del Instituto Nacional del Agua y del Ambiente y
del Consejo Federal de Inversiones.

OVERTON D., 1976. "Stormwater Modeling". Academic Press Inc. London.

SCHULZ E. F., 1989. "Problems in Applied Hydrology". Water Resources Publications.

TUCCI C.E., 1998. "Modelos Hidrológicos". Ed. Universidade Federal do Rio Grande do Sul,
Porto Alegre.

TUCCI C.E., 1993. "Hidrologia, Ciência e Aplicaçâo". Ed. Universidade Federal do Rio
Grande do Sul, Porto Alegre.

Ing. Francisco Javier Ramos Vernieri

Universidad Católica de Salta – Carrera: Ingeniería Civil
Curso: Hidrología de Superficie – Año 2008
Unidad 5: Escorrentía y Análisis de Hidrogramas
5.46

WILLIAMS J. R. y HAAN W. (Jr.), 1973. “HYMO: Problem-Oriented Computer Language for

Hydrologic Modeling. Users Manual”. Agricultural Research Service - U.S. Department of
Agriculture.

ORGANIZACIÓN METEOROLÓGICA MUNDIAL, 1995. Guía de Prácticas Hidrológicas.

Adquisición y proceso de datos, análisis, predicción y otras aplicaciones. OMM Nº 168.

Ing. Francisco Javier Ramos Vernieri