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.
해양에서 발생하는 파랑은 여러 가지 종류가 있으며, 이중 바람에 의해서 발생한 파랑은 선박 및 해양구조물의 주요 설계 인자 중의 하나이다. 바람에 의해서 발생된 파랑은 불규칙하고 이러한 불규칙파를 해석하기 위해서 통계적인 해석 방법을 사용하고 있다. 최근에 전형적인 해양파의 최대 파고의 범위를 넘어서는 매우 높은 파의 특성에 대한 연구가 많이 진행되고 있으며, 이런 파도를 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) 분포를 사용하였다.
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에
계측된 파형의 시계열자료에서 연속적인 두 개의 영점상향통과(Zero up-crossing) 사이의 구간을 한 개의 파도로 정의하였다. 계측주기(Sampling period)가 0.469 초 이어서 개별 파도 주기의 오차가 ±0.235초가 된다. 개별 파도주기의 오차를 줄이기 위해서 영점을 통과하는 구간에서는 보간법을 사용하여 정확한 영점통과 시점을 적용하였다.
파형
여기서
파형의 분포확률을 살펴보기 위해서, Gaussian 분포와, 비정규 확률분포에 사용되는 Gram-Chalier A 급수를 사용하였고, Gram- Chalier A 급수는 다음과 같이 정의 된다.
Table 1에 파형의 모멘트를 표시하였다. 비대칭도와 첨도가 각각 0와 3을 벗어난 값을 나타내고 있어, 정규분포에서 벗어나고 있는 것을 볼 수 있다. D1420의 경우는 비대칭도가 음의 값을 나타내고 있으며, D1520의 경우 비대칭도와 첨도가 가장 큰 값을 가지는 것을 볼 수 있다. 첨도는 |
[Table 1] Various statistical characteristics of surface elevations
Various statistical characteristics of surface elevations
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) 등이 최대파의 주기가 평균주기 근처에서 발생한다고 언급한 결과와 일치를 한다.
[Table 2] Summary of various wave statistical characteristics from time series of surface elevation.
Summary of various wave statistical characteristics from time series of surface elevation.
[Table 3] Probability of wave height distribution p(H) for D1420, D1520 and D1620.
Probability of wave height distribution p(H) for D1420, D1520 and D1620.
파고의 확률분포를 Rayleigh 분포로 표시할 수 있으며, 파고의 초과확률(Probability of exceedance)은 다음과 같이 주어진다(Longuet-Higgins, 1952).
비선형분포에서 파고의 초과확률분포는 Edgeworth-Rayleigh(ER) 분포로 표시할 수 있으며, 다음과 같이 표시할 수 있다. (Mori et al., 2002; Mori, 2004).
파랑스펙트럼의 대역폭(Bandwidth)을 나타내는 파라미터
여기서
Modified Pierson spectrum (또는 Bretschneider spectrum)은 다음과 같이 정의된다.
JONSWAP spectrum은 다음과 같이 표시된다(Tucker, 1991).
파형을 고속 푸리에 변환(Fast fourier transform, FFT)하고, 고조파(Harmonics)를 평균하여 최종 주파수 분해능(Frequency resolution)을 0.005Hz로 계산한 파랑스펙트럼의 특성들을 Table 4에 표시하였다. 그리고 Fig. 7~Fig. 9에 파랑자료에서 얻은 스펙트럼과 Modified Pierson spectrum 그리고 JONSWAP spectrum과 비교를 하였다. JONSWAP spectrum이 실제 해상자료와 더 일치하는 것을 볼 수 있다.
[Table 4] Various characteristics of wave spectrum
Various characteristics of wave spectrum
Freak wave는 일반적으로 이상지수
여기서
Draupner 파랑자료를 이용해 Rayleigh(
Fig. 13~Fig. 15에서는 파고와 파주기의 상관관계를 보여주고 있다. 최대파의 주기는 13초~14초 사이의 값을 보이는 것을 볼 수 있다. 특히 Draupner 1520에서는 최대파의 주기는 12초 정도로 분포된 파주기의 평균과 매우 근접한 값을 보이고 있다. 따라서 최대파는 평균파주기와 비슷한 파주기를 가지는 것을 볼 수 있다. 이것은 최대파의 발생은 파랑의 중첩에 의해 발생할 수 있다는 것을 보여준다. Freak wave발생을 수치 시뮬레이션한 결과에 의하면 평균
Freak wave의 발생 원인과 발생 확률에 대한 연구는 여러 연구자에 의해 연구되고 있다. 그러나 발생 확률에 대해서는, 해상상태의 불규칙성과 비선형 현상 등으로 인해서 해석하기가 매우 복잡한 상황이다. 여러 가지 인자 중에서 파형의 첨도 (Kurtosis)의 값이 Freak wave 발생과 밀접한 관계가 있으며 첨도가 큰 경우 Freak wave 발생확률이 높아진다(Mori, 2004; Stansell, 2004; Shin et al., 2007). Table 5에 파랑자료에서 계산한 첨도(Kurtosis)와
[Table 5] H1/3/?s vs. κ kurtosis of surface elevation
H1/3/?s vs. κ kurtosis of surface elevation
Table 6에서는 유의파첨도와 대역폭 파라미터를 다시 표시하였다. 여기서 D1520은 가장 큰 유의파첨도를 나타내고 있으며, 실제 해상상태에서 발생할 수 있는 최대 유의파첨도인 1/13(=0.077)에 근접하는 값을 보이고 있다.
[Table 6] Significant steepness and spectrum bandwidth
Significant steepness and spectrum bandwidth
이 논문에서는 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이지만 중첩의 원리 범위 안에 있는 이상지수의 값을 가지고 있다.