基于條件模擬理論的波浪數(shù)值模擬

張 夢(mèng)，柳淑學(xué)，李金宣，張昊宸

(大連理工大學(xué) 海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，大連 116024)

條件模擬是地質(zhì)統(tǒng)計(jì)學(xué)中的一個(gè)概念。在地質(zhì)統(tǒng)計(jì)學(xué)中，因研究對(duì)象是一種具有結(jié)構(gòu)特征的隨機(jī)變量，對(duì)它進(jìn)行數(shù)學(xué)模擬時(shí)要求保持一定的空間相關(guān)性，即保持實(shí)驗(yàn)信息具有相同的協(xié)方差函數(shù)或變差函數(shù)，在地質(zhì)統(tǒng)計(jì)學(xué)中稱非條件模擬(unconditional simulation)；如果再增加一個(gè)更嚴(yán)格的約束條件，將模擬條件化，即令在各觀測(cè)點(diǎn)處的模擬值等于該點(diǎn)的實(shí)測(cè)值，稱其為條件模擬(conditional simulation)[1]。

近年來，由于一些海上工程規(guī)模的龐大和工程問題的復(fù)雜性，其對(duì)波浪荷載精確程度的要求也越來越高。目前在實(shí)驗(yàn)室中進(jìn)行物理模擬主要是通過模擬頻譜，不能模擬波形。很多時(shí)候研究人員常常需要研究某一特定短波列或波群對(duì)建筑物產(chǎn)生的作用，且這段短波列包含于滿足某一特定波譜的隨機(jī)波列中，這往往需要進(jìn)行逐波分析以找到滿足此條件的隨機(jī)波列，費(fèi)時(shí)費(fèi)力。因此如何有效地產(chǎn)生這樣的特殊波列對(duì)研究波浪對(duì)于復(fù)雜工程作用問題具有重要的意義。

起初，波浪模擬主要都是通過模擬頻譜的方式來實(shí)現(xiàn)的。Borgman[2]建立了線性波浪疊加法和線性過濾法，過濾法主要利用白噪聲濾波的方式對(duì)波浪進(jìn)行模擬，線性疊加法中初相位隨機(jī)分布，因此也稱為隨機(jī)相位譜法。后來，人們發(fā)現(xiàn)有時(shí)是某種特殊的波列對(duì)建筑物產(chǎn)生主要的破壞作用[3]，因此需要模擬一些特殊波列。劉思等基于改進(jìn)的波包譜[4]的經(jīng)驗(yàn)公式，建議了單向不規(guī)則波群的數(shù)值模擬方法[5]，但只是模擬了同時(shí)滿足頻譜及群性要求的波浪，并不是嚴(yán)格地模擬了不規(guī)則波列。Borgman[6]提出波浪的條件模擬概念，即在滿足波浪的統(tǒng)計(jì)特征情況下，人為加入一些約束，闡述了如何利用正態(tài)分布的條件概率理論，在時(shí)域上產(chǎn)生一段既包含特殊短波列，又屬于某一波譜的目標(biāo)波列。但當(dāng)計(jì)算波列的長(zhǎng)度增大時(shí)，此方法暴露出一個(gè)很大的缺點(diǎn)：自相關(guān)函數(shù)矩陣越來越病態(tài)導(dǎo)致計(jì)算結(jié)果出現(xiàn)嚴(yán)重偏差。Spanos[7]建立的ARMA波浪模擬模型在計(jì)算時(shí)遇到同樣的病態(tài)Toeplitz矩陣。Medina和Sanchez-Carratala[8]對(duì)數(shù)值波浪模擬中的許多計(jì)算模型進(jìn)行了比較，通過在波浪中加入白噪聲的方式很好地解決了這個(gè)問題。Gime′nez等[9]從上述方法得到啟示，并將其運(yùn)用到條件波浪模擬中，解決了算法中出現(xiàn)的病態(tài)矩陣問題。Hudspeth等[10]在二維水槽中對(duì)目標(biāo)波列進(jìn)行了重現(xiàn)，同時(shí)從能量的角度討論了嵌入波列的參數(shù)對(duì)模擬產(chǎn)生的影響。

天然海浪波列中存在的極限波浪，是造成海洋工程建筑物破壞的主要因素之一，而波浪聚焦是目前實(shí)驗(yàn)室產(chǎn)生極限波浪的主要方法，即讓不同方向、不同頻率的波浪在傳播過程中互相調(diào)制，使波能發(fā)生集中，形成極端大波。但是這樣模擬波浪的過程中畸形波的產(chǎn)生具有隨機(jī)性，因此如何實(shí)現(xiàn)在不規(guī)則波的波列中嵌入已知的畸形波過程，在特定地點(diǎn)和時(shí)間產(chǎn)生包含已知畸形波列的確定性模擬，亦是進(jìn)一步提高物理模型試驗(yàn)水平的需要。Oggiano等[11]基于數(shù)值二維水槽，針對(duì)一極限波浪過程，建立可以復(fù)演該過程的造波邊界條件，但其只是針對(duì)孤立的極限波浪過程，并非針對(duì)整個(gè)不規(guī)則波過程。Buldakov等[12]則對(duì)于給定的二維單向特定的極限波浪及其破碎型式進(jìn)行了確定性的模擬研究，但是其針對(duì)的是短序列的波群序列，而非一完整的包含畸形波的不規(guī)則波序列。Pierella等[13]對(duì)于二維極限波浪，基于流函數(shù)理論，將其嵌入一不規(guī)則波波列，實(shí)現(xiàn)其在物理水槽中的確定性的模擬。

本文基于Borgman的條件模擬波浪理論，通過加入白噪聲的方式，解決了計(jì)算中因病態(tài)矩陣導(dǎo)致的數(shù)值不穩(wěn)定問題，進(jìn)而建立了二維數(shù)值波浪模擬程序，可在滿足某一特定波譜的隨機(jī)波列中，在任意合理時(shí)間點(diǎn)，嵌入某一特定波列后，產(chǎn)生的新的隨機(jī)波列仍然滿足上述特定波譜。分別以規(guī)則波和聚焦波為嵌入波列，研究了嵌入波的長(zhǎng)度、周期、波高和嵌入點(diǎn)等對(duì)模擬波浪產(chǎn)生的影響。

1 基于條件模擬的波浪模擬方法

所謂條件波浪模擬，首先通過基于給定波浪參數(shù)和波浪譜，模擬生成一初始的隨機(jī)波列ηu(t)，可稱之為非條件模擬；之后在指定的時(shí)間點(diǎn)嵌入一段特殊的波列ηe(t)，從而形成新的同時(shí)仍滿足給定的波浪譜和統(tǒng)計(jì)特征參數(shù)的隨機(jī)波列ηc(t)。

對(duì)于初始的隨機(jī)波列ηu(t)，可采用非條件波浪模擬得到。即假定海浪可看作一平穩(wěn)隨機(jī)過程，由無數(shù)多個(gè)不同的余弦波隨機(jī)疊加而成[3]

(1)

(2)

把代表M個(gè)區(qū)間內(nèi)波能的各余弦波動(dòng)疊加起來，即得不規(guī)則波的波面

(3)

為了在初始波列中某一時(shí)間點(diǎn)嵌入某一特定的波列，可將初始波列的波面ηu(t)前后分為3段，表示為{ηu-,vu,ηu+}T，其中vu代表將要被嵌入的短波列替代的部分，而目標(biāo)波列可表示為{ηc-,vu,ηc+}T；若ve表示將要嵌入的波列，即要求

vc=ve

(4)

接下來就是如何求解ηc-、ηc+，使得所形成的新的波列既滿足所要求的波浪譜，同時(shí)嵌入已知的某特定波列。這里假設(shè)波浪的波面滿足正態(tài)分布，則由多向正態(tài)分布的條件概率公式可得

(5)

式中：C12T、C11-1為初始波列自相關(guān)函數(shù)矩陣的分塊矩陣，其表達(dá)式可參見Hudspeth等[10]，詳細(xì)證明可見Borgman[6]。

2 數(shù)值模擬及結(jié)果分析

2.1 嵌入規(guī)則波數(shù)值模擬結(jié)果分析

表1 初始波列參數(shù)Tab.1 Unconditional simulation wave parameters

表2 嵌入規(guī)則波參數(shù)Tab.2 Embedded regular wave parameters

(6)

式中：ηu為初始波列的過程線；ηc為目標(biāo)波列的過程線；N為序列長(zhǎng)度；Hmax為初始波列的最大波高。需要說明的是，式(6)中誤差累加計(jì)算不包括嵌入波列的部分。

1-a 組次R1 1-b 組次R2

圖2 初始波列的自相關(guān)函數(shù)曲線圖 圖3 波面誤差隨嵌入點(diǎn)的變化(組次R2)Fig.2 Auto-correlation of the unconditional simulation wave Fig.3 Variation of RMS error coefficient with the embedded time point

圖4 目標(biāo)波列與初始波列特征波高的比值隨嵌入點(diǎn)的變化(組次R2) 圖5 目標(biāo)波列與初始波列特征周期的比值隨嵌入點(diǎn)的變化(組次R2)Fig.4 Comparison of the variation of the ratios between the unconditional simulation and conditional simulation characteristic wave height with the embedded time point Fig.5 Comparison of the ratios between the variation of the unconditional simulation and conditional simulation characteristic wave period with the embedded time point

圖6為各組次目標(biāo)波列與初始波列的波譜比較圖，表3為嵌入規(guī)則波后波列各統(tǒng)計(jì)參數(shù)計(jì)算結(jié)果。具體來說，從組次R1、R2的統(tǒng)計(jì)參數(shù)和譜對(duì)比來看，進(jìn)一步說明嵌入點(diǎn)的改變對(duì)模擬產(chǎn)生的波浪分析結(jié)果影響不大；同樣，從組次R1、R3的結(jié)果對(duì)比來看，當(dāng)嵌入規(guī)則波列的長(zhǎng)度增大了9倍，但其統(tǒng)計(jì)參數(shù)和譜的變化亦較小，因此可以認(rèn)為當(dāng)嵌入規(guī)則波列的周期和波高與初始波列的有效周期和波高相差不大時(shí)，嵌入規(guī)則波列的長(zhǎng)度變化不會(huì)對(duì)模擬波浪產(chǎn)生明顯的影響；從組次R1、R4的結(jié)果對(duì)比來看，嵌入規(guī)則波列的波高增大了1倍，而其統(tǒng)計(jì)有效波高參數(shù)和譜亦沒有產(chǎn)生很大的變化，但是如果嵌入的規(guī)則波波高較大或大于初始波列的最大波高時(shí)，新的波列的最大波高會(huì)大于初始波列的最大波高。

而從組次R1、R5的結(jié)果對(duì)比來看，嵌入波列的周期會(huì)對(duì)模擬得到的新的波列產(chǎn)生較大的影響。從圖1組次R5可以明顯看出，目標(biāo)波列與初始波列相比，其在嵌入波列周圍產(chǎn)生了較大的特異波浪，從表3中也可以看到各統(tǒng)計(jì)波高明顯增大，從圖6組次R1、R5的對(duì)比可以看到，嵌入規(guī)則波后的波浪譜整體明顯變大。也就是說，嵌入規(guī)則波的周期與初始波列譜峰周期偏離越大，嵌入特定規(guī)則波后波浪的特征與初始波列特征相差就越大。

6-a 組次R1、R2 6-b 組次R1、R3

表3 嵌入規(guī)則波時(shí)初始波列和目標(biāo)波列的統(tǒng)計(jì)參數(shù)表Tab.3 Statistical parameters of unconditional simulation and conditional simulation after embedding regular waves

綜合上述討論可以看出，當(dāng)嵌入規(guī)則波列時(shí)，嵌入點(diǎn)、嵌入波列的長(zhǎng)度和波高變化對(duì)初始波列的統(tǒng)計(jì)特征參數(shù)和譜的影響不大，可以認(rèn)為條件模擬是有效的，但是當(dāng)嵌入規(guī)則波的波高較大時(shí)，嵌入規(guī)則波后的波列的最大波高會(huì)增大。事實(shí)上，從某種程度上看，嵌入波列的周期變化也沒有對(duì)統(tǒng)計(jì)參數(shù)和譜產(chǎn)生較大的影響，但是由于嵌入波列的周期偏離了初始波列的有效周期，統(tǒng)計(jì)特性的不同使得初始波列產(chǎn)生類似“排斥”效應(yīng)，從而在嵌入波浪附近產(chǎn)生了特異波浪。這種含有遠(yuǎn)遠(yuǎn)大于隨機(jī)波列有效波高的特異波的長(zhǎng)波列不滿足常規(guī)波浪特征，因此嵌入波列的周期應(yīng)與初始波列的周期基本一致。

2.2 病態(tài)矩陣問題

在上述計(jì)算過程中，有時(shí)會(huì)出現(xiàn)如圖7所示情況。圖7中7-a、7-b、7-c分別代表嵌入一個(gè)波高為0.1 m，周期分別為2.16 s、2.18 s、2.2 s，長(zhǎng)度為一個(gè)周期的規(guī)則波。在周期改變很小的情況下，圖7-b的模擬結(jié)果卻出現(xiàn)了很大的波動(dòng)。其原因是因?yàn)橛?jì)算公式(5)中的矩陣C11有時(shí)為病態(tài)的，使得計(jì)算結(jié)果不穩(wěn)定。Spanos[7]通過ARMA算法進(jìn)行波浪模擬時(shí)遇到了相似的問題。根據(jù)Spanos和Mignolet[14]的研究，發(fā)現(xiàn)其原因是譜在低頻率時(shí)波浪譜密度值為0引起的。根據(jù)此特點(diǎn)，Medina和Sanchez-Carratala[8]通過在波浪譜中人為地加入白噪聲，從而很好地解決了這個(gè)問題。假定加入的白噪聲占總能量的ε，則新的譜為

7-a Te=2.16 s 7-b Te=2.18 s

(7)

式中：fL、fH分別表示頻率的上下界，其取決于模擬波浪時(shí)在低頻和高頻處略去的波浪的能量占波浪總能量的比例[3]，本文取1%。雖然此種方法改變了原有譜形，但是Gime′nez等[9]認(rèn)為無論是實(shí)驗(yàn)室產(chǎn)生的波浪還是海面上記錄的波浪都存在著噪聲，因此當(dāng)ε取得足夠小時(shí)，S*可近似等于S。圖7-d為取ε=0.000 5所計(jì)算得到的新結(jié)果，與圖7-b相比有了很大改善。本文所涉及到的譜的計(jì)算全部采用S*。

2.3 嵌入聚焦波數(shù)值模擬結(jié)果分析

2.3.1 二維聚焦波的產(chǎn)生方法(相速度法)[15]

聚焦波可通過模擬頻譜的方法產(chǎn)生。對(duì)于二維波浪，任一點(diǎn)處波浪自由表面可以表示為不同頻率和不同幅值的余弦波線性疊加后的結(jié)果，即

(8)

式中：aj為第j個(gè)組成波的幅值；kj為第j個(gè)組成波的波數(shù)；ωj為波角頻率；φj為j個(gè)組成波的初相位；N為組成波的總個(gè)數(shù)。kj和ωj滿足波浪的色散關(guān)系

(9)

式中：g和d分別為重力加速度和水深。若假定波浪在指定時(shí)刻t=tb聚焦于x=xb位置處，即各組成波在x=xb處波峰疊加，那么需滿足

cos(kjx-ωjt-φj)=1

(10)

則各組成波的初相位應(yīng)滿足下式

φj=kjxb-ωjtb

(11)

此時(shí)，波浪的波面表達(dá)式為

(12)

由式(12)可以看出，對(duì)于給定的聚焦時(shí)間tb，聚焦點(diǎn)的坐標(biāo)xb，頻率數(shù)N，影響波浪聚焦面參數(shù)主要為各組成波的幅值aj。通常把A定義為聚焦點(diǎn)處的波浪幅值，則有

(13)

則生成聚焦波浪組成波的幅值可以表示為

(14)

式中：頻譜S(f)為本文采用的JONSWAP譜。

2.3.2 數(shù)值模擬及結(jié)果分析

以聚焦點(diǎn)為中點(diǎn)向兩邊對(duì)稱截取足夠長(zhǎng)度得到欲嵌入的聚焦波列，表4為嵌入的聚焦波參數(shù)。同理，分別研究嵌入點(diǎn)、嵌入聚焦波的歷時(shí)長(zhǎng)度、幅值和譜峰周期對(duì)模擬波浪產(chǎn)生的影響。

8-a 組次fc1 8-b 組次fc2

9-a 組次fc1、fc2 9-b 組次fc1、fc3

表4 嵌入聚焦波參數(shù)Tab.4 Embedded focusing wave parameters

圖8為嵌入聚焦波后各組次目標(biāo)波列與初始波列的波面比較圖，從結(jié)果來看，其主要特點(diǎn)與嵌入規(guī)則波時(shí)類似，即嵌入聚焦波列后目標(biāo)波列波形只在嵌入波列周圍發(fā)生變化。圖9為嵌入聚焦波后各組次目標(biāo)波列與初始波列的波譜比較圖，表5為嵌入聚焦波后波列各統(tǒng)計(jì)參數(shù)變化表。從各組次對(duì)比結(jié)果來看，嵌入聚焦波時(shí)，其影響因素也與嵌入規(guī)則波時(shí)相似，即嵌入點(diǎn)、嵌入聚焦波的長(zhǎng)度、幅值不會(huì)對(duì)模擬產(chǎn)生明顯的影響，而嵌入聚焦波的譜峰周期會(huì)對(duì)模擬波浪產(chǎn)生較大的影響。

表5 嵌入聚焦波時(shí)初始波列和目標(biāo)波列的統(tǒng)計(jì)參數(shù)表Tab.5 Statistical parameters of unconditional simulation and conditional simulation after embedding focusing waves

這里也可以看到與嵌入規(guī)則波不同的地方。從fc1和fc3組成結(jié)果的對(duì)比可以看到，聚焦波列的長(zhǎng)度變大反而使得其統(tǒng)計(jì)參數(shù)和譜變化更小了，這是因?yàn)榫劢共械拈L(zhǎng)度越長(zhǎng)，統(tǒng)計(jì)特性更加與初始波列趨于一致，初始波列產(chǎn)生的“排斥”因此就越小。從周期的影響上來看，聚焦波譜峰周期的變化比規(guī)則波周期的變化產(chǎn)生的影響更大，同樣的嵌入聚焦波的周期為1.5倍的初始波列有效周期，嵌入聚焦波時(shí)產(chǎn)生的特異波是嵌入規(guī)則波時(shí)的1.7倍。

3 結(jié)論

本文基于Borgman的條件波浪模擬理論，建立了相應(yīng)的數(shù)值波浪模擬程序，可實(shí)現(xiàn)模擬生成既滿足某一特定波浪譜，又包含了某一特殊短波列的隨機(jī)波浪，并通過加入白噪聲的方式解決了計(jì)算中產(chǎn)生的數(shù)值不穩(wěn)定問題。

嵌入波浪參數(shù)對(duì)于嵌入波浪后模擬波浪的影響研究結(jié)果表明，在進(jìn)行條件波浪模擬時(shí)，嵌入波列的周期是主要影響因素，其與初始波列的有效周期越接近，則模擬效果越好，因此需要限定嵌入波列周期在一定區(qū)間范圍內(nèi)。而在限定周期范圍內(nèi)，嵌入點(diǎn)、嵌入波列的波高和長(zhǎng)度的變化不會(huì)對(duì)模擬產(chǎn)生明顯的影響，但是如果嵌入波列的波高較大，嵌入波列后波浪的最大波高會(huì)增大。