Modelacion de Campos
Modelacion de Campos
Modelacion de Campos
INFORME
_____________________________________
Fecha: ______________________________________
TABLA DE CONTENIDOS
1. INTRODUCCIN .................................................................................................. 4
2. OBJETIVO............................................................................................................... 7
3. APLICACIONES GENERALES............................................................................ 8
4. DESCRIPCIN GENERAL DE LA ZONA DE ESTUDIO ............................... 9
4.1. EL LUGAR ........................................................................................................... 9
4.2. CONDICIONES OCEANOGRFICAS GENERALES ................................................ 9
4.2.1. Temperatura del Agua de Mar .................................................................... 9
4.2.2. Salinidad.................................................................................................... 13
4.2.3. Oxgeno Disuelto (OD)............................................................................. 13
4.2.4. Slidos en Suspensin, Turbidez y Transparencia .................................... 14
4.2.5. Mareas ....................................................................................................... 19
4.2.6. Corrientes .................................................................................................. 22
4.2.7. Vientos ...................................................................................................... 23
4.2.8. Olas ........................................................................................................... 24
5. INFORMACIN UTILIZADA PARA LA MODELACIN MATEMATICA
..................................................................................................................................... 26
5.1. MAREAS ........................................................................................................... 26
5.2. CORRIENTES ..................................................................................................... 30
5.2.1. Mediciones directas eulerianas ............................................................... 31
5.2.2. Corridas de flotadores lagrangeanas....................................................... 38
5.2.3. Modelos globales de corrientes................................................................... 45
5.2.4. Estado actual del conocimiento.................................................................. 47
5.3. VIENTOS ........................................................................................................... 48
5.4. BATIMETRA ..................................................................................................... 50
6. EL MODELO NUMRICO................................................................................. 51
6.1. GENERALIDADES .............................................................................................. 51
6.2. ANTECEDENTES. MODELACIN DEL CAMPO DE CORRIENTES EN LA ZONA . 52
6.3. MODELACIN 2D VS. 3D ................................................................................ 52
6.4. SISTEMA DE COORDENADAS Y ORIENTACIN DE LA GRILLA ........................ 54
6.5. REFERENCIA TEMPORAL Y MANEJO DEL TIEMPO ........................................... 55
6.6. PERODO DE ESTABILIZACIN ......................................................................... 55
6.7. NMERO DE COURANT.................................................................................... 56
6.8. BATIMETRA ..................................................................................................... 56
6.9. FORZANTES ...................................................................................................... 56
6.9.1. Bordes ........................................................................................................ 56
6.9.2. Mareas ....................................................................................................... 57
6.9.3. Corrientes .................................................................................................. 58
6.9.4. Vientos ...................................................................................................... 58
6.9.5. Forzantes combinados................................................................................ 59
6.10. VISCOSIDAD TURBULENTA ............................................................................ 59
6.11. RESISTENCIA................................................................................................... 59
6.12. CALIBRACIN DEL MODELO.......................................................................... 60
7. ESCENARIOS DE MODELACIN Y RESULTADOS.................................... 67
7.1. ESCENARIOS DE MODELACIN........................................................................ 67
7.2. RESULTADOS .................................................................................................... 68
CONCLUSIONES..................................................................................................... 88
1. INTRODUCCIN
provientes
de
la
combinacin
de
las
fuerzas
citadas
Las corridas del modelo, esto es, series de tiempo para cada punto del
dominio de modelacin, quedan almacenadas en soporte magntico
para posterior uso en otras aplicaciones.
2. OBJETIVO
3. APLICACIONES GENERALES
Los modelos, una vez calibrados pueden ser aplicados a problemas del
proyecto, por ambientalistas, diseadores, autoridades, contratistas,
etc. Podrn ser consultados (corridos) y actualizados en cualquier
momento.
4.1. El Lugar
La zona de estudio corresponde a la Baha de Paracas, ubicada en la
Provincia de Pisco, Departamento de Ica, Regin Los LibertadoresWari. La baha de Paracas y alrededores, se denomina "zona de Pisco",
por las caractersticas morfolgicas del litoral y la batimetra.
los
hundimientos
movimientos
9
advectivos.
Las
10
11
Tabla 4.2.1.1. Temperatura del agua de mar en las Bahias de Pisco y Paracas
Punto de muestreo
Profundidad (m)
Temperatura (C)
T1(0)
T1(5)
T1(5)
T1(10)
T1(10)
T2(0)
T2(5)
T2(5)
T2(10)
T2(10)
T3(0)
T3(5)
T3(5)
T3(10)
T3(10)
T4(0)
T4(5)
T4(5)
T4(10)
T4(10)
T5(0)
T5(5)
T5(5)
T5(10)
T5(10)
T5(15)
T5(15)
T6(0)
T6(5)
T6(5)
T6(10)
T6(10)
T7(0)
T7(5)
T7(5)
T7(10)
T7(10)
T8(0)
T8(5)
T8(5)
T8(10)
T8(10)
T9(0)
T9(5)
T9(5)
T9(10)
T9(10)
T10(0)
T10(5)
T10(5)
T10(10)
T10(10)
T11(0)
T11(5)
T11(5)
T11(10)
T11(10)
T12(0)
T12(5)
T12(5)
T12(10)
T12(10)
T13(0)
T13(5)
T13(5)
T13(10)
T13(10)
T14(0)
T14(5)
T14(5)
T14(10)
T14(10)
Superf
Superf
4.7
Superf
10.0
Superf
Superf
4.8
Superf
9.8
Superf
Superf
5.5
Superf
9.5
Superf
Superf
5.1
Superf
9.4
Superf
Superf
5.6
Superf
9.9
Superf
15.5
Superf
Superf
5.1
Superf
9.5
Superf
Superf
4.9
Superf
9.7
Superf
Superf
4.5
Superf
6.8
Superf
Superf
4.7
Superf
6.7
Superf
Superf
5.2
Superf
8.1
Superf
Superf
3.0
Superf
4.2
Superf
Superf
4.3
Superf
8.6
Superf
Superf
5.2
Superf
9.2
Superf
Superf
4.4
Superf
9.1
23.0
23.3
22.3
21.4
19.8
24.1
24.6
22.4
23.9
20.0
26.7
24.7
23.1
23.8
18.2
29.1
25.2
22.8
24.7
19.1
25.2
24.6
23.5
23.8
21.5
23.7
18.6
26.6
25.7
23.8
25.6
23.5
26.5
26.3
24.3
26.5
19.6
29.5
26.4
26.2
26.3
21.8
29.5
26.3
25.6
26.2
23.6
26.4
25.1
22.5
25.6
21.3
28.5
26.7
26.6
26.5
25.2
21.8
21.3
16.7
21.1
16.1
20.8
22.4
19.6
21.1
15.5
16.5
17.6
17.4
19.1
17.6
23
23.3
24.18
23
25.04
24.64
26.04
26.24
24.18
26.7
19.4
19.8
17.64
12
4.2.2. Salinidad
Los valores de salinidad medidos, oscilaron entre un mnimo de 3.02%
en superficie en el transecto 5 (cota de 10 m), y un mximo de 3.24% a
10 m de profundidad en el transecto 1 (cota de 10 m), con un promedio
de 3.16% (EIA, Lnea de Base Ambiental, Oceanografa, ERM, 2002). La
Tabla 4.2.2.1. resume los valores obtenidos en el campo.
Tabla 4.2.2.1. Salinidad del agua de mar en las Bahias de Pisco y Paracas
Punto de muestreo
T1(10)
T1(10)
T2(10)
T3(10)
T4(10)
T5(10)
T5(10)
T6(10)
T7(10)
T8(10)
T9(10)
T10(10)
T11(10)
T12(10)
T13(10)
T14(10)
Profundidad (m)
Superf
10.0
Superf
Superf
9.4
Superf
9.9
Superf
Superf
Superf
6.7
Superf
Superf
Superf
Superf
Superf
Salinidad (%)
3.23
3.24
3.18
3.17
3.23
3.02
3.18
3.12
3.11
3.11
3.21
3.13
3.15
3.19
3.22
3.21
13
14
Tabla 4.2.3.1. Niveles de Oxgeno Disuelto en el agua de mar en las Bahias de Pisco y
Paracas
Punto de muestreo
Profundidad (m)
T1(0)
T1(5)
T1(5)
T1(10)
T1(10)
T2(0)
T2(5)
T2(5)
T2(10)
T2(10)
T3(0)
T3(5)
T3(5)
T3(10)
T3(10)
T4(0)
T4(5)
T4(5)
T4(10)
T4(10)
T5(0)
T5(5)
T5(5)
T5(10)
T5(10)
T5(15)
T5(15)
T6(0)
T6(5)
T6(5)
T6(10)
T6(10)
T7(0)
T7(5)
T7(5)
T7(10)
T7(10)
T8(0)
T8(5)
T8(5)
T8(10)
T8(10)
T9(0)
T9(5)
T9(5)
T9(10)
T9(10)
T10(0)
T10(5)
T10(5)
T10(10)
T10(10)
T11(0)
T11(5)
T11(5)
T11(10)
T11(10)
T12(0)
T12(5)
T12(5)
T12(10)
T12(10)
T13(0)
T13(5)
T13(5)
T13(10)
T13(10)
T14(0)
T14(5)
T14(5)
T14(10)
T14(10)
Superf
Superf
4.7
Superf
10.0
Superf
Superf
4.8
Superf
9.8
Superf
Superf
5.5
Superf
9.5
Superf
Superf
5.1
Superf
9.4
Superf
Superf
5.6
Superf
9.9
Superf
15.5
Superf
Superf
5.1
Superf
9.5
Superf
Superf
4.9
Superf
9.7
Superf
Superf
4.5
Superf
6.8
Superf
Superf
4.7
Superf
6.7
Superf
Superf
5.2
Superf
8.1
Superf
Superf
3.0
Superf
4.2
Superf
Superf
4.3
Superf
8.6
Superf
Superf
5.2
Superf
9.2
Superf
Superf
4.4
Superf
9.1
7.7
7.5
6.5
10.3
5.9
5.8
6.3
4.5
9.1
5.1
14.2
7.7
5.7
8.5
3.8
18.5
8.6
4.2
9.0
4.4
7.2
6.9
6.6
8.2
5.7
8.5
6.5
11.0
7.2
4.1
8.7
5.8
6.1
8.1
6.1
10.8
7.2
10.2
10.3
10.2
9.7
2.8
12.8
8.7
5.5
8.6
6.1
6.9
6.9
5.5
8.6
5.5
11.5
9.9
9.5
10.2
5.8
6.4
4.3
3.1
6.1
4.3
6.0
6.3
5.1
6.6
4.5
4.1
4.7
4.4
5.4
4.9
15
6.16
7.98
8.94
7.08
7.36
7.66
8.64
8.34
6.68
9.38
4.84
5.7
4.7
Tabla 4.2.4.1. Turbidez del agua de mar en las Bahias de Pisco y Paracas
Punto de muestreo
Profundidad (m)
Turbiedad (NTU)
T1(0)
T1(5)
T1(5)
T1(10)
T1(10)
T2(0)
T2(5)
T2(5)
T2(10)
T2(10)
T3(0)
T3(5)
T3(5)
T3(10)
T3(10)
T4(0)
T4(5)
T4(5)
T4(10)
T4(10)
T5(0)
T5(5)
T5(5)
T5(10)
T5(10)
T5(15)
T5(15)
T6(0)
T6(5)
T6(5)
T6(10)
T6(10)
T7(0)
T7(5)
T7(5)
T7(10)
T7(10)
T8(0)
T8(5)
T8(5)
T8(10)
T8(10)
T9(0)
T9(5)
T9(5)
T9(10)
T9(10)
T10(0)
T10(5)
T10(5)
T10(10)
T10(10)
T11(0)
T11(5)
T11(5)
T11(10)
T11(10)
T12(0)
T12(5)
T12(5)
T12(10)
T12(10)
T13(0)
T13(5)
T13(5)
T13(10)
T13(10)
T14(0)
T14(5)
T14(5)
T14(10)
T14(10)
Superf
Superf
4.7
Superf
10.0
Superf
Superf
4.8
Superf
9.8
Superf
Superf
5.5
Superf
9.5
Superf
Superf
5.1
Superf
9.4
Superf
Superf
5.6
Superf
9.9
Superf
15.5
Superf
Superf
5.1
Superf
9.5
Superf
Superf
4.9
Superf
9.7
Superf
Superf
4.5
Superf
6.8
Superf
Superf
4.7
Superf
6.7
Superf
Superf
5.2
Superf
8.1
Superf
Superf
3.0
Superf
4.2
Superf
Superf
4.3
Superf
8.6
Superf
Superf
5.2
Superf
9.2
Superf
Superf
4.4
Superf
9.1
19.82
16.20
12.77
2.55
6.33
15.33
11.74
12.48
4.49
7.01
7.00
9.66
8.31
2.96
3.12
1.31
6.61
4.47
6.00
13.25
2.68
1.51
0.65
1.97
2.76
2.59
3.59
0.52
0.75
0.77
1.28
3.59
0.22
1.06
2.57
0.91
1.92
0.74
1.20
2.08
2.48
13.00
1.25
1.53
2.22
1.65
2.79
0.48
0.46
1.24
0.97
4.23
3.81
1.34
1.36
1.41
1.98
1.77
1.02
5.48
0.70
6.14
1.36
0.14
0.64
0.65
1.22
4.11
0.18
0.44
0.14
0.17
10.2
6.2
6.3
2.3
1.4
1.3
3.9
1.8
1.5
0.8
16
Se realiz un estudio (Off Shore Site Investigation Required for Sub Sea
Pipe LPG Pisco Terminal, Camisea Proyect, Pisco. Task 5000 Turbidity;
Golder Associates Per S.A., December 2002.) complementario de niveles
de turbidez y slidos totales disueltos a lo largo de la traza de la obra.
Los resultados de las mediciones se pueden observar en las Figura
4.2.4.1 y Tabla 4.2.4.2
17
18
4.2.5. Mareas
19
Durante la presencia del fenmeno de "El Nio", los valores del nivel
medio del mar se incrementan entre 30 a 40 centmetros.
Las mareas que llegan a las costas del Per procedentes del norte, es
decir que si una pleamar pasa por un determinado punto, despus de
un tiempo pasar por otro punto ms al sur. Por lo que la hora de las
pleamares y las bajamares va a ser diferente a todo lo largo del litoral.
20
Figura 4.2.5.1. Posiciones de las estaciones de muestreo a lo largo de la traza de la obra, Playa Lobera, Pisco
21
Puerto de Pisco
latitud:
13 43.0S
longitud:
76 14.0W
0.58 m,
Puerto El Callao
latitud:
12 03.0S
longitud:
77 09.0W
0.55 m,
22
Gran parte del flujo de agua en la franja comprendida entre 100 y 200 Km de
la costa peruana es predominantemente hacia el sur, contrariamente a los
vientos imperantes y a la circulacin asociada al sistema de Per-Chile.
Las corrientes costeras son dominadas por las mareas y los vientos.
4.2.7. Vientos
23
Tabla 4.2.7.1. Relacion velocidad del viento y direccion y sufrecuencia de ocurrencia, para el
Aeropuerto de Pisco.
4.2.8. Olas
costa
afuera,
provenientes
de
la
ejecucin
de
modelos
meteorolgicos globales.
24
Tabla 4.2.8.1. Relacin entre la altura significativa y periodo de las olas y su nmero de
ocurrencia, para el sitio de implantacin de la obra.
Figura 4.2.8.2. Relacion altura significativa y direccion de las olas y su numero de ocurrencia,
para el sitio de implantacin de la obra.
25
5.1. Mareas
En la Tabla 5.1.1 se presentan los valores correspondientes a las amplitudes de
las cuatro componentes de marea principales en Pisco y Callao obtenidas de
las Admiralty Tide Tables. Puede observarse que en la regin estudiada la
componente de marea predominante es la semidiurna principal (M2)
siguindole en importancia la diurna (K1).
Tabla 5.1.1. Amplitudes (en metros) de las principales componentes armnicas de la marea en
Pisco y El Callao. Periodo de la componente en horas: T.
Puerto
Pisco
Callao
M2
T=12.42060122
H=0.22
H=0.23
S2
T=12.00000000
H=0.07
H=0.08
K1
T=23.93446966
H=0.15
H=0.14
O1
T=25.81934166
H=0.07
H=0.07
y, con los valores de la Tabla 5.1.1 se obtiene que es mayor que 0.6 en ambos
casos, lo cual indica que la marea en este sector de la costa peruana es mixta
preponderantemente semidiurna (0.25 F 1.50).
Esto indica que durante un mes lunar, habr perodos con caractersticas
semidiurnas (dos pleamares y dos bajamares diarias) y otros con
caractersticas diurnas (una pleamar y una bajamar diaria).
26
Camisea Project Pisco Per Task 1000 - Ocean Currents, marzo de 2003) entre el 5
y el 15 de marzo de 2003. Por ejemplo, se observa que el 7 de marzo los niveles
del mar presentan caractersticas semidiurnas (dos pleamares y dos bajamares
por da) mientras que el 11 de marzo la marea presenta un claro
comportamiento diurno (una pleamar y una bajamar diaria), como indica la
Figura 5.1.1.
60
50
40
30
20
10
03/3/14 1:01
03/3/14 13:01
03/3/13 1:01
03/3/13 13:01
03/3/12 1:01
03/3/12 13:01
03/3/11 1:01
03/3/11 13:01
03/3/10 13:01
03/3/10 1:01
03/3/9 1:01
03/3/9 13:01
03/3/8 1:01
03/3/8 13:01
03/3/7 1:01
03/3/7 13:01
03/3/6 1:01
03/3/6 13:01
03/3/5 13:01
03/3/5 1:01
Fecha y Hora
Figura 5.1.1. Mediciones de niveles del mar en Pisco. Alguna singularidad en la curva se
asocia a un error en el procesamiento o dficit en el funcionamiento del sistema de obtencin de
datos.
27
Comparacin entre el nivel del mar predicho por la Tabla de Marea en El Callao
y los resultados del Modelo Global de Mareas TOPEX
1.2
1.1
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
02/6/1 0:00
02/6/1 9:00
02/6/1 18:00
02/6/2 3:00
02/6/2 12:00
02/6/2 21:00
02/6/3 6:00
02/6/3 15:00
02/6/4 0:00
02/6/4 9:00
02/6/4 18:00
02/6/5 3:00
02/6/5 12:00
02/6/5 21:00
02/6/6 6:00
02/6/6 15:00
02/6/7 0:00
02/6/7 9:00
02/6/7 18:00
02/6/8 3:00
02/6/8 12:00
02/6/8 21:00
02/6/9 6:00
02/6/9 15:00
02/6/10 0:00
02/6/10 9:00
02/6/10 18:00
02/6/11 3:00
02/6/11 12:00
02/6/11 21:00
02/6/12 6:00
02/6/12 15:00
02/6/13 0:00
02/6/13 9:00
02/6/13 18:00
02/6/14 3:00
02/6/14 12:00
02/6/14 21:00
02/6/15 6:00
02/6/15 15:00
0.0
Fecha y Hora
Figura 5.1.2. Mediciones de niveles del mar en El Callao y el clculo del Modelo Global Topex
para el mismo lugar. El modelo sigue la prediccin de tabla, basada en mediciones del lugar, en
amplitud y fase, correctamente.
28
Tanto las alturas como los instantes de las pleamares y bajamares publicadas
en la Tabla de Marea que edita la Direccin de Hidrografa y Navegacin de la
Marina del Per, se utilizaron para evaluar la calidad de los datos de marea
del modelo global. Las Figuras 5.1.2. y 5.1.3. muestran las salidas del TOPEX
para los puertos de El Callao y Pisco y los valores predichos por la Tabla de
Marea para esos lugares para un mismo perodo.
Comparacin entre el nivel del mar predicho por la Tabla de Marea en Pisco
y los resultados del Modelo Global de Mareas TOPEX
1.0
0.9
0.6
0.5
0.4
0.3
0.2
0.1
02/6/7 9:00
02/6/7 21:00
02/6/7 15:00
02/6/7 3:00
02/6/6 21:00
02/6/6 9:00
02/6/6 15:00
02/6/6 3:00
02/6/5 21:00
02/6/5 9:00
02/6/5 3:00
02/6/5 15:00
02/6/4 21:00
02/6/4 9:00
02/6/4 15:00
02/6/4 3:00
02/6/3 21:00
02/6/3 9:00
02/6/3 15:00
02/6/3 3:00
02/6/2 9:00
02/6/2 21:00
02/6/2 15:00
02/6/2 3:00
02/6/1 21:00
02/6/1 9:00
02/6/1 15:00
02/6/1 3:00
02/5/31 21:00
02/5/31 9:00
02/5/31 15:00
0.0
Fecha y Hora
Figura 5.1.3. Mediciones de niveles del mar en Pisco y el clculo del Modelo Global Topex para
el mismo lugar. El modelo sigue la prediccin de tabla, basada en mediciones del lugar, en
amplitud y fase, correctamente.
Debe aclararse que el modelo computa el nivel del mar en las vecindades del
punto donde regularmente se hacen las mediciones de mareas costeras. A
esto, en parte, se asocian las diferencias entre los resultados del modelo y la
prediccin.
29
5.2. Corrientes
Los modelos hidrodinmicos requieren de datos de corrientes para calibrarlos
y validarlos. Es decir, disponiendo de las series temporales de corrientes
(intensidad y direccin) en algunos puntos de inters obtenidas mediante
mediciones directas, pueden realizarse simulaciones numricas y comparar
los resultados obtenidos con las mediciones.
Propietario
Pluspetrol
Peru
Corporation
S.A.
Autor
H&O
Contenido
Descripcin de la corriente
medida en dos sitios. Los
datos son representados en
rosas de corrientes
Abril
de 2002
Pluspetrol
Peru
Corporation
S.A.
Sandwell
Engineering
Inc.
Estudio de Impacto
Ambiental - Oceanografa
(Lnea de Base)
2002
Dic. de
2002
Pluspetrol
Peru
Corporation
S.A.
Pluspetrol
Peru
Corporation
S.A.
Marzo
de 2003
Pluspetrol
Peru
Corporation
S.A.
30
Environmental
Resources
Management,
Per
Golder
Estudio de las corrientes
Associates Per
en el sitio de inters
S.A.
mediante el seguimiento de
flotadores. Presentan
tablas y grficos de las
trayectorias medidas.
Golder
Estudio de las corrientes
Associates Per
en el sitio de inters
S.A.
mediante el seguimiento de
flotadores. Presentan
tablas y grficos de las
trayectorias medidas.
Estudio complementario al
de diciembre de 2002.
Estacin E2:
Estacin E1:
Latitud:
13 46 12.10 S
Longitud:
76 14 44.66 W,
Latitud:
13 46 03.99 S
Longitud:
76 15 57.06 W.
Del anlisis de los datos, los autores sugieren que el ingreso de la corriente es
por ambos lados de la baha, principalmente durante las primeras horas del
da, y la salida se da por el fondo, en medio de la baha sobre la zona ms
ENVIRONMENTAL RESOURCES MANAGEMENT
31
32
Velocidad
(cm/s)
0-5
5-10
10-15
15-20
20-25
25-30
30-40
>40
TOTAL %
N
0,33
10,71
5,96
0,22
0,00
0,00
0,00
0,00
17,22
NE
0,33
18,10
11,37
1,21
0,22
0,00
0,00
0,00
31,24
E
0,11
7,17
2,43
0,00
0,00
0,00
0,00
0,00
9,71
Direccion
SE
S
0,55
0,33
2,54
1,99
0,55
4,86
0,00
4,30
0,00
0,11
0,00
0,00
0,00
0,00
0,00
0,00
3,64
11,59
SW
0,00
3,09
5,63
4,97
1,66
0,11
0,00
0,00
15,45
W
0,00
4,08
2,21
0,00
0,00
0,00
0,00
0,00
6,29
NW
0,11
3,31
1,43
0,00
0,00
0,00
0,00
0,00
4,86
Total
(%)
1,77
50,99
34,44
10,71
1,99
0,11
0,00
0,00
100,00
N
2,34
9,26
4,68
3,15
0,71
0,00
0,00
0,00
20,1
NE
2,24
13,33
7,22
2,75
0,81
0,00
0,10
0,00
26,4
E
2,03
4,27
1,42
0,31
0,10
0,00
0,00
0,00
8,1
Direccion
SE
S
1,63
1,73
3,87
13,53
0,81
0,71
0,00
0,00
0,00
0,00
0,00
0,00
0,00
0,00
0,00
0,00
6,3
16,0
SW
1,53
9,66
4,37
0,10
0,10
0,00
0,00
0,00
15,8
W
0,81
0,92
0,61
0,00
0,00
0,00
0,00
0,00
2,3
NW
1,02
2,64
1,22
0,00
0,00
0,00
0,00
0,00
4,9
Total
(%)
13,33
57,48
21,06
6,31
1,73
0,00
0,10
0,00
100,0
N
1,6
10,3
2,5
0,1
0,0
0,0
0,0
0,0
14,5
NE
0,2
18,0
13,5
0,5
0,1
0,0
0,0
0,0
32,3
E
0,0
5,1
5,5
0,1
0,0
0,0
0,0
0,0
10,7
Direccion
SE
S
0,2
1,4
5,1
5,9
4,6
3,0
0,6
0,4
0,0
0,0
0,0
0,0
0,0
0,0
0,0
0,0
10,5
10,7
33
SW
0,7
4,0
2,6
0,4
0,0
0,0
0,0
0,0
7,7
W
0,3
2,1
2,7
0,0
0,0
0,0
0,0
0,0
5,1
NW
1,7
5,6
1,1
0,0
0,0
0,0
0,0
0,0
8,4
Total
(%)
6,10
56,10
35,50
2,10
0,10
0,00
0,00
0,00
100
0.20
0.15
0.10
0.05
0.00
-0.05
-0.10
-0.15
-0.20
-0.25
-0.25
-0.20
-0.15
-0.10
-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.25
0.20
0.15
0.10
0.05
0.00
-0.05
-0.10
-0.15
-0.20
-0.20
-0.15
-0.10
-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
34
0.25
0.20
0.15
0.10
0.05
0.00
-0.05
-0.10
-0.15
-0.20
-0.20
-0.15
-0.10
-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.5
0.3
0.1
-0.1
-0.3
-0.5
-20
20
60
100
140
180
Figura 5.2.1.4. Curva de marea calculada por el MIKE 21 para el perodo del 16-22 de mayo 2002.
Puede observarse que la componente M2 (demidiurna de 12.4 h de perodo) no es importante y s
lo es la componente K1 (diurna de perdo 24.8 h). La separacin entre los picos ms altos
corresponde a la componente K1. Al final del registro la componente semidiurna comienza a
parecerse en magnitud a la diurna.
35
0.06
0.04
0.02
0.00
-0.02
-0.04
-0.06
-100
100
300
500
700
900
1100
Figura 5.2.1.5. Resultado del anlisis armnico de la serie de corrientes para el perodo 16-22
de mayo de 2002. Se muestra como nica componente significativa la onda K1 en la
componente N-S, en concordancia con la onda de marea. La separacin entre picos es de 24.8 h.
0.4
0.3
0.2
0.1
0.0
-0.1
-0.2
-0.3
-0.4
-20
20
60
100
140
180
Figura 5.2.1.6. Curva de marea calculada por el MIKE 21 para el perodo del 31 de mayo-7 de
junio de 2002. Puede observarse que la componente M2 (demidiurna de 12.4 h de perodo) es
importante como la componente K1 (diurna de perdo 24.8 h). La separacin entre los picos ms
altos corresponde a la componente K1.
36
0.04
0.02
0.00
-0.02
-0.04
-0.06
-0.08
0
42
84
126 168 210 252 294 336 378 420 462 504 546 588 630 672 714 756 798 840 882 924 966
Figura 5.2.1.7. Resultado del anlisis armnico de la componente E-W del registro de
corrientes del perodo 31 de mayo al 7 de junio de 2002. Puede observarse que las componentes
M2 y K1 son importantes del mismo modo que resultaron en la marea, para el mismo perodo.
0.03
0.02
0.01
0.00
-0.01
-0.02
-0.03
-0.04
0
42
84
126 168 210 252 294 336 378 420 462 504 546 588 630 672 714 756 798 840 882 924 966
Figura 5.2.1.8. Resultado del anlisis armnico de la componente N-S del registro de corrientes
del perodo 31 de mayo al 7 de junio de 2002. Puede observarse que las componentes M2 y K1
son importantes del mismo modo que resultaron en la marea, para el mismo perodo.
37
38
Tabla 5.2.2.1. Condiciones generales y caractersticas de las corridas de flotadores. Notas: N nmero de la
corrida, h.ini hora inicial, h.fin hora final.
N
DD/M
M/AA
h.ini
h. fin
BH
-4
06/03/2
003
8.30
15.00
BH
-6
06/03/2
003
8.30
20.30
BH
-8
12/11/2
002
5.00
17.00
BH
-10
Viento
(descrip.
gral.)
Marea
(descrip.
gral.)
Leves del
W
Semidiurna
Amplitud:
48 cm.
13/11/2
002
5.00
17.00
Leves del
W por la
maana.
Regulares
del SW por
la tarde.
Viento S
por la
maana.
Viento del
N por la
tarde.
A las
12.00
comienza a
soplar
fuerte
viento del
S. Se
levanta el
flotador.
Viento leve
del W por
la maana.
Regulares
del SW por
la tarde.
Calma por
la maana.
Ventoso
por la
tarde.
19/11/2
002
4.00
17.00
Viento del
S. Leve por
la maana,
fuertes por
la tarde.
Semidiurna.
Amplitud:
70 cm.
Primeras
horas
vientos
muy leves
del S.
Luego
leves del
W/SW.
Semidiurna.
Amplitud:
30 cm.
14/11/2
002
4.00
14.00
07/03/2
003
9.00
19.00
BH
-12
08/03/2
003
6.00
19.00
Dos corridas.
Se mueven
hacia NE.
Se varan en
costa.
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
Se mueve
lentamente
hacia el N.
NO
NO
NO
Se mueve
hacia el N.
Luego
cambio
brusco al E y
al N.
Se mueve
hacia el
NNE.
Se mueve
hacia el NW.
Cambio
brusco al
SW y N.
Se mueven
hacia el N.
Se mueve
lentamente
hacia la
costa.
Presenta
caracterstic
as mixtas.
Amplitud:
20 / 25 cm.
Se mueve
hacia el
N/NW. Al
final giro
abrupto al E
Presenta
caracterstic
as mixtas.
Amplitud:
60 cm.
Se mueve
hacia el W.
Luego hacia
el NNE.
Se mueve
hacia el N.
Cambia
abruptamente
hacia S.
Se mueve
hacia el W.
Luego hacia
el N.
Semidiurna.
Amplitud:
43 cm.
Se mueve
hacia el SE.
Luego hacia
el NNE.
Se mueve
hacia el
SSW. Luego
hacia el
NNE.
Se mueve
hacia el
SSW. Luego
hacia el ESE.
Presenta
caracterstic
as mixtas.
Amplitud:
64 / 25 cm.
Se mueve
hacia el
WNW. Giro
gradual
hacia SE.
Se mueve
hacia el
WNW.
Luego hacia
SSW y ESE.
Se mueve
hacia el
WNW. Giro
brusco hacia
el SSW.
NO
NO
NO
Se mueve
hacia el
NNW. Giro
hacia E y
ENE.
Se mueve
hacia el
NNE.
Se mueve
hacia el
NNE.
Semidiurna
Amplitud:
55 cm.
12
Se mueve
hacia el N.
Cambia
hacia el E.
Se mueve
hacia el W.
Luego hacia
el NE.
Se mueve
hacia el
NNE.
NOTA. Los casilleros con NO indican que no se realizaron corridas con veletas ubicadas a
esas profundidades. Se indican, adems, las amplitudes de marea correspondiente a los perodos
de medicin. Para los lapsos con desigualdades diurnas se indican la mayor y menor amplitud
registrada.
_______________________________________
39
DD/M
M/AA
h.ini
h. fin
18/11/2
002
4.00
16.00
Viento
(descrip.
gral.)
Vientos
leves del
N.
BH
-16
20/11/2
002
4.00
15.00
Predomina
ntes del S.
Marea
(descrip.
gral.)
12
NO
NO
NO
NO
NO
Se mueve
hacia el
NNE.
Se mueve
hacia el
NNE.
Se mueve
hacia el
NNE.
Se mueve
hacia el
NNE.
09/03/2
003
5.00
17.00
Calma por
la maana. Desigualdad
Leves del
es diurnas.
W por la
Amplitud:
tarde.
44 / 26 cm.
Se mueve
hacia el NW,
luego hacia
el N/NNE.
Se mueve
hacia el
NNE.
Se mueve
hacia el
NNE.
P200
10/03/2
003
4.00
17.00
Vientos
leves.
Predomina
ntes del W
/ SW.
Desigualdad
es diurnas.
Amplitud:
46 / 21 cm.
Se mueve
hacia el NW,
luego hacia
el NNE.
Se mueve
hacia el NW,
luego hacia
el NNE.
Se mueve
hacia el NE.
P500
10/03/2
003
5.00
17.00
Vientos
leves.
Predomina
ntes del W
/ SW.
Desigualdad
es diurnas.
Amplitud:
46 / 21 cm.
Se mueve
Se mueve
hacia el ESE, hacia el ESE,
luego hacia
luego hacia
el NE.
el NE.
Se mueve
Se mueve
Se mueve
Se mueve
Se mueve
hacia el NW, hacia el ESE, hacia el ESE, hacia el ESE, hacia el ESE,
luego hacia
luego hacia
luego hacia
luego hacia
luego hacia
el NNE.
el NE.
el NE.
el N.
el N.
NOTA. Los casilleros con NO indican que no se realizaron corridas con veletas ubicadas a
esas profundidades. Se indican, adems, las amplitudes de marea correspondiente a los perodos
de medicin. Para los lapsos con desigualdades diurnas se indican la mayor y menor amplitud
registrada.
_______________________________________
40
En Tambo de Mora las corrientes fluyen hacia el Sur con baja intensidad.
Por otro lado, las corrientes de marea se encontraban en una epoca de bajas
velocidades debido a la combinacin de fases, desfavorable para la intensidad
de las corrientes, entre las componentes astronmicas M2 (semidiurna) y K1
(diurna).
Figura 5.2.2.1. Corridas de flotadores durante marzo de 2003 con el objeto de obtener
informacin de corrientes para la calibracin del MIKE 21. La figura cubre el dominio de
modelacin.
ENVIRONMENTAL RESOURCES MANAGEMENT
41
Figura 5.2.2.2. Corridas de flotadores durante marzo de 2003 para la zona de la Transecta A.
Figura 5.2.2.3. Corridas de flotadores durante marzo de 2003 para la zona de la Tambo de
Mora
Figura 5.2.2.4. Corridas de flotadores durante marzo de 2003 para la zona de la Transecta B.
42
-0.05
-0.10
-0.15
-0.20
-0.25
-0.30
-0.35
-0.40
-0.45
-0.50
-0.50
-0.45
-0.40
-0.35
-0.30
-0.25
-0.20
-0.15
-0.10
-0.05
0.00
0.05
0.10
0.00
-0.05
-0.10
-0.15
-0.20
-0.25
-0.30
-0.35
-0.40
-0.40
-0.35
-0.30
-0.25
-0.20
-0.15
-0.10
-0.05
0.00
43
0.15
0.10
0.05
0.00
-0.05
-0.10
-0.15
-0.15
-0.10
-0.05
0.00
0.05
0.10
0.15
0.20
Como estas mediciones son las nicas disponibles costa afuera, sern
utilizados para calibrar y validar las simulaciones realizadas con el MIKE 21,
dada la necesidad que el modelo tiene de ellas.
-0.04
-0.25
0.08
-0.43
-0.22
-0.09
Vientos del S
Calma y vientos
Calma y vientos
del SW
del W
Las corridas fueron suspendidas antes de las 13 h de medicin
por condiciones adversas de tiempo y oleaje
44
45
El modelo de alta resolucin OCCAM est basado en el modelo de Bryan-CoxSemtner pero, a diferencia de ste, no tiene tapa rgida y, por lo tanto,
permite intercambios atmosfricos de masa y de cantidad de movimiento.
Las salidas del modelo (componentes Norte y Este de la corriente para todos
los niveles correspondientes a cada nodo del dominio) se encuentran
disponibles en Internet con un discretizacin temporal de 15 das. Para
mayores detalles sobre el modelo y sobre cmo descargar datos, puede
accederse a:
http://www.soc.soton.ac.uk/jrd/occam/welcome.html
46
Sobre
la
base
de
las
publicaciones
informes
tcnicos
citados
47
5.3. Vientos
Informacin histrica de vientos en Pisco fue introducida anteriormente.
Durante el perodo de mediciones de corrientes con flotadores, se midieron
vientos en Pisco, cuyas caractersticas se incluyen en las Figuras 5.3.1. a 5.3.2.
Los datos muestran que los vientos medidos tienen periodicidad del orden de
las 24 h. Es decir, repiten un comportamiento diurno, al que eventualmente se
le suman valores medios, como se observa en las figuras mencionadas
anteriormente.
7
6
5
4
3
2
1
0
-1
03/3/14 0:00
03/3/14 12:00
03/3/13 0:00
03/3/13 12:00
03/3/12 12:00
03/3/12 0:00
03/3/11 0:00
03/3/11 12:00
03/3/10 0:00
03/3/10 12:00
03/3/9 0:00
03/3/9 12:00
03/3/8 0:00
03/3/8 12:00
03/3/7 0:00
03/3/7 12:00
03/3/6 0:00
03/3/6 12:00
-2
Fecha y Hora
Figura 5.3.1. Mediciones de viento en Pisco durante las corridas de flotadores. Componente EW. La componente tiene una marcada periodicidad del orden de las 24 h. La velocidad media de
esta componente durante el perodo de medicin fue de 2.3 m/s.
48
4
3
2
1
0
-1
-2
-3
-4
-5
03/3/14 0:00
03/3/14 12:00
03/3/13 12:00
03/3/13 0:00
03/3/12 0:00
03/3/12 12:00
03/3/11 0:00
03/3/11 12:00
03/3/10 0:00
03/3/10 12:00
03/3/9 0:00
03/3/9 12:00
03/3/8 0:00
03/3/8 12:00
03/3/7 0:00
03/3/7 12:00
03/3/6 0:00
03/3/6 12:00
-6
Fecha y Hora
Figura 5.3.2. Mediciones de viento en Pisco durante las corridas de flotadores. Componente NS. La componente tiene una marcada periodicidad del orden de las 24 h. La velocidad media de
esta componente durante el perodo de medicin fue de 1.8 m/s.
700
600
600
500
500
400
400
300
300
200
200
100
100
0
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
0.45
0
0.50
Frequencia (1/hora)
Figura 5.3.3. Mediciones de viento en Pisco durante las corridas de flotadores. Representacin
espectral de la componente E-W.
Vientos en Pisco. Representacion espectral de la componente N-S del viento.
700
600
600
500
500
400
400
300
300
200
200
100
100
0
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
0.45
0
0.50
Frequencia (1/hora)
Figura 5.3.4. Mediciones de viento en Pisco durante las corridas de flotadores. Representacin
espectral de la componente N-S.
49
7
6
5
4
3
2
1
0
45
90
135
180
225
270
315
360
Direccin
Figura 5.3.5. Diagrama de dispersin de las componentes E-W y N-S del viento.
5.4. Batimetra
Los datos de profundidad fueron extrados de las cartas nuticas Per
Hidronav Nro. 226 y 2262. Las profundidades de las cartas estn referidas al
nivel medio de bajamares de sicigias ordinaria.
50
6. EL MODELO NUMRICO
6.1. Generalidades
El modelo hidrodinmico MIKE 21, desarrollado por el Danish Hydraulic
Institute, es un sistema de modelado para flujos bidimensionales (2D) a
superficie
libre.
El
MIKE
21
es
aplicable
para
simulaciones
El MIKE 21 puede ser utilizado para simular una amplia gama de situaciones
las cuales incluyen condiciones asociadas a flujos de marea, ondas de
tormenta y calidad de agua, entre otras. En el APNDICE A se brinda una
breve descripcin del modelo hidrodinmico y en el APNDICE B algunas
aplicaciones del modelo y en el APNDICE C la licencia del DHI para su uso.
plantearn
escenarios
caracterizados
por
forzantes
mareolgicos,
51
implementado
un
modelo
numrico
hidrodinmico
Del anlisis comparatorio entre los datos registrados en Pisco y los resultantes
de la modelacin surgen diferencias del nivel de mar comprendidos entre 2 y
15 cm. Asociadas probablemente al efecto meteorolgico. Los resultados
indican corrientes internas entre la pennsula y la isla San Galln. Algunos
ensayos realizados incorporando el efecto del viento como forzante indican
que el sistema responde claramente al efecto atmosfrico.
52
53
estudiado que permitan elaborar una condicin inicial realista para un modelo
3D.
Los puntos de la grilla adoptada son 114 en X (E-W) y 147 en Y (N-S). Dado
que las celdas de la grilla son cuadradas y de 500 m por 500m, la longitud del
dominio de modelacin en X es de 57 km y en Y 73.5 km.
54
simulaciones
fue
entre
2880
4320.
Consecuentemente,
55
6.8. Batimetra
Con
los
datos
de
profundidad
obtenidos
de
las
cartas
indicadas
6.9. Forzantes
6.9.1. Bordes
56
Contorno superior del dominio, Y = 146 (lmite Norte del rea estudiada,
aproximadamente en la latitud 13 22 S).
6.9.2. Mareas
57
6.9.3. Corrientes
6.9.4. Vientos
Los datos de viento disponibles, estadsticos y en forma de series de tiempo,
sern aplicados a la modelacin matemtica, por medio de la configuracin de
escenarios realistas. Los valores de viento fueron introducidos anteriormente
en este informe.
58
6.11. Resistencia
El coeficiente de friccin de fondo fue utilizado para calibrar el modelo. En
este trabajo se utiliza un coeficiente contante para toda el rea modelada. Se
realizaron
algunas
simulaciones
preliminares
utilizando
diferentes
59
60
Sin embargo, es muy poco lo que se sabe sobre las corrientes fuera de las
bahas, muy necesarias para la calibracin del modelo, dado que el dominio de
modelacin se extiende varias decenas de kilmentros costa afuera, como fue
indicado.
Como resultado de esto, se dispuso medir costa afuera como fue mencionado
anteriormente. Las mediciones, en profundidades entre 50 m y 60 m,
indicaron velocidades de corrientes con componentes N-S importantes a lo
largo de unas 12 h, aun en condiciones de viento en contra. Por otro lado, las
corrientes de marea se encontraban en una epoca de bajas velocidades debido
a la combinacin de fases, desfavorable para la intensidad de las corrientes,
entre las componentes astronmicas M2 (semidiurna) y K1 (diurna).
El prximo paso fue establecer las condiciones de contorno en los lmites del
dominio de modelacin para esta contracorriente. Por un lado, el modelo
OCCAM dispone de valores de corrientes variables en la profundidad
(modelo multicapa) solamente en 3 puntos del dominio de modelacin
ubicados sobre el borde oeste del mismo, por lo que no se dispuso de una
distribucin de estas corrientes hasta la costa, tanto al sur como al norte del
dominio. Se decidi, en consecuencia, utilizar las mediciones de corrientes en
Tambo de Mora (cercanas al borde norte del dominio) para ajustar los valores
del mdulo de la velocidad provista por el OCCAM, respetando las
direcciones del mismo.
61
Con respecto a las mareas, las variaciones del nivel del mar modeladas por el
MIKE 21 coinciden con las amplitudes y fases medidas y predichas por la
Tabla de Marea para Pisco. Las bondades del TOPEX para calcular la marea de
El Callao y Pisco fueron presentadas anteriormente.
De los resultados obtenidos surge que los niveles del mar son simulados
satisfactoriamente, logrndose un alto grado de ajuste entre los valores de
tablas y los calculados para las vecindades de en Pisco por medio de las
simulaciones numricas realizadas.
Comparacin entre los niveles del mar predichos por la Tabla de Mareas y
los calculados por el MIKE21 en base a los resultados del modelo global Topex
1.2
Niveles de Tabla
1.0
Niveles de MIKE21
0.8
0.6
0.4
0.2
0.0
02/5/16 0:00
02/5/16 4:01
02/5/16 7:58
02/5/16 12:00
02/5/16 16:01
02/5/16 19:58
02/5/17 0:00
02/5/17 4:01
02/5/17 7:58
02/5/17 12:00
02/5/17 16:01
02/5/17 19:58
02/5/18 0:00
02/5/18 4:01
02/5/18 7:58
02/5/18 12:00
02/5/18 16:01
02/5/18 19:58
02/5/19 0:00
02/5/19 4:01
02/5/19 7:58
02/5/19 12:00
02/5/19 16:01
02/5/19 19:58
02/5/20 0:00
02/5/20 4:01
02/5/20 7:58
02/5/20 12:00
02/5/20 16:01
02/5/20 19:58
02/5/21 0:00
02/5/21 4:01
02/5/21 7:58
02/5/21 12:00
02/5/21 16:01
02/5/21 19:58
02/5/22 0:00
02/5/22 4:01
02/5/22 7:58
02/5/22 12:00
02/5/22 16:01
02/5/22 19:58
-0.2
Fecha y Hora
Figura 6.12.1. Comparacin entre los valores del nivel del mar, predichos por la Tabla de
Mareas (solamente pleamares y bajamares), y los correspondientes calculados por el MIKE 21
alimentado por los resultados del modelo global TOPEX.
62
Comparacin entre el nivel del mar calculado por el Modelo Global de Mareas
Topex y el correspondiente calculado por el MIKE 21, cerca de Pisco
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.0
-0.1
-0.2
-0.3
-0.4
-0.5
MIKE 21
TOPEX
-0.6
03/3/28 1:58
03/3/28 18:00
03/3/27 10:01
03/3/26 1:58
03/3/26 18:00
03/3/25 10:01
03/3/24 1:58
03/3/24 18:00
03/3/23 10:01
03/3/22 1:58
03/3/22 18:00
03/3/21 10:01
03/3/20 1:58
03/3/20 18:00
03/3/19 10:01
03/3/18 1:58
03/3/18 18:00
03/3/17 10:01
03/3/16 1:58
03/3/16 18:00
03/3/15 10:01
03/3/14 1:58
03/3/14 18:00
03/3/13 10:01
03/3/12 1:58
03/3/12 18:00
03/3/11 10:01
-0.7
Fecha y Hora
Figura 6.12.2. Comparacin entre los valores horarios del nivel del mar calculados por el
TOPEX y los correspondientes calculados por el MIKE 21 en las vecindades de Pisco,
alimentado, en sus bordes, por los resultados del modelo global TOPEX.
63
Comparacin entre el nivel del mar calculado por el Modelo Global de Mareas
Topex y el correspondiente calculado por el MIKE 21, costa afuera
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.0
-0.1
-0.2
-0.3
-0.4
MIKE 21
TOPEX
-0.5
-0.6
03/3/28 1:58
03/3/28 18:00
03/3/27 10:01
03/3/26 1:58
03/3/26 18:00
03/3/25 10:01
03/3/24 1:58
03/3/24 18:00
03/3/23 10:01
03/3/22 1:58
03/3/22 18:00
03/3/21 10:01
03/3/20 1:58
03/3/20 18:00
03/3/19 10:01
03/3/18 1:58
03/3/18 18:00
03/3/17 10:01
03/3/16 1:58
03/3/16 18:00
03/3/15 10:01
03/3/14 1:58
03/3/14 18:00
03/3/13 10:01
03/3/12 1:58
03/3/12 18:00
03/3/11 10:01
-0.7
Fecha y Hora
Figura 6.12.3. Comparacin entre los valores horarios del nivel del mar calculados por el
TOPEX y los correspondientes calculados por el MIKE 21 costa afuera, alimentado, en sus
bordes, por los resultados del modelo global TOPEX
64
-0.3
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1.0
Figura 6.12.4. Distribucin de velocidades en la vertical por medio del anlisis de las corridas
de flotadores. Nota: z es la distancia de la superficie del mar hasta la posicin de la pala de
arrastre y h es la profundidad local
Tabla 6.12.1. Comparacin de los valores medios de corrientes medidos por los flotadores y
los correspondientes calculados por el MIKE 21 para Transecta A, Transecta B y Tambo de
Mora
Componentes
Transecta A
Transecta B
Tambo de Mora
de la
Velocidad
MK21
Esperado
MK21
Esperado
MK21
Esperado
E-W (m/s)
-0.05
-0.04
-0.16
-0.22
0.01
N-S (m/s)
-0.39
-0.38
-0.20
-0.17
-0.11
-0.08
Nota: (*) la calibracin se realiz con la componente N-S de la Contracorriente del Per, sin
vientos. Durante la corrida se produjo viento del W que hizo derivar al flotador hacia la costa.
65
Tabla 6.12.2. Comparacin de las corridas de flotadores para el mes de marzo de 2003 con los
resultados del MIKE 21.
Corrida
Comentarios
Observacin (cm/s)
MIKE 21 (cm/s)
BH-4
Simulacin
6
6
cualitativamente bien.
Subestima la velocidad.
BH-6
Simula bien.
5-7
8
BH-10
No reproduce el
7-9
6
movimiento inicial hacia
el sur.
El desplazamiento hacia
el N/NNE es bien
simulado.
BH-12
Simula bien.
5-6
7
P-200
Salvo la trayectoria en
3
6
superficie (z=0/3m), el
resto es bien simulado.
Subestima la velocidad.
P-500
Salvo la trayectoria, el
3
6
resto es bien simulada.
Subestima la velocidad.
66
Direccin
Ocurrencia
para todos de la
Velocidad (m/s)
direccin
considerada
4
2.1
NE
1.1
2.0
SE
1.9
SW
SW
SW
12
10
11
10
12
NW
13
NW
13.9
28.0
22.3
7.8
14
15
16
MAREA SOLAMENTE
67
7.2. Resultados
De acuerdo a las mediciones disponibles de vientos, mareas y corrientes,
utilizadas para calibrar el MIKE 21 y los resultados suministrados por el
propio modelo, las corrientes en las Bahas de Pisco y Paracas son producto de
los siguientes factores:
un
comportamiento
peridico
(marzo
2003,
presentados
enteriormente) alrededor de un valor medio en ambas componentes E-W y NS. Las corrientes muestran un comportamiento similar, alcanzando un
mximo de velocidad de 0.13 m/s en la componente N-S, con velocidades
mximas del viento de 9 m/s, esto significa que se obtienen valores de
corriente (promediada en la profundidad) de 1.5% del valor de la velocidad
del viento. La componente E-W es de escasa magnitud dada la presencia de la
costa. La Figura 7.2.1. muestra los resultados de la modelacin para un punto
sobre la traza del proyecto de Terminal Marima.
68
0.08
0.04
0.00
-0.04
-0.08
03/3/10 19:58
03/3/10 16:01
03/3/10 7:58
03/3/10 12:00
03/3/10 4:01
03/3/10 0:00
03/3/9 19:58
03/3/9 16:01
03/3/9 7:58
03/3/9 12:00
03/3/9 4:01
03/3/9 0:00
03/3/8 19:58
03/3/8 16:01
03/3/8 7:58
03/3/8 12:00
03/3/8 4:01
03/3/8 0:00
03/3/7 19:58
03/3/7 16:01
03/3/7 7:58
03/3/7 12:00
03/3/7 4:01
-0.12
03/3/7 0:00
0.12
Fecha y Hora
Figura 7.2.1. Corrientes generadas por un viento variable en el tiempo, para un punto ubicado
sobre la traza proyectada para la Terminal Martima, entre la costa y el extremo oeste de la
obra. El valor medio de la componente E-W es menor a 0.01 m/s, mientras que el
correspondiente a la componente N-S es de 0.04 m/s.
69
70
71
72
73
74
75
76
77
78
79
80
81
82
peridica diurna
83
84
85
Tabla 7.1.2. Resumen con los resultados de los escenarios de modelacin seleccionados
Viento
Escenario
Direccin
del
Velocidad
(m/s)
Ocurrencia %
para todos de la
direccin
considerada
2.1
NE
1.1
2.0
SE
1.9
SW
SW
13.9
28.0
9
SW
12
10
11
10
12
NW
13
NW
14
22.3
7.8
86
Comentarios
Nota: B.P. Baha de Paracas, P.P. Pennsula
de Paracas. Sitio: extremo W de la Terminal
Martima
Corrientes hacia el S disminuyendo hacia la
B.P. El agua entra por la margen E de la
baha y sale bordeando la P.P. por el W.
Dentro de la B.P. se forma una celda de
circulacin debido a su batimetra.
Velocidades en el Sitio de unos 0.1 m/s y la
mitad de ellas dentro de la B.P.
Patrn de circulacin semejante al anterior.
La aguas empujadas hacia la costa se alinean
con esta en direccin a la B.P. Velocidades
en el Sitio de unos 0.1 m/s y la mitad de ellas
dentro de la B.P.
Patrn de circulacin semejante al anterior
con velocidades ms bajas. La aguas
empujadas hacia la costa se alinean con esta
en direccin a la B.P. Velocidades en el Sitio
de unos 0.05 m/s y semejantes dentro de la
B.P.
Corrientes hacia el N. dentro de las bahas.
El agua entra por la margen W de la baha y
sale bordeando la P.P. por el E. Dentro de la
B.P. se forman celdas de circulacin debido a
su batimetra. Velocidades en el Sitio
menores de 0.1 m/s y menores de0.05 m/s
dentro de la B.P.
Con velocidades del viento de 4m/s, las
corrientes van hacia el N con velocidades en
el Sitio menores de 0.1 m/s y menores de
0.05 m/s dentro de la B.P. Con velocidades
del viento de 8m/s, las corrientes duplican su
intensidad.
Con velocidades del viento de 4m/s, las
corrientes van hacia el N con velocidades en
el Sitio menores de 0.1 m/s y menores de
0.05 m/s dentro de la B.P. Con velocidades
del viento de 8m/s, las corrientes duplican su
intensidad. Con velocidades del viento de
12m/s, las corrientes exceden 0.25 m/s fuera
de la B.P. y son menores de 0.2 m/s dentro
de la misma.
Corrientes hacia el S, dbiles. Entran por el
lado W de la B.P. y el balance se produce por
la salida hacia en N por el lado E.
Corrientes hacia el S. Entran por el lado E de
la B.P. y el balance se produce por la salida
hacia en N por el lado W.
Las partculas se desplazan y se separan del
mismo modo que lo hacen los flotadores.
Siguen la direccin del viento a velocidad
media disminuyendo y aumentando su
velocidad segn la periodicidad del viento.
En las vecindades del borde E de la P.P.
algunas partculas entran a la baha y otras
se alejan hacia el N, an habiendo sido
lanzadas a corta distancia entre s.
Tabla 7.1.2. Continuacin. Resumen con los resultados de los escenarios de modelacin
seleccionados
Escenario
Viento
15
16
MAREA SOLAMENTE
Comentarios
Nota: B.P. Baha de Paracas, P.P. Pennsula
de Paracas. Sitio: extremo W de la Terminal
Martima
La contracorriente, importante fuera de las
bahas, no influye significativamente dentro
de ellas, an sin vientos ni marea.
Durante las bajantes las corrientes vacian la
B.P. y lo contrario ocurre durante las
crecientes. Las velocidades son muy bajas
dependiendo de la condicin de marea de la
poca
87
CONCLUSIONES
Dada la complejidad que presenta la Baha de Paracas desde el punto de vista
ambiental y la necesidad de conocer mejor la circulacin para tratar aspectos
vinculados a los sedimentos en suspensin y su desplazamiento y eventuales
derrames, se ha revisado la informacin disponible relevante a la circulacin,
se ha calibrado un modelo matemtico y se lo ha utilizado para el clculo del
campo de corrientes para toda la regin de influencia del proyecto.
Sobre los Datos Disponibles y su Anlisis
88
Los datos muestran que los vientos medidos tienen periodicidad del
orden de las 24 h. Es decir, repiten un comportamiento diurno, al que
eventualmente se le suman valores medios. Durante las mediciones, de
corta duracin, se encontraron vientos medios del SW con velocidades
medias de unos 3 m/s y oscilaciones superpuestas con hasta unos 6 m/s.
El anlisis espectral de los datos muestra un pico de energa significativo
en una frecuencia algo menor a 0.05 1/h (algo mayor a 20 h de perodo).
89
Las velocidades obtenidas por medio de las corridas costa afuera, entre la
Isla San Guilln y la Pennsula de Paracas, alcanzaron valores medios,
durante el perodo de medicin, de unos 0.45 m/s, siempre hacia el S, a
pesar de soplar viento en contra del S. Situaciones similares, aunque con
velocidades menores, se registraron en los otros dos lugares. Durante la
poca, las corrientes de marea se localizaban en un mnimo, por lo que no
se encontr otra explicacin a las corrientes medidas que no fuera la
Contracorriente del Per.
90
La comparacin entre los niveles del mar calculados por el modelo MIKE
21 -utilizado para la modelacin- fue inmediatamente satisfactoria en
todos los casos, con diferencias tpicas de la modelacin matemtica, por
lo que se dio por calibrado respecto al desplazamiento de la superficie del
mar.
91
para
entender
completamente
su
funcionamiento,
es
por
escenarios
Con vientos del W las corrientes se dirigen hacia el S, pero muy dbiles y
modificables por las circunstancias de la marea. Entran por el lado W de
la Baha de Paracas y el balance se produce por la salida hacia en N por el
lado E de la misma.
92
Los vientos del S, SE y SW, con una frecuencia de ocurrencia del 44%
aproximadamente, generan corrientes hacia el N. Las corrientes
generadas por los vientos del SW resultaron las ms fuertes y del orden
de los 0.25 m/s para velocidades de viento de 12 m/s.
Con vientos del S, las corrientes van hacia el N con velocidades en el sitio
menores de 0.2 m/s y menores de 0.1 m/s dentro de la Baha de Paracas.
93