安蒙華，蔣勤

（1.中交第三航務(wù)工程勘察設(shè)計(jì)院有限公司，上海200032；2.河海大學(xué)，江蘇南京210029）

超標(biāo)準(zhǔn)潮波作用下越浪量及越浪流特性研究

安蒙華1，蔣勤2

（1.中交第三航務(wù)工程勘察設(shè)計(jì)院有限公司，上海200032；2.河海大學(xué)，江蘇南京210029）

為探求臺(tái)風(fēng)等引起的超標(biāo)準(zhǔn)風(fēng)暴潮過境時(shí)海堤的破壞機(jī)理，通過物理模型試驗(yàn)對(duì)超標(biāo)準(zhǔn)水動(dòng)力條件下帶胸墻斜坡堤的越浪量及越浪流特性進(jìn)行了研究?；谝?guī)則波的試驗(yàn)數(shù)據(jù)，建立了相對(duì)堤頂超高、相對(duì)墻頂超高以及相對(duì)胸墻高度與越浪量之間的經(jīng)驗(yàn)關(guān)系公式。同時(shí)提出了在研究堤頂前端水舌厚度與越浪量關(guān)系時(shí)應(yīng)采用的無因次參數(shù)形式以及在研究越浪流堤頂沿程變化時(shí)應(yīng)采用的相對(duì)位置參數(shù)。

超標(biāo)準(zhǔn)潮位；超標(biāo)準(zhǔn)波高；斜坡式海堤；胸墻；越浪量；越浪流

0 引言

我國(guó)地處太平洋西岸，由于受到形成于西北太平洋的熱帶風(fēng)暴的影響，夏季東南沿海經(jīng)常發(fā)生風(fēng)暴潮災(zāi)害。海堤作為海岸防護(hù)的工程設(shè)施，是抵御風(fēng)暴災(zāi)害的重要屏障。近年來，東南沿海地區(qū)發(fā)生了多次超過海堤原有設(shè)計(jì)標(biāo)準(zhǔn)的風(fēng)暴潮災(zāi)害，使得海堤因波浪越頂而破壞，甚至潰決。但要保證海堤不發(fā)生越浪是不現(xiàn)實(shí)的，在綜合考慮安全性及經(jīng)濟(jì)性的基礎(chǔ)上，很多國(guó)家包括我國(guó)東南沿海的部分省份都采用了允許部分越浪的海堤設(shè)計(jì)標(biāo)準(zhǔn)。對(duì)于這一類海堤，平均越浪量就成為衡量其安全性的一個(gè)重要設(shè)計(jì)指標(biāo)。

自20世紀(jì)初，國(guó)內(nèi)外學(xué)者對(duì)海堤越浪量課題開展了深入的研究。譬如，T.Saville[1]、合田良實(shí)[2]、Owen[3]、Van Der Meer[4-5]等先后對(duì)海堤的越浪量問題進(jìn)行了大量的試驗(yàn)研究和分析工作。特別是Van Der Meer提出的越浪量計(jì)算公式，被列入許多歐洲國(guó)家的海堤設(shè)計(jì)規(guī)范。在我國(guó)也有大量針對(duì)斜坡堤越浪的研究報(bào)告及相關(guān)的經(jīng)驗(yàn)公式。虞克等[6]（1992）、賀朝敖等[7]（1995）、周益人[8]（2007）、陳國(guó)平[9]（2010）都各自提出了斜坡堤平均越浪量計(jì)算公式。隨著對(duì)海堤越浪問題研究的不斷深入，學(xué)者們開始關(guān)注越浪流水動(dòng)力特性及其對(duì)海堤安全穩(wěn)定性的影響。Holger Schüttrumpf[10]（2001）及Van-Gent[11]（2002）、代英男[12]（2011）、朱偉娜[13]（2012）相繼對(duì)越浪流進(jìn)行了研究，探討了越浪流及相關(guān)參數(shù)的計(jì)算方法。

盡管目前已經(jīng)有許多關(guān)于越浪量及越浪流的試驗(yàn)成果及相應(yīng)的經(jīng)驗(yàn)計(jì)算公式，但這些計(jì)算公式適用的海況及海堤的結(jié)構(gòu)形式不同，幾乎每個(gè)經(jīng)驗(yàn)公式都存在適用性的限制。目前針對(duì)風(fēng)暴過境時(shí)產(chǎn)生的超過設(shè)計(jì)標(biāo)準(zhǔn)的水動(dòng)力條件對(duì)斜坡式海堤影響的研究尚不成熟。因此本文通過物理模型試驗(yàn)方法，針對(duì)超標(biāo)準(zhǔn)潮位及波浪條件，對(duì)帶胸墻斜坡堤的越浪量及越浪流的變化特性進(jìn)行了研究，基于規(guī)則波的試驗(yàn)結(jié)果，提出了越浪量值與胸墻相關(guān)參數(shù)的經(jīng)驗(yàn)公式，并對(duì)堤頂前端的水舌厚度及其在堤頂?shù)难爻套兓匦赃M(jìn)行了分析。本研究提出的分析方法對(duì)把握臺(tái)風(fēng)對(duì)海堤的破壞機(jī)理，評(píng)估超標(biāo)準(zhǔn)潮位和波高條件下斜坡式海堤越浪量及越浪流特性都具有一定的參考意義和實(shí)用價(jià)值。

1 模型試驗(yàn)

本次試驗(yàn)選取了設(shè)計(jì)標(biāo)準(zhǔn)為30 a一遇的海堤斷面。試驗(yàn)采用規(guī)則波。選取的潮位和波浪條件分別為30 a，50 a，100 a及200 a一遇，并進(jìn)行潮位和波浪的不同組合，試驗(yàn)工況均屬于超設(shè)計(jì)標(biāo)準(zhǔn)的水動(dòng)力條件。

1.1 儀器設(shè)備

物理模型試驗(yàn)水槽長(zhǎng)80 m、寬1.0 m、高1.5 m。其中一端安裝由計(jì)算機(jī)自動(dòng)控制的不規(guī)則造波機(jī)，生成所需要模擬的波浪要素，另一端設(shè)置消波系統(tǒng)。試驗(yàn)時(shí)將水槽分割為寬0.5 m的兩部分，一部分鋪設(shè)試驗(yàn)斷面，另一部分用以消除波浪的二次反射。

平均越浪量的測(cè)量采用圖1所示的接水裝置，其中導(dǎo)水槽寬10 cm。

圖1 越浪量測(cè)量方法示意圖Fig.1Metering method of overtopping rate

越浪流的測(cè)定采用電容式浪高儀及DJ800波浪采集系統(tǒng)，采樣頻率取100 Hz，采樣時(shí)長(zhǎng)取60 s。

在設(shè)置波高儀時(shí)，首先在光滑平板上挖出與波高儀底部絕緣體相同尺寸的空隙，再將波高儀置入其中，以消除波高儀震動(dòng)對(duì)波高數(shù)據(jù)測(cè)量精度的影響。

1.2 試驗(yàn)方法

本試驗(yàn)采用正態(tài)模型，模型幾何比尺取1∶15，根據(jù)Froude數(shù)相似率，各物理量模型比尺如下：

此外，采用規(guī)則波造波，每組試驗(yàn)的持續(xù)造波時(shí)間為120 s，并重復(fù)3次。越浪量取3次量測(cè)結(jié)果的平均值。

1.3 試驗(yàn)斷面

試驗(yàn)斷面如圖2所示，試驗(yàn)中取斷面前坡坡度為1∶2.5，后坡坡度為1∶2，堤頂高度為46 cm（原型6.9 m），堤頂寬度為53.3 cm（原型8 m）。采用4種不同的堤前水深，在34.9～41.7 cm之間（原型5.24～6.25 m），直立式胸墻選取3種不同的高度，在3.3～10 cm之間（原型0.5～1.5 m）。

圖2 試驗(yàn)斷面圖Fig.2Dike profile of experiment

模型制作過程中要保證其重量和幾何相似，重量誤差控制在3%以內(nèi)，幾何誤差控制在1%以內(nèi)。

1.4 試驗(yàn)組次

不同的潮位和規(guī)則波組合如表1所示，分別對(duì)胸墻高度P為0 cm、3.3 cm、6.7 cm、10 cm（原型0 m、0.5 m、1.0 m、1.5 m）的4種海堤斷面進(jìn)行了64組工況的越浪量試驗(yàn)。

2 試驗(yàn)結(jié)果分析

2.1 越浪量

影響斜坡堤越浪量的因素包括水動(dòng)力要素和海堤結(jié)構(gòu)形式，即：

式中：H為堤前波高；T為波周期；Hc為胸墻頂?shù)届o水面的距離；Rc為堤頂?shù)届o水面的距離；P為胸墻高度；B為堤頂寬度；d為堤前水深；m為前坡坡度；mb為堤前底坡；k△為糙滲系數(shù)。

表1 海堤越浪量試驗(yàn)各工況中潮位波浪條件Table 1Tide and wave conditions for overtopping rate experiment in different conditions

為考察胸墻高度對(duì)越浪量的影響，本次試驗(yàn)中取B=8 m，m=2.5，mb=+∞，k△=1，并取H/L=0.49為固定值，保持各試驗(yàn)工況下的波陡基本相同。鑒于Hc=Rc+P，將堤前波高H作為基本變量，可建立如式（2）所示的海堤越浪量Q與相對(duì)墻頂超高Hc/H，相對(duì)堤頂超高Rc/H及相對(duì)胸墻高度P/H之間的定量關(guān)系，進(jìn)而考察這3個(gè)無因次變量對(duì)海堤越浪量的影響。

式中：Rc/H為堤頂相對(duì)高度；Hc/H為墻頂相對(duì)高度；P/H為胸墻相對(duì)高度。

2.1.1 越浪量與墻頂超高的關(guān)系

賀朝敖[7]（1995）基于水力學(xué)中流量系數(shù)的概念，將相對(duì)墻頂超高Hc/H作為基本變量，利用最小二乘法原理對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行線性回歸，建立了海堤越浪量與各影響參數(shù)的經(jīng)驗(yàn)關(guān)系式，指出流量系數(shù)2πQ/HL與Hc/H具有指數(shù)關(guān)系。

參照賀朝敖[7]的方法，本研究以相對(duì)墻頂超高Hc/H和相對(duì)越浪量Q/（gH3）0.5作為基本變量，對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行整理，得到圖3所示的關(guān)系曲線。

圖3 相對(duì)越浪量Q/(gH3)0.5與Hc/H之間的關(guān)系Fig.3Relation between dimensionless overtopping rate Q/(gH3)0.5and Hc/H

相對(duì)越浪量與相對(duì)堤頂超高間具有明顯的指數(shù)關(guān)系。在Hc/H=0.3～2.0之間，越浪量隨堤頂超高以指數(shù)關(guān)系減少，這與賀朝敖[7]（1995）對(duì)一般外海潮波條件下越浪量與堤頂超高間的研究結(jié)論相吻合。

2.1.2 越浪量與堤頂超高的關(guān)系

對(duì)于相對(duì)堤頂超高Rc/H和相對(duì)越浪量Q/（gH3）0.5，試驗(yàn)數(shù)據(jù)顯示兩者并沒有特定的相關(guān)關(guān)系，數(shù)據(jù)較散亂。如圖4所示。

通過擬合分析，得到置信度為95%的超標(biāo)準(zhǔn)風(fēng)暴潮條件下相對(duì)越浪量Q/（gH3）0.5與相對(duì)墻頂超高Hc/H的經(jīng)驗(yàn)關(guān)系公式如下：

圖4 相對(duì)越浪量Q/(gH3)0.5與Rc/H之間的關(guān)系Fig.4Relation between dimensionless overtopping rate Q/(gH3)0.5與Rc/H

2.1.3 越浪量與胸墻高度的關(guān)系

周益人[9]（2008）基于物理模型試驗(yàn)指出胸墻高度P對(duì)堤頂越浪量具有重要影響，并提出了考慮胸墻高度的斜坡式海堤越浪量計(jì)算公式。參照其方法，本節(jié)以相對(duì)胸墻高度P/H和相對(duì)越浪量Q/（gH3）0.5為基本參數(shù)，對(duì)胸墻高度P>0工況的試驗(yàn)數(shù)據(jù)進(jìn)行整理。通過擬合分析，發(fā)現(xiàn)二者間具有指數(shù)關(guān)系，其經(jīng)驗(yàn)關(guān)系式為：

從圖5可見，在P/H=0.2～0.95之間，盡管試驗(yàn)值的分布相對(duì)分散，但是相對(duì)越浪量值隨相對(duì)胸墻高度的變化的總體趨勢(shì)與擬合得到的指數(shù)曲線一致，表明當(dāng)相對(duì)胸墻高度增加時(shí)，相對(duì)越浪量值隨之成指數(shù)關(guān)系減少。

圖5 相對(duì)越浪量Q/(gH3)0.5與P/H之間的關(guān)系Fig.5Relation between dimensionless overtopping rate Q/(gH3)0.5and P/H

2.2 堤頂越浪流

當(dāng)發(fā)生越浪時(shí)，無胸墻結(jié)構(gòu)的斜坡式海堤，會(huì)形成連續(xù)的堤頂越浪流。這一水體會(huì)引起堤體的滲透及侵蝕，是造成海堤破壞的主要原因之一。在實(shí)際海堤設(shè)計(jì)中，不僅要評(píng)估海堤越浪量大小，還要對(duì)堤頂越浪流等因素進(jìn)行定性及定量分析。本節(jié)針對(duì)超標(biāo)準(zhǔn)風(fēng)暴潮作用下，堤頂越浪流水舌厚度與越浪量之間的關(guān)系及其沿程變化特性進(jìn)行研究。

2.2.1 堤頂前端水舌厚度與越浪量的關(guān)系

Schüttrumpf及Van Gent以堤頂前端水舌厚度hc（xc=0）作為主要參數(shù)，分別給出了堤頂越浪流的流速及水舌厚度的計(jì)算公式。再利用堤頂末端流速與水舌厚度的乘積來計(jì)算堤頂越浪量。參照Schüttrumpf及Van Gent的方法，本節(jié)對(duì)堤頂前端水舌厚度與越浪量之間的關(guān)系進(jìn)行了研究。

試驗(yàn)結(jié)果表明，相對(duì)越浪量Q/（gH3）0.5與以波高為基本參數(shù)描述的相對(duì)堤頂水舌厚度hc/H之間雖具有一定的線性關(guān)系，但是試驗(yàn)數(shù)據(jù)的分布相對(duì)散亂。而相對(duì)越浪量與以堤頂超高為基本參數(shù)描述的相對(duì)堤頂水舌厚度hc/Rc具有良好的指數(shù)關(guān)系。本次擬合區(qū)間為hc/Rc=0.03～0.09，結(jié)果如圖6所示。

圖6 相對(duì)水舌厚度hc(xc=0)/Rc與Q/(gH3)0.5的關(guān)系圖Fig.6Relation between dimensionless overtopping flowthichness hc(xc=0)/Rcand Q/(gH3)0.5

以堤前波高H和堤頂相對(duì)超高Rc為基本參數(shù)描述的相對(duì)水舌厚度，后者較前者與相對(duì)越浪量之間的關(guān)系更加明顯。建議在對(duì)堤頂水舌厚度與越浪量關(guān)系的研究中采用hc（xc=0）/Rc的形式描述相對(duì)水舌厚度。

2.2.2 堤頂水舌厚度沿堤頂位置xc/B的變化趨勢(shì)

Schüttrumpf指出堤頂相對(duì)水舌厚度hc（xc）/ hc（xc=0）只與計(jì)算點(diǎn)在堤頂?shù)南鄬?duì)位置xc/B有關(guān)。并以超出靜水面的波浪虛擬爬高值作為計(jì)算堤頂越浪流水舌厚度的邊界條件，計(jì)算越浪水體在堤頂?shù)乃嗪穸取?/p>

依據(jù)Schüttrumpf的研究思路，本節(jié)對(duì)hc（xc）/ hc（xc=0）隨xc/B的變化規(guī)律進(jìn)行探討。圖7為在30 a、50 a、100 a、200 a一遇的波浪條件下堤頂相對(duì)水舌厚度隨位置的變化規(guī)律。其中不同類型的曲線代表不同重現(xiàn)期的水位條件。圖中的虛線為由Schüttrumpf的理論公式得出的計(jì)算結(jié)果。

選用hc/Rc作為基本變量來描述相對(duì)水舌厚度與相對(duì)越浪量之間的經(jīng)驗(yàn)關(guān)系公式為：

圖7 不同波浪條件下hc/hc(xc=0)隨xc/B的變化規(guī)律Fig.7Relation between hc/hc(xc=0)and xc/B withdifferent wave conditions

與Schüttrumpf的公式計(jì)算值相比，在30 a及50 a一遇的波浪條件下，實(shí)測(cè)水舌厚度在2號(hào)與3號(hào)測(cè)點(diǎn)存在一定的波動(dòng)。從4號(hào)測(cè)點(diǎn)開始接近Schüttrumpf的公式計(jì)算值。在100 a一遇的波浪條件下，試驗(yàn)數(shù)據(jù)大于理論值，但差距較小。從5號(hào)測(cè)點(diǎn)處開始，貼近理論值。在200 a一遇波浪條件下，試驗(yàn)值大于理論值，在6號(hào)測(cè)點(diǎn)處才開始貼近理論值。

2.2.3 堤頂水舌厚度隨xc/L0的變化趨勢(shì)

Bosman[14]（2007）在對(duì)于堤頂越浪流的研究中引入無因次參數(shù)xc/L0（L0為深水波長(zhǎng)）。參照其研究方法，以不同重現(xiàn)期水位分類，得到圖8所示的堤頂相對(duì)水舌厚度與xc/L0之間的關(guān)系。其中，不同類型的曲線代表不同重現(xiàn)期的波浪條件。

圖8 不同波浪條件下hc/hc(xc=0)隨xc/L0的變化規(guī)律Fig.8Relation between hc/hc(xc=0)and xc/L0with different wave conditions

由圖8可見，相同水位條件下堤頂中前部水舌厚度的規(guī)律性較之相同波浪條件下的情況要差。但水體在堤頂中后部達(dá)到穩(wěn)定后所表現(xiàn)出來的規(guī)律性要明顯優(yōu)于前者。特別是在30 a及50 a一遇的工況下，堤頂中后部的數(shù)據(jù)基本在一條直線上。

綜上所述，研究堤頂前端水體運(yùn)動(dòng)對(duì)于水舌厚度的影響，宜以xc/B作為相對(duì)位置參數(shù)進(jìn)行考察。而對(duì)于研究堤頂中后部水舌厚度的變化趨勢(shì)，宜選取xc/L0作為相對(duì)位置參數(shù)。

3 結(jié)語

本文采用物理模型試驗(yàn)方法，針對(duì)超標(biāo)準(zhǔn)潮位及波浪條件下帶胸墻斜坡堤的越浪過程進(jìn)行了研究。通過試驗(yàn)數(shù)據(jù)分析了斜坡式海堤相對(duì)越浪量與相對(duì)堤頂超高、相對(duì)墻頂超高以及相對(duì)胸墻高度之間的相關(guān)關(guān)系。同時(shí)，基于無胸墻工況下的越浪流試驗(yàn)結(jié)果，提出了在研究堤頂前端水舌厚度與越浪量關(guān)系時(shí)應(yīng)采用的無因次參數(shù)形式以及在研究越浪流堤頂沿程變化時(shí)應(yīng)采用的相對(duì)位置參數(shù)。

[1]SAVILLE T J R.Laboratory date on wave run-up and overtopping on shore structures[R].Washington D C:Beach Erosion Board,US. Army Corps of Engineers,1955.

[2]合田良實(shí).港口建筑物的防浪設(shè)計(jì)[M].劉大中，孫巨才，譯.北京：海洋出版社，1984.

GODA.Wave protection design of harbor structures[M].LIU Da-zhong,SUN Ju-cai,translate.Beijing:China Ocean Press, 1984.

[3]OWEN M W.Design of seawalls allowing for overtopping[R]. Wallingford:Hydraulics Research Station,1980.

[4]DE WAAL J P,VAN DEER MEER J W.Wave run-up and overtopping on coastal structures[C]//Proceedings of the twenty-third international conference on coastal engineer.Venice,Italy,1992.

[5]VAN DEER MEER J W.Wave run-up and wave overtopping at dikes[R].Delft,Netherland:Technical Advisory Committee on Flood Defence,2002.

[6]虞克，余廣明.斜坡堤越浪實(shí)驗(yàn)研究[J].水利水運(yùn)科學(xué)研究. 1992（3）：211-219.

YU Ke,YU Guang-ming.Experimental study of overtopping volume on sloping breakwater[J].Hydro-Science and Engineering, 1992(3):211-219.

[7]賀朝敖，任佐皋.帶胸墻斜坡堤越浪量的試驗(yàn)研究[J].海洋工程，1995，5（2）：62-70.

HE Zhao-ao,REN Zuo-gao.Experimental study on wave overtopping of rubble mound breakwater with vertical wall[J].The Ocean Engineering,1995,5(2):62-70.

[8]周益人.波浪作用下堤壩防護(hù)問題試驗(yàn)研究[D].南京：河海大學(xué)，2008.

ZHOU Yi-ren.Research on dike protection issues under wave effect[D].Nanjing:Hohai University,2008.

[9]陳國(guó)平，周益人，嚴(yán)士常.不規(guī)則波作用下海堤越浪量試驗(yàn)研究[J].水運(yùn)工程，2010（3）：1-6.

CHEN Guo-ping,ZHOU Yi-ren,YAN Shi-chang.Test study on wave overtopping under irregular wave action[J].Port&Waterway Engineering,2010(3):1-6.

[10]SCHüTTRUMPF H.Wave overtopping flow on sea dikes-experimental and theoretical investigation[D].Braunschweig:Technical University of Braunschweig,2001.

[11]VAN GENT M R A.Low-exceedance wave overtopping events: measurements of velocities and the thickness of water-layers on the crest and inner slope of dikes[R].Delft:Report DC030202/H3803, 2002.

[12]代英男.斜坡堤越浪流厚度及堤后次生波的試驗(yàn)研究[D].大連：大連理工大學(xué)，2011.

DAI Ying-nan.Experimental research on layer thickness of wave overtopping flow and secondary wave conditions over sloping breakwater[D].Dalian:Dalian University of Technology,2011.

[13]朱偉娜.規(guī)則波作用下海堤越浪流特征的實(shí)驗(yàn)研究[D].上海：上海交通大學(xué)，2012.

ZHU Wei-na.An experimental study on characteristics of overtopping flow against seadikes in regular wave[D].Shanghai:Shanghai Jiaotong University,2012.

[14]BOSMAN G.Velocity and flow depth variations during wave overtopping[D].Delft:Civil Engineering&Geosciences,2007.

Characteristics of overtopping rate and overtopping flow under the effect of superstandard tide and wave

AN Meng-hua1,JIANG Qin2
（1.CCCC Third Harbor Consultant Co.,Ltd.,Shanghai 200032,China; 2.Hohai University,Nanjing,Jiangsu 210029,China）

In order to study the destructive mechanism of superstandard storm surge caused by typhoon,we studied the overtopping rate and overtopping flow of the sloped seadike with parapet under superstandard tide and wave through physical experiments.The empirical formulas between the relative crest ultra-height,relative wall-top ultra-height,relative parapet height and the overtopping rate were established based on experimental data of regular waves.In addition,the appropriate dimensionless parameter was proposed for the research on the relation between overtopping flow thickness and overtopping rate. The proper relative position parameter was also defined to study the overtopping flow thickness on the dike crest.

superstandard tide;superstandard wave height;sloped seadike;parapet;overtopping rate;overtopping flow

U656.2；P731.22

A

2095-7874（2017）06-0011-06

10.7640/zggwjs201706003

2016-10-27

2016-12-26

安蒙華（1989—），男，浙江舟山人，碩士，助理工程師，主要從事港口及航道設(shè)計(jì)工作。E-mail：anmh@theidi.com