隨機(jī)波列中極端波浪波形和波高分布的數(shù)值試驗(yàn)研究

夏維達(dá), 馬玉祥, 董國海

(大連理工大學(xué) 海岸與近海工程國家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116024)

極端波浪是波高顯著高于背景波浪的大波,由極端波浪引發(fā)的海洋事故每年都會(huì)造成巨大的經(jīng)濟(jì)損失和嚴(yán)重的人員傷亡。極端波浪具有波高大、非線性強(qiáng)等特點(diǎn),并且出現(xiàn)之前無明顯征兆,因此會(huì)對(duì)海洋結(jié)構(gòu)物和船舶的航行安全造成巨大威脅[1-3]。對(duì)其的研究雖然已有幾十年,但是目前對(duì)極端波浪生成和特性的研究仍然是學(xué)界的研究熱點(diǎn)[4-9];波浪波高分布和極端波浪波形的研究對(duì)工程設(shè)計(jì)也有重要的意義。在一般的海洋結(jié)構(gòu)物設(shè)計(jì)中,通常采用斯托克斯波為設(shè)計(jì)波浪,波高根據(jù)長期的水文數(shù)據(jù)確定;但是實(shí)際的海洋環(huán)境復(fù)雜,波浪為隨機(jī)過程,所以為了更精確地計(jì)算結(jié)構(gòu)物響應(yīng),通常需要長時(shí)間的模擬。Boccotti[10]和Tromans等[11]經(jīng)過嚴(yán)密的理論推導(dǎo),給出了隨機(jī)波列中最可能出現(xiàn)的極端波波群的波面表達(dá)式,也叫新波模型(new wave model),為海洋結(jié)構(gòu)物設(shè)計(jì)計(jì)算提供了新的思路。該模型的有效性也得到了驗(yàn)證[12-15]。之后Walker等[12]在新波模型的基礎(chǔ)上加入了非線性項(xiàng),并驗(yàn)證表明非線性修正后的模型能夠較好地符合海洋中實(shí)測(cè)極端波浪的波形。Santo等[14]也研究了墨西哥灣颶風(fēng)下觀測(cè)到的極端波浪,并與新波模型進(jìn)行了對(duì)比,發(fā)現(xiàn)實(shí)測(cè)的波浪具有更加明顯的垂直不對(duì)稱性并且具有明顯的三維傳播特征。Whittaker等[13]后來也通過實(shí)測(cè)數(shù)據(jù)驗(yàn)證了包含二階非線性的新波模型在淺水條件的適用性。極端波浪的波形特征在波浪與結(jié)構(gòu)物相互作用的研究中也至關(guān)重要。Slunyaev等[16]運(yùn)用數(shù)值方法,模擬了隨機(jī)波浪的傳播,將波浪形態(tài)的分為前峰波、后峰波、前谷波和后谷波4類,發(fā)現(xiàn)大部分極端波浪的波峰后的波谷相對(duì)峰前的波谷要更深。Sergeeva等[17]也發(fā)現(xiàn)了類似的結(jié)果。Ma等[18]通過模型實(shí)驗(yàn)研究了逆流對(duì)規(guī)則波列中由調(diào)制不穩(wěn)定性導(dǎo)致的畸形波的波陡、波形的不對(duì)稱性的影響,發(fā)現(xiàn)逆流的存在會(huì)改變畸形波的不對(duì)稱性。

隨機(jī)波浪的波高分布在工程設(shè)計(jì)中具有重要的作用。Longuet-Higgins[19]對(duì)有效波高、平均波高和極大波高的關(guān)系做了系統(tǒng)性的研究,并給出了窄譜下的最大波高分布。Longuet-Higgins[20]接著在線性波浪模型下推導(dǎo)出了波浪波高服從Rayleigh分布;Kinsman[21]在現(xiàn)場觀測(cè)中發(fā)現(xiàn)偏態(tài)會(huì)導(dǎo)致波面分布不再服從正態(tài)分布,于是他擬合出一種更符合實(shí)測(cè)數(shù)據(jù)的波面分布。Tayfun等[22]和Forristall[23]在考慮了二階非線性后各自提出了二階的波峰分布。但是二階Tayfun分布和Forristall分布只考慮了約束波對(duì)波峰分布的影響,并沒有考慮自由波共振作用的影響。Mori等[6,24]在窄譜的假定下,將波面高程概率密度函數(shù)用Edgeworth級(jí)數(shù)展開后,提出了MER分布。Tayfun等[25]同樣考慮了四階近似共振對(duì)波峰分布的影響,推導(dǎo)出了三階Tayfun分布。這些波高分布模型都不是完美或者適用于所有工況的,需要根據(jù)實(shí)際情況和實(shí)際觀測(cè)數(shù)據(jù)選擇合適的波高分布模型。本文通過數(shù)值模擬多組不同工況下隨機(jī)波的演化,對(duì)極端波浪的平均波形做了研究,驗(yàn)證了新波模型的有效性;最后對(duì)比了實(shí)驗(yàn)數(shù)據(jù)的波高,波峰分布與各種理論分布(Rayleigh分布、MER分布、Forristall分布和Tayfun分布等),探討各種分布的適用性。本文中極端波浪的標(biāo)準(zhǔn)為波峰ηc/Hs>1.25的波浪或波高H/Hs>2.0的波浪。

1 數(shù)值模型及波浪參數(shù)

1.1 數(shù)值模型介紹

本文數(shù)值模擬采用的模型是高階譜模型,由Dommermuth等[26]和West等[27]提出,能夠滿足完全非線性自由表面條件。該模型可以模擬任意階數(shù)的波浪非線性傳播和相互作用。高階譜方法的優(yōu)點(diǎn)是收斂性好、計(jì)算速度快、精度高;并且不受水深、譜寬等的限制。對(duì)于不可壓縮的無粘無旋理想流體,引入自由水面處的速度勢(shì)函數(shù):Φs(x,y,t) =Φ(x,y,η,t),將其代入到自由表面邊界條件中,則自由表面條件可表達(dá)為:

(1)

(2)

(3)

式中:w(m)和Φ(m)分別表示m階的垂向速度和速度勢(shì);Φ(m)可以通過偽譜方法求解速度勢(shì)的初邊值問題得到。將式(3)代入式(1)和(2)便得到了波面升高η和表面速度勢(shì)Φs的演化方程。演化方程可以通過四階龍格庫塔法進(jìn)行求解。為了模擬造波邊界,同時(shí)引入了附加速度勢(shì)模擬造波板的運(yùn)動(dòng),造波邊界條件精確到二階非線性。同樣也加入了吸收邊界來減少波浪的反射,吸收邊界的實(shí)現(xiàn)采用對(duì)式(1)加入壓力項(xiàng)[28],使得波面升高衰減。

1.2 波浪參數(shù)

表1 波浪參數(shù)

2 結(jié)果與討論

2.1 極端波浪平均波形特征與新波模型的對(duì)比

新波模型基于線性理論,沒有譜寬的限制,可以描述極端波浪的平均波形;并且理論上平均波形也是最有可能出現(xiàn)的極端波浪波形。新波模型波面函數(shù)正比于隨機(jī)過程的自相關(guān)函數(shù),其表達(dá)式為:

(4)

式中:α為波峰高度;S(ω)為波浪場的譜密度函數(shù);σ為波浪的標(biāo)準(zhǔn)差。Lindgren推導(dǎo)出了極端波浪波峰附近波面的方差,也稱為Lindgren方差[30],其表達(dá)式為:

(5)

Walker[12]基于斯托克斯波展開式,在線性新波理論的基礎(chǔ)上加入了非線性項(xiàng),三階非線性新波模型表達(dá)式為:

(6)

式中:ηLH為希爾伯特變換后的線性新波波面;h為水深;S22和S33取值參考Walker[12]中的方法。

為了研究極端波浪的波形特征,分別從6組數(shù)值模擬工況中提取出所有波峰高度滿足ηc/Hs>1.25的波浪,并對(duì)波面進(jìn)行平均,然后與三階新波模型進(jìn)行對(duì)比。圖1為有效波高相同(Hs=0.03 m),譜寬依次減小的4組工況C01(γ=1.0),C02(γ=3.3),C03(γ=5.0)和C04(γ=7.0)模擬得到的極端波浪平均波形與三階修正新波模型的對(duì)比圖;圖中實(shí)線為極端波浪波面平均后的結(jié)果,短線為新波模型波面;灰色區(qū)間為新波模型理論平均波面95%置信區(qū)間的范圍;置信區(qū)間范圍的上下限為η±2σE;其中σE為波面均值分布的標(biāo)準(zhǔn)差,也就是均值的標(biāo)準(zhǔn)誤差;σE=(σL/N)0.5;N為極端波浪個(gè)數(shù),每個(gè)工況的極端波浪個(gè)數(shù)在200左右。圖中時(shí)間橫坐標(biāo)和空間縱坐標(biāo)分別以譜峰周期Tp和有效波高Hs做了無因次化處理,并且都將極端波浪波峰時(shí)刻作為原點(diǎn)。從對(duì)比結(jié)果可以看到新波模型與數(shù)值模擬結(jié)果吻合較好,并且平均波形都在三階新波模型95%置信區(qū)間的范圍內(nèi);但是對(duì)于譜寬最寬的工況C01,波峰兩側(cè)的波谷位置處,三階新波模型預(yù)測(cè)的波谷值較大;而對(duì)于其余工況,誤差主要出現(xiàn)在極端波浪生成之前,波浪周期相對(duì)新波模型更短,波高也更小。從圖中可以看到4組不同譜寬工況下三階新波模型均能夠較好地描述極端波浪的波形特征,但是對(duì)于窄譜,誤差相對(duì)更小。從圖1中還可以看到隨著譜寬的減小,無因次波峰高度明顯增大;4組工況平均后的波峰高度與有效波高的比值ηc/Hs分別為:1.38、1.39、1.43、1.44、并且波谷也隨著譜寬的減小而增大;表明譜寬會(huì)影響極端波浪的波高,譜寬越窄,極端波浪波高越大。除此之外,隨著譜寬的減小,極端波浪水平的不對(duì)稱性也更加明顯。極端波浪前端的波浪周期明顯小于譜峰周期,尤其是最大波前第2個(gè)波浪;而極端波后的波浪周期基本等于譜峰周期。由于極端波浪前的周期小于波浪后的周期,根據(jù)線性色散關(guān)系,極端波浪前的群速也小于極端波浪后的群速,導(dǎo)致傳遞的能量在波峰附近聚集,這說明了線性色散關(guān)系對(duì)極端波浪的生成有著重要的作用,Lighthill[31]也給出了類似的解釋,此時(shí)色散關(guān)系是極端波浪形成的重要原因。

圖1 C01、C02、C03和C04工況極端波浪平均波形與三階新波模型對(duì)比

圖2為譜寬相同,有效波高不同的4組工況模擬得到的極端波浪平均波形與三階修正新波模型的對(duì)比圖;從對(duì)比結(jié)果可以看到在波陡較大時(shí),三階修正的新波模型明顯高估最大波峰兩側(cè)的波谷深度和波峰高度,超出三階新波模型95%置信區(qū)間的范圍;并且在有效波高較大的工況C06和C07中,數(shù)值模擬得到的極端波在波峰兩側(cè)更加陡峭;產(chǎn)生誤差的原因有可能是因?yàn)楫?dāng)波陡更大時(shí),波浪的非線性更強(qiáng);另一個(gè)原因可能是因?yàn)槿A新波模型的推導(dǎo)是基于斯托克斯展開,只考慮了約束波的影響,并沒有考慮自由波的共振對(duì)波高的影響。進(jìn)一步觀察模擬得到的極端波浪波形還可以發(fā)現(xiàn):隨著有效波高的增大,平均后的波峰高度與有效波高的比值分別為:ηc/Hs=1.40,1.44,1.46,1.45,波谷深度與有效波高比值分別為ηt/Hs=0.91,0.76,0.67,0.62;波峰波谷比值為ηc/ηt=1.54,1.89,2.18,2.34;結(jié)果表明當(dāng)有效波高增大時(shí),極端波浪的垂直不對(duì)稱性顯著增大。然而4組工況的極端波浪的無量綱波高分別為H/Hs=2.31,2.20,2.13,2.06,也就是說極端波浪的無量綱波高隨著非線性的增加反而減小。在波陡較大的3組工況C05、C06和C07中基本沒有出現(xiàn)水平周期上的不對(duì)稱性,這可能是由于隨著有效波高的增大,色散關(guān)系不再為極端波浪形成的主要原因,非線性已經(jīng)成為主導(dǎo),共振導(dǎo)致局部波譜變化。

所以總體看來,修正的三階新波模型基本可以較好地描述極端波浪的平均形狀;但是對(duì)于窄譜和波陡較小的工況,誤差更小;在波陡較大時(shí),實(shí)際波浪更為陡峭,新波模型會(huì)高估最大波峰兩側(cè)的波谷和波峰高度。極端波浪的波高隨著譜寬的減小和波陡的增大而增大;而無量綱波高隨著非線性的增加反而減小。在有效波高較小時(shí),色散作用時(shí)極端波浪形成的重要原因;而當(dāng)有效波高增大時(shí),非線性占據(jù)主導(dǎo)作用,此時(shí)水平不對(duì)稱性降低,垂直不對(duì)稱性增強(qiáng)。

2.2 波高分布對(duì)比分析

Rayleigh分布是線性窄譜下的波浪理論分布,并沒有考慮非線性的影響。MER分布考慮了峰度的影響,峰度值越大,說明該數(shù)據(jù)系列中的極端值越多。高斯分布的峰度值為3.0,若實(shí)際分布的峰度值大于3.0,則表示比高斯分布存在更多的極端值。由于波浪的非線性作用,實(shí)際海浪的峰度值大于3.0。

圖3為有效波高相同,但譜寬不同的4組工況模擬得到的波高分布;從圖中可以看到在H/Hs<2.0時(shí),4種工況Rayleigh分布能夠較好地預(yù)測(cè)波高分布,基本上和模擬的波高分布重合;當(dāng)H/Hs>2.0時(shí),Rayleigh分布明顯低估了極端波浪的發(fā)生概率。2.02.5時(shí),MER分布只在工況C02時(shí)與實(shí)驗(yàn)數(shù)據(jù)吻合較好,這是由于譜寬較窄時(shí)波浪的共振非線性作用增強(qiáng)(BFI增大),極端波浪出現(xiàn)概率增加,MER分布也不能準(zhǔn)確描述極端波浪分布。

圖3 C01、C02、C03和C04工況的累積波高概率分布與Rayleigh和MER分布的對(duì)比

圖4為譜寬相同,有效波高不同的4組工況模擬得到的波高分布與理論波高分布的對(duì)比圖。同樣的Rayleigh分布低估了極端波浪發(fā)生概率,雖然MER分布比Rayleigh分布較好反映極端波,但仍會(huì)低估極端大浪的出現(xiàn)概率。

圖4 C02、C05、C06和C07工況的累積波高概率分布與Rayleigh和MER分布的對(duì)比

圖5給出了波峰分布的對(duì)比,可以看出Rayleigh分布、Forristall分布和二階Tayfun分布都低估了大波峰出現(xiàn)的概率;Forristall分布和二階Tayfun差別并不大,因?yàn)閮烧叨际强紤]了二階非線性;而三階Tayfun分布在ηc/Hs<1.25范圍內(nèi)與實(shí)驗(yàn)波峰分布吻合較好,而當(dāng)波峰更大時(shí),三階Tayfun分布也會(huì)低估極端波浪發(fā)生概率。雖然三階Tayfun也是在窄譜假定下得到的理論分布,但是從結(jié)果看來在寬譜下的結(jié)果也很好。

圖5 C01、C02、C03和C04工況的波峰概率分布與理論分布的對(duì)比

圖6也可以發(fā)現(xiàn)三階Tayfun分布在ηc/Hs<1.25時(shí)基本與實(shí)驗(yàn)分布一致,誤差主要發(fā)生在ηc/Hs>1.25時(shí);波陡越大,高階非線性作用增強(qiáng),三階Tayfun分布與實(shí)際結(jié)果相差越大;但是在有效波高最大的工況誤差反而減小了。

圖6 C02、C05、C06和C07工況波峰概率分布與理論分布的對(duì)比

3 結(jié)論

1)修正三階新波模型能夠較好地描繪極端波浪的波形;但是在波況的波陡較大時(shí),實(shí)際波浪更為陡峭,并且具有明顯的高波峰淺波谷的特點(diǎn);修正的三階新波模型會(huì)高估最大波峰兩側(cè)的波谷和波峰高度。極端波浪的波高會(huì)隨著譜寬的減小和波陡的增大而增大,但是無量綱波高反而隨著波陡的增加而減小。

2)當(dāng)H/Hs>2.0時(shí)MER分布與實(shí)際波浪分布吻合較好;當(dāng)波浪更大時(shí),對(duì)于H/Hs>2.5時(shí),MER分布只在部分工況與實(shí)驗(yàn)數(shù)據(jù)吻合較好,有效波高增大也導(dǎo)致MER誤差增大。與MER分布類似,在波峰分布中,三階Tayfun分布在ηc/Hs<1.25時(shí)基本與實(shí)驗(yàn)分布一致;譜寬越窄誤差越大,有效波高增大誤差也會(huì)先增大后回落。