Modelo Determinístico Estocástico de Lutz Scholz (Scholz, 1980)

GENERALIDADES
Este modelo hidrológico, es combinado por que cuenta con una estructura
determinística para el cálculo de los caudales mensuales para el año promedio
(Balance Hídrico – Modelo determinístico); y una estructura estocástica para la
generación de series extendidas de caudal (Proceso markoviano – Modelo
Estocástico). Fue desarrollado por el experto Lutz Scholz para cuencas de la sierra
peruana, entre los años 1979-1980, en el marco de Cooperación Técnica de la
República de Alemania a través del Plan Meris II.

Determinado el hecho de la ausencia de registros de caudal en la sierra peruana, el

modelo se desarrolló tomando en consideración parámetros físicos y meteorológicos
de las cuencas, que puedan ser obtenidos a través de mediciones cartográficas y de
campo. Los parámetros más importantes del modelo son los coeficientes para la
determinación de la Precipitación Efectiva, déficit de escurrimiento, retención y
agotamiento de las cuencas. Los procedimientos que se han seguido en la
implementación del modelo son:

1. Cálculo de los parámetros necesarios para la descripción de los fenómenos de

escorrentía promedio.

2. Establecimiento de un conjunto de modelos parciales de los parámetros para el

cálculo de caudales en cuencas sin información hidrométrica. En base a lo
anterior se realiza el cálculo de los caudales necesarios.
3. Calibración del modelo y generación de caudales extendidos por un proceso
Markoviano combinado de precipitación efectiva del mes con el caudal del mes
anterior.

Este modelo fue implementado con fines de pronosticar caudales a escala mensual,
teniendo una utilización inicial en estudios de proyectos de riego y posteriormente
extendiéndose el uso del mismo a estudios hidrológicos con prácticamente cualquier
finalidad (abastecimiento de agua, hidroelectricidad etc). Los resultados de la
aplicación del modelo a las cuencas de la sierra peruana, han producido una
correspondencia satisfactoria respecto a los valores medidos.
ECUACIÓN DEL BALANCE HÍDRICO
La ecuación fundamental que describe el balance hídrico mensual en mm/mes es la
siguiente:
𝐶𝑀𝑖 = 𝑃𝑖 − 𝐷𝑖 + 𝐺𝑖 − 𝐴𝑖 (7)
Dónde:
𝐶𝑀𝑖 = Caudal mensual (mm/mes)
𝑃𝑖 = Precipitación mensual sobre la cuenca (mm/mes)
𝐷𝑖 = Déficit de escurrimiento (mm/mes)
𝐺𝑖 = Gasto de la retención de la cuenca (mm/mes)
𝐴𝑖 = Abastecimiento de la retención (mm/mes)

Asumiendo:

1. Que para períodos largos (en este caso 1 año) el Gasto y Abastecimiento de la
retención tienen el mismo valor es decir Gi = Ai, y

2. Que para el año promedio una parte de la precipitación retorna a la atmósfera

por evaporación.

Reemplazando (P-D) por (CP), y tomando en cuenta la transformación de unidades

(mm/mes a m3/seg) la ecuación (7) se convierte en:

𝑄 = 𝑐 ′ 𝐶𝑃(𝐴𝑅) (8)

Que es la expresión básica del método racional.

Dónde:
Q = Caudal (m3/s)
c’ = coeficiente de conversión del tiempo (mes/seg)
C = coeficiente de escurrimiento
P = Precipitación total mensual (mm/mes)
AR = Área de la cuenca (m2)
COEFICIENTE DE ESCURRIMIENTO
Se ha considerado el uso de la fórmula propuesta por L. Turc:

𝑃−𝐷
𝐶= (9)
𝑃

Los coeficientes parciales de las curvas se encuentran con las expresiones.

𝐶𝑃−𝑃𝐸𝐼𝐼
𝐶𝐼 = (10.1)
𝑃𝐸𝐼−𝑃𝐸𝐼𝐼

𝐶𝑃−𝑃𝐸𝐼𝐼𝐼
𝐶𝐼𝐼 = (10.2)
𝑃𝐸𝐼𝐼−𝑃𝐸𝐼𝐼𝐼

𝐶𝑃−𝑃𝐸𝐼𝐼
𝐶𝐼𝐼𝐼 = (10.3)
𝑃𝐸𝐼𝐼𝐼−𝑃𝐸𝐼𝐼
Dónde:
C = Coeficiente de escurrimiento
CI, CII y CIII= Coeficientes para cada grupo de curvas
P = Precipitación Total anual (mm/año)
PEI, PEII, PEIII= Precipitación efectiva para cada grupo de curvas (mm)
D = Déficit de escurrimiento (mm/año)

Para la determinación de D se utiliza la expresión:

1
𝐿 = 300 + 25(𝑡 ) + 0.05(𝑇)3 = 646.70, 𝐷=𝑃 1 (11)
𝑃2 2
(0.9+ 2 )
𝐿

Siendo:
L = Coeficiente de Temperatura
T = Temperatura media anual (°C)

Dado que no se ha podido obtener una ecuación general del coeficiente de escorrentía
para la toda la sierra, se ha desarrollado la fórmula siguiente, que es válida para la
región sur:

C  3.16E12( P 0.571)( EP 3.686 );    r  0.96 (12)

D  1380  0.872( P)  1.032( EP );      r  0.96 (13)
Dónde:
C = Coeficiente de escurrimiento
D = Déficit de escurrimiento (mm/año)
P = Precipitación total anual (mm/año)
EP = Evapotranspiración anual según Hargreaves (mm/año)
r = Coeficiente de correlación

La evapotranspiración potencial, se ha determinado por la fórmula de Hargreaves:

𝐸𝑃 = 0.0075(𝑅𝑆𝑀)(𝑇𝐹)(𝐹𝐴) (14)

𝑛
𝑅𝑆𝑀 = 0.075(𝑅𝐴)√ (15)
𝑁

𝐹𝐴 = 1 + 0.06(𝐴𝐿) (16)

Dónde:
RSM = Radiación solar media
TF = Componente de temperatura
FA = Coeficiente de corrección por elevación
TF = Temperatura media anual (°F)
RA = Radiación extraterrestre (mm de agua / año)
(n/N) = Relación entre insolación actual y posible (%)
50 % (estimación en base a los registros)
AL = Elevación media de la cuenca (Km)

Para determinar la temperatura anual se toma en cuenta el valor de los registros de

las estaciones y el gradiente de temperatura de -5.3 °C 1/ 1000 m, determinado para
la sierra.

PRECIPITACIÓN EFECTIVA

Para el cálculo de la Precipitación Efectiva, se supone que los caudales promedio

observados en la cuenca pertenecen a un estado de equilibrio entre gasto y
abastecimiento de la retención. La precipitación efectiva se calculó para el coeficiente
de escurrimiento promedio, de tal forma que la relación entre precipitación efectiva y
precipitación total resulta igual al coeficiente de escorrentía.
Para fines hidrológicos se toma como precipitación efectiva la parte de la precipitación
total mensual, que corresponde al déficit según el método del USBR (precipitación
efectiva hidrológica es el antítesis de la precipitación efectiva para los cultivos).
A fin de facilitar el cálculo de la precipitación efectiva se ha determinado el polinomio
de quinto grado:

𝑃𝐸 = 𝑎0 + 𝑎1 𝑃 + 𝑎2 𝑃 2 + 𝑎3 𝑃 3 + 𝑎4 𝑃 4 + 𝑎5 𝑃 5 (17)

Dónde:
PE = Precipitación efectiva (mm/mes)
P = Precipitación total mensual (mm/mes)
𝑎 i = Coeficiente del polinomio

A continuación se muestra los valores límite de la precipitación efectiva y la tabla 01

muestra los tres juegos de coeficientes, ai, que permiten alcanzar por interpolación
valores de C, comprendidos entre 0.15 y 0.45.

Límite superior para Precipitación Efectiva

Curva I: PE=P-120.6 para P> 177.8 mm/mes

Curva II: PE=P-86.4 para P>152.4 mm/mes

Curva III: PE=P-59.7 para P>127.0 mm/mes

Tabla 01: Coeficientes para el cálculo de precipitación efectiva

Coef. I II III
a0 -0.018 -0.021 -0.028
a1 -0.019 0.1358 0.2756
a2 0.0011 -0.002 -0.004
a3 -1E-05 4E-05 6E-05
a4 1E-07 -9E-08 1E-07
a5 -3E-10 -9E-11 -1E-09

De esta forma es posible llegar a la relación entre la precipitación efectiva y

precipitación total
 𝑄 𝑃𝐸𝑖
𝐶= = ∑12
𝑖=1 (18)
𝑃 𝑃

Dónde:

C = Coeficiente de escurrimiento
Q = Caudal anual
P = Precipitación Total anual

∑12
𝑖=1 𝑃𝐸𝑖 = suma de la precipitación efectiva

RETENCIÓN DE LA CUENCA

Bajo la suposición de que exista un equilibrio entre el gasto y el abastecimiento de la

reserva de la cuenca y además que el caudal total sea igual a la precipitación efectiva
anual, la contribución de la reserva hídrica al caudal se puede calcular según las
fórmulas:

𝑅𝑖 = 𝐶𝑀𝑖 − 𝑃𝑖 (19.1)

𝐶𝑀𝑖 = 𝑃𝐸𝑖 + 𝐺𝑖 − 𝐴𝑖 (19.2)

Dónde:
𝐶𝑀𝑖 = Caudal mensual (mm/mes)
𝑃𝐸𝑖 = Precipitación Efectiva Mensual (mm/mes)
𝑅𝑖 = Retención de la cuenca (mm/mes)
𝐺𝑖 = Gasto de la retención (mm/mes)
𝐴𝑖 = Abastecimiento de la retención (mm/mes)
𝑅𝑗 =𝐺𝑖 para valores mayores que cero (mm/mes)
𝑅𝑗 = 𝐴𝑖 para valores menores que cero (mm/mes)

Sumando los valores de G o A respectivamente, se halla la retención total de la cuenca

para el año promedio, que para el caso de las cuencas de la sierra varía de 43 a 188
(mm/año).

Relación entre descargas y retención

Durante la estación seca, el gasto de la retención alimenta los ríos, constituyendo el

caudal o descarga básica. La reserva o retención de la cuenca se agota al final de la
estación seca; durante esta estación la descarga se puede calcular en base a la
ecuación:
 𝑄𝑡 = 𝑄0 𝑒 −𝑎(𝑡) (20)

Dónde:
𝑄𝑡 = Descarga en el tiempo t
𝑄0 = Descarga inicial
𝑎 = Coeficiente de agotamiento
t = Tiempo
Al principio de la estación lluviosa, el proceso de agotamiento de la reserva termina,
comenzando a su vez el abastecimiento de los almacenes hídricos. Este proceso está
descrito por un déficit entre la precipitación efectiva y el caudal real. En base a los
hidrogramas se ha determinado que el abastecimiento es más fuerte al principio de la
estación lluviosa continuando de forma progresiva pero menos pronunciada, hasta el
final de dicha estación

Coeficiente de agotamiento

Mediante la ecuación (100.1) se puede calcular el coeficiente de agotamiento “a”, en

base a datos hidrométricos. Este coeficiente no es constante durante toda la estación
seca, ya que va disminuyendo gradualmente.
Con fines prácticos se puede despreciar la variación del coeficiente “a” durante la
estación seca empleando un valor promedio.
El coeficiente de agotamiento de la cuenca tiene una dependencia logarítmica del área
de la cuenca.

𝑎 = 𝑓(𝐿𝑛𝐴𝑅) (21.1)

𝑎 = 3.1249𝐸67(𝐴𝑅)−0.1144 (𝐸𝑃)−19.336 (𝑇)−3.369 (𝑅)−1.429 (21.2)

𝑟 = 0.86

El análisis de las observaciones disponibles muestra, además cierta influencia del

clima, la geología y la cobertura vegetal. Se ha desarrollado una ecuación empírica
para la sierra peruana:
En principio, es posible determinar el coeficiente de agotamiento real mediante aforos
sucesivos en el río durante la estación seca; sin embargo cuando no sea posible ello,
se puede recurrir a las ecuaciones desarrolladas para la determinación del coeficiente
“a” para cuatro clases de cuencas:

- 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):
 𝑎 = −0.00252(𝐿𝑛𝐴𝑅) + 0.034 (21.3)

- Cuencas con agotamiento rápido. Retención entre 50 y 80 mm/año y vegetación

poco desarrollada (puna):
𝑎 = −0.00252(𝐿𝑛𝐴𝑅) + 0.030 (21.4)

- Cuencas con agotamiento mediano. Retención mediana (80 mm/año) y vegetación

mezclada (pastos, bosques y terrenos cultivados):
𝑎 = −0.00252(𝐿𝑛𝐴𝑅) + 0.026 (21.5)
- Cuencas con agotamiento reducido. Debido a la alta retención (> 100 mm/año) y
vegetación mezclada:
𝑎 = −0.00252(𝐿𝑛𝐴𝑅) + 0.023 (21.6)

Dónde:
𝑎 = coeficiente de agotamiento por día
AR = área de la cuenca (km2)
EP = evapotranspiración potencial anual (mm/año)
T = duración de la temporada seca (días)
R = retención total de la cuenca (mm/año)

ALMACENAMIENTO HÍDRICO

Tres tipos de almacenes hídricos naturales que inciden en la retención de la cuenca

son considerados:
 Acuíferos
 Lagunas y pantanos
 Nevados

La determinación de la lámina “L” que almacena cada tipo de estos almacenes está
dado por:

- Acuíferos:
𝐿𝐴 = −750(𝐼 ) + 315 (𝑚𝑚/𝑎ñ𝑜) (22.1)

Siendo:
LA = lámina específica de acuíferos
I = pendiente de desagüe: I <= 15 %

𝐿𝐿 = 500 (𝑚𝑚/𝑎ñ𝑜) (22.2)

- Lagunas y Pantanos
Siendo:
LL = Lámina específica de lagunas y pantanos

- Nevados
𝐿𝑁 = 500(𝑚𝑚/𝑎ñ𝑜) (22.3)

Siendo:
LN = lámina específica de nevados

Las respectivas extensiones o áreas son determinadas de los mapas o aerofotografías.

Los almacenamientos de corto plazo no son considerados para este caso, estando los
mismos incluidos en las ecuaciones de la precipitación efectiva.

ABASTECIMIENTO DE LA RETENCIÓN

El abastecimiento durante la estación lluviosa es uniforme para cuencas ubicadas en

la misma región climática. En la región del Cusco el abastecimiento comienza en el
mes de noviembre con 5%, alcanzando hasta enero el valor del 80 % del volumen final.
Las precipitaciones altas del mes de febrero completan el 20 % restante, y las
precipitaciones efectivas del mes de marzo escurren directamente sin contribuir a la
retención. Los coeficientes mensuales expresados en porcentaje del almacenamiento
total anual se muestran en la tabla 02.

Tabla 02: Almacenamiento hídrico durante la época de lluvias (%)

Región Oct Nov Dic Ene Feb Mar Total

Cusco 0 5 35 40 20 0 100
Huancavelica 10 0 35 30 20 5 100
Junín 10 0 25 30 30 5 100
Cajamarca 25 -5 0 20 25 35 100

La lámina de agua Ai que entra en la reserva de la cuenca se muestra en forma de

déficit mensual de la Precipitación Efectiva PEi. Se calcula mediante la ecuación:
𝑅
𝐴𝑖 = 𝑎𝑖 ( ) (23)
100

Siendo:
𝐴𝑖 = abastecimiento mensual déficit de la precipitación efectiva (mm/mes)
𝑎𝑖 = coeficiente de abastecimiento (%)
R = retención de la cuenca (mm/año)

DETERMINACIÓN DEL CAUDAL MENSUAL PARA EL AÑO PROMEDIO

Está basado en la ecuación fundamental que describe el balance hídrico mensual a

partir de los componentes descritos anteriormente:

𝐶𝑀𝑖 = 𝑃𝐸𝑖 + 𝐺𝑖 − 𝐴𝑖 (24)

Dónde:
𝐶𝑀𝑖 = Caudal del mes i (mm/mes)
𝑃𝐸𝑖 = Precipitación efectiva del mes i (mm/mes)
𝐺𝑖 = Gasto de la retención del mes i (mm/mes)
𝐴𝑖 = abastecimiento del mes i (mm/mes)

1.2.2.11.- GENERACIÓN DE CAUDALES PARA PERIODOS EXTENDIDOS

A fin de generar una serie sintética de caudales para períodos extendidos, se ha

implementado un modelo estocástico que consiste en una combinación de un Proceso
Markoviano de primer orden, según la ecuación (25) con una variable de impulso, que
en este caso es la precipitación efectiva en la ecuación (26):

𝑄𝑡 = 𝑓 (𝑄𝑡=1 ) (25)

𝑄 = 𝑔(𝑃𝐸𝑡 ) (26)

Con la finalidad de aumentar el rango de valores generados y obtener una óptima

aproximación a la realidad, se utiliza además una variable aleatoria.

𝑍 = 𝑧(𝑆)√(1 − 𝑟 2 ) (27)

𝑄𝑡 = 𝐵1 + 𝐵2 (𝑄𝑡−1 ) + 𝐵3 (𝑃𝐸𝑡 ) + 𝑧(𝑆)√1 − 𝑟 2 (28)

La ecuación integral para la generación de caudales mensuales es:

Dónde:

Qt = Caudal del mes t

Q t-1 = Caudal del mes anterior
 PE t = Precipitación efectiva del mes
B1 = Factor constante o caudal básico.

Se calcula los parámetros B1, B2, B3, r y S sobre la base de los resultados del modelo
para el año promedio por un cálculo de regresión con Qt como valor dependiente y Qt-
1 y PEt, como valores independientes. Para el cálculo se recomienda el uso de software

comercial (hojas electrónicas) o de uso específico (programas elaborados tales como

el SIH).

El proceso de generación requiere de un valor inicial, el cual puede ser obtenido en

una de las siguientes formas:
 Empezar el cálculo en el mes para el cual se dispone de un aforo
 Tomar como valor inicial el caudal promedio de cualquier mes,
 Empezar con un caudal cero, calcular un año y tomar el último valor como
valor Qo sin considerar estos valores en el cálculo de los parámetros
estadísticos del período generado.

TEST ESTADÍSTICOS

Para determinar la calidad de la coincidencia de los caudales generados con los

observados, se desarrolla la comparación de los promedios y desviaciones tipo de los
valores históricos y los generados.

Para probar si los promedios salen de la misma población, se utiliza el test de Student
(Prueba "t"). Esta prueba debe ser desarrollada para cada mes. Se compara el valor
de t con el valor límite tp,n que indica el límite superior que, con una probabilidad de
error del P%, permite decir que ambos promedios pertenecen a la misma población.

La comparación estadística de promedios se realiza mediante el test de Fischer

(Prueba "F") que se compara con el valor límite Fp/2 (%) , (n1,n2)

RESTRICCIONES DEL MODELO

El modelo presenta ciertas restricciones de uso o aplicación tales como:

 El uso de los modelos parciales, únicamente dentro del rango de calibración

establecido.
 Su uso es únicamente para el cálculo de caudales mensuales promedio.
  Los registros generados en el período de secas presentan una mayor
confiabilidad que los valores generados para la época lluviosa.
 La aplicación del modelo se restringe a las cuencas en las que se ha calibrado
sus parámetros (sierra peruana: Cusco, Huancavelica, Junín, Cajamarca)

Es importante tener en cuenta las mencionadas restricciones a fin de garantizar una

buena performance del modelo.