UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.

1
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

TITULAR: Ing. Diego Taus

JTP: Ing. David Hardoy

FACULTAD TECNOLOGICA NACIONAL

REGIONAL AVELLANEDA

CATEDRA DE HIDROLOGIA
Y
OBRAS HIDRAULICAS

UNIDAD TEMATICA Nº2

TRATAMIENTO PROBABILÍSTICO DE LA INFORMACIÓN HIDROLÓGICA
CURSO 2020

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
2
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

TRATAMIENTO PROBABILÍSTICO DE LA INFORMACIÓN HIDROLÓGICA

Los procesos hidrológicos (tormentas, lluvias, crecidas, sequías, etc.) evolucionan en el

espacio y en el tiempo en forma parcialmente predecible y en forma parcialmente
aleatoria, a estos procesos se los conoce como estocásticos.
En algunos casos la componente aleatoria del proceso es tan grande, comparada con la
variables deterministicas que se justifica que se trate el tema como puramente aleatorio.

ALGUNOS CONCEPTOS DE PROBABILIDAD

Variable aleatoria, sus valores se rigen por leyes de azar y debe utilizarse el concepto de
probabilidad.
La variable aleatoria puede ser:

• Discretas: toman valores solo en puntos determinados

Xt Xt

0 2 4 6 N t 0 24 6 N t

Variable aleatoria discreta

• Continuas: toman valores en cualquier punto del eje o entre dos puntos extremos

Xt

0 24 6 N t

Variable aleatoria continua

Funciones de probabilidad: son las funciones que definen la probabilidad asociada a la

variable aleatoria
Se utilizan de dos tipos:

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
3
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

• Funciones de densidad
• Funciones de distribución

La función de densidad f(x) se utiliza para obtener las probabilidades por integración, por
ejemplo, la probabilidad de que X se halle entre x1 y x2 se obtiene según:

Funcion de densidad
x2

P( x1
X
x2)= ∫ f(x) dx

x1

La función de distribución F(x) representa la probabilidad acumulada desde -∞a X1 es decir:

X1

F(x) = ∫ f(x) dx
-∞

f(x)

F(x)

X1 X

F(x)
1.0

F(x)

Una variable aleatoria puede ser la precipitación, descripta por una distribución de
probabilidad, que determina la posibilidad de que una observación p de la variable caiga en
un rango especificado de P.

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
4
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Variables aleatorias conjuntas

Puede suceder que existan variables aleatorias conjuntas, en este caso se debe
considerar dos o más variables aleatorias simultaneamente. Se utilizan las funciones
conjuntas de distribución y densidad.
En el caso de dos variables:
x1 y1

F(x) = ∫ ∫ f(x,y) dx dy
-∞ -∞

Variables aleatorias independientes:

son aquellas en que:

f(x,y)= f(x).f(y) ∧ F(x,y)=F(x). F(y)

Muestra : en un conjunto de observaciones X1,X2 , X3 ........Xn de la variable aleatoria ,

que son sacadas de una población infinita , que posee propiedades estadísticas constantes.

La muestra constituye una estimación del fenómeno real y son los datos disponibles sobre la
población a partir de los cuales habrán de inferirse las propiedades de esta.
La población es el proceso real, que nunca es del todo conocido.

Espacio muestral: es el conjunto de todas las muestras que puedan extraerse de una
población.

Evento: es un subconjunto del espacio muestral.

Por ejemplo el espacio muestral para los caudales medios anuales de un río es
teóricamente el rango de cero a infinito positivo (con limite superior obviamente mas
cercanos que esto) y un evento x puede ser la ocurrencia de un caudal medio anual
superior que un Q de 300 m3/s.
La probabilidad de un evento, es la frecuencia relativa del mismo cuando el tamaño de la
muestra aumenta.

NA
P(A)= lim
n⇒∞ n

NA: numero de observaciones en el rango del evento A

n: numero de observaciones de la muestra

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
5
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Principos de las probabilidades

1- Probabilidad total
Si U es el conjunto universal , y U = A1 +A 2+A3+.......+An
donde Ai = subconjuntos de U
P(A1) +P(A2) + ........+P(An) = P(U)=1

2-Complementariedad
Si A´ es el complemento de A ⇒ P(A´) = P(U) -P(A) = 1-P(A)

3 -Probabilidad Condicional
P (B/A)=P(A ∩ B) / P(A) A y B son independientes P(B∩A) = P(B).P(A)

U
A B

A∩ B

Una de las funciones de densidad de probabilidad más conocidas es la distribución normal

(-½(x-µ)2/σ2)
1 exp
f(x) = rango de la variable -∞<x<+∞
√ 2π σ

si utilizamos una variable reducida z , tendremos

1 exp (-z2/2)
z = (x- µ) / σ ⇒ f(z) = rango de la variable -∞
z
+∞
√ 2π

ésta es la conocida distribución normal estándar, la función de distribución de probabilidad

normal estándar será :

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
6
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

1 exp (-u2/2)
F(z) = .du donde u es una variable de integración
√ 2π

Los valores de F(z) se pueden aproximar mediante los polinomio de Abramowitz y Stegun
F(z) = B para z<0
F(z) = 1- B para z≥0

B= ½ (1+0.196854 . z  + 0.115194 . z 2 + 0.000344 . z 3 + 0.019527. z 4 )-4

PARÁMETROS ESTADÍSTICOS

El objetivo de la estadística es extraer la información esencial de un conjunto de datos, a

través de las muestras de la variable aleatoria, que son los datos disponibles de una
población infinita, la cual posee propiedades estadísticas constantes, se trata de inferir y
estimar las propiedades del fenómeno real (población) que nunca es del todo conocido,
existiendo un cierto grado de incertidumbre.

Parámetro de la población Estadística de la muestra

- Media aritmética.
x2 ∑ Xi
µ= E(x) = ∫ x f(x) dx X=
x1 n

Es el promedio de los valores de la muestra

- Mediana.

F(x) F(x)
1.0 1.0

0.5 0.5
F(x)=0.5
X
x x

x tal que F(x)=0.5 Valor de la informacion al 50% del percentil

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
7
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

- Media geométrica.

(1/n)
Anti log E(logx) ∏ xi

- Varianza.

∑ (Xi - X)2
σ2 = E (x- µ)2 S2 =
(n - 1)

- Desviación estandar

∑ (Xi - X)2
σ= E (x- µ)2 S=
(n - 1)

- Coeficiente de variación

σ s
CV = CV=
µ X

- Coeficiente de asimetría

n ∑ (Xi - X)3
E (x- µ)3 Cs =
γ= (n - 1)(n-2)S3
σ3

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
8
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

DISTRIBUCIÓN DE PROBABILIDAD PARA LAS VARIABLES HIDROLÓGICAS

Distribución Normal Rango (-∞ , +∞ )

La precipitación anual, calculada como la suma de los efectos de muchos eventos

independientes, tienden a seguir una distribución normal.
Principal limitación: variables hidrológicas no son negativas, información hidrológica tiende a
ser asimétrica.

(-½(x-µ)2/σ2)
1 exp
f(x) = rango de la variable -∞<x<+∞
√ 2π σ donde µ = X σ = s

Distribución Lognormal

Describe la distribución de la conductividad hidráulica en un medio poroso, la distribución de

tamaño de gotas de lluvia en una tormenta.
Principal limitación: tiene solamente dos parámetros, requiere que los logaritmos de los
datos sean simétricos alrededor de su media

(-½(y-µy)2/σy2)
1 exp
f(x) = rango de la variable 0<x<+∞
x σ √ 2π donde y =log x µy = y σy = sy

Distribución Poisson

Modelo de probabilidad para embalse lineal donde λ =1/k, donde k es la constante de

almacenamiento del embalse lineal.

1
f(x)= λ e λ=
-λx
donde x0
X

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
9
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Distribución Gamma

Se la ha aplicado a la distribución de la precipitación (milímetros de lluvia caído) en

tormentas.

λβ ⋅ X (β −1) ⋅ e −λx x X2 1
f (x) = donde x ≥ 0 λ = β= =
Γ(β) Sx 2 Sx 2 CV

Γ( )= funcion Gamma

Cabe recordar que Γ(β)=(β -1)!, para β entero y positivo, sino se calcula mediante.

Γ(β) = ∫ U(β-1) e-u du

0

Distribución Pearson Tipo III

Se la utiliza para describir la probabilidad de picos de crecientes máximos anuales.

λβ ⋅ ( X − ε)(β −1) ⋅ e −λ( x − ε ) Sx 4

f (x) = donde x ≥ ε λ = β=
Γ(β) β Cs 2

ε = X − Sx β
Γ( )= funcion Gamma

Distribución Log-Pearson Tipo III

Se la utiliza en [Link]. como distribución estándar para análisis de frecuencia de crecientes

máximas anuales.

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
10
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

λβ ⋅ ( y − ε)(β −1) ⋅ e −λ( y − ε ) Sy 4

f (x) = donde log x ≥ ε λ = β=
x ⋅ Γ(β) β Csy 2

ε = y − Sy β
Γ( )= funcion Gamma

Distribución Valor Extremo

Se demuestra que las distribuciones de valores extremos seleccionados de conjuntos de
muestras de cualquier distribución de probabilidad convergen en una de las tres formas de
distribución de valor extremo. Esta distribución tiene 3 tipos, la Tipo1, es muy conocida
como distribución de Gumbel.

Valor extremo Tipo I

1 (x-u) (x-u)
f(x) = exp - exp
α α α

donde el rango de la variable -∞<x<+∞

6
α = Sx
π

U= X –0.5772α

REGRESIÓN LINEAL SIMPLE

El análisis por regresión lineal simple permite relacionar la linealidad entre dos variables de
un suceso. Este análisis es utilizado en ocasiones para evaluar la linealidad que existe entre
la lluvia y la escorrentía que la misma genera sobre una cuenca. El procedimiento es el
siguiente: Sea la recta
Se trata de conocer los parámetros α y β que me indique la existencia de linealidad entre las
variables x e y.

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
11
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

(el cuadrado del error tiene que ver con que este puede ser en mas o en menos al valor real)

Donde:

Reordenando:

Despejando:

Para comprobar si existe una linealidad entre variables se calcula el coeficiente de

correlación mediante la expresión:

Cuanto mas cercano a 1 ó -1 mas eficiente es la linealidad entre variables

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
12
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

ANÁLISIS DE FRECUENCIA

Los sistemas hidrológicos son afectados por eventos extremos, tormentas muy fuertes,
crecientes, sequías, precipitaciones máximas, etc. La magnitud de un evento extremo esta
inversamente relacionado con su frecuencia de ocurrencia, es decir eventos muy severos
ocurren con menor frecuencia que eventos mas moderados.

El periodo de retorno T de un evento X≥ xτ puede definirse como el valor esperado de τ

E(τ) es el valor promedio medido sobre un numero de ocurrencias suficientemente grande.
Entonces el periodo de retorno de un evento de magnitud dada es el intervalo de recurrencia
promedio entre eventos que igualan o exceden una magnitud especificada.

Si p es la probabilidad de que X≥xτ , para cada observación existen dos posibilidades:

1- de éxito es decir X≥xτ de probabilidad p

2- de falla es decir X<xτ de probabilidad 1-p

Dado que las observaciones son independientes, la probabilidad de un intervalo de

recurrencia de duración τ es el producto de las probabilidades de τ-1 fallas seguidas por un
exito∴(1-p) τ-1.p

E (ττ) = T=Στ.(1-p) τ-1.p= 1/p

τ=1

Conclusión
La probabilidad de ocurrencia de un evento en cualquier observación es el inverso de su
periodo de retorno.

≥xττ) = 1/ T
P(X≥

Si bien se ha mencionado la reciprocidad entre el tiempo de recurrencia y la frecuenta de

ocurrencia en base a un concepto empírico existe un desarrollo analítico que permite
justificar la expresión:

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
13
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Dado que se trata de una

INTERVALO O PERIODO DE RECURRENCIA “T”

Para aclarar el concepto de T desde el punto de vista hidrológico aplicamos los conceptos
de probabilidad por ejemplo al cálculo de las crecidas máximas.

Se define el intervalo o periodo de recurrencia como el periodo promedio en años dentro

del cual una crecida determinada será igualada o superada solo una vez.

Si se ordenan en forma decreciente los caudales máximos de cada año de una serie de
caudales, siendo N el numero de años de observación, asignando a cada uno el numero de
orden m, el caudal mayor tendrá un numero de orden m=1 y el ultimo de la serie tendrá un
numero de orden m=N, observar que el ultimo termino de la serie es N solo si se ordena un
caudal máximo por año.

Asi definido el intervalo de recurrencia T en años de una crecida de orden m, será:T=N/m

Si N = 30, el periodo de recurrencia de la crecida de orden 3 será de 10años y la crecida

máxima tendrá un periodo de recurrencia de 30 años (m=1), es decir una probabilidad de
ocurrencia del orden del 3.33%.

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
14
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

RELACION ENTRE LA VIDA UTIL DE LA OBRA ; T DEL EVENTO DE DISEÑO Y RIESGO

DE OCURRECIA DE UN EVENTO DE T SUPERIOR

Entendido el concepto de periodo de retorno se calculara cual es la probabilidad de que un

evento con periodo de retorno de T años ocurra al menos una vez en N años.

1. Primero consideramos la situación que no ocurra el evento de T años en N años,

Esto requeriría una secuencia de N fallas sucesivas, recordando que p es la probabilidad de
que X≥≥xττ , tendremos que la probabilidad de que X<<xττ es 1-p ⇒ P(X<xτ cada año) = (1-p)
Como cada evento es independiente del otro, es decir que ocurra una precipitación, caudal
máximo, o cualquier otra variable hidrológica, no esta condicionada a las que ya han
sucedido, tendremos que la probabilidad para N fallas sucesivas:

P(X<xτ cada año durante N años) = (1-p)(1-p) ........(1-p)= (1-p)N

El complemento de esta situación es lo requerido, es decir, si el conjunto Universal esta

formado por la unión de los siguientes subconjuntos

P(X≥xτ al menos una vez en N años) =1- P(X<xτ cada año durante N años)
P(X≥xτ al menos una vez en N años) = 1-(1-p)N

pero p=1/T ⇒ ≥xττ al menos una vez en N años) = 1-(1-1/ T)N

P(X≥

Si llamamos N a la vida útil esperada de la obra, Tr el periodo de retorno del evento

adoptado para el diseño del proyecto, P es la probabilidad de que un evento de recurrencia
mayor ocurra al menos una vez durante la VU, por lo tanto P es el riesgo aceptado de falla.

Aplicando la formula tendremos los siguientes datos.

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
15
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

PERIODO DE RETORNO NECESARIO PARA UNA PROBABILIDAD DE OCURRENCIA

ACEPTADA PÁRA LA VIDA ÚTIL DE LA OBRA

Por ejemplo la probabilidad P que un evento ocurra en un periodo de recurrencia de 1000

años, en una obra de 50 años de vida útil será:

P= 1-(1-1/T)N= 1-(1-1/1000)50 ≈ 5%

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
16
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

TIEMPO DE RECURRENCIA EN FUNCIÓN DE LAS CARACTERÍSTICAS DE LA OBRA

El tiempo de recurrencia esta relacionado con las dimensiones y la importancia de la obra a
ejecutar. En la siguiente planilla se dan estimaciones de los Tr para los distintos tipos de
estructuras

Es muy común el desarrollo de un análisis económico del impacto de la obra sobre los
daños que esta sea capaz de originar en función de su capacidad de servicio.
Por lo tanto es de utilidad la determinación del un tiempo de recurrencia óptimo que permita
contemplar los costos de la obra frente a los costos de los daños que se puedan generar si
la obra se ve superada en su capacidad. Para ello se realiza un análisis hidroeconómico.

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
17
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

ANALISIS HIDROECONOMICO

Para un proyecto determinado, el tiempo de recurrencia (Tr) optimo del evento hidrológico,
para el diseño de la obra, puede surgir de un análisis económico.
Como es de esperar, a medida que el Tr de diseño aumenta, los costos de la obra se
incrementan, pero los daños esperados que pueden suceder disminuyen, debido a una
mayor previsión (mayor diámetro de las cañerías de desagües, mayor capacidad de
conducción de canales, etc.).
Por lo tanto sumando al costo de la obra, los costos de los daños esperados anualmente,
el Tr óptimo de diseño surgirá como el valor que hace mínimo al costo total.

Costo de los daños

0 2 5 10 25 50 100 200 (Tr del evento hidrologico)

0 0.5 0.2 0.1 0.04 0.02 0.01 0.005 (Prob. anual de excedencia del e.h. (1/Tr)

Este grafico representa el costo del daño que resultaría si ocurriese un evento hidrológico,
con un Tr especificado
Si adoptamos para diseñar una obra un tiempo de recurrencia Tr, para nuestro evento
hidrológico (e.h.), del que surge una magnitud Xtr, para dicho evento.

Cualquier e.h. de magnitud inferior no provocaría daños, dado que la estructura de la obra
los evitaría, pero no podría evitar daños para e.h. tal que X >Xt

si X ≤ Xtr , la obra evita los daños

⇒
si X >Xt , la obra no puede evitar los daños

Por consiguiente el costo del daño esperado (CDE), se encuentra multiplicando la

probabilidad de que un e.h. de magnitud X ocurra en un año dado , por el daño que
resultaria de este evento D(x) , integrando para X >Xt

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
18
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Costo Anual Esperado = ∑ D(x) f(x) dx

Xtr
Para cuantificar los daños causados por e.h., se debe determinar el área, la duración, la
profundidad y estimar los costos en perdidas materiales, costos sociales (de muy difícil
estimación), costos de restauración, costos de limpieza, etc. Para cada Tr del e.h.
analizado.

DETERMINACION DEL Tr OPTIMO

El costo total de la obra surge de la suma del:

• costo de inversión para la construcción
• costo del riesgo asumido , para eventos de Tr superior que puedan ocurrir

Costo de los daños

El area bajo la curva

El area bajo la curva
representa el costo
representa el costo anual
anual del daño evitado
del daño esperado o costo
del riesgo

0 2 5 10 25 50 100 200 (Tr del evento hidrologico)

Costos

CT

CO
CRD

0 2 5 10 25 50 100 200 (Tr del evento hidrologico)

CO: costo de la obra para un evento hidrologico de diseño de Tr

CDO: costo del riesgo de daño

CT =costo total =CO+CRD

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
19
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Para CT mínimo, el valor de Tr es el óptimo, el cual se deberia adoptar, para el e.h. de

diseño. Demás esta decir que no siempre se puede adoptar éste Tr óptimo, dado que
pueden existir otras restricciones a la hora de formular el proyecto, como puede ser, entre
otras, partidas de dinero asignadas con anterioridad que fijan el monto total de las obras.
Se ajunta un ejemplo práctico del análisis hidroeconómico en el anexo de la unidad

METODOS DE ESTIMACION DE CUANTILES

METODOS NO PARAMETRICOS
Las metodologías empleadas en la obtención de X-T (donde X podrá ser Qmaximos ,
Precipitaciones máximas , Q mínimos , etc.), admiten una primera clasificación, según se utilicen
modelos paramétricos o modelos no paramétricos.

Los métodos no paramétricos no modelan la serie de datos mediante formulas algebraicas

basadas en uno o mas parámetros que es preciso ajustar a partir de la muestra, sino que
utilizan directamente dicha muestra. El método no paramétrico mas simple es el ajuste
gráfico de una curva a los datos de la muestra representados en papel de probabilidad, con
la que se establece una relación X-T.

Las series de datos a utilizar (de caudales, precipitaciones máximas, etc. ) pueden ser de
dos tipos :

1-series de valores extremos :Aquí los datos son las crecidas máximas de cada año, estas
series tienen la ventaja de ajustarse mejor a la condición de independencia de crecida de un
año a otro.

2-series de duración parcial : aquí se toman los valores de caudal superiores a un cierto
caudal base.

GRÁFICAS DE PROBABILIDAD
Como una verificación de que la distribución de probabilidad se ajusta a un conjunto de
datos hidrológicos, estos se pueden graficar en un papel de probabilidad que linealice la
función distribución, luego los datos graficados se ajustan por medio de una recta con
propósitos de interpolación y extrapolación. El ajuste gráfico realizado de forma manual
permite introducir de forma cualitativa tanto información histórica como regional, pero
presenta el problema de la subjetividad del método. La aplicación de este método exige
asignar a cada uno de los datos una determinada probabilidad muestral que permita
representarlos en un determinado papel de probabilidad, siendo habitual el uso del
denominado “papel Gumbel”. Para la asignación de la recurrencia (inversa de la
probabilidad) se ordenan en forma decreciente los valores máximos de cada año de la serie,
siendo N el numero de años de observación, y asignando a cada uno el numero de orden m.
Realizado por el Ing. DIEGO TAUS 2020 UT Nº2
REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
20
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Se podrá asignar, el intervalo de recurrencia T de un dato de orden m, mediante la

aplicación de formulas del tipo:

1 α
N+1 - 2α
T (xj) = =
p m-α

Siendo:
Xj= valor que ocupa el lugar j de la serie ordenada de mayor a menor
N = numero total de datos
α = parámetro de la formula

Distintos valores de α conducen a diferentes expresiones, siendo las más empleadas las
siguientes:

• T= (N+1)/m α=0 Expresión de Weibull

Aquí se indica un periodo de retorno un año mayor que el periodo de retorno del registro del
valor máximo

• T= N/(m-0.5) α = 0.5 Formula de Hazen (1930)

• T= (N+0.4)/ (m-0.3) α= 0.3 Formula propuesta por Chegoday

Usada en Rusia y en los países de Europa Oriental

• T= (N+0.12)/ (m-0.44) α= 0.44 Formula propuesta por Gringorten

Otras formulas propuestas que no involucran el parámetro α son:

• T=N/m Conocida como la Formula de California

Aquí para N=m se le asigna una probabilidad del 100% , que puede ser difícil de graficar en
una escala de probabilidad , como un ajuste de ésta se utiliza generalmente alguna de las
formulas anteriores

• T=(3N+1)/(3m-1) Conocida como la Formula de Tukey

Utilizada para asignar tiempos de recurrencias en precipitaciones máximas
Realizado por el Ing. DIEGO TAUS 2020 UT Nº2
REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
21
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

• T=(N -¼)/(m-3/8) Conocida como la Formula de Blom

Diversos trabajos realizados por Cunnane (1978), indican que la formula de Weibull,
conduce a resultados sesgados, siendo preferible la de Gringorten si la población
subyacente es Gumbel (o exponencial) o la Hazen para poblaciones con un sesgo mayor.
Para datos normalmente distribuidos, encontró que la formula de Blom esta bastante
cercana a no ser sesgada. Con independencia del empleo o no del método gráfico, es
siempre recomendable comparar gráficamente los resultados obtenidos por otros
procedimientos, con los datos debidamente representados.
Una vez que los datos se han identificados y ordenado y las posiciones de graficacion
calculados puede elaborarse un gráfico de magnitud x-T para ajustar gráficamente una
distribución, para ello se elabora un cuadro de la siguiente manera:

DISTRIBUCIÓN DE GUMBEL

El problema estadístico a resolver es el ajustar alguna ley de distribución de probabilidad a

la variable aleatoria X , caudal de crecida de cada año , definida a partir de una serie de N
muestras , que contiene cada una 365 valores extremos de una serie (caudales máximos
diarios). Si aplicamos la distribución de Gumbel tendremos:
P(x)= 1-exp(-exp(-y)) donde P es la probabilidad que un caudal máximo anual X sea
igualado o excedido.
P= exp(-exp(-y)) donde P es la probabilidad que un caudal máximo anual X no sea
superado.
y es la variable reducida ⇒ y= (X-u)/α α
siendo u y α parámetros cuya estimación a partir de los datos de la muestra puede hacerse
por los métodos de máxima verosimilitud (muy complicado en este caso ) o por el método de
los momentos, aplicando esto ultimo resulta:

α
u = XMEDIA - 0.5772.α u :es la moda de la distribución, punto de mayor densidad de
probabilidad

α= 6/π σ =0.7796968 σ ≈0.78 σ

(X - XMEDIA + 0.45.σ) (1)
reemplazando se obtiene ⇒ y =
0.78.σ

X: es la magnitud del evento cuya probabilidad de ocurrencia es P

XMEDIA : es la media aritmética de los eventos de la serie
σ: es la desviación estándar de la muestra

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
22
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Otra forma de calcular “y” (variable auxiliar) es la siguiente:

si P(x)= 1-exp(-exp(-y))= 1/T

⇒ exp(-exp(-y))= (T-1)
⇒ y = -ln(ln ((T-1)/T) ( 2)

combinando la (1) y la (2) resulta

αy
X= XMEDIA +α

Nota: Hay que tener en cuenta que las expresiones de u = XMEDIA - 0.577α α y α= 0.78 σ
son validas para series de datos extensas, si no lo son es decir para muestras pequeñas, se
deben corregir resultando las siguientes expresiones:
u = XMEDIA - yn .σ/ σn
α= σ/σσn
Donde yn y σn fueron calculados por Gumbel para distintos valores de n (longitud de la
serie):

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
23
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

DISTRIBUCIÓN TEÓRICA DE CRECIENTES

Debido a la longitud de los registros es normalmente corta no es posible determinar la

distribución de frecuencias mas apropiada para ser usada para analizar las probabilidades
asociadas con crecientes. La distribución Log-Pearson Tipo III ha sido recomendada para
las agencias federales de los [Link]., también la primera distribución asintótica de los
valores extremos (Distribución de Valor Extremo Tipo I) comúnmente conocida como
distribución Gumbel también tiene amplia aceptación.

DISTRIBUCIÓN LOG-PEARSON III

El procedimiento recomendado para el uso de esta distribución consiste en convertir la serie

de datos a sus logaritmos y luego calcular los siguientes parámetros

Media Media logX = 1/n . ∑log X

Desviación estándar σlogx = ∑ (log X- Media logX)2

n-1

n ∑ (log X- Media logX)3

Coeficiente de asimetría G=
(n-1)(n-2)(σlogx)3

El valor de X para cualquier nivel de probabilidad se puede calcular a partir de la ecuación:

log X= Media logX +K σlogx

donde K es un valor tabulado (ver tabla Nº1)

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
24
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

OTRAS DISTRIBUCIONES DE INTERÉS HIDROLÓGICO SON LAS SIGUIENTES:

CONCLUSIÓN
De lo visto en la unidad se ponen en evidencia los distintos recursos utilizados para la
determinación del tiempo de recurrencia mas adecuado para cada obra:

Como dato obtenido del análisis estadístico de los registros disponibles.

Según el tipo de Obra.
En función de la magnitud de la Obra.
En función de la vida útil de la Obra.
Como dato obtenido del análisis de un balance hidroeconomico de la obra.
Etc.

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
25
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

RESUMEN DE LAS FUNCIONES ANALIZADAS EN ESTA UNIDAD

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
26
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
27
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
28
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

TABLA DE AJUSTE DE LA MEDIA Y DESVIACIÓN SEGÚN GUMBEL

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
29
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

TablaNº1 - VALORESDEKPARALADISTRIBUCIÓNLOG-PEARSONTIPOIII

Coeficiente PERIODODERETORNOENAÑOS
De 1.0101 1.25 2 5 10 25 50 100
Asimetría NIVELDEPROBABILIDADEN%
G 99 80 50 20 10 4 2 1
3.0 -0.667 -0.636 -0.396 0.420 1.180 2.278 3.152 4.051
2.8 -0.714 -0.666 -0.384 0.460 1.210 2.275 3.114 3.973
2.6 -0.769 -0.696 -0.368 0.499 1.238 2.267 3.071 3.889
2.4 -0.832 -0.725 -0.351 0.537 1.262 2.256 3.023 3.800
2.2 -0.905 -0.752 -0.330 0.574 1.284 2.240 2.970 3.705
2.0 -0.990 -0.777 -0.307 0.609 1.302 2.219 2.912 3.605
1.8 -1.087 -0.799 -0.282 0.643 1.318 2.193 2.848 3.499
1.6 -1.197 -0.817 -0.254 0.675 1.329 2.163 2.780 3.388
1.4 -1.318 -0.832 -0.225 0.705 1.337 2.128 2.706 3.271
1.2 -1.449 -0.844 -0.195 0.732 1.340 2.087 2.626 3.149
1.0 -1.588 -0.852 -0.164 0.758 1.340 2.043 2.542 3.022
0.8 -1.733 -0.856 -0.132 0.780 1.336 1.993 2.453 2.891
0.6 -1.880 -0.857 -0.099 0.800 1.328 1.939 2.359 2.755
0.4 -2.029 -0.855 -0.066 0.816 1.317 1.880 2.261 2.615
0.2 -2.178 -0.850 -0.033 0.830 1.301 1.818 2.159 2.472
0 -2.326 -0.842 0 0.842 1.282 1.751 2.054 2.326
-0.2 -2.472 -0.830 0.033 0.850 1.258 1.680 1.945 2.178
-0.4 -2.615 -0.816 0.066 0.855 1.231 1.606 1.834 2.029
-0.6 -2.755 -0.800 0.099 0.857 1.200 1.528 1.720 1.880
-0.8 -2.891 -0.780 0.132 0.856 1.166 1.448 1.606 1.733
-1.0 -3.022 -0.758 0.164 0.852 1.128 1.366 1.492 1.588
-1.2 -3.149 -0.732 0.195 0.844 1.086 1.282 1.379 1.449
-1.4 -3.271 -0.705 0.225 0.832 1.041 1.198 1.270 1.318
-1.6 -3.388 -0.675 0.254 0.817 0.994 1.116 1.166 1.197
-1.8 -3.499 -0.643 0.282 0.799 0.945 1.035 1.069 1.087
-2.0 -3.605 -0.609 0.307 0.777 0.895 0.959 0.980 0.990
-2.2 -3.705 -0.574 0.330 0.752 0.844 0.888 0.900 0.905
-2.4 -3.800 -0.537 0.351 0.725 0.795 0.823 0.830 0.832
-2.6 -3.889 -0.499 0.368 0.696 0.747 0.764 0.768 0.769
-2.8 -3.973 -0.460 0.384 0.666 0.702 0.712 0.714 0.714
-3.0 -4.051 -0.420 0.396 0.636 0.660 0.666 0.666 0.667

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
30
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

TABLA Nº2 - VALORES DE Kt PARA LA DISTRIBUCIÓN PEARSON TIPO III

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
31
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

EJEMPLO DE APLICACIÓN
Se analizará el siguiente registro de caudales aplicando los métodos vistos
Años Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Nº DATOS MAX ANUAL

1966 7.7 7.9 13.8 10.5 5.8 4.2 4.5 4.1 3.9 8.3 8.2 5.8 1 13.8

1967 8.9 19.5 17.2 11.1 6.7 8.7 10.8 9.2 7.1 11.7 7.4 5.5 2 19.5

1968 10.3 12.1 25.6 16 6 4.3 7.6 5.3 8.5 11.3 7.8 6.2 3 25.6

1969 10.1 17.4 12.3 17.9 10.4 6.3 3.8 10.1 10 9.2 11.2 14.4 4 17.9

1970 20.3 21.2 16.8 8.8 9.1 12.1 5.7 2.7 5.6 11.3 11.5 9.4 5 21.2

1971 11.2 15.8 21.3 13.7 6.6 8.3 7.7 6.2 6.4 11 8.1 10.2 6 21.3

1972 12.6 15.5 22.3 16.6 8.8 6.4 8.7 4.8 4.6 8.7 13.7 13 7 22.3

1973 11.5 20.5 20 16.8 9.8 8.9 8.5 7 6.6 9.5 8 10.6 8 20.5

1974 9.2 19.5 24.7 11.6 4.3 6.3 9.8 6.1 5.7 12.5 15.2 14.4 9 24.7

1975 10.8 20.8 16.3 12.7 9.2 10.1 9.4 8.9 9.2 10.9 9 6.8 10 20.8

1976 12.6 19.8 18.2 18 9.1 7.3 7.1 8.1 4.8 4.1 7.1 9.3 11 19.8

1977 10.7 12.5 13.1 15.6 7.6 7.1 4.8 5.5 8.7 11.7 7.7 7.6 12 15.6

1978 8.1 9.7 20.2 15.2 10.5 7.5 5.6 4.9 5.2 5.5 4.6 7 13 20.2

1979 7.4 14 24.6 19.5 8.4 2.8 1.8 5.5 6.4 5.1 3.5 6.1 14 24.6

1980 13.7 21.6 22 22.8 10.8 5.7 3.9 3.2 3.7 10.4 12.2 12.5 15 22.8

1981 11.8 13.6 21.5 15 6.3 3.4 3.7 3.5 2.1 10 9.9 13.5 16 21.5

1982 14.6 14.8 18.6 18.7 12.9 4.9 3.4 3.3 4.8 12.2 11.2 16.3 17 18.7

1983 19.6 19.4 23.2 17.6 11.1 5.1 4.2 2.8 4.5 10.9 8.3 15.2 18 23.2

1984 11.9 21.9 25.6 19.6 12 8 7.8 4.8 5.1 11 12.8 8.4 19 25.6

1985 9.2 15.3 17.2 13 13.5 7.2 5.1 6.4 6.3 11.2 10.3 13.3 20 17.2

1986 9.8 16.6 19.5 17.3 8.4 3.4 3.6 3.5 5.5 8.5 8.9 10.3 21 19.5

1987 9.3 13.2 15.4 15.1 10 4 6.7 6.2 5 9.1 11.7 9 22 15.4

1988 12.7 17.9 13.1 16.4 11 5.6 5.7 5.1 5.4 9.5 11.2 9.9 23 17.9

1989 18.3 20.2 25.6 16.8 9 6.2 4.3 3.9 5.6 13.4 7.2 5.3 24 25.6

1990 10.8 19.5 18.5 18.7 11.7 8.3 6.7 4.5 3.7 9.3 11.9 12.5 25 19.5

1991 13.7 16.3 34.7 18.4 9.9 7.6 4.9 5.2 4.8 7 7.4 10.1 26 34.7

1992 10.3 16.8 16 13.2 9.8 7.9 5.9 3.6 5.8 5.9 9.8 9.7 27 16.8

1993 11.3 20.2 34.5 21.2 7 4.6 5.6 3.9 5.3 10.2 9 21.1 28 34.5

1994 24.4 17.8 28.6 22.6 9.6 6.2 6.7 7.6 6.3 5.7 7.8 11.9 29 28.6

1995 8.3 13.5 15.2 12.7 10.6 4.8 4.2 3 3 7.4 16.1 17.6 30 17.6

1996 15.4 16.7 22.1 27.4 14.2 7.3 5.7 4 4.1 8.3 5.3 4.6 31 27.4

1997 9.8 14 14.7 10.8 6.9 3.8 4.7 4.4 4.8 7.2 11 12.9 32 14.7

1998 9.4 14 19.5 15.1 9.2 6.6 5.4 4.3 4 13 15.7 10.3 33 19.5

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
32
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

METODO GUMBEL
Del registro de datos se realiza el siguiente desarrollo:
Años Caudales
1966 13.8
1967 19.5 X media 21.47
1968 25.6 Desviación 4.98
1969 17.9
1970 21.2
1971 21.3
1972 22.3
1973 20.5
1974 24.7
1975 20.8
1976 19.8
n = 33
1977 15.6 σn = 1.12
1978 20.2
1979 24.6
yn = 0.54
1980 22.8
1981 21.5
u = X media - yn. σ/σn = 19.1
1982 18.7
1983 23.2 α = σ/σn = 4.4
1984 25.6
1985 17.2
1986 19.5
1987 15.4
1988 17.9
1989 25.6 Luego Y será igual a:
1990 19.5
1991 34.7
1992 16.8 Y = (X - Xmedia + yn . σ/σn) σn/σ
1993 34.5
1994 28.6
1995 17.6 Luego X = 4.4 Y + 19.1
1996 27.4
1997 14.7
1998 19.5

Conocida la recta se puede calcular el caudal para los distintos Tr

Tr (Tr-1)/Tr ln((Tr-1)/Tr) Y = -ln(-ln((Tr-1)/Tr)) X = 4.44Y + 19.07

2 0.5 -0.693 0.367 20.7

5 0.8 -0.223 1.500 25.7
10 0.9 -0.105 2.250 29.1
50 0.98 -0.020 3.902 36.4
100 0.99 -0.010 4.600 39.5
200 0.995 -0.005 5.296 42.6

Resumen:
Caudal máximo 21 m3/seg para Tr = 2 años
Caudal máximo 29 m3/seg para Tr = 10 años
Caudal máximo 36 m3/seg para Tr = 50 años
Caudal máximo 39 m3/seg para Tr = 100 años
Caudal máximo 43 m3/seg para Tr = 200 años

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
33
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

METODO PEARSON TIII

3
Nº Orden Años caudales Log X (X - Xmedia)
1 1966 13.8 1.140 -451.2
2 1967 19.5 1.290 -7.6
3 1968 25.6 1.408 70.5
4 1969 17.9 1.253 -45.5
5 1970 21.2 1.326 0.0
6 1971 21.3 1.328 0.0
7 1972 22.3 1.348 0.6
8 1973 20.5 1.312 -0.9
9 1974 24.7 1.393 33.7
10 1975 20.8 1.318 -0.3
11 1976 19.8 1.297 -4.7
12 1977 15.6 1.193 -202.2
13 1978 20.2 1.305 -2.0
14 1979 24.6 1.391 30.7
15 1980 22.8 1.358 2.4
16 1981 21.5 1.332 0.0
17 1982 18.7 1.272 -21.2
18 1983 23.2 1.365 5.2
19 1984 25.6 1.408 70.5
20 1985 17.2 1.236 -77.8
21 1986 19.5 1.290 -7.6
22 1987 15.4 1.188 -223.6
23 1988 17.9 1.253 -45.5
24 1989 25.6 1.408 70.5
25 1990 19.5 1.290 -7.6
26 1991 34.7 1.540 2315.8
27 1992 16.8 1.225 -101.8
28 1993 34.5 1.538 2212.4
29 1994 28.6 1.456 362.5
30 1995 17.6 1.246 -57.9
31 1996 27.4 1.438 208.6
32 1997 14.7 1.167 -310.2
33 1998 19.5 1.290 -7.6
Total = 43.603 3807.6

Cantidad elementos 33
Media logX = 1.321
X media = 21.470
σlogX = 0.096
Desviación = 4.979

Coef. de asimetría G = n . Σ ( X - Xmedia )3 = 1.03

(n-1) . (n-2) . ( σ )
3

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
34
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Tr K K σlogX Media logX + K σlogX X

2 -0.168 -0.016 1.305 20.2

5 0.750 0.072 1.393 24.7
10 1.340 0.128 1.450 28.2
50 1.800 0.172 1.494 31.2
100 3.050 0.292 1.613 41.1
200 3.370 0.323 1.644 44.1

Resumen
Caudal máximo 20 m3/seg para Tr = 2 años
Caudal máximo 25 m3/seg para Tr = 5 años
Caudal máximo 28 m3/seg para Tr = 10 años
Caudal máximo 31 m3/seg para Tr = 50 años
Caudal máximo 41 m3/seg para Tr = 100 años
Caudal máximo 44 m3/seg para Tr = 200 años

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
35
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

METODO LOG-PEARSON TIII

3
Nº Orden Años Caudales Log X (logX - MedialogX)
1 1966 13.8 1.140 -0.0060
2 1967 19.5 1.290 0.0000
3 1968 25.6 1.408 0.0007
4 1969 17.9 1.253 -0.0003
5 1970 21.2 1.326 0.0000
6 1971 21.3 1.328 0.0000
7 1972 22.3 1.348 0.0000
8 1973 20.5 1.312 0.0000
9 1974 24.7 1.393 0.0004
10 1975 20.8 1.318 0.0000
11 1976 19.8 1.297 0.0000
12 1977 15.6 1.193 -0.0021
13 1978 20.2 1.305 0.0000
14 1979 24.6 1.391 0.0003
15 1980 22.8 1.358 0.0000
16 1981 21.5 1.332 0.0000
17 1982 18.7 1.272 -0.0001
18 1983 23.2 1.365 0.0001
19 1984 25.6 1.408 0.0007
20 1985 17.2 1.236 -0.0006
21 1986 19.5 1.290 0.0000
22 1987 15.4 1.188 -0.0024
23 1988 17.9 1.253 -0.0003
24 1989 25.6 1.408 0.0007
25 1990 19.5 1.290 0.0000
26 1991 34.7 1.540 0.0105
27 1992 16.8 1.225 -0.0009
28 1993 34.5 1.538 0.0101
29 1994 28.6 1.456 0.0025
30 1995 17.6 1.246 -0.0004
31 1996 27.4 1.438 0.0016
32 1997 14.7 1.167 -0.0037
33 1998 19.5 1.290 0.0000
Total = 43.603 0.0105

Cantidad de elementos 33
Media logX = 1.321
σlogX = 0.096

Coef. de asimetría G = n . Σ ( logX - Media logX )3 = 0.40

(n-1) . (n-2) . ( σlogX )
3

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
36
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Tr K K σlogX Media logX + K σlogX X

2 -0.066 -0.006 1.315 20.7
5 0.816 0.078 1.399 25.1
10 1.317 0.126 1.447 28.0
50 2.261 0.217 1.538 34.5
100 2.615 0.250 1.572 37.3
200 2.949 0.282 1.604 40.2

Resumen
Caudal máximo 21 m3/seg para Tr = 2 años
Caudal máximo 25 m3/seg para Tr = 5 años
Caudal máximo 28 m3/seg para Tr = 10 años
Caudal máximo 35 m3/seg para Tr = 50 años
Caudal máximo 37 m3/seg para Tr = 100 años
Caudal máximo 40 m3/seg para Tr = 200 años

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
37
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

METODOS GRAFICOS

Se realiza el reordenamiento de los caudales asignándole un numero de orden:

Años Caudales Nº orden Orden dec

1966 13.8 1 34.7
1967 19.5 2 34.5
1968 25.6 3 28.6
1969 17.9 4 27.4
1970 21.2 5 25.6
1971 21.3 6 25.6
1972 22.3 7 25.6
1973 20.5 8 24.7
1974 24.7 9 24.6
1975 20.8 10 23.2
1976 19.8 11 22.8
1977 15.6 12 22.3
1978 20.2 13 21.5
1979 24.6 14 21.3
1980 22.8 15 21.2
1981 21.5 16 20.8
1982 18.7 17 20.5
1983 23.2 18 20.2
1984 25.6 19 19.8
1985 17.2 20 19.5
1986 19.5 21 19.5
1987 15.4 22 19.5
1988 17.9 23 19.5
1989 25.6 24 18.7
1990 19.5 25 17.9
1991 34.7 26 17.9
1992 16.8 27 17.6
1993 34.5 28 17.2
1994 28.6 29 16.8
1995 17.6 30 15.6
1996 27.4 31 15.4
1997 14.7 32 14.7
1998 19.5 33 13.8

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
38
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Se calculan las frecuencias relativas según las expresiones de cada autor:

Según California : Según Tukey : Según Blom :

Fr = m Fr = 3.m - 1 Fr = m - 3/8
N 3.N + 1 N - 1/4

Nº Orden Años Caudales nº orden s/California s/Tukey s/Bloom

ord dec. Fr Tr = 1/Fr Fr Tr = 1/Fr Fr Tr = 1/Fr
1 1966 34.7 1 0.03 33.00 0.02 50.00 0.02 52.40
2 1967 34.5 2 0.06 16.50 0.05 20.00 0.05 20.15
3 1968 28.6 3 0.09 11.00 0.08 12.50 0.08 12.48
4 1969 27.4 4 0.12 8.25 0.11 9.09 0.11 9.03
5 1970 25.6 5 0.15 6.60 0.14 7.14 0.14 7.08
6 1971 25.6 6 0.18 5.50 0.17 5.88 0.17 5.82
7 1972 25.6 7 0.21 4.71 0.20 5.00 0.20 4.94
8 1973 24.7 8 0.24 4.13 0.23 4.35 0.23 4.30
9 1974 24.6 9 0.27 3.67 0.26 3.85 0.26 3.80
10 1975 23.2 10 0.30 3.30 0.29 3.45 0.29 3.40
11 1976 22.8 11 0.33 3.00 0.32 3.13 0.32 3.08
12 1977 22.3 12 0.36 2.75 0.35 2.86 0.35 2.82
13 1978 21.5 13 0.39 2.54 0.38 2.63 0.39 2.59
14 1979 21.3 14 0.42 2.36 0.41 2.44 0.42 2.40
15 1980 21.2 15 0.45 2.20 0.44 2.27 0.45 2.24
16 1981 20.8 16 0.48 2.06 0.47 2.13 0.48 2.10
17 1982 20.5 17 0.52 1.94 0.50 2.00 0.51 1.97
18 1983 20.2 18 0.55 1.83 0.53 1.89 0.54 1.86
19 1984 19.8 19 0.58 1.74 0.56 1.79 0.57 1.76
20 1985 19.5 20 0.61 1.65 0.59 1.69 0.60 1.67
21 1986 19.5 21 0.64 1.57 0.62 1.61 0.63 1.59
22 1987 19.5 22 0.67 1.50 0.65 1.54 0.66 1.51
23 1988 19.5 23 0.70 1.43 0.68 1.47 0.69 1.45
24 1989 18.7 24 0.73 1.38 0.71 1.41 0.72 1.39
25 1990 17.9 25 0.76 1.32 0.74 1.35 0.75 1.33
26 1991 17.9 26 0.79 1.27 0.77 1.30 0.78 1.28
27 1992 17.6 27 0.82 1.22 0.80 1.25 0.81 1.23
28 1993 17.2 28 0.85 1.18 0.83 1.20 0.84 1.19
29 1994 16.8 29 0.88 1.14 0.86 1.16 0.87 1.14
30 1995 15.6 30 0.91 1.10 0.89 1.12 0.90 1.11
31 1996 15.4 31 0.94 1.06 0.92 1.09 0.94 1.07
32 1997 14.7 32 0.97 1.03 0.95 1.05 0.97 1.04
33 1998 13.8 33 1.00 1.00 0.98 1.02 1.00 1.00

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
39
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Graficas de California, Tukey y Blom respectivamente:

80
70
60
50
y = 5.690ln(x) + 16.24
40
30
20
10
0
1.00 10.00 100.00 1000.00 10000.00

70
60
50
40 y = 5.281ln(x) + 16.32
30
20
10
0
1.00 10.00 100.00 1000.00 10000.00

70
60
50
40 y = 5.222ln(x) + 16.44
30
20
10
0
1.00 10.00 100.00 1000.00 10000.00

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
40
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

EJEMPLO PRÁCTICO CON EL USO DE LA REGRESIÓN LINEAL

Se trata de establecer la linealidad entre la precipitación anual sobre una cuenca (variable X)
y la escorrentía que esta genera sobre la misma (Variable Y):

n año lluvia anual m Escurr anual m

X Y
1 1966 1.637 0.557
2 1967 1.655 0.436
3 1968 1.884 0.562
4 1969 2.018 0.776
5 1970 1.536 0.600
6 1971 1.349 0.391
7 1972 1.327 0.359
8 1973 1.555 0.567
9 1974 1.405 0.635
10 1975 1.335 0.495
11 1976 1.391 0.580
12 1977 1.133 0.325

Se calculan los parámetros para hallar los términos de la recta

-0.026 + 0.362 X

Factor de correlacion : 0.71 La correlación es fiable

Luego si quisiera conocer cual sería la escorrentía para una precipitación de 2.5m = 0.9m

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
41
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

EJEMPLO PRÁCTICO CON EL USO DE LA REGRESIÓN NO LINEAL

Se utiliza en los casos donde existe una correspondencia entre variables del tipo no lineal.

Ecuacion No Lineal
β
y = α.x
ln y = lnα + β.lnx
y = e^ (lnα + β.lnx)

Donde

a=
Σ y '⋅Σ ( x ' ) 2
− Σ ( x '⋅ y ' ) ⋅ Σ x '
n ⋅ Σ ( x ' ) 2 − (Σx ' ) 2
n ⋅ Σ ( y '⋅ x ' ) − Σ y '⋅ Σ x ' ρ = β . σx /σy
β
n ⋅ Σ ( x ' ) 2 − (Σ x ' ) 2
Coeficiente de correlacion

Con :
ln y = y'
ln α = a
ln x = x'

Luego se calculan los parámetros:

Entonces
y' = a + β.x'

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]
 UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
42
TRTATAMIENTO PROBABILISTICO DE LA INFORMACION UT2 de
HIDROLOGICA 42

HIDROLOGÍA Y OBRAS HIDRÁULICAS

Se trata entonces de establecer la correlación entre la curva de caudales de un río y la altura

del mismo.

n Caudal m3/s Altura m

X Y
1 531.0 2.5
2 294.0 1.5
3 288.0 1.5
4 159.0 0.8
5 1635.0 5.8
6 1705.0 6.0
7 1089.0 4.2
8 1560.0 5.6
9 937.0 3.8
10 1013.0 4.1
11 616.0 2.6
12 210.0 1.1
13 201.0 1.0
14 146.0 0.7
15 120.0 0.5
16 111.0 0.5
17 81.0 0.5
18 449.0 2.0
19 369.0 1.7
20 422.0 1.9
21 266.0 1.4
22 247.0 1.3
23 280.0 1.4

En este caso la ecuación a construir resolver es la siguiente:

n = 23.0
Σ y' = 12.4
Σ( x' )
2
= 826.5
Σ( x' y' ) = 88.9 a -4.586
Σ x' = 136.4 α=e
a
0.010
( Σ x' ) β
2
= 18599.2 0.864

Luego aplicando la ecuación

Se tiene que para un caudal de 5000 m3/seg la altura será de 16 m

Realizado por el Ing. DIEGO TAUS 2020 UT Nº2

REV. Colaboración Ing DAVID HARDOY AÑO
Archivo: [Link]