ESTIMACIÓN DE LA RECARGA MEDIANTE EL
PROGRAMA HBV (Introducción a los Modelos)
Jorge Jódar Bermúdez
Científico Titular del IGME-CSIC en Zaragoza
Mayo de 2025
Índice
1. Introducción a los Modelos
2. Modelos Agregados
2
1. Introducción a los Modelos
3
1. Introducción a los Modelos
Modelo Conceptual Modelo Numérico
Modelo Conceptual:
Una representación simplificada del funcionamiento
del Acuífero
4
1. Introducción a los Modelos
Modelo Conceptual Modelo Numérico
• Una representación simplificada del
funcionamiento del Acuífero
• Ha de incluir todos los procesos relevantes en el
balance de masas de Agua
… y ha de incluir una estimación previa de la contribución al
balance global de cada uno de esos procesos
5
1. Introducción a los Modelos
Modelo Conceptual Modelo Numérico
• Una representación simplificada del
funcionamiento del Acuífero
• Ha de incluir todos los procesos relevantes en el
balance de masas de Agua
Modelo Numérico
… y ha de incluir una estimación previa de la contribución al
balance global de cada uno de esos procesos, que el modelo
numérico ha de reproducir
6
1. Introducción a los Modelos
Modelo Conceptual Modelo Numérico
El modelo conceptual es un “ente vivo”, que evoluciona
de la mano de la información disponible de su
implementación numérica 7
1. Introducción a los Modelos
Modelo Conceptual → Ecuación de Balance
La variación de almacenamiento ( ∆S ) de agua en el
acuífero es igual a las entradas (𝑄𝑖𝑛 ) menos las salidas
(𝑄𝑖𝑛 ) de agua del acuífero
∆S = 𝑄𝑖𝑛 − 𝑄𝑜𝑢𝑡
Estado estacionario ➔ ∆S = 0 ∀𝑡 ➔ 𝑄𝑖𝑛 = 𝑄𝑜𝑢𝑡
El nivel de agua (h) en
el acuífero no varía con
h el tiempo, aunque salga
agua del acuífero!
8
1. Introducción a los Modelos
Modelo Conceptual → Ecuación de Balance
La variación de almacenamiento ( ∆S ) de agua en el
acuífero es igual a las entradas (𝑄𝑖𝑛 ) menos las salidas
(𝑄𝑖𝑛 ) de agua del acuífero
∆S = 𝑄𝑖𝑛 − 𝑄𝑜𝑢𝑡
Estado estacionario ➔ ∆S = 0 ∀𝑡 ➔ 𝑄𝑖𝑛 = 𝑄𝑜𝑢𝑡
Estado transitorio ➔ ∆S = ∆S 𝑡
∆S = 0 ➔ 𝑄𝑖𝑛 = 𝑄𝑜𝑢𝑡 ➔ el nivel del acuífero no varía
∆S > 0 ➔ 𝑄𝑖𝑛 > 𝑄𝑜𝑢𝑡 ➔ el nivel del acuífero sube
el acuífero se llena
∆S < 0 ➔ 𝑄𝑖𝑛 < 𝑄𝑜𝑢𝑡 ➔ el nivel del acuífero baja
el acuífero se vacía
1. Introducción a los Modelos
Escala Temporal de los Modelos
La elección de la escala espacial del modelo depende del sistema
hidrológico, y los procesos asociados que queremos simular
10
1. Introducción a los Modelos
Discretización Temporal
∆S 𝑥, 𝑡 = 𝑄𝑖𝑛 − 𝑄𝑜𝑢𝑡
Modelo estacionario ➔ ∆S 𝑥 = 0 ∀𝑡 ➔ 𝑄𝑖𝑛 = 𝑄𝑜𝑢𝑡
Modelo transitorio ➔ ∆S = ∆S 𝑥, 𝑡
Integrar temporalmente la ecuación de balance
considerando discretización diaria, mensual, anual, etc…
11
1. Introducción a los Modelos
Escala Espacial de los Modelos
La elección de la escala espacial del modelo depende del sistema
hidrológico a caracterizar, y del ámbito espacial de los procesos
hidrológicos que queremos considerar y su influencia en el balance
12
1. Introducción a los Modelos
Discretización Espacial ➔ Parametrización
Información necesaria y/o deseable:
Geometría + Condiciones de Contorno y Condiciones Iniciales
La parametrización depende de la información disponible, y del
objeto del modelo ➔ La discretización espacial será 3D, 2D, 1D, o
Agregada (0D) 13
1. Introducción a los Modelos
Discretización Espacial
Integración de la ecuación de Flujo ∆S 𝑥, 𝑡 = 𝑄𝑖𝑛 − 𝑄𝑜𝑢𝑡
𝜕ℎ
2D ➔ 𝑆 = ∇ 𝑇∇ℎ + 𝑟
𝜕𝑡
𝜕ℎ
3D ➔ 𝑆𝑠 = ∇ 𝐾∇ℎ + 𝑟
𝜕𝑡
Diferencias Finitas Volúmenes Finitos Elementos Finitos Elementos de Contorno
El método de integración es “invisible” al modelador, pero
viene condicionado por la plataforma en la que el modelador
implementa el modelo numérico 14
1. Introducción a los Modelos
Escala Espacial y Escala Temporal
[Link]
[Link] [Link]/21/3859/2017/[Link]
Parecen independientes, pero no…
15
2. Modelos Agregados (Modelos Conceptuales)
Modelos Conceptuales (i.e., No modelos de IA)
Son modelos que utilizan simplificaciones con sentido físico, preservando la
física del proceso hidrológico, como, por ejemplo:
• Áreas homogéneas (celdas o grupos de celdas)
• Escorrentía superficial hortoniana
• Los flujos horizontales en el interior del suelo
no son tan importantes (i.e., percolación
vertical de la infiltración)
Realizan la transferencia de la escorrentía hacia aguas abajo (traslación o
tránsito, tiempos de concentración, etc.) mediante esquemas simplificados
Existe equifinalidad (i.e., no unicidad) en los parámetros hidrológicos
considerados: Diferentes sets de parámetros generan la misma respuesta del
modelo numérico. ➔ "Houston, we have a problem“ (Jack Swigert dixit)
16
2. Modelos Agregados (Modelos Conceptuales)
Modelos Físicamente Basados (3D, 2D y 1D)
Hydraulic Conductivity Zones Recharge Zones
Zone-1 Zone-3 Zone-5 Zone-1 Zone-3 Zone-5
Zone-2 Zone-4 Zone-6 Zone-2 Zone-4 Zone-6
Modelos con una cierta zonificación paramétrica correspondiente
al ámbito espacial de aplicación de las diferentes características
hidrogeológicas consideradas e el modelo
17
2. Modelos Agregados (Modelos Conceptuales)
Modelos Agregados (0D)
Integración Ec. de Flujo sobre todo el dominio
𝜕ℎ
∆S 𝑡 = න 𝑆 𝑑𝐴 = න ∇ 𝑇∇ℎ + 𝑟 𝑑𝐴 = 𝑄𝑖𝑛 − 𝑄𝑜𝑢𝑡
𝜕𝑡
Ω Ω
El hecho de no tener un dominio espacial no implica que sean
sencillos. Pueden considerar multitud de procesos hidrológicos
con formulaciones. 18
2. Modelos Agregados (Modelos Conceptuales)
Modelos Agregados (0D)
Modelos-0D
∄ Parametrización espacial (a priori)
∄ Condiciones de contorno
∃ Parámetros de leyes constitutivas
∃ Condiciones iniciales (e.g., precalentamiento)
19
2. Modelos Agregados (Modelos Conceptuales)
Modelos Agregados (0D)
Los principales procesos hidrológicos que se incluyen en la mayoría
de los modelos son:
• Precipitación
• Perdidas Iniciales
• Evapotranspiración
• Infiltración (Suelo)
• Percolación (ZNS, Recarga)
• Escorrentía superficial (Directa, Q0)
• Escorrentía subsuperficial (Interflujo, Q1)
• Escorrentía Subterránea (Flujo Base, Q2)
• Escorrentía Total (Qtot= Q0 + Q1 + Q2)
• Traslación flujo; según escala temporal
2. Modelos Agregados (Modelos Conceptuales)
Modelos Agregados (0D)
Simplifican la cuenca a modo de flujos y tanques (almacenamientos)
Témez (Jódar et al., 2017) HBV (Jódar et al., 2014) AWBM (Boughton, 2004)
2. Modelos Agregados (Modelos Conceptuales)
Modelos Agregados (0D)
2. Modelos Agregados (Modelos Conceptuales)
Modelos Agregados (0D)
25