llidioloi>a Eslndistica pgina 156

l l i i l r o l o s E.stadslica p g i n a 157

6.4 DISTRIBUCION G U M B E L

con lo cual, la funcin densidad reducida Gumbel es:

La distribucin Gumbel, en una de las distribuciones de valor extremo, es llamada tanibiOn

-Y

Valor Extremo Tipo I , Fisher Tippett tipo I o distribucin doble exponencial.

-Y - c
g(Y) = EXP(-Y - EXP(-Y)) = e

1. FUNCIN A C U M U L A D A

y la fmcin acumulada reducida Gumbel, es;

La funcin de distribucin acumulada de la distribucin Gumbel, tiene la forma:

F(x) = EXP(-EXT(-(x - p) / a))

... (73)
,' :

-e

... (68)

G(Y) = EXP(-EXP(-Y)) = e

(X

- p)

F(x) = e ^

Los valores correspondientes x e Y, estn relacionados por:

y la relacin:

... (69)

para: - x < x < x

F(x) = G ( Y )
x- n
Y=
x= n +

donde:
O<a < x

... (74)

...

aY

3. ESTIM.ACIN D E P A R M E T R O S . M T O D O D E M O M E N T O S

, es el parmetro de escala

Utilizando el mtodo de momentos, se obtieneri las siguientes relaciones:

-X < p < X

, es el parmetro de posicin, llamado tambin valor central o moda

- Moda:
- Media:

2 . FUNCIN DENSIDAD
Derivando la funcin de distribucin acumulada, ec.(68), con respecto a x, se obtiene la
funcin densidad de probabilidad, es decir:
dF(x)
f(x) =
dx
I
f ( x ) = EXP(-(x - p)/a - EXP(-(x - pya))
a

xmoda = n
E(x) = X = n + a C

donde C, es la constante de Euler, cuyo valor es:

1
I
I
C = lim [ 1 + + + .... +
2
3
n
n ->x

..

(70)

- Inn]

C = 0.5772156649
por lo tanto:
X = 1.1 +0.57721 a

1
f(x)= a

-(X

-(x - p)/a
- p)/a - e

n2a2
. (71)

para:

. Varianza:

E[(x - E(x))2] = s2 =

(77)

de donde se obtiene:

-X < X < X

La variable aleatoria reducida Gumbel, se define como:

a =

S = 0.78 S

(78)

X- p

Y=

.... (76)

... (72)

H = X - 0.57721 a = X - 0.45 S

.... (79)

llldriiliiul E t t i i d l i l l i * |I>(||IH | I

l l i i l i l o g i a Estadstica | i u i i i a 158

- Calculo de los parmetros de la serie de caudales;

4. APLICACIN EN H I D R O L O G A
La ley de Gumbel o ley de valores extremos, se utiliza generalmente para ajustar a una

Media:

X = 957.59

expresin matemtica, las distribuciones empricas de frecuencia de caudales mximos

Desviacin estndar:

S = 682.72

anuales, precipitaciones mximas anuales, etc.

- Clculo de los parmetros de la distribucin Gumbel:
El programa calcula estos parmetros utilizando las ees. (78) y (79), obteniendo:
Es importante verificar antes de aplicar esta distribucin de probabilidad que los coeficientes de asimetra y curtosis de la distribucin emprica sean del mismo orden que los

parmetro! de escala :

valores poblacionales.

parmetro de posicin . \i = 650.3237

Uno de los inconvenientes del uso de esta distribucin, es que en una distribucin doble

a = 523.3182

- Prueba de bondad de ajuste

exponencial, la variable puede tomar cualquier valor, por lo que se puede asignar proba-

Los resultados que se obtienen del programa scm:

lidades a valores negativos de la variable aleatoria, cuestin que resta significacin fsica a la

A =0.1454
A Q = 0.2525

aplicacin, debido a que las variables hidrolgicas toman solamente valores positivos o cero.

, para un nivel de significacin del 5 %

-Decisin

5. EJEMPLO

Siendo:

A = 0.1454 < A Q = O 2525

Se tiene el registro de caudales mximos de 29 aos, para la estacin 9-3 Angostura, como

se concluye que los datos se ajustan a una distribucin Gumbel, con un nivel de

se muestra en la tabla.

significacin del 5 %

En este ro se desea construir una presa de almacenamiento, calcular el caudal de diseo

para el venedor de demasas, con un periodo de retorno de 50 aos. Usar la distribucin
Gumbel.

1660
618
876
563
824
557

917
683
740
520
824
818

3800
934
1120
360
1230
1030

1410
779
610
367
522
418

2280
921
1150
658
581

2. Clculo del caudal de diseo, para un perodo de retorno de 50 aos.

1;
F(Q = q) = P ( Q < q ) = l
T
1
F(Q = q) = 1
50
F(Q = q) = 0.98 = F(y)
-e-y
F(y)=e

=0.98

SOLUCION:

-e-y = ln 0.98

1. Ajuste de los datos de la serie a la distribucin Gumbel

e-y = 0.020202707

1.1 Para crear el archivo con la serie de datos, se usa el programa del listado 6.

-y = In (0.020202707)

1.2 Para realizar la prueba de bonda de ajuste, se utiliza el programa del listado 9. Este

y = 3.9019

programa calcula adems los parmetros de la serie y los parmetros de la distribucin,

algunos clculos parciales que se obtienen, es comoje muestra^ _

pero de la ec. (72):

Q-p
y=

=3.9019

I l i i l i o l o B a E s l a d i s l i c a ngina 160

H i d r o l o g a E.s(adslca p y i n a Ifi!

Q = H + 3.9019xa
Q = 650.3237 + 3.9019 x 532.3182
Q = 2,727.38 m^/s

LISTADO

9. Programa en QuicIdSASIC, para realizar la prueba de SniirnovKolmogorov, para la distribucin Gumbel

y^^^^^^i^^^iii******^*****

'
'

* * * * * * *

* * * * * * * * * * *

Prueba de bondad de ajuste

a una Distribucin GUMBEL

y * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *

DEFDBL P-T
D I M X(60), P(60), F(60), Z(60), D(60), FL2(60), ZL2(60), FL3(60), ZL3(60)
D I M Y(60), FG(60), GY(60)
GOSUB ingresodatos
CLS
LOCATE 12, 20: PRINT "Por favor espere, realizando clculos"
FOR m = 1 TO n
P(m) = m / ( n + l )
NEXTm
'**** Clculo de los parmetros de la Distribucin Gumbel ****
PI# = 3.141592654#
ALFA = SQR(6) / PIs DESV
M U = X M - .5772156649# * ALFA
' * * Clculo de la Funcin Gumbel Acumulada ****
FOR i = 1 T O n
Y(i) = (X(i) - MU) / ALFA
FG(i) = EXP(-EXP(-Y(i)))
NEXTi
'Clculo del estadstico Smirnov-Kolmogorov
FOR i = 1 T O n
D(i) = ABS(P(i) - FG(i))
NEXTi
' Salida por impresora
F= 1
LPRINT
LPRINT TAB(23); ' T A B L A DE VALORES DE P(x), F(x) y D"
LPRINT
LPRINT TAB(15), STRING$(56, "=")

LPR1NTTAB(15),
m * X(ni) P{x)
LPRINT " Y ( l ) * F(x) * D(F-P)
LPRINT TAB(15); STR1NG$(56. "=")
FOR m = 1 TO n
LPRINT TAB( 15);"* ";
'
LPRINT USING "###*"; m;
LPRINT USING 'W#U#M * X(m);
LPRINT USING
* "; P(m);
LPRINT USING "#.## * "; Y(m);
LPRINT USING "#M## * "; FG(m);
LPRINT USING "#.####"; D(m);
LPRINT.." *" "
NEXTm
LPRINT T A B ( 15); STRING$(56,"=")
LPRINT : LPRINT

'

LPRINT U S I N G " parmetro de escala a:

UUMMMM
"; ALFA
LPRINT USING " parmetro de posicin p : #######.####"; NRJ
LPRINT
' Clculo del delta mximo
MAX=D(1)
FOR i = 1 TO n
IF D(i) < M A X THEN
ELSE
M A X = D(i)
DMAX = MAX
ENDIF
NEXTi
tablavalor:
CLS
LOCATE 2, 15: PRINT "Tabla de valores de significacin ("; STR!XGS(1, 224);")"
LOCATE 8, 24: PRINT "NIVEL DE SIGNIFICACION"
LOCATE 10, 19: PRINT "NS= 0.20
(1)";
PRINT SPACES(IO); "NS= O 10 (2)"
LOCATE 12, 19: PRINT "NS= 0,05
(3)";
PRINT SPACE$(10), "NS= 0.01
(4)"
LOCATE 19, 50: PRINT "Opcin : _ "
inkey:
0P$ = INKEYS: IF OPS = "" THEN GOTO inkey
IF OPS = " 1" THEN DO = 1.07 / SQR(n): GOTO estimacin
IF OPS = "2" THEN DO = I 22 / SQR(n): GOTO estimacin
IF OPS = "3" THEN DO = 1.36 / SQR(n): GOTO estimacin
IF OPS = "4" THEN DO = 1.63 / SQR(n): GOTO estimacin
CLS
LOCATE 12, 30: PRINT "OPCION NO V . ^ I D A "
LOCATE 17. 23 I N f UT "PULSE R E T U R N PARA CONTINUAR"; a$

Hidrolo|>i'a E s t a d s t i c a piigina 165

lliilioloRia Esladislica p g i n a 164

6.5 PROBLEMAS PROPUESTOS

En el caso que se ajusten a las distribuciones tericas mencionadas, calcular para c? 'a
una de las distribuciones:

L En la siguiente tabla, se muestra el registro de descargas medias diarias mxiiiiii'.

correspondientes a 28 aos, en m^/s.

La probabilidad de que un caudal anual, est comprendido entre 350.m3/s y 800

m3/s.

6597

9231

5153

10165

4813

8551

6229

4757

La probabilidad de que un caudal anual, sea menor o igual que 850 m^/s.

4898

El caudal anual que ocurre con un periodo de retorno de 10 aos.

7645

6654

4870

5219

4360

4672

4474

3794

5719

3437

7399

4870

3029

3851

4983

8155

8126

6456

4898

3. A partir de un registro de 70 aos, de volmenes anuales en M M 3 (millones de metros

cbicos), se desea realizar el ajuste a una distribucin log-normal de tres parmetros,

Realizar la prueba de bondad de ajuste, para ver si se ajustan a una distribucim gamm.i

usando la prueba Smirnov-Kolmogorov.

de dos parmetros.
En este sentido se calcul la funcin de distribucin emprica, utilizando la frmula de

Si se ajustan a esa distribucin, determinar la probabilidad de que una descarga media

Weibull, de los cuales en la siguiente tabla, se muestran 12 puntos:

mxima diaria, est comprendida entre 6500 in-^/s y 9300 m-'/s.

Q
(M.M3)

P(x) = m/(n+l)

Realizar la prueba de bondad de ajuste, para ver si se ajustan a una distribucin gamma
de tres parmetros, si el ajuste es bueno, determinar la descarga media diaria mxima

292.696

0056

33.036

para un periodo de retorno de 100 aos.

382 376
493.536

0.15S
0 197

644.736

0310

709.546

0.409

805.966

0.507

765

923.456

0.606

420

1058.286

0.704

530

1132.276

0.803
0.901
0.944

2. Dada la serie de caudales medios anuales, en m^/s , correspondiantes a un registro de 40

aos, los mismos que se muestran en la tabla que sigue

200
555
540

480
370
400

430
695
250

8S0
690
420

470
730
750

350

3;o
780

630
800
360

0.098 .

580,

520

550

610

590

536

690

548

1386.224

530

570

450

650

320

680

610

290

1700.716

Realizar la prueba de bondad de ajuste, para ver si los datos de la serie se ajustan a la

Adems, mediante un proceso computacional y utilizando el mtodo de mxima verosimi-

distribuciones:

litud, se estimaron sus parmetros, resultando:

parmetro de posicin : xo = -997.603

Normal

media de In (x - xo) = 7.4865

Log-normal de 2 parmetros

varianza de In (x - xo) = 0.0553

Log-normal de 3 parmetros

Ilidrutii|>ia EaluillMira

(tituliit ](\y

Ilidiolosiii Esludistica pKnn

con estos datos, se pide:

BIBLIOGR.AFIA

Indicar si la serie de voli'inicncs anuales, se ajusta a la distribucin log-normal de 3

Aliaga, Segundo "Conceptos Bsicos de Estadstica y Probabilidades". Deparatmento de
Ingeniera Agrcola. Instituto Tecnolgico de Costa Rica, Cartago - Costa Rica,
1982.

parmetros.

CONSULTADA

Hallar en forma grfica (en papel de distribucin log-normal), el valor del caudal que
Aliaga, Segundo . "Tratamiento de Datos Hidrometeorolgicos". Editorial Horizonte
Latinoamericano, Lima - Per, 1983.

satisface la siguiente relacin:

P ( Q < q ) = 0.9

4. En la siguiente tabla, se muestra el registro de los cudales mximos instantneos anuales,

en m^/s, de la estacin Guardia, de la cuenca del Ro Tempisque:

Ao

Caudal mximo anual

m^/s

Ao

Caudal mximo anual

m^/s

1962

273

1977

734

1963

1045

1978

134

1964

1875

1979

1706

1965

2876

19S0

2409

1966

1967

19S1

551

1967

1643

1982

597

1968

133

1983

80

1969

186

1984

606

1970

883

19S5

616

197J

3687

1986

287

1972

756

19S7

121

1973

186

19S8

350

1974

616

1989

315

1975

224

1990

294

1976

122

1991

3384

Determinar estadsticamente si stos caudales mximo anuales, se distribuyen segn la

distribucin Gumbel. Si as fuera, determinar el caudal mximo anual, con un perodo de
retorno de 50 aos.

Aliaga, Segundo. "Hidrologa Estadstica". Editorial Horizonte Latinoamericano, Lima Per, 1985. '
Amisial, Roger. "Distribuciones Tericas Comunmente Utilizadas". Centro Interamericano
de Desarrollo Integral de Aguas y Tierra (CIDIAT), Mrida - Venezuela, 1978.
Haan, Ch. "Statistical Methods in Hydrology". lowa, State University Press, U.S.A.,
1977.
Kreyszig, Erwin. "Introduccin a la Estadstica Matemtica: Principios y Mtodos".
' Editorial Limusa, Mxico, D.F., 1987.
Salas Jos. "Computer Workshop in Statistical Hydrology". Hydrology and Water
Resources Program, Colorado State University, Fort Collins, Colorado, U.S.A.,
1978.
Villn, Mximo. "Diseo de Capacidad de Embalses por el Mtodo Experimental - Teoria
del Range", Tesis para optar el grado de Magister Scientiae en Ingeniera de
Recursos de Agua y Tierra, Universidad Nacional Agraria "La Molina", Lima Per, 1983.
Villn, Mximo. "Distribuciones Usadas en Hidrologa". Departamento de Ingenieria
Agrcola, Instituto Tecnolgico de Costa Rica, 1985.
Yevjevich, Vujica. "Probabily and Statistics in Hydrology". Water Resources Publications,
Fort Collins, Colorado, U.S.A., 1972.
Zamanillo, Eduardo. "Aplicacin de Distribuciones Estadsticas a Variables Hidrolgicas".
Universidad Nacional del Litoral, Departamento de Hidrologa General y Aplicada,
Argentina, 1980.