Advanced
On the Statistical Characteristics of the New Year Wave
On the Statistical Characteristics of the New Year Wave
Journal of Ocean Engineering and Technology. 2013. Feb, 27(1): 102-108
Copyright © 2013, Korean Society of Ocean Engineers
  • Received : January 18, 2013
  • Accepted : February 20, 2013
  • Published : February 28, 2013
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
도영 김

Abstract
In this paper time series wave data, which were measured at the Draupner platform in the North Sea on 1995, are used to investigate statistical characteristics of nonlinear wave. Various statistical properties based on time and frequency domain are examined. The Gram-Chalier distribution fits the probability of wave elevation better than the Gaussian distribution. The skewness of wave profile is 0.393 and the kurtosis is 4.037 when the freak wave is occurred. The nonlinearity of D1520 data is higher than two adjacent wave data. AI index of the New Year Wave is 2.11 and the wave height is 25.6m. The zero crossing wave period of the New Year Wave is 12.5s which is compared to the average zero up-crossing period 11.3s. The significant steepness of wave data is 0.077 when the freak wave was occurred. H1/3/ƞ s does not increases as the kurtosis increases and the values is close to 4. The New Year Wave belongs to highly nonlinear wave data packet but the AI index is within linear focusing range.
Keywords
1. 서 론
해양에서 발생하는 파랑은 여러 가지 종류가 있으며, 이중 바람에 의해서 발생한 파랑은 선박 및 해양구조물의 주요 설계 인자 중의 하나이다. 바람에 의해서 발생된 파랑은 불규칙하고 이러한 불규칙파를 해석하기 위해서 통계적인 해석 방법을 사용하고 있다. 최근에 전형적인 해양파의 최대 파고의 범위를 넘어서는 매우 높은 파의 특성에 대한 연구가 많이 진행되고 있으며, 이런 파도를 Freak wave (또는 Rougue wave)로 부르고 있다. 일반적으로 Freak wave는 파고가 유의파고의 2배가 넘는 큰 파도로 정의되고 있다. 그러나 이 정의에 따른다면 파고의 절대 크기가 크지 않더라도 Freak wave로 분류될 수 있기 때문에 Freak wave의 정의에서 파도의 절대 크기가 10m이상이 되는 조건을 추가로 포함하기도 한다.
이전에는 Freak wave는 아주 드물게 발생하는 것으로 생각되었다. 그러나 해상에서 습득한 파랑자료를 통하여 매우 큰 파도의 발생이 실제로 일어나고 있으며 생각보다 많이 발생하는 것이 발견되었다( Mori, 1997 ; Liu and Pinho, 2004 ). 1991년에는 파고 29m의 파도가 계측되었다( Holliday et al., 2006 ). 1995년에는 25m 크기의 파도가 계측되었다( Haver, 2003 ). 이러한 거대 파도들은 선박이나 해양구조물에 심각한 손상을 초래할 수 있다. 실례로 1982년에 Newfoundland의 Grand bank에서 운영중이던 Mobil사의 시추선을 거대한 파도가 조종석을 강타했고, 이 여파로 시추선이 전복되어 침몰되었고 승선한 승무원 84명 전원이 사망하였고, 1969년과 1999년 사이에 돌연히 발생한 Freak wave에 22척의 초대형선이 전손 또는 심각한 손상을 당하였고 그 결과 542명의 인명손실이 발생되었다( Lawton, 2001 ).
본 논문에서는 1995년 1월 1일 북해에 설치된 Draupner platform에서 계측된 Freak wave가 포함된 파랑자료를 분석하여 각종 통계적인 특성을 검토하여 Freak wave와 같은 거대한 비선형 파랑의 특성을 살펴보고자 한다. 파형의 비선형성의 특성을 살펴보기 위해 비정규분포함수인 Gram-Chalier A분포함수를 사용하였고, 파고의 확률분포에는 Edgeworth-Rayleigh(ER) 분포를 사용하였다.
2. 파랑자료의 개요
Draupner platform은 북해 노르웨이 쪽 수심 70m인 지역에 설치되어 있다. 이 Platform에는 수면 아래 방향으로 레이저가 설치되어 있어 수면과 센서 사이의 순간 거리의 시계열 자료를 계측할 수 있다. 이 논문에서는 1995년 1월 1일 14:20부터 1시간 간격으로 각 20분간 계측된 3개의 파랑자료에 대해서 살펴보기로 한다. 각각의 자료는 Draupner 1420, Draupner 1520 그리고 Draupner 1620 또는 간단히 D1420, D1520 그리고 D 1620으로 표시하기로 한다. 여기서 뒤의 숫자는 계측이 시간된 시간을 나타낸다. 즉 1420은 14시 20분에 계측된 것을 표시한다. 계측자료의 계측주가(Sampling period)는 0.4687573s(= 2.1333Hz)이다. Fig. 1 ~ Fig. 3 에 D1420, D1520, D1620 데이터에 의한 수면형상의 시계열을 표시하였다. Fig. 2 t =280초 부근에 파정이 매우 큰 파도가 포함되어 있는 것을 볼 수 있으며, 이 파도의 크기는 25.6m이다. 이 거대 파도는 1월 1일에 계측되어서 New year wave라고도 부른다.
PPT Slide
Lager Image
Time series of surface elevation for the Draupner 1420
PPT Slide
Lager Image
Time series of surface elevation for the Draupner 1520
PPT Slide
Lager Image
Time series of surface elevation for the Draupner 1620
3. 시계열 파형의 분석
계측된 파형의 시계열자료에서 연속적인 두 개의 영점상향통과(Zero up-crossing) 사이의 구간을 한 개의 파도로 정의하였다. 계측주기(Sampling period)가 0.469 초 이어서 개별 파도 주기의 오차가 ±0.235초가 된다. 개별 파도주기의 오차를 줄이기 위해서 영점을 통과하는 구간에서는 보간법을 사용하여 정확한 영점통과 시점을 적용하였다.
파형 ƞ 의 모멘트는 파형의 확률분포 p 를 이용하여 다음과 같이 정의 된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 μ 는 평균(Mean), σ 2 는 분산(Variance), 𝛾 는 비대칭도(Skewness), κ 는 첨도(Kurtosis)이다. 정규분포인 경우 κ 는 3.0다.
파형의 분포확률을 살펴보기 위해서, Gaussian 분포와, 비정규 확률분포에 사용되는 Gram-Chalier A 급수를 사용하였고, Gram- Chalier A 급수는 다음과 같이 정의 된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
Table 1 에 파형의 모멘트를 표시하였다. 비대칭도와 첨도가 각각 0와 3을 벗어난 값을 나타내고 있어, 정규분포에서 벗어나고 있는 것을 볼 수 있다. D1420의 경우는 비대칭도가 음의 값을 나타내고 있으며, D1520의 경우 비대칭도와 첨도가 가장 큰 값을 가지는 것을 볼 수 있다. 첨도는 | ƞ |가 큰 경우 즉 분포의 꼬리 부분에서 분포확률이 정규분포와 비교해서 얼마나 많거나 부족한 지와 직접적으로 관계되어 진다. 따라서 첨도가 3에서 벗어난 정도는 파정의 최대치나 파저의 최소치의 확률분포가 정규분포로부터 얼마나 벗어난 것과 직접적인 연관이 있게 된다. 계측된 파형의 확률분포를 Fig. 4 ~ Fig. 6 에서 정규분포와 Gram-Chalier A분포와 비교를 하였다. Gram-Chalier A분포가 정규분포보다 최대치의 위치와 크기를 잘 따르는 것을 볼 수 있다. 세 개의 자료 중에서 D1520이 정규분포에서 가장 많이 벗어나는 것을 볼 수 있다.
Various statistical characteristics of surface elevations
PPT Slide
Lager Image
Various statistical characteristics of surface elevations
PPT Slide
Lager Image
Probability distribution of surface elevation for the Draupner 1420
PPT Slide
Lager Image
Probability distribution of surface elevation for the Draupner 1520
PPT Slide
Lager Image
Probability distribution of surface elevation for the Draupner 1620
4. 파고의 분포
Table 2 에는 영점상향통과법(Zero upcrossing method)에 의해서 계산된 파랑의 각종 통계특성을 정리하여 표시를 하였다. 평균영점상향주기는 D1520이 가장 짧아서 포함된 파도의 수가 가장 많은 것을 볼 수 있다. D1520에 포함된 최대파고는 25.6m로 D1420, D1620에 포함된 최대파고보다 약 5m정도 큰 것을 볼 수 있다. 그러나 D1520의 평균파고 및 유의파고의 크기, 제곱근평균제곱 파고(Root mean square, RMS) 파고는 오히려다른 두 개에 비해 작은 것을 볼 수 있다. 최대파고에 대응하는 파도의 주기는 15.6초, 12.5초 그리고 14.3초로 나타났으며, D1520에 포함된 최대파의 주기가 가장 짧은 것으로 나타났다. 파고의 확률분포는 Table 3 에 표시를 하였다. New year wave의 파고는 25.6m, 파주기는 12.5초를 나타내었다. 이것은 Draupner 1520의 평균 영점상향 파주기 11.3초 그리고 파랑스텍트럼에서 얻어지는 평균주기 11.9초와 비슷한 값을 보이는 것을 볼 수 있다. 이것은 Longuet-Higins(1980) , Tucker(1991) 등이 최대파의 주기가 평균주기 근처에서 발생한다고 언급한 결과와 일치를 한다.
Summary of various wave statistical characteristics from time series of surface elevation.
PPT Slide
Lager Image
Summary of various wave statistical characteristics from time series of surface elevation.
Probability of wave height distribution p(H) for D1420, D1520 and D1620.
PPT Slide
Lager Image
Probability of wave height distribution p(H) for D1420, D1520 and D1620.
파고의 확률분포를 Rayleigh 분포로 표시할 수 있으며, 파고의 초과확률(Probability of exceedance)은 다음과 같이 주어진다( Longuet-Higgins, 1952 ).
PPT Slide
Lager Image
비선형분포에서 파고의 초과확률분포는 Edgeworth-Rayleigh(ER) 분포로 표시할 수 있으며, 다음과 같이 표시할 수 있다. ( Mori et al., 2002 ; Mori, 2004 ).
PPT Slide
Lager Image
PPT Slide
Lager Image
5. 파랑스펙트럼
파랑스펙트럼의 대역폭(Bandwidth)을 나타내는 파라미터 v ( Longuet-Higgins, 1980 ), 평균 영점상향 통과주기 Tzs , 유의파첨도(Significant steepness) SS 및 스펙트럼법에서 정의된 유의파고 Hs 는 다음과 같이 정의된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 mn 은 파랑스펙트럼의 n차 모멘트를 나타낸다.
Modified Pierson spectrum (또는 Bretschneider spectrum)은 다음과 같이 정의된다.
PPT Slide
Lager Image
JONSWAP spectrum은 다음과 같이 표시된다( Tucker, 1991 ).
PPT Slide
Lager Image
파형을 고속 푸리에 변환(Fast fourier transform, FFT)하고, 고조파(Harmonics)를 평균하여 최종 주파수 분해능(Frequency resolution)을 0.005Hz로 계산한 파랑스펙트럼의 특성들을 Table 4 에 표시하였다. 그리고 Fig. 7 ~ Fig. 9 에 파랑자료에서 얻은 스펙트럼과 Modified Pierson spectrum 그리고 JONSWAP spectrum과 비교를 하였다. JONSWAP spectrum이 실제 해상자료와 더 일치하는 것을 볼 수 있다.
Various characteristics of wave spectrum
PPT Slide
Lager Image
Various characteristics of wave spectrum
PPT Slide
Lager Image
Comparison of wave Spectrum of D1420
PPT Slide
Lager Image
Comparison of wave Spectrum of D1520
PPT Slide
Lager Image
Comparison of wave Spectrum of D1620
6. 파랑자료의 분석
Freak wave는 일반적으로 이상지수 AI (Abnormality index)를 이용해 다음과 같이 정의 된다( Kalif et al., 2009 ).
PPT Slide
Lager Image
여기서 H 는 최대 파고 H 1/3 은 파랑의 시계열 자료에 포함된 파고 중 상위 1/3의 평균인 유의파고이다. 그러나 종종 스펙트럼법에서 정의된 유의파고를 사용하기도 한다. AI ≥2 만을 Freak wave 발생 기준으로 한다면, 파고의 절대 크기와 관계없이 낮은 해상상태에서도 Freak wave가 정의되고 발생할 수 있다. 따라서 Freak wave의 정의에서 파도의 절대 크기가 10m이상이 되는 조건을 추가로 포함하기도 한다( Stansell, 2004 ). 또 어떤 경우는 Freak wave의 정의에 파정의 크기에 대한 기준을 추가하기도 한다. 유의파고는 영점상향통과법에 의한 H 1/3 과 스펙트럼법 정의에 의한 Hs 가 사용되고 있는데. 일반적으로 H 1/3 Hs 보다 5%~10% 작은 값을 가진다( Forristal, 1978 ; Tucker, 1991 ; Kim, 2007 ). 따라서 파랑자료를 비교할 때 AI 가 어떻게 정의되었는지를 확인해야 한다.
Draupner 파랑자료를 이용해 Rayleigh( Hrms ) 분포함수와 ER분포함수로 계산한 파고의 초과확률분포 결과를 Fig. 10 ~ Fig. 12 에 표시를 하였다. 이 그림에서 파고의 전체적인 초과확률 분포는 Hs/H 1/3 이 1.4까지는 파랑자료와 ER, Rayleigh( Hrms ) 분포가 잘 일치를 하고 있으나 파고가 높아지면 자료에서 벗어나 확률을 과소평가하는 것을 볼 수 있다.
PPT Slide
Lager Image
Probability of wave height exceeding the given value of H for the Draupner 1420
PPT Slide
Lager Image
Probability of wave height exceeding the given value of H for the Draupner 1520
PPT Slide
Lager Image
Probability of wave height exceeding the given value of H for the Draupner 1620
Fig. 13 ~ Fig. 15 에서는 파고와 파주기의 상관관계를 보여주고 있다. 최대파의 주기는 13초~14초 사이의 값을 보이는 것을 볼 수 있다. 특히 Draupner 1520에서는 최대파의 주기는 12초 정도로 분포된 파주기의 평균과 매우 근접한 값을 보이고 있다. 따라서 최대파는 평균파주기와 비슷한 파주기를 가지는 것을 볼 수 있다. 이것은 최대파의 발생은 파랑의 중첩에 의해 발생할 수 있다는 것을 보여준다. Freak wave발생을 수치 시뮬레이션한 결과에 의하면 평균 AI 의 값은 2.11을 가지며 파랑의 중첩에 의해서도 AI 의 값은 최대 2.5~3.0 정도까지 가질 수 있다( Kim, 2010 ). 따라서 Draupner 1520에 포함된 New year wave는 중첩의 원리 범위 안에 있는 이상지수의 값을 가진다.
PPT Slide
Lager Image
Wave period and wave height relation for the Draupner 1420
PPT Slide
Lager Image
Wave period and wave height relation for the Draupner 1520
PPT Slide
Lager Image
Wave period and wave height relation for the Draupner 1620
Freak wave의 발생 원인과 발생 확률에 대한 연구는 여러 연구자에 의해 연구되고 있다. 그러나 발생 확률에 대해서는, 해상상태의 불규칙성과 비선형 현상 등으로 인해서 해석하기가 매우 복잡한 상황이다. 여러 가지 인자 중에서 파형의 첨도 (Kurtosis)의 값이 Freak wave 발생과 밀접한 관계가 있으며 첨도가 큰 경우 Freak wave 발생확률이 높아진다( Mori, 2004 ; Stansell, 2004 ; Shin et al., 2007 ). Table 5 에 파랑자료에서 계산한 첨도(Kurtosis)와 H 1/3의 관계를 표시하였다. Fig. 16 에서는 Table 5 에서 표시한 값과 ER급수에 의한 추정치와 Mori의 실험한 결과와 비교를 하였다. 두 경우 모두 첨도가 커지면 H 1/3 / ƞs 값이 커지는 경향을 나타낸다. 파형이 정규분포인 경우에 첨도의 값은 3을 가지게 되며 H 1/3 / ƞs 값은 4를 나타낸다. 그리고 파고의 분포가 Rayleigh(mo)분포를 따르는 경우는 H 1/3 / ƞs 가 첨도에 무관하며 그 크기는 4로 일정한 값을 가진다. New year wave가 포함된 Draupner 1520에서는 첨도가 커졌지만 오히려 H 1/3 / ƞs 가 작아진 것을 볼 수 있다. 그리고 그 값의 크기는 Mori의 실험결과보다 ER급수를 사용해 예측된 결과와 정성적으로 일치를 한다.
H1/3/ƞsvs.κkurtosis of surface elevation
PPT Slide
Lager Image
H1/3/ƞs vs. κ kurtosis of surface elevation
PPT Slide
Lager Image
H1/3/ƞs vs. kurtosis of Draupner data
Table 6 에서는 유의파첨도와 대역폭 파라미터를 다시 표시하였다. 여기서 D1520은 가장 큰 유의파첨도를 나타내고 있으며, 실제 해상상태에서 발생할 수 있는 최대 유의파첨도인 1/13(=0.077)에 근접하는 값을 보이고 있다.
Significant steepness and spectrum bandwidth
PPT Slide
Lager Image
Significant steepness and spectrum bandwidth
7. 결 언
이 논문에서는 Draupner 1420, 1520, 1620 시계열 파랑자료를 비교해서 New year wave의 통계적 특성을 살펴보았다. Freak wave가 발생한 D1520의 경우 파형의 비대칭도는 0.393와 첨도(Kutosis)는 4.037로 다른 두 개에 비해 가장 비선형성이 많이 나타나고 상대적으로 정규분포에서 가장 많이 벗어나는 것을 볼 수 있다. 파형의 분포는 정규분포를 벗어나며 Gram-Chalier A분포로 표시할 수 있다. New year wave의 이상지수 AI는 2.21이며 파고는 25.6m, 파주기는 12.5초를 나타내었다. 이것은 Draupner 1520의 평균 영점상향 파주기 11.3초 그리고 파랑스펙트럼에서 얻어지는 평균주기 11.9초와 비슷한 값을 가진다. 이것은 최대파의 주기가 평균주기 근처에서 발생하는 것을 볼 수 있다. New year wave가 포함된 D1520은 실제 해상상태에서 발생할 수 있는 최대 유의파첨도에 근접하는 매우 큰 유의파첨도의 값을 보이고 있다. JONSWAP spectrum이 Modified Pierson spectrum보다 실제 해상자료와 더 일치하는 것을 볼 수 있다. New year wave는 극한 해상상태에서 발생한 거대한 freak wave이지만 중첩의 원리 범위 안에 있는 이상지수의 값을 가지고 있다.
Acknowledgements
본 연구는 2010년 홍익대학교 학술연구진흥비의 지원으로 수행된 연구결과 중 일부입니다.The author would like to thank Dr. Sverre Haver of Statoil for kindly providing the Draupner wave data.
References
Forristall G.Z. 1978 On the Statical Dstribution of Wave Heights in a Storm J. Geophys. Res. 83 2353 - 2358    DOI : 10.1029/JC083iC05p02353
Haver S. 2003 Freak Wave Event at Draupner Jacket January 1 1995 PTT-KU-MA Staoil
Holliday N.P. , Yelland M.Y. , Pascal R. 2006 Were Extreme Waves in the Rockall Trough the Largest Ever Recorded Geophys. Res. Lett. 33 L05613 -    DOI : 10.1029/2005GL025238
Kalif C. , Pleinovsky E. , Slunyaev A. 2009 Rogue Waves in the Ocean Springer
Kim D.Y. 2007 On Individual Wave Height Distribution of Ocean Waves Journal of Ocean Engineering and Technology 21 (2) 6 - 11
Kim D.Y. 2011 On the Characteristics of Freak Wave Occurrence J. Korean Soc. for Marine Environmental Eng. 14 (2) 138 - 145    DOI : 10.7846/JKOSMEE.2011.14.2.138
Liu P.C. , Pinho U.F. 2004 Freak Waves-More Frequent than Rare! Annales Geophysics 22 1839 - 1842    DOI : 10.5194/angeo-22-1839-2004
Lawton G. 2001 Monsters of the Deep (The Perfect Wave) New Scientist 170 (2297) 28 - 32
Longuet-Higgins M.S. 1952 On the Statistical Distribution of the Heights of Sea Waves J. Marine Res. 11 245 - 266
Longuet-Higgins M.S. 1980 On the Distribution of the Heights of Sea Waves: Some Effects of Non-linearity and Finite Band-width J. Geophy. Res. 85 (C3) 1519 - 1523    DOI : 10.1029/JC085iC03p01519
Mori N. 1997 Occurrence Properties of Giant Freak Waves in Sea Area around Japan J. waterway, port, coastal Engineering (July/Aug) 209 - 213
Mori N. 2004 Occurrence Probability of a Freak Wave in a Nonlinear Wave Field Ocean Engineering 31 165 - 175    DOI : 10.1016/S0029-8018(03)00119-7
Mori N. , Paul C.L. , Yasuda T. 2002 Analysis of Freak Wave Measurements in the Sea of Japan Ocean Engineering 29 1399 - 1414    DOI : 10.1016/S0029-8018(01)00073-7
Shin S.H. , Hong K.Y. , Moon J.S. 2007 Wave Data Analysis for Investigation of Freak Wave Characteristics J. Lorean Navigation and Port Research 31 (6) 471 - 478    DOI : 10.5394/KINPR.2007.31.6.471
Stansell P. 2004 Distributions of Extreme Wave, Crest and Trough Heights Measured in the North Sea Ocean Engineering 32 1015 - 1036
Tucker M.J. 1991 Waves in Ocean Engineering Ellis Horwood