HYDROIILA (Manual) Autores:
Helms Kennedy Rodríguez Venegas.
[email protected] Cristhian Rogger Esquen Lozano.
[email protected] Asesores:
Ing. Juan Pablo García Rivera
[email protected] Doc. Hans Sánchez Tueros
[email protected]
Autores: Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
1. INSTALACIÓN
De los archivos contenidos en la carpeta del programa, deberá ejecutarse el que tiene el nombre asignado de “HydroIILA”. Se visualizará la ventana que se muestra a continuación en la cual se deberá hacer clic en “siguiente ” :
Posteriormente se mostrará la siguiente siguiente ventana en la cual se deberá hacer clic en siguiente.
Nota: Deberá instalarse en la carpeta destinada por defecto, caso contrario ciertos archivos de “ayuda” no podrán visualizarse.
Autores: Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Posteriormente se visualizará otra ventana en la cual se deberá hacer clic en “siguiente”, de esta manera se dará inicio al proceso de instalación del programa, tal como figura en la siguiente imagen:
En su escritorio se creará un acceso directo del programa, tal como se presenta en la siguiente imagen:
2. EJECUCIÓN
El programa será ejecutado como cualquier otro software, al hacerlo, se mostrará una ventana de bienvenida, tal como se muestra a continuación:
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
3. INTERFAZ
Se presenta imagen con la interfaz de trabajo de HydroIILA.
4. PROCESAMIENTO DE DATOS 4.1. Duración de la Precipitación:
La duración de la precipitación es la primera opción que se tiene que determinar en HydroIILA para el cálculo de la intensidad de la lluvia. Existen dos fórmulas para ello, una cuando la duración (t) es menor a tres horas y otra para “t” mayor o igual a tres.
Para t < 3 horas:
Para t >= 3 horas: Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
4.2. Datos de Entrada: a) Cuando t < 3h , se activará la opción “zona” y aparecerá un cuadro desplegable con las tres zonas establecidas por el método IILA, Costa Centro y Sur, Sierra y Costa Norte y Selva, tal como se detalla en la imagen siguiente:
Una vez elegida la zona en estudio, seguidamente se activará el recuadro de Subzona, donde se deberá determinar aquella donde se encuentre el proyecto en estudio. La imagen siguiente lo detalla:
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
b) Cuando t >= 3h , se activará directamente la opción “Subzona” y aparecerá un cuadro desplegable con las 38 subzonas establecidas por el método IILA, tal como se detalla en la imagen siguiente:
Una vez fijada la zona y Subzona en estudio, según sea el caso, se activarán los casilleros de los parámetros faltantes u otros cargados por defectos. Se deberá digitar la información solicitada por HydroIILA. Por ejemplo, seleccionaremos t > 3h y l a Subzona 123_1, tal como se muestra en la imagen siguiente, y los parámetros faltantes serán: cota, n, T y el intervalo.
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Definici ón d e los datos de Entrada:
b, Eg, kg, a y n: Son factores de la zona y Subzona en análisis. Cota: Elevación sobre el nivel del mar ( metros). Dc: Distancia a la cordilla (metros). Dm: Distancia al mar (metros). T: Es el periodo de retorno (años). t: Duración de la tormenta (minutos). Intervalo: Es el intervalo de tiempo en que se dividirá t.
Digitalizados los datos faltantes se hará clic en “ejecutar” para realizar los cálculos respectivos, tal como se puede observar en la imagen presente:
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Luego de ello se nos mostrará la matriz de resultados, tal como se presenta a continuación:
4.3. Gráficos
Esta opción se activará luego de “ejecutar” la i nformación ingresada a HydroIILA y permitirá visualizar los gráficos de hietrogramas y las curvas de Intensidad-duración y frecuencia (I-D-F), tal como se muestra a continuación:
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
4.4. Exportar
Exportar permitirá llevar la matriz de información a Excel, bastará con hacer un solo clic para realizar la tarea.
4.5. Limpiar
Esta herramienta limpiará todos los datos ingresados a HydroIILA y permitirá un nuevo procesamiento de datos.
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
4.6. Calculadora
Al hacer clic sobre esta herramienta nos mostrará la calculadora de Windows.
4.7. Salir
Esta opción permitirá salir por definitiva de HydroIILA.
5. AYUDAS
Dentro de ayudas podemos encontrar el mapa de las subzonas del Perú y la descripción del programa:
5.1. Mapa de subzonas o
El mapa de las subzonas en que se encuentra dividido el Perú se encuentra en formato JPG y KMZ.
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Mapa subzonas formato jpg.
Mapa subzonas formato kmz.
6. Acerca del programa
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
MÉTODO DEL INSTITUTO ITALOLATINOAMERICANO, SENAMHI Y UNI EN LA OBTENCIÓN DE PRECIPITACIONES E INTENSIDADES MÁXIMAS ANUALES
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
CONCEPTUALIZACIONES: En el año de 1982 se publicó el “Estudio de la hidrología en el Perú”, que nació de la cooperación técnica entre el Instituto Italo-Latinoamericano (IILA), el Servicio Nacional de Meteorología e Hidrología (SENAMHI) y la universidad nacional de Ingeniería (UNI), el cual entre sus múltiples alcances, permite determinar los máximos anuales de altura e intensidades de lluvia en todo el Perú, para lo cual subdividió al país en 7 zonas y 38 sub zonas pluviométricas teniendo como fundamento de la investigación la climatología en el país, las alturas de lluvia anual, así como las alturas de lluvia máxima diaria, en periodos de observaciones más o menos largos otorgadas por las 1377 estaciones meteorológicas que en ese entonces se tenía instaladas en todo el país. La fórmula de intensidad general, dada por el estudio es:
, = (1 + ´ ) − La cual se cumple para duraciones de lluvia comprendidas entre: 3 horas ≤ t ≤ 24 horas. Para duraciones de lluvia menores a tres horas, se tiene la siguiente ecuación:
, = (1 + ´ )( + )− El cálculo de la altura de lluvia o precipitación se puede calcular con la relación:
, =
ℎ,
Dónde:
, : Intensidad de la lluvia en “t” horas, para un determinado periodo de retorno “T”.
ℎ, : Altura de lluvia en “t” horas, para un determinado periodo de retorno “T”. T: periodo de retorno en años. t: duración de la lluvia en horas. b: factor de aumento de duración de lluvia en horas, b: 0.5 horas (costa centro y sur), b: 0.4 horas (sierra), b: 0.2 horas (costa norte y selva). Usar este factor solo para duraciones de lluvia menores a tres horas. K´: Constante definida por las características de cada zona pluviométrica, asume un valor único en todos los puntos de la zona seleccionada (ver cuadro 1).
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
a: Constante definida por las características de cada sub zona pluviométrica (ver cuadro 2). n: Constante que junto con el factor “a” definen cada sub zona pluviométrica, la diferencia es que “n” toma un valor único por sub zona mientras “a” es un factor variable (ver cuadro 2) El valor de “a”, es calculado mediante la ecuación:
() =
()
=
() 15.2
Dónde:
(): Ley de variación propia de cada sub zona pluviometrica, varía en función de “z”: el cual puede ser: altitud, distancia al mar o distancia a la cordillera (ver cuadro 1)
: valor “t” equivalente al día, el cual deducido para diferentes estaciones varía en torno a un valor medio de 15.2 horas
Nota: Los criterios de estimación de “n” y “a” y de la ley () presuponen que en la zona exista algún pluviógrafo. Lo cual no ocurre para l a mayor parte de las sub zonas, por lo que no queda otra cosa que:
En las sub zonas donde se carece del dato “n” (ver cuadro 2), asumir para este un valor comprendido entre aquellos que caracterizan las subzonas vecinas.
Fijado “n”, deducir la constante “a” con la relación ya mostrada, =
15.2
.
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
=
Cuadro 1. Subdivisión del territorio en zonas y sub zonas pluviométricas y valores de los parámetros Kg y , que definen las distribuciones de probabilidades de hg (máxima altura de lluvia anual) en cada punto de estas.
Zona
123
4
5a
5b
K´g
K´g = 0.553
K´g = 0.861
K´g = 11* εg^(-0.85)
K´g = 130* εg^(-1.4)
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Sub zona
εg
123-1
εg = 85
123-2
εg = 75
123-3
εg = 100 - 0.022
123-4
εg = 70 - 0.019
123-5
εg = 24
123-6
εg = 30.5
123-7
εg = -2 + 0.006 y
123-8
εg = 26.6
123-9
εg = 23.3
123-10
εg = 6 + 0.005 y
123-11
εg = 1 + 0.005 y
123-12
εg = 75
123-13
εg = 70
4.- 1
εg = 20
5 - a1
εg = -7.6 + 0.006 y (y> 2300)
5 -a2
εg = 32 + 0.177 Dc
5 - a3
εg = -13 + 0.010 y (y>2300)
5- a4
εg = 3.8 + 0.0053 y (y>1500)
5 - a5
εg = -6 + 0.007 y (y>2300)
5 -a6
εg = 1.4+ 0.0067
5 - a7
εg = -2 + 0.007 y (y>2000)
5 - a8
εg = 24+ 0.0025 y
5 - a9
εg = 9.4+ 0.0067 y
5 - a10
εg = 18.8+ 0.0028 y
5 - a 11
εg = 32.4+ 0.004 y
5 - a12
εg = 19.0+ 0.005 y
5 - a13
εg = 23.0+ 0.0143 y
5 - a14
εg = 4.0+ 0.010 y
5 - b1
εg = 4 + 0.010 y (y>1000)
5 -b2
εg = 41
5 - b3
εg = 23.0+ 0.143 y
5 -b4
εg = 32.4+ 0.004 y
5 - b5
εg = 9.4+ 0.0067 y
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
6 9
10
K´g = 5.4* εg^(-0.6)
K´g = 22.5* εg^(-0.85)
K´g = 1.45
6. - 1
εg = 30.0 - 0.50 Dc
9.- 1
εg = 61.5
9.- 2
εg = -4.5 + 0.323 Dm (30≤ Dm ≤110)
9.- 3
εg = 31+ 0.475 (Dm - 110) (Dm < 110) εg = 12.5 + 0.95 Dm
10.- 1
Cuadro 2. Parámetro “a” y “n” que junto con K´ definen las curvas de probabilidad pluviométrica en cada punto de las sub zonas. Sub zona
número total de estaciones
valor de n
valor de a
123-1 123-3 123-13 123-5 123-6 123-8 123-9 123-10 123-11 5 - a2 5 - a5 5 - a10
2
0.357
32.2
3
0.405
a= 37.85 - 0.0083 y
2
0.432
2
0.353
9.2
4
0.38
11
9
0.232
14
1
0.242
12.1
14
0.254
a= 3.01 + 0.0025 Y
5
0.286
a= 0.46 + 0.0023 y
2
0.301
a= 14.1 - 0.78 Dc
1
0.303
a= -2.6 + 0.0031 y
1
0.434
a= 5.8 + 0.0009 y
Cuadro 3.ESTACIONES METEREOLOGICAS: número, estado y entidad operadora. ENTIDAD SENAMHI Sector Agricultura Sector Energía y Minas Diversos TOTAL %
En Func.
Susp.
Total
1011
70
1081
31
77
108
106
49
155
11
22
33
1159
218
1377
84.2%
15.8%
100.0%
De las cuales solo se consideró en el estudio solo las estaciones de medida que disponen de un período de observaciones de al menos 6 años, la sub división según este criterio se muestra en el siguiente cuadro:
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Cuadro 4. Subdivisión de las estaciones pluviométricas en función de l a duración ng , en años, del periodo de observación. ng
Número de Estaciones
6 - 10 11 - 16 17 - 25 > 25
98 270 223 27
Delimit ación de las Zonas y sub zonas pluvio métricas Para la delimitación de las sub zonas, se ha hecho particular referencia al mapa ecológico del Perú a escala 1:1.000.000, publicado por la desaparecida ONERN (Oficina Nacional de Evaluación de Recursos Naturales), también se tomó como base de criterio para la subdivisión, el mapa climático de Trewartha para el Perú, del cual se tomó como referencia las situaciones meteorológicas que pueden verificarse dentro del territorio peruano, y se pudo distinguir: -
-
Un clima tropical A, que afecta a la selva amazónica (clima Ar , tropical húmedo) y la parte más al norte de la faja litoral del pacífico (clima Aw tropical húmedo y seco); Un clima B, desértico, que afecta a la parte restante de la faja litoral del pacífico; Un clima H, de las tierras altas, que afecta a todo el relieve andino.
De lo cual se pudieron distinguir 10 zonas pluviométricas: La Zona 1, en la selva amazónica, con el clima Ab, al sur y con el clima Ha, al norte. La zona 2, en la vertiente amazónica del relieve andino, con clima Ha. Las zonas 3 y 4, con la zona 4 que interrumpe la continuidad de la zona 3, en el altiplano con clima Hb. La zona 5 y 6, con la zona 6 que interrumpe la continuidad de la zona 5, en la vertiente del pacífico y la vertiente del norte (cuenca del marañón) del relieve andino con clima Hc. La zona 7, en la faja litoral del pacífico situada al sur 10° de lati tud sur con clima B. Las zonas 8 y 9, en la faja del pacífico inmediatamente al norte de la zona 7, con clima Ac. La zona 10, en el extremo norte de la faja litoral del Pacífico con clima Ad. Finalmente en base al valor K´ característico de cada zona, determinado por los registros pluviométricos, se definió K´= constante, soló para las zonas 1, 2, 3 ,4 y 10. Además desde el punto de vista estadístico se determinó que la diferencia
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
de los valores de K´ para las zonas 1, 2 y 3 no eran significativas, por lo que se les agrupo en una sola zona denominada con el símbolo 123. La zona 5 en referencia al valor K´ fue subdividida en las zonas 5a y 5b. Todo esto permitió al estudio concluir que el territorio puede ser subdividido en 7 zonas pluviométricas: 123, 4, 5a, 5b, 6, 9 y 10 dejando fuera las zonas 7 y 8 de la faja litoral, en las cuales la distribución de probabilidad de h g (altura de la lluvia diaria máxima) y ht (altura de la lluvia máxima en “t” horas) no tiene interés.
Fig1. Mapa de la división de las 7 zonas y 38 subzonas en todo el Perú.
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Ejemplo d e Aplicación: Se cuenta con los datos de registros históricos de lluvia máxima anuales en 24 horas de la estación Salpo, se desea conocer el porcentaje de variación de las precipitaciones e intensidades máximas en 24 horas, para diferentes periodos de retorno, los cuales serán de T: 2, 5, 10, 25, 50, 100 y 500 años. Obtenidas al comparar los resultados otorgados por el método IILA y por las funciones de distribución.
ESTACIÓN SAL PO La estación Salpo ubicada a 3412 m.s.n.m., con una Latitud de 8° 00’ 31.6”, Longitud de 78° 36’ 33.4” localizada en el distrito de Salpo, provincia de Otuzco, departamento de la Libertad, esta estación cuenta con un registro histórico de 46 años.
Método IILA-SENAMHI- UNI Para el cálculo de la intensidad, Se tiene la ecuación:
, = (1 + ´ ) − Además se debe cumplir que la duración de la lluvia, este comprendida entre: 3 ℎ ≤ ≤ 24 ℎ. Nuestra duración de lluvia es de 24 horas con lo que cumple con la ecuación para hallar nuestra intensidad. Conociendo que nuestra estación se encuentra dentro de la subzona 5a – 7, y con una elevación de 3412 m.s.n.m se pueden conocer los parámetros:
ɛ : 2 + 0.007(y) : 2 + 0.007 ∗ 3412 ∶ . : 11 ∗ ɛ −.: .
n: 0.33, se optó por tomar este valor, en relación del promedio entre las dos subzonas vecinas (5a 5 y 5a 10), que cuentan que con este valor “n” definido. :
ɛg :. 15,2
La ecuación de intensidad por el método IILA, para la estación Quiruvilca es:
, = 8.9149(1 + 0.7985 ∗ ) .− De igual manera se puede conocer la altura de precipitación por medio de la relación:
, =
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
ℎ
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Cuadro resumen de precipitaciones e intensidades máximas en 24 horas obtenidas por el método IILA: Probabilidad de excedencia
Periodo de intensidad máx precipitación máx retorno (años) 24 hrs (mm/hr) 24 hrs (mm)
0.5
2
1.32
31.56
0.2
5
1.65
39.65
0.1
10
1.91
45.76
0.04
25
2.24
53.85
0.02
50
2.50
59.96
0.01
100
2.75
66.08
0.002
500
3.35
80.28
Funciones de distribució n: Del registro histórico de 46 años de la estación Salpo, se tiene la siguiente tabla:
n
AÑO DE REGISTRO
PRECIP MAX EN 24 HRS (mm)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
1964
16.3 17.6 22.3 38.5 16.8 27 31.6 22.2 21 21.2 14.4 28.4 20.2 31.6 19.6 69.5 20.2 35.3 30.6 48.3 28.6 26.3 23.3 25 20 26.5
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46
21.5 20.6 22 25 40 23 35.5 26 28.5 31.8 24 29 29 29 20.8 16.4 33 20.4 43 30.5
1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009
Las funciones de distribución analizadas usadas en hidrología, fueron: la Función Normal, Función Log Normal, Función Gumbel, Función Log Gumbel, Función Pearson 3 o Gamma 3 parámetros y la Función Log Pearson 3. A las cuales se les aplico un ajuste de distribución gráfica y la prueba de ajuste de Kolmogorovsmirnov, el resumen de estas pruebas de ajuste se pueden observar en la siguiente tabla: Función de distribución Normal Log Normal Gumbel Log Gumbel Pearson 3 Log pearson 3
kolmogorov
Tabla
0.128
0.2
0.070
0.2
0.087
0.2
0.088
0.2
Inconsistente
0.2
0.065
0.2
De la cual se puede concluir que la función que mejor ajuste tiene por tener un delta calculado menor a los demás y por lo mismo también menor al delta de tabla, es la función Log Pearson 3. A continuación se detalla la prueba de ajuste Kolmogorov- Smirnov, para la función Log Pearson3 , los parámetros de esta función y la prueba de ajuste gráfico a la que también fue sometido.
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Distribució n de Probabilidad Log Pearson 3 – prueba de ajuste smirnov n
año
Pmax 24 (mm)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
1974
14.4
2.66722821
0.00406265
0.0212766
0.017213948
1964
16.3
2.79116511
0.03431222
0.04255319
0.008240967
2005
16.4
2.79728133
0.03706044
0.06382979
0.026769344
1968
16.8
2.82137889
0.04928486
0.08510638
0.035821519
1965
17.6
2.8678989
0.0794406
0.10638298
0.026942374
1978
19.6
2.97552957
0.18158487
0.12765957
0.053925299
1988
20
2.99573227
0.2051318
0.14893617
0.05619563
1976
20.2
3.0056826
0.21713854
0.17021277
0.046925778
1980
20.2
3.0056826
0.21713854
0.19148936
0.025649182
2007
20.4
3.0155349
0.22927047
0.21276596
0.016504511
1991
20.6
3.02529108
0.24150572
0.23404255
0.007463163
2004
20.8
3.03495299
0.25382323
0.25531915
0.001495917
1972
21
3.04452244
0.26620281
0.27659574
0.010392931
1973
21.2
3.05400118
0.27862515
0.29787234
0.019247191
1990
21.5
3.06805294
0.29729885
0.31914894
0.021850091
1992
22
3.09104245
0.32838834
0.34042553
0.012037193
1971
22.2
3.10009229
0.34076744
0.36170213
0.020934692
1966
22.3
3.10458668
0.34693776
0.38297872
0.036040963
1995
23
3.13549422
0.38962616
0.40425532
0.014629161
1986
23.3
3.14845336
0.40757559
0.42553191
0.017956326
2000
24
3.17805383
0.44842427
0.44680851
0.001615757
1987
25
3.21887582
0.50379779
0.46808511
0.035712683
1993
25
3.21887582
0.50379779
0.4893617
0.014436087
1997
26
3.25809654
0.55517365
0.5106383
0.044535356
1985
26.3
3.26956894
0.56976143
0.53191489
0.037846535
1989
26.5
3.27714473
0.57927156
0.55319149
0.026080068
1969
27
3.29583687
0.60229212
0.57446809
0.027824035
1975
28.4
3.34638915
0.66107134
0.59574468
0.065326658
1998
28.5
3.34990409
0.66495637
0.61702128
0.047935097
1984
28.6
3.35340672
0.66880063
0.63829787
0.030502754
2001
29
3.36729583
0.68377446
0.65957447
0.024199992
2002
29
3.36729583
0.68377446
0.68085106
0.002923397
2003
29
3.36729583
0.68377446
0.70212766
0.018353199
2009
30.5
3.41772668
0.73442907
0.72340426
0.011024811
1982
30.6
3.42100001
0.73751251
0.74468085
0.007168341
1970
31.6
3.45315712
0.76647285
0.76595745
0.000515399
1977
31.6
3.45315712
0.76647285
0.78723404
0.020761197
1999
31.8
3.45946629
0.77187214
0.80851064
0.036638494
2006
33
3.49650756
0.8017265
0.82978723
0.028060735
1981
35.3
3.56388296
0.84824322
0.85106383
0.002820611
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Z=LNX
f(Y)
F(M)
Kolmogorov . Smirnov
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
41 42 43 44 45 46
1996
35.5
3.5695327
0.85170964
0.87234043
0.020630787
1967
38.5
3.65065824
0.89474758
0.89361702
0.001130554
1994
40
3.68887945
0.91105135
0.91489362
0.003842263
2008
43
3.76120012
0.93602022
0.93617021
0.000149995
1983
48.3
3.87743156
0.96335528
0.95744681
0.005908473
1979
69.5
4.24132675
0.99466266
0.9787234
0.015939258
media des estan
3.25372979 7.20595557
coef asim
0.74504792 0.74504792
0.3061955
max calculado =
alfa beta
2.6843911
Z0
7.20595557 7.20595557 0.11406516 2.43178132
kOLMOGOROV (tabla) =
0.06532666
1
0.2
Análisis Gráf ico d e ajuste
Gráfico de Probabilidad 1.2 1 0.8 0.6 0.4 0.2 0 0
0.2
0.4
0.6 F logpearson3
0.8
1
1.2
Recta
Se concluye que los datos del registro histórico se ajustan a la distribución Log Pearson 3 y es además la función con mejor ajuste.
Precipit aciones e intensidades máximas en 24 horas por la función Log Pearson 3 Para los periodos de retorno antes mencionados, se tiene el siguiente cuadro: Probabilidad Periodo de de excedencia retorno (años)
intensidad máx 24 hrs (mm/hr)
precipitación máx 24 hrs (mm)
0.500
2
1.04
24.93
0.200
5
1.37
32.93
0.100
10
1.62
38.95
0.040
25
1.98
47.45
0.020
50
2.27
54.44
0.010
100
2.58
62.03
0.002
500
3.44
82.46
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Cuadro comparativo de precipitación LOG PEARSON vs IILA: Periodo de retorno (años) 2 5 10 25 50 100 500
precipitación máx IILA (mm)
precipitación máx Log Pearson 3 (mm)
Diferencia en ppton (mm)
% variación
31.56
24.93
6.63
21.01%
39.65
32.93
6.72
16.95%
45.76
38.95
6.81
14.88%
53.85
47.45
6.40
11.89%
59.96
54.44
5.52
9.21%
66.08
62.03
4.05
6.13%
80.28
82.46
-2.18
-2.71%
Cuadro comparativo de Intensidad LOG PEARSON vs IILA: Periodo de retorno (años) 2 5 10 25 50 100 500
Intensidad máx IILA (mm)
Intensidad máx Log Pearson 3 (mm)
Diferencia en Intensid (mm)
% variación
1.32
1.04
0.28
21.01%
1.65
1.37
0.28
16.95%
1.91
1.62
0.28
14.88%
2.24
1.98
0.27
11.89%
2.50
2.27
0.23
9.21%
2.75
2.58
0.17
6.13%
3.35
3.44
-0.09
-2.71%
Gráficos comparativos IILA vs Log Pearson3:
Precipitación (mm) 90 80 70 60 50 40 30 20 10 0 2
5
10
25 logper3
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
50
100
500
IILA
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros
Intensidad (mm/hr) 4 3.5 3 2.5 2 1.5 1 0.5 0 2
5
10
25 logper3
50
100
500
IILA
Se concluye que las precipitaciones e intensidades máximas halladas por el método IILA y la función Log Pearson 3, presentan poco porcentaje de variación para periodos de retorno superiores a 25 años, además se nota cierta tendencia de la función log Pearson a tener mayores valores de precipitaciones que el IILA conforme aumenta el periodo de retorno véase T: 500 años, se nota además una constante en la diferencia entre las precipitaciones e intensidades de ambos métodos, sobre todo para periodos de retorno comprendidos entre 2 y 100 años.
Autores: Helms Kennedy Rodríguez Venegas
Cristhian Rogger Esquen Lozano.
Asesores: Ing. Juan P. García Rivera
Doc. Hans Sánchez Tueros