Art Alca Tara
Art Alca Tara
Art Alca Tara
case of the flow rates, we used the Yonan station with a data of 34 years. Prior to the
model, was analyzed by the SAMS and HYDRACCESS programs.
The results of calibration of the model show a high efficiency with a Nash of 83 %,
after optimization of the parameters X1 =5.88 and X2 =0.72. The coefficient of
determination is highly reliable with a value of 0.83. In the validation of the model, it also
won a high efficiency of 81.3 % and a R2 coefficient of 0.86. In the sensitivity of the
model, it was determined that the most sensitive parameter is X2, therefore, has the
greatest influence on the results of the model.
The simulated flow rates for the stations of Porcon and Asuncin, determine that
there is a difference in flow rates during the course of the year. In the northern part of the
bedside of the basin flows are always present; on the other hand, in the southern part of the
bedside of the basin, the flows are a function of precipitation, droughts exist in the absence
of rain.
KEY WORDS: rainfall-runoff model, basin, calibration, simulation.
I.- INTRODUCCIN
El presente trabajo, est vinculado con la aplicacin del modelo GR2M de lluviaescorrenta en la cuenca del ro Jequetepeque, cuenca perteneciente a la vertiente del
pacfico, ubicado en la costa norte del Per.
Los datos en anlisis, pertenecen a estaciones distribuidas espacialmente con la mayor
cantidad de data (precipitacin y temperatura), y la estacin Yonan, con datos de caudales,
cuyo periodo corresponde a los aos 1965 1998.
Dada la problemtica de la deficiencia de datos de caudales, lo que se espera obtener en el
presente estudio es la aplicacin eficiente del modelo a la cuenca del ro Jequetepeque; el
problema se extiende a la mayor cantidad de cuencas y micro cuencas de la parte alta,
lugares en los cuales se presenta como informacin valiosa los datos de precipitacin y
temperatura.
El programa SAMS, ha permitido el anlisis de los datos y la generacin de caudales
sintticos con proyeccin de fenmenos del nio; as mismo, el programa Hydraccess ha
permitido distribuir la influencia espacial de las estaciones y delimitacin de influencia de
la intensidad de precipitacin espacial a lo largo de la cuenca.
Durante el proceso de utilizacin del modelo GR2M, se ha realizado los procesos de
calibracin para un periodo de 240 meses, y el periodo de validacin de 120 meses, como
parte del planteamiento de los objetivos para el estudio de la cuenca del ro Jequetepeque.
Con la finalidad de determinar el parmetro que tiene mayor influencia en los resultados
del modelo, se plante el proceso de sensibilidad del modelo con variaciones porcentuales
del 5% de incrementos. Base de estos procesos, se discute los resultados que arroja el
modelo mediante cuadros y figuras.
II.-ANTECEDENTES Y MARCO TERICO
2.1 Cuenca del ro Jequetepeque
El rea de estudio se encuentra en la parte norte del Per, ubicada entre las regiones de La
Libertad (provincias: Pacasmayo y Chepn) y Cajamarca hacia el este (Provincias de
Cajamarca, Contumaz, San Pablo y San Miguel). Hidrolgicamente pertenece a la
vertiente del Pacfico, con una superficie aproximada de 4443.63 km2, (INADE, 2002).
Donde
Q = cauda, en m3 /s.
A = rea en m2.
V = velocidad en m/s.
Figura N 2. Medida de la seccin de un ro.
(1)
P
(1)
evaporation
(2)
(3)
Production X1
store
(2)
(3)
S2
(4)
S
(4)
P2
(5)
P1
P3
(5)
(6)
Outside of
the basin
S X 1
S
1
X1
P1 P S S1
S1
X2
(7)
(8)
S1 1
S
1 1 1
X
1
X1
with
tanh
with
tanh
X1
P2 S 2 S
S2
1/ 3
S 3
1 2
X 1
P3 P1 P2
R1 R P3
60 mm
(6)
Routing
store
(7)
(8)
R2 X 5 .R1
R22
Q
R2 60
R R2 Q
4.0
5.0
Los datos de caudales, se obtuvo de la estacin hidrolgica Yonan, con datos histricos de
34 aos, desde el ao 1965 al ao 1998. Con la finalidad de estandarizar la data de las
estaciones que no contaban con datos desde el ao 1965, se realiz proyecciones y
simulaciones de precipitaciones sinttica a travs del programa SAMS.
Para el anlisis de lluvia escorrenta, se analiz los datos correspondientes a 4 estaciones
tambin por un periodo de 34 aos (1965-1998). Los datos a lo largo de la cuenca, han
sido trabajados previamente en el programa Hydraccess para determinar aplicar mtodos
de interpolacin y distribucin espacial de las estaciones. As mismo, el modelo GR2M, es
utilizado para la generacin de caudales mensuales a partir de la lluvia y
evapotranspiracin, (Tabla N1).
Tabla N 1. Distribucin de las estaciones en la cuenca Jequetepeque.
Nombre de la
Estacin
Llapa
Porcn
Asuncin
Contumaz
Yonan
Tipo de
estacin
Meteorolgica
Meteorolgica
Meteorolgica
Meteorolgica
Hidrolgica
Provincia
San Miguel
Cajamarca
Cajamarca
Contumaz
Contumaz
Distrito
San Miguel
Cajamarca
Asuncin
Contumaz
Yonan
Longitud
Latitud
-78.83028 -6.999722
-78.62666 -7.033889
-785155
-7.3121
-78.83028 -7.355833
-79.10028 -7.250278
Altitud
msnm
2900
3261
2194
2610
600
120
90
Y=X
90
60
60
30
30
0
.1
.5
10
30
50
70
90
98
99.5
99.9
.1
.5
10
30
50
70
90
98
99.5
99.9
120
Y = ln(X + a)
90
60
3.5
3
2.5
30
2
1.5
0
.1
.5
10
30
50
70
90
98
99.5
99.9
.1
.5
10
30
50
70
90
98
99.5
99.9
Las muestras estadsticas que arrojan el programa SAMS de las series originales y series
transformadas de la estacin en anlisis, se presentan la Tabla 2.
Tabla N 2. Resultados estadsticos del anlisis de caudales en la estacin Yonan.
Station Yonan
Parmetros
Mean
Series originales 310.5
Series
5.609
Transformadas
StDev CV
Skewness Min
Max acf(1)
168.1 0.5414 1.6553
79.40 970.3 -0.0166
0.5239 0.0934 -0.2398
4.387 6.879 0.0045
acf(2)
-0.0474
0.0045
Por otro lado, de acuerdo a los anlisis de los modelos Auto regresivos de Media Mvil,
que se presentan en la Tabla 3, cuyo principio es generar el presente o futuro en funcin de
lo que ha ocurrido en el pasado, los valores que nos arroja el programa para la estacin en
estudio, nos muestra el valor referencial de trabajo es el ARMA (1,0), por ser el menor
valor, dado que minimiza el error de pronstico de la matriz de suma de errores.
Tabla N 3. Resultados de los modelos autoregresivos.
Mean:
Variance:
AICC:
SIC:
Estacin 1
ARMA(1,0)
ARMA(2,0)
5.54386
5.54386
0.374783
0.374783
38.361
40.222
37.501
40.474
ARMA(3,0)
5.54386
0.374783
39.337
40.537
Season
1
2
3
4
5
6
7
8
9
10
11
12
Mean
3.09629
3.64282
4.0474
3.82215
3.13472
2.28819
1.5303
1.15403
0.830742
2.03154
2.52541
2.45322
Variance
0.662601
0.834269
0.838159
0.618956
0.340606
0.250824
0.365754
0.241532
0.517606
0.479121
0.303288
0.764442
PARMA(1,0)
AICC
SIC
27.017
28.217
25.233
26.433
17.224
18.423
7.880
9.080
4.199
5.399
-36.170
-34.970
-27.538
-26.338
-12.454
-11.254
-15.765
-14.566
29.659
30.859
17.792
18.991
20.932
22.132
Usando la base de datos en log-normal y con la ARMA (1,0), para la estacin en estudio,
se determina la generacin de series sintticas, para 100 aos, corriendo el modelo,
obtenemos los siguientes datos del grfico. El color azul, corresponde a la data de series
histricas y el de color rojo, son las series sintticas o series generadas. As mismo se
puede observar 4 grandes fenmenos del nio hasta el ao 2094, siendo el ms intenso en
los prximos 10 aos. Tabla 4 y Figura 6.
Figura N 6. Generacin de series sintticas de la estacin Yonan para 100 aos.
Annual Series
Historical
900
Generated
Flow
700
500
300
100
10
40
70
100
130
Year
10
11
29.7
44.9
41.9
174.8
Los resultados que nos arroja el modelo, en el proceso de calibracin con la optimizacin
del programa Solver, son ms aceptables, dado que el valor del Nash obtenido en este
proceso, tiene una eficiencia del 83% aplicado a un periodo de 240 meses, como se
observa en la Figura 9; al mismo tiempo que el valor de los parmetros se modifican
(parmetros ptimos): cuyos valores son:
X1 = 5.88
X2 = 0.72
Nash = 83
Figura N 9. Valor de los parmetros, calibrando el modelo.
Stationname
Catchmentarea (km)
Modelparameters
x1: Capacity of the production store (mm)
x2: Waterexchangecoefficient (mm)
Period
Length of the warm-up period (months)
Length of the test period (months)
Start date
End date
Efficiencycriteria (%)
Nash(Q)
Nash(VQ)
Nash(ln(Q))
Bias
Cuenca Jequtepeque
3642
Transformed
Real
5.88
357.81
0.72
0.72
12
240
01/1966
01/1986
83.0
83.8
83.0
95.6
As mismo, en forma grfica (Figura 10), se observa que los flujos simulados, respecto a
los observados, guarda una buena correlacin, sin embargo se observa algunos resultados
no muy satisfactorios en seis periodos, siendo tres de ellos muy pronunciados: mar-78
12
(Qsim = 34.62 y Qobs = 7.402), abr-78 (Qsim = 39.07 y Qobs = 14.44) y ene-83 (Qsim =
92.25 y Qobs = 60.78). En el resto de los periodos, explica su proceso de generacin de
lluvia en caudales, como se observa en la Figura 10.
Figura N 10. Hidrograma de lluvia, caudales observados y simulados.
En la Figura 11, se considera la distribucin a partir de 1965, donde observa los cambios
en la velocidad de llenado de la produccin de tanques, por lo que se observa poca
estabilidad en el modelo, donde las variaciones S, al principio del periodo en evaluacin,
estn influenciadas generalmente por las Precipitaciones (ejemplo enero del 67, 68, 69 y
70), presentando dos periodos de tendencia marcados de los aos 68 72 y de los aos 80
al 83.
Figura N 11. Evolucin del nivel S del reservorio de produccin.
13
14
Feb-95
Feb-96
Feb-94
Feb-93
Feb-92
Feb-91
Feb-90
Feb-89
Feb-88
Feb-87
Feb-86
Ene-85
Ene-86
Ene-84
Ene-83
Ene-82
Ene-81
Ene-80
Ene-79
Ene-78
1/3
Ene-77
Ene-76
Ene-75
Ene-74
Ene-73
Ene-72
Ene-71
Ene-70
Ene-69
Ene-68
Ene-66
Ene-67
2/3
Para este proceso, el periodo de validacin, corresponde a 120 meses, es decir 1/3 de la
data en estudio, y como se observa en la figura 14, el resultado del Nash, en este proceso,
result con una eficiencia alta del 81.3%, mantenindose los parmetros X1 =5.88 y el
X2=0.72.
Figura N 14. Periodo de validacin.
Cuenca Jequtepeque
Stationname
3642
Catchmentarea (km)
Modelparameters
Transformed
Real
x1: Capacity of the production store (mm)
5.88
359.26
x2: Waterexchangecoefficient (mm)
0.72
0.72
Period
Length of the warm-up period (months)
12
Length of the test period (months)
120
Start date
02/1986
End date
02/1996
Efficiencycriteria (%)
Nash(Q)
Nash(VQ)
Nash(ln(Q))
Bias
81.3
80.5
78.4
98.6
As mismo, como se puede observar en la Figura 15, los resultados son bastante aceptables,
dado que los caudales simulados, se ajustan a los caudales observados en todo este periodo,
esto es un reflejo de los periodos de la generacin de lluvias, por lo cual el modelo es
vlido y aplicable para esta cuenca.
15
Los resultados estadsticos de la lnea de tendencia, con la ecuacin lineal de los caudales
simulados, en funcin de los valores observados, mantienen una alta correlacin de 0.86,
por lo que nos indica que si existe validez del modelo, como se observa en la Figura 16.
Figura N 16. Lnea de tendencia y coeficiente R2, periodo de validacin.
sensible (Tabla 6), por lo que es el parametro que tiene mayor influencia en los resultados
del modelo para nuestra cuenca; para ello variamos cada uno de los parametros
manteniendo el otro fijo.
Este proceso ha sido trabajado con los caudales generados durante todo el periodo de
calibracion del modelo. En forma grfica, se observa en las Figuras 17, donde se grafica la
sensibilidad para cada uno de los parametros en estudio; la tendencia de uno de ellos
aumenta con mayor rapidez a medida que se incrementan los porcentajes.
Tabla N 6. Datos del proceso de sensibilidad de parametros
%
0.75
0.50
0.25
0.20
0.15
0.10
0.05
0.00
-0.05
-0.10
-0.15
-0.20
-0.25
-0.50
-0.75
Qsim1
35.5
22.6
17.9
17.3
17.0
16.8
16.8
17.0
17.4
18.0
18.8
19.9
21.4
35.6
43.8
Qsim2
47.3
35.9
25.7
23.6
22.1
20.2
18.8
17.0
15.3
14.1
12.5
11.4
10.0
4.7
1.3
Qsim1%
108.54
33.04
5.06
2.01
-0.16
-1.21
-1.18
0.00
2.35
5.91
10.66
17.34
26.15
109.54
157.67
Qsim2%
178.07
110.91
51.04
38.88
30.05
18.69
10.47
0.00
-10.13
-17.36
-26.55
-33.13
-41.45
-72.56
-92.65
existen glaciares en esta parte de la cuenca, sin embargo existe escorrenta aun en pocas
de estiaje. El rea de influencia de estos caudales son los distritos de Llapa, Cochn y
Cajamarca.
Figura N 18. Simulacin de caudales a partir de la precipitacin. Est. Porcn.
En la Figura 19, observamos que los caudales simulados en las pocas secas, no presentan
corrientes de agua, esto estara relacionado con las caractersticas del suelo y los procesos
de intemperismo, dado que las precipitaciones tambin son menores respecto a la estacin
anterior.
Figura N 19. Simulacin de caudales a partir de la precipitacin. Est. Asuncin.
La simulacin de caudales de la estacin Asuncin, corresponden para los ros San Juan,
Paigal y Quinuas, cuya rea de influencia es el distrito de San Juan, Jess, Asuncin, y
parte de Contumaz.
18
V.-CONCLUSIONES
1. La aplicacin del modelo GR2M, en la cuenca del rio Jequetepeque, presenta una alta
eficiencia, tanto en el periodo de calibracin (83%), como en el periodo de validacin del
modelo (81.3%), por lo que se ajusta a la aplicacin de este modelo satisfactoriamente a
nuestra cuenca en estudio.
2. Los parmetros de calibracin, adquieren valores de X1= 5.88 y X2=0,72. As mismo, el
anlisis de sensibilidad, nos muestra que el parmetro que tiene mayor influencia en los
resultados del modelo es el X2 (coeficiente de intercambio de aguas subterrneas), pues es
el ms sensible, frente a la variacin de pequeos cambios en los parmetros del modelo.
3. La aplicacin del modelo, dado su alta eficiencia para la cuenca del rio Jequetepeque, se
adecua a la generacin de caudales a partir de la informacin de precipitacin bsicamente
en la parte de las micro cuencas (parte alta), donde se adolece de informacin de caudales.
4. A nivel de cabecera de cuenca, se observa que existe una variacin de los caudales en la
parte norte de la cabecera de cuenca, respecto a una disminucin en los caudales para la
parte sur de la cabecera de cuenca del ro Jequetepeque, simulado en dos estaciones
ubicados en las partes ms altas de la cuenca en estudio.
VI.-AGRADECIMIENTOS
El agradecimiento muy cordial a mis docentes, quienes me orientaron en este maravilloso
mundo de la investigacin. As mismo el agradecimiento al CONCYTEC, por su apoyo en
mis estudios y las investigaciones sobre uno de los recursos base como son los recursos
hdricos.
VII.-REFERENCIAS BIBLIOGRFICAS.
1.
2.
3.
4.
5.
6.
7.
8.
9.
Sveinsson, O., Lee, T., Salas, J., Lane, W. y Frevert, D. (2009). Stochastic Analysis
Modeling and Simulation. Colorado State University. USA.
20