涂俊豪，陶愛峰，蘇俊瑋，徐 嘯，張海明

（1.海岸災(zāi)害及防護教育部重點實驗室（河海大學(xué)），江蘇 南京 210024；2.河海大學(xué) 港口海岸與近海工程學(xué)院，江蘇 南京 210024）

0 引言

由于能源短缺問題日益嚴重，大力發(fā)展可再生能源成為世界各國的共識。波浪能具有環(huán)保可再生、能量密度高等優(yōu)勢，開發(fā)利用波浪能，成為解決能源問題的重要突破口。人們對波浪能利用的研究已有百年歷史，大多數(shù)學(xué)者通過對波能發(fā)電裝置的改進和陣列布局優(yōu)化來提高發(fā)電裝置的發(fā) 電 效 率[1]～[5]。近 年 來，有 學(xué) 者 發(fā) 現(xiàn)，可 以 通 過 布 拉格共振來提高波浪能量密度，從而來提高發(fā)電裝置的發(fā)電效率。

1915年，布拉格父子發(fā)現(xiàn)，當X射線的波長為晶體間距的兩倍時，X射線的反射最強，這種波列的相關(guān)干涉，被稱為布拉格共振[6]。Davies A G[7]發(fā)現(xiàn)，特定波長的表面重力波會與連續(xù)的正弦周期地形發(fā)生共振，即使沙波波高不大，表面重力波的振幅也會在發(fā)生共振以后大幅增加。Elandt R B[8]發(fā)現(xiàn)，可以借助布拉格共振驅(qū)動雙層密度流體的表面波能量向界面波傳遞，從而保護浮式建筑物；同時，Elandt R B提出，可以利用布拉格共振，讓表面重力波的能量匯聚，從而提高波浪能量密度和發(fā)電裝置的發(fā)電效率。Tao A[9]通過波浪水槽試驗發(fā)現(xiàn)，利用布拉格共振來提高波浪能量的密度，可以將發(fā)電裝置的發(fā)電功率提升10倍以上。蘇俊瑋[10]基于高階譜方法發(fā)現(xiàn)，當?shù)匦纬蔞形布局的連續(xù)正弦時，波浪與地形發(fā)生布拉格共振反射后，地形前方會產(chǎn)生聚焦效應(yīng)，波浪振幅增幅較二維地形顯著增大。但是，該研究只考慮了共振最強狀態(tài)對應(yīng)的入射波布拉格共振引發(fā)的波能聚焦特性，而真實的波浪由不同頻率的波組成。

因此，本文采用文獻[10]中的數(shù)學(xué)模型進一步研究了在二維布拉格共振的共振主頻附近，V形布局地形在不同頻率入射波下 （保持地形波陡與入射波波陡不變）的布拉格共振聚焦特性，研究的成果可為布拉格共振發(fā)電的實際應(yīng)用提供理論參考。

1 波浪與地形的耦合模型及驗證

利用不同頻率的規(guī)則波進行正向入射和斜向入射（以 θ=19.47°為例，其中 θ為地形法線和入射波的夾角），以復(fù)演文獻[11]中的兩組布拉格共振數(shù)值模擬試驗，并求得不同頻率下的反射系數(shù)R，將模擬結(jié)果與文獻[11]的攝動理論結(jié)果進行對比，結(jié)果見圖1。從圖1可以看出，數(shù)值模擬結(jié)果與理論曲線吻合良好，證實了模型計算結(jié)果的可靠性。

圖1 不同頻率的布拉格共振反射系數(shù)Fig.1 Bragg resonance reflection coefficients at different frequencies

2 數(shù)值模擬

本節(jié)將具體介紹本次模擬的基本參數(shù)、入射波參數(shù)、地形參數(shù)和模擬組次的設(shè)置情況。下文中：kx為入射波x方向上的波數(shù)；kb為地形波數(shù)；T，A，L和k分別為入射波的周期、振幅、波長和波數(shù)；d為地形振幅；h為水深。

2.1 模型基本參數(shù)

模型計算區(qū)域的大小為32 m×32 m，共512×512個計算節(jié)點。按每周期T/Δt=64的時間步長進行模擬，共模擬波浪傳播TS/T=20個周期。相關(guān)參數(shù)按照文獻[11]進行設(shè)置，非線性階數(shù)M=3。

2.2 入射波設(shè)置

在本文中，入射波的初始波面和勢函數(shù)按微幅波理論進行計算，同時保證各組頻率下入射波的波陡相等，即kA=0.05。按長度轉(zhuǎn)換無量綱參數(shù)L′/L=2π/32，將 所 有 參 數(shù) 無 量 綱 化。

2.3 地形設(shè)置

根據(jù)文獻[10]的研究，各組頻率入射波對應(yīng)的V形 布 局 地 形 的 夾 角 α為90～180°。d，h，kb滿 足d/h=0.31，kbd=0.16。圖2為V形地形在空間布局的立面圖（總體以y=16 m對稱，以V形地形的頂點為x方向的坐標原點，并向著來波方向增大，ξ為地形高度）。

圖2 起伏地形空間布局立面圖Fig.2 Stereoscopic schematic diagram of rippled bottom spatial layout

本研究通過改變2kx/kb來改變?nèi)肷洳ǖ念l率，考慮到大部分海浪譜為單峰譜，且譜寬較窄，所 以 當2kx/kb＞1.1或2kx/kb＜0.9時，布 拉 格 共 振 的波浪反射率 （反射率的大小等于反射波高與入射波高的比值）雖然會增大但仍較小，而且對應(yīng)頻率的入射波能量較小。因此，本文研究的頻率為2kx/kb=0.9～1.1。

2.4 模擬組次

綜合前文討論，本文研究的入射波頻率為2kx/kb=0.9～1.1，從中均勻取17個頻率，研究的V形布局地形的夾角α為90～180°，從中取15個角度，共構(gòu)成255組模擬組次。

3 結(jié)果與分析

入射波會與連續(xù)正弦地形發(fā)生共振反射，反射波經(jīng)由V形布局地形匯聚，在地形前方與入射波疊加，形成不完全立波。振幅顯著增大的點組成的區(qū)域，被稱為波能聚焦區(qū)。20個周期內(nèi)有效波高Hs與初始波高H比值（Hs/H）最大的點被稱為波能聚焦點。通過對比各模擬組次波能聚焦點的Hs/H，找出各頻率入射波對應(yīng)V形布局地形的最優(yōu)夾角，并探究最優(yōu)V形布局地形夾角下，各頻率入射波的波能聚焦區(qū)分布及波能聚焦點的運動特性。

3.1 入射波對應(yīng)V形布局地形的最優(yōu)夾角

根據(jù)模擬結(jié)果，計算出模擬計算區(qū)域內(nèi)512×512個計算點20個周期內(nèi)的有效波高，各模擬組次波能聚焦點的Hs/H如圖3所示（為了清晰地表示波能聚焦點的Hs/H與頻率和V形布局地形夾角 的 關(guān) 系，以2kx/kb=0.9，0.95，1，1.05，1.1的5組頻率入射波的模擬組次結(jié)果為例畫圖）。從圖3可以看出：隨著V形布局地形夾角 α的逐漸增大，波能聚焦點的Hs/H先增大后減??；當2kx/kb偏離1較少時，波能聚焦點的Hs/H最大值對應(yīng)的 α=162.24°；當2kx/kb偏離1較多時，波能聚焦點的Hs/H最大值對應(yīng)的 α不等于162.24°。

圖3 波能聚焦點的Hs/H隨α的變化曲線Fig.3 The variation of Hs/H at wave energy focusing point as a function ofα

不同 α下，波能聚焦點的Hs/H隨入射波頻率的變化曲線如圖4所示。從圖4可以看出：當α=90°時，入射波的反射率為零，各頻率入射波的波能聚焦點的Hs/H均很??；在不同 α下，波能聚焦點的Hs/H均隨著頻率的增大而呈現(xiàn)出先增大后減小的變化趨勢。結(jié)合圖1和文獻[10]的研究可知，當2kx/kb偏離1較大或 α接近90°時，波浪的反射率較小，布拉格共振作用較弱，波能聚焦點的Hs/H均較小，模擬結(jié)果易受到水深變化引起的波浪非線性作用等因素的影響，因此，模擬結(jié)果不能較好地反映布拉格共振的特性。由此可以得到結(jié)論，在本文的研究頻率范圍內(nèi)，各頻率入射波對應(yīng)V形布局地形的最優(yōu)夾角均相等，且α=162.24°。因此，本文將進一步研究最優(yōu)V形布局地形夾角下，不同頻率入射波的布拉格共振特性。

圖4 波能聚焦點的Hs/H隨入射波頻率的變化曲線Fig.4 The variation of Hs/H at wave energy focusing point as a function of incident wave frequency

3.2 入射波的波能聚焦區(qū)的Hs/H分布

通過統(tǒng)計計算區(qū)域內(nèi)各計算點20個周期內(nèi)的波動數(shù)據(jù)，得到各計算點的有效波高，將各點的有效波高與初始波高進行對比，可知各點的Hs/H如圖5所示。從圖5可以看出：當2kx/kb偏離1較大時，波高增大區(qū)域主要集中在地形上方；當2kx/kb接近1時，布拉格共振作用較強，大量反射波由V形布局地形聚焦后，在地形前方形成波能聚焦區(qū)。由于地形對稱布置，發(fā)生布拉格共振反射時，有效波高分布也以y=16 m為對稱軸對稱分布。波浪由V形布局地形反射以后，能量向?qū)ΨQ軸聚焦，因此對稱軸上的波能聚焦區(qū)的Hs/H最大，并向著兩側(cè)逐漸降低，波能聚焦點的位置在對稱軸上。因此，本文將繼續(xù)研究對稱軸上的波能聚焦點的Hs/H分布。

圖5 不同頻率入射波的波能聚焦區(qū)的Hs/H分布Fig.5 Spatial distribution of Hs/H within the wave energy focusing region at different incident wave frequencies

在對稱軸上-2.5～22.5 m的范圍內(nèi)，各計算點的有效波高如圖6所示（以地形結(jié)束位置為0，地形前方為正）。從圖6可以看出：當2kx/kb接近1時，布拉格共振反射率較大，地形前方有一系列波高明顯增大且Hs/H相近的點，這些點分布在對稱軸上5 m（對稱軸上地形起始邊界位置）后的第1～15個波峰點之間；當2kx/kb偏離1較大時，布拉格共振強度和反射率均較低，地形前方的Hs/H沒有明顯增大，整個計算區(qū)域內(nèi)Hs/H最大點在地形上方，這是因為水深變化而產(chǎn)生的波浪非線性作用對波高的影響較大，故本文不作分析。由此推測，若入射波是不規(guī)則波，則波能聚焦點在V形布局對稱軸上x=5～12 m處。

圖6 對稱軸上-2.5～25 m之間的計算點的Hs/H分布Fig.6 Hs/H distribution on symmetrical axis between-2.5～25 m

3.3 波能聚焦點的位置及運動特性

當α=162.24°時，不同頻率入射波的波能聚焦點的位置如圖7所示（以地形結(jié)束位置為0，地形前方為正）。

圖7 最優(yōu)V形布局地形夾角下，不同頻率入射波的波能聚焦點的位置Fig.7 The positions of wave energy focusing point at different incident wave frequencies under the condition of optimal V-shaped bottom angle

從圖7可以看出：當2kx/kb＜1時，隨著2kx/kb逐漸增大，波能聚焦點的Hs/H逐漸增大，波能聚焦點的位置逐漸從地形上方向地形前方移動；當2kx/kb＞1時，隨著2kx/kb逐漸增大，波能聚焦點的Hs/H逐漸減小，波能聚焦點的位置逐漸靠近地形；當0.925≤2kx/kb＜1.05時，波能聚焦點在對稱軸上x=10 m附近，此時波能聚焦點的Hs/H最大，且波能聚焦點的位置在地形前方相差不大。因此，下文將繼續(xù)討論該頻域內(nèi)波能聚焦點的運動特性（以2kx/kb=0.95，1和1.05時 為 例）。

圖8為波能聚焦點的波面時間序列圖（η為波面高程）。當波能聚焦點處于穩(wěn)定狀態(tài)時，每相鄰周期的波峰值不相等。因此，本文選取20個周期內(nèi)各計算點的有效波高來反映各點的能量大小。波能聚焦點達到穩(wěn)定狀態(tài)所需時間與穩(wěn)定狀態(tài)的Hs/H值成正比，穩(wěn)定狀態(tài)的波高增幅越大，達到穩(wěn)定所需的時間越長。

圖8 波能聚焦點的波面時間序列Fig.8 Time series of free-surface elevation at wave energy focusing point

4 結(jié)論

本文基于高階譜方法，討論了不同V形布局地形對不同頻率入射波的影響特性，找出了各頻率入射波對應(yīng)的最優(yōu)V形布局地形夾角，并探究了最優(yōu)V形布局地形夾角下，各頻率入射波的波能聚焦區(qū)分布及波能聚焦點的運動特性，得到以下結(jié)論。

①對于不同頻率的入射波，隨著V形布局地形夾角 α的增大，波能聚焦點的Hs/H呈現(xiàn)出先增大后減小的變化趨勢，并且當 α=162.24°時，各頻率入射波的波能聚焦點的Hs/H均達到最大，因此，各頻率入射波的最優(yōu)V形布局地形夾角為162.24°。

②在最優(yōu)V形布局地形夾角下，各頻率入射波的波能聚焦區(qū)均對稱分布，且對稱軸上的波能聚焦區(qū)計算點的Hs/H最大。當2kx/kb偏離1較大時，波能聚焦區(qū)主要分布在地形上方；當2kx/kb接近1時，波能聚焦區(qū)主要分布在地形前方。

③當2kx/kb接近1時，在波能聚焦區(qū)對稱軸上存在一系列增幅明顯增大且增大倍數(shù)相近的計算點，這些計算點分布在15個立波波長范圍內(nèi)，且波能聚焦點位于x=10 m附近；當2kx/kb偏離1較大時，波能聚焦點位于地形上方。

本文只研究了V形布局地形下不同頻率規(guī)則波的布拉格共振特性，不規(guī)則波可以看成由多組不同頻率規(guī)則波按不同海浪譜組合而成。當入射波為不規(guī)則波時，V形布局地形上的布拉格共振特性，還有待探討。