Simulación Precipitación-Descarga Usando El Modelo MHC
Simulación Precipitación-Descarga Usando El Modelo MHC
Simulación Precipitación-Descarga Usando El Modelo MHC
1. Resumen
En los proyectos hidráulicos se requiere estimar los caudales disponibles en los puntos de
interés, tarea que se ve dificultada por la falta de información o porque la longitud del
registro es corta. Por ello, se utilizan, entre otros procedimientos, modelos de simulación
hidrológica que permiten obtener los caudales superficiales a la salida de la cuenca a
partir de los datos de precipitación.
El modelo que se presenta, llamado modelo MHC, fue desarrollado en base al modelo
Stanford por el Dr. Edilberto Guevara, Ing. Agrícola de la Universidad Nacional Agraria La
Molina, Director General de Postgrado de la Universidad de Carabobo, Venezuela. Este
nuevo modelo ofrece mayor versatilidad para ingresar los datos, además resuelve el
problema de flujo superficial y en canal de forma más sencilla. El modelo MHC es un
modelo tipo determinístico conceptual, que involucra gran parte de los procesos
hidrológicos que tienen lugar en la atmósfera y en las capas del suelo. Por esta razón se
espera que los resultados obtenidos con este modelo sean de mucha utilidad para la
correcta estimación del recurso hídrico en una cuenca con escasa información de
caudales. Para ello, es requisito indispensable lograr la calibración del modelo en la
cuenca utilizando información histórica de buena calidad.
El modelo MHC está dividido en dos fases: la fase suelo y la fase canal, donde la fase
suelo alimenta a la fase canal. En el Gráfico 2.1 se ve la interacción entre estas fases, así
como la relación entre los almacenamientos de humedad de la fase suelo.
Evaporación Precipitación
Detención Infiltración
Superficial
Evapotranspiración
Evapo-
transpiración
Almacenamiento
Almacenamiento
Zona Percolación
Agua
Inferior
Subterránea
Canal
(Tránsito)
La Subcuenca media es de interés para los fines de aplicación del modelo, está
controlada por la estación hidrométrica Santo Domingo y a su vez se ha delimitado en 9
subcuencas. El Cuadro 3.1 presenta la lista de las subcuencas de interés, indicando el
área de la cuenca de drenaje y la altitud media correspondiente.
4. Información hidrometeorológica
Altitud
Estación Periodo Resolución Latitud Longitud
(msnm)
Altitud
Estación Periodo Resolución Latitud Longitud
(msnm)
Santo Domingo 1921-1999 Diaria 11° 23' 77° 03' 697
5. Tratamiento de la información
Posteriormente, dado que se necesita que todas las estaciones pluviométricas tengan
información con resolución diaria en el periodo seleccionado para la calibración y
validación, fue necesario evaluar la correlación simple a nivel mensual entre las
estaciones que tienen sólo información mensual vs aquellas estaciones que tienen
información diaria. Las correlaciones se hicieron mes a mes, de este proceso se
Donde:
Cada una de las 9 subcuencas se modeló en la fase suelo para calcular las escorrentías
producidas por cada una de ellas. En el Cuadro 6.1 se muestra algunos de los valores de
las escorrentías calculas para las diferentes subcuencas.
Escorrentías Carac Vichaycocha Añasmayo Media1 Media2 Media3 Media4 Baños Huataya
d001 2.07 2.07 2.07 2.07 2.07 2.07 2.07 2.07 2.07
d002 1.83 1.83 1.83 1.83 1.83 1.83 1.83 1.83 1.83
d003 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75
d004 2.92 3.50 2.14 2.73 2.42 2.09 1.84 3.02 2.52
d005 2.56 3.31 1.63 1.91 1.91 1.60 1.61 2.60 2.06
d006 1.64 1.58 1.54 1.83 1.61 1.53 1.55 1.69 1.61
d007 1.52 1.51 1.46 1.51 1.50 1.47 1.48 1.52 1.49
d008 1.89 2.94 1.41 1.68 1.44 1.41 1.42 2.16 1.51
d009 2.53 3.12 1.41 1.78 1.80 1.35 1.36 2.51 2.01
d010 2.34 3.60 1.35 2.23 1.52 1.30 1.31 2.71 1.70
d011 2.76 3.29 2.37 2.30 2.50 2.36 1.65 2.75 2.46
d012 1.51 1.30 1.23 1.36 1.41 1.21 1.20 1.42 1.45
d013 1.24 1.25 1.16 1.22 1.20 1.17 1.15 1.26 1.21
d014 1.19 1.21 1.15 2.00 1.15 1.12 1.11 1.45 1.16
d015 1.95 2.25 1.13 1.70 1.44 1.08 1.07 1.99 1.57
d016 3.71 5.41 3.21 2.93 3.18 3.23 1.88 3.92 3.06
d017 1.46 1.94 1.13 1.39 1.20 1.12 1.04 1.61 1.23
d018 1.71 1.89 1.18 1.50 1.39 1.15 1.03 1.73 1.45
d019 1.24 1.05 1.14 1.12 1.24 1.14 1.00 1.17 1.21
d020 1.00 1.01 1.27 1.03 1.16 1.30 1.04 1.03 1.05
Donde:
Cuando el número de Nash se acerca a uno, significa que los valores simulados se
acercan bastante a los datos observados.
Para lograr la calibración se probó con los valores de los parámetros, haciéndolos variar
uno a uno a fin de lograr representar los valores observados en la estación Santo
Domingo. Para esta aplicación particular, debido a la inexistencia de datos de
evapotranspiración y temperatura en la cuenca se decidió considerar también la
evapotranspiración con un factor de corrección de evapotranspiración calculada con el
método de Thornwaithe, a evapotranspiración real. Este factor de corrección fue objeto
de calibración.
Para el caso en estudio, se muestra en el Gráfico 6.1 muestra los caudales simulados vs
observados de la calibración, el número de Nash obtenido es 0.47. Este número
representa la mejor calibración obtenida para el modelo en este caso.
90
80
70
60
Caudal (m3/s)
50
Observado
40 Simulado
30
20
10
0
01-ene-97 01-mar-97 01-may-97 01-jul-97 01-sep-97 01-nov-97
El Gráfico 6.2 muestras los caudales simulados vs los observados de la validación, para
este caso se trabajó con los parámetros determinados en la fase de calibración. El
número de Nash obtenido es 0.43.
140
120
100
Caudal (m3/s)
80
Observado
Simulado
60
40
20
0
01-ene-98 01-mar-98 01-may-98 01-jul-98 01-sep-98 01-nov-98
7. Conclusiones y Recomendaciones
Se ha realizado la calibración y modelación del modelo MHC con el objetivo de simular
adecuadamente los caudales en los periodos de estiaje; sin embargo en las épocas de
avenidas los caudales simulados se acercan a lo que se ha observado.
Para mejorar la calidad de los resultados, se recomienda trabajar con el método de los
polígonos de Thiessen modificados, de esta forma los valores obtenidos de precipitación
serian más representativos.