Hidrología Estadística en Ingeniería
Hidrología Estadística en Ingeniería
1.1 Introducción
Los estudios hidrológicos requieren del análisis de información hidrometeorológica, esta información
puede ser de datos de precipitación, caudales, temperatura, evaporación, infiltración, etc.
2. PROBABILIDADES
2.1Definición
( )=
1. 0 ≤ P(A) ≤ 1, para todo AϵS (para todo evento A su probabilidad es positiva y cero si el evento es
imposible).
2. P(S)=1
3. P(A1UA2UA3U…UAN)=P(A1+A2+A3+….+AN)=P(A1)+P(A2)+P(A3)+….P(AN).Si
A1+A2+A3+…+AN, es una serie de sucesos mutuamente excluyentes.
Una de las formas de representar las probabilidades de las variables hidrológicas son las funciones de
probabilidad (funciones de densidad), y las funciones de probabilidad acumuladas que a
continuación se mencionan.
Cuando el número n de valores que puede tomar una variable aleatoria X es finito, se dice que la
variable aleatoria X es discreta
A la función y gráfica que asocia una probabilidad a dicha variable aleatoria X se denomina
función de probabilidad discreta f(xi)
Cuando el número de valores n que puede tomar una variable aleatoria X es infinito, se dice que la
variable aleatoria X es continua. Este tipo de variables es más frecuente en hidrología.
La función que asocia una probabilidad a dicha variable se denomina función de probabilidad continua
o función de densidad f(xi). Esta función representa la probabilidad que toma una variable aleatoria X,
la representación gráfica se muestra en la Figura 3
Si X es una variable aleatoria discreta o continua, se define la función de distribución acumulada F(x),
como la probabilidad de que la variable aleatoria X tome cualquier valor menor o igual a x y se designa
por:
F(x)=P(X≤x)
Tal que:
P(X ≤ x) + P(X ≥ x) = 1
En hidrología la variable más frecuente es una variable continua, se analizara la función de distribución
acumulada de esta variable, que está representada por:
( )= ( ≤ )= ( )
( ≤ ≤ ) = ( )− ( ) = ( )
Lo que significa que la probabilidad de un evento a≤x≤b, es igual al área que hay bajo la curva de la
función de densidad f(xi) entre x=a y x=b, ver Figura No 5
0≤F(x)≤1
Es decir que la función de distribución acumulada está en el rango de cero y la unidad o 100%,
dependiendo si se trabaja en porcentajes o decimales.
La Figura No 7 nos permite ver el porcentaje de las observaciones que están por encima (Fxi) o
debajo (1-Fxi) del valor xi con respecto al total.
El Periodo de Retorno T, se define como el tiempo o lapso promedio entre la ocurrencia de un evento
igual o mayor a una magnitud dada, dicho de otra forma, es el intervalo de recurrencia promedio para
un cierto [Link]ísticamente el Periodo de Retorno es la inversa de la probabilidad de excedencia,
es decir:
1
=
( ≥ )
1
=
1− ( ≥ )
Considerar por ejemplo la variable “caudal máximo del año, Q max” para n años.
La gráfica correspondiente para una serie de 41 años será:
La media histórica de esta serie de 41 años resulta 14.9 m3/s. Ahora considerar por ejemplo el valor 20
m3/s. Trazar una recta a 20 m3/s en el gráfico. Realizar el conteo de años transcurridos entre eventos
mayores a 20 m3/s:
Una vez que se presentó el evento “Q>20 m3/s” en el segundo año, transcurrieron 2 años antes de que
se volviera a presentar dicho evento. Luego transcurrieron 5 años, luego 2 años, etc.
Considerando varias centenas de años, el periodo de retorno T será el valor esperado de esos lapsos
de tiempo. Entonces en el ejemplo descrito T puede ser estimado como sigue:
2+5+2+5+6+2+2+1+8+5
= = 3.80 ñ
10
Lo que significa:
Considerando varias centenas de años, el valor de 20 m3/s es excedido en promedio una vez cada 3.8
años, es decir, el periodo de retorno del valor de 20 m3/s es de 3.8 años.
Con otras palabras, en el transcurso de un año cualquiera se tiene una probabilidad de uno en 3.8 (o
sea 26%) de que Q max sea igual o mayor a 20 m3/s.
El periodo de retorno a adoptar para el diseño de una estructura hidráulica debería ser el resultado del
análisis costo-beneficio. A mayor periodo de retorno mayor la obra y en consecuencia más cara y el
beneficio también podría ser más grande. Sin embargo la evaluación de los beneficios es
frecuentemente muy difícil de utilizar, por lo que en la práctica se adoptaran periodos de retorno en
base a la práctica usual. En la Tabla 1, se muestra periodos de retorno recomendados para el
cálculo de caudales de diseño de estructuras menores.
Se dan a conocer otras tablas presentando periodos de retornos recomendados para diferentes tipos
de estructuras civiles: La Tabla No 2 es de carácter general e incluye diversas obras, la Tabla 3 es
exclusivo para obras hidráulicas en carreteras, la Tabla 4 está en función al tipo de área a proteger y
la Tabla 5 en para el diseño de vertederos de embalses.
Por lo común el ingeniero diseña una obra para resistir una avenida de cierta magnitud.
Se define el riesgo de fallo R de un diseño como la probabilidad de que la avenida para la cual se
diseña la obra sea excedida en el transcurso de N años, esto es consider ado como una situación de
riesgo, pues la obra se diseña para soportar cierta avenida máxima, y crecientes mayores podrían
hacerle daño o incluso destruirla, poniendo en riesgo vidas humanas e infraestructuras que están aguas
abajo.
De forma más sencilla se entiende por riesgo de fallo a la probabilidad de que un evento con un
periodo de retorno de T años ocurra al menos una vez en N años.
= ( ≥ ñ )=1− 1− ( ≥ )
1
= ( ≥ ñ )=1− 1−
Dónde:
T = Periodo de Retorno;
N = Años
P( X≥ x) = Probabilidad de excedencia
R = Riesgo de fallo o probabilidad de que un evento con periodo de retorno T años ocurra al menos una
vez en N años.
De la misma manera se puede definir la confiabilidad que viene a ser el complemento del riesgo de fallo,
que se define como la probabilidad de que un evento con periodo de retorno de T años no ocurra en N
años, la confiabilidad se puede expresar de la siguiente manera:
ñ ñ =1− ( ≥ )
1
ñ ñ = 1− = ( )
También es posible calcular el periodo de retorno a partir del riesgo de fallo y del númerode años,
como sigue a continuación:
1
=
ln 1 −
1−
Dónde:
m : Numero de orden.
Una vez seleccionada la fórmula empírica a utilizar, se procede a ordenar los datos de la muestra de
menor a mayor, después se les asigna la probabilidad empírica, que es la probabilidad de no
excedencia. Si se ordena de mayor a menor, la probabilidad asignada será la probabilidad de
excedencia.
b. Ploteo.
El propósito del papel de probabilidad es el de linealizar la relación de probabilidad de tal manera que
los datos graficados se acomoden a una recta, generalmente con fines de comparación. Es una forma
de determinar si una serie de datos está siendo representada de mejor manera por una distribución de
probabilidades en comparación con otras distribuciones de probabilidades teóricas.
Este procedimiento es conocido como la prueba de bondad de ajuste gráfico, que nos sirve para poder determinar
si los datos se ajustan a la distribución representada por el papel de probabilidades. Más adelante se presentarán
las pruebas de bondad de ajuste estadístico
Distribución Normal
Distribución Gumbel
a. Distribución Normal.
También denominada distribución gausiana. Se dice que una variable aleatoria X tiene una distribución
normal, cuando su función de densidad de probabilidad es:
1 1 −
( )= −
√2 2
Dónde:
Para su aplicación lo más fácil es la utilización de una tabla que relacione Z versus f(Z) para lo cual
se ha definido la variable estandarizada como:
−
=
1
( )= −
√2 2
Una característica importante de la distribución normal estándar es que tiene la media cero y la varianza
igual a uno.
1
( )=
√2
O su equivalente:
1
( )= ( )=
√2
Para el cálculo de la función de distribución acumulada se recurre a la tabla de la ley normal que
está en función de la variable estandarizada Z.
La distribución normal es de gran utilidad en hidrología, siendo algunas de sus principales aplicaciones:
Como referencia para comparar varias distribuciones teóricas de ajuste con una distribución empírica.
Para realizar el ajuste se utiliza el papel de probabilidades de la ley normal junto a su recta trazada
analíticamente.
Se dice que una variable aleatoria X tiene una distribución log-normal, cuando su función de densidad
de probabilidad se define como:
Para 0<x<∞
1 1 ln −
( )= −
√2 2
Para -∞<y<∞
1 1 −
( )= −
√2 2
Dónde:
1 1 ln −
( )= ( )= −
√2 2
O su equivalente:
1 1 −
( )= ( )= −
√2 2
Si:
− ln −
= =
1
( )=
√2
Figura No 10 Función de densidad de la distribución Log Normal
Una vez realizada la transformación con la variable estandarizada Z, utilizar las tablas de la ley normal
para el cálculo de la probabilidad o la función acumulada.
Como referencia para comparar varias distribuciones teóricas de ajuste con una
distribución empírica.
Este es una de las distribuciones más utilizadas en hidrología, se dice que una variable aleatoria X,
tiene una distribución Gama o Pearson tipo III, si su función de densidad de probabilidad es:
( − ) ( )
( )=
( )
( − ) ( )
( )=
( )
Dónde:
f(x) : Función de densidad de la variable x. F(x) : Función de distribución acumulada.
x : Variable aleatoria.
x0 : Origen de la variable x, parámetro de posición.
β : Parámetro de escala.
ɣ : Parámetro de forma.
Γ(ɣ): Función gama completa.
Figura No 11 Función de densidad de la distribución Pearson Tipo III
= +
Dónde:
X : Variable analizada, con una probabilidad dada.
Media
∑
=
Desviación Estándar
∑( − )
=
−1
Coeficiente de Asimetría
∑(− )
= =
( − 1)( − 2)
Para determinar el factor de frecuencia, es necesaria la utilización de tablas, para lo cual es necesario
calcular el coeficiente de asimetría y la probabilidad o período de retorno respectivo para la variable
analizada.
En el caso de la distribución log-Pearson tipo III, el procedimiento es el mismo, lo único que cambia
es que se deberá trabajar con los logaritmos de las variables, y se utilizará la misma tabla para
determinar el factor de frecuencia.
La distribución Pearson tipo III es de gran utilidad en hidrología, siendo algunas de sus principales
aplicaciones:
Como referencia para comparar varias distribuciones teóricas de ajuste con una distribución
empírica.
Análisis de errores aleatorios en las observaciones o mediciones hidrológicas.
Para aplicar inferencia estadística.
Para realizar ajustes de distribución empírica de variables hidrológicas de precipitación, caudales,
temperatura, etc., tales como valores anuales, mensuales o valores acumulados anuales, mensuales.
La distribución Gumbel es también llamada distribución de Valores Extremos Tipo I o distribución doble
exponencial. Se dice que una variable aleatoria X tiene una distribución
1
( )=
Dónde:
f(x) : Función de densidad de Gumbel de la variable x.
x : Variable independiente.
α : Parámetro de escala.
μ : Parámetro de posición, llamado moda.
e : Base de logaritmo neperiano.
( )=
Donde F(x) es la función de distribución acumulada de la ley Gumbel. Una forma de calcular α y
μ es con las ecuaciones 10.33 y 10.34 respectivamente, y están en función de los parámetros de la
media ( ) y la desviación estándar (S) de la muestra.
√6
= ≈ 0.78
= − 0.57721
Realizar ajustes de distribución empíricas de variables hidrológicas tales como valores de caudales
máximos anuales, mensuales o precipitaciones máximas anuales, entre otros.
Como referencia para comparar varias distribuciones teóricas de ajuste con una distribución
empírica.
Para efectuar inferencias estadísticas.
Las pruebas estadísticas, tienen por objeto calificar el hecho de suponer que una variable aleatoria, se
distribuya según una cierta función de probabilidades.
Las pruebas de bondad de ajuste grafico más utilizado en hidrología se mencionó en el acápite 1.3 con
la ayuda del papel de probabilidades.
A continuación se detallarán las pruebas de bondad de ajuste estadístico más utilizadas en hidrología
que son:
Chi –Cuadrado
Smirnov – Kolmogorov
Esta prueba es comúnmente usada para verificar la bondad de ajuste de la distribución empírica a
una distribución teórica conocida, fue propuesta por Karl Pearson en 1900.
−
=
= =
Dónde:
θi : Numero de valores observados en el intervalo de clase i. ei : Numero de valores esperados
en el intervalo de clase i.
: Valor calculado de Chi Cuadrado, a partir de los datos.
−
=
Dónde:
Ni : Número de observaciones que caen dentro de los límites de clases ajustadas del intervalo i
N : Tamaño muestral i
Pi : Probabilidad igual para todos los intervalos de clases
Pi =1/k o ei=PiN
Simplificando la última ecuación se obtiene la fórmula computacional desarrollada por Markovic.
= −
TABLA No 7 DISTRIBUCIONES DE X2
Si el Chi cuadrado calculado es menor o igual que el valor tabular, es decir ≤ en tonces se acepta la
hipótesis que el ajjuste es bueno al nivel de significación seleccionado.
Si el Chi cuadrado calculado es mayor que le valor tabular, es decir: > , entonces el ajuste es malo y
se rechaza la hipótesis, siendo necesario probar con otra distribución teórica.
Esta prueba es de fácil aplicación, es válida sólo para ajustes a la distribución normal, en la práctica se
usa para cualquier modelo de ajuste.
b. Prueba de Smirnov-Kolmogorov.
= | ( ) − ( )|
Dónde:
D : Estadístico de Smirnov-Kolmogorov, cuyo valor es igual a la diferencia máxima existente entre
la probabilidad ajustada y la probabilidad empírica
F(x) : Probabilidad de la distribución teórica
P(x) : Probabilidad experimental o empírica de los datos
Si D0 es un valor crítico para un nivel de significación α , se tiene que:
[ | (x) − P(x)| ≥ 0] = ≥ =
También:
< =1−
1. Calcular la probabilidad empírica o experimental Px de los datos, para esto se puede utilizar las
formulas de la Tabla 10.6, de estos el mas recomendado es la fórmula de Weibull, que se indica a
continuación:
=
+1
Dónde:
Px : Probabilidad empírica o experimental. m : Numero de orden.
N : Numero de datos.
2. Calcular la probabilidad teórica F(x), utilizando la ecuación de la función acumulada
F(x) de los modelos teóricos o tablas elaboradas para tal fin.
3. Calcular las diferencias F(x)-P(x)
4. Seleccionar la máxima diferencia: D máx. Ι F(x)-P(x) Ι
5. Calcular el valor crítico del estadístico D, es decir D0, para un nivel de significancia α=0.05 y N
igual al número de datos, la relación para calcular Do es la siguiente:
=
6. Comparar el valor del estadístico D, con el valor crítico D0, con los siguientes criterios de decisión:
Si:
D<D0 El ajuste es bueno, al nivel de significación seleccionado
D≥D0 El ajuste no es bueno, al nivel de significación seleccionado, siendo necesario probar con
otra distribución.
Esta prueba de ajuste no requiere del conocimiento a priori de la función de distribución teórica, es
aplicable a distribuciones de datos no agrupados y de cualquier distribución teórica.
Comparándola con la prueba Chi-cuadrado, no requiere que la frecuencia absoluta de cada clase sea
igual o mayor que 5, esta no es una prueba exacta, sino una prueba aproximada.
Análisis de consistencia de datos
Se propone utilizar el modelo de Water Resources council para la realización de ajustes de series
dudosas, que se alejan significativamente de la tendencia de la información recomendada para detectar
los putnos que se separan de la tendencia de los valores máximos por encima o por debajo y facilitar la
toma de decisiones de retención o eliminación de datos que se puedan afectar significativamente la
magnitud de los parámetros estadísticos con un nivel de significancia del 10%.
Utilizando la ecuación de frecuencia se puede hallar los datos dudosos altos con la siguiente expresión:
= ̅+ ∗
Dónde:
̅ y S variables estadísticas para
XH es el umbral dato dudoso alto en unidades logarítmicas
un tamaño de la muestra. Una ecuación similar se utiliza para datos dudosos
bajos:
= ̅− ∗
PH 10 xH
Precipitación mínima aceptada
PL 10 xL
Utilizando los valores Kn (que contienen los valores para la prueba de datos dudosos de la tabla de U.S.
de Water Resources), teniendo en cuenta el tipo de distribución de los datos analizados. Los picos de
crecientes considerados como datos dudosos bajos o altos se eliminan del registro y se repite el análisis
de frecuencia.
La hidrología es una ciencia de la tierra que estudia las aguas de la tierra: ocurrencia, circulación,
distribución, propiedades físicas químicas y su relación con los seres vivientes.
La tarea de estudiar la hidrología es la de especificar la oferta de agua y sus constituyentes a un
sistemas de recursos de agua, como ser un río, lago, o un sistema de acuíferos y trazar el movimiento
del agua y constituyentes paso por el sistema: uso de agua, control de agua, control de contaminación
de agua.
3.1 Oferta de agua
Según el tipo de aprovechamiento que se desee realizar, los estudios para esta determinación son
diferentes, aunque todos se basan en el estudio de la cuenca hidrográfica y en su balance hídrico.
La hidrología es la ciencia que estudia a la cuenca hidrográfica y los cambios que experimenta en su
balance hídrico a lo largo del tiempo (ciclo hidrológico).
El ciclo hidrológico es el conjunto de procesos involucrados en el intercambio de humedad del mar hacia
la tierra y de ésta nuevamente de regreso al mar, que experimenta el agua en la naturaleza, tanto en su
estado (sólido, líquido o gaseoso) como en su ubicación (ríos, lagos, agua subterránea, nieve, etc).
Entre las fuentes de agua que pueden aprovecharse con fines de riego, se tienen las siguientes:
ríos permanentes
Agua superficial ríos intermitentes
arroyos y quebradas
pozos perforados
pozos excavados
Agua subterránea
manantiales
lagos y lagunas naturales
embalses artificiales
Embalses de agua
estanques de agua (atajados)
i. Cartografía
Imágenes satelitales, fotos aéreas, cartas geográficas IGM escala 1:50.000 y 1:250.000, para la
determinación de la cuenca de aporte.
Se utilizan para conocer la ubicación geográfica del área analizada, su extensión y algunas de sus
principales características físicas. Las cartas 1:250,000 permiten conocer principalmente la ubicación y
vías de comunicación. En proyectos de aprovechamiento es recomendable utilizar cartas 1:50,000.
Una cuenca es una zona de la superficie terrestre en donde (si fuera permeable) las gotas de lluvia que
caen sobe tiende a ser drenadas por el sistema de corrientes hacia un mismo punto de salida
El parte aguas es una línea imaginaria formada por los puntos de mayor nivel topográfico y que separa
la cuenca de las cuencas vecinas
Además de los factores directos que intervienen en el fenómeno de las crecidas, hay que considerar
otros indirectos que influyen notablemente, como son: las dimensiones y la forma de la cuenca vertiente,
el trazado de la red de drenaje, la topografía, el estado de la superficie (suelos, vegetación, superficies
de aguas libres, etc.).
3.6.1 Cuenca
Es el área que contribuye al escurrimiento y que proporciona parte o todo el flujo de la corriente principal
y sus afluentes. Está delimitada por el parte aguas o línea divisoria que define el contorno de la cuenca
hidrográfica y que la separa de las cuencas adyacentes.
Es una línea imaginaria formada por los puntos de mayor nivel topográfico y que separa la cuenca de las
cuencas vecinas.
El área de una cuenca es el área plana en proyección horizontal encerrada por su divisoria.
Hay clasificaciones muy diversas para las cuencas en función de su tamaño. Es importante tener una
clasificación de las cuencas en cuanto a su tamaño pues hay métodos y técnicas de cálculo aplicables
solamente a cuencas pequeñas, y otros, en cuencas grandes. Además el comportamiento del
escurrimiento es diferente (Campos, 1993). En una cuenca pequeña el escurrimiento es influenciado por
las condiciones físicas del suelo y la vegetación. En cambio, para cuencas grandes el efecto del
almacenamiento en el cauce llega a ser pronunciado y hay que darle mayor atención a la hidrología del
río principal.
El tamaño de la cuenca es importante, por ejemplo: Si consideramos dos cuencas de forma similar, una
más grande que la otra. La gota de lluvia que se desplaza del punto más lejos aguas arriba de la cuenca
más grande tiene que recorrer una mayor distancia y, por lo tanto, tarda más en alcanzar la salida dela
cuenca; en cambio la gota de lluvia que provienen del punto más lejos aguas arriba de la cuenca as
pequeña, lógicamente tardará menos tiempo en llegar a la salida.
Es el sistema de cauces por el que fluyen los escurrimientos de la cuenca, de manera temporal o
permanente. Las corrientes que la forman pueden ser clasificadas para reflejar el grado de ramificación
o bifurcación dentro de la cuenca. La clasificación más utilizada es la propuesta por Horton.
Se determina como se muestra en la figura. Una corriente de orden 1 es un tributario sin ramificaciones,
una de orden 2 tiene sólo tributarios de primer orden, etc. Dos corrientes de orden 1 forman una de
orden 2, dos corrientes de orden 3 forman una de orden 4, etc., pero, por ejemplo, una corriente de
orden 2 y una de orden 3 forman otra de orden 3. El orden de una cuenca es el mismo que el de la
corriente principal en su salida; así, por ejemplo, el orden de la cuenca de la figura es 4. Nótese que el
orden de una cuenca depende en mucho de la escala del plano utilizado para su determinación; en este
sentido, las comparaciones entre una cuenca y otra deben hacerse con cuidado, especialmente cuando
los planos correspondientes no están a la misma escala o están editados por diferentes organismos.
Figura N18. Orden de corrientes
Existen varios criterios para determinar la pendiente de la cuenca en función de las curvas de nivel que
pasan por ella, (Springall, 1970), (Campos, 1993). Uno de los más usados es el que considera la
pendiente igual a la longitud total de las curas de nivel dentro de ella, multiplicada por el desnivel
constante entre el tamaño de la cuenca.
a *h i i
Hm i 1
i n
a
i 1
i
El perfil longitudinal del cauce es un gráfico de los valores de sus recorridos horizontales (abscisas)
contra sus cambios de elevaciones respectivas (ordenadas).
La pendiente del cauce principal se estima por diferentes procedimientos, simplificados o no.
Los más utilizados son: los basados en el perfil del río, y el criterio simplificado de dividir el desnivel del
río entre su longitud).
3.6.11 Geología, suelos y cubierta vegetal en la cuenca
El estudio de la geología, rocas, suelos, de la cuenca debe estar encaminado a clasificar en términos
generales “la permeabilidad de los suelos (textura), localización de las aguas subterráneas y áreas de
afluencias y efluencias subterráneas.
Con relación a los suelos será conveniente estudiar la variabilidad y predominio de los suelos en toda la
cuenca, su capacidad de infiltración, su textura y características hidrofísicas de interés para evaluar si
influencia en las pérdidas iniciales de una tormenta.
En cuanto a la cobertura vegetal, es necesario conocer los cultivos predominantes, áreas desérticas y de
bosques, ya que la vegetación controla la acción y movimiento del agua.
Es conveniente tener estudios en la cuenca, de donde obtener el porcentaje de área de cuenca cubierta
por cada tipo de suelo, vegetación y zonas geológicas.
3.7.1 Lluvias
- horarias
- diarias
- mensuales
La Organización Meteorológica Mundial (OMM) recomienda que un pluviómetro cubra como mínimo las
siguientes superficies:
Curvas de descarga
Hidrogramas
- Pozos: pruebas de bombeo
Embalses:
Levantamientos topográficos - batimetría
Relación altura – área y altura – volumen.
3.7.3 Otros datos climatológicos
Geología
Transporte de sedimentos
- Relleno de información
- Generación de datos
Para todos estos estudios se pueden aplicar diferentes metodologías, pero la garantía de los resultados
es una función directa de los datos disponibles. Períodos de registro de 20 años recién se consideran
representativos y no deberían usarse períodos menores a 10 años.
4. MODELACION DE CAUDALES MENSUALES
Se compara la cuenca con otras de características semejantes que cuenten con datos de aforo y de
lluvia.
P = Q + I + E
P = precipitación
Q = escorrentía superficial
I = infiltración subsuperficial y subterránea
E = evaporación directa y evapotranspiración
El parámetro de mayor dificultad de estimación es la infiltración.
Es el método más utilizado para el estudio de cuencas pequeñas en nuestro país, pese a su inexactitud
y grandes limitaciones.
V = C * P * A
Aunque la fórmula parece resolver de manera sencilla el problema de obtener la escorrentía directa V
correspondiente a una lluvia P, la determinación del valor apropiado del coeficiente de escorrentía C es
sumamente compleja.
El valor C varía según las características físicas y topográficas de la cuenca, tipos de suelos, tipo de
cobertura vegetal e incluso la distribución de las lluvias en el tiempo. Los valores que se encuentran en
la literatura deben ser tomados con mucho criterio y en lo posible validados con observaciones directas,
a cargo de ingenieros con experiencia en este campo.
A manera indicativa, se dan a continuación algunos valores de C, tomados del libro “Hidrología para
estudiantes de ingeniería civil”, de W. Chereque (Lima, 1989), aplicables al caso de tormentas y áreas
muy pequeñas.
En este método se usan tres variables para determinar el escurrimiento: la precipitación, la humedad
antecedente y el complejo hidrológico suelo-vegetación.
( P 0 .2 S ) 2
Q
P 0 .8 S
En la que
1000
Número de Curva CN
10 S
Tabla N8. Números de curva
Grupo A, tiene bajo potencial de escorrentía
Consiste en comparar la cuenca en estudio con cuencas con datos de aforos. En México se han
obtenido valores límites del coeficiente de escurrimiento anual en función de la superficie de la cuenca,
de la precipitación y de la cobertura vegetal, los cuales se dan a continuación.
Hasta 800 mm Ce = 0% a 5%
De 800 a 1200 mm Ce = 5% a 15%
De 1200 a 1500 mm Ce = 15% a 35%
Mayor a 1500 mm Ce = 35% a 50%
b. Por fórmulas
Se basan en el tipo de suelo, uso del suelo y la precipitación anual, en forma similar, aunque más
simplificada, al método del número de curva del SCS. Las fórmulas se consideran válidas para valores
de precipitación anual variables entre 350 y 2250 mm.
Suelos muy permeables, tales como arenas profundas y loes poco compactos
Suelos medianamente permeables, tales como arenas de mediana profundidad, loes algo más
compactos que los correspondientes a suelos tipo A y terrenos migajosos.
Suelos casi impermeables, tales como arenas o loes muy delgados sobre una capa impermeable, o bien
arcillas.
Según el tipo de suelo y el uso o cobertura, se obtienen valores de K de la siguiente tabla N9
Tipo de Suelo
Uso o cobertura del suelo
A B C
Cultivos:
En hilera 0.24 0.27 0.30
Legumbres o rotación en pradera 0.24 0.27 0.30
Granos pequeños 0.24 0.27 0.30
Bosque:
Cubierto más del 75% 0.07 0.16 0.24
Cubierto del 50% al 75% 0.12 0.22 0.26
Cubierto del 25% al 50% 0.17 0.26 0.28
Cubierto menos del 25% 0.22 0.28 0.30
Fórmulas:
( P 250)
Ce K *
- Para K < 0.15: 2000
( P 250) ( K 0.15)
Ce K *
- Para K > 0.15 2000 1 .5
Dónde:
Se refiere al uso de aguas de la cuenca que quiere ser aprovechada para otros fines, tales como agua
potable, generación de energía eléctrica o uso en minería.
En este grupo se incluyen también los derechos de terceros, es decir de pobladores de otras áreas
agrícolas que también utilizan las aguas de la cuenca y que deben ser establecidos con exactitud para
evitar conflictos posteriores.
Por último se debe citar la necesidad de dejar un porcentaje de caudal natural de un río como caudal
ecológico. El Ministerio del Medio Ambiente y Agua indica que un 20% del caudal natural del río en
época de estiaje debe ser considerado caudal ecológico.
“Una avenida es un incremento rápido del caudal de un río, que lleva aparejados incrementos anormales
del nivel de aguas, provocando la inundación de terrenos próximos” (Francés, 1994)
“Una avenida es el paso por el río de caudales extraordinarios de gran magnitud” (MOPU, 1987)
“Crecida o avenida de un río es la elevación rápida, generalmente breve, del nivel de las aguas hasta un
máximo a partir del cual desciende más lentamente” (Heras, 1990)
En conclusión, el elemento natural, imprescindible para que se produzca una crecida, es la presencia de
una lluvia máxima o tormenta (Chávez, 1994). De ahí la importancia de estudiar el carácter y el origen
de las lluvias máximas en la zona de estudio.
Las obras que se realizan con el fin de modificar o mitigar los extremos del ciclo hidrológico pueden
resumirse de forma práctica como:
Obras de aprovechamiento
Son los proyectos principalmente constituidos por los embalses que se realizan para modificar el ciclo
hidrológico. Almacenan el agua en periodos húmedos ara poder disponer de ella en los periodos de
déficit.
Obras de control
Son los proyectos cuyo objetivo es brindar protección contra las inundaciones y están constituidos por
embalses, cauces de alivio, diques, bordos de protección, etc.
Obras de drenaje
Son las obras que se construyen para mitigar los efectos de las crecidas de los ríos y que dan paso a las
avenidas y están entre ellas los sistemas de drenaje urbanos y rurales, puentes, alcantarillas, etc.
Para los estudios de obras hidráulicas y prevención de avenidas tales como los aliviaderos de las
presas, puentes, alcantarillas, canales de drenaje, diques, etc., es necesario definir la avenida de diseño
o avenidas de proyecto.
La selección del caudal máximo a admitir para un proyecto hidráulico implica un balance de los
problemas a resolver, los riesgos y vulnerabilidad de la zona a las inundaciones, los daños en vidas
humanas y pérdidas económicas debido a las crecidas, y los costos de las obras, entre otros.
Es por eso que la determinación de la tormenta y avenida de diseño se relaciona con una probabilidad
determinada de ocurrir estos eventos durante un periodo dado y es considerada la probabilidad de
diseño de la obra.
Existen varios procedimientos para seleccionar la probabilidad de diseño de diferentes obras. Hay
métodos analíticos, gráficos y normas técnicas de selección. Usualmente según el tipo y categoría de la
obra, su importancia, los riesgos, daños y consecuencias de un mal funcionamiento y aspectos
económicos entre otros, se seleccionan las probabilidades de diseño. (ver tabla 10)
Drenajes viales 5 – 20 20 – 5
Drenaje urbano 4 – 20 25 – 5
Drenaje agrícola 10 – 25 10 – 4
Diques 1 – 10 100 – 10
Además de los factores directos que intervienen en el fenómeno de las crecidas, hay que considerar
otros indirectos que influyen notablemente, como son: las dimensiones y la forma de la cuenca vertiente,
el trazado de la red de drenaje, la topografía, el estado de la superficie (suelos, vegetación, superficies
de aguas libres, etc.).
La tormenta de diseño se estima por varios procedimientos. Uno de ellos es considerando la tormenta
máxima meteorológicamente posible para una cuenca y un intervalo de tiempo dados, sin tomar en
cuenta los posibles cambios climáticos a largo plazo. Esta sería la tormenta máxima probable TMP, la
cual se utiliza en el diseño de las obras de desagüe de las grandes presas, donde un fallo de la cortina
ocasionaría pérdidas humanas y una gran catástrofe económica.
Para presas pequeñas y otros tipos de obras de menor categoría (canales de desvío, diques, drenaje
urbano y agrícola, etc.) el concepto de tormenta de diseño está relacionado con la probabilidad de
diseño de la obra o probabilidad que se selecciona para definir la tormenta que es capaz de evacuar
dicha obra sin destruirse.
La avenida de diseño se define como la avenida máxima que puede resistir la obra sin dañarse. Y puede
obtenerse probabilísticamente, fijando un determinado nivel de riesgo, o determinísticamente.
Hay otras definiciones para las avenidas o crecidas, como por ejemplo, las dadas por Heras (1980):
Así en cada país se definen normas Técnicas o regulaciones para definir según la obra la selección de la
Crecida de Diseño.
Estas avenidas se pueden caracterizar por toda o una parte de los elementos siguientes:
Estos elementos se definen por diferentes métodos hidrológicos, teniendo en cuenta el origen y el
comportamiento de las crecidas, además de las características físicas de la cuenca y la disponibilidad de
información meteorológica e hidrológica disponible.
Sus resultados se expresan mediante algunas de las formas siguientes: mapas isoyeticos, curvas IDF,
hietogramas o pluviogramas reales o de un periodo de retorno, etc.
Los datos faltantes en los registro de lluvias diarias son bastante frecuentes y se deben a diversas
causas. Para reconstruir estos datos se utilizan varias metodologías que pueden ser revisadas en varias
bibliografías, por lo general se estiman los datos basándose en los registros de estaciones cercanas.
Se recomienda utilizar los métodos de regresión lineal y el de US National Weather Service para deducir
el dato faltante de la lámina de lluvia diaria de una estación a partir de los registros de ese día en las
estaciones pluviométricas vecinas.
El control de la calidad del dato hidrológico consistirá en detectar valores erróneos o dudosos ya sea por
exceso o por defecto. La existencia de los valores dudosos se podrá confirmar por medio de
comparaciones entre los registros de estaciones vecinas. Estos valores de confirman, se corrigen o bien
se eliminan.
Los resultados del procesamiento de las observaciones del régimen de lluvias y escurrimientos de una
cuenca, resultan válidos sólo en el caso en que se empleen series homogéneas. Cualquier estudio
hidrológico debe comenzar con la comprobación de la homogeneidad y aleatoriedad de las series de
observaciones que pretende utilizar.
Existen varios procedimientos para comprobar la calidad de datos de mediciones, los más
recomendados para el análisis de los datos hidrológicos máximos son:
Uno de los objetivos del procesamiento de las series de eventos máximos es cuantificarlos para diversos
periodos de retorno T.
Para la conversión de lluvia diaria a lluvia en 24 horas se utiliza un factor de 1,13 para ajustarlo por
intervalo fijo y único de observación, pues los registro de lluvias maximas diarias, se toman de 8:00 am a
8:00 am del día siguiente y entonces con tal corrección la lluvia representativa se convierte en la lluvia
máxima de 24 horas.
Para estimar probabilidades empíricas, se recomiendan las fórmulas de Weibull y la de Hazen para las
series de valores extremos:
Weibull
,
= Hazen
Dónde:
Las distribuciones de probabilidad adecuadas a los registros de lluvias máximas diarias son básicamente
tres:
Para comprobar cuán bien se adapta una serie empírica a una ley de probabilidades de distribución
deben realizarse las pruebas estadísticas que miden la bondad de ajustes teóricos, entre ellas se tiene
Chi cuadrado y de Smirnoff-Kolmogorof.
5.4 Lluvia máxima
En los estudios hidrológicos de crecidas para definir los parámetros de la avenida de diseño es
necesario definir todos o algunos de los parámetros siguientes:
Para la determinación de la crecida máxima de diseño existen varios métodos entre ellos:
a) Métodos directos
Se utilizan cuando no hay registros de caudales máximos en la cuenca. Existe una gran variedad de
métodos. Tienen como objetivo principal la ventaja de la facilidad de cálculo y como limitaciones,
fundamentalmente, lo simplificado que asumen el proceso de formación de las avenidas y que sus
parámetros están deducidos para cuencas específicas. La mayoría de los métodos relacionan el caudal
máximo con el área de la cuenca y muchos otros con las intensidades máximas de las precipitaciones y
características físicas geográficas de la cuenca. Dan como resultado Qmaxp y en algunos métodos los
hidrogramas sintéticos de avenidas
Este método fue desarrollado por el Soil Conservation Service (SCS) de los estados Unidos; se aplica a
cuencas medianas como pequeñas. El parámetro de mayor importancia de la lluvia generadora, es la
altura de esta, pasando su intensidad a un segundo plano. Su principal aplicación es la estimación de las
cantidades de escurrimiento tanto en el estudio de avenidas máximas, como en el caso del cálculo de
aportaciones líquidas.
d) Métodos estadísticos
Se utilizan en cuencas con registros, adecuados de gastos máximos instantáneos. Estos registros
permiten formar series estadísticas, a las cuales se controla su calidad, y se les hace un análisis de
probabilidad para definir la distribución teórica de probabilidades de mejor ajuste, para evaluar y
determinar el parámetro de diseño Qmaxp.
e) Métodos Hidrológicos
Son los que mejor representan el fenómeno de las crecidas. Necesitan como datos básicos: hietogramas
de las lluvias máximas, características físicas morfométricas de la cuenca y en algunos casos,
hidrogramas de avenidas. Se obtienen como resultados de estos métodos: el hidrograma de la venida y
el Qmax, dependiendo del tipo de lluvia que los provoca. Estos métodos son los que generalmente son
utilizados en la modelación hidrológica determinística o conceptual de las crecidas.
En el presente documento se desarrollara el método del número de curva SCS (Soil Conservation
Service).
5.5.2 Método del número de curva del SCS (Soil Conservation Service)
Para desarrollar el método del número de curva primeramente se pasará a definir algunos conceptos
como:
Se denomina tiempo de concentración tc, al tiempo transcurrido, desde que una gota de agua cae, en el
punto más alejado de la cuenca hasta que llega a la salida de esta. Este tiempo es función de ciertas
características geográficas y topográficas de la cuenca.
El tiempo de concentración debe incluir los escurrimientos sobre terrenos, canales, cunetas y los
recorridos sobre la misma estructura que se diseña.
Todas aquellas características de la cuenca tributaria, tales como dimensiones, pendientes, vegetación,
y otras en menor grado, hacen variar el tiempo de concentración.
A) Fórmulas empíricas
Dónde:
= . .
Dónde:
Ecuación desarrollada por el SCS (Soil Conservation Service) a partir de información de cuencas de uso
agrícola; ha sido adaptada a pequeñas cuencas urbanas con áreas inferiores a 2000 acres.
.
.
−
= .
Dónde:
El nombre del Número de Curva deriva de una serie de curvas, cada una de las cuales lleva el número
de curva CN, que varía de 1 a 100. Un número de curva de 100, indica que toda la lluvia escurre, y un
número de 1, indica que toda la lluvia se infiltra, por los números de curvas representan coeficientes de
escorrentía.
Este método es utilizado para estimar la escorrentía total a partir de datos de precipitación y otros
parámetros de las cuencas de drenaje. El método fue desarrollado utilizando datos de un gran número
de cuencas experimentales, y se basa en la siguiente relación:
Dónde:
= −
−
=
Dónde:
=
+
Reemplazando Pe, se tiene
( − )
=
− +
Los autores del método, por datos experimentales obtuvieron una relación entre Ia y S, la cual es:
( − . )
=
+ .
Dónde:
Esta es la ecuación principal del método. Se debe tener presente que es esta ecuación, P y S deben
tener las mismas unidades y el Qobtenido, también tendrá esas mismas unidades.
El SCS después de estudiar un gran número de cuencas estableció una relación para estimar S a partir
del número de Curva CN, mediante la siguiente ecuación:
=
+
O también:
= −
Dónde:
Figura 29. Duración para el exceso de precipitación para 6 horas de precipitación y para CN 40-100
Figura N30. Familias del Hidrograma
El SCS presenta la tabla de numero de curva para diferentes prácticas agrícolas, diferentes condiciones
hidrológicas y grupo hidrológico de suelos, La tabla fue elaborada para una relación Ia=0.2S y para una
condición de humedad antecedente promedio (CHA II)
a) Condición hidrológica
Grupo de Descripción
suelos
El proceso de cálculo del caudal máximo utilizando la metodología del SCS, es el siguiente:
Paso 2 Se calculan las lluvias de duración 6 horas y periodos de retorno de acuerdo a las avenidas del
proyecto.
Paso 3 Con base en el Número de Curva de la cuenca, se calcula la escorrentía para cada una de las
lluvias determinadas en el paso anterior con la ecuación:
( − . )
=
+ .
Paso 4 En función de la magnitud del tiempo de concentración se determina el valor del gasto unitario
(q), interpolando linealmente si es necesario con los ábacos que se tiene del SCS.
Paso 5 Por último, se multiplican el gasto unitario (q), la escorrentía (Q) y el área de la cuenca (A), para
obtener el gasto máximo (Max) en m3/seg.
Figura N31. Hidrograma tabulado SCS
El método de Muskingum fue presente por McCarthy (1938) y maneja relaciones caudal almacenamiento
variables. Este método modela el almacenamiento en un cauce mediante la combinación de dos tipos de
almacenamientos.
Durante el avance de la avenida el caudal de entrada es mayor que el de la salida y se forma lo que se
denomina cuña positiva y durante la recesión el caudal de entrada es menor al caudal de salida,
formándose una cuña negativa.
= + +
Dónde:
∆
− + 2
=
∆
(1 − ) +
2
∆
+
= 2
∆
(1 − ) +
2
∆
(1 − ) −
= 2
∆
(1 − ) +
2
Donde
.
= 0.18 .
Figura N33 Gráficos para aproximar el valor de X
Es una de las técnicas más recomendadas, está clasificada como un método hidrológico pero
proporciona valores comparables a los métodos hidráulicos.
Ecuaciones:
= + +
∆
−2
=
∆
2(1 − ) +
∆
+2
=
∆
2(1 − ) +
∆
2(1 − ) −
=
∆
2(1 − ) +
Dónde:
Dentro los sistemas automáticos de información hidrológica, están los SIG, los sistemas de información
geográfica con el cual podemos extraer toda la información física de la cuenca en estudio, como ser:
delimitación de la cuenca en estudio, área de la cuenca, longitud de los cursos de agua, pendientes,
tiempo de concentración.
7. APLICACIONES INFORMATICAS
Con toda la información básica obtenida mediante los Sistemas de Información Geográfica,
posteriormente podemos conocer los caudales de máximas crecidas haciendo uso del programa
computacional HEC HMS.
Para realizar estos cálculos parte de los datos de un hietograma de precipitación y a través de una serie
de modelos se calcula el exceso de precipitación (o precipitación neta), restando las pérdidas por
infiltración y detención, sobre la base de una función de velocidad de infiltración del agua en los suelos.
Los excesos de precipitación resultantes son luego seguidos mediante hidrogramas unitarios o técnicas
de la onda cinemática a la salida de la subcuenca generando un hidrograma.
La aplicación tiene una interfaz gráfica, que permite al usuario introducir la información necesaria para
una simulación, manejar los componentes de análisis hidrológicos, a través de módulos integrados y
obtener respuestas gráficas o tabuladas de fácil comprensión e impresión.
Establecer dimensiones preliminares de alcantarillas viales, luz en los puentes y sistemas de drenaje
pluvial.
Estudiar el impacto ambiental de las crecientes, causado por la transformación del uso de los suelos
en cuencas rurales y urbanas.