INGENIERÍA HIDRÁULICA Y AMBIENTAL, VOL. XXXII, No. 2, May-Ago 2011, p.

03-09

Modelación hidrológica del río Levisa

con el modelo HEC-HMS

INTRODUCCIÓN Resumen / Abstract

El desarrollo de la Voluntad Hidráulica en Cuba, procla-
mada por el Comandante en Jefe de la Revolución Fidel Dada la gran importancia del estudio de avenidas para
Castro Ruz, alcanza su máxima expresión en el momen- protección de la población y preservación de inversiones
to actual con la construcción del Trasvase Este-Oeste en las obras del trasvase Este-Oeste, se realizó la
(conocido como la obra hidráulica del siglo en Cuba) para modelación matemática de las avenidas para la presa Levisa
el aprovechamiento futuro del potencial hidráulico del no- con el modelo HEC-HMS. Para suplir la falta de informa-
reste montañoso oriental en zonas deficitarias del oeste ción de lluvias máximas, se estableció una elación Lluvia
y cuyas presas presentan también objetivo de protección Máxima del 1% contra Altura cuyo resultado se transformó
contra inundaciones en los ríos del macizo Nipe-Sagua- en hietograma con la Curva Regional Adimensional de
Baracoa. Una de las presas perteneciente a éste conjunto Acrecentamiento de la Lámina y el programa TableCurve2D.
es la presa Levisa en el río homónimo. Se modeló para el 1% de probabilidad en 48 subvariantes.
Los resultados se compararon con la relación regional
En el trópico húmedo expuesto a huracanes es muy Escurrimiento Máximo del 1% contra Área Tributaria esta-
frecuente que los ríos montañosos como éste, con pen- blecida por la autora para la región Nipe-Sagua-Baracoa,
dientes importantes, presenten grandes avenidas como resultando definitiva la combinación de submodelos: Pér-
consecuencia de precipitaciones extremas (huracanes o didas NC, Hietograma Especificado e Hidrógrafo Unitario
bajas presiones estacionarias) que provocan pérdidas de de Clark.
vidas, destruyen poblaciones y objetivos económicos. Por Palabras clave: modelación de avenidas, modelación
ello adquiere gran importancia el estudio más preciso hidrológica, modelo HEC-HMS.
posible de dichas avenidas para el diseño de la presa
Levisa con la estabilidad requerida, de manera que se Given the great importance of floods definition for
proteja a la población aguas abajo de inundaciones natu- population protection and investment preservation in the
rales o por rotura de presa y se preserve la inversión a construction works of the East-West Diversion, mathematical
ejecutar en dicha obra y los beneficios a producir por la modeling of floods was carried out for Levisa dam with
misma sin sobrevalorar su costo económico. HEC-HMS model. Lacking maximum rainfall data, a 1%
Rainfall vs. Height relationship was established and its
La modelación matemática es una de las manifesta- results were transformed into a hyetogram by using the
ciones de la revolución científico-técnica contemporánea Regional Nondimensional Depth Increase Curve and soft-
ampliamente generalizada en países desarrollados. Las ware TableCurve2D. It was modeled for 1% probability
metodologías empleadas en las investigaciones with 48 subvariants. The results were compared with the
hidrológicas han evolucionado en los últimos años gra- 1% Maximum Flow vs. Tributary Area regional relationship,
cias a la generalización del uso de los modelos matemá- established by the author for Nipe-Sagua-Baracoa region,
ticos. El Sistema de Modelación Hidrológica del Centro being definitive the following submodel combination: CN
de Ingeniería Hidrológica del ejército de los EE. UU.: HEC- Losses, Specified Hyetogram and Clark Unitary Hydrograph.
HMS (por sus siglas en Inglés) es un programa de simula- Keywords: flood modeling, hydrological modeling, HEC-HMS
model.

Valentina Estrada Sifontes, Empresa de Investigaciones y Proyectos Hidráulicos de Holguín

e-mail: [email protected]
Adriel Reyes Gómez, Empresa de Investigaciones y Proyectos Hidráulicos de Holguín, e-mail: adriel.re [email protected]

3
Modelación hidrológica del río Levisa con el modelo HEC-HMS

ción hidrológica tipo evento, lineal y semidistribuido, de- Baracoa, y que se caracterizan por la heterogeneidad y
sarrollado para estimar los hidrógrafos de salida en una extrema cortedad en varios casos, de sus períodos de
cuenca o varias subcuencas (caudales máximos y tiem- observación (oscilan entre los 7 y 44 años).
pos al pico) a partir de condiciones extremas de lluvia,
aplicando para ello algunos de los métodos de cálculo de Se contó además con documentos técnicos relaciona-
hidrogramas de diseño, pérdidas por infiltración, flujo base dos con estudios anteriores en el territorio. Como base
y conversión en escorrentía directa (USACE 2000). Ha cartográfica se utilizaron los planos E 1: 25 000 y 50
alcanzado gran popularidad en la práctica hidrológica in- 000.
ternacional pero aún es escasa su aplicación en el país.

Tratamiento de la información hidrometeorológica

CARACTERÍSTICAS DE LA CUENCA E Para la caracterización de la lluvia máxima diaria anual
INFORMACIÓN HIDROMETEORO LÓGICA del 1% de probabilidad (Estrada et al. 2010), requerida
Características de la cuenca para la determinación de la avenida de diseño, primera-
La cuenca del río Levisa se encuentra al este de la mente fueron definidos todos los valores de la lluvia máxi-
provincia de Holguín, en la vertiente norte de la Sierra ma diaria para cada año de observación de los pluviómetros
Cristal, corre de sur a norte y desemboca en el océano señalados. A partir de estas series de observación se
Atlántico, en la bahía homónima, limita al oeste con la definieron los fenómenos extremos más significativos de
cuenca del río Cabonico, al sur con el río Mícara y las la zona por sus valores y extensión territorial, destacán-
cabeceras del río Mayarí y al este con la cuenca del río dose el correspondiente a 1963 (Ciclón Flora), con valo-
Arroyo Blanco y afluentes derechos del río Mayarí. El área res entre 238 y 331 mm y que abarcó todo el territorio,
de la cuenca correspondiente al cierre en estudio es de seguido en importancia por los años 2000, 1977, 94 y 60.
91 km2 con una altura media de 500 m SNM, tiene forma
ovalada, con pendientes bastante desarrolladas y laderas Dada la heterogeneidad en las longitudes de las series
muy abruptas, en ocasiones casi verticales. y, sobre todo, en la observación por los diferentes
pluviómetros de los fenómenos extremos más significati-
Las mayores alturas se localizan en el parteaguas don- vos, fue necesario establecer la homogeneidad estadísti-
de sobresale el Pico Cristal con sus 1320 m SNM. El ca de la zona, siendo restablecidos por el método de las
cauce de estiaje es poco profundo, tiene un ancho de 10 a distancias los valores correspondientes al ciclón Flora en
20 m, seguidos en ambos márgenes de terrazas muy los equipos que no poseen observación del mismo dada
pronunciadas y abruptas debido a la proximidad de la se- su significativa influencia en las series de lluvia máxima.
rranía. Los suelos son típicos de la región, formados por Se realizó el análisis para los 11 equipos más cercanos al
arcillas rojas de origen laterítico. La vegetación es unifor- área de estudio determinándose que las series de obser-
me, conformada por bosque de pinos en las terrazas y en vación se agrupan en una zona estadísticamente homo-
las márgenes crecen arbustos y matorrales. En la cuenca génea (ZEH) con respecto a las medias y a las varianzas,
existen zonas muy afectadas por la explotación minera. excluyendo 2 equipos por su cortedad y mala calidad.

Se realizó el tratamiento estadístico de las series de

Disponibilidad de información lluvia máxima aplicándoseles 13 curvas teóricas de pro-
En la cuenca en estudio no existe información babilidad a través de los programas PQMAX y AFMULTY
pluviométrica, en general en la zona se cuenta con 18 con utilización del coeficiente de Mamedov, determinán-
pluviómetros. Las series de observación de las precipita- dose la curva de mejor ajuste gráfico en cada caso, re-
ciones están caracterizadas por la heterogeneidad de sultando en la mayoría la Exponencial. El coeficiente de
sus períodos (oscilan entre los 10 y 58 años) y de su Mamedov (k) permite corregir la probabilidad de los fenó-
distribución territorial, quedando grandes fajas altimétricas menos extraordinarios al asignar a los valores atípicos
sin observación, lo que dificulta la definición de los una probabilidad de ocurrencia más cercana a la realidad
estadígrafos de la presa Levisa. que la correspondiente a la que obtienen a partir de las
series de observación, por lo general cortas.
Solamente se cuenta con observaciones del
escurrimiento medio en la cuenca Levisa, aguas abajo del
cierre de estudio, las mismas son extremadamente cor-
X
k (1)
tas, existiendo sólo un año completo de observaciones: Xm
1965, y un período con algunas observaciones mensua-
les y no se cuenta con observaciones del escurrimiento
máximo. En consecuencia se utilizaron otras 7 estacio-
nes de la región físico-geográfica del macizo Nipe-Sagua-
INGENIERÍA HIDRÁULICA Y AMBIENTAL, VOL. XXXII, No. 2, May-Ago 2011
4
 Valentina Estrada Sifontes y Adriel Reyes Gómez

Esta lámina máxima fue transformada en el hietograma

m de 1 día de duración con intervalos de 1 h por la Curva
p * 100 (2)
N 1 k3 Adimensional de Acrecentamiento de la Lámina del 1%
de probabilidad para la región montañosa de Oriente, de-
donde: finida por Trusov (1986) y utilizando el programa Tabla-
X- Valor extremo Curve 2D para el ajuste de curvas. La intensidad máxima
Xm- Promedio de la serie incluyendo X igual a 127 mm en una hora obtenida en el mismo fue
p- probabilidad de ocurrencia de X variada de posición desde el inicio al centro, 2/3 y final del
m- número de orden de X en la serie ordenada hietograma para manejar en la modelación diferentes va-
descendentemente riantes de su ocurrencia. Ver figura 2.
N- Número de valores observados en la serie incluyen- Al no poderse desarrollar la calibración y validación del
do X modelo en la cuenca por falta de observaciones
hidrométricas, sus resultados tendrán un carácter aproxi-
Se estableció una relación lluvia máxima del 1% de mado. A modo de comparación fue establecida una Rela-
probabilidad contra su altura con todos los equipos proce- ción Regional de dependencia del escurrimiento máximo
sados (ver figura 1) y se determinó el valor correspondien- del 1% de probabilidad con el área de la cuenca para las
te a la altura media de la cuenca igual a 440 mm. estaciones hidrométricas del macizo Nipe-Sagua-Baracoa

Figura 1- Relación Lluvia Máxima 1% - Altura

Figura 2- Variantes de posición de la intensidad máxima

INGENIERÍA HIDRÁULICA Y AMBIENTAL, VOL. XXXII, No. 2, May-Ago 2011
5
Modelación hidrológica del río Levisa con el modelo HEC-HMS

Figura 3- Relación Regional Escurrimiento Máximo1%-Área Tributaria.

Región Nipe-Sagua-Baracoa

y por la cual se obtuvo un gasto del 1% de probabilidad culo, que tiene un carácter lineal y semidistribuido; y posi-
para la presa Levisa igual a 1475 m3/s (Estrada et al. bilita estimar los hidrógrafos de salida en una cuenca o
2010). Ver Figura 3 varias subcuencas a partir de condiciones extremas de
lluvias, aplicando para ello algunos de los métodos de
Para ello fue necesario, a partir de los valores del gas- cálculo de hidrogramas de diseño, pérdidas por infiltra-
to máximo del ciclón Flora determinado por las huellas de ción, flujo base y conversión en escorrentía directa
la avenida para algunas estaciones (Bulat 1964), definir el (USACE 2000).
valor aproximado correspondiente al resto en función del
área de la cuenca y agregar el mismo a las series de El mismo cuenta con numerosas ventajas: se basa en
observación. En el tratamiento estadístico de las series métodos hidrológicos para gastos máximos en cuencas
de escurrimiento máximo anual de dichas estaciones, con de cualquier tamaño, permite la modelación de diferentes
aplicación del coeficiente de Mamedov (1), se utilizaron situaciones, es capaz de procesar una voluminosa infor-
los programas AFMULTY y PQMAX, con un total de 15 mación, posee buenas funciones de pre y post procesa-
curvas probabilísticas para valores extremos, de las cua- miento por lo que los datos son más fáciles de editar,
les fueron escogidas la de mejor ajuste o el promedio de modificar y visualizar en pantalla, tiene capacidad para
un grupo de ellas para cada estación. interactuar con modelos digitales, es un programa valida-
do por la práctica internacional, es relativamente sencillo
En el caso de las series más cortas (7 y 8 años) se y fácil de utilizar y es de libre adquisición.
siguieron las recomendaciones de Semionov (1969) para
el tratamiento estadístico de valores mensuales para se- Presenta algunas desventajas relacionadas con las
ries cortas de observación a partir de la definición de los condiciones de su aplicación, como son la insuficiente
períodos de acuosidad interanual y cuya fórmula de posi- información de precipitaciones en la cuenca en estudio y
ción probabilística depende del número de meses por pe- que existe poca experiencia de su aplicación en cuencas
ríodos húmedo y seco. del país.

El programa permite introducir la información necesa-

MODELACIÓN HIDROLÓGICA ria para una simulación, manejar los componentes de aná-
La modelación de la avenida de diseño se realizó con lisis hidrológico a través de módulos integrados y obtener
el modelo HEC-HMS, modelo de simulación hidrológica respuestas gráficas o tabuladas de fácil comprensión e
de amplio uso internacional y poderoso algoritmo de cál- impresión. Trabaja con tres módulos básicos que definen
INGENIERÍA HIDRÁULICA Y AMBIENTAL, VOL. XXXII, No. 2, May-Ago 2011
6
 Valentina Estrada Sifontes y Adriel Reyes Gómez

en su conjunto el proyecto de simulación de la cuenca, micas innecesarias. La fórmula de Kirpich fue calibrada
los cuales deben definirse completamente antes de ini- con base en los datos registrados en 7 cuencas rurales
ciar la corrida de la solución: Módulo Meteorológico, Mó- de Tennesse, en Estados Unidos.
dulo de Cuenca y Módulo de Control.
Las cuencas consideradas presentaban una superficie
Una simulación calcula la transformación de lluvia-cau- inferior a 50 hectáreas y cauces bien definidos, con pen-
dal en el modelo de la cuenca, dada la entrada del modelo dientes que oscilaban entre 3 % y 10 %. Para su utiliza-
meteorológico y las especificaciones de control que defi- ción en cuencas con características diferentes a las del
nen el período de tiempo durante el cual se realizará la proceso de calibrado (Ferrer Polo 1993) es necesario te-
simulación y el intervalo de tiempo a utilizar. Los compo- ner en cuenta lo indicado al respecto por D. A. Chin en su
nentes de los datos de entrada, tales como las series libro "Water Resources Engineering", en el que se inclu-
temporales, tablas y datos son requeridos como yen algunos coeficientes correctores. En el caso de cuen-
parámetros o condiciones de contorno tanto en el modelo cas naturales el mismo es igual a 2.0. Sólo al aplicar di-
de la cuenca como en el meteorológico. cha corrección a la fórmula de Kirpich (lo que comúnmen-
te no se hace) se obtienen valores semejantes al de la
El modelo de la cuenca representa la cuenca física, se fórmula de Temez utilizada en el presente estudio.
desarrolla por el usuario incluyendo y conectando elemen-
tos hidrológicos. Estos últimos usan modelos matemáti- Resultados de la modelación
cos para describir los procesos físicos que se producen La modelación, en dependencia de las variantes utili-
en la cuenca. En la presente aplicación el modelo fue zadas, arrojó valores de gastos máximos del 1 % de pro-
concebido como una cuenca unitaria debido a que no exis- babilidad que varían en 480 m3/s, con un intervalo entre
ten tributarios importantes comparables al río principal ni 1200 y 1680 m3/s y un promedio de 1440 m3/s, luego de
cambios significativos en los regímenes geológico y desechar algunas variantes con valores muy bajos para
climático, la cuenca es abarcable en su totalidad por las las características de humedad y pendientes de la cuen-
lluvias ciclónicas y no existen objetivos internos con in- ca en estudio y que se apartaban del resultado de la Rela-
terés de modelación, de tal forma que para el cálculo de ción Regional Escurrimiento Máximo1%-Área Tributaria
los gastos máximos no es imprescindible descomponer en un 20% o más.
la cuenca. Tomando un 5% de error con respecto al resultado de
la Relación Regional utilizada como comparación, los
La modelación se realizó en 48 variantes y subvariantes mejores resultados se obtuvieron por las siguientes com-
a partir de la combinación de los diferentes submodelos binaciones de submodelos:
que integran el HEC-HMS y en dependencia de la infor-
mación existente en la cuenca, considerando que no existe 1. Modelo de pérdidas: Iniciales y Constantes; modelo
flujo base: meteorológico: Hietograma Especificado con la posición
de la intensidad máxima a 2/3 y al centro del hietograma;
Modelo de pérdidas (considerando una alta humedad y modelo de transformación lluvia-escurrimiento: Hidrógrafo
antecedente): Iniciales y constantes; y Número de la Cur- Unitario de Clark.
va (NC por sus siglas en Inglés) del Servicio de Conser-
vación de Suelos (SCS por sus siglas en Inglés). de EE. 2. Modelo de pérdidas: NC del SCS; modelo meteoro-
UU. lógico: Hietograma Especificado con la posición de la
intensidad máxima a 2/3 y al centro del hietograma; y
Modelo meteorológico: Hietograma Especificado (con modelo de transformación lluvia-escurrimiento: Hidrógrafo
4 subvariantes según la posición de la intensidad máxima Unitario de Clark.
al principio, 2/3, medio o final del hietograma); y Tormenta
Tipo del SCS (en 4 subvariantes según el tipo de tormen- 3. Modelo de pérdidas: NC del SCS; modelo meteoro-
ta). lógico: Tormenta del SCS Tipo 3; y modelo de transforma-
ción lluvia-escurrimiento: Hidrógrafo Unitario del SCS.
 Modelo de transformación lluvia-escurrimiento:
Hidrógrafo Unitario de Clark e Hidrógrafo Unitario del SCS. Es necesario destacar que, para la cuenca en estudio,
la posición de la intensidad máxima a 2/3 o al centro del
Para la determinación del tiempo de concentración hietograma especificado no provoca diferencias significa-
solicitado por el modelo se utilizó la ecuación de Temez tivas en el resultado de la modelación. Como resultado
(Ferrer Polo 1993). Comúnmente en el país se determina definitivo de la modelación se adoptó la tercera combina-
por la fórmula de Kirpich 1940, o de California, que arrojan ción de submodelos cuyo gasto máximo coincide con el
resultados semejantes entre sí, por lo general bajos, lo resultado de la Relación Regional, obteniendo así el
que exacerba el gasto pico y provoca inversiones econó- hidrógrafo de la avenida del 1% de probabilidad: distribu-
INGENIERÍA HIDRÁULICA Y AMBIENTAL, VOL. XXXII, No. 2, May-Ago 2011
7
Modelación hidrológica del río Levisa con el modelo HEC-HMS

Figura 4- Hidrógrafo resultante de la modelación

ción temporal de los gastos, valor extremo y volumen de Recomendaciones

la avenida. Ver figura 4. Utilizar los resultados obtenidos en el diseño de la
presa Levisa.

CONCLUSIONES Y RECOMENDACIONES Utilizar las experiencias obtenidas de la aplicación

Conclusiones del modelo en la modelación hidrológica de otras cuencas
La aplicación del modelo HEC-HMS permitió ampliar nacionales.
la caracterización hidrológica de la cuenca en estudio lo-
grando resultados importantes tales como el gasto máxi- Desarrollar, a partir del modelo confeccionado, el
mo y el hidrograma de la avenida de diseño. estudio de escenarios futuros relacionados con el cambio
climático o transformaciones particulares en la cuenca
Esta simulación, aunque constituye una aproxima- Levisa, así como la ulterior ampliación del propio modelo
ción a la realidad del fenómeno modelado dado que en la para la transformación de las avenidas por el embalse en
cuenca en estudio no existen observaciones para la cali- diseño.
bración y validación del modelo, se considera aceptable
luego de ser comparada con la Relación Regional
Escurrimiento Máximo1%-Área Tributaria para la región REFERENCIAS
Nipe-Sagua-Baracoa. Bulat, B. G. (1964). "Avenidas catastróficas en las pro-
vincias orientales de Cuba del 4 al 8 de octubre de
El trabajo presenta como novedoso la aplicación, en 1963". EIPH. Holguín
una cuenca montañosa de la región oriental del país, de
un modelo hidrológico de amplio uso internacional, la utili- Estrada, V.; Moreno, I. y Cruz, S. (2010). "Trasvase Este-
zación de la ecuación de Temez para la determinación del Oeste. Tramo Sagua - Melones. Presas Levisa,
tiempo de concentración y la confección del modelo Cabonico y Grande. Estudio hidrológico". EIPH. Holguín.
hidrológico de la cuenca tributaria de la presa Levisa, lo
que posibilita el estudio de escenarios futuros. Ferrer Polo, F. J. (1993). "Recomendaciones para el cál

INGENIERÍA HIDRÁULICA Y AMBIENTAL, VOL. XXXII, No. 2, May-Ago 2011

8
 Valentina Estrada Sifontes y Adriel Reyes Gómez

culo hidrometeorológico de avenidas". Centro de Estudios USACE (2000). "Hydrologic Modeling System HEC-HMS".
y Experimentación de Obras Públicas. Madrid. Hydrologic Engineering Center.

Semionov, E. (1969). "Cálculo de los caudales máximos Trusov I. (1986). "Intensidades máximas de las precipita-
de agua por series cortas de observación". INRH. La ciones en las regiones naturales de Cuba". INRH. La
Habana. Habana.

Recibido: enero del 2011

Aprobado: febrero del 2011

INGENIERÍA HIDRÁULICA Y AMBIENTAL, VOL. XXXII, No. 2, May-Ago 2011

9