Modelos hidrológicos

semi-distribuidos
(SOCONT, SAC-SMA, HBV, GR4J,…)
Javier GARCÍA HERNÁNDEZ

SEMINARIO DE FORMACIÓN EN MODELIZACION

HIDROLOGICA E HIDRAULICA CON RS MINERVE

Cusco, 2 de octubre de 2016

1. Introducción
2. Modelos hidrológicos
3. División de sub-cuencas en bandas de altura
4. Estación meteorológica virtual

2
1. Introducción

3
 1. Introducción

Modelo hidrológico semi-distribuido

El modelo semi-distribuido divide la cuenca en diferentes sub-

cuencas.

A continuación, el caudal fluvial de cada una de estas sub-

cuencas se calcula a partir de los flujos totales (infiltración,
escorrentía,…) obtenidos por medio de los modelos
hidrológicos u otros métodos.

Finalmente, estos volúmenes de escorrentía se propagan río

abajo para obtener un cálculo estimado del caudal en la salida
de la cuenca más amplia.

4
2. Modelos hidrológicos

5
 2. Modelos hidrológicos

Snow: “Snow model”

SWMM: “Storm Water Management Model”

GSM: “Glacier Snow Melting”
SOCONT: “SOil CONTribution”
GR4J: “Génie Rural à 4 paramètres Journalier”
HBV: “Hydrologiska Byråns Vattenbalansavdelning»
SAC-SMA: “SACramento Soil Moisture Accounting”

6
 2. Modelos hidrológicos

Snow: “Snow model”

• Precipitación
• Temperatura

• Precipitación equivalente
• Altura de nieve
• Theta

• 6 parámetros
• 2 variables de estado

7
 2. Modelos hidrológicos

Object Name Units Description Regular Range

Asn mm/°C/day Degree-day snowmelt coefficient 0.5 to 20
Critical relative water content of the snow
ThetaCri - pack 0.1
bp s/m Melt coefficient due to liquid precipitation 0.0125
Minimum critical temperature for liquid
Snow Tcp1 °C precipitation 0
Maximum critical temperature for solid
Tcp2 °C precipitation 4
Tcf °C Critical snowmelt temperature 0
HsnowIni m Initial snow height -
ThetaIni - Initial relative water content in the snow pack -

8
 2. Modelos hidrológicos

GSM: “Glacier Snow Melting”

• Precipitación
• Temperatura

• Caudal de fusión de nieve

• Caudal de fusión glaciar
• Caudal total
• Precipitación equivalente
• Altura de nieve

• 11 parámetros
• 4 variables de estado

9
 2. Modelos hidrológicos

Object Name Units Description Regular Range

A m2 Surface of infiltration >0
Asn mm/°C/day Degree-day snowmelt coefficient 0.5 to 20
Critical relative water content of the snow
ThetaCri - pack 0.1
bp s/m Melt coefficient due to liquid precipitation 0.0125
Agl mm/°C/day Degree-day icemelt coefficient 0.5 to 20
Kgl 1/day Release coefficient of icemelt reservoir 0.1 to 5
Ksn 1/day Release coefficient of snowmelt reservoir 0.1 to 5
Minimum critical temperature for liquid
GSM Tcp1 °C precipitation 0
Maximum critical temperature for solid
Tcp2 °C precipitation 4
Tcf °C Critical snowmelt temperature 0
0
Tcg °C Critical glacier melt temperature

HsnowIni m Initial snow height -

QsnowIni m3/s Initial outflow of linear snow reservoir -
QglacierIni m3/s Initial outflow of linear glacier reservoir -
Initial relative water content in the snow
ThetaIni - pack -
 2. Modelos hidrológicos

SWMM: “Storm Water Management Model”

• Precipitación

• Caudal de escorrentía
• Altura de escorrentía

• 3 parámetros
• 1 variable de estado

11
 2. Modelos hidrológicos

Object Name Units Description Regular Range

A m2 Surface of runoff >0
L m Length of the plane >0
SWMM J0 - Runoff slope >0
K m1/3/s Strickler coefficient 0.1 to 90
HIni m Initial water level downstream of the surface -

12
 2. Modelos hidrológicos

SOCONT: “SOil CONTribution”

• Precipitación
• Temperatura
• ETP

• Caudal total
• Caudal de escorrentía
• Caudal de base (infiltración)
• ETR
• Precipitación equivalente
• Altura de nieve
• Saturación
• Intensidad neta
• Intensidad de infiltración
• Altura de escorrentía

• 12 parámetros
• 4 variables de estado

13
 2. Modelos hidrológicos

Regular
Object Name Units Description
Range
A m2 Surface of infiltration >0
Asn mm/°C/day Degree-day snowmelt coefficient 0.5 to 20
ThetaCri - Critical relative water content of the snow pack 0.1
bp s/m Melt coefficient due to liquid precipitation 0.0125
Minimum critical temperature for liquid
Tcp1 °C precipitation 0
Maximum critical temperature for solid
Tcp2 °C precipitation 4
Tcf °C Critical snowmelt temperature 0
HGR3Max m Maximum height of infiltration reservoir 0 to 2
SOCONT 0.00025 to
KGR3 1/s Release coefficient of infiltration reservoir
0.1
L m Length of the plane >0
J0 - Runoff slope >0
Kr m1/3/s Strickler coefficient 0.1 to 90
HsnowIni m Initial snow height -
HGR3Ini m Initial level in infiltration reservoir -
Initial runoff water level downstream of the
HrIni m surface -
ThetaIni - Initial relative water content in the snow pack -
 2. Modelos hidrológicos

HBV: “Hydrologiska Byråns Vattenbalansavdelning”

• Precipitación
• Temperatura
• ETP
• Caudal total
• Caudal de escorrentía
• Caudal du reservorio superior Caudal de
• Caudal du reservorio inferior infiltración
• ETR
• Precipitación equivalente
• Altura de nieve
• Humedad del suelo
• Intensidad de recarga
• Altura de agua del tanque superior
• Intensidad de percolación
• Altura de agua del tanque inferior

• 15 parámetros
• 5 variables de estado
15
 2. Modelos hidrológicos

HBV: “Hydrologiska Byråns Vattenbalansavdelning”

Object Name Units Description Regular Range
A m2 Surface of the basin >0
CFMax mm/°C/day Melting factor 0.5 to 20
CFR - Refreezing factor 0.05
CWH - Critical relative water content of the snow pack 0.1
TT °C Threshold temperature of rain/snow 0 to 3
TTInt °C Temperature interval for rain/snow mixing 0 to 3
TTSM °C Threshold temperature for snow melt 0
Beta - Model parameter (shape coefficient) 1 to 5
FC m Maximum soil storage capacity 0.050 to 0.65
PWP m Soil permanent wilting point 0.030 to 0.65
HBV
SUMax m Upper reservoir water level threshold 0 to 0.10
Kr 1/day Near surface flow storage coefficient 0.05 to 0.5
Ku 1/day Interflow storage coefficient 0.01 to 0.4
Kl 1/day Baseflow storage coefficient 0 to 0.15
Kperc 1/day Percolation storage coefficient 0 to 2
HsnowIni m Initial snow height -
WHIni - Initial relative water content in the snow pack -
HIni m Initial humidity -
SUIni m Initial upper reservoir water level -
SLIni m Initial lower reservoir water level -
 2. Modelos hidrológicos

GR4J: “Génie Rural à 4 paramètres Journalier”

• Precipitación
• ETP

• Caudal total
• Caudal de escorrentía
• Caudal de base
• ETR
• Intensidad de percolación
• Intensidad hacia el tanque de propagación
• Intensidad del intercambio subterráneo de agua
• Altura de agua en el tanque de propagación
• Altura de agua en el tanque de producción

• 5 parámetros
• 2 variables de estado

17
 2. Modelos hidrológicos

GR4J: “Génie Rural à 4 paramètres Journalier”

Object Name Units Description Regular Range

A m2 Surface of the basin >0
X1 m Capacity of production store 0.1 to 1.2
X2 m Water exchange coefficient -0.005 to 0.003
GR4J X3 m Capacity of routing store 0.02 to 0.3
X4 day UH time base 1.1 to 2.9
SIni m Initial water content in the production reservoir -
RIni m Initial water level in the routing reservoir -

18
 2. Modelos hidrológicos

SAC-SMA: “SACramento Soil Moisture Accounting”

• Precipitación
• ETP

• Caudal total
• Caudal de superficie (de
escorrentía e interflujo)
• Caudal subterráneo (principal y
suplementario)
• Caudal de percolación profunda
• ETR
• Altura de agua de la superficie
impermeable, de la zona superior en
tensión, zona inferior en tensión, zona
superior libre, zona inferior libre, zona
primaria de agua libre

• 17 parámetros
• 6 variables de estado
19
 2. Modelos hidrológicos

SAC-SMA: “SACramento Soil Moisture Accounting”

Object Name Units Description Regular Range
A m2 Surface of the basin >0
Adimp - Maximum fraction of an additional impervious area due to saturation 0 to 0.2
Pctim - Permanent impervious area fraction 0 to 0.05
Riva - Riparian vegetarian area fraction 0 to 0.2
UztwMax m Upper zone tension water capacity 0.01 to 0.15
UzfwMax m Upper zone free water capacity 0.005 to 0.10
Uzk 1/day Interflow depletion rate from the upper zone free water storage 0.10 to 0.75
Zperc - Ratio of maximum and minimum percolation rates 10 to 350
Rexp - Shape parameter of the percolation curve 1 to 4
Pfree - Percolation fraction that goes directly to the lower zone free water storages 0 to0.6
SAC- LztwMax m The lower zone tension water capacity 0.05 to 0.40
LzfpMax m The lower zone primary free water capacity 0.03 to 0.80
SMA
LzfsMax m The lower zone supplemental free water capacity 0.01 to 0.40
Rserv - Fraction of lower zone free water not transferable to lower zone 0 to 1
Lzpk 1/day Depletion rate of the lower zone primary free water storage 0.001 to 0.03
Lzsk 1/day Depletion rate of the lower zone supplemental free water storage 0.02 to 0.3
Side - Ratio of deep percolation from lower zone free water storages 0 to 0.5
AdimIni m Initial tension water content of the Adimp area -
UztwIni m Initial upper zone tension water content -
UzfwIni m Initial upper zone free water content -
LztwIni m Initial lower zone tension water content -
LzfpIni m Initial lower zone free supplemental content -
LzfsIni m Initial lower zone free primary content -
20
 2. Modelos hidrológicos

Modelos hidrológicos existentes en RS MINERVE

Objects Parameters Initial Conditions

SWMM 4 1
GSM 11 4
GR4J 5 2
SOCONT 12 4
HBV 15 5
SAC-SMA 17 6

21
 2. Modelos hidrológicos

GR4J + Snow

Posibilidad de combinar
un modelo GR4J y
un modelo Snow

• 5 + 5 parámetros
• 2 + 2 variables de estado

22
 2. Modelos hidrológicos

SAC-SMA + Snow

Posibilidad de combinar
un modelo SAC y
un modelo Snow

• 17 + 6 parámetros
• 8 + 2 variables de estado

23
 2. Modelos hidrológicos

Modelos hidrológicos o combinaciones comúnmente utilizadas

en RS MINERVE

Objects Parameters Initial Conditions

SWMM 4 1
GSM 11 4
GR4J+Snow 11 4
SOCONT 12 4
HBV 15 5
SAC-SMA + Snow 23 8

24
3. División de sub-cuencas en bandas de altura

25
 3. División de sub-cuencas en bandas de altura

Concepto semi-distribuido
División en sub-cuencas

0 50 km
26
 3. División de sub-cuencas en bandas de altura

Concepto semi-distribuido
División de cada sub-cuenca en bandas altitudinales

27
 3. División de sub-cuencas en bandas de altura

Concepto semi-distribuido
Cálculo del centro de gravedad (X, Y, Z) y de la superficie de
cada banda de altura y cálculo de las características de los ríos
GSM-Glacier
SOCONT
Kinematic river
Junction

28
 3. División de sub-cuencas en bandas de altura

Concepto semi-distribuido
Cálculo del centro de gravedad (X, Y, Z) y de la superficie de
cada banda de altura y cálculo de las características de los ríos
GSM-Glacier
SOCONT
Kinematic river
Junction

29
 3. División de sub-cuencas en bandas de altura

Concepto semi-distribuido
Agrupación de las bandas de altura de cada sub-cuenca en
sub-modelos y creación de las confluencias

Sub-model
Kinematic river
Junction

30
 3. División de sub-cuencas en bandas de altura

Concepto semi-distribuido
Modelo hidrológico listo para ser creado en el programa
de simulación hidrológica

Sub-model
Kinematic river
Junction

31
4. Estación meteorológica virtual

32
 4. Estación meteorológica virtual

Estación Meteorológica Virtual

Centro de
gravedad

33
 4. Estación meteorológica virtual

Estación Meteorológica Virtual

Polígonos de Thiessen
X X
dPs
PS  CoeffPs  ( Pk 
dz
 ( z s  zi )) X
dTs
TS  CoeffT s  (Tk   ( z s  zi ))
dz

dETPs Centro de
ETPS  CoeffETPs  ( ETPk   ( z s  zi ))
dz gravedad

• 9 parámetros
• 0-6 parámetros de calibración

X 34
 4. Estación meteorológica virtual

Estación Meteorológica Virtual

Shepard
Inverso distancia2 X

Centro de
gravedad

• 9 parámetros
• 0-7 parámetros de calibración

X 35
 4. Estación meteorológica virtual

Estación Meteorológica Virtual

Shepard P
n
dPs 1
 Pi 
i 1 dz
 ( z s  zi )  2
d i ,s
PS  CoeffPs  n
1
 2
i 1 d i , s

Centro de
gravedad

• 9 parámetros
• 0-7 parámetros de calibración

36
 4. Estación meteorológica virtual

Estación Meteorológica Virtual

Shepard T
n
dTs 1
T 
i 1
i
dz
 ( z s  zi )  2
d i ,s T
TS  CoeffT s  n
1
 2
i 1 d i , s

Centro de
gravedad

• 9 parámetros
• 0-7 parámetros de calibración

37
 4. Estación meteorológica virtual

Estación Meteorológica Virtual

Shepard ETP
n
dETPs 1
 ETPi 
i 1 dz
 ( z s  zi )  2
d i ,s
ETPS  CoeffETPs  n
1
 2
i 1 d i , s

Centro de
gravedad

• 9 parámetros
• 0-7 parámetros de calibración

38
 4. Estación meteorológica virtual

Estación Meteorológica Virtual

Object Name Units Description Regular Range

X, Y, Z m Coordinates of the virtual station -
Search radius m Search radius of weather stations >0
Minimal number of stations used
No. min. of
- for interpolation (higher priority ≥1
stations
than “Search Radius”)
Gradient P m/s/m Precipitation gradient -
Station Gradient T °C/m Temperature gradient -0.007 to -0.004
Gradient ETP m/s/m Evapotranspiration gradient -
Coeff P - Multiplying correction coefficient 0.5 to 2

CoeffT °C Adding correction coefficient -2 to 2

Coeff ETP - Multiplying correction coefficient 0.5 to 2

39
 4. Estación meteorológica virtual

Evapotranspiración potencial
Diferentes métodos de calculo teórico
Turc
McGuinness
Oudin
Thornthwaite
Penman-Monteith
Wendling
Hamon
Hande
Hargreaves-Samani
Makkink
Priestley-Taylor
Primault

40
 4. Estación meteorológica virtual

Evapotranspiración potencial
Diferentes métodos de calculo teórico
Turc
McGuinness
Oudin
Thornthwaite
Penman-Monteith
Wendling
Hamon
Hande
Hargreaves-Samani
Makkink
Priestley-Taylor
Primault

41
 4. Estación meteorológica virtual

Evapotranspiración potencial
Turc
T
ETP  CoeffETP  K   ( Rg  50) if T  0 K  0.4 if Month  February
T  15
ETP  0 if T  0 K  0.37 if Month February

McGuinness
Rg Ta  5
ETP  CoeffETP   if T  5
   68
ETP  0 if T  5

Oudin
Re T  5
ETP  CoeffETP   if T  5
   100
ETP  0 if T  5

42
 “Remember that all models are wrong; the practical
question is how wrong do they have to be to not be useful”
(Box and Draper, 1987)

43
Gracias por
su atención