Presas Calculo
Presas Calculo
Presas Calculo
ISSN: 1390-4272
1
Doctor en Ciencias Técnicas, especialista en Obras Hidrotécnicas, docente a tiempo completo de las
materias de Principios de la Hidráulica y Diseño de Obras Hidrotécnicas. Consultor en diseño y
construcción de proyectos hidráulicos.
2
Ingeniero Geólogo, especialista en geotermia y geotecnia de obras hidráulicas. Docente de Ingeniería
Geológica y Manejo de Cuencas Hidrográficas. Consultor en ingeniería geológica y geotecnia para el
diseño y construcción de Obras Hidráulicas.
Resumen
De manera general, en el Ecuador existe una sensible carencia de información hidrométrica,
indispensable para los diseños de obras hidráulicas. Con miras a paliar esta situación desventajosa, se
presentan algunas ecuaciones lógicas, resultantes de un análisis de masas, para calcular los caudales
mínimos, medios y máximos que, junto con otras, permiten evaluar el coeficiente y el módulo
específico de escorrentía. Lo anterior se ilustra con aplicaciones para varias cuencas del País.
Palabras clave: Caudal mínimo, caudal medio, caudal máximo, caudal ecológico, coeficiente de
escorrentía, Módulo específico de escorrentía.
Abstract
In most cases, in Ecuador there is a noticeable lack of hydrometric information, essential for the
design of hydraulic projects. In order to overcome this disadvantage, have been developed some
logical equations resulting from mass analysis to calculate the minimum, average and maximum
flows, well as others that to assess the runoff coefficient and specific modulus. This is illustrated with
some apps to several watersheds of the country.
1. INTRODUCCIÓN
El diseño de ingeniería de los proyectos de aprovechamiento hídrico requiere de información
hidrológica para su adecuado dimensionamiento. Desafortunadamente, en el Ecuador se dispone
únicamente de datos de precipitaciones y, muy escasamente, de caudales registrados, en vista de lo
cual se dificulta determinar la escorrentía y evaluar la producción hídrica de las cuencas. Para paliar
esta situación, en la práctica, se recurre a modelos empíricos y semi-empíricos desarrollados en otros
países, de acuerdo con sus propias condiciones climáticas, que no son necesariamente aplicables al
territorio ecuatoriano.
1
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
En la literatura científica son muy conocidas las ecuaciones de Creager, Témez, Verni-King, para la
evaluación de caudales de crecida con diferentes períodos de retorno, que no siempre se ajustan a las
condiciones climáticas del Ecuador. El presente trabajo asume que la precipitación media es el
principal parámetro a considerarse para evaluar la escorrentía de una cuenca.
2. DESARROLLO
Determinación de Caudales
Caudal Medio, Coeficiente de Escorrentía y Módulo Específico. A partir de una combinación de
las ecuaciones de Voscresiensky [1], Zhelezniakov, G. Negovskaya, T. y Ovcharov, E. [2] y
Zhivotovsky [3], y considerando un equilibrio de masas, o volúmenes, se propone una ecuación lógica
para evaluar los principales parámetros hidrológicos para el diseño de obras de aprovechamiento
hídrico.
El volumen promedio de precipitaciones anuales en cualquier cuenca, expresado en m3/s, es igual a:
Y el Coeficiente de Escorrentía:
En el Ecuador es más probable encontrar información sobre precipitaciones, por lo que no es difícil
determinar a través de cualquiera de los métodos existentes en la bibliografía, o a partir del Mapa de
Isoyetas elaborado por el Instituto Nacional de Meteorología e Hidrología (INAMHI). También existe
un Mapa de los Módulos Específicos de Escorrentía, elaborado por Pourrut (1995) [4], Figura 1.
En el caso de cuencas con presencia de glaciares o humedales, la fórmula antes propuesta para el
coeficiente de escorrentía debe incorporar la columna de agua, resultante del deshielo o del aporte de
los humedales. Si es que no se toma en cuenta esta columna, el valor del coeficiente de escorrentía
puede resultar mayor a la unidad. Para calcular los caudales por el Método Racional se recomienda
2
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
Caudal Mínimo y Caudal Ecológico. El caudal mínimo de una cuenca se lo utiliza generalmente
como referencia para determinar el caudal ecológico, por lo que se recomienda considerarlo igual a:
3
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
Si se asume que en cada cuenca, la vida animal y vegetal tiene que estar adaptada a las condiciones
naturales del sitio, y que esa adaptación se relaciona directamente con el volumen de agua disponible,
incluyendo las variaciones extremas, se observa que la anterior ecuación propuesta corresponde a las
condiciones mínimas de equilibrio. De acuerdo a lo enunciado, resultaría factible tomar este valor
como el caudal ecológico, sin dejar de lado las variaciones resultantes de condiciones específicas, que
se las determine mediante estudios especializados, o características particulares del uso del agua,
determinadas a partir de factores paisajísticos, de recreación u otros. El caudal mínimo propuesto
corresponde aproximadamente al 97% de probabilidad de ocurrencia de los caudales en la cuenca; en
cualquier caso, el caudal ecológico no supera un valor igual a 2,5 veces este caudal mínimo.
Caudal Máximo. Se asume como premisa que los caudales máximos están siempre relacionados con
un determinado período de retorno en años (T), en vista de lo cual cualquier fórmula empírica debe
involucrar este factor. En segundo lugar, se tiene que el exponente f del área es variable, de acuerdo
con su tamaño.
Para cuencas pequeñas, Voscresiensky (1956) demostró que el exponente f es igual a 1 y, para grandes
cuencas, igual a 0,5, aunque puede variar entre 0,4 y 0,7 [8]. Como ejemplo de lo expuesto, para el
caso de los caudales máximos, se tiene la ecuación de Témez, citada por Robredro (s.f.) [9] o la
envolvente de Creaguer, citada por Martínez et al. [10]. A partir de la poca información disponible en
4
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
el Ecuador sobre los caudales máximos, se pudo establecer que el punto de inflexión de la variación
del exponente f se lo encuentra en cuencas con áreas comprendidas entre 40 y 50 km2, por lo que se
puede establecer que:
A partir del análisis de la información recopilada en el presente trabajo para varias cuencas del
Ecuador, se confirmó, con buen grado de aproximación que; para cuencas con un área A ≤ 45 km2, el
exponente f= 1 y, para cuencas de A > 45 km2, f= 0,5.
Por lo anotado, se confirma que no se debe utilizar una misma ecuación para cuencas de cualquier
área. Esta misma premisa se aplica para la Fórmula Racional, respecto a la cual la mayoría de autores
recomiendan aplicarla solo en cuencas con áreas menores de 50 a 200 km2 y para aplicaciones fuera de
este rango, la ecuación necesita algunas modificaciones como las citadas por Jianpeng et al. [11].
Para el caso del Ecuador, para cuencas con A>45 km2, se propone utilizar la siguiente ecuación para el
cálculo de los caudales máximos correspondientes a diferentes períodos de retorno T:
En la que: a es un coeficiente variable entre 0,01 y 0,1. Con la poca información disponible, que
consta en [12] [13] [14] [15] [16] [17] [18] [19], [20], [21], se determinó que existe una relación
variable entre la precipitación media y el coeficiente a, Figura 2.
Se nota, de manera clara, que la función a=F(P) tiene un mínimo para el valor de 2500 mm. Con el
fin de obtener una mejor correlación, se ha dividido a esta función en dos partes: la primera
corresponde a cuencas con precipitación entre 500 y 2500 mm (para zonas con precipitaciones
menores a 1000, por condiciones específicas de la cuenca, puede variar notablemente el valor de a); y
una segunda parte, para zonas con precipitaciones entre 2500 y 4000 mm, que se la puede ampliar a
5000 mm.
En el caso que exista información de otras cuencas similares (igual franja de precipitaciones), se
recomienda determinar el coeficiente a por analogía, a través de las ecuaciones propuestas.
Se recomienda utilizar a1 =1,8, aunque puede variar entre 1,6 y 2,0 para cuencas con áreas cercanas a
45 km2, como se observa en la Figura 3.
Los principales afluentes son los ríos: Platanales, Sánchez, Cañitas y Coñaque, que, en más del 75%
de su longitud, se desarrollan entre las cotas 100 y 25 msnm. Esto significa que en la cuenca existen
profundos encañonados en medio de laderas muy inclinadas. Las nacientes del Río Grande se ubican
en la cota 460 y, apenas ha recorrido 2 km ya se encuentra en la cota 175, por lo que la pendiente
promedio inicial es i =0,1425; después de 5 km de recorrido se encuentra en la cota 75, (i =0,02). En
los 14,8 km que faltan hasta el sitio de cierre, la pendiente es i =0,0034. La pendiente media de la
cuenca im = 0,0257, con un ancho promedio b=157,12/21,8 =7,207 km.
6
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
Para establecer las precipitaciones sobre la cuenca se utilizaron los datos procesados por el INAMHI,
en el “Estudio Hidrológico de Inundaciones en la Cuenca Alta del Río Chone” elaborado por
Gutiérrez, Góngora y Melo (2008) [22], que es bastante completo y actualizado, porque cubre un
periodo de 42 años (1964-2005), dentro del cual se incluye información de varios eventos El Niño
(ENOS) y de cuyos resultados se consideró el valor de =1309 mm para el análisis del presente
trabajo, que corresponde al obtenido por el método de Kriging.
Coeficiente de escorrentía. A partir del mapa elaborado por Pourrut (1995) [23] se obtuvo el valor
promedio del módulo de escorrentía para la cuenca, que corresponde a M0 = 25 l/s/km2= 0,025
m/s/km2, con lo que resulta que C = 0,602, Q0 = 3,92 m3/s y Qmin = Qeco = 0,206 m3/s.
De acuerdo con los datos presentados por ACOLIT, Cia. Ltda. [24] el volumen medio anual disponible
es de 99,76 hm3, de donde se puede obtener que el caudal medio es 3,16 m3/s, y el caudal ecológico se
estableció en Qeco= 0,200 m3/s.
Existen varias ecuaciones empíricas, muchas de las cuales están citadas por Pérez y Rodríguez [25],
Martínez [26], Krochin [27] y Voscresiensky (1956) [28], de cuya aplicación el rango de caudales
resultantes es desde 88 a 1884 m3/s, dejando mayor incertidumbre que certeza.
Métodos semiempíricos
Se consideran semiempíricos los modelos que, a partir de cierta información experimental, llegan a
inferir resultados como, por ejemplo, los que utilizan las huellas dejadas por las grandes crecidas
históricas, o la información recabada de los habitantes de una determinada área. Según el testimonio
de los habitantes del sector “La Caída”, en el invierno de 1987 la crecida del Río Grande llegó a la
7
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
cota 32 msnm que, a partir de estimaciones realizadas por los autores representa un caudal Q = 695,1
m3/s.
Con base en un análisis expost de la crecida del 4 de marzo de 2013, que se realizó en la sección de
cierre de la presa, los autores evaluaron un caudal de crecida en Q=703 m3/s. Para el caso de zonas
con la presencia de construcciones y cultivos aguas abajo, Martínez, Fernández. y Salas (s.f.) [29], de
la Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación de México
(SAGARPA), recomiendan considerar que los caudales históricos representan el 50% del caudal
máximo, a partir de lo cual se tendría que:
Qmax=1406 m3/s con el valor determinado para la crecida del 4 de marzo de 2013.
Los datos atinentes a cuencas cercanas constituyen un buen referente para la estimación de caudales.
Para el presente caso se dispone de un informe elaborado por Moncayo (1979) [30], en el que se
determina el caudal máximo del Río Jama Qmax= 3847 m3/s, con un área de la cuenca A= 1094 km2. A
partir de este dato se estableció el valor de la constante de la Envolvente de Creaguer en C=66,52; por
lo cual, el caudal máximo de la cuenca del Río Grande sería:
A partir de los datos hidrológicos preliminares del aprovechamiento hidráulico de la cuenca del Río
Jama, elaborados por INASSA (2013) [31], para la cuenca del Río Yescas, afluente del Jama que
colinda con la del Río Grande, con un área A=187,9 km2, se conocen los caudales que aparecen en la
Tabla 3, con los que se puede inferir los del Río Grande, mediante una simple relación de áreas.
Caudal R.
Grande 298.6 431.2 500.0 613.1 828.6 1034.1 1308.8 1430.4
Tabla 3. Caudales para las Cuencas de los Ríos Yescas y Grande
Se destaca que, en los cuatro casos anteriores, el caudal para un período de retorno de 10.000 años es
cercano a 1.400 m3/s.
El valor de k depende de las unidades que se utilice para expresar la Intensidad I (16,67 para mm/min
y 0,278 para mm/h). Para determinar las intensidades se requiere conocer el tiempo de concentración
de la cuenca, para lo cual se pueden utilizar algunas fórmulas citadas por Torres et al. (2012) [32],
Ibáñez et al. (s/f) [33], Garrido (2008) [34] y el Ministerio de Transporte y Comunicación de Perú
(2008) [35], Tabla 4:
8
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
GIANDOTTI
4,48
(1934)
TEMEZ
6,31
(1977)
GOROSHKOV
4,33
(1979)
Al ser notoria la divergencia que presentan los valores del tiempo de concentración, se utilizó el
método de la velocidad del flujo de Goroshkov (1979) [36], como definitorio. Con base en algunos
aforos realizados en 2013 en la sección del Campamento Río Grande, se determinó que v = 1,3 m/s y
al considerar que, en las secciones aguas arriba, la velocidad debe ser ligeramente mayor, se asumió v
=1,4 m/s, a pesar de que para zonas montañosas la bibliografía recomienda tomar valores más altos
que el propuesto. El valor así obtenido resulta bastante cercano al que determina la ecuación de
Giandotti (4,48 horas).
Con base en lo anterior, se asume que el tiempo de concentración es igual a 260 min o 4,33 horas.
En la siguiente Tabla 5 se presentan los valores calculados de las intensidades en minutos y horas.
AUTOR\Periodo
de Retorno 10 25 50 100 500 1000 5000 10000
NADAL 152.4 206.9 260.6 328.3 561.4 707.3 1209.5 1523.8
RAWS 110.9 150.6 189.7 239 408.6 514.8 880.3 1109.1
S/N 154.9 210.2 264.9 333.7 570.6 718.9 1229.3 1548.8
KESTLIN -
NIKOLAY 137.1 186 234.4 295.3 504.9 636.1 1087.7 1370.3
SOKOLOVSKY 123.5 167.6 211.1 266 454.9 573.1 979.9 1234.6
VEL. DE
CRECIDA 148.2 201.1 253.4 319.2 545.9 687.7 1176 1481.6
Tabla 6. Caudales máximos determinados con la Fórmula Racional
Se observa que, entre los autores citados, existe una amplia variación de resultados que se hace más
notoria a medida que se incrementa el período de retorno. Existen otras fórmulas para la evaluación de
caudales, en las que se utiliza como parámetro de cálculo las precipitaciones, como las propuestas por
los siguientes autores:
Fórmula modificada de Verni-King, citada por la compañía ENACON S.A. [37] y el Ministerio de
Obras Públicas de Chile [38].
Donde, P24 – es la precipitación máxima diaria en mm; C(T) – coeficiente de escorrentía que depende
de la zona y del período de retorno. Un coeficiente parecido utilizaba el ex Instituto Ecuatoriano de
Recursos Hídricos (INERHI) para la determinación de caudales, Krochin (1986) [39], que actualizado,
se lo puede estimar como C(T) = 0,3*T0,19 .
Fórmula de Témez
Donde, P24 – es la precipitación máxima diaria en mm; Fs –el factor de reducción por simultaneidad de
la lluvia y es igual a:
A las anteriores se suman las propuestas en el presente trabajo. Como la precipitación media anual de
la Cuenca del Río Grande es de 1309 mm, el caudal se determinó con la ecuación para A> 45 km2 y
el coeficiente a con:
10
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
AUTOR\Periodo de
Retorno 10 25 50 100 500 1000 5000 10000
VERNI-KING 128.6 190.9 249.3 322.8 556.4 695.4 1141.5 1401.8
TÉMEZ 177.4 296.4 401.7 523.6 856.7 1025.1 1472.9 1690.3
SANDOVAL&AGUILERA 162.2 326.8 451.3 575.9 865.1 989.6 1278.8 1403.3
Tabla 7. Caudales máximos con ecuaciones que utilizan precipitaciones
El cálculo de los caudales máximos se lo puede realizar por subcuencas; en este caso, la ecuación a
utilizarse es la desarrollada para áreas menores a 45 km2. Los resultados, para diferentes períodos de
retorno, aparecen en la Tabla 8, que al compararlos con los de la Tabla 7 tienen un buen grado de
aproximación.
Método estadístico
El informe “Actualización de los estudios definitivos del Proyecto de Propósito Múltiple Chone”,
Tomo III, Vol. 1, año 2008, presenta una serie de caudales medidos aforados entre los años 1971 a
1984, que aparecen en la Tabla 9, excepto el de la fila 1, que corresponde al caudal asumido que pasó
en el año 1998 y que alcanzó la cota 32 en el sitio “La Caída”.
A pesar que se trata de una serie corta, puede procesársela para obtener el Coeficiente de Variación
Cv, y establecer un Coeficiente de Asimetría Cs, como aparece en la Tabla 9.
N° 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Caudales,
695 375 265 202 168 150 140 126 118 117 89 84 83 79 72 35
m3/s
Caudal
165 *** Cv 0,99 Cs =2*Cv =1,97=2
Medio
Tabla 9. Procesamiento de los caudales del Río Grande
11
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
Utilizando la Distribución de Pearson III, se establece la curva de Caudales que consta en la Figura 5,
cuyos valores se presentan en la Tabla 10.
Como punto de referencia, en la Tabla 11 se presenta los caudales, establecidos por ACOLIT (2008)
en el diseño final del Proyecto de Propósito Múltiple Chone.
3. CONCLUSIONES
El presente trabajo propone ecuaciones aplicables a cuencas con escasa información hidrométrica,
para determinar los caudales; máximo, medio y mínimo, fundamentados en el equilibrio volumétrico
de las masas de agua. Los caudales máximos corresponden a diferentes períodos de retorno.
Las fórmulas para caudales máximos propuestas se separan para cuencas con áreas menores o iguales
a 45 km2 y mayores de 45 km2, coincidiendo con las recomendaciones para la aplicación de la
Formula Racional. La ecuación de caudales máximos para cuencas pequeñas (A45 km2) es
directamente proporcional al área y a la precipitación media de la cuenca.
Para cuencas grandes se determinó que el coeficiente de proporcionalidad de caudales máximos es
variable en función de la precipitación, con un valor mínimo correspondiente a la precipitación media
12
Revista Ciencia UNEMI-Ecuador. 2014, Vol. 7, N° 2, pp. 100-110. ISSN: 1390-4272
de 2500 mm. Por esta razón otras fórmulas, de otros autores, no son aplicables cuando se obtienen
para otro rango de precipitaciones.
4. REFERENCIAS BILBIOGRÁFICAS
14