Tema - 4 - Caracterización Del Oleaje PDF
Tema - 4 - Caracterización Del Oleaje PDF
Tema - 4 - Caracterización Del Oleaje PDF
de Ingeniería
Puertos y Costas I
1 Introducción
1 Introducción
Necesidad de caracterización
El diseño de actuaciones así como la realización de diversidad de estudios en una
determinada zona litoral requieren la definición de las características del oleaje real
existente en dicha zona costera o en una ubicación específica dentro de ella.
Oleaje regular
El oleaje regular o monocromático (formado por una onda con un único período o
frecuencia) puede ser caracterizado en una ubicación de interés específica por medio de
sólo tres variables fundamentales: altura de ola, H, período, T y ángulo de incidencia, .
Oleaje real
El oleaje real no se corresponde con un oleaje regular, caracterizándose por presentar
cierta variabilidad en los parámetros mencionados, denominándose de este modo oleaje
irregular.
A largo plazo: análisis de las condiciones de oleaje para un período temporal extenso
(e.g., varios años), recurriendo para ello al empleo de distribuciones estadísticas.
Para llevar a cabo cualquiera estos análisis es preciso disponer de datos de partida fiables!!
1 Introducción
Estado de mar
Situaciones a corto plazo durante las cuales las características del oleaje son estables
estadísticamente y energéticamente (duraciones del orden de decenas de minutos)
Estado de mar
Situaciones a corto plazo durante las cuales las características del oleaje son estables
estadísticamente y energéticamente (duraciones del orden de las decenas de minutos).
E1 E2 E3 E4 E5
Los estados de mar se pueden caracterizar por unos pocos parámetros: H, T, , etc.
Técnicas caracterización corto plazo estados: i) análisis en el tiempo ii) análisis espectral
Fundamentos
Estudio del oleaje asumiendo que éste se encuentra formado por una sucesión aleatoria
de ondas diferentes que son analizadas de forma individual.
Ondas irregulares
No son olas
H T
H1 T1
H2 T2
H1 > H2 > H3 >…> HN
H3 T3
… …
HN TN
Altura de ola cuadrática media, Hrms, que se calcula como la raíz cuadrada de la media
de los cuadrados de las altura de ola identificadas:
0.5
1 N
2
H rms
N H i
i 1
1 N /n 1 N /n
N /n
H 1 /n H i T 1/n Ti
i 1 N /n i 1
En el caso de n =3, se analiza el tercio de las olas con altura más grande,
denominándose los parámetros obtenidos como altura de ola significante, Hs , y
período significante, Ts :
1 N /3 1 N /3
H s H 1/3
N / 3 i 1
H i T s T 1/3
N / 3
Ti
i 1
Altura de ola media del n % de las ondas más grandes, Hn10-3, (e.g., altura de ola
media superada por el 5% de las olas más grandes, H0.05).
Altura de ola máxima, Hmax , que se corresponde con la mayor Hi del registro, y el
período de onda máximo, Tmax , que se obtiene como el período de la onda con Hmax .
2H H
2
H H
2
p H 2 exp p H exp
H rms H rms 2 Hm
2
4 Hm
Es común representar la densidad de probabilidad para la altura de ola normalizada
X=H /Hrms o X =H /Hm , de modo que se tenga una distribución más suavizada, e.g.,:
2H H 2
p H / H rms p H H rms exp
H rms H rms
Una forma más intuitiva para determinar la función de densidad de la altura de ola
normalizada es definir la variable X = H /Hrms o X =H /Hm de modo que, eg., X = H /Hrms:
2X X 2
p X exp
X rms
2
X rms
2H H 2 2H H
2
p H 2 exp p H dH 2 exp dH
H rms H rms H rms H rms
Para determinar la FDA la expresión anterior debe ser integrada. De este modo, la
probabilidad de que la altura de ola de una onda individual, H, no supere/supere un
determinado valor, H ′, viene dada por:
H' H ' 2
P (H H ' ) p (H )dH 1 exp
0
H rms
H' H ' 2
P (H H ' ) 1 p (H )dH exp
0
H rms
X' H ' 2
P (X X ' ) 1 p (X )dX exp
0
H rms
p (X )
P (X )
P (X ) P (X X ' )
X X
N
H ' HH ’
rms ln
n
La altura de ola 1/n, H1/n , esto es la altura de ola media de las N /n olas con mayor
altura se determina como:
H 1/n p (H )dH Hp (H )dH
H' H'
H 1/
H
Hp (H )dH
n Hp (H )dH
n
'
H'
H p (H )dH
'
donde H′ se corresponde con la altura de ola excedida por las 1/n olas más grandes. Si las
olas consideradas se corresponden con el tercio de altura mayor, la altura de ola obtenida
se corresponde con Hs de modo que:
H p (H )dH 1 / n 1 / 3
'
Hm
Hp (H )dH
0
p (H )dH
0
De este modo, se pueden obtener las diferentes relaciones entre las variables
fundamentales:
n H1 /n /Hrms H1 /n /Hm
Las más importantes:
H m 0.63H s
1
H rms H s 0.707H s
2
H 1/10 1.27H s
H max 1
1/2
1/2
Hmax más probable para número de olas N: ln N (2.887ln N )
H 1/3 1.416
(Hmax /H1/3 = 1,5 - 2)
Puertos y Costas I Tema4 Caracterización del oleaje 27
2 Caracterización a corto plazo
Análisis en el tiempo
Tm (0.7 1.1)T s
T 1/ 10 (0.6 1.3)T s
CCL
¿Por qué crees que las noticias que se muestran para las mismas fechas difieren
tanto?
CCT
Calcula:
CCT y CCL
Se acompaña un registro de elevaciones del nivel del mar de un minuto de duración
medido en una boya en aguas profundas. Aplicando un análisis en el dominio del tiempo
mediante el criterio de paso ascendente por cero:
a) Estima los siguientes parámetros descriptores del estado de mar: Hmax , Tmax , Hm , Tm ,
Hrms , H1 /3 y T1 /3 .
b) Representa en el gráfico que se acompaña la distribución de probabilidad de las olas
identificadas en el registro. ¿Por qué crees que difiere tanto de la distribución de Rayleigh?
CCT y CCL
Fundamentos
Espectro de la luz (Newton, 1671): Cómo varía la intensidad de la luz (su color) en función
de la longitud de onda.
Fundamentos
El análisis espectral del oleaje asume Ola regular 1 η1 (a1 , ω1, δ1 , θ1 )
que un oleaje irregular está formado por
+ +
la superposición de numerosos
Ola regular 2 η2 (a2 , ω2 , δ2 , θ2 )
componentes, siendo cada componente
+ +
una onda regular con una determinada
amplitud y frecuencia. Ola regular 3 η3 (a3 , ω3 , δ3 , θ3 )
= =
N
(t ) ai cos (i t i )
i 1
N
(t ) ai cos (2 f i t i )
i 1
t
siendo ai la amplitud, i la frecuencia angular, fi la frecuencia lineal y i la fase
correspondientes a la componente i de cada uno de los N componentes considerados.
Fundamentos
Ejemplo obtención oleaje irregular superposición 3 olas regulares
η(t)
a? ω? δ?
CCL
Representa el espectro de
amplitudes del oleaje presentado
formado por tres componentes
individuales a través de una
gráfica a (m) Vs f (Hz).
N
Varianza x : x 2
Espectros del oleaje. Espectro de varianza i 1 CP Física
1 2 a2
a (f )cos (2 ft )d (2 ft )
2 2 2 2
(t ) acos (2 ft )
2 0 2
ai2
i i
2 2
2
Determinación directa a partir del espectro de amplitudes.
1 1 ai2
E i gH i g 2ai g
2
2
g i2 g 2i
8 8 2
a 2 (f ) / 2 E v (f ) S (f )
E v (f )
f
Ev (f )[m2Hz-1]
a 2/2 [m2]
f [Hz] f [Hz]
Ev (f )[m2Hz-1]
Ei ai2 / 2
E f (f ) g gE v (f )
f f
Energía (ud superf.) en cada banda centrada en f
con ancho, f :
E f ,f gE v (f )f f [Hz]
i
Uso de forma indistinta como E (f )
0
i 1
Momentos espectrales
Se define como momento espectral de orden n como: mn f n E v (f )df
0
m0 E v (f )df
0
Ev (f )[m2Hz-1]
Representa, el área total del espectro, permitiendo así
determinar la energía de un determinado estado de mar:
E t g E v (f )df E t g E v (f )df gm0
0 0
f [Hz]
Hm0
H'
Hp (H )dH
3 ' Hp (H )dH 4.04... m0 4 m0
H
H p (H )dH
'
CP Oleaje regular N
1 1 i
H 2
1
Ei gH i2 Et g i 1
gH rms
2
8 8 N 8
1 1 1
gH m2 0 gH rms
2
H rms Hm0
16 8 2
1 m0
• Período correspondiente a la frecuencia media, Tm01 : Tm 01
fm m1
• Período energético, Te , que se corresponde con el período de una onda regular con
la misma energía que la onda irregular representada en un determinado espectro:
m 1
Te
m0
Puertos y Costas I Tema4 Caracterización del oleaje 45
2 Caracterización a corto plazo
Análisis espectral
El espectro direccional
El oleaje real se caracteriza porque cada una de sus componentes posee una dirección de
propagación determinada. Para su caracterización espectral, es preciso recurrir al
denominado espectro direccional que considera la distribución de la energía de modo
conjunto para las diferentes frecuencias y direcciones:
E v (f , ) E v (f )D (f , )
donde D es la función de distribución direccional: D (f , )d 1
E v (f )
Espectros paramétricos
Los espectros de diferentes estados de mar presentan parecidos muy notables. Esta
similitud da lugar a que sea posible definir de modo aproximado un espectro de un
estado de mar con base a únicamente un número limitado de parámetros.
Pierson-Moskotwitz
Ev (f )[m2Hz-1]
f [Hz]
Puertos y Costas I Tema4 Caracterización del oleaje 47
2 Caracterización a corto plazo
Análisis espectral
Espectros paramétricos
(f f )2
JONSWAP:
4
exp
p)
5 f 2 f 2
2
E vJ (f ) g (2 ) f exp
2 4 5
p
4 fp
a 0.07paraf f p
b 0.09paraf f p
siendo un parámetro de escala dependiente de las condiciones de generación del
oleaje y el factor de pico (peak enhancement factor). Espectro estrictamente válido
para oleajes no totalmente desarrollados*.
Espectros paramétricos
Espectro JONSWAP de un estado de mar con diferentes valor de factor de pico
Ev (f )[m2Hz-1]
f [Hz]
2
cos 2
( ) para
D ( ) 2
0 para
2
CCT y CCL
Determina los principales
parámetros espectrales (Hm0 , Tp y
Tm ) del espectro de densidad de
varianza de un estado de mar en
una localización en la costa de
Galicia (Figura). Emplea para ello
métodos gráficos aproximados.
1 Introducción
Definición
La caracterización del oleaje a largo plazo proporciona una descripción de la ocurrencia de
los descriptores a corto plazo (parámetros estados de mar) para varios años, décadas, etc.
Procedimiento
Se analizan de forma individual o conjunta la ocurrencia de la altura de ola de los estados
de mar, por ser el parámetro más relevante, el período y/o la dirección del oleaje.
Análisis preliminar
Análisis de las series temporales de los datos disponibles en una determinada localización
a través de representaciones gráficas o tablas: curvas de estados de mar, tablas e
histogramas de distribución de altura de ola y período, rosas de oleaje, etc., que cubren
usualmente períodos concretos de tiempo, e.g. anuales, estacionales, mensuales, etc.
Oleaje de diseño
Oleaje referente a condiciones extremales para el cual se diseña una determinada
actuación litoral. La altura de ola asociada a dicho oleaje, o altura de ola de diseño, Hd ,
puede ser expresada a través de diversas variables (e.g. Hmax , Hs , H1/10 , etc.).
Estudios en el litoral
La caracterización del régimen extremal del oleaje puede ser de gran interés para la
realización de estudios relacionados con la aparición de eventos poco frecuentes.
Selección de temporales
Para una adecuada definición del régimen extremal únicamente se deben considerar los
valores de riesgo o extremos, que normalmente están asociados a temporales.
El Método de los Máximos Anuales que consiste en analizar únicamente los máximos
anuales de los registros. El problema de este método reside en que de forma general
existe un número limitado de años de registros en una determinada ubicación, por lo que
este método no es válido en la mayoría de las ocasiones.
no independientes
Sin embargo, las excedencias no son un valor fijo sino aleatorio, por lo que es preciso
combinar la GPD con la distribución del número aleatorio de excedencias que se
asume que sigue una Distribución de Poisson.
exp 1 k sik 0
P (X x )
x
exp 1 sik 0
donde:
x B
P (X x ) exp exp
A
Distribución de Weibull (Fisher-Tippett tipo III):
x B C
P (X x ) 1 exp
A
Donde:
x C
P (X x ) 1 exp
A
Período de retorno
Período de retorno, Tr : duración promedio durante la cual los eventos extremales exceden
un determinado valor, i.e., el número de años que en término medio transcurre entre dos
temporales que superan un determinado valor de altura de ola. Dicha altura de ola se
denomina de este modo altura de ola de retorno, HTr (o Hr). La relación entre la
probabilidad de superación de HTr y Tr viene dada por la expresión:
1 1
P (H HTr )¿? P (H HTr )¿? 1
donde: ¿? Tr ¿? Tr
Período de retorno
Vida útil: tiempo durante el cual se pretende que dicha actuación sea capaz de
realizar la función para la cual se diseña.
Se puede obtener fácilmente una relación entre los tres términos (Tr , L , E ) considerando:
Riesgo admisible (E )
Probabilidad de que un evento extremo se presente
por lo menos una vez a lo largo de la vida útil.
L
1
E 11
T r
Distribución de Gumbel:
HTr B 1 ¿? HTr B
P (H HTr ) exp exp ln ln
A P (H H ¿? )
Tr A
1 Tr
HTr B A ln ln B A ln ln
P (H H )
Tr
rT 1
Distribución de Weibull:
H B C ¿? H B C
P (H HTr ) exp Tr ln P (H HTr ) ln exp
Tr
¿?
A A
1/C 1/C
HTr B A ln P (H HTr ) B A ln(Tr ) (dist. triparamétrica)
1/C 1/C
HTr A ln P (H HTr ) A ln(Tr ) (dist. biparamétrica)
E (T ) aH c
donde E (T ) es el valor estimado del período para un determinado temporal con altura de
ola H. Es usual realizar dicho análisis en términos de Tp y Hs de modo que:
E (T p ) aH sc
realizándose el ajuste de esta expresión para los datos de temporales identificados.
Ajuste de distribuciones.
Método de máxima verosimilitud (Maximum likelihood estimation, MLE)
En el LSE se realizan estimaciones de los mínimos cuadrados calculadas a partir del ajuste
de una recta de regresión al conjunto de datos poseyendo ésta un error de mínimos
cuadrados (suma mínima de las desviaciones elevada al cuadrado).
y mx b
Esta ecuación se ajusta para que represente lo mejor posible los datos de la muestra
dados por los pares de valores [x, P (H ≤ x )] o [x, P (H x )]. A partir de los valores de m
y b obtenidos del ajuste se determinan los parámetros A y B de la función de distribución.
Para ello, los temporales independientes disponibles en los registros son identificados a
través de uno de los métodos descritos, usualmente el método POT, y las probabilidades
de no excedencia y excedencia de un determinado valor de altura de ola obtenidas a
través de la asignación de cada temporal a intervalos con una resolución específica
dependiente del número de temporales disponible.
Puertos y Costas I Tema4 Caracterización del oleaje 71
3 Caracterización a largo plazo
Régimen extremal
x HTr
2. Se toman como variables reducidas: ¿?
y ln P (¿?
H HTr )
1/C
4. Por tanto y ¿?
1
x
B
A¿? A
y mx¿? b
5. Si se comparan las ecuaciones: 1 ¿? B
1 ¿? B m b
y x A ¿? A
A ¿? A
Nota:
A = β; B = α ; C = γ
2º. Cálculo de Tr
3º. Se obtiene, mediante régimen extremal, la altura de ola de retorno (e.g. altura
significante) correspondiente a ese Tr (determinación de Hd a través de Rayleigh)
Ajuste de distribuciones método mínimos cuadrados (Seminarios Interactivos*)
4º. El T (e.g. Tp) asociado a esa altura (e.g. Hs) se puede obtener mediante fórmula
ajustada
¿DATOS EN EMPLAZAMIENTO?
En este caso, las condiciones de oleaje obtenidas a través del procedimiento expuesto
deben ser propagadas hasta la ubicación o zona de interés, considerando los procesos
litorales presentados (Transformación de ondas) haciendo uso de técnicas numéricas.
CCT
En un determinado punto del mar, el régimen extremal viene dado por la
función de distribución de Weibull triparamértrica con parámetros: A = 0.048; B
= 2.992; C = 0.330. Dichos parámetros han sido obtenidos a través del método
POT a partir de una muestra con λ = 1.08 (temporales al año). Determina:
El oleaje de diseño en dicha ubicación para una actuación litoral con un riesgo
admisible E =0.3 y una vida útil L = 25 años.
Solución: Hd = 7,1 m
Tipos de análisis
Análisis por medio de gráficas (todos los estados de mar durante período de tiempo).
Distribuciones
Distribución de Weibull triparamétrica o biparamétrica.
Si se analizan todos los datos sin tener en consideración la direccionalidad del oleaje
se obtiene el denominado régimen escalar del oleaje; en caso de emplear los datos
correspondientes a un determinado sector direccional, se obtiene el régimen
direccional referente a dicho sector. La probabilidad de no excedencia de una altura
de ola condicionada a que un oleaje provenga de un determinado sector se obtiene
como: P (H x ) P (A B )
P (H x / ) P (A / B )
P ( ) P (B ) CP Matemáticas
Aplicaciones
El estudio del régimen medio es de particular interés en ciertas aplicaciones en Ingeniería
de Costas y Puertos.
CCT
Dada una ubicación costera en donde se dispone de una caracterización del régimen
medio a través de una distribución de Weibull triparamétrica con parámetros: A=1.89;
B=0.58; C =1.31, estima:
a) Probabilidad de que Hs sea menor a 5 m
b) Hs12
En este sentido, existe un estado límite normal de operación a partir del cual la
explotación y operatividad se ven afectadas por causas externas a la propia
infraestructura, no existiendo existir en ningún caso daño estructural en ella. Este
estado está determinado fundamentalmente por diferentes agentes y procesos
litorales, tales como el viento, las corrientes, el oleaje, etc., siendo de forma usual el
oleaje el factor más limitante.
http://www.youtube.com/watch?v=QyeA8eqF14c
CCT
Se proporcionan (Tabla) los ajustes obtenidos para un distribución Weibull triparamétrica de los
regímenes medios direccionales en aguas profundas frente a la costa de Cartagena para las tres
principales direcciones del oleaje que afectan a la dársena de Escombreras. Del mismo modo se
proporciona la frecuencia de presentación de dichas direcciones.
AP A B C Frec. (%)
SSW 0.61 0.36 1.4 10
SW 0.98 0.34 1.6 18.6
WSW 0.99 0.35 1.62 7.7
Para poder transferir los oleajes desde aguas profundas hasta la bocana de la dársena se ha realizado
un estudio de transformación del oleaje del cual se obtiene:
Estima la operatividad de la dársena en horas anuales sabiendo que la altura significante límite
operativa es de Hs=0.5 m.
En este caso, las condiciones de oleaje obtenidas a través del procedimiento expuesto
deben ser propagadas hasta la ubicación de interés, considerando los procesos litorales
presentados (Transformación del oleaje), y haciendo uso de técnicas numéricas.
1 Introducción
Tipos de modelos
De modo general los modelos de transformación del oleaje se pueden clasificar en
modelos espectrales y no espectrales.
Modelos espectrales
Modelos no espectrales:
• De rayos
• De Boussinesq
• No hidrostáticos
Los modelos más empleados son los modelos espectrales. No obstante en muchos casos
su aplicación no es recomendable (e.g. presencia de obstáculos).
1 Introducción
https://www.youtube.com/watch?v=mTIspbiHers