趙宇

摘 要：本文主要就海塘工程波要素計算方法進(jìn)行探討，其中設(shè)計波要素的深水波要素均采用風(fēng)推浪的計算方法，按莆田公式計算確定。淺水變形采用數(shù)模計算和經(jīng)驗公式推求，并考慮海塘前保灘順壩對波浪的消減作用；結(jié)合周邊已建工程的設(shè)計波要素，并比較數(shù)值模擬與經(jīng)驗公式的推求結(jié)果，分析本工程堤前波要素計算方法的合理性。

關(guān)鍵詞：海塘工程 波要素 數(shù)值模擬 經(jīng)驗公式

1.計算標(biāo)準(zhǔn)

某海塘達(dá)標(biāo)工程位于上海市長興島，防洪標(biāo)準(zhǔn)參照《上海市海塘規(guī)劃（2011～2020）》執(zhí)行。本工程建于1997年，原設(shè)計抗御標(biāo)準(zhǔn)為50年一遇，本次海塘達(dá)標(biāo)標(biāo)準(zhǔn)提高至200年一遇高潮位+12級風(fēng)下限（32.7m/s）。

2.計算工況

主海塘的不利風(fēng)向取N、NNE、NE、ENE、E五個方向，其中E風(fēng)向的深水波要素較其他幾個風(fēng)向的最大。

3.數(shù)值模擬計算

模型采用的是丹麥水利研究所DHI開發(fā)的MIKE21軟件中波譜模型SW進(jìn)行淺水區(qū)波浪數(shù)值計算，SW為無結(jié)構(gòu)網(wǎng)格模型，能較為高效和精確地模擬工程中岸線曲折多邊和島嶼眾多的特殊地形。

3.1計算原理

模型根據(jù)波作用守恒原理進(jìn)行波浪傳播變形計算，基本方程如下：

式中：E （σ，θ）—能量密度譜，σ—相對頻率，θ—波向角。

3.2計算范圍及網(wǎng)格

模型網(wǎng)格外海網(wǎng)格較粗，網(wǎng)格邊長約4km；在工程區(qū)對網(wǎng)格進(jìn)行加密至10m。

3.3模型參數(shù)設(shè)置

根據(jù)經(jīng)驗，長江口杭州灣地區(qū)底摩擦系數(shù)（Nikuradse Roughness）取0.001，white capping cdis和DELTAdis分別取3和0.25。

3.4模型邊界條件

E向和NE向模型邊界條件根據(jù)各計算水位的深水波要素給；N向堤前波浪受崇明島和崇明東灘阻擋作用，北邊波浪對堤前影響很小，故可不考慮波能傳入，采用側(cè)邊界。

3.5累積率計算原理

累積率計算：

在順堤起點、東側(cè)堤起點和東側(cè)堤終點各設(shè)置一個計算點經(jīng)過模型計算得到灘前的波要素，計算見表1。

4.經(jīng)驗公式計算

深水波要素計算根據(jù)《海堤工程設(shè)計規(guī)范》（GB/T51015-2014）所推薦莆田公式。由《上海風(fēng)浪實測資料統(tǒng)計分析及風(fēng)浪譜分析研究》的研究成果表明：上海水域（以大治河測站和中港測站的資料驗證）風(fēng)浪以莆田方法計算的波高、波周期與實測值吻合最好。

4.1深水波要素計算原理

根據(jù)工程水域情況，采用莆田試驗站公式進(jìn)行波浪要素計算：

（1）莆田試驗站公式：

5.結(jié)果分析與結(jié)論

由于數(shù)值模型法只能計算到灘地的波高，因此本文針對兩種方法計算的淺水波要素成果進(jìn)行對比和分析。由計算結(jié)果可知，數(shù)值模擬計算得到的淺水平均波高略小于經(jīng)驗公式計算結(jié)果，但差異不大，絕對差值在0.1m之內(nèi)，波周期差值在0.5s之內(nèi)，兩者計算結(jié)果基本一致。對比相鄰工程《長興北沿海塘達(dá)標(biāo)保灘工程初步設(shè)計》中的波要素計算成果，本工程考慮潛堤消浪后的平均波高最大值為1.37m，參照長興北沿海塘達(dá)標(biāo)保灘工程中北側(cè)堤考慮潛堤消浪后的平均波高，可以認(rèn)為本節(jié)中的波要素計算成果是合理的。

參考文獻(xiàn)：

[1]黃朝煊，袁文喜.海堤風(fēng)浪爬高計算方法探討[J].人民長江，2012，43（19）：34-37.

[2]賀海洪，陳善民，胡海忠.浙江省海塘工程波浪要素計算分析與比較[J].浙江水利科技，2005（4）：13-14.

[3]唐巨山，袁文喜，曾甄.海堤堤頂高程確定中若干問題的探討[J].水力發(fā)電，2008，34（2）： 39-42.

[4]張從聯(lián)，趙吉國，江洧.海堤設(shè)計波浪計算有關(guān)問題探討[J].水利技術(shù)監(jiān)督，2005，13（3）： 40-42.