Primer Parcial Hidrología
Abigail Gual
Balance Hidrológico Criterio de Horton
𝑑𝑆
Cambio en el almacenamiento 𝑑𝑡 = 𝐼 (𝑡) − 𝑄(𝑡)
𝑑𝑆
Para un sistema cerrado 𝑑𝑡 = 0
Balance Global
𝑃 + 𝑅 + 𝐺 − 𝐸 − 𝑇 = ∆𝑆
𝑃 + 𝑅 + 𝐺 − 𝐸𝑇 = ∆𝑆
𝑅 = 𝑅𝑒 − 𝑅𝑠
*Manejar siempre las mismas unidades.
Ecuación de continuidad aplicada a un volumen de
fluido de fase única. Tiempo discreto.
𝑗 = 𝑖𝑛𝑑𝑖𝑐𝑒 𝑑𝑒 𝑡𝑖𝑒𝑚𝑝𝑜 𝑑𝑖𝑠𝑐𝑟𝑒𝑡𝑜
𝑆𝑜 = 𝐴𝑙𝑚𝑎𝑐𝑒𝑛𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑖𝑛𝑖𝑐𝑖𝑎𝑙
∆𝑠1 = 𝑎𝑙𝑚𝑎𝑐𝑒𝑛𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑖𝑛𝑡𝑒𝑟𝑛𝑜
Cantidad de agua almacenada 𝑆𝑗
Caudal incremento (∆𝑄) = se saca el promedio de Q,
se multiplica por el tiempo y se divide entre Ac.
Almacenamiento incremento (∆𝑆) = 𝑃𝑟𝑒𝑐𝑖𝑝. −∆𝑄
Almacenamiento Acumulado (𝑆𝑗) = Se inicia con el
que indica el ej.
Cuenca hidrológica Elevación de la cuenca
Leyes de Horton. Relación de bifurcación (Rb)
𝑎 = á𝑟𝑒𝑎 𝑒𝑛𝑡𝑟𝑒 2 𝑐𝑢𝑟𝑣𝑎𝑠
Ni = Número de cauces de orden i 𝐴 = á𝑟𝑒𝑎 𝑡𝑜𝑡𝑎𝑙
𝐸 = 𝑒𝑙𝑒𝑣𝑎𝑐𝑖ó𝑛 𝑚𝑒𝑑𝑖𝑎
Índice de compacidad 𝑒 = 𝑒𝑙𝑒𝑣. 𝑚𝑒𝑑𝑖𝑎 𝑒𝑛𝑡𝑟𝑒 2 𝑐𝑢𝑟𝑣𝑎𝑠 𝑑𝑒 𝑛𝑖𝑣𝑒𝑙
Densidad de drenaje
P = perímetro de la cuenca
S = área del círculo
Si Ic = 1, la cuenca es perfectamente circular, si se 𝐿 = 𝑙𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑙𝑎𝑠 𝑐𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒𝑠
𝐴 = á𝑟𝑒𝑎 𝑑𝑒 𝑑𝑟𝑒𝑛𝑎𝑗𝑒
aleja de 1, entonces es una cuenca alargada.
Criterio de Alvord
Sc = Inclinación que
presenta la cuenca.
D = Desnivel entre curvas Corriente de masividad
de nivel (km).
A = Área de la cuenca
(km2).
L = Longitud total de las curvas de niveles (km).
ℎ = 𝑒𝑙𝑒𝑣. 𝑚𝑒𝑑𝑖𝑎 𝑑𝑒 𝑙𝑎 𝑐.
Sumatoria de todas las 𝐴 = á𝑟𝑒𝑎 𝑑𝑒 𝑑𝑟𝑒𝑛𝑎𝑗𝑒.
longitudes.
Primer Parcial Hidrología
Abigail Gual
Precipitación Dobles acumulaciones
Intensidad Se basa en el graficado de dos series anuales de la
precipitación registrada en dos estaciones.
Curva masa: Son las precipitaciones acum. Y de ahí
se sacan las intensidades.
Hietograma: P(mm) contra t, horas.
Precipitación media en una cuenca
Media aritmética
Se grafican los datos de la columna (6) en abscisas contra
los datos de la columna (8) en ordenadas y se verifica la
Poligonos de Thiessen necesidad o no de efectuar una corrección. En caso
afirmativo, deben corregirse los valores erróneos de la
Pn: precipitación media de
columna (7) y presentarse en la columna (9).
cada poligono.
Aforo
An: área de cada poligono
AT: área total de la cuenca.
Isoyetas
Pn: precipitación media entre 2 Se saca la semisuma de las velocidades, luego la
isoyetas (semisuma). semisuma de las profundidades, la diferencia de
An: área comprendida entre 2 distancia, se multiplican las columnas y se realiza
isoyetas. la suma de todas las filas.
Datos faltantes Curva de calibración
Tomar 3 estaciones cercanas que posean datos
confiables.
Si difiere en menos de 10% - media aritmetica,
Más del 10% :
Se grafica altura contra caudal.
Tiempo de concentración
U.S National Weather Service
Primer Parcial Hidrología
Abigail Gual
Vd es la sumatoria de los caudales directos
𝑄𝑇 = 𝑄𝑏 + 𝑄𝑑
Cálculo de la lluvia efectiva. Se cuentan con los
datos de la tormenta y los caudales generados por
la lluvia.
Índice Phi. Es la tasa constante de abstracciones o
pérdidas en mm/hr.
Separación del escurrimiento base y esc. Directo.
Método de la línea recta
Estimación del índice Phi.
Método de flujo base fijo
Si la diferencia es negativa quiere decir que no hay
escurrimiento, por lo tanto, lo tomamos como cero.
Coeficiente de escurrimiento
Es la relación entre la tasa pico de escorrentía
directa y la intensidad promedio de precipitación
de una tormenta. Puede usarse para la
precipitación y escurrimiento de una tormenta.
Método de pendiente variable
Hidrograma Unitario
Pasos para tipo 1
1. Se hallan Qb y [Link] Qb es una recta solo se
toman en cuenta los valores entre lo dos puntos.
2. Sumatoria de Qd (Vd), y se calcula He (Vd/Ac)
3. HU cada ordenada de Qd entre He.
4. Graficar Qt, Qb y Qd.
Pasos para tipo 2.
1. Se tiene t, Q, y HU de una duración conocida (ej.
3 hrs), para cambiar la duración de HU, se crean
columnas de HU, por ej. Si quiero una duración de
*Los puntos de inflexión nos indican cuantos tipos de 6 hrs, se crean 6 columnas de HU. Posteriormente
flujos existen. Si tengo 2, indica que tengo 3 tipos de vamos copiando la columna de HU(3hrs) y vamos a
flujo. Se grafica Tiempo vs lnQ. desfasar 3 espacios en cada columna de manera
Precipitación en exceso o precipitación efectiva
sucesiva. (El número de espacios depende de la
duración de nuestra tormenta, si es de 4h, se dejan
4 espacios).
2. Se crea una columna de curva S correspondiente
a mi duración inicial (3h).
Primer Parcial Hidrología
Abigail Gual
Y esa columna será la suma de la fila, que va desde Estimación de K para el modelo de Horton
mi HU inicial hasta mi HU de 6 hrs (7 columnas). Este Se tiene la fc y fo, el volumen de infiltración f(t) y la
paso de hace para cada fila. duración.
3. Para obtener la curva S nueva (en este caso de 6h), Se integra la ecuación sustituyendo los valores.
se toma la columna de la curva S de 3h, se copia y se
desfasa 6 espacios.
4. Se crea una columna de diferencias de curvas S y
se hace la resta de Curva de 3 h menos curva de 6 h.
(Se excluyen los negativos)
5. Se determina la relación N
𝐿𝑎𝑠 ℎ𝑜𝑟𝑎𝑠 𝑞𝑢𝑒 𝑑𝑒𝑠𝑒𝑜 6
𝑁= = =2 Se tiene las tasas de precipitación I (mm/h) y la de
𝐿𝑎𝑠 ℎ𝑜𝑟𝑎𝑠 𝑞𝑢𝑒 𝑡𝑒𝑛𝑔𝑜 3
Y se crea una columna de relación N. infiltración f (mm/h).
6. Se determina la nueva HU, que es la dif. De curvas La precipitación total de la tormenta es la suma de
entre la relación N. las intensidades por delta T.
7. Se grafica T contra HU (toda la tabla).
Pasos para tipo 3.
1. Se tiene el tiempo, la HU de 3hrs, la ec. De Qb y
La infiltración acumulada es la suma de las
nos dan el escurrimiento total por una lluvia en
filtraciones por delta T.
exceso de 36 mm (he).
2. Para obtener nuestra Qd, será HU*He.
Y la altura de la precipitación efectiva o altura de
3. Se saca la Qb con la ec. escurrimiento (He) es la diferencia de la
4. Se obtiene QT que es la suma de Qb y Qd. precipitación total menos la infiltración acumulada.
Índice de infiltración Phi.
Infiltración
Modelo de Horton
Fc: capacidad de
infiltración después del Para estimar el índice Phi
tiempo prolongado de Se tiene el tiempo, las precipitaciones (I mm/h), y
humedad (mm/h) o el porcentaje de la precipitación efectiva con
volumen de infiltración respecto a la lluvia total, se suma todas las
acumulado. intensidades (Precipitación total)
Fo: Capacidad de infiltración inicial.
K: constante que representa la tasa de descenso de
f.
Profundidad acumulada de infiltración para un Se saca la precipitación efectiva (He) que es la
punto dado. sumatoria de las intensidades por el porcentaje
proporcionado.
Se crea la columna de I-Phi*deltaT, dándole un
valor de Phi. Se excluyen los negativos asignándoles
el valor de 0.
La suma de esa columna debe ser igual a He para
que se acepte Phi.
Primer Parcial Hidrología
Abigail Gual
Infiltrómetro
Se quiere calcular el agua infiltrada e infiltración, se
tiene las horas, el enrase y las lecturas.
Se crea una columna delta T, comenzando por cero,
es el intervalo de tiempo.
Se crea una columna de lámina parcial que será la
diferencia del enrase menos la lectura, el siguiente Las es se hallan de la tabla 3.1 calor latente.
valor será la lectura menos la siguiente lectura y así Para calcular ea (buscamos con la temperatura el
sucesivamente. medio).
Se crea la columna de tiempo acumulado. E con la es del agua.
Y otra de lámina acumulada. El Vol. Total perdido es la E*Ac (m3/d)
Por último, se obtiene la capacidad de infiltración La evaporación. Meyer.
que será lamina acum/tiempo acum expresado en
mm/h. Se dejó la primera casilla en blanco
Para graficar la lamina infiltrada es tiempo acum
contra lamina acum y para la gráfica de capacidad
de infiltración es tiempo acum contra capacidad de
infiltración.
Ajuste de curva
Para estimar la k del modelo de Horton con una f(t)
de prueba.
Fo y Fc se sacan de la tabla proporcionada, el primer
valor se determina observando el comportamiento
de la grafica de manera que intercepte con el eje Y.
Se crea una columna de f(t) Horton y se calculan los
datos utilizando la ecuación y una k de prueba.
Se saca una columna de RMS que será la diferencia
de la f(t) prueba menos la f(t) Horton y se eleva al
W2 la altura dada en el ejercicio
cuadrado. W1 a 2 m.
Se saca la sumatoria de esa ultima columna y se Medición de la evaporación con estanques
divide entre el número de datos.
Se grafica la f(t) de prueba y la f(t) Horton, y se
variará la k hasta que las dos graficas coincidan. De
igual manera la curva ajustada es aquella cuyo valor
de RMS es el mínimo.
Evaporación y transpiración
Humedad relativa (r)
Evapotranspiración (Doorenbos and Pruit)
La evaporación. Métodos de transferencia de
masas. Herbeck.
Solución
Obtener Kc dependiendo del cultivo (promedio).
Depende de la humedad (r), velocidad (W) e
insolación (I) se saca a y b (donde coincidan).
Primer Parcial Hidrología
Abigail Gual
De acuerdo a la latitud (checar si es N o S) tomar el
mes y latitud para hallar p.
Se calcula f (°C) y posteriormente se calcula la ETP
(mm/d)
Thornwaite and Wilm
Métodos probabilísticos y estadísticos de eventos
extremos.
Funciones de densidad de probabilidad
Se dan los años y P (mm), se crea una columna de n
(número de datos), se copia la columna de P y se
ordena de menor a mayor, por ultimo se crea una
columna de n/N*100, donde n corresponde a la fila
y N el numero total de datos. Se grafica P(mm)
ordenados contra la última columna. El gráfico nos
muestra la probabilidad.
Solución Periodo de retorno
Se calcula T^1.514 con las temperaturas dadas y se
suman, excluyendo los negativos (se eliminan).
Posteriormente se calcula J.
Se calcula C que es la parte del exponente de la
fórmula. Periodo de ocurrencia
Para calcular los km se halla los valores de la tabla
3.5 dependiendo de la latitud, si los valores no son
directos se interpola.
Se calcula la ETP (mm/mes).
Evapotranspiración Potencial (Penman)
Probabilidad de no excedencia
Primer Parcial Hidrología
Abigail Gual
Método de la variable reducida Para Gumbel
Paso 1. Encontrar el valor medio y la desviación Se repiten los mismo pasos a excepción de
estándar de las lluvias máximas anuales. En Excel
PROMEDIO(se selecciona toda la tabla) y
DESVEST.M().
Paso 2. Utilizando el valor medio y la desv. Est.
Calculadas se calcula los momentos de la fdp. El valor de Pi en excel es PI().
Se despejan a y c.
Paso 3. Calcular la variable reducida (YT) para los
periodos de retorno (Tr) de T años. Funcionamiento de los vasos
Paso 4. Se calcula la precipitación máxima anual
para los periodos de retorno (Tr) de T años.
Para encontrar el volumen útil y NAMINO.
Métodos de Momentos Se saca la columna de Xi-Di, luego la columna de (Xi-
Se dejan los valores máximos por cada año. Di)acum.
Se hallan las intensidades de cada uno, las P(mm) Se saca el valor de S de la columan de los
entre la duración y se reporta en mm/hr. acumulados. (Se identifica primero el máximo
Se crea una columna de orden que va desde 1 hasta positivo y luego el máximo negativo)
N. S = P1 -P2
Se saca el periodo de retorno Se calcula la columna de Vol. Que es S+(Xi-Di) para
el primer valor. Posteriormente es el recien calulado
menos el siguiente valor de (Xi-Di). Siguiendo una
forma escalonada. Los valores que superan la S, se
quedan como el valor de S. Ej. En excel
=SI((G18+E19)>1020,1020,G18+E19)
*Si la suma es mayor que 1020, arroja el valor de
Se copia toda la tabla de intensidades y se 1020, sino, reporta la suma.
selecciona ordenar datos de mayor a menor por Se calcula la columna de derrame, es el mismo
cada columna. algoritmo con la diferencia de que se reporten los
Debajo de esa misma tabla se saca el promedio y la valores excedentes.
varianza de cada columna. PROMEDIO() y VAR.S() =SI((G18+E19)<1020,"",((G18+E19)-1020))
Se saca a y b, a es la raíz de la varianza *Si la suma es menor que 1020, no des ningun valor,
sino, reporta la diferencia.
Por ultimo, crear una columna que reporte el estado
del vaso, en excel:
Se crea una tabla colocando esta formula para cada =SI((G19=0),"Vacio",SI((G19)<1020,"Con
columna de I. agua","Lleno"))
*Si el vol. Es igual a cero, reporta vacio, si es menor
a 1020, reporta con agua, si no (solo queda la
Se grafica la fila de duración (min) contra la fila de posibilidad de que sea mayor a 1020), reporta lleno.
intensidad (tomados de la última tabla) para una Tr
de x años.
Primer Parcial Hidrología
Abigail Gual
Método de Muskingum
Se crean las columnas de constantes, colocando en la primera fila los valores de las constantes
0.5∆𝑡 − 𝑘𝑥
𝐶0 =
𝑘(1 − 𝑥) + 0.5∆𝑡
𝑘𝑥 + 0.5∆𝑡
𝐶1 =
𝑘(1 − 𝑥) + 0.5∆𝑡
𝑘(1 − 𝑥) − 0.5∆𝑡
𝐶2 =
𝑘(1 − 𝑥) + 0.5∆𝑡
Comprobar las restricciones
𝐶0 + 𝐶1 + 𝐶2 = 1 (𝐶𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑑𝑎𝑑)
𝑘
≤ ∆𝑇 ≤ 𝑘 (𝐸𝑠𝑡𝑎𝑏𝑖𝑙𝑖𝑑𝑎𝑑)
3
Por ejem. Si dT=6
11
≤ 6 ≤ 11 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒
3
Crear la ulltima columna aplicando la fórmula de
𝑂𝑖+1 = 𝐶0 𝐼𝑖+1 + 𝐶1 𝐼𝑖 + 𝐶2 𝑂𝑖
Se grafica T contra Li y T contra Oi
Para hallar los valores de x y k, se tiene las Ii y las Oi, primero se crean las columnas de
(𝐼𝑖 + 𝐼𝑖+1 ) (𝑂 + 𝑂𝑖+1 )
I̅ = ̅= 𝑖
O
2 2
Posteriormente se crea la columna de Si+1
𝑆𝑖+1 = (I̅ − O̅ )∆𝑡 + 𝑆𝑖
Iniciando con cero en la fila de las celdas vacías.
Se crea la sig. Columna de
[𝑥𝐼𝑖+1 + (1 − 𝑥)𝑂𝑖+1 ]
Utilizando la primera casilla para asignarle un valor de x
Se grafica Si+1 contra la ultima columna. Variar el valor de x hasta que se pueda apreciar una
linea recta.
K será la pendiente de la recta, tomando dos puntos de ella y se reporta en hrs. (La pendiente se
multiplica por 3600).
Primer Parcial Hidrología
Abigail Gual
Primer Parcial Hidrología
Abigail Gual
Primer Parcial Hidrología
Abigail Gual
