張景新,王 穎,王 偉,寇 軍,王佳飛,劉 樺

(1.上海交通大學船舶海洋與建筑工程學院,上海 200240;2.上?；瘜W工業(yè)區(qū)物業(yè)管理有限公司,上海 201507)

沿海地區(qū)夏季頻繁的臺風對海岸防護工程有諸多危害,涉及越浪、護岸結(jié)構(gòu)的破壞及堤身的穩(wěn)定等。臺風引起的風浪在灘面上傳播,伴隨著波浪的破碎,引起灘面上大量的床沙起動、懸浮及輸移。灘面沖淤變化不僅影響著波浪的傳播、爬高及越浪,同時也影響著堤身結(jié)構(gòu)的穩(wěn)定,甚至引起海岸災害。

鑒于海岸防護工程的重要性,波浪作用下灘面沖淤演變的研究一直以來備受關(guān)注,國內(nèi)外近年來已有大量的研究成果[1-4],涉及試驗研究及數(shù)值模擬。試驗研究借助于大型水槽設(shè)備的開發(fā)及測量手段的不斷完善,已取得了豐碩的研究成果[4]。數(shù)值模擬手段在海灘輸沙的計算中也發(fā)揮著重要的作用。近年來,復雜流體運動模型或水沙兩相流模型在海岸帶潮流場、波浪場及泥沙運動的研究中逐漸得以發(fā)展,在水沙運動的機理性研究中發(fā)揮著重要的作用[3,5]。但對于實際海岸工程而言,大型數(shù)值模型的開發(fā)及應(yīng)用還難以達到工程實際應(yīng)用的程度。早期發(fā)展的海岸帶橫向輸沙的一維數(shù)學模型形式相對簡單,計算工作量較小,借助于實測資料及合理的經(jīng)驗系數(shù)的選取,該類模型在實際工程應(yīng)用中具有一定的適用性。一維數(shù)學模型應(yīng)用中的首要問題是確定橫向輸沙率,目前已有眾多輸沙率計算方法[6-7]。依據(jù)輸沙的動力因素不同,研究者將泥沙在空間上分為波浪上涌帶和波浪淺化帶及破碎帶分別加以考慮[7]。泥沙輸沙率計算公式繁多,如何選擇或發(fā)展適合特定水動力條件下的輸沙率公式非常重要。一維數(shù)學模型涉及的輸沙率計算公式多為經(jīng)驗型公式,但借助于適當?shù)膶崪y資料及經(jīng)驗系數(shù)的選取,往往可獲得較滿意的計算結(jié)果。

筆者建立了波浪作用下灘面演變的一維數(shù)學模型,波浪場模擬采用波浪能量守恒方程,該方程計入了波浪破碎效應(yīng)。波浪模型給出波參數(shù)的空間分布,應(yīng)用于輸沙率的計算。橫向輸沙率計算公式參考了眾多研究成果,計入了波浪上涌帶泥沙輸沙率、波浪淺化帶及破碎帶輸沙率及波生流對泥沙輸運的作用。將該模型應(yīng)用于杭州灣某海防工程堤前灘面的沖淤演變計算,模擬了某次臺風過境期間風浪所引起的灘面沖淤變化。通過數(shù)值模擬驗證了該模型的精度,分析了若干影響因素及該模型的適用性。

1 數(shù)學模型

1.1 波浪場計算

文中波浪場的計算基于規(guī)則波,通過波浪能量(或波作用量)守恒方程的求解,獲得波浪要素的空間分布。近岸水域波浪運動受地形等因素的影響,運動形式復雜,淺化、折射及繞射現(xiàn)象并存,筆者僅考察入射波浪在橫向的傳播演化及其所引起的橫向輸沙,將控制方程簡化為橫向的一維模型。近岸帶水體運動受到諸多因素的影響,如床面摩阻作用隨水深的減小逐漸顯著。而波浪破碎是該水域波浪運動的顯著特征,本文模型僅概化了破波帶內(nèi)波浪運動的破碎效應(yīng),控制方程如下:

式中:Cg(x)為波浪群速度;E(x)為單位波長的波能;Hw(x)為波高的空間分布;Dw(x)為由波浪破碎引起的波能耗散項;W為破波因子函數(shù)(W建立了空間波高與當?shù)厮畹年P(guān)系,通過設(shè)定破波指標可將波浪破碎這一物理過程模式化,大量應(yīng)用成果表明該處理方法具有相當廣泛的適用性),本文破波因子函數(shù)計算式來源于文獻[8];K為計算系數(shù);d為當?shù)厮?γ為破波指標,γ=0.73。

利用有限差分法求解上述波浪能量守恒方程,求得波浪參數(shù)的空間分布,進而可計算床面應(yīng)力等物理量,并進一步用于計算橫向輸沙率。

1.2 床面演化控制方程

床面橫向演化的一線方程給出了局部床面高程的時間變化與輸沙率橫向變化的關(guān)系,表達式如下:

式中:z為床面高程;ε為床沙孔隙率;ρs和 ρ分別為泥沙的密度和水的密度;g為重力加速度;q為輸沙率;qR為波浪上涌帶的輸沙率;qw為波浪淺化帶及破碎帶內(nèi)的輸沙率。

海岸帶橫向輸沙率的空間分布如圖1所示。圖中區(qū)域1,2,3分別代表破波帶外、破波帶和波浪上涌帶;xB,xR,xM分別為破波點的橫坐標、上涌帶的輸沙率峰值點的橫坐標和波浪最大爬升高度處的橫坐標;dB,dR分別為xB,xR處的水深。觀察圖1可知,破波帶(區(qū)域2)內(nèi)的輸沙強度最大,波浪上涌帶的輸沙強度隨著波浪的爬升逐漸減弱。利用迎風格式的有限差分法數(shù)值求解式(2),而式(2)的求解首先需要給出橫向輸沙率的計算表達式。

圖1 近岸帶橫向輸沙強度空間分布[7]

1.3 海岸帶橫向輸沙率

海岸帶橫向輸沙率按空間分布可分為2部分,即波浪上涌帶的輸沙率qR和波浪淺化帶及破碎帶內(nèi)的輸沙率qw。參考已有的研究成果[7-8],分別給出各輸沙率的計算表達式。海浪上涌帶的qR計算公式如下:

式中:βeq為床沙平衡角;β為床面坡度;c1為系數(shù),c1=10-3～10-1;R為波浪爬升高度;Hs0為入射波的有效波高;L0為深水波波長。

上涌帶外側(cè)qR值通過式(4)計算得到:

式中:c3為系數(shù),c3=0.2～0.3。

淺化帶及破波帶內(nèi)的波浪輸沙率qw可分為推移質(zhì)輸沙率qb和懸移質(zhì)輸沙率qs,參考已有的研究成果[3,7],分別給出各自的計算表達式:

式中:εb為推移質(zhì)輸沙效率 ,εb=0.1;εs為懸移質(zhì)輸沙效率,εs=0.01;Φ為床沙顆粒的內(nèi)摩擦角,tan Φ=0.6;fw為波浪摩阻系數(shù);Ud為水深平均的離岸流速;F和B分別為床面摩阻和波浪破碎引起的動能耗散率;T為波周期;Uw為波浪邊界層處水質(zhì)點運動速度幅值;z0為床面粗糙尺度;ωs為泥沙沉速;d50為床沙中值粒徑。

離岸流體運動主要分布于下層水體,該流動與向岸的波浪場質(zhì)量輸移運動和破碎波引起的上層水體的向岸運動相平衡,Ud可由式(7)計算:

式中:c為波浪相速度;E為波能。

懸移質(zhì)輸沙率計算表達式中F及B的計算參考相關(guān)研究成果[7]:

2 模型應(yīng)用及討論

將模型應(yīng)用于杭州灣某海防工程堤前灘面沖淤演變的數(shù)值模擬研究中,該海防工程位于杭州灣北岸,工程位置及地形測量斷面位置見圖2(圖中符號G代表斷面)。2005年“麥莎”臺風前后,分別開展了堤前水下地形的測量工作,2次測量前后間隔大約3個月。2次局部地形的變化主要由此次臺風過境期間的風浪所引起。首先根據(jù)現(xiàn)場實測資料確定計算水文條件,臺風風向ESE,取入射波的有效波高Hs0=2.7m,相應(yīng)周期T=6.5s,堤前平均潮位 η0=2.0m,波向由實測資料確定為正向入射。模型計算波參數(shù)取相應(yīng)的有效波參數(shù),計算歷時24h(實際臺風過境歷時約72h)。堤前灘面床沙中值粒徑d50=0.03mm,將2005年5月地形作為初始地形,采用所設(shè)定的水文參數(shù),模型給出了此次臺風過境期間風浪所引起的沖淤變化量。該海防工程堤前建有混凝土護岸的緩坡,故在輸沙率的計算中未考慮波浪上涌帶的泥沙輸沙率。

測量資料給出了沿岸4個固定斷面的地形資料,將各斷面的模擬結(jié)果與實測資料比較,如圖3所示。圖3分別給出了4個固定斷面2005年5月和8月的實測地形值和臺風后灘面地形的模擬值,x軸負向為外海方向,正向為堤身方向。計算結(jié)果顯示在堤腳的鄰近區(qū)域各斷面普遍出現(xiàn)灘面沖刷,而離岸方向普遍出現(xiàn)淤積。實測結(jié)果顯示在鄰近堤腳的區(qū)域內(nèi)灘面沖刷明顯,離岸方向沖淤幅度減弱。計算值與實測值在鄰近堤腳的區(qū)域吻合較好。分析流場動力因素,隨著水深的減小,入射波波高逐漸增大,直至破碎。破波帶內(nèi)波浪對泥沙的輸運起著非常重要的作用,而在遠離堤腳的水域,潮流作為泥沙輸運的動力因素逐漸增強。分析該水域水文實測資料,潮流運動表現(xiàn)為沿岸往復流動,與波向近似正交,其輸沙方向為沿岸。本文模型僅關(guān)注了波浪作用下的泥沙橫向輸運,未計及潮流對泥沙輸運的影響。以G1斷面為例,模擬結(jié)果顯示該斷面內(nèi)床沙的沖淤總量平衡,而實測值顯示該斷面的總沖刷量大于淤積量,意味著部分起懸泥沙可能為潮流所輸運,表現(xiàn)為沿岸輸沙。除潮流影響外,破碎波引起的沿岸流動也是重要的動力因素之一。沿岸流在淺水近岸處受地形等因素影響較弱,泥沙輸沙方向基本維持在斷面內(nèi)。遠離堤腳,隨著水深的增加,地形等影響作用減弱,沿岸流動增強,沿岸輸沙強度漸增,這一因素未計入本模型,也是影響遠離堤腳處計算偏差漸增的可能原因之一。

圖2 工程位置圖及地形測量斷面位置

圖3 灘面地形沖淤變化

3 結(jié) 語

海岸帶橫向泥沙輸運的一維數(shù)學模型以其形式簡單、計算高效等特點,具有較高的實用性。模型成功應(yīng)用的前提是橫向輸沙率的準確計算。海岸帶流動條件復雜,波浪的淺化、破碎及上涌都對泥沙運動起著重要的作用,不同輸沙形式是由不同的動力因素決定的。面對種類繁多的輸沙率公式,還有待分析總結(jié)及發(fā)展更適合的計算模型。深入研究該水域的流體運動,量化各流動物理量是改進模型的前提。此外,如文中提到的這類工程設(shè)施,不僅需要考慮風浪,還需考慮潮流作用,這方面的工作還有待深入。并且,文中模型應(yīng)用于實際海岸的床面沖淤,針對規(guī)則波,而實際風浪為隨機波列。雖然對于近岸床面的沖淤模擬結(jié)果較好,但仍需進一步研究隨機波作用下的灘面沖刷。

:

[1]RAKHA K A,DEIGAARD R,BROKER I.A phase-resolving cross shore sediment transport model forbeach profile evolution[J].Coastal Engineering,1997,31(4):231-261.

[2]SPIELMANN K,ASTRUC D,THUAL O.Analysis of some key parametrizations in a beach profile morphodynamical model[J].Coastal Engineering,2004,51(10):1021-1049.

[3]LEONT'YEV I O.Short-term shoreline changes due to crossshore structures:a one-line numerical model[J].Coastal Engineering,1997,31(1):59-75.

[4]van RIJN L C.Prediction of dune erosion due to storms[J].Coastal Engineering,2008,56(4):441-457.

[5]CHRISTENSEN E D,WALSTRA D J,EMERAT N.Vertical variation of the flow across the surf zone[J].Coastal Engineering,2002,45(3):169-198.

[6]錢寧,萬兆惠.泥沙運動力學[M].北京:科學出版社,1983:270-359.

[7]LEONT'YEV I O.Numericalmodeling of beach erosion during storm event[J].Coastal Engineering,1996,29(2):187-200.

[8]KIRBY J T,DALRYMPLE R A.Modeling waves in surf-zones and around islands[J].Journal of Waterway,Port,Coastal and Ocean Engineering,1986,112(1):78-93.