UNIVERSIDAD NACIONAL DE INGENIERÍA
1
POSGRADO-FIC UNI
MAESTRIA EN CIENCIAS EN INGENIERIA CIVIL CON MENCIÓN EN
HIDRAULICA
• MODELO SANS – Simulación estocástica
TRABAJO N° 8
DOCENTE:
MSc. WALTER OBANDO LICERA
ESTUDIANTES:
MUÑOZ HERMITAÑO, BETSY
CURSO: HIDROLOGIA AVANZADA Lima, 09 de Agosto del 2020
SUB CUENCA SHULLCAS
2
Stochastic
Stochastikos: Palabra Griega que hace referencia a algo que contiene una
variable aleatoria. En hidrología se usa para describir un sistema que tiene en el,
un elemento aleatorio.
INTRODUCCIÓN
A pesar de la existencia de paquetes sofisticados para el análisis matemático y
estadístico como STATGRAPHICS, MINITAB, SAS, MATLAB y otros; la simulación
de series hidrológicas requiere software especializado por dos razones fundamentales:
Por la naturaleza periódica de los procesos hidrológicos y porque algunas series
hidrológicas incluyen características complejas como dependencia en el tiempo y
memoria.
En ese contexto, durante los últimos 30 años se han desarrollado innumerables
programas de computadora para el modelamiemto de
series de tiempo hidrológicas y en 1996 elUSBR liberó, al mercado, el paquete SAMS,
escrito en lenguaje C y Fortran que opera en ambiente WINDOWS, y desde entonces,
se editaron las versiones SAMS 2000 y SAMS 2009 que incorporan las correcciones y
modificaciones sugeridas por los usuarios
PRINCIPALES CAPACIDADES Y LIMITACIONES DEL SAMS
1. Analiza datos anuales y estacionales (hasta 12 estaciones por año).
2. Incluye 3 tipos de opciones de normalización o transformación de datos 3
3. Permite trabajar con hasta 40 sitios y para propósitos de desagregación
multivariado con grupos de 10 sitios
4. Incluye esquemas para modelar sistemas complejos de redes
5. Máximo número de años de datos de entrada 600
6. El número de nuestras a generarse es ilimitado
7. El número de años a generarse es ilimitado
BREVE DESCRIPCION DEL SAMS
SAMS consta de muchas opciones de menú que permite al usuario eligir
varias alternativas de análisis disponibles. Consiste de tres módulos:
1.Análisis Estadístico de Datos
2. Fijación del Modelo Estocástico
3.Generación de Series Sintéticas
Análisis Estadístico de Datos “Statistical Analysis of Data”
• Ploteo de datos: ayuda a detectar tendencias, cambios, salidas, o errores en los
datos. 4
• Verificación de la normalidad y Transformación de datos: uando técnicas de
transformación Logarítmica, Potencial y de Box-Cox .
• Determinación de los parámetros estadísticos: media, desviación estándar,
coeficiente de asimetría, correlación serial (datos anuales), correlación de
estación a estación (datos estacionales), correlaciones cruzadas; y estadísticas
relativas a sequías, excedencias y requerimiento de almacenamiento. Módulo:
Fijación del Modelo Estocá
MANUAL DE SAMS 2007
PASO 1: Creación de datos
5
Para comenzar a trabajar con el SAMS se debe tener una serie de tiempo
completa, puede ser anuales o estacionales (mensual, estacional, etc), estos
datos se pueden extraer de una tabla en formato EXCEL.
ESTACION:ANGASMAYO
N° Años ene feb mar abr may jun jul ago set oct nov dic
1 1963 18 35.7 43.4 26 11.4 7.5 6.1 6.1 5.4 6 8 15.8
2 1964 7.1 17 39.5 20.7 9.5 6 4 3.8 4 4 5.7 4.2
3 1965 8.2 41.7 38.5 12.6 6.5 5 4 3.7 3.6 4.1 4.3 5.7
4 1966 12.7 11.3 19.1 6.7 5.2 3.7 3 3.3 4.1 7.4 11.1 33.6
5 1967 23.8 63 71.7 29.2 10.3 6.7 4.9 4.8 7.5 13.1 5.1 7.7
6 1968 24.4 18.9 37.1 13.3 5.5 4 4.5 6.7 4.5 3.7 10.4 11.4
7 1969 6.9 19 21.9 20.2 6.6 4.1 4.8 4.7 4.7 4.9 4.2 19.9
8 1970 48.8 43.7 25.5 23.3 13.1 7 5.8 7.3 7.5 5.7 5 12.4
9 1971 21.8 53 52.5 21.5 8.4 5.8 4.6 4.2 3.6 5.6 6 16.2
Para crear un
10 1972 35.8 47.2 72.7 55.2 19.3 9.1 6.9 6.1 9.1 8.8 5.3 11.3
11 1973 46.8 77.7 77.3 42.6 18.1 10.5 8.1 7.2 7.5 10 10.3 23.3
12 1974 51.8 95.6 47.9 27.4 11.5 7.5 5.6 5.2 4.7 4.7 4.8 4.8
13
14
1975
1976
6.8
36.4
18.7 54.6 18.5 12.9 6.1
57.8 46 20.8 8.9 7.5
4.6
9.9
4.1
14
4.1
7.4
4.9
4.2
5.5
4
13.5
4.4
archivo .dat podemos
15
16
17
1977
1978
1979
11.7
30.4
8.6
26.2 50 15.4 11.5 5.4
43.1 30.7 20 7.5 5.3
29.3 42.7 24.2 8.6 5.2
4.1
8.7
4.7
3.5
13
4.1
6.2
9
3.7
7.8
4.4
4.8
14.6
8
6
10
10.4
4.7
utilizar softwares
18
19
1980
1981
5.7
28.3
10.9 34.4 12.5 1.4 1.1
87.2 45.5 19.5 6.1 2.7
0.9
2.8
1.7
5.5
3.4
4.1
10.2
4.6
6.1
10.6
12.5
23.9
diversos.
20 1982 41.3 60 38.8 25.7 8.1 3.4 2.6 2.5 3.1 7.9 19.6 18.7
21 1983 16 11.7 26.9 14.8 5.5 2.5 1.8 1.7 1.7 2.2 1.5 2
22 1984 20.2 72.7 63.4 36.3 11.8 6 3.5 2.7 2.6 7.7 13.6 32.5
23 1985 19.8 37.1 38.1 33.2 14.8 10 6.5 6.8 9.1 6.9 6.5 9.9
24 1986 40.6 63.8 62.5 36.4 27.6 25.5 23.5 21.5 19.5 17.4 15.4 13.4
25 1987 30.8 27.8 19.5 8 3.8 2.6 2.8 2.3 2.6 2.9 3.9 7.3
26 1988 29.1 36.6 24.5 27.3 8.5 4.2 3.8 3.8 4.7 7.2 5 4.5
27 1989 26.9 38.1 42.9 25.1 9.3 7.6 5.7 4.2 3.3 4.4 4 3.7
28 1990 9.9 7.7 10.4 5.4 4.5 5.5 3.8 3.2 3.2 4.7 10.2 10.6
29 1991 12.3 8.9 25.2 11.8 8.2 4.5 4.8 3.7 4.3 4.4 4.7 4.3
30 1992 4.5 5 11.4 4.4 3 3 2.8 2.8 2.2 3 2.6 2.3
31 1993 7.8 22.6 37 22.1 11.9 5.5 4.1 3.6 4 4 12 27.5
32 1994 27.3 56.3 36.5 35.4 11.7 6.5 4.8 4.4 4.4 4.2 3.9 4.5
33 1995 9.5 12.1 32 13.7 5.4 3.6 3.5 3.1 3.3 3.9 4.5 5.1
34 1996 19 49.4 33.9 23.1 7.7 4.8 4.2 3.8 3.9 3.7 3.7 5.6
35 1997 14.4 53 29.8 8.2 5.3 3.9 3.4 3.2 3.2 3.6 5.9 13.5
36 1998 26.9 46.7 33.7 23.4 5.8 4.4 3.9 3.4 3.1 4.3 4.7 5.7
37 1999 11.7 56 58.4 32.4 13.4 5.8 4.7 4.1 5.1 6 4.9 7.7
38 2000 22.9 54.3 54.4 21.9 9.8 5.9 5.2 4.9 4.7 8.7 5.8 11.7
39 2001 65.9 43.5 75 28.4 10.4 7.2 6.2 5.7 6.2 6.2 7.3 8.9
40 2002 9.8 38.2 39.6 19.4 8.7 5.5 5.1 4.5 5.4 6.6 12.8 27.7
En nuestro caso ya tenemos los datos en el archivo .dat
PASO 2: importar datos
Abrimos el program Sams 2007
Abrimos nuevo proyecto
Importar el archivo de datos
PASO 3: ANALISIS DE LOS DATOS
El ploteo de la serie de tiempo puede ayudar a detectar de manera visual las 8
tendencias, saltos, outliers y errores en la data, para lograr esto se selecciona lo
siguiente:
El ploteo de la serie de tiempo se puede hacer de manera anual o para un mes especifico,
esto se define en la pestaña ANNUAL/SEASONAL, cada “Season”
9
El ploteo de la serie a nivel anual. Se observa un pico extremo en el año 1986
10
En la grafica se observa que los valores de autocorrelacion (ACF) están dentro del
intervalo de confianza (las dos líneas paralelas entrecortadas) lo cual garantiza que este
mes (Season 1) tiene observaciones con independencia serial, se debe realizar el mismo
grafico para todos los meses y verificar este comportamiento, en caso de ver que los
valores del ACF exceden el intervalo de confianza se debe revisar los datos, corregirlos y
explicar este comportamiento.
Los datos de entrada deben ser estacionarios, es decir libres de tendencias,
para garantizar esta característica se debe realizar el proceso de normalización
o estandarización de los datos, para ello se procede asi: 11
Cuando los datos son mensuales se debe seleccionar la opción Seasonal en el cuadro
de la izquieda (Transf. Bar), este proceso se debe hacer para todos los meses, es decir,
se repite 12 veces ya que la generación de series sintéticas a “seasonal” es decir, a
nivel mensual, se genera resultados y datos para cada mes independientes los unos de
los otros.
12
A nivel “seasonal” (mensual) seleccionar las opciones en la Transf. Bar tal
como se muestra en el gráfico:
13
14
El proceso de estandarización fue correcto, todos los meses pasaron los test de
confiabilidad como se muestra en el gráfico, en caso de tener algún resultado
“reject” se debe analizar la calidad de los datos, caso contrario se pueden buscar 15
otras opciones de normalización o en su defecto aceptar la normalización pero
explicando los motivos de este resultado y bajo su propio riesgo.
PASO 3: SELECCIONANDO EL MODELO ESTOCASTICO
Para proceder a elegir que modelo es el más adecuado para generar las series 16
sintéticas se debe primero ajustar los datos mensuales a las funciones que hemos
definido previamente, para ello se procede:
Se observa que ahora se variaron los parámetros
PAR y PMA (que ahora ya no es
0 sino 1), luego colocar Fit Model:
17
4.-GENERACIÓN DE SERIES SINTÉTICAS
Se procede primero a seleccionar el modelo más adecuado para el mes 18
elegido, la generación de series sintéticas se hace por mes de la siguiente
manera:
Luego se indica cuantos datos se desea generar, se recomienda un numero menor a la
longitud total de los datos, para este caso, como se tienen 50 años, se deberán generar
a lo mucho 50 datos.
19
Como este modelo SOLO ES PARA EL MES DE ENERO, se eligen las
siguientes opciones en el cuadro que aparece:
20
5) RESULTADOS
Las series sintéticas generadas se presentan en formato .txt en el directio de trabajo 21
(en la carpeta donde guardo el proyecto SAMS.
22
GRACIAS 23
