0387-B3
Roberto Pizarro T.; Manuel Soto B.; César Farías D. 1
Se calibraron dos modelos de simulación integral de cuenca, Modelos T y Témez, para las cuencas del río Purapel, y del río Achibueno, en la VII Región del Maule, Chile. La calibración consideró el planteamiento original de los modelos, pero diferenciando parámetros para distintos meses o grupos de meses. La validación, se analizó en base a pruebas estadísticas no paramétricas. Así, los resultados muestran que, en general, un 82% y un 61% de la aportación de las cuencas de los ríos Purapel y Achibueno es explicada por los modelos Témez y T, respectivamente.
Se recomienda el modelo Témez de 8 parámetros para la cuenca del río Purapel y el modelo T de 24 parámetros para la cuenca del río Achibueno, como herramientas para la evaluación del recurso hídrico, ya que incluyen los principales procesos del ciclo hidrológico en su simulación y presentan un número de parámetros manejable.
Two Rainfall-Runoff models, Témez and T Model, were calibrated, for Purapel river basin and Achibueno river basin in the 7° region of Chile. The models were calibrated based in original parameters, however some modifications were included, clusting months for new parameters. Models were proved with non parametrics fitness of fit. So, results show 82% and 61% of runoff level in Purapel river basin and Achibueno river basin are explained by Témez model and T model respectively.
Finally, is suggested the use of Témez model for Purapel and T Model for Achibueno.
En la VII Región del Maule, zona donde habitualmente se presenta la problemática del recurso hídrico, surge la necesidad de llevar a cabo una gestión adecuada de éste, considerando el crecimiento demográfico y la oferta natural de agua. Para ello, se hace necesaria la evaluación de las cuencas de la Región y la comprensión de los fenómenos hídricos al interior de éstas. En dicho proceso, los modelos matemáticos son herramientas fundamentales.
Dentro de los modelos matemáticos, los modelos de simulación integral de cuenca, son capaces de generar series de datos de aportación, a partir de datos meteorológicos. Dos conocidos modelos de simulación de este tipo, son el Modelo T (Carmona et al, 1990), y el Modelo Témez (Estrela, 1988 y Ruiz et al, 1996).
Así, esta investigación pretendió contribuir al conocimiento del Ciclo Hidrológico en las cuencas de los ríos Purapel y Achibueno, a partir de la calibración y validación, a nivel mensual, de los modelos precipitación escorrentía T y Témez.
Dentro de las fases metodológicas desarrolladas en la presente investigación, destacan la recopilación de la información (bibliográfica, cartográfica, meteorológica y fluviométrica), el proceso de la información (completación de datos faltantes y obtención de precipitaciones areales y aportaciones específicas) y principalmente, la Modelación matemática que contempló:
Calibración: Se programaron y calibraron los modelos sobre una planilla electrónica, efectuándose el ajuste de los parámetros minimizando la sumatoria del cuadrado de las diferencias entre las aportaciones reales y las aportaciones calculadas, empleando 2/3 de la serie de datos disponibles. El problema de optimización, resultó ser el siguiente:
Sea
Ac depende del valor de los parámetros dada la función:
Luego, Función Objetivo:
Sujeto a:
Validación: En esta fase se evaluó el ajuste brindado por los parámetros. Se observó que la variable aportación no mostró una distribución Normal, y dado que los resultados de pruebas que exigen este supuesto fundamental pierden toda consistencia estadística, se emplearon pruebas no paramétricas: Prueba U de Mann-Withney (Mason, 1993), y R2 (Llamas, 1993), Prueba de Kruskal-Wallis K-W (Taucher, 1997), en la cuenca del río Achibueno y Prueba Kolmogorov-Smirnov K-S, (Stell et al, 1980), en la cuenca del río Purapel.
Variantes introducidas a los modelos: Se elaboraron modificaciones al planteamiento original diferenciándose los parámetros para distintos meses o grupos de meses, flexibilizando así la capacidad de ajuste de los modelos.
Modelo T: El modelo T fue propuesto originalmente en 1955; asume en el suelo de la cuenca una capacidad de almacenamiento de humedad (_). Dada una precipitación (Pi) y una Evapotranspiración Potencial (ETPi) el almacenamiento de humedad al final del periodo viene dado por:
Se producirá escorrentía (Qi) cuando la humedad total en el suelo exceda su capacidad de almacenamiento, siendo este excedente (T), el que se expresa por:
Parte de este excedente constituye escorrentía durante el periodo siendo modelado por un parámetro (Φ), el cual expresa la proporción de agua que permanece en la cuenca. Así, la escorrentía para el período es:
y la humedad excedente al final del mes es:
El modelo así planteado contempla el ajuste de dos parámetros Φ y λ.
Modelo Témez: Este modelo, creado por Témez en 1977, supone dividido el perfil del suelo en una zona superior insaturada, y otra inferior que se halla completamente saturada, asemejando su comportamiento a un embalse subterráneo que desagua en la red superficial. El agua precipitada (P), se divide en evapotranspiración (ET) y excedente (T), cuyo cálculo se efectúa según la siguiente ley:
donde:
La humedad en el suelo al final del mes es,
con una evapotranspiración real igual a:
La ley de infiltración al acuífero es:
Esta infiltración se convierte en recarga (R) al acuífero, mientras la otra parte de T va a formar parte de la escorrentía superficial o directa.
Las ecuaciones que modelan el comportamiento del acuífero son:
;
El volumen de agua almacenado en el acuífero es:
La aportación total se calcula como la aportación superficial más la aportación subterránea. Así, el modelo contempla el ajuste de 4 parámetros Hmax, C, Imáx y alfa.
Los cuadros 1 y 2 muestran los resultados obtenidos en la calibración de los modelos Témez y T, respectivamente, para las cuencas del río Purapel y Achibueno. Durante la calibración se encontró el problema de existencia de abundantes mínimos locales, viéndose afectado el algoritmo de calibración por lo que se conoce como regiones de atracción múltiple (Estrela, 1988; Thyer et al, 1999). Así, se optó por emplear numerosas combinaciones de valores iniciales para los parámetros, seleccionando aquella cuyo valor en su función objetivo resultó mínima.
Cuadro 1. Parámetros Modelo Témez, Cuenca Río Purapel.
Parámetro |
Hmaxh |
Hmaxs |
Ch |
Cs |
Valor |
339,2 |
293,0 |
0,01 |
0,67 |
Parámetro |
Imaxh |
Imaxs |
αh |
αs |
Valor |
1000 |
20,5 |
0,01354 |
0,00891 |
Cuadro 2. Parámetros Modelo T, Cuenca Río Achibueno.
Parámetro |
Abr |
May |
Jun |
Jul |
Φ |
317,17 |
103,04 |
228,41 |
421,75 |
λ |
0,0001 |
0,6641 |
0,4631 |
0,3333 |
Parámetro |
Agoo |
Sep |
Oct |
Nov |
Φ |
254,97 |
126,05 |
1,0 |
33,26 |
λ |
0,5666 |
0,6684 |
0,6883 |
0,6050 |
Parámetro |
Dic |
Ene |
Feb |
Mar |
Φ |
41,59 |
600 |
188,93 |
600 |
λ |
0,4778 |
0,4410 |
0,4861 |
0,4780 |
Las figuras 1 y 2 muestran las gráficas de aportaciones reales y simuladas para las cuencas de los ríos Purapel y Achibueno, construidas en cada caso para el modelo que presentó el mejor ajuste. Allí se aprecia la relación que guardan las aportaciones reales respecto de las simuladas y como estas últimas se asemejan presentando una fluctuación estacional.
Figura 1. Aportaciones específicas reales y simuladas, Modelo Témez, cuenca Río Purapel.
Figura 2. Aportaciones específicas reales y simuladas, Modelo T, cuenca Río Achibueno.
El cuadro 3, se muestra los resultados de las pruebas estadísticas sobre el modelo Témez en la cuenca del río Purapel. Se aprecia que, en términos generales, este modelo explica un 82% de la variabilidad mensual de la aportación, además de entregar resultados satisfactorios dada la aceptación de las pruebas de U de Mann-Whitney y de K-S.
Cuadro 3. Test de Bondad de Ajuste, Modelo Témez, Cuenca del río Purapel.
Coef, Det, (R2) |
U de Mann - Whitney |
Kolmogorov - Smirnov | ||||||
N° Datos |
Calib, |
Valid, |
Estadístico |
Ut |
Conclusión |
Estadístico |
K-S d(95%) |
Conclusión |
136 |
0,77 |
0,82 |
0,13 |
±1,96 |
Se Acepta Ho |
0,10 |
0,16 |
Se Acepta Ho |
En el cuadro 4, se muestran los resultados de las pruebas de bondad de ajuste del modelo T sobre la cuenca del río Achibueno. Tanto el test de comparación de medias U de Mann-Whitney, como el test de comparación de varianzas Kruskal-Wallis, arrojan como conclusión que no existen diferencias significativas entre la aportación real y la simulada. Por su parte, el valor de R2 expresa que sólo un 61% de la variabilidad de la aportación simulada, estaría siendo explicada por el modelo.
Cuadro 4. Test de Bondad de Ajuste, Modelo T, Cuenca del río Achibueno.
Coef, Det, (R2) |
U de Mann - Whitney |
Kruskal - Wallis | ||||||
N° Datos |
Calib, |
Valid, |
Estadístico |
Ut |
Conclusión |
Estadístico |
X t(0,05) |
Conclusión |
192 |
0,71 |
0,61 |
0,48 |
±1,96 |
Se Acepta Ho |
1,43 |
3,84 |
Se Acepta Ho |
Los autores agradecen al proyecto FONDECYT 1010590, la posibilidad de llevar a cabo esta investigación.
Carmona A., Ayuso J. y Ayuso J.L.(1990). "Estudio de algunos Modelos Determinísticos Lluvia-Escorrentía. Aplicación a una cuenca Real". Centro de Estudios Hidrográficos, Ministerio de Obras Públicas Transportes y Medio Ambiente. Ingeniería civil Nº 77. Madrid, España. pp. 13-28.
Dirección General De Aguas. (1998). "Apuntes del curso de capacitación en modelación hidrológica", Ministerio de Obras Públicas, Dirección General de Aguas. Santiago, Chile. 120 p.
Estrela, T. (1988). " Los Modelos de Simulación Integral de Cuenca y su utilización en estudios de Recursos Hídricos". Centro de Estudios Hidrográficos, Ministerio de Obras Públicas Transportes y Medio Ambiente. Ingeniería Civil Nº 72. Madrid, España. pp.83-95.
Llamas, J. (1993). "Hidrología General, Principios y Aplicaciones"., Servicio Editorial Universidad del País Vasco. España. 635 p.
Mason R. y Lind D. (1995). "Estadística para administración y economía". The University of Toledo, Ohio. U.S.A. Editorial Alfaomega. México. 911 p.
Ruiz, J. M., Estrela, T. y Quintas, L. (1996). "Modelización hidrológica distribuida en el proyecto Guadiana. El Modelo SIMPA 1.0 (Sistema integrado modelación Precipitación-Aportación)".
Stell R. and Torrie J. (1980). "Principles and procedures of statistics". McGraw-Hill Book Company. New York. U.S.A. 633 p.
Taucher, E. (1997). "Bioestadística". Vicerrectoría de Asuntos Académicos y estudiantiles. Universidad de Chile. Editorial Universitaria. Santiago, Chile. 310 p.
Thyer M., Kuczera G. and Bates B. (1999). "Probabilistic optimization for conceptual rainfall-runoff models: A comparison of the shuffled complex evolution and simulated abbealing algorithms". Water resources research. Vol 35, Nº3. Marzo 1999. Pag. 767-773.
1 Depto. de Gestión Forestal y Ambiental, Universidad de Talca. Casilla 721, Talca, Chile. E-mail: [email protected], [email protected]