WEAP

Water Evaluation And Planning System

Calibración
Básica
De Caudales
Un tutorial para
Creando el Modelo Basico.............................................................2
Corriendo el Modelo y Cambiando Parámetros Fisicos...........................................................13

April 2014

Stockholm Environment Institute Abril 2014

Creando el Modelo Basico

1. ABRIENDO EL MODELO Y ESCOGIENDO EL ÁREA DE INTERES

Abra el modelo llamada “Piura_Paramos” y una vista de todo el mundo debe aparecer.
Haga clic derecho sobre el mapa y escoja “Set Area Boundaries”. Definir un rectángulo (que
aparece de color verde) sobre el área general de Alto Piura. Puede usar el “zoom” mueva
hacia la izquierda de la pantalla para obtener resolución más clara:

Haga clic OK cuando el área está definido. Vamos a enfocar en los partes altos de los
ríos Quiroz y Chipillico.

Ahora agregue una capa de SIG al modelo (Add Vector Layer):

2
 *Haga clic derecha*

En el pantalla que aparece, encuentre la carpeta del modelo “Piura_Paramos” dentro WEAP Areas la
carpeta “GIS”, donde existen 3 capas de SIG. Escoja el archivo llamada “Paramos-Quiroz.shp”. Se puede
modificar los aspectos de la capa como el color y borde si usted prefiere:

3
Haga clic OK y repita los pasos con la capa dentro la misma carpeta llamada “Rios-Quiroz_Chipillico” con
líneas azules para esta capa de ríos. El área debe aparecer como se observa en la imagen:

Las polígonas amarillos significan las áreas de paramos en la cuenca Quiroz. Hacer “zoom” al área de los
páramos porque vamos a modelar los caudales en la cabecera del rio Quiroz que fluyen de los páramos.

2. CREANDO UN RIO Y UN “CATCHMENT”

Una vez modificado la vista del modelo para “zoom” al área de los páramos, con las
herramientas que se muestran al lado izquierdo “Río” construya un río que empieza por
debajo de los páramos y continúa sobre la línea de referencia del río Quiroz. Asigne un
nombre “Alto Río Quiroz” al río y presione OK, como se muestra en la imagen:

4
Importante: El río debe fluir de sur a norte-oeste.

Ahora cree un “Catchment” dentro de los páramos que significa un parte del área de páramos que
contribuye a los caudales en la cabecera del río Quiroz. Arrastre un “Catchment” (punto verde) en el
sistema y póngalo dentro del área de páramos al lado de la cabecera del Río Quiroz y llámelo “Paramo”.

Luego ponga una vía de corriente del catchment al río. Para ello arrastre una línea de
“Runoff/Infiltration” que empiece sobre el catchment “Paramos” y termine en el punto inicio del río:

5
3. DEFINIENDO LOS PARAMETROS FISICOS DEL CATCHMENT “PARAMO”
Cambie a la vista “Data” (impresione el botón “Data” debajo del “Schematic” en el lado
izquierda del WEAP) y cuando aparece un menú que dice “Select Catchment Method”,
escoja “Rainfall Runoff (soil moisture method)” y haga clic OK. Esto significa el tipo de
cálculo que WEAP va a usar para producir los valores de corriente y caudal del catchment.
Esta opción modela dos capas de suelo, cada uno con parámetros únicos.

Revise que los parámetros generales de año y paso de tiempo para modelar, que se
encuentran en el menú “General / Years and Time Steps”, tengan la siguiente información:

Current Accounts Year (Año Inicio): 2009

Last Year of Scenarios (Año Final): 2010

Time Steps per Year (Numero de pasos de tiempo cada año): 365

Water Year Start (Inicio año hidrológico): Setiembre

En el árbol de opciones del lado izquierdo, vaya al “Demand Sites and Catchments” y
encuentre el catchment “Paramo”. Haga clic sobre este y mire los de parámetros que se
debe definir al lado derecha. Empiece con la pestaña “Land Use” y “Area” para definir el
área del páramo que se va modelar. Primero, debe cambiar las unidades (“Units”) a “Área
 Hectares” (ha):

6
Una vez definido las unidades, ponga el valor 604 en el espacio al lado “Paramo”:

7
Siga estos pasos para cada parámetro en “Land Use”, definiendo las unidades
apropiadas y valores que se muestran en la tabla:

“Land Use” Parámetro Valor Unidades

Área 604 Ha
Kc 2
Soil Water Capacity (Humedad Máxima del Suelo—Capa 300 mm
Superior)
Deep Water Capacity (Humedad Máxima del Suelo – Capa 210 mm
Inferior)
Runoff Resistance Factor (Factor de Resistencia del Suelo a 2
Corriente)
Root Zone Conductivity (Conductividad del suelo en la zona 240 mm/day
de raíces)
Deep Conductivity (Conductividad del suelo en la zona más 360 mm/day
bajo)

8
Preferred Flow Direction (Dirección de Flujo Preferida) If(PrevTSValue(Interflo
w[m^3])/Area[M^2] *
1000 > 15, 0.75, 0.25)
Initial Z1 (Agua inicio en esta capa del suelo) 15 %
Initial Z2 (Agua inicio en esta capa del suelo) 10 %

4. DEFINIENDO LOS DATOS CLIMATICOS DEL CATCHMENT “PARAMO”

Ahora vaya a la pestaña al lado derecho de “Land Use” llamada “Climate” para ingresar los
datos históricos de clima del catchment “Paramo”. Los datos se encuentran dentro de las
tablas en Excel dentro de la carpeta llamada “Piura_Paramos” en “WEAP Areas” y carpeta
“Datos”. Existen varias maneras de ingresar los datos, en este ejercicio vamos usar las
funciones que hay dentro WEAP para ingresar los datos automáticamente, como la función
“ReadFromFile…”. Ponga estas expresiones para cada variable climático:

Variable Valor
Precipitation (Precipitación) ReadFromFile(datos\NCI_Precip.csv, 4)

Temperature (Temperatura) ReadFromFile(Datos\Temperatura.csv)

Humidity (Humedad del Aire) ReadFromFile(Datos\Humedad_Relativa.csv)

Wind (Viento) ReadFromFile(Datos\Velocidad_Viento.csv)

Cloudiness Fraction (Nubosidad) 0.8

Latitude (Latitud) -4.75

Cada variable debe graficarse de la siguiente manera (Ej. Precipitación)

9
5. CREANDO UN ESTACIÓN HIDROLOGICA E IMPORTANDO DATOS
OBSERVADOS DE CAUDALES

Regrese a la vista “Schematic” y “zoom” otra vez al área del río para creer una estación
hidrológica con datos históricos de caudales. Arrastre un “Streamflow Gauge” sobre el río
en cualquier lugar debajo del ingreso de flujo del catchment “Paramo” y llámelo “Caudales
Históricos”:

10
El “Streamflow gauges” es un punto de referencia para comparar los caudales que el
modelo está produciendo contra los datos reales que fueron observados. Con esta
comparación, vamos a determinar los mejores valores de los parámetros físicos del
catchment “Paramos” para producir caudales más similares a los datos históricas.

Haga clic derecha sobre el Estación (punto azul) para ingresar los datos de caudales
históricas. Vaya al “Edit Data  Streamflow Data”:

11
Para ingresar los datos en tablas en Excel de los caudales históricos ponga esta expresión:

ReadFromFile(datos\NCI_Caudal.csv)/1000

Los caudales deben mostrarse de esta manera:

12
Corriendo el Modelo y Cambiando
Parámetros Físicos

1. ANALISIS DE RESULTADOS PRIMEROS

Para correr el modelo, haga clic el botón llamada “Results” y cuando aparece un mensaje
que dice “Results are out of date, do you want to recalculate now?”, haga clic “Yes” para
empezar a correr el modelo. Cuando el modelo ha corrido encuentre los resultados de
caudal del río Quiroz. Abra el menú de resultados en el parte superior de la pantalla y
busque Streamflow”. También lo puedes ver en “Supply and Resources  River 
Streamflow”. Haga clic “Streamflow”, mostrara los valores de caudales diarios del río
Quiroz periodo 2009-2010.

En el lado derecha, hay una leyenda con las series de tiempo que corresponden a los valores
del modelo para el río Quiroz y también los valores históricos. Queremos comparar el
punto del río en el modelo que corresponde con el punto de la estación con los valores
históricos. Para hacer eso, abra el menú sobre la leyenda y seleccione el segundo opción
llamada “selected Alto Quiroz Nodes and Reaches” y luego seleccione: “1\Catchment
Inflow Node 1” (valores del modelo) y “1\Caudales Historicos” (valores históricos).

13
Haga clic OK y seleccione el opción sobre el grafico que dice “Daily Average” para mostrar
valores diarios promedios por los 2 años de registro.

El gráfico debe mostrar lo siguiente:

14
En este ejemplo, la línea verde significa caudales históricos y la línea roja significa caudales
modelados. Note que las líneas son similares en algunos partes pero no acuerdan en la
mayoría del tiempo.

Otra manera de visualizar la correspondencia entre los dos caudales del modelo es el
gráfico de porciento excedencia que se puede ver con el botón “Percent of Time Exceeded”
como se muestra en la siguiente imagen:

15
Este gráfico muestra el caudal (eje horizontal) y el porcentaje de la diferencia entre la serie
de valores de ambos caudales en el (eje vertical). Note que los valores modelados son
menor que los históricos en los dos vistas de los resultados.

Con la finalidad de lograr un modelo que produce valores de caudal muy similares a los
caudales históricos, se hace necesario modificar algunos parámetros físicos del catchment
“Paramo” hasta lograr una mejor similitud entre ambos caudales, a este proceso se le
denomina “Calibración” que lo veremos a continuación.

2. CAMBIANDO LOS PARAMETROS FISICOS PARA CALIBRAR EL MODELO

Cada cuenca en el mundo es diferente y depende de varias condiciones físicos locales. Por
lo tanto, cada modelo necesita una calibración única con valores de los parámetros físicos
que son razonables y que producen valores dentro de los modelos similares a condiciones
históricos. Sin un modelo calibrado, no se tiene la plena confianza que los resultados
obtenidos estén representando condiciones de la cuenca.

En este caso, necesitamos aumentar los valores modelados del caudal y también lograr un
patrón de caudales que correlacione mejor con los caudales históricos. Con varios
parámetros físicos, muchas veces es muy difícil identificar que parámetro variar, pues para
ello es importante entender el estado de la cuenca y conocer bien la función que cumple
cada parámetro en el modelo.

16
En este ejemplo, vamos a cambiar dos parámetros para lograr un modelo calibrado: Kc
(coeficiente de evapotranspiración de la cobertura vegetal) y Soil Water Capacity
(capacidad de agua en el capa superior del suelo).

1) Kc

Este parámetro es un valor que significa la magnitud relativa de

evapotranspiración de la cobertura vegetal. Un valor alto significa mayor
consumo de agua que sale del sistema por evapotranspiración, y no entra el rio
como corriente (Escurrimiento). Los primeros resultados de modelo están por
debajo de los resultados históricos, por tanto es probable que el valor para este
parámetro sea demasiado alto.

Regrese a la vista “Data” y cambie el valor de Paramo  Land Use  Kc de 2 a

1.5.

Corre el modelo y revise los resultados.

Activar y desactivar el botón “Percent of Time Exceeded” para ver los dos tipos
de resultados que comparen los valores modelados e históricos.

Note que los valores modelados han aumentado y están más cerca a los valores
históricos:

17
Pruebe con y otro valor de Kc menor. Regrese a la vista “Data” y modificar el valor de Kc
de 1.5 a 1.2. Corre el modelo otra vez y analizar los resultados:

18
Ahora la magnitud de los caudales modelados están más cercanas en magnitud a los
valores históricos, pero hay diferencias en el patrón de los caudales modelados. Hay
mucho más variabilidad en los valores históricos que los valores modelados. Para mejorar
eso, necesitamos cambiar un otro parámetro físico.

2) “Soil Water Capacity” (Capacidad de agua dentro la capa superior del suelo)

Este parámetro tiene efectos en la sensibilidad de caudales, porque significa la

cantidad del agua que puede existir dentro el suelo. Un gran valor significa una
gran cantidad de agua que el suelo puede almacenar, como una esponja. Un
valor menor va a producir más variabilidad en los caudales porque los flujos del
rio son afectados directamente por los eventos de lluvia. Para lograr más
variabilidad del caudal, necesitamos disminuir el valor original a este parámetro.

Regrese al modelo y cambiar “Paramo”  Land Use  Soil Water Capacity de

300 a 100. Corre el modelo y analizar los resultados:

19
Estés resultados tienen mayor similitud en magnitud y variabilidad de los resultados
históricos y el modelado, ahora podemos decir que ha mejorado la calibración. Si usted
quiere una mejor similitud y tratar de representar aun mejor el modelo, puede continuar
cambiando y probando otros valores, así como otros parámetros, para este ejercicio lo
dejamos con los valores actuales.

20
3. CREANDO ESCENARIOS DE CAMBIO CLIMATICO
Una vez calibrado el modelo, vamos a producir algunos escenarios de cambio climático
para ver los efectos de un aumento de temperatura en 2010 de 2 °C y de cambios en
precipitación dentro el catchment “Paramos”.

Regrese al “Data” y presione el botón “Manage Scenarios” en el parte superior de la

pantalla para crear nuevos escenarios. Dentro de este nuevo menú, haga clic “Add” para
agregar un nuevo escenario y llamelo “Temperatura +2 C” y haga clic OK.

Note el nuevo escenario que aparece debajo de “Reference”. Ahora haga clic sobre
“Current Accounts” y luego haga clic nuevamente “Add” para crear un segundo nuevo
escenario. Llamelo “Reduccion de Lluvia 15%”. Finalmente, haga clic sobre “Current
Accounts” y “Add” un tercero escenario se llama “Aumento de Temperatura y Reduccion
de Lluvia”. La pantalla de escenarios debe aparecer de esta manera:

21
Asegurar que cada nuevo escenario se basado en “Current Accounts”. Este se indicado en
el menú al lado la lista de escenarios. Por ejemplo “Aumento de Temperatura y Reduccion
de Lluvia is based on: Current Accounts”. Asegurar que cada escenario tiene “Current
Accounts” al lado derecho.

Es muy importante notar que los nuevos escenarios solo aplican en los años después del
año inicial (“Current Accounts Year”), en este caso, todos los años después 2009. Porque
solo son dos años de modelamiento, los nuevos escenarios están aplicando solo a 2010, y

22
cualquier cambio realizamos dentro estos escenarios solo van a afectar los resultados de
2010.

Haga clic “Close” y cambiar al escenario al “Temperatura +2 C” con el menú al lado

izquierda del botón “Manage Scenarios”. Dentro el escenario “Temperatura +2 C”, vamos a
modificar los valores de temperatura, para ello incluimos un aumento de temperatura
absoluta de 2 centígrados cada día. Encuentre Paramo  Climate  Temperature y cambie
el expression a:

ReadFromFile(Datos\Temperatura.csv) + 2

Note el efecto que esta expresión ha generado en el patrón de temperatura en el grafico que
se muestra.

Ahora cambie al escenario llamada “Reducción de Lluvia 15%” y encuentre Paramo 

Climate  Precipitation. Modifique este al:

ReadFromFile(datos\NCI_Precip.csv, 4)*0.85

23
Esta expresión significa que todos los lluvias durante 2010 solo tienen 85% los valores
originales.

Ahora cambie al escenario llamada “Aumento de Temperatura y Reducción de Lluvia” y

en el catchment Paramo  Climate  Temperature y cambie la expresión a:

ReadFromFile(Datos\Temperatura.csv) + 2

También en el catchment Paramo  Climate  Precipitation. Modifique este al:

ReadFromFile(datos\NCI_Precip.csv, 4)*0.85

4. CORRIENDO LOS ESCENARIOS DE CAMBIO CLIMATICO

Corre el modelo y cuando los resultados aparecen, cambie la leyenda para mostrar los
varios escenarios en lugar de los valores modelados e históricos. Del menú sobre la
leyenda, escoge “All Scenarios”. Debe mostrar lo siguiente:

24
El escenario “Reference” significa el escenario original. Para ver los resultados solo del
2010, necesite cambiar los años, que se encuentran en la parte inferior de la pantalla. Haga
clic al lado de “All Years” y escoge del menú que aparece “Selected Years”. De la lista de
años, escoge solo 2010 y haga clic OK:

Ahora los resultados solo están mostrando días de 2010 para cada escenario. Se puede
activar o desactivar las varias líneas que corresponden a los escenarios diferentes para ver
los efectos de los cambios en temperatura y precipitación.

25
Los efectos del aumento de temperatura en caudal son muy pequeñas porque no hay nieve
en esta cuenca y solo hay un poco más de evapotranspiración de la vegetación. Para ver las
diferencias exactas, puede ver los resultados en tabla en el botón “Table” al lado de “Chart”
sobre el grafico:

Los cambios en precipitación tienen muchos efectos en caudal porque la oferta del agua ha
disminuido en 15%.

26
Otro resultado muy útil son las sumas anuales de caudal, que se puede ver con el botón
“Annual Total”:

Gracias por su atención.

27