MODELO DE GENERACION DE

DESCARGAS MEDIAS
MENSUALES

BALANCE HIDROLOGICO
DE LUTZ SCHOLZ

ING. ALBERTO QUESQUEN R

AGOSTO 2016
 MODELO GENERACION DE DESCARGAS
BALANCE HIDRICO - LUTZ SCHOLZ
Fischer, ecuación fundamental del BH en mm/mes: Cmi = Pi - Di + Gi - Ai
Lutz, señala: ABASTECIMIENTO Pi
Ai = ai (R/100)

RETENCION (R) : Gi - Ai

GASTO
Qt/Qo=e^(-a.t)
bo = e^(-a.t)

Cmi = Pi - Di + Gi - Ai
Ri = Gi - Ai
Turc: Ci=(Pi-Di)/Pi , Ci.Pi=(Pi- Di) ESCURRIMIENTO
USBR: PE Hidrológico es la antitésis de
la PE de los cultivos.
En situación de equilibrio:
Cmi = PEi + Ri
Cmi = PEi + Gi -Ai
Descripción de los Elementos
del Balance Hidrológico
Donde:

CMi = Caudal mensual (mm/mes)

Pi = Precipitación mensual sobre la cuenca (mm/mes)
Di = Déficit de escurrimiento (mm/mes)
Gi = Gasto de la retención de la cuenca (mm/mes)
Ai = Abastecimiento de la retención (mm/mes)
C = coeficiente de escurrimiento
Ri = Retención de la cuenca (mm/mes)
PE = Precipitación efectiva (mm/mes)
Qt = descarga en el tiempo t
Qo = descarga inicial
a = Coeficiente de agotamiento
t = tiempo
bo = Relación entre descarga mes actual y mes anterior y
que corresponde al primer mes de estiaje.
 ELEMENTOS DEL BALANCE HIDROLÓGICO
1.0 PRECIPITACION (P): Referida a la precipitación media de la cuenca, en mm

2.0 COEF. ESCURRIMIENTO (C): Adimencional

3.0 PRECIPITACION EFECTIVA (PE) : En el sentido hidrológico, es una parte de la
precipitación total mensual que produce escorrentía.
Utiliza la antitesis del método del USBR, es decir que de la relación de lluvia una porción es
aprovechada por las plantas y lo que queda produce escorrentía a la que llamó PE.

PRECIPITACION EFECTIVA
METODO DEL USBR
Precipitación Total Porcentaje del Precipitación Efectiva
Mensual Aumento Acumulada
(Rango) (Rango) (Rango)
mm % mm
0.0 - 25.4 90 - 100 22.9 - 25.4
25.4 - 50.8 85 - 95 44.5 - 49.5
50.8 - 76.2 75 - 90 63.5 - 72.4
76.2 - 101.6 50 - 80 76.2 - 92.7
101.6 - 127.0 30 - 60 83.8 - 102.9
127.0 - 152.4 10 - 40 86.4 - 118.1
> - 152.4 0 - 10 118.1 - 120.6
Fuente: Generación de Caudales para la Sierra del Perú. Misión Técnica Alemana 1980
Ex-Programa Nacional de Pequeñas y Medianas Irrigaciones
 PRECIPITACION EFECTIVA
USBR MODIFICADO

Precipitación PORCION DE LA PRECIPITACION

Total Mensual Aprovechable por las Plantas Déficit ó Escorrentía
Límite Superior CI C II C III PE I PE II PE III
mm. mm. mm. mm. mm. mm. mm.

25.4 25.4 22.9 20.3 0.0 2.5 5.1

50.8 49.5 44.5 39.4 1.3 6.3 11.4
76.2 72.4 63.5 54.6 3.8 12.7 21.6
101.6 92.7 76.2 59.7 8.9 25.4 41.9
127.0 107.9 83.8 59.7 19.1 43.2 67.3
152.4 108.1 86.4 59.7 34.3 56.0 92.7
177.8 120.6 86.4 59.7 57.2 91.4 118.1
Fuente: Generación de Caudales para la Sierra del Perú. Misión Técnica Alemana 1980

Ex-Programa Nacional de Pequeñas y Medianas Irrigaciones

grado:
PE = ao + a1P1 + a2P2 + a3P3 + a4P4 +
a5P5

Los límites de la PE son:

CPEI : PE=P-120.6 para P=177.8
mm/mes
CPEII :PE=P-86.4 para P=152.4
mm/mes
CPEIII :PE=P-59.7 para P=127.6
mm/mes
COEFICIENTES CURVA I CURVA II CURVA III

Coeficientes
a0
del
-0.018
Polinomio
-0.021
para el Cálculo
-0.028
PE a1 -0.0185 0.1358 0.2756
a2 0.001105 -0.002296 -0.004103
a3 -1204 E-8 4349 E-8 5543 E-8
a4 -1440 E-9 -890 E-9 1240 E-9
a5 -2850 E-12 -879 E-13 -1420 E-11

E : Simbolo que indica que la cifra siguiente vale como Exponente a la base 10
FUENTE : Generación de Caudales para la Sierra del Perú. Misión Técnica Alemana 1980
Ex-Programa Nacional de Pequeñas y Medianas Irrigaciones
Como se ha mencionado anteriormente la PE se calcula
con el polinomio de quinto grado, teniendo en cuenta que
la suma de los valores de Pei (de los 12 meses del año) se
igual al producto: Precipitación total anual (P) por el
Coeficiente de escurrimiento (C ).
𝒏=𝟏𝟐

∑ 𝑷𝑬𝒊=𝑪 . 𝑷
𝒊=𝟎

CPEI para C = 0.15

CPEII para C = 0.30
CPCIII para C = 0.50
 Se debe ajustar las curvas de la PEi, es decir si la cuenca
tiene un coeficiente C = 0.40, entonces se determinará la
PEi de las CPEII y CPEIII, mediante la interpolación.
 1) k1.PEII + k2.PEIII = PE
 2) siendo k1 + k2 =1
Donde:
k1 y k2 : Coeficientes de interpolación
PEII y PEIII: PE según las curvas II y III.
PE : Precipitación efectiva.
4.0 RETENCION DE LA CUENCA (R): Depende directamente de las características
geomorfológicas de la cuenca, principalmente de los almacenamientos hídricos naturales.
Se distinguen tres tipod de mayor importancia:
o Acuíferos
o Lagunas y Pantanos
o Nevados

LAMINA DE AGUA ACUMULADA

EN LOS TRES TIPOS DE ALMACENAMIENTOS
Tipo Pendiente REi Ai
Almacenamientos (Areas Parciales)
% mm/año km2

2 300 Ai
Acuíferos 8 250 Ai
Potenciales 15 200 Ai

Lagunas 500 Ai

Nevados 500 Ai

Retención = ΣREi*Ai
ACuenca
REi = Retención de cada tipo
Ai= Area parcial de cada tipo
Fuente: Generación de Caudales para la Sierra del Perú. Misión Técnica Alemana 1980
Ex-Programa Nacional de Pequeñas y Medianas Irrigaciones
La retención total se transforma en una lámina de
agua sobre la cuenca para el año promedio en
mm/año, que por definición es igual al Gasto (G)
durante la estación seca e igual al Abasto (A) durante
la estación lluviosa, transformándose el Balance
Hídrico en :

 Balance Estación Seca :

Q = PE + G
 Balance Hídrico Estación Lluviosa :
Q = PE – A
La secuencia del Proceso es la siguiente:
 4a. GASTO DE LA RETENCION DE LA CUENCA (G):

Se presenta en la estación seca afectada por el Coeficiente de Agotamiento,

para su estimación se utiliza aforos o en base a las relaciones siguientes:

- Cuencas con agotamiento muy rápido. Debido a temperaturas elevadas (>10°C)

y retención que va de reducida
(50 mm/año) a mediana (80 mm/año):
a = - 0.00252.LN(AR) + 0.034
- Cuencas con agotamiento rápido. Retención entre 50 y 80 mm/año y vegetación
poco desarrollada (puna):
a = - 0.00252.LN(AR) + 0.030
- Cuencas con agotamiento mediano. Retención mediana (80 mm/año) y vegetación
mezclada (pastos, bosques
y terrenos cultivados):
a = - 0.00252.LN(AR) + 0.026
- Cuencas con agotamiento reducido. Debido a la alta retención (> 100 mm/año) y
vegetación mezclada:
a = - 0.00252.LN(AR) + 0.023
donde:
a = coeficiente de agotamiento por día
AR = área de la cuenca (km2)
LN = Logaritmo Natural
Durante la estación seca el Gasto (+) de la Retención alimenta a los ríos y
mantiene el caudal básico. Esta reserva se agota al final de la estación seca y
durante este tiempo se puede representar matematicamente la curva de
agotamiento a través de la formula de Maillet:
Qt = Qo. e^(-a.t) O también Qt/Qo = e^(-30.a) = bo^(i)

Donde:
Qt : Descarga en el tiempo t
Qo : Descarga Inicial
a : Coeficiente de agotamiento
t : Tiempo (30 días)
bo : Coeficiente del gasto del mes de inicio
i : Valor en el primer mes es =1 y se incrementa en una unidad
en forma sucesiva hasta el final del último mes de estiaje.
Por lo que la contribución mensual del gasto de la retención
en la época seca es:
Gi = R (bo)^i / Suma(bo)^i
Donde:
Gi : Gasto mensual de la retención en mm/mes
R : Retención de la Cuenca en mm/mes
 4b. ABASTECIMIENTO DE LA RETENCION DE LA CUENCA (A):

Corresponde a la estación lluviosa y su descrpción no se puede precisar

mediante una fórmula sencilla por la complejidad de las características
geomorfológicas y climáticas de la cuenca, teniendo en consideración
que la suma de los valores del los coeficientes ai=100%, en función de la
Retención de la cuenca.
La lámina de agua Ai, que se adiciona como reserva de la cuenca se
muestra como un déficit mensual de la PE y se calcula mediante la
ecuación:
Ai = ai . R /100
Donde:
Ai : Abastecimiento mensual, déficit de la PE, en mm/mes
ai : Coeficiente de abastecimiento en %
R : Retención de la cuenca en mm/mes.
 APLICACION MICROCUENCA QDA. UCANAN

Esta microcuenca drena sus aguas por la margen derecha del río
Chillón; esta quebrada se ubica en la parte alta, entre las cotas
3340 a 4670 msnm; en un área de drenaje de 5.379 km2. Tiene una
longitud de río principal de 6.461 km y una pendiente media de
20.58 % y su altitud media es de 3976 msnm.
 ELEMENTOS DEL MODELO

1. PRECIPITACION : CUADRO_01

2. COEFICIENTE DE ESCURRIMIENTO: C=0.48

3. PRECIPITACION EFECTIVA : CUADRO_02

Se calcula entre la curva PEII y PEIII, obteniendo K1=-1.08018 y K2= 2.08018
4. RETENCION : 32.56 mm/año.
5. COEFICIENTE AGOTAMIENTO : Formula Agotamiento Rápido = 0.0298
Abr. b1 = 0.4090
May. b2 = 0.1673
Jun. b3 = 0.0684
Jul. b4 = 0.0280
Ago. b5 = 0.0114
6. COEFICIENTE ABASTECIMIENTO:

ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC
16.6 29.4 35.3 5.1 1.0 0.9 2.7 3.1 5.9 = 100 %
 Retención de la Microcuenca Ucanan
Tipo Pendiente REi Ai (A. Parciales) Km2
Ucanan
% mm/año 5.379

Acuíferos 2 300 0.03318

Potenciales 8 250 0.11618
15 200 0.68081

Lagunas 500

Nevados 500
Retención = ΣREi*Ai Retención en mm/año
Atotal Cuenca 32.56

REi = Retención de cada tipo

Ai= Area parcial de cada tipo
Area Total en km2

BALANCE HIDROLOGICO PARA EL AÑO PROMEDIO

Precipitación en la Microcuenca Contribución de la Retención Caudales
Mes Total Efectiva Gasto Abastecimiento Generados
P PE II PEIII PE BI GI ABA AI
mm/mes mm/mes mm/mes mm/mes mm/mes mm/mes mm/mes l/s

ENE. 74.02 12.24 20.93 30.32 0.000 0.00 0.166 5.41 24.9 50.01
FEB. 97.17 22.72 37.57 53.6 0.000 0.00 0.294 9.57 44.03 97.9
MAR. 105.64 27.75 45.32 64.31 0.000 0.00 0.353 11.49 52.82 106.08
ABR. 29.42 3.05 6.03 9.25 0.409 19.47 0.051 1.65 27.06 56.16
MAY. 4.31 0.55 1.12 1.73 0.1673 7.96 0.01 0.31 9.38 18.85
JUN. 0.03 0 0.01 0.01 0.0684 3.26 0.00 0.00 3.27 6.78
JUL. 0.07 0.01 0.02 0.03 0.028 1.33 0.00 0.01 1.36 2.73
AGO. 0.12 0.02 0.03 0.05 0.0114 0.54 0.00 0.01 0.59 1.18
SET. 4.06 0.52 1.06 1.64 0.000 0.00 0.009 0.29 1.35 2.79
OCT. 13.84 1.55 3.18 4.94 0.000 0.00 0.027 0.88 4.05 8.14
NOV. 16.06 1.76 3.6 5.59 0.000 0.00 0.031 1.00 4.59 9.53
DIC. 35.02 3.67 7.11 10.82 0.000 0.00 0.059 1.93 8.89 17.85

AÑO 379.76 73.84 125.98 182.29 0.6841 32.56 1.000 32.55 182.29 31.5
 PARAMETROS DEL MODELO DE GENERACION

MES Q(T) Q(T-1) PE(T)

mm mm mm

ENE. 24.90 8.89 30.32

FEB. 44.03 24.90 53.60
MAR. 52.82 44.03 64.31
ABR. 27.06 52.82 9.25
MAY. 9.38 27.06 1.73
JUN. 3.27 9.38 0.01
JUL. 1.36 3.27 0.03
AGO. 0.59 1.36 0.05
SET. 1.35 0.59 1.64
OCT. 4.05 1.35 4.94
NOV. 4.59 4.05 5.59
DIC. 8.89 4.59 10.82

Parámetros para la Generación:

B1 = 0.437
B2 = 0.354
B3 = 0.617
Desviación Estándar Total = 1.498
Coeficente de Correlación Múltiple = 0.996

MODELO DE GENERACION DE CAUDALES MEDIOS MENSUALES

Q(I,J)=B1+B2 * Q(I,J-1) + B3 * PE(I,J) + S(1-R²)½*Z(J)

Donde:
Q(I,J) = caudal generado del año I, mes J.
Q(I,J-1) = caudal del año I, mes (J-1).
PE(I,J) = precipitación efectiva del año I, mes J.
S = error estándar.
R = coeficiente de correlación múltiple.
Z(J) = número aleatorio normal (0,1) del mes J.
B1, B2 y B3 = coeficientes de regresión.

Q(I,J)= 0.437+0.354 * Q(I,J-1) + 0.617 * PE(I,J) + 1.498*(1-0.996²)½*Z(J)

Con r = 0.996
 Calibración del Modelo
Para la calibración del modelo se tiene que tener datos de
aforos históricos del rio/quebrada. Generalmente se tiene
carencia de esta información, por lo que se contrastó las
descargas medias generadas con las del río Chillón
registradas en la estación Puente Magdalena con un área
de drenaje de 1010 km2. En vista que la microcuenca
tiene un área muy pequeña, se usó las descargas
específicas (m3/s/km2)
El periodo de calibración de los caudales históricos del río
Chillón en la estación Puente Magdalena es de 30 años
1948/1978
El modelo se ejecutó con los parámetros encontrados
inicialmente y comparados con los caudales históricos del
río Chillón, a fin de encontrar series sintéticas que
tengan la tendencia de los caudales históricos. Por tal
motivo se realizó una serie de ejecuciones del modelo a
fin de calibrar sus parámetros. En este caso se cambiaron
los coeficientes de abastecimiento en cada corrida.

DESCRIPCION ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC
Q_Histórico Río Chillón Area=1010 km2 (1948/1978)
QPromedio Histórico m3/s 11.160 19.840 23.940 12.120 3.590 2.130 1.490 1.170 0.790 1.120 1.790 5.100
Qespecífico m3/s/km2 0.011 0.020 0.024 0.012 0.004 0.002 0.001 0.001 0.001 0.001 0.002 0.005
Q_Generado Qda. Ucanan Area=5.379 km2 (1964/2013)
QPromedio Generado m3/s 0.060 0.106 0.125 0.059 0.023 0.009 0.004 0.002 0.004 0.008 0.011 0.021
Qespecífico m3/s/km2 0.011 0.020 0.023 0.011 0.004 0.002 0.001 0.000 0.001 0.001 0.002 0.004

0.03
CAUDALES m3/s/km2

0.02

0.02

0.01

0.01

0.00
ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC
Serie_Generada 0.011 0.020 0.023 0.011 0.004 0.002 0.001 0.000 0.001 0.001 0.002 0.004
Serie Histórica 0.011 0.020 0.024 0.012 0.004 0.002 0.001 0.001 0.001 0.001 0.002 0.005

TIEMPO EN MESES
 Evaluacion Estadística de los Resultados

Para evaluar los resultados de la calibración

y validación Se han utilizado parámetros
estadísticos, recomendados por Moriasi et
al, 2007; Debels, 2010, los cuales se
describen a continuación:
• Indice de eficiencia de Nash-Sutcliffe-NSE
(1970) es una estadístico normalizado que
determina la magnitud relativa de la
varianza residual en comparación con la
variación de los datos de medición. Este
índice fue empleado como criterio para
evaluar el poder reproductivo del modelo
hidrológico (Debels, 2010), adoptando la
interpretación recomendada por la
literatura (Moriasi et. al, 2007), según
Tabla de Interpretación.
El Indice de eficiencia de Nash se
representa en la siguiente fórmula:

N
2
(Q
Observado
i
Simulado
Qi )
NSE = 1
i=1
N
2
i=1
(Q Observado
i Q
Promedio
)
Otros índices de error son comúnmente
utilizados en la evaluación del modelo
(Moriasi et al, 2007; Debels, 2010); más los
utilizados se describen a continuación:
• Razón RMSE/ Desviación Estándar de las
observaciones (RSR) es una versión
estandarizada del RMSE (Legates &
McCabe, 1999); se calcula como el
cociente entre el RMSE y la desviación
estándar de los datos observados, como se
muestra en la siguiente ecuación y se
interpreta basados en los criterios
establecidos por Moriasi et. al. (2007), que
se resumen en la según Tabla
Interpretación.

N
2
RMSE (Q
Observado
i
Simulado
Qi )
RSR = = i=1

DESVESTObs
N

(Q Observado Promedio
)
2
i Q
i=1
• Sesgo porcentual (Percent Bias-PBIAS) es la
tendencia promedio de los datos simulados
de ser mayores o menores que los datos
observados (Gupta et al., 1999). Se adoptan
los criterios recomendados por Moriasi et.
al. (2007) para su interpretación (según
Tabla de Interpretación). La expresión para
el cálculo es la siguiente ecuación:

(Q Observado
i
Simulado
Qi )
PBIAS = i=1
N
X 100

i=1
(Q Observado
i )

TABLA DE INTERPRETACION DE LOS PARAMETROS ESTADISTICOS PARA

CALIBRACION Y VALIDADCION
 RESULTADOS DE LOS PARAMETROS ESTADISTICOS
DE LA CALIBRACION

MICROCUENCA PARAMETROS ESTADISTICOS

PBIAS NSE RSR
PBIAS<+-10 0.76<NSE<=1.0 0.00<RSR<=0.50

INTERPRETACION Muy Bueno Muy Bueno Muy Bueno

QDA.UCANAN 3.663 0.993 0.083
DETERMINACION DE LA PENDIENTE EN UNA CUENCA
MODELO BALANCE HIDROLOGICO – LUTZ
ARCGIS 10.0

1 ARCTOOLBOX
CURVAS DE NIVEL

3D ANALYST
CUENCA CORTADA
TIN MANAGEMENT
2 CREACION TIN
ABRIR TABLA ATRIBUTO
CONVERSION 1. AÑADIR CAMPO NUMERICO
OF TIN PARA EL AREA Y SE CALCULA
3 CREACION RASTER 2. AÑADIR CAMPO TXT PENDIENTE
Cellsize 5m
RASTER SUFACE
SLOPE Seleccionar por Atributos
4 Calculadora\clip\GridCode:
%
"GridCode" Escribir: = 1
RASTER RECLASS APLICAR
CLIP DERECHO COLUMNA PEND.
RECLASSIFY
FIELD CALCULATOR:
5 RECLASIFICACION ESCRIBIR: " 2 %"
2%, 8%, 15%, >15 Repetir para los otras pendientes
CONVERSION TOOLS
FROM RASTER CLIP DERECHO COLUMNA PEND.
SUMARIZE
6 RASTER TO POLIGONO ABRIR \+ AREA
CLIP: SUM
ANALISYS TOOLS
EXTRACT \ CLIP SE OBTIENE
7 1. ARCHIVO SHAPE POLIGONO TABLA CON LOS
2. ARCHIVO SHAPE LIMITE CUENCA RESULTADOS

http://glcfapp.glcf.umd.edu:8080/esdi/
 CUADRO 01_PRECIPITACION TOTAL MENSUA
Microcuenca Qda. Ucanan

Año Ene Feb Mar Abr May Jun Jul Ago

1964 65.4 117.5 65.5 48.7 19.3 0.0 0.0 0.0

1965 99.3 76.4 145.7 31.5 9.4 0.0 0.0 0.1
1966 19.0 89.4 151.5 30.9 0.1 0.0 0.0 0.0
1967 173.6 177.7 79.6 46.5 22.9 0.0 0.0 0.0
1968 34.7 21.3 16.9 9.9 0.3 0.0 0.0 0.0
1969 11.9 72.1 112.0 32.0 0.2 0.0 0.0 0.0
1970 210.5 47.1 85.3 52.4 0.4 0.0 0.0 0.0
1971 79.6 138.3 160.0 45.6 4.7 0.0 0.0 0.0
1972 80.4 97.7 160.9 46.0 21.4 0.0 0.0 0.0
1973 111.4 77.0 11.6 56.4 0.5 0.0 0.0 0.1
1974 95.7 93.4 116.5 41.3 5.8 0.2 0.1 0.0
1975 54.9 113.2 136.7 23.0 0.0 0.3 0.2 0.3
1976 104.1 156.7 7.9 17.0 0.4 0.0 0.0 0.0
1977 41.8 93.8 154.2 31.4 0.2 0.0 0.0 0.0
1978 74.9 85.2 160.8 32.4 0.1 0.0 0.0 0.0
1979 28.8 77.8 157.5 36.3 0.0 0.0 0.0 0.0
1980 34.6 44.6 65.8 4.5 1.2 0.0 0.0 0.0
1981 47.2 150.2 148.3 10.5 0.0 0.0 0.0 0.0
1982 99.6 78.3 18.7 18.8 0.0 0.0 0.0 0.1
1983 49.0 101.7 160.9 30.7 0.1 0.0 0.1 0.0
1984 114.6 192.7 192.5 45.0 18.6 0.0 0.0 0.0
1985 11.7 102.9 97.9 8.4 0.0 0.0 0.0 0.0
1986 57.1 62.6 131.8 79.8 13.4 0.0 0.0 0.0
1987 128.7 111.2 0.0 0.0 0.0 0.0 0.0 0.0
1988 56.0 52.8 37.9 31.9 3.5 0.0 0.0 0.0
1989 201.3 178.0 209.5 41.7 13.2 0.0 0.0 0.0
GRACIAS POR LA
ATENCION