Modelacion de Campos

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 93

INFORME

Pluspetrol Per Corporation


Modelacin Matemtica del
Campo de Corrientes en las
Bahas de Pisco y Paracas,
Per.
Abril 2003
Ref. PLU_03_510

ERM Per S.A.


Environmental Resources Management
Vctor A. Belande 147, Va Principal 140 Ed. Real 6,
P6
Lima 27 Per
Telfono 51-1-211-2677
Fax 51-1-211-2526
Web Site: www.erm.com

INFORME

Pluspetrol Per Corporation


Modelacin Matemtica del
Campo de Corrientes en las
Bahas de Pisco y Paracas,
Per.
Abril 2003
Ref. PLU_03_510
En nombre y por cuenta de ERM Per
Aprobado por: Alberto Sambartolome
Firma: ______________________________________
Cargo:

_____________________________________

Fecha: ______________________________________

Este documento ha sido elaborado por ERM Per con la debida


competencia, diligencia y cuidado con arreglo a los trminos del
contrato estipulado con el Cliente y nuestras condiciones generales
de suministro, utilizando los recursos concertados.
ERM Per declina toda responsabilidad ante el cliente o terceros por
cualquier cuestin que no est relacionada con lo anteriormente
expuesto.
Este documento tiene carcter reservado para el Cliente. ERM Per
no asume ninguna responsabilidad ante terceros que lleguen a
conocer este informe o parte de l.

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

ENVIRONMENTAL RESOURCES MANAGEMENT

1. INTRODUCCIN

Durante los estudios de evaluacin de impactos ambientales para la


construccin de la Terminal Martima en Playa Lobera, han sido
identificados algunos potenciales problemas que, su tratamiento y
evaluacin dependen, en parte, del conocimiento de la circulacin
costera en el rea.

Uno de estos problemas, lo constituye la puesta en suspensin de


sedimentos finos depositados en el fondo por acciones tpicas de la
construccin, tales como el armado de trincheras para la instalacin de
ductos, las corrientes inducidas por hlices de remolcadores y buques,
arrastres de anclas, movimientos de tuberas sobre el fondo y otros.
Otro de los problemas identificados lo constituye un potencial derrame
de algn producto establecido en los ductos, para su exportacin.

En ambos casos, la evolucin de las "manchas" de sedimentos o


producto, sus concentraciones, superficies ocupadas, velocidades de
traslacin y destinos geogrficos, constituyen aspectos de importancia
a lo largo de las evaluaciones y acciones ambientales relativas a
proyectos como la terminal martima en cuestin.

El seguimiento y prediccin de estos fenmenos dependen muy


significativamente de las corrientes (circulacin) dentro de la zona de
influencia de la terminal y la baha.

Las mediciones disponibles son de corta duracin, puntuales y dentro


de la bahia, por lo que se desconoce la circulacin general dentro de la
baha y vecindades. En estas condiciones, es difcil la interpretacin
fsica de las fuerzas que generan las corrientes en un ambiente con
notorios accidentes geograficos y batimetra variada.

En general, y en base a la informacin disponible, las corrientes en el


rea son la resultante de :
ENVIRONMENTAL RESOURCES MANAGEMENT

Corrientes puras de marea, caracterizadas por tener perodos


cclicos, similares a los perodos de las componentes armnicas
dominantes para la marea del lugar.

Corrientes inducidas por el viento. La friccin viscosa del viento


sobre la superficie del mar induce una fuerza de arrastre que
genera una corriente, no necesariamente alineada con la direccin
del viento.

Corrientes permanentes debidas a patrones de circulacin ocenica


general. Estas corrientes influyen la circulacin dentro de la baha,
debido a su cercana a aguas profundas.

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, el presente informe incluye
una revisin de la informacin disponible relevante a la circulacin y el
clculo del campo de corrientes para toda la regin de influencia del
proyecto.

Se han analizado las corrientes, considerando su

variabilidad en el espacio y el tiempo, para escenarios frecuentes en el


lugar

provientes

de

la

combinacin

de

las

fuerzas

citadas

anteriormente. Esta tarea no es posible hacerla de modo artesanal, por


lo que se ha elegido un modelo matemtico.

Para estudiar el campo de corrientes en la zona de Baha Pisco se


implement el modelo hidrodinmico (HD) bidimensional (2D) MIKE
21 del DHI (Danish Hydraulic Institute de Dinamarca) que consiste en
una herramienta matemtica, esto es un modelo matemtico, para el
calculo de velocidades y direcciones de la corriente y desplazamientos
del nivel del mar, debidos a fuerzas tales como mareas y vientos.

Este modelo es aceptado en aplicaciones de ingenieria y medio


ambiente (evaluaciones ambientales, prediccion y monitoreo) y muy
generalizado su uso en todo el mundo, incluyendo aplicaciones a obras
del tipo de la terminal martima en la Bahia de Pisco.

ENVIRONMENTAL RESOURCES MANAGEMENT

El dominio considerado para la modelacin numrica se extiende


aproximadamente desde 13 22 S hasta 14 02 S, en la direccin N-S y
desde 76 40 W hasta la costa, en la direccin W-E. Se adopt una
resolucin espacial de 500 m por 500 m con la cual quedaron
convenientemente representados los accidentes costeros y la batimetra
de la zona en estudio. La calibracin y validacin del modelo se
realizaron utilizando datos disponibles de marea y corrientes.

Con el objeto de estudiar las corrientes en el sitio de inters se


disearon algunos escenarios tpicos considerando la importancia de
los forzantes de las corrientes. De esta manera se obtuvieron y
analizaron cartas de corrientes para la zona de estudio. Se incluyen
grficos y tablas para facilitar la interpretacin de los resultados
obtenidos en este trabajo.

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.

El modelo queda en condiciones de ser revalidado en caso de


disponerse de ms mediciones de corrientes, especialmente.

ENVIRONMENTAL RESOURCES MANAGEMENT

2. OBJETIVO

El objetivo principal del presente trabajo es determinar la circulacin


dentro de la baha de Pisco, rea de inters de este proyecto, para cada
instante, y para situaciones representativas (condiciones de viento y
marea, principalmente), usando las herramientas de la simulacin
matemtica, calibrados mediante las mediciones disponibles a la fecha.

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

Los resultados pueden ser aplicados a distintos aspectos ambientales


del proyecto tales como:

Evaluaciones de impactos ambientales ms slidas y


convincentes

Optimizacin de las medidas de mitigacin

Optimizacin del plan de monitoreo

Prediccin y seguimiento de las consecuencias ambientales


de la obra

ENVIRONMENTAL RESOURCES MANAGEMENT

4. DESCRIPCIN GENERAL DE LA ZONA DE ESTUDIO

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.

La zona comprendida entre Caucato y la desembocadura del ro Pisco,


y el fondo de la Baha de Paracas se encuentran casi al nivel del mar,
elevndose entre el entrante de Sequin y Punta Pejerrey, donde se
inicia la zona de acantilados. El relieve del fondo submarino frente al
borde costero de Pisco, es de suave pendiente, hacindose ms
marcado frente a la desembocadura del ro Pisco.

Plano 4.1.1. Bahias de Pisco, Paracas y Sistema de Islas.

4.2. Condiciones Oceanogrficas Generales


4.2.1. Temperatura del Agua de Mar

Las variaciones de la temperatura superficial del mar estn


relacionadas fundamentalmente con la absorcin de la radiacin sobre
la circulacin atmosfrica, las corrientes superficiales, los afloramientos
costeros,

los

hundimientos

ENVIRONMENTAL RESOURCES MANAGEMENT

movimientos
9

advectivos.

Las

temperaturas promedio mensuales muestran, a todo lo largo de la


costa peruana, un comportamiento cclico anual, con valores
relativamente altos durante los meses de verano y bajos en invierno
(EIA, Lnea de Base Ambiental, Oceanografa, ERM, 2002)

Los valores promedios de las transectas 14, 13 y 12 son representativos


de la fra corriente que baa la costa del Per, aunque en la zona de
Pisco, que incluye la baha de Paracas, debido a la proteccin de la
pennsula del mismo nombre de los sistemas de circulacin ocenica
que proceden del sur, y la formacin de Eddys, ocasionan que las
temperaturas en el rea mencionada, sean mayores a las zonas
adyacentes, tanto hacia el sur como hacia el norte. Su variabilidad
estacional es relativamente baja, y slo se altera con la ocurrencia del
fenmeno de "El Nio", en que aparecen aguas clidas.

La variacin de la temperatura con respecto a la profundidad, mostr


una capa de mezcla localizada en los 2 primeros metros de
profundidad. En las estaciones de poca profundidad, se observ un
gradiente negativo desde la superficie, mostrando una dbil
termoclina.

Esta situacin nos muestra la existencia de 2 capas de agua de


diferentes densidad, comportndose como una barrera de distribucin
que delimita 2 masas de agua de diferentes caractersticas fsicas y
dinmicas, lo que permiti la seleccin de las profundidades de
medicin de corrientes. En tal sentido, las corrientes superficiales se
tomaron por encima de los 2 metros de profundidad, mientras que las
de fondo o subsuperficiales, se registraron por debajo de los 2 metros
de profundidad y a 1 metro sobre el fondo (EIA, Lnea de Base
Ambiental, Oceanografa, ERM, 2002)

En la Tabla 4.2.1.1. se muestran los resultados de las mediciones de


temperatura del agua de mar por transectos distribuidos a lo largo de
las Bahias de Pisco y Paracas y para diferentes profundidades. La
ubicacin de los transectos se indica en la Figura 4.2.1.1.

ENVIRONMENTAL RESOURCES MANAGEMENT

10

Figura 4.2.1.1. Bahias de Pisco, Paracas y ubicacin de las transectas.

ENVIRONMENTAL RESOURCES MANAGEMENT

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

Valores promedio por


transecto
21.96

23

23.3

24.18

23

25.04

24.64

26.04

26.24

24.18

26.7

19.4

19.8

17.64

El valor mnimo registrado fue de 15.5 C en el transecto 13 (cota de 10


m) a 9.2 m de profundidad.
ENVIRONMENTAL RESOURCES MANAGEMENT

12

El valor mximo registrado fue de 29.5 C en los transectos 8 y 9 (cota


de o m), en superficie. El promedio de temperaturas resgistrado a lo
largo de todos los transectos fue de 23.2 C.

Respecto a la distribucin de la temperatura promedio por transectos, se pudo


observar un incremento gradual desde el transecto 1 (21.96 C) hasta un
mximo en el transecto 9 (26.24 C), donde luego se observa una distribucin
aleatoria, hasta llegar al mnimo registrado en el transecto 14 (17.64 C).

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

4.2.3. Oxgeno Disuelto (OD)

El oxgeno disuelto a nivel superficial, present concentraciones entre


5,41 y 4,25 mg/l con un promedio de 4,73 mg/l. La estructura vertical
del oxgeno, de acuerdo a los valores encontrados en superficie y cerca
del fondo, muestran un comportamiento similar a la estructura
trmica, es decir con un gradiente negativo desde la superficie. Los
valores cerca del fondo fluctuaron entre 5,13 a 1,19 mg/l, con un
ENVIRONMENTAL RESOURCES MANAGEMENT

13

promedio de 2,49 mg/l. (EIA, Lnea de Base Ambiental, Oceanografa,


ERM, 2002)

En cuanto a la distribucin promedio de los niveles de OD a lo largo de


los transectos, hubo una mxima concentracin en el transecto 11, con
un valor de 9.38 mg/l y una mnima de 4.7 mg/l en el transecto 14. El
valor medio de concentracin fue de 7.2 mg/l. .La Tabla 4.2.3.1.
muestra los niveles medidos de OD.

4.2.4. Slidos en Suspensin, Turbidez y Transparencia

Los slidos en suspensin se distribuyeron con concentraciones en


superficie, entre 7,5 y 2 mg/l, con un promedio de 5,5 mg/l, y en el
fondo, con valores entre 11 y 4 mg/l, con un valor promedio de 6,9
mg/l. Las concentraciones encontradas muestran valores por debajo
de los lmites permisibles.

Concentraciones mayores de 25 mg/l, pueden producir asfixia,


principalmente a las poblaciones bentnicas que viven adheridas al
fondo o tienen poca movilidad. (EIA, Lnea de Base Ambiental,
Oceanografa, ERM, 2002)

En el documento se incluyen mediciones de turbidez entre la zona del


Ro Pisco (Transecto 1), a lo largo de la Baha de Paracas, hasta las
vecindades de Puerto San Martn (Transecto 14). Los valores de
turbidez obtenidos, ordenados por transectas, se presentan en la Tabla
4.2.4.1.

Desde el punto de vista de la Turbidez, la zona se puede dividir en


dos: al norte de Playa Lobera, con valores mayores a 6 NTU y al sur de
ese lugar con valores menores a 6 NTU.

ENVIRONMENTAL RESOURCES MANAGEMENT

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)

Oxigeno disuelto (mg/l)

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

ENVIRONMENTAL RESOURCES MANAGEMENT

15

Valores promedio por


Transecto
7.58

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

Valores promnedio por


Transecto
11.5

10.2

6.2

6.3

2.3

1.4

1.3

3.9

1.8

1.5

0.8

Los valores mximos registrados se encontraron inmediatamente al sur


del Ro Pisco con valores medios mayores a 10 NTU.
ENVIRONMENTAL RESOURCES MANAGEMENT

16

Los valores mnimos registrados se encontraron en los transectos 13 y


14, con valores medios menores a 1 NTU.
En cuanto a la transparencia del agua, la misma vari entre 3,8 a 3,2
metros de profundidad, con un promedio de 3,30 metros. Las
estaciones en donde se presentaron las ms bajas transparencias fueron
las ubicadas ms lejos de costa. En general, estos valores muestran
condiciones normales (EIA, Lnea de Base Ambiental, Oceanografa, ERM,
2002).

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

ENVIRONMENTAL RESOURCES MANAGEMENT

17

Figura 4.2.4.1. Resultados de las mediciones de turbidez a lo largo de la traza de la


obra.
Playa Lobera, Pisco

El valor mximo registrado corresponde al Bore Hole 10 (BH 10) con 15


NTU, medido a una profundidad de 11 m. El mnimo corresponde al
BH 12 con 5 NTU, medido a una profundidad de 6 m. El valor
promedio total durante el relevamiento fue de 8.4 NTU.

ENVIRONMENTAL RESOURCES MANAGEMENT

18

Tabla 4.2.4.2. Resultados de las mediciones de turbidez, slidos disueltos totales y


Disco Secchi a lo largo de la traza de la obra.

4.2.5. Mareas

Para determinar las caractersticas mareales de la zona en estudio, se


ha utilizado la Tabla de Marea que edita la Direccin de Hidrografa y
Navegacin de la Marina para el Puerto de Pisco, observndose que las
mareas son del tipo semi-diurno, es decir que se presentan dos
pleamares y dos bajamares en un da mareal (24 horas 50 minutos). La
amplitud media es de 0,58 metros, mientras que la amplitud en sicigias
(luna llena y/o nueva) es de 0,76 metros.

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

20

Figura 4.2.5.1. Posiciones de las estaciones de muestreo a lo largo de la traza de la obra, Playa Lobera, Pisco

ENVIRONMENTAL RESOURCES MANAGEMENT

21

Los instantes de ocurrencia y las alturas correspondientes a las pleamares y


bajamares para el Puerto de Pisco y El Callao puede obtenerse de la tabla
mencionada. A continuacin se presentan los valores correspondientes a las
amplitudes de sicigia y medias para cada uno de los sitios mencionados.

Puerto de Pisco
latitud:

13 43.0S

longitud:

76 14.0W

amplitud en sicigia: 0.76 m,


amplitud media:

0.58 m,

Puerto El Callao
latitud:

12 03.0S

longitud:

77 09.0W

amplitud en sicigia: 0.73m,


amplitud media:

0.55 m,

Comparando los instantes de ocurrencia de pleamares correspondientes a


distintas localidades del litoral peruano puede concluirse que la marea se
propaga de Norte a Sur.
4.2.6. Corrientes
El sistema de corrientes ocenicas costa afuera del Per tiene una direccin
predominante hacia el Norte y difiere de otro sistema adyacente a la
plataforma o costero, que presenta gran variabilidad espacio-temporal debido
principalmente a los accidentes geogrficos y otros factores como ser las
mareas, la topografa del fondo y los vientos locales.

Debido a la presencia de puntas o promontorios rocosos que sobresalen del


litoral, es frecuente la formacin de sistemas de circulacin rotatorios,
ENVIRONMENTAL RESOURCES MANAGEMENT

22

producto del rozamiento forzado de la masa de agua en su desplazamiento.


Esta situacin, similar a la que se produce en la baha del Callao, origina un
cambio localizado de la direccin de la corriente.

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.

Esta corriente hacia el sur, llamada tambin Contracorriente, parece


extenderse sobre la plataforma continental hasta donde comienza a
manifestarse una zona de surgencia (upwelling) costera asociada a la
presencia de los vientos alisios, confinada esta ltima a profundidades
inferiores a 50 m. Antecedentes recientes de esta contracorriente (Strub et al.,
Altimeter Observations of the Peru-Chile Countercurrent, Geophysical Research
Letters, 1995) indica que diferentes autores notaron un alto grado de
variabilidad tanto espacial como temporal de esta contracorriente y, por lo
tanto, todas las conclusiones relativas a la estructura espacial y variabilidad
estacional son hasta el momento tentativas.

Las corrientes costeras son dominadas por las mareas y los vientos.
4.2.7. Vientos

La zona de estudio est ubicada en una franja subecuatorial en la que


imperaran los vientos alisios con direccin SE. La regularidad en intensidad y
direccin de los alisios, que predomina sobre todo en regiones ocenicas, se ve
alterada por contrastes trmicos entre el continente y el ocano y por factores
orogrficos. Este es el caso de la regin costera de Pisco. En la zona estudiada
la direccin del viento cambia durante el da: durante la maana los vientos
soplan predominantemente del SW cambiando hacia la direccin N durante la
tarde. Por otra parte, ciertos das durante la tarde, se manifiestan vientos que
alcanzan intensidades de hasta 30 nudos, conocido regionalmente como
Paracas. Estos vientos ocurren fundamentalmente de julio a septiembre.

ENVIRONMENTAL RESOURCES MANAGEMENT

23

Los vientos histricos en Pisco, agrupados estadsticamente, se presentan en la


Tabla 4.2.7.1.

Tabla 4.2.7.1. Relacion velocidad del viento y direccion y sufrecuencia de ocurrencia, para el
Aeropuerto de Pisco.

4.2.8. Olas

En el rea de estudio no se tienen mediciones de olas disponibles. En el


informe (Sandwell - Chevron Texaco. NAWC LNG Receiving Terminal 142552,
Metocean Report, April 2002, Draft Report) se incluye un estudio del clima de
olas en el sitio de la obra.

Para la concrecin de este estudio fue necesario recurrir a datos de vientos


historicos

costa

afuera,

provenientes

de

la

ejecucin

de

modelos

meteorolgicos globales.

Una vez determinados estos vientos se procedio a modelar el clima de olas


(GROW2000) resultante de ese campo de vientos, obtenindose un campo
estadstico de olas costa afuera.

El campo estadstico de olas dentro de la Bahia de Pisco, en el sitio de la obra,


se calculo por medio de la utilizacin de un modelo matemtico de
transformacin espectral de olas llamado SWAN.

El clima de olas resultante, en las Bahia de Pisco, se incluye en las siguientes


tablas

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

25

5. INFORMACIN UTILIZADA PARA LA MODELACIN MATEMATICA

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

Con los valores de las constantes presentadas en la tabla anterior puede


calcularse el factor F (Dronkers, J.J. Tidal Computations, North Holland Publishing
Co, msterdam, 1965) el cual caracteriza el rgimen de la marea. Dicho
parmetro est definido de la siguiente manera:

F = (HK1 + HO1) / (HM2 + HS2)

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).

Esto se observa claramente en la curva de niveles del mar medidos (Report on


Phase II Off-Shore Site Investigation Required For Sub Sea Pipe Lpg Pisco Terminal ENVIRONMENTAL RESOURCES MANAGEMENT

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.

Mediciones del nivel del mar en Pisco.


Perodo: 5/03-14/03-2003
90
80
70

Nivel del Mar (cm)

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.

Todos los modelos hidrodinmicos requieren condiciones en los contornos o


bordes abiertos del dominio de cmputo. Para el caso de simulaciones con
forzantes mareolgicos esto siempre requiere de un tratamiento particular ya
que rara vez se dispone de mediciones directas de niveles del mar en aguas
abiertas. En estos casos, por lo general, los resultados de modelos globales son
de utilidad para representar las variaciones espacio-temporales de los niveles
del agua en los bordes del dominio.

Actualmente, el modelo global de marea TPXO.6 o TOPEX es el que mejor


ajusta (en el sentido de mnimos cuadrados) la ecuacin de marea de Laplace
con los datos satelitales resultantes de la misin TOPEX/Poseidon.

ENVIRONMENTAL RESOURCES MANAGEMENT

27

En este modelo global, la marea es representada mediante cuatro


componentes armnicas semidiurnas (M2, S2, N2 y K2), cuatro diurnas (K1, O1,
P1 y Q1) y dos de largo perodo (Mf y Mm) sobre una grilla global de 1440 x 721
nodos, con una resolucin espacial de 0.25, tanto en latitud como en longitud.

Los resultados del modelo, con intervalo de muestreo horario, se encuentran


disponibles en Internet. Una completa informacin sobre el modelo y un
instructivo acerca de cmo descargar datos puede obtenerse en:
http://www.oce.orst.edu/po/research/tide/global.html
En este sitio puede accederse, entre otros, a software especfico para
predecir alturas y corrientes de marea para lapsos y lugares requeridos.

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

Modelo Global TOPEX


Tabla de Marea en El Callao

0.9
0.8

Nivel del Mar (m)

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

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

Tabla de Marea en Pisco_PI

0.9

Modelo Global TOPEX


0.8
0.7

Nivel del Mar (m)

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.

De la comparacin surge la validacin del TOPEX para ser utilizado como


condicin de borde para el MIKE 21.

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

No obstante, dichas simulaciones deben realizarse considerando a todos los


forzantes presentes, siendo los ms importantes, la marea, el viento y, dado
este caso particular, los sistemas de corrientes marinas predominantes costa
afuera.

Para la zona de inters se dispone de informes tcnicos preliminares (ver


Tabla 5.2.1) que brindan informacin relativa a la circulacin de las aguas,
describen mediciones directas de corrientes de corta duracin obtenidas con
instrumental oceanogrfico y seguimiento de flotadores.

Tabla 5.2.1. Informes tcnicos de corrientes correspondientes a la zona en estudio


Ttulo
Ao
Estudio de corrientes marinas Junio
para la construccin de un
de 2002
muelle en la playa Loberia Pisco

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

Metocean Report: Pisco LPG


Terminal

Abril
de 2002

Pluspetrol
Peru
Corporation
S.A.

Sandwell
Engineering
Inc.

Estudio de Impacto
Ambiental - Oceanografa
(Lnea de Base)

2002

Off-Shore Site investigation


required for sub sea pipe
LPG Pisco Terminal
Camisea Project
Pisco Per
Task 4000 Ocean Currents
Phase II
Off-Shore Site nvestigation
required for sub sea pipe
LPG Pisco Terminal Camisea Project
Pisco Per
Task 1000 Ocean Currents

Dic. de
2002

Pluspetrol
Peru
Corporation
S.A.
Pluspetrol
Peru
Corporation
S.A.

No incluye datos. Breve


descripcin de la corriente
basada en las Pilots Charts
y en cartas nuticas
inglesas.
Descripcin general de las
corrientes de la zona

Marzo
de 2003

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

5.2.1. Mediciones directas eulerianas

Se realizaron mediciones de corrientes (Estudio de corrientes marinas para la


construccin de un muelle en la playa Loberia Pisco, Proyecto del Gas de Camisea, H &
O, junio de 2002.), a diferentes niveles de la columna de agua, en dos estaciones
de muestreo, utilizando para tal efecto un correntmetro porttil marca Valeport
Model 105 & 106 Self Recording, programado para obtener registros de
velocidad y direccin de la corriente a intervalos de 10 minutos y durante siete
das en cada estacin.

Las mediciones se realizaron, en la Estacin 1 (E1) a 5 y a 10 m sobre el fondo


marino donde la profundidad local es de 15.6 m, y en la Estacin 2 (E2) a 2 m
sobre el lecho marino, de profundidad local igual a 5.4 m.

La posicin de las estaciones fueron las siguientes:

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.

Los resultados de la medicin en forma estadstica se incluyen en las Tablas


5.2.1.1. a 5.2.1.3. y las Figuras 5.2.1.1. a 5.2.1.3, muestran los diagramas de
dispersin en trminos de las componentes N-S y E-W de las corrientes medidas.

Las velocidades de las corrientes superficiales ( a 5 m de la superficie del mar),


fueron menores de 0.30 m, a 5 m del fondo menores de 0.25 m/s , para una
altura de la columna de agua de 15 m.

El profundidades de 5 m y a 2 m del fondo las corrientes alcanzaron un


mximo de 0.25 m/s.

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

profunda. Adems, indican que la presencia de viento norte durante las


maanas produce un mayor ingreso de agua a la baha, originando cambios
en la circulacin.

La serie de tiempo de corrientes presenta caractersticas de ser afectada por


vientos y mareas. Sin embargo, no es evidente la contribucin de cada una de
esas fuerzas al valor de las corrientes. Para la calibracin del modelo, este
tema es crucial, debido a que es necesario calibrarlo con cada fuerza en
particular.

Se decidi utilizar un procedimiento de filtrado para tratar de separar, del


registro de corrientes, las contribuciones de la marea a las propias corrientes.
Para ello se tomaron dos series de mediciones de corrientes: una del 16 al 22
de mayo y la otra desde el 31 de mayo al 7 de junio de de 2002.

Durante el primer perodo, se calcularon las alturas con al MIKE 21 como se


muestra en la Figura 5.2.1.4. Puede observarse que la componente M2
(demidiurna de 12.4 h de perodo) no es importante. Esto se refleja en el
anlisis armnico de la serie de corrientes, el cual muestra como nica
componente significativa la onda K1 en la componente N-S, como muestra la
Figura 5.2.1.5. Cabe destacar que la energa de la componente K1 puede
explicar solamente el 11.1% de toda la energa del registro de corrientes, en
este perodo. Esto es, la contribucin de la marea al campo de corrientes no es
importante en este caso.

Seguidamente, se examin el segundo registro de corrientes. Aqu,


nuevamente se calcularon las alturas de marea por medio del MIKE 21, como
muestra la Figura 5.2.1.6. Puede observarse que las componentes M2 y K1 son
importantes en este caso. Esto se refleja en el anlisis armnico de la serie de
corrientes, donde ambas componentes (E-W y N-S) resultan significativas,
como lo muestran las Figuras 5.2.17. y 5.2.1.8.

La energa de las componentes M2 y K1 (y una dbil M3) pueden explicar el


8.1% de toda la energa del registro de corrientes en la componente E-W y el
25.9% en la componente N-S.

ENVIRONMENTAL RESOURCES MANAGEMENT

32

Tabla 5.2.1.1. Estadistica de Corrientes Marinas en la Bahia de Pisco (E1)


Profundidad : 10 metros (16 al 22 de mayo), Frecuencia (%) de incidencia de
corrientes

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

Tabla 5.2.1.2. Estadistica de Corrientes Marinas en la Bahia de Pisco (E1)


Profundidad : 5 metros (24 al 31 de mayo), Frecuencia (%) de incidencia de corrientes
Velocidad
(cm/s)
0-5
5-10
10-15
15-20
20-25
25-30
30-40
>40
TOTAL %

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

Tabla 5.2.1.3. Estadstica de Corrientes Marinas en la Bahia de Pisco (E2)


Profundidad : 2m (31 de mayo al 7 de junio), Frecuencia (%) de incidencia de
corrientes
Velocidad
(cm/s)
0-5
5-10
10-15
15-20
20-25
25-30
30-40
>40
TOTAL %

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

ENVIRONMENTAL RESOURCES MANAGEMENT

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

Mediciones de corrientes con instrumento.


Perodo: 16/05-22/05/2002, Intervalo de muestreo: 10 minutos
0.25

0.20

Componente N-S de la Velocidad (m/s)

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

Componente E-W de la Velocidad (m/s)

Figura 5.2.1.1. Mediciones de corrientes frente a Playa Lobera. Estacin E1, 15 m de


profundidad local, instrumento ubicado a 10 m del fondo del mar. Diagrama de dispersin
componentes E-W y N-S de las corrientes.

Mediciones de corrientes con instrumento.


Perodo: 24/05-31/05/2002, Intervalo de muestreo: 10 minutos
0.40
0.35
0.30

Componente N-S de la Velocidad (m/s)

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

Componente E-W de la Velocidad (m/s)

Figura 5.2.1.2. Mediciones de corrientes frente a Playa Lobera. Estacin E1, 15 m de


profundidad local, instrumento ubicado a 5 m del fondo del mar Diagrama de dispersin
componentes E-W y N-S de las corrientes.

ENVIRONMENTAL RESOURCES MANAGEMENT

34

Mediciones de corrientes con instrumento.


Perodo: 31/05-07/06/2002, Intervalo de muestreo: 10 minutos

0.25

0.20

Componente N-S de la Velocidad (m/s)

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

Componente E-W de la Velocidad (m/s)

Figura 5.2.1.3. Mediciones de corrientes frente a Playa Lobera. Estacin E2, 5 m de


profundidad local, instrumento ubicado a 2 m del fondo del mar Diagrama de dispersin
componentes E-W y N-S de las corrientes.

Nivel del mar para el perodo 16 al 22 de Mayo de 2002


0.7

0.5

Nivel del Mar (m)

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

35

Estimacin de corrientes de marea - componente N-S


Porcentaje de varianza explicado de la serie original: 11.1% (K1)
Perodo 16 al 22 de Mayo de 2002

Componente N-S de corriente de marea (cm/s)

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.

Nivel del mar para el perodo 31 de Mayo al 7 de Junio de 2002


0.5

0.4

0.3

Nivel del Mar (m)

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

36

Estimacin de corrientes de marea - componente E-W


Porcentaje de varianza explicado de la serie original: 8.1% (K1, M2, M3)
Perodo: 31 de Mayo al 7 de Junio de 2002
0.06

Componente E-W de corriente de marea (cm/s)

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.

Estimacin de corrientes de marea - componente N-S


Porcentaje de varianza explicado de la serie original: 25.9% (K1,M2,M3)
Perodo: 31 de Mayo al 7 de Junio de 2002
0.04

Componente N-S de corriente de marea (cm/s)

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

37

5.2.2. Corridas de flotadores lagrangeanas

A continuacin, se sintetizan los resultados obtenidos mediante corridas de


flotadores realizadas en baha Pisco en inmediaciones de Lobera y Tambo de
(Report on Phase II Off-Shore Site Investigation Required For Sub Sea Pipe Lpg Pisco
Terminal - Camisea Project Pisco Per Task 1000 - Ocean Currents, marzo de
2003). La Tabla 5.2.2.1. resume las condiciones generales y caractersticas de
las corridas de Flotadores.

De la informacin disponible provista por las corridas de flotadores dentro de


la baha, se puede concluir lo siguiente:

La velocidad y direccin de la boya est altamente influenciada por el


efecto del viento.

Por lo general la boya con veleta superficial se mueve ms rpidamente


que las que tienen veleta en profundidad. En consecuencia se concluye que
las aguas superficiales se mueven ms rpidamente que las profundas
(esto es compatible con las observaciones eulerianas).

Ciertos das, dentro de la baha, se manifiestan corrientes no asociadas al


viento. Estas corrientes provienen del Sur y entran en la baha bordeando
la pennsula Pejerey.

Las boyas arrojadas cerca de la costa muestran un comportamiento ms


errtico comparativamente a las de aguas profundas. Por lo general, las
boyas con veletas profundas presentan una tendencia a moverse
predominantemente hacia la lnea de costa.

Con el objeto de obtener mediciones fuera de la baha para contribuir a la


calibracin del modelo, se planearon y ejecutaron (Report on Phase II Off-Shore
Site Investigation Required For Sub Sea Pipe Lpg Pisco Terminal - Camisea Project
Pisco Per Task 1000 - Ocean Currents, marzo de 2003) 3 corridas de flotadores
en los siguientes lugares: Transecta A, Transecta B y Tambo de Mora, como
muestra la Figura 5.2.2.1. Las Figuras 5.2.2.2. a 5.2.2.4. muestran los resultados
de las corridas en forma vectorial.

ENVIRONMENTAL RESOURCES MANAGEMENT

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

Profundidad de la veleta (metros)


3
6
9

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.
_______________________________________

ENVIRONMENTAL RESOURCES MANAGEMENT

39

Tabla 5.2.2.1. CONTINUACIN. 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
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.)

Profundidad de la veleta (metros)


3
6
9

Preponderan Predominant Predominant Predominant


temente
e-mente
e-mente
e-mente
semidiurna.
hacia el S.
hacia el S.
hacia el S.
Amplitud:
60 cm.
Desigualdad
es diurnas.
Amplitud:
37 / 85 cm.

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.
_______________________________________

Los resultados en forma de diagrama de dispersin para las componentes E-W


y N-S, se presentan en las Figuras 5.2.2.5. a 5.2.2.7.

De la informacin disponible provista por las corridas de flotadores fuera de


la baha, se puede concluir lo siguiente:

Las mediciones realizadas en la transecta A (Boquern) muestra que las


corrientes van hacia el Sur aun cuando el viento viene del sur, es decir, en
contra de la corriente.

Las mediciones desarrolladas sobre la transecta B muestran que las


corrientes fluyen principalmente hacia el Suroeste.

ENVIRONMENTAL RESOURCES MANAGEMENT

40

En Tambo de Mora las corrientes fluyen hacia el Sur con baja intensidad.

Estas mediciones, en profundidades entre 50 m y 60 m, indicaron velocidades


de corrientes con componentes N-S muy 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).

En consecuencia, las velocidades medidas por los flotadores tienen explicacin


solamente en la Contracorriente del Per, por lo menos costa afuera en los
lugares de medicin Transectas A y B y Tambo de Mora.

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

42

Mediciones de corrientes con flotadores. Transecta A.


Perodo: 06:20-16:30, 12/03/2003
0.10
0.05
0.00

Componente N-S de la Velocidad (m/s)

-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

Componente E-W de la Velocidad (m/s)

Figura 5.2.2.5. Mediciones de corrientes con flotadores en la Transecta A. Diagrama de


dispersin componentes E-W y N-S de las corrientes.

Mediciones de corrientes con flotadores. Transecta B.


Perodo: 05:20-16:10, 13/03/2003

0.00

-0.05

Componente N-S de la Velocidad (m/s)

-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

Componente E-W de la Velocidad (m/s)

Figura 5.2.2.6. Mediciones de corrientes con flotadores en la Transecta B. Diagrama de


dispersin componentes E-W y N-S de las corrientes.

ENVIRONMENTAL RESOURCES MANAGEMENT

43

Mediciones de corrientes con flotadores. Transecta Tambo de Mora.


Perodo: 05:50-16:30, 14/03/2003
0.20

0.15

Componente N-S de la Velocidad (m/s)

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

Componente E-W de la Velocidad (m/s)

Figura 5.2.2.7. Mediciones de corrientes con flotadores en la Transecta Tambo de Mora.


Diagrama de dispersin componentes E-W y N-S de las corrientes

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.

El objetivo de la calibracin y validacin es reproducir las corrientes


observadas, con intensidades y direcciones en el rango de los valores
esperados. Dado que las corrientes muestran valores sostenidos de direccin
en el tiempo de medicin con algunas variaciones de intensidad durante ese
perodo (unas 12 h), se ha decidido representarlas por medio de sus valores
medios, como indica la Tabla 5.2.2.2.
Tabla 5.2.2.2. Valores medios de la velocidad durante el perodo de modelacin, del orden
de 12 h, para los sitios de corridas de flotadores en Transecta A, Transecta B y Tambo de
Mora.
Valores Medios de la
Transecta A
Transecta B
Tambo de Mora
Velocidad
Componente E-W de la
Velocidad (m/s)
Componente N-S de la
Velocidad (m/s)
Notas

-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

ENVIRONMENTAL RESOURCES MANAGEMENT

44

Un ejemplo de corrida de flotadores para la zona de la Pennsula de Paracas se


puede ver en la Figura 5.2.2.8.

Figura 5.2.2.8. Ejemplo de corrida de flotadores en el la zona proyectada para la Terminal


Marima. Los flotadores se dirigen hacia en N con distintas profundidades de sus palas de
arrastre.

5.2.3. Modelos globales de corrientes

Dada la falta de datos de corrientes costa afuera y la necesidad del modelo de


condiciones de contorno en sus borde, fue necesario recurrir a un modelo
global de corrientes. Algunos detalles del mismo se explican a continuacin.

La circulacin de las aguas de Ocano Pacfico Sur, adyacente a las costas


peruanas, es representada por el modelo OCCAM con un relativamente alto
grado de detalle. El Proyecto Ocean Circulation and Climate Advanced
Modelling (OCCAM) est siendo desarrollado por la Natural Environment
Research Council (NERC) de United Kindom y muchos de los resultados son
pblicos y estn disponibles en Internet.

Este proyecto involucra a investigadores de universidades e institutos


cientficos de UK, estando el grupo de cientficos principales en el
Southampton Oceanography Centre. Uno de los objetivos de este proyecto
ENVIRONMENTAL RESOURCES MANAGEMENT

45

es el desarrollo de modelos ocenicos de alta resolucin para contribuir, por


un lado, al estudio de la circulacin ocenica de gran escala y, por otro,
interpretar los datos del World Ocean Circulation Experiment (WOCE).

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.

El modelo considera que la velocidad es cero en todos los contornos slidos y


que los gradientes de temperatura potencial y salinidad normales a dichos
contornos (incluyendo al fondo) son tambin cero. En la superficie libre se
utilizan campos atmosfricos para forzar al ocano. Estos incluyen tanto a los
campos de viento y de presin atmosfrica, como as tambin, a los flujos de
calor y agua a travs de la superficie del mar. Sobre el lecho ocenico el
esfuerzo de fondo horizontal depende de la corriente correspondiente a la
capa de fluido adyacente al fondo.

El modelo OCCAM tiene una resolucin espacial de 0.25 tanto en latitud


como en longitud y posee treinta y seis niveles en la vertical. Los espesores de
las capas varan entre 20 m cerca de la superficie hasta 255 m a la profundidad
de 5500 m.

El esfuerzo del viento utilizado como forzante corresponde a interpolaciones


lineales temporales de los promedios mensuales climatolgicos . Los flujos de
calor y de agua en superficie fueron calculados mediante la aplicacin de
tcnicas de relajacin numrica sobre la capa superior del modelo.

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

ENVIRONMENTAL RESOURCES MANAGEMENT

46

5.2.4. Estado actual del conocimiento

Sobre

la

base

de

las

publicaciones

informes

tcnicos

citados

precedentemente se concluye que las corrientes en la zona de estudio son


variables e inducidas principalmente por el viento, la marea y las corrientes de
contorno.

Dado que no se dispone de registros de corrientes de al menos un mes de


duracin, que permitan discriminar todos los efectos de la corriente marea
mediante anlisis armnico, la cuantificacin de los efectos de cada uno de los
mecanismos mencionados se hace altamente dificultosa. En consecuencia, es a
priori aventurado hablar de patrones de circulacin caractersticos para la
zona en estudio en base a las mediciones disponibles. No obstante, los datos
de campo obtenidos permiten cuantificar las intensidades de corrientes
esperadas para la zona.

Cabe sealar el papel preponderante que juegan los accidentes geogrficos


costeros en la zona de inters. Es de prever que la Pennsula de Paracas
produzca un giro de las aguas tanto para corrientes predominantes del S como
as tambin del N. Este efecto hidrodinmico es vastamente conocido por los
oceangrafos costeros y su interpretacin es an ms compleja en presencia de
vientos variables y bajo flujos rotatorios como los producidos por la marea.

La aplicacin de una herramienta objetiva como MIKE 21 produce cartas de


corrientes asociadas a forzantes especficos y, de esta manera, puede evaluarse
la preponderancia de cada uno de ellos.

Para ello, se plantearn escenarios caracterizados por forzantes mareolgicos,


atmosfricos y de corrientes de contorno con el objeto de simular las corrientes
asociadas a cada uno de ellos. A los efectos de validar los resultados, las
corrientes obtenidas se compararn con las observaciones de campo
disponibles.

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

Vientos en Pisco. Serie de tiempo de la componente E-W del viento.


Periodo: 6/3 al 14/3/2003, Numero de casos:N=216, Intervalo de muestreo: 1 hora
10
9
8

Componente E-W del Viento (m/s)

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.

Dada la marcada periodicidad de los vientos medidos, se decidi investigar


espectralmente estos vientos.

Las Figuras 5.3.3. y 5.3.4. muestran la representacin espectral de las


componentes de la velocidad del viento, para el perodo considerado. Se
puede observar claramente un aporte importante de la energa ubicado en
frecuencias algo menores de 0.05 1/h (esto es perodos algo mayores a 20 h)
ENVIRONMENTAL RESOURCES MANAGEMENT

48

Vientos en Pisco. Serie de tiempo de la componente N-S del viento.


Periodo: 6/3 al 14/3/2003, Numero de casos:N=216, Intervalo de muestreo: 1 hora
8
7
6
5

Componente N-S del Viento (m/s)

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.

Vientos en Pisco. Representacion espectral de la componente E-W del viento.

Energia *N/2 (m2/s2)

Periodo: 6/3 al 14/3/2003, Numero de casos: N=216, Intervalo de muestreo: 1 hora


700

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.

Energia *N/2 (m2/s2)

Periodo: 6/3 al 14/3/2003, Numero de casos:N=216, Intervalo de muestreo: 1 hora


700

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

49

Estas caractersticas peridicas en el comportamiento del viento son


potencialmente significativas al momento de la modelacin matemtica. La
Figura 5.3.5. Muestra un diagrama de dispersin de las componentes E-W y
N-S del viento , para el perodo de medicin.

Vientos en Pisco. Diagrama de dispesin direccin-velocidad del viento.


Periodo: 6/3 al 14/3/2003, Numero de casos:N=216, Intervalo de muestreo: 1 hora
10
9
8

Velocidad del Viento (m/s)

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

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

hidrodinmicas/ambientales en lagos, estuarios, bahas, reas costeras y


mares.

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.

El mdulo hidrodinmico (HD) es el mdulo bsico en el paquete de


programas del MIKE 21. Con este mdulo se pueden simular fluctuaciones del
nivel del agua y flujos asociados a una gran variedad de forzantes como, por
ejemplo, marea, viento y corrientes. El modelo resuelve las ecuaciones de
conservacin de masa y de cantidad de movimiento verticalmente integradas,
las cuales incluyen los trminos de friccin de fondo, esfuerzos del viento,
gradientes de presin baromtrica, fuerza de Coriolis y mezcla turbulenta,
entre otros.

Con el modelo MIKE 21 pueden obtenerse cartas de corrientes asociadas a


forzantes especficos las cuales permiten evaluar el efecto de cada uno de ellos.
Se

plantearn

escenarios

caracterizados

por

forzantes

mareolgicos,

atmosfricos y de corrientes de contorno con el objeto de simular la circulacin


de las aguas en la zona de la baha de Pisco. A los efectos de validar los
resultados, las corrientes obtenidas se compararn con las observaciones de
campo disponibles. Finalmente, se discutirn los principales patrones de
circulacin obtenidos.

ENVIRONMENTAL RESOURCES MANAGEMENT

51

6.2. Antecedentes. Modelacin del Campo de Corrientes en la Zona


La Universidad Nacional Mayor de San Marcos (Yzocupe, V. Modelo Numrico
de Flujo Bidimensional de la Baha de Paracas, Laboratorio de Fluidodinmica,
Departamento de Ingeniera Mecnica de Fluidos, Facultad de Ciencias Ficas) ha
construido

implementado

un

modelo

numrico

hidrodinmico

bidimensional de la Baha de Paracas.

En este trabajo los autores aplican el modelo con el objeto de estudiar la


circulacin de las aguas ante forzantes mareolgicos y meteorolgicos. Se
utiliz una discretizacin espacial de 926 m, un parmetro de Courant de 0.5,
un incremento temporal de 7.09 s. y valores variables para un coeficiente de
rugosidad de Manning comprendidos entre 0.02 y 0.05.

En la frontera oeste se utiliz una condicin de borde resultante del anlisis


armnico de dos series de una semana de valores horarios del nivel del mar
registrados en la estacin Pisco. Las ondas entran al dominio con un frente de
76. En los contornos norte y sur se utilizaron condiciones de radiacin.

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.

6.3. Modelacin 2D vs. 3D


Los modelos hidrodinmicos bidimensionales, verticalmente integrados,
conocidos comnmente como 2D, son de relativamente fcil calibracin,
econmicos (desde el punto de vista de los tiempos de ejecucin), brindan
resultados de interpretacin directa y son robustos (es decir, pueden
adaptarse casi sin restricciones a cualquier rea somera) por lo cual son
altamente recomendados para estudios de circulacin costera. Sin embargo, no
son aplicables para la simulacin de, por ejemplo, celdas de circulacin
ENVIRONMENTAL RESOURCES MANAGEMENT

52

vertical, surgencias (upwelling) o sitios que presenten cambios importantes


de la corriente con la profundidad, en trminos de velocidad pero ms
importante en la direccin.

Los modelos hidrodinmicos tridimensionales, multicapas conocidos como


modelos 3D permiten estudiar la variacin de la corriente con la profundidad
y analizar los campos de corrientes verticales brindando un panorama
completo de la circulacin de las aguas. Sin embargo, tienen una mayor
cantidad de variables de ajuste por lo cual son de compleja calibracin (deben
ajustarse las salidas del modelo con perfiles de velocidades medidos en el
campo), debe gastarse un esfuerzo adicional para la elaboracin y
determinacin de las condiciones de borde (estn generalmente anidados a
modelos 3D regionales que los alimentan) y son ms propensos a
inestabilidades numricas cuando se incorporan los campos trmicos y salinos
con las estratificaciones verticales correspondientes.

Considerando todos los requerimientos que necesita la calibracin de los


modelos 3D estos frecuentemente se aplican a sitios para los cuales se dispone
de batimetras de detalle y se posee de un conocimiento cuantitativo y
cualitativo de la circulacin general bidimensional (es decir, se conocen los
patrones bsicos de circulacin y las caractersticas y variabilidad de los
forzantes principales)

Para el caso de la zona de inters no se dispone de mediciones directas del


perfil de la corriente en ningn sitio. Por otro lado, aunque a priori se deba
considerar a la marea, a las corrientes de borde y al viento como no se conoce
exactamente el efecto que cada uno de ellos produce sobre la circulacin de las
aguas en la zona de baha Pisco. Tal es el caso que, al no disponerse de
mediciones prolongadas (ms de un mes) de corrientes, la discriminacin de
los efectos de cada uno de los mecanismos mencionados se hace altamente
dificultosa.

Adems, aunque los antecedentes disponibles indiquen que en algunos sitios


particulares la columna de agua presenta algun tipo de estratificacin, no se
dispone de mediciones exhaustivas del campo trmico y salino en el dominio

ENVIRONMENTAL RESOURCES MANAGEMENT

53

estudiado que permitan elaborar una condicin inicial realista para un modelo
3D.

Como en el presente estudio se pretende determinar los patrones principales


de circulacin costera, de escasa profundidad, para evaluar los posibles
comportamientos de hipotticos derrames de combustible o dispersin de
partculas de sedimento producido por dragado y, de alguna manera,
establecer cuales son los forzantes dominantes de dicha circulacin, la
aplicacin de un modelo HD 2D parece ser la herramienta ideal a utilizarse en
este trabajo.

6.4. Sistema de Coordenadas y Orientacin de la Grilla


El modelo trabaja sobre un sistema de coordenadas rectangulares planas x-y
(grilla). Dado que la costa presenta accidentes geogrficos significativos como
por ejemplo la pennsula Paracas, no parecer existir una orientacin particular
del sistema de referencia que posibilite una mejor representacin de los
accidentes geogrficos costeros. En consecuencia, se adopt un sistema de
referencia con orientacin tradicional: con el eje x apuntando hacia el Este y el
eje y hacia el Norte.

El dominio considerado para la modelacin numrica se extiende


aproximadamente desde 13 22 S hasta 14 02 S, en la direccin N-S y desde
76 40 W hasta la costa, en la direccin W-E.

El origen del sistema de referencia corresponde exactamente al vrtice SW del


dominio del modelo (14 01 42.6S, 76 40 W).

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

54

6.5. Referencia Temporal y Manejo del Tiempo


El modelo MIKE 21 necesita la siguiente informacin bsica para el correcto
manejo del tiempo:

Cantidad de pasos temporales necesarios para completar la simulacin


numrica.

El paso temporal (es decir, el incremento temporal entre pasos


sucesivos). En este estudio se utiliz un paso temporal igual a 10 y15 s

Sistema de referencia temporal. Dado que en este trabajo se utiliza


informacin mareolgica y de corrientes provenientes de los modelos
globales TOPEX y OCCAM, respectivamente, los cuales se manejan
con hora GMT, y observaciones de campo variadas, se adopt un
sistema de referencia temporal nico el cual corresponde al de la hora
oficial peruana para sincronizar toda la informacin disponible.

6.6. Perodo de Estabilizacin


El perodo de estabilizacin (frecuentemente denominado warm-up period)
es el nmero de pasos temporales sobre los cuales el forzante se incrementa
gradualmente desde cero hasta alcanzar el 100% de su valor real.

El nmero de pasos adoptado para el perodo de estabilizacin en las


presentes

simulaciones

fue

entre

2880

4320.

Consecuentemente,

considerando uno de los pasos temporales utilizados (15 s), a modo de


ejemplo, el tiempo de estabilizacin resultante fue de 12 h. Esta porcin de las
simulaciones fueron descartadas del anlisis de los resultados.

Para la realizacin de las simulaciones numricas se adopt la modalidad de


arranque en fro (conocida como cold start mode) para la cual se utiliza un
campo de corrientes iniciales nulo y un nivel de la superficie libre constante
para todo el dominio modelado.

ENVIRONMENTAL RESOURCES MANAGEMENT

55

6.7. Nmero de Courant


El sistema computa el nmero de Courant mximo resultante de la batimetra
y de la resolucin temporal adoptada, resultando para el caso estudiado un
valor igual a 2.4 y 1.6 (intervalo de tiempo igual 15 s y 10 s respectivamente).

El DHI recomienda reducir el paso temporal si el nmero de Courant excede


el valor 8-10. Considerando la batimetra utilizada y la discretizacin temporal
propuesta el modelo result estable en todas las simulaciones realizadas.

6.8. Batimetra
Con

los

datos

de

profundidad

obtenidos

de

las

cartas

indicadas

anteriormente, se gener (en el dominio de modelacin) una grilla batimtrica


aplicando una tcnica de interpolacin espacial provista por software SURFER
7.0. la Figura 6.8.1. muestra el dominio y la batimetra.

6.9. Forzantes
6.9.1. Bordes

El modelo MIKE 21 necesita condiciones especficas de niveles del mar o de


flujo de agua a travs de todos los nodos correspondientes a los bordes
abiertos del recinto estudiado.

En el presente estudio se consideraron tres contornos abiertos: dos dispuestos


a lo largo de paralelos (lmites Norte y Sur) y uno a lo largo de un meridiano
(lmite Oeste). A continuacin se indican las posiciones y coordenadas de los
tres contornos abiertos.

Contorno inferior del dominio, Y = 0 (lmite Sur del rea estudiada,


aproximadamente en la latitud 14 02 S).

ENVIRONMENTAL RESOURCES MANAGEMENT

56

Figura 6.8.1 Dominio de modelacin y batimetra. El punto al norte de la Baha de Paracas


corresponde al extremo oeste de la obra de la Terminal Martima proyectada. Puede observada
la discretizacin de los contornos con el caracterstico dentado.

Contorno superior del dominio, Y = 146 (lmite Norte del rea estudiada,
aproximadamente en la latitud 13 22 S).

Contorno lateral del dominio, X = 0 (lmite Oeste del rea estudiada, 76


40 W).

6.9.2. Mareas

Para la generacin de las condiciones de borde de marea a lo largo de los tres


contornos abiertos del recinto se utilizaron series de tiempo de valores
horarios del nivel del mar generadas con diez constantes armnicas de marea
obtenidas del modelo global TOPEX.

ENVIRONMENTAL RESOURCES MANAGEMENT

57

Se utilizaron para la generacin de la condicin de borde cuatros series


temporales de niveles del mar: dos correspondientes a los vrtices NW y SW
del dominio (mar afuera) y dos en puntos costeros (intersecciones entre los
contornos inferior y superior con la lnea de costa).

La determinacin de los niveles del mar en los puntos intermedios de los


contornos se realiz mediante la interpolacin lineal (provista por el MIKE
21), ponderada por la distancia, utilizando las series disponibles en los vrtices
y la costa.

6.9.3. Corrientes

Para la generacin de las condiciones de borde de corrientes en los nodos


correspondientes a los tres contornos abiertos del recinto se utilizaron
transportes o caudales especficos (en m3s-1/m) quincenales resultantes del
modelo global ocenico OCCAM, dada la ausencia de mediciones directas de
corrientes.

Como la resolucin de este modelo global (0.25 de latitud y longitud,


aproximadamente 27 km) es mucho menor que la del dominio de este estudio
(500 m), solamente 3 nodos del OCCAM coincidieron dentro del recinto de
modelacin. En consecuencia, fue necesario extrapolar hasta la costa e
interpolar entre los nodos disponibles del OCCAM.

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

58

6.9.5. Forzantes combinados


A la hora de calibrar el modelo numrico se deben considerar e incluir los
forzantes de significacin, simultneamente.

Los campos de velocidades resultantes de estas simulaciones deben ajustarse


lo mejor posible a las observaciones de campo. Al proceso de comparacin
entre las salidas del modelo y observaciones se lo denomina calibracin. El
modelo una vez calibrado se transforma en una herramienta eficaz para el
anlisis de eventos o situaciones particulares y puede ser utilizado como un
instrumento para toma de decisiones y de manejo costero.

6.10. Viscosidad Turbulenta


El coeficiente de viscosidad turbulenta (o eddy) fue considerado variable
siguiendo una de las opciones que suministra el MIKE 21. El valor de la
viscosidad resulta de los valores de los gradientes locales de velocidad.

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

coeficientes de friccin adoptndose finalmente el que mejor reproduca a las


observaciones de campo. El DHI recomienda utilizar valores comprendidos
entre 20 y 40 m0.33/s argumentando que en general los mejores resultados se
obtienen utilizando valores comprendidos en ese rango. Luego de realizar
diversos ensayos se adopt un valor igual a 25, constante para todo el recinto
modelado.

ENVIRONMENTAL RESOURCES MANAGEMENT

59

6.12. Calibracin del Modelo


Con el objeto de obtener resultados realistas del MIKE 21, se procedi a la
calibracin del mismo. Este procedimiento consiste en correr el modelo por
primera vez con valores aproximados de sus parmetros y condiciones de
contorno.

Una vez corrido el modelo y obtenidos sus resultados, se procedi a


compararlos con mediciones directas de corrientes y marea y con los
productos de los modelos globales. Cuando la coincidencia fue pobre se
procedi a ajustar los parmetros del modelo, correrlo nuevamente y as
sucesivamente hasta que se encontr un compromiso entre los resultados del
modelo y las mediciones y modelos globales.

Ciertamente, en este tipo de procedimiento, la calidad de la calibracin es


mayor cuando se dispone de mediciones de largo perodo en puntos
distribuidos a lo largo y ancho del dominio de modelacin. En la mayora de
las aplicaciones en todo el mundo, esto no es posible. Se dispone, en general,
de mediciones costeras y por perodos cortos de tiempo, principalmente de
mareas y vientos, pero no de corrientes.

Para la calibracin del modelo se consideraron:

Datos de alturas de plemares y bajamares obtenidas de Tablas de Marea


para Pisco, correspondientes a mayo y junio de 2002.

Mediciones de niveles del mar en Pisco obtenidas en marzo de 2003.

Mediciones de corrientes obtenidas con correntgrafo (2002).

Mediciones de trayectorias de corrientes medidas con flotadores (2002-03).

En este caso, durante el 2002 y 2003, se realizaron una serie de mediciones de


corrientes en la zona elegida para la construccin de la Terminal Martima y
unas mediciones costa afuera de corta duracin. Se dispuso de datos de
mareas provenientes de la estacin de Pisco y vientos para esos perodos de
medicin.

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

Se sabe de la existencia de la Contracorriente del Per por trabajos cientficos


publicados y de los resultados del Modelo Global OCCAM. No fue posible
encontrar mediciones directas de esta contracorriente aunque se le asocian
velocidades en superficie de hasta 0.5 m/s, a unos 100-300 km de la costa.

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).

En consecuencia, las velocidades medidas por los flotadores tienen explicacin


solamente en la Contracorriente del Per, por lo menos costa afuera en los
lugares de medicin y en el perodo de medicin.

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.

Una vez calibrado y corrido el modelo con esta contracorriente, se pudo


observar que no tiene un efecto de significacin dentro de las bahas, donde
las corrientes son dominadas por otras fuerzas como las mareas y el viento,
ENVIRONMENTAL RESOURCES MANAGEMENT

61

especialmente este ltimo. Se observ tambin que la contracorriente influye


en la generacin de celdas (fuera de la Baha de Paracas) de circulacin con
determinadas condiciones de viento.

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.

Se decidi analizar la capacidad del MIKE21 en simular la marea astronmica.


Para ello se forz al modelo con niveles del mar extrados del modelo global
TPX0.6 y se simul el lapso de tiempo correspondiente.

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

Nivel del Mar (m)

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

62

En la Figura 6.12.1. se presentan los resultados del MIKE 21 para el nodo


correspondiente a la zona de Pisco. En dicha figura se comparan los resultados
del modelo con los valores tabulados. Puede apreciarse que el modelo sigue
los valores predichos con una alta concordancia tanto en la fase como as
tambin en las alturas. Se concluye, en consecuencia, que el MIKE 21
reproduce satisfactoriamente las variaciones del nivel del mar asociadas a la
marea.

Se realiz una prueba comparando los resultados del modelo TOPEX y el


MIKE 21 para Pisco. La Figura 6.12.2. muestra los resultados. Otra prueba para
un punto costa afuera se realiz para el objeto de comparacin, como lo indica
la Figura 6.12.3. Los resultados son satisfactorios en ambos casos.

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

Nivel del Mar (m)

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

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

Nivel del Mar (m)

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

Se realiz una comparacin entre las mediciones de corrientes en Transecta A


y B y Tambo de Mora con los resultados del MIKE 21. Dado que los flotadores
dispusieron sus palas de arrastre a unos 6 m de profundidad, las corrientes
medidas por los mismos son mayores que las velocidades medias calculadas
por el MIKE 21. En base a las corridas de flotadores, con palas de arrastre a
distintas profundidades, se redujeron en 10%, para la comparacin. Este valor
fue determinado a partir de la estimacin de la distribucin vertical de
velocidades provistas por las corridas de flotadores. Una de esas
distribuciones ms regulares se muestra en la Figura 6.12.4.

Debe aclararse que este mtodo para el cmputo de la distribucin de


velocidades no es el ms idneo, pero el nico disponible dado las mediciones
realizadas. Los flotadores se mueven permanentemente y por lo tanto el perfil
de velocidades instantneo en un punto no es posible calcularlo. La Tabla
6.12.1. presenta los resultados de la calibracin.

Complementariamente se realiz una comparacin entre los resultados de las


corridas de flotadores para ciertas condiciones y se las compar con los
resultados provistos por el MIKE 21. La Tabla 6.12.2. muestra los resultados.

ENVIRONMENTAL RESOURCES MANAGEMENT

64

Perfil de velocidades en la isobata de 10 metros. Corrida BH-10.


0.0
-0.1
-0.2

Profundidad Normalizada (z/h)

-0.3
-0.4
-0.5
-0.6
-0.7
-0.8
-0.9
-1.0

Velocidad de la Corriente (cm/s)

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.

Las velocidades simuladas con el MIKE 21 (forzado con corrientes de borde y


viento) son comparadas con las intensidades de la corrientes medias
estimadas de las corridas de flotadores.

Otras comparaciones con flotadores, en corridas con datos confiables de


vientos, se realizaron a pesar que el MIKE 21 calcula las corrientes para puntos
determinados de la grilla y no sigue partculas como son los flotadores. Los
resultados se presentan en la Tabla 6.12.2.

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

66

7. ESCENARIOS DE MODELACIN Y RESULTADOS


7.1. Escenarios de Modelacin
Los escenarios de modelacin fueron seleccionados con el objeto de disponer
de un campo de velocidades realista aun cuando la frecuencia de ocurrencia
de los vientos sea muy baja, dado que estudios posteriores como el
desplazamiento de sedimentos en suspensin, requerirn de todos los
escenarios capaces de comprometer la biota, por ms remotos que sean.

Los escenarios seleccionados, habindose reducido el nmero de casos de


magnitud de la velocidad, son los siguientes:

Tabla 7.1.1. Escenarios de modelacin seleccionados


Viento
Escenario

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

VARIABLE SEGN MEDICIONES

15

SIN MAREA NI VIENTO

16

MAREA SOLAMENTE

ENVIRONMENTAL RESOURCES MANAGEMENT

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:

Los vientos principalmente, tanto aquellos sostenidos en el tiempo como


aquellos peridicos con periodicidad diurna manifiesta;

Las mareas, que influyen moderadamente en el campo de corrientes con


componentes peridicas semidiurnas y diurnas principalmente.

La Contracorriente del Per, con un aporte apenas perceptible en la


regin costera pero con influencia en el patrn de circulacin cercano a
ese lugar y dominante en aguas ms profundas.

Los accidentes geogrficos y batimetra asociada que estimulan patrones


de circulacin cerrada o eddy manifiestos aisladamente a lo largo de la
costa y especialmente en la Baha de Paracas, con condiciones particulares
de vientos.

Se ha investigado el campo de corrientes en la zona proyectada para la


Terminal Martima con una serie de tiempo de vientos medidos, durante
marzo de 2003.

Las corrientes han sido calculadas por el MIKE 21 en presencia de la


Contracorriente del Per y vientos como se ha mencionado. Los vientos
presentan

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

68

Serie de tiempo de corrientes en un lugar equidistante entre la costa y


el extremo oeste de la Terminal Martima
Terminal Martima. Corrientes forzadas por una serie de tiempo de viento medido.

0.08

0.04

0.00

-0.04

Compnente E-W, X=91, Y=57


Componente N-S, X=91, Y=57

-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

Velocidad de las Componentes N-S y E-W (m/s)

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.

Durante el perodo de medicin de vientos, su velocidad media fue de 2.9 m/s


con direccin del SW, aproximadamente. Bajo estas condiciones, la velocidad
media de la componente N-S de la corriente (promediada en la profundidad)
result ser de 0.04 m/s, esto es 1.4% de la velocidad del viento. La velocidad
media de la componente E-W de la corriente es apenas perceptible.

Los resultados, para los escenarios seleccionados, se presentan en forma


grfica en formato de velocidades sobre la batimetra.

A modo de resumen de lo visto en los grficos anteriores se presenta en la


Tabla 7.2.1. un detalle de los escenarios corridos con comentarios.

ENVIRONMENTAL RESOURCES MANAGEMENT

69

Resultados de la Modelacin Matemtica


Escenario 1. Viento del N. Velocidad del viento : 4 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del N de velocidad 4 m/s. Frecuencia de los vientos del N para
todas las intensidades: 2.1%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

70

Resultados de la Modelacin Matemtica


Escenario 2. Viento del NE. Velocidad del viento : 4 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del NE de velocidad 4 m/s. Frecuencia de los vientos del NE para
todas las intensidades: 1.1%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

71

Resultados de la Modelacin Matemtica


Escenario 3. Viento del E. Velocidad del viento : 4 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del E de velocidad 4 m/s. Frecuencia de los vientos del E para
todas las intensidades: 2.0%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

72

Resultados de la Modelacin Matemtica


Escenario 4. Viento del SE. Velocidad del viento : 4 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del SE de velocidad 4 m/s. Frecuencia de los vientos del SE para
todas las intensidades: 1.9%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

73

Resultados de la Modelacin Matemtica


Escenario 5. Viento del S. Velocidad del viento : 4 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del S de velocidad 4 m/s. Frecuencia de los vientos del S para
todas las intensidades: 13.9%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

74

Resultados de la Modelacin Matemtica


Escenario 6. Viento del S. Velocidad del viento : 8 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del S de velocidad 8 m/s. Frecuencia de los vientos del S para
todas las intensidades: 13.9%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

75

Resultados de la Modelacin Matemtica


Escenario 7. Viento del SW. Velocidad del viento : 4 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del SW de velocidad 4 m/s. Frecuencia de los vientos del SW
para todas las intensidades: 28.0%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

76

Resultados de la Modelacin Matemtica


Escenario 8. Viento del SW. Velocidad del viento : 8 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del SW de velocidad 8 m/s. Frecuencia de los vientos del SW
para todas las intensidades: 28.0%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

77

Resultados de la Modelacin Matemtica


Escenario 9. Viento del SW. Velocidad del viento : 12 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del SW de velocidad 12 m/s. Frecuencia de los vientos del SW
para todas las intensidades: 28.0%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

78

Resultados de la Modelacin Matemtica


Escenario 10. Viento del W. Velocidad del viento : 5 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del W de velocidad 5 m/s. Frecuencia de los vientos del W para
todas las intensidades: 22.3%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

79

Resultados de la Modelacin Matemtica


Escenario 11. Viento del W. Velocidad del viento : 10 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del W de velocidad 10 m/s. Frecuencia de los vientos del W para
todas las intensidades: 22.3%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

80

Resultados de la Modelacin Matemtica


Escenario 12. Viento del NW. Velocidad del viento : 4 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del NW de velocidad 4 m/s. Frecuencia de los vientos del NW
para todas las intensidades: 7.8%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

81

Resultados de la Modelacin Matemtica


Escenario 13. Viento del NW. Velocidad del viento : 8 m/s

Campo de velocidades resultante de la superposicin de la Contracorriente del Per y


corrientes debidas a un viento del NW de velocidad 8 m/s. Frecuencia de los vientos del N
W para todas las intensidades: 7.8%.
Nota: el punto al norte de la Baha de Paracas corresponde al extremo oeste del proyecto.
Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

82

Resultados de la Modelacin Matemtica


Escenario 14. Viento variable sobre registro real. Viento medio del SW con variacin

peridica diurna

Campo de velocidades en trminos de partculas seleccionadas. Inicio de la corrida en el S.


Los vientos tuvieron una velocidad media de 2.9 m/s, durante el perdo de modelacin.

ENVIRONMENTAL RESOURCES MANAGEMENT

83

Resultados de la Modelacin Matemtica


Escenario 15. Contracorriente del Peru. Sin viento ni marea.

Campo de velocidades resultante de la Contracorriente del Per solamente.


Abajo: detalle del campo de corrientes en la Baha de Paracas.

ENVIRONMENTAL RESOURCES MANAGEMENT

84

Resultados de la Modelacin Matemtica


Escenario 16. Marea solamente.

Corrientes bajo condiciones de marea bajante (Arriba) y creciente (Abajo)

ENVIRONMENTAL RESOURCES MANAGEMENT

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

VARIABLE SEGN MEDICIONES

ENVIRONMENTAL RESOURCES MANAGEMENT

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

SIN MAREA NI VIENTO

16

MAREA SOLAMENTE

ENVIRONMENTAL RESOURCES MANAGEMENT

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

El anlisis de los datos disponibles de mareas, vientos y corrientes ha


permitido concluir que:

En la regin estudiada la componente de marea predominante es la


semidiurna principal (M2) siguindole en importancia la diurna (K1), en
una relacin aproximada de 3 a 1, con una amplitud de la componente M2
de 0.22 m. 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). Este
fenmeno tiene influencia sobre las corrientes que, a lo largo de un mes
lunar, pasarn por pocas de baja intensidad (una pleamar por da) y alta
intensidad (dos pleamares por da) , de acuerdo a las caractersticas de la
marea. Las mareas en la zona de Pisco son bien conocidas y confiable su
prediccin por tabla.

Se ha utilizado el Modelo Global de Mareas TPXO.6 o TOPEX que es un


modelo terico que corre sobre batimetra real, calibrado con los datos
satelitales resultantes de la misin TOPEX/Poseidn y estaciones de
marea costera. Las comparaciones de los resultados de este modelo con
datos de las estaciones de marea de El Callao y Pisco es excelente y
recomendable para su uso en otros puntos de la costa peruana, por lo
menos en el segmento costero comprendido entre esos dos lugares. El
sitio de Internet para acceder al modelo se indica en el interior de este
trabajo.

ENVIRONMENTAL RESOURCES MANAGEMENT

88

Estadsticamente, los vientos dominantes son del SW (con frecuencia de


ocurrencia del 28%), W (22.3%), S(13.9%) y NW(7.8%). Las calmas tienen
un porcentaje de ocurrencia de 20.9% y el 7.1% restante se distribuye
entre las otras direcciones en porcentajes semejantes. Las velocidades ms
importantes se han registrado con vientos de SW con hasta unos 15 m/s.
Vientos con velocidades de hasta 8 m/s ocupan el 96.1%.

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).

A partir de la informacin disponible, las corrientes medidas en


profundidad de 15 m en la zona proyectada para la Terminal Martima,
alcanzan velocidades mximas de unos 0.30 m/s (del SW). Las
direcciones ms frecuentes de las corrientes son hacia el NE y N que
suman un porcentaje de ocurrencia algo menor del 50%. Le siguen en
importancia las corrientes hacia el SW y S que suman un 27%,
aproximadamente.

En vista de la dificultad presentada para obtener, por inspeccin, el aporte


de la marea al campo total de corrientes, se realiz el anlisis armnico de
las series de velocidades, de una semana de duracin cada una.

El resultado sugiere que el comportamiento de las corrientes sigue a la


poca de la marea. Es decir, durante mareas diurnas dominantes, se
observa en las corrientes la componente K1 como la ms importante. La
energa de esta componente puede explicar hasta el 11% de la energa
total del registro de corrientes. El resto se asocia al viento principalmente,
dadas las caractersticas de los residuos que resultan de restar de las
velocidades totales las velocidades asociadas a la componente K1

ENVIRONMENTAL RESOURCES MANAGEMENT

89

El otro caso analizado corresponde a una poca de la marea donde no


solamente se encuentra presente, con significacin la componente K1 sino
tambin la componente semidiurna M2. El anlisis armnico de la serie de
corrientes, muestra a ambas componentes de la velocidad (E-W y N-S)
influenciadas por la armnicas M2 y K1 , que pueden explicar el 8.1% de
toda la energa del registro de corrientes en la componente E-W y el 25.9%
en la componente N-S. Dada la escasa longitud de los registros, es
recomendable considerar estas cifras como una estimacin.

El aporte de las mareas al campo de corrientes de la zona del proyecto no


es despreciable pero es ciertamente menor que el aporte del viento.

El campo de corrientes costa afuera es desconocido y solo estimables sus


caractersticas por medio de las mediciones con flotadores de marzo del
2003, de unas 12 h de duracin, entre la Isla San Guilln y la Pennsula de
Paracas, frente al Cerro Talpo y frente a Tambo de Mora, en aguas de 5060 m de profundidad, la literatura cientfica y la utilizacin del modelo
global OCCAM.

La corriente prevaleciente costa afuera es la llamada Contracorriente del


Per, que corre de N a S, con velocidades mximas de unos 0.5 m/s, en la
franja limitada por los 100 y 300 km de las costas peruanas. Ms cercana a
la costa corre con velocidades menores.

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

90

Sobre la Calibracin del Modelo

La calibracin del modelo se realiz considerando datos de alturas de


pleamares y bajamares predichas y medidas en Pisco, mediciones de
corrientes obtenidas con correntgrafo y mediciones de trayectorias de
corrientes medidas con flotadores.

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.

La comparacin con las corrientes medidas con flotadores costa afuera


utilizada para establecer la Contracorriente del Per en el modelo- result
variada segn el lugar de medicin y las condiciones meteorolgicas. Para
el establecimiento de las condiciones de contorno en los lmites del
dominio de modelacin para esta contracorriente, se utiliz el modelo
OCCAM que 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. Una vez calibrado y
corrido el modelo con esta contracorriente, se pudo observar que el MIKE
21 representaba correctamente el campo de corrientes en los otros dos
puntos costa afuera, entre la Isla San Guilln y la Pennsula de Paracas y
frente al Cerro Talpo, por lo que se lo dio por calibrado en cuanto a la
Contracorriente del Per.

La calibracin con las corrientes dentro de la baha, se realiz por medio


de la modelacin en condiciones de vientos y las velocidades medias de
los flotadores, que arrojaron resultados aceptables. Se calcul, adems, las
magnitudes de las corrientes de marea por medio del MIKE 21 y se

ENVIRONMENTAL RESOURCES MANAGEMENT

91

encontraron comparables con las deducidas por medio del anlisis


armnico de los registros de corrientes. Dado que las mareas y las
corrientes de marea requieren de mediciones de por lo menos un mes de
duracin

para

entender

completamente

su

funcionamiento,

es

recomendable obtener esa longitud de registro por medio de mediciones


complementarias y revalidar el modelo. En todo caso, esta situacin no
representa invalidar los resultados del modelo dado que las corrientes de
marea no son dominantes, como se ha demostrado.
Sobre los Resultados de la Modelacin

Los resultados de la modelacin se ordenaron

por

escenarios

seleccionados por las condiciones de viento, principalmente. Se ha tratado


de disponer de un campo de velocidades realista aun cuando la
frecuencia de ocurrencia de los vientos sea muy baja, dado que estudios
posteriores como el desplazamiento de sedimentos en suspensin,
requerirn de todos los escenarios capaces de comprometer la biota, por
ms remotos que sean.

En lneas generales y con variada intesidad, los vientos de N, NE, E, W y


NW generan corrientes hacia el S y tienen una frecuencia de ocurrencia
sumada de 35%, aproximadamente. Sin embargo, los vientos del W
(22.3%) generan corrientes muy dbiles.

Con viento N las corrientes se dirigen al S disminuyendo hacia la Baha


de Paracas. El agua entra por la margen E de la baha y sale bordeando la
Pennsula de Paracas por el W. Dentro de la Baha de Paracas se forma
una celda de circulacin tpica disparada por la batimetra. Las
velocidades en el sitio elegido para el proyecto son de unos 0.1 m/s y la
mitad de ellas dentro de la baha. Con vientos del NE, E y NW el patrn
de circulacin es semejante.

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.

ENVIRONMENTAL RESOURCES MANAGEMENT

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.

Con vientos del SW las corrientes van hacia el N con velocidades en el


sitio que exceden 0.25 m/s fuera de la Baha de Paracas y son menores de
0.2 m/s dentro de la misma.

Las velocidades mximas dentro de la Baha de Paracas, si bien muy


variables por el patrn de circulacin, son consistentemente menores que
las calculadas para el sitio del proyecto. El agua circula a baja velocidaddentro de esta baha, entrando por su margen W y saliendo por su
margen E o viceversa.

Finalmente, se ha generado una importante base de datos de corrientes para la


regin, en base a distintos escenarios que representan la relacin de fuerzas
entre los fenmenos que generan las corrientes caractersticas del lugar.

ENVIRONMENTAL RESOURCES MANAGEMENT

93

También podría gustarte