程永舟，李青峰，潘昀，韓二品，王晶

（1.長沙理工大學 水利工程學院，湖南 長沙 410004；2.水沙科學與水災害防治湖南省重點實驗室，湖南 長沙 410004）

1 引言

波浪循環(huán)荷載對海床床面形態(tài)以及孔隙水壓力產(chǎn)生影響，破碎波對海床的影響更大，研究破碎波作用下海床響應這一問題對于海岸工程具有重要的理論價值以及實際意義。

波浪由深海向淺海傳播過程中自身非對稱性逐漸增強，當局部能量過于集中時會導致破碎，關于波浪破碎的判定，目前最具有代表性的是合田指標［1］與Nelson指標［2］。波浪的破碎形態(tài)可以區(qū)分為崩破波、卷破波、塌破波以及激破波［3］。波浪作用下的泥沙運動以及海床響應問題在早期就引起國內(nèi)外學者的重視，到現(xiàn)在已經(jīng)日趨成熟。蔣昌波等［4］對海嘯波作用下的岸灘剖面演變規(guī)律進行了研究，Young等［5］通過水槽實驗，研究了孤立波破碎帶床面形態(tài)的演變以及海床中的孔隙水壓力。Seelam等［6］研究了未破碎孤立波作用下的床面剪切力。鐘佳玉等［7］研究了規(guī)則波與不規(guī)則波作用時海床孔隙水壓力的變化規(guī)律，得出了規(guī)則波作用下不同波浪基本要素對孔隙水壓力的影響。程永舟等［8—9］通過物理模型實驗，研究了破碎波作用下床面形態(tài)變化以及海床中孔隙水壓力響應問題，同時分析了波浪參數(shù)的變化對孔隙水壓力的影響。破碎帶水體運動機理十分復雜，本文著重對不同波浪作用下的床面形態(tài)以及孔隙水壓力進行分析，同時探討不同破碎波類型對海床響應的影響大小。

2 實驗設計

2.1 實驗裝置

本次實驗在長沙理工大學港航中心的波浪水槽中進行，水槽兩側為透明玻璃，總長40m，寬0.5m，高0.8m。試驗中造波機采用大連理工大學研制的造波控制系統(tǒng)，通過攝像機對波浪以及床面形態(tài)進行采集，孔隙水壓力量測采用天津水運工程科學研究所研制的SG2008型水工實驗數(shù)據(jù)采集處理系統(tǒng)，浪高量測采用WG－50浪高儀和優(yōu)泰克采集系統(tǒng)。

2.2 模型設計

試驗中采用中值粒徑為0.47mm的泥沙在水槽中鋪設坡度為1∶20的海床模型。首先根據(jù)預備實驗確定波浪的破碎區(qū)，然后根據(jù)波浪破碎位置確定浪高儀、孔隙水壓力傳感器位置，不同波況根據(jù)破碎位置的不同而調(diào)整上述采集儀器的位置。圖1所示的實驗布置圖是根據(jù)入射波高為10cm，周期為2s，水深為35cm的規(guī)則波情況布置的，其中布置7個浪高儀，在6＃斷面埋置1組孔隙水壓力傳感器，該組垂向布置4個孔隙水壓力傳感器，其中第一個壓力傳感器置于斜坡表面，豎直向下每隔5cm布置一個，壓力傳感器從上至下標號為P1、P2、P3、P4。

圖1 實驗布置示意圖

2.3 波浪參數(shù)

試驗中，不同的波浪破碎時有可能產(chǎn)生不同類型的破碎波，不同類型的破碎波對海床中孔隙水壓力的影響不同，即使相同波況的波浪在破碎時也可能導致不同類型的破碎。波浪破碎類型可以由下式確定：

式中，β為斜坡坡度；H0為深水波高；L0為深水波長，其中，

式中，T為周期，k為波數(shù)，h為水深。

當0.46＜ξ0＜3.3時，波浪破碎類型為卷破波，當ξ0＜0.46時，波浪破碎類型為崩破波，當ξ0＞3.3時，波浪破碎類型為激破波。試驗中的波浪采用孤立波和規(guī)則波，選取部分波況見表1。

表1 波況參數(shù)

3 實驗結果分析

3.1 孤立波和規(guī)則波對海床的影響

破碎波作用下床面形態(tài)變化劇烈，然而孤立波和規(guī)則波破碎時對床面形態(tài)的影響完全不同，對比分析了水深35cm、波高10cm的孤立波和水深35cm、波高10cm、周期2s的規(guī)則波兩組波況破碎時對床面形態(tài)的影響，其波形歷時曲線見圖2、3，從圖中可以看出，波浪從深海向淺海傳播過程中，波峰逐漸尖銳，波高逐漸增大至破碎時達到最大值，破碎后波高變得很小。

3.1.1 床面形態(tài)分析

圖4描述的是孤立波和規(guī)則波造波20次后破碎帶地形，從圖中可以得知，孤立波和規(guī)則波對床面形態(tài)影響不同，孤立波作用下，床面在波浪破碎前形成沙壩，在破碎處形成沙坑，而在規(guī)則波作用下，床面在破碎處形成沙壩，在破碎后形成沙坑。

圖2 孤立波歷時曲線

圖3 規(guī)則波歷時曲線

圖4 破碎帶床面形態(tài)

通過影像采集分析了孤立波造波次數(shù)以及規(guī)則波造波時間對床面形態(tài)的影響，從圖5、6可以得知。隨著造波次數(shù)的增加，沙層堆積和掏刷程度都會有所增加，然而與規(guī)則波相比，孤立波破碎段較長，導致了床面形態(tài)受孤立破碎波影響范圍更大，且形成的沙壩和沙坑尺度較大。

3.1.2 孔隙水壓力分析

孤立波和規(guī)則波在向岸傳播過程中發(fā)生破碎，但波浪破碎時對海床孔隙水壓力的影響不同，見圖7，將相同水深和波高的孤立波和規(guī)則波破碎斷面的孔隙水壓力進行對比分析，發(fā)現(xiàn)孤立波破碎處不同深度處的孔隙水壓力相比規(guī)則波破碎處的孔隙水壓力較大，且規(guī)則波破碎處孔隙水壓力沿深度方向衰減較快。另外將表層P1值相同的兩組波況進行對比分析，從圖8中可以得知，兩者在床面表層的孔隙水壓力大小相同，但是隨著深度的增加，孤立波破碎處孔隙水壓力幅值衰減較慢，而規(guī)則波破碎處孔隙水壓力幅值衰減較快。

3.2 波面變化對孔隙水壓力的影響

規(guī)則波在破碎前，波形歷時變化曲線很光滑，而在破碎時波面變化不穩(wěn)定，為了研究波面變化對海床不同深度孔隙水壓力的影響，對比分析了相同水深和波高的孤立波和規(guī)則波破碎時的波面歷時曲線和不同深度處的孔隙水壓力歷時曲線。通過圖9可以得知，不同海床深度的孔隙水壓力隨著波面的升高而增大，隨著波面的降低而減小，波面變化的歷時曲線與孔隙水壓力歷時曲線相似。然而，孤立波在爬坡和回落階段（圖9a）以及規(guī)則波箭頭所示波谷段（圖9b），波面變化較小，而不同深度處的孔隙水壓力卻持續(xù)下降，表明波面的變化不能完全反應孔隙水壓力的變化規(guī)律。另外，波面的變化對不同海床深度的孔隙水壓力影響程度不同，波面變化對淺層的孔隙水壓力變化規(guī)律影響程度更大，波面變化歷時曲線與淺層孔隙水壓力變化歷時曲線最相似，與深層的孔隙水壓力歷時曲線相似度較低。

圖5 孤立波作用下床面變化

圖6 規(guī)則波作用下床面變化

圖7 相同波高時孔隙水壓力幅值

圖8 相同P1值時孔隙水壓力幅值

圖9 破碎波波面與孔隙水壓力歷時曲線

同時研究了波浪破碎時波面變化形態(tài)對孔隙水壓力梯度的影響，將P1－P2、P2－P3、P3－P4所得數(shù)據(jù)的歷時曲線與波面歷時曲線進行對比分析，見圖10。從圖中可以看出，在波峰段，海床淺層的孔隙水壓力梯度值最大，沿著深度方向，孔隙水壓力梯度值逐漸減小。在孤立波爬坡和回落階段（圖10a）以及規(guī)則波波谷段（圖10b），當波面變化不大時，不同深度的孔隙水壓力梯度相差不大。另外，波面變化歷時曲線與孔隙水壓力梯度歷時曲線極其相似，其相似度比波面變化與孔隙水壓力歷時曲線之間的相似度更高。這表明波面的變化更能反映一種壓力梯度的變化。

圖10 破碎波波面與孔隙水壓力梯度歷時曲線

3.3 不同破碎波類型對孔隙水壓力的影響

分析了入射波高為10cm、水深35cm、周期2s的規(guī)則波，通過式（1）計算得知，ξ0為0.295，即波浪開始破碎類型為崩破波，從圖11中可以看出，床面形態(tài)未發(fā)生變化時，波浪破碎類型為崩破波，隨著破碎帶地形發(fā)生變化，形成沙壩及沙坑，波浪破碎類型由崩破波轉變?yōu)榫砥撇?。隨著波浪破碎類型的變化，土體中的孔隙水壓力隨之發(fā)生變化，圖12描述的造波1min前的孔隙水壓力幅值以及造波1～2 min的孔隙水壓力幅值，從實驗結果可以得知，當波浪破碎為崩破波時，土體中孔隙水壓力幅值相對較小，當波浪破碎為卷破波時，土體中孔隙水壓力幅值相對較大。

圖11 地形變化及波浪破碎類型變化

圖12 不同類型破碎波作用下孔壓幅值

4 結論

通過對不同破碎波作用下床面形態(tài)及海床響應的實驗研究，發(fā)現(xiàn)孤立波與規(guī)則波對海床床面形態(tài)變化的影響不同，與規(guī)則波相比，孤立波破碎區(qū)較長，導致了床面形態(tài)受孤立破碎波影響范圍更大，孤立波破碎前床面形成沙壩，破碎處床面形成沙坑，且形成的沙壩和沙坑尺度較大。而規(guī)則波破碎處床面形成沙壩，破碎后床面形成沙坑，且形成的沙壩和沙坑尺度較小。對比分析波面歷時變化曲線與孔隙水壓力歷時變化曲線，發(fā)現(xiàn)兩者變化規(guī)律相似，進一步對波面歷時曲線與孔隙水壓力梯度歷時曲線進行對比分析，發(fā)現(xiàn)波面歷時變化曲線與孔隙水壓力梯度歷時變化曲線更為相似，說明波面變化更能反映海床內(nèi)部孔隙水壓力梯度的變化。破碎波沖擊海床，造成床面形態(tài)發(fā)生變化，然而床面形態(tài)的變化同時對波浪破碎類型產(chǎn)生影響，當床面形成沙壩和沙坑時，破碎類型可由崩破波轉化為卷破波，且卷破波作用下的孔隙水壓力幅值較崩破波作用下的孔隙水壓力幅值大。

［1］ Goda Y.A synthesis of breaker indices［J］.Transactions of Japan Society of Civil Engineers（part 2），1970，2：227－229.

［2］ Nelson R C.Design wave height on very mild slope－An Experimental study［J］.Civil Engineering，1987，29（3）：157－161.

［3］ Galvin C J.Breaker type classification on three laboratory beaches［J］.Journal of Geophysical Research，1968，73（12）：3651－3659.

［4］ 蔣昌波，陳杰，程永舟，等.海嘯波作用下泥沙運動：Ⅰ.岸灘剖面變化分析［J］.水科學進展，2012，23（5）：665－672.

［5］ Young Y L，Xiao H，Maddux T.Hydro－and morpho－dynamic modeling of breaking solitary waves over a fine sand beach，PartⅠ：Experimental study［J］.Marine Geology，2010，269（3/4）：107－118.

［6］ Seelam J K，Gurard P A，Baldock T E.Measurement and modeling of bed shear stress under solitary waves［J］.Coastal Engineering，2011，58（9）：937－947.

［7］ 鐘佳玉，鄭永來，倪寅.波浪作用下砂質(zhì)海床孔隙水壓力的響應規(guī)律試驗研究［J］.巖土力學，2009，30（10）：3188－3193.

［8］ Cheng Y Z，Jiang C B，Zhao L P，et al.Experimental study of pore water pressure and bed profile change under regular breaking waves［J］.China Ocean Engineering，2012，26（3）：457－468.