孫家文，房克照，何棟彬，張振偉

(1.國(guó)家海洋環(huán)境監(jiān)測(cè)中心，大連 116023；2.大連理工大學(xué) 海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室， 大連 116024；3.國(guó)家海洋局海島研究中心，平潭 350400)

隨著人類對(duì)自然資源需求的增長(zhǎng)以及社會(huì)經(jīng)濟(jì)發(fā)展的需要，人類生產(chǎn)活動(dòng)的范圍不斷地向海洋擴(kuò)張，在這一現(xiàn)實(shí)背景下，能為生產(chǎn)作業(yè)提供陸地依托的島礁，無(wú)疑將扮演越來越重要的角色。同時(shí)，島礁本身也具有豐富的自然資源，比如生物資源、礦物資源、海洋能資源。我國(guó)是海洋大國(guó)，海島眾多，面積超過500 m2的海島有6 500多個(gè)，在500 m2以下的海島約有上萬(wàn)個(gè)[1]。

深海島礁地處開敞性海域，自然環(huán)境惡劣，浪、潮和流等水動(dòng)力作用顯著，是島礁工程要考慮的主要環(huán)境荷載。而且島礁工程遠(yuǎn)離大陸，屬遠(yuǎn)海孤島建設(shè)，建設(shè)和維護(hù)成本昂貴[2]，更要充分考慮上述水動(dòng)力因素，如設(shè)計(jì)波要素的確定、波浪循環(huán)荷載對(duì)結(jié)構(gòu)穩(wěn)定性和地基承載力的影響等[3-4]，以保證工程的安全性和持久性。此外，需要明確工程建設(shè)對(duì)水動(dòng)力的影響以防惡化島礁水動(dòng)力環(huán)境，導(dǎo)致島礁的侵蝕和消失。在生態(tài)環(huán)境保護(hù)方面，島礁面積小、承載力和生態(tài)系統(tǒng)脆弱。我國(guó)深水島礁多為珊瑚島礁，珊瑚以及附屬生物群落對(duì)生長(zhǎng)環(huán)境要求苛刻。而人類活動(dòng)的增加和深海資源開采不可避免帶來環(huán)境污染，需要重視上述水動(dòng)力因素作用下物質(zhì)輸運(yùn)問題以加強(qiáng)對(duì)深海島礁生態(tài)環(huán)境的保護(hù)。此外，生物群落間營(yíng)養(yǎng)物質(zhì)輸運(yùn)的研究也依賴于對(duì)水動(dòng)力的了解?？梢灶A(yù)見，在將來很長(zhǎng)一段時(shí)間內(nèi)，島礁工程和生態(tài)環(huán)境保護(hù)所涉及的島礁水動(dòng)力問題會(huì)日益突出，亟需開展相關(guān)研究以提供理論和技術(shù)支持。

本文結(jié)合國(guó)內(nèi)外學(xué)者的相關(guān)研究，具體闡述島礁水動(dòng)力環(huán)境的復(fù)雜之處和目前研究關(guān)注的熱點(diǎn)與難點(diǎn)。同時(shí)，對(duì)現(xiàn)有的現(xiàn)場(chǎng)觀測(cè)實(shí)驗(yàn)、物理模型實(shí)驗(yàn)和數(shù)值模型進(jìn)行總結(jié)歸納和介紹，并指出當(dāng)前研究存在的不足。

1 島礁地貌特征

島礁附近海底條件復(fù)雜，當(dāng)波浪從深遠(yuǎn)海傳播到島礁附近，由于其復(fù)雜的地形地貌，波浪將發(fā)生折射、繞射和反射等，當(dāng)波浪變陡到一定限度后，便產(chǎn)生破碎?？梢哉f，島礁地貌的特殊性決定了其周圍水動(dòng)力環(huán)境的復(fù)雜性，而島礁地貌的特殊性主要體現(xiàn)在下述兩個(gè)方面。

圖1 島礁地形示意圖Fig.1 Schematic of a reef topography

一是具有極其陡峭的礁前斜坡以及多級(jí)階梯狀地形。典型的深水島礁地形見圖1，通常由以下幾部分組成，前礁、礁冠、礁坪、瀉湖和沙島。前礁一端連接外海，一端連接礁冠，水深急劇變化，形成陡峭的礁前斜坡。另外，島礁存在的多級(jí)階梯狀地形也是其地貌特殊性的一個(gè)體現(xiàn)。在水深600～850 m范圍內(nèi)普遍發(fā)育水下階地，水下階地多達(dá)5級(jí)，最常見的坡度為0.5°～3°和9°～11°[5-6]。

二是發(fā)育有粗糙度極大而其水深很淺的礁坪。礁坪上水深很淺，低潮時(shí)甚至露出水面，對(duì)于岸礁，礁坪分布較廣，可延伸至沙島(或沙洲)，對(duì)于堡礁和環(huán)礁，礁坪寬度變窄，其后方為瀉湖。礁坪—瀉湖受到礁冠的掩護(hù)作用，波能相對(duì)較低，往往有生物群落發(fā)育，生物成固的珊瑚礁巖土具有易破碎、多空隙、高壓縮等特性，珊瑚礁體表層主要是由鈣質(zhì)沙形成的松散沙礫層，鈣質(zhì)沙孔隙比范圍為0.54～2.97，遠(yuǎn)高于石英砂的通常范圍(0.4～0.9)[3],這相當(dāng)于存在粗糙度極大且可滲透的底床。

上述兩方面共同構(gòu)成了島礁地形的復(fù)雜性。島礁地形上的水動(dòng)力問題是眾多學(xué)者和工程師共同關(guān)注的焦點(diǎn)并取得了許多研究成果，下面將介紹國(guó)內(nèi)外學(xué)者在該領(lǐng)域的研究方向及相關(guān)工作，包括島礁水動(dòng)力的物理實(shí)驗(yàn)?zāi)Ｐ秃蛿?shù)值計(jì)算模型研究情況。

2 主要研究方向

波浪在深海島礁獨(dú)特的地形上傳播時(shí)，發(fā)生的變形以及相關(guān)水動(dòng)力現(xiàn)象極其復(fù)雜。國(guó)際上的研究相對(duì)較多，研究?jī)?nèi)容更為廣泛，主要集中在以下幾個(gè)方面：

圖2 Nwogu和Demirbilek實(shí)驗(yàn)設(shè)置圖[9]Fig.2 Experiment setup of Nwogu and Demirbilek[9]

一是對(duì)波浪能量損耗和礁坪增水的研究。外海波浪入射至島礁區(qū)域，礁前斜坡和礁冠處水深急劇變淺，波浪發(fā)生劇烈破碎，礁坪上水深較淺且粗糙度極高，進(jìn)一步增強(qiáng)了對(duì)波浪的衰減作用。Lugo等[7]在Margarita礁的現(xiàn)場(chǎng)觀測(cè)結(jié)果表明，通過礁坪后波浪有效波高減小了82%，能量衰減達(dá)96%。Ferrario等[8]分析了包括亞特蘭大、太平洋和印度洋等區(qū)域的珊瑚礁的監(jiān)測(cè)數(shù)據(jù)，指出珊瑚礁對(duì)波能的削弱平均達(dá)到了97%，而其中礁冠貢獻(xiàn)最大，消耗了86%的能量。島礁上波高的急劇衰減，導(dǎo)致輻射應(yīng)力的梯度出現(xiàn)變化，引起礁坪平均水面的上升——礁坪增水。由于礁坪水深很淺，增水不容忽視，其值甚至與水深相當(dāng)。波浪能量衰減和礁坪增水是對(duì)立統(tǒng)一的兩個(gè)方面，從產(chǎn)生機(jī)制上看，前者是后者產(chǎn)生的原因，從作用效果上看，前者對(duì)島礁工程有利而后者不利，均需要引起重視。對(duì)這兩個(gè)問題的理論與數(shù)值分析主要基于輻射應(yīng)力理論、波能損耗方程和相關(guān)的波浪計(jì)算模型[9-10]。也有一些波高衰減的估算公式提出，其中Yao等[10]基于一維完全非線性Boussinesq方程，針對(duì)波浪通過岸礁引發(fā)的增減水和波高變化問題建立了相應(yīng)的數(shù)值模型。物理模型實(shí)驗(yàn)則多圍繞簡(jiǎn)單島礁地形上的一維波浪傳播問題展開[11-12]。Nwogu和Demirbilek[9]參照關(guān)島東南側(cè)島礁地形的典型剖面設(shè)計(jì)了一個(gè)二維斷面模型對(duì)不規(guī)則波進(jìn)行實(shí)驗(yàn)(圖2)，對(duì)波浪的能量譜進(jìn)行了分析。Quiroga和Cheung[13]開展了一系列二維水槽實(shí)驗(yàn)來分析底床粗糙度對(duì)孤立波傳播的影響，并與普通混凝土床面的實(shí)驗(yàn)結(jié)果做了比對(duì)。目前的研究主要探討波高衰減和礁坪增水與島礁地形特征參數(shù)的關(guān)系。

二是對(duì)低頻波浪運(yùn)動(dòng)的研究。根據(jù)Van Dongeren[14]的工作，島礁上劇烈的波浪破碎一方面導(dǎo)致短波(周期為5～25 s)能量的損耗，另一方面波浪破碎點(diǎn)的移動(dòng)誘發(fā)低頻波浪運(yùn)動(dòng)(周期為25 s～10 min)。此外，在礁冠處未發(fā)生破碎的、波高較小的涌浪在礁坪和瀉湖傳播時(shí)，波—波非線性作用也導(dǎo)致低頻波浪的產(chǎn)生[9,15]。礁坪—瀉湖上的低頻波浪運(yùn)動(dòng)對(duì)島礁水動(dòng)力影響顯著，主要包括：1)容易導(dǎo)致波浪爬高的增大，島礁上發(fā)生的極端波浪災(zāi)害事件與低頻波浪相關(guān)[9]；2)低頻波動(dòng)頻率與礁坪—瀉湖固有頻率一致時(shí)，易誘發(fā)水體共振[9]，是島礁碼頭建設(shè)需要考慮的問題；3)低頻波浪的存在顯著改變礁坪—瀉湖流體波動(dòng)和流場(chǎng)形態(tài)[16]。相比于波浪能量損耗和礁坪增水，近幾年才開始重視對(duì)島礁低頻波浪運(yùn)動(dòng)的研究，相關(guān)的物理模型實(shí)驗(yàn)[9-10]和數(shù)值計(jì)算結(jié)果[9-10,14-15,17]較少。

圖3 礁坪—瀉湖體系水循環(huán)[22]Fig.3 Wave-driven circulation of a reef-lagoon system[22]

三是對(duì)波生流的研究。傳統(tǒng)意義上的波生流是波浪在近岸破碎后引起的時(shí)均水流，如沿岸流和裂流，而島礁上的波生流又增加了一個(gè)特殊的組成部分，即由礁坪增水引起的流向礁坪—瀉湖體系的水流運(yùn)動(dòng)?；谳椛鋺?yīng)力理論的波生流模式表明輻射應(yīng)力變化和平均水面變化是導(dǎo)致波生流的原因[11]，而礁坪上波高變化不大，波生流的產(chǎn)生原因主要為礁坪增水。礁坪增水引起的水流向岸運(yùn)動(dòng)是礁坪—瀉湖體系水循環(huán)(圖3)的重要來源，對(duì)礁坪—瀉湖體系的生態(tài)系統(tǒng)有重要的意義，因此得到了廣泛的關(guān)注[18]，但這些研究偏重于礁坪—瀉湖體系的流量問題。對(duì)島礁波生流流速和流場(chǎng)特征的研究罕有報(bào)道[14,19]。

國(guó)內(nèi)對(duì)深海島礁水動(dòng)力的研究比較薄弱，處于起步階段。這主要是因?yàn)槲覈?guó)過去的海岸和海洋工程僅局限于海岸和近海區(qū)域，很少涉及到深海島礁。目前的研究結(jié)果主要是對(duì)島礁地形上波能衰減和波浪破碎指標(biāo)進(jìn)行研究。黎滿球等[20]使用珊瑚礁坪上的海浪測(cè)量數(shù)據(jù)，應(yīng)用統(tǒng)計(jì)方法得到波高、波能衰減系數(shù)與相對(duì)水深關(guān)系的擬合方程。物理模型實(shí)驗(yàn)[21-23]一般都是用簡(jiǎn)單斜坡地形模擬島礁地形，分析波高、波能的衰減和波浪破碎情況，其中，姚宇等[22]考慮了波浪水流共存情況。柳淑學(xué)等[23]采用斷面物理模型實(shí)驗(yàn)，對(duì)比分析了規(guī)則波和不規(guī)則波在島礁地形上的傳播過程、破碎指標(biāo)、破碎后波高的衰減，在此基礎(chǔ)上，柳淑學(xué)等[24]進(jìn)行了三維波浪在島礁地形上傳播破碎特性試驗(yàn)研究。上述物理實(shí)驗(yàn)都是在光滑底床條件下開展的。

數(shù)值計(jì)算方面，目前國(guó)內(nèi)多采用相位識(shí)別模型，如淺水方程模型、緩坡方程模型、Boussinesq模型，其中以Boussinesq類波浪模型最為典型，它既包含色散性又包含非線性、能夠描述近岸區(qū)域波浪傳播變形、破碎、爬坡、波流共存、強(qiáng)非線性流體運(yùn)動(dòng)等諸多動(dòng)力因素。劉海清等[25]利用基于NS方程的數(shù)值水槽將對(duì)島礁地形簡(jiǎn)化為臺(tái)階[26]對(duì)非破碎規(guī)則波進(jìn)行了簡(jiǎn)單計(jì)算，田超[4]等利用緩坡方程計(jì)算島礁周圍波浪場(chǎng)，姚宇等[27]采用OpenFOAM數(shù)值模擬珊瑚島礁附近孤立波傳播變形，方亞冰等[28]采用自適應(yīng)有限元求解緩坡方程所建立的數(shù)值模型對(duì)島礁地形上波浪破碎進(jìn)行了模擬研究。房克照等[29-30]建立了基于高階Boussinesq水波方程的波浪傳播數(shù)學(xué)模型，模擬了孤立波和隨機(jī)波浪在潛礁地形上傳播變形，數(shù)值結(jié)果同實(shí)驗(yàn)數(shù)據(jù)吻合較好，尤其是實(shí)現(xiàn)了對(duì)海岸動(dòng)邊界和波浪破碎的處理。對(duì)于如何在數(shù)值模型中考慮底摩擦的影響，目前常見的做法是在動(dòng)量方程中添加底摩擦作用項(xiàng)，該項(xiàng)一般與流速相聯(lián)系，同時(shí)引入底摩擦系數(shù)等可調(diào)系數(shù)。

3 當(dāng)前研究的局限和不足

縱觀上述研究，對(duì)島礁海域的近岸水動(dòng)力研究主要有三種手段：現(xiàn)場(chǎng)實(shí)驗(yàn)觀測(cè)、物理模型實(shí)驗(yàn)和數(shù)值模型計(jì)算。圍繞這三種主要手段，國(guó)內(nèi)外都可以看到相應(yīng)工作的報(bào)道，已經(jīng)取得了一些成果，但仍然存在著局限和不足，主要集中在以下幾個(gè)方面：

圖4 三維島礁示意圖和測(cè)點(diǎn)分布[39]Fig.4 Schematic of a three-dimensional and measurement locations[39]

一是實(shí)驗(yàn)方面，現(xiàn)有的物理實(shí)驗(yàn)?zāi)Ｐ筒荒芊从吵鰨u礁地形的復(fù)雜性，同時(shí)缺少現(xiàn)場(chǎng)實(shí)驗(yàn)觀測(cè)數(shù)據(jù)。島礁地形具有陡峭的前坡，水下多級(jí)臺(tái)階和水深極淺的礁坪，有按照實(shí)際島礁地形變化設(shè)置了多級(jí)斜坡的物理模型[9,30]，但島礁水深變化大，由于比尺限制，很難在實(shí)驗(yàn)室內(nèi)開展實(shí)驗(yàn)，目前多采用簡(jiǎn)化的辦法，使用簡(jiǎn)單斜坡或潛礁地形來進(jìn)行實(shí)驗(yàn)，礁坪段也簡(jiǎn)化為平底，沒有高程變化，類似的簡(jiǎn)化地形難以模擬上述島礁地形的復(fù)雜變化，同時(shí)缺少三維島礁地形的模型實(shí)驗(yàn)，可考慮參考實(shí)際島礁地形，建立對(duì)應(yīng)的三維地形來開展研究。現(xiàn)場(chǎng)觀測(cè)實(shí)驗(yàn)方面，與物理模型實(shí)驗(yàn)不同的是，現(xiàn)場(chǎng)觀測(cè)實(shí)驗(yàn)缺少相應(yīng)的規(guī)范指導(dǎo)，另外，其需要面對(duì)復(fù)雜多變的實(shí)際環(huán)境，觀測(cè)結(jié)果受多重因素(地形、水文和氣候等)影響，同時(shí)受不可控制因素(如風(fēng)、浪入射條件)影響大，這些對(duì)現(xiàn)場(chǎng)觀測(cè)實(shí)驗(yàn)的實(shí)驗(yàn)設(shè)計(jì)、數(shù)據(jù)采集和分析處理帶來很大的難度。目前的現(xiàn)場(chǎng)觀測(cè)實(shí)驗(yàn)有關(guān)注島礁海域近岸水動(dòng)力特性的[16,20,32-33]，也有對(duì)底床粗糙度進(jìn)行觀測(cè)分析的[34-36]，但更為多見的是大范圍海域的觀測(cè)記錄，對(duì)于島礁近岸流域，相應(yīng)的現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)仍比較少，尤其是國(guó)內(nèi)更是很少見到相應(yīng)的報(bào)道。所以，這方面也有待于對(duì)島礁海域觀測(cè)實(shí)驗(yàn)的開展，一方面可以從觀測(cè)數(shù)據(jù)入手分析流場(chǎng)，另一方面也為島礁水動(dòng)力數(shù)值模型的驗(yàn)證提供數(shù)據(jù)積累。

圖5 三維島礁上孤立波傳播的瞬時(shí)波面[39]Fig.5 Snapshots of computed water surface for solitary wave transformation over a three-dimensional reef[39]

二是海岸工程中常用的計(jì)算模型難以勝任深海島礁水動(dòng)力的計(jì)算。在地形變化劇烈的島礁區(qū)域，水動(dòng)力具有水深跨度大、不同尺度波動(dòng)(如潮流、短波和底頻波浪)共存、波浪破碎劇烈、水—陸動(dòng)邊界問題突出和非線性作用強(qiáng)等特點(diǎn)，其數(shù)值計(jì)算具有挑戰(zhàn)性。海岸工程中常用的兩類模型為相位平均模型和相位識(shí)別模型，前者通過輻射應(yīng)力考慮短波貢獻(xiàn)，不能充分考慮短波非線性以及長(zhǎng)短波間相互作用，如SWAN模型不能準(zhǔn)確模擬島礁低頻波浪的產(chǎn)生和運(yùn)動(dòng)[37]。相位識(shí)別模型在海岸工程應(yīng)用廣泛，以Boussinesq類模型最為成功，尤其是采用近年發(fā)展的具備間斷捕捉能力的格式后，此類模型在計(jì)算島礁波浪傳播問題時(shí)取得了良好的效果[38-39]。具備間斷捕捉能力的Boussinesq類模型的主要特征是通過控制色散項(xiàng)來捕捉波浪破碎其中，如Fang等[30]通過求解水平二維完全非線性Boussinesq方程，引入兩個(gè)波浪破碎指標(biāo)，建立了具有間斷捕捉能力的數(shù)值模型，并將其用于二維和三維島礁(圖4和圖5)上孤立波傳播的計(jì)算。但必須指出的是，上述采用相位識(shí)別模型的計(jì)算多針對(duì)物理模型實(shí)驗(yàn)尺度的問題展開，用于實(shí)際工程時(shí)，仍會(huì)存在應(yīng)用水深和緩坡假定的限制[14,17]。近年發(fā)展的新型Boussinesq水波方程在一定程度上能彌補(bǔ)上述缺陷，但方程表達(dá)形式復(fù)雜、涉及三次及以上高階導(dǎo)數(shù)，其數(shù)值離散并非易事。此外，Boussinesq類模型本質(zhì)上屬于水平二維模型，不能準(zhǔn)確刻畫流場(chǎng)三維特征[40]；絕大多數(shù)Boussinesq 類模型都建立在無(wú)旋假定的基礎(chǔ)上，不能有效處理破波帶內(nèi)顯著存在的有旋流體運(yùn)動(dòng)[41]。Navier-Stokes 類模型雖然能給出較為詳細(xì)的流場(chǎng)信息，但由于計(jì)算效率的問題當(dāng)前很難用于深海島礁水動(dòng)力問題的研究，相關(guān)的計(jì)算結(jié)果未見報(bào)道。

近年興起的非靜壓波浪模型[17,42]，雖然具備深海島礁水動(dòng)力計(jì)算的潛質(zhì)，但缺乏系統(tǒng)的理論分析和適當(dāng)?shù)挠?jì)算方法。該模式求解的仍為 Navier-Stokes方程(或歐拉方程)，但其假定自由表面為單值函數(shù)，可以方便地通過求解沿水深積分的連續(xù)方程和運(yùn)動(dòng)學(xué)邊界條件得到，比常用的方法(如VOF方法，LEVEL SET方法等)執(zhí)行簡(jiǎn)單、計(jì)算效率高。最近發(fā)展的動(dòng)力學(xué)邊界條件施加方法使得在沿水體方向劃分較少層數(shù)就能獲得較好的色散性能，進(jìn)一步提高了模型的計(jì)算效率[44]。此類模型的性能，包括色散性和非線性是沿水體分層方式，包括層數(shù)和劃分位置和變量布置方式的函數(shù)(和數(shù)值離散方式相關(guān))，僅Bai和Cheung[45]和Zhu等[46]給出最多至三層的非靜壓模型的理論分析結(jié)果，方程性能缺乏系統(tǒng)的理論分析，難以支持這類模型的合理應(yīng)用和計(jì)算結(jié)果的準(zhǔn)確分析。 發(fā)展具備間斷捕捉能力的計(jì)算方法是建立此類模型的另外一個(gè)主要任務(wù)。復(fù)雜的島礁地形導(dǎo)致波浪的運(yùn)動(dòng)異常復(fù)雜，島礁波浪破碎劇烈且具有極其陡峭的波前、水—陸動(dòng)邊界問題突出、礁坪上流動(dòng)狀態(tài)轉(zhuǎn)化劇烈，容易形成類似水躍(或者激波)的潮涌[38]，上述水動(dòng)力現(xiàn)象的準(zhǔn)確計(jì)算對(duì)任何一種波浪計(jì)算模型而言都極具挑戰(zhàn)性，要求數(shù)值格式具備良好的守恒性，能捕捉間斷，但具備此性能的非靜壓計(jì)算模型非常有限[17,43-44]。通常采用具有動(dòng)量守恒性質(zhì)的差分格式，這類格式建立在交錯(cuò)網(wǎng)格上，不便于方程離散和程序編制。將有限體積方法用于非靜壓模型的求解是近年興起的比較新穎的思路[17,42-43,47]，但當(dāng)前此類模型多采用Roe格式或者HLL格式計(jì)算數(shù)值通量，每個(gè)時(shí)間步都需要求解所有網(wǎng)格界面的黎曼問題，計(jì)算效率較低，對(duì)于三維問題更為突出。如何構(gòu)建準(zhǔn)確、高效的數(shù)值通量算法和提高計(jì)算效率，仍是當(dāng)前主要的研究方向。

三是當(dāng)前研究對(duì)島礁海床的大粗糙度特征考慮不夠充分?，F(xiàn)場(chǎng)觀測(cè)表明，島礁床面粗糙度極大，對(duì)波浪產(chǎn)生的阻力相當(dāng)可觀，底摩擦系數(shù)較普通的沙質(zhì)和淤泥質(zhì)海岸大1～2個(gè)量級(jí)。除表面摩擦力外，島礁通常發(fā)育有生物群落，導(dǎo)致床面高地起伏，這一方面會(huì)降低海床上的有效水深，另一方面流體會(huì)在生物群落間流動(dòng)，上述兩個(gè)因素均直接影響流場(chǎng)形態(tài)[15,22]。大粗糙島，礁海床附近流體運(yùn)動(dòng)非常復(fù)雜，涉及冠層流、邊界層和孔隙滲流等[18]。而現(xiàn)有的物理模型實(shí)驗(yàn)多是在較光滑的床面上進(jìn)行的，不能充分考慮大粗糙度海床對(duì)波浪運(yùn)動(dòng)的影響。Lowe 等[48]將大粗糙度海床處理為剛性植被，在物理模型實(shí)驗(yàn)中以小圓柱代替，其對(duì)流場(chǎng)的作用表現(xiàn)為拖曳力，如何確定剛性植被的拖曳力系數(shù)當(dāng)前仍無(wú)統(tǒng)一認(rèn)識(shí)，一般而言為過水面積和植被密度等的函數(shù)，對(duì)于實(shí)際島礁海床難以概化出這幾個(gè)參量[48]。Quiroga和Cheung[13]通過在光滑床面上間隔安置具有一定高度的小木條增加粗糙度研究孤立波傳播問題(圖6和圖7)，提供了一種探索大粗糙海床對(duì)波浪影響的途徑，但沒有考慮實(shí)際島礁海床的滲透性。

圖6 水槽實(shí)驗(yàn)布置圖(Quiroga和Cheung[13])Fig.6 Schematic of wave flume experiments(Quiroga and Cheung[13])圖7 用木條模擬粗糙度實(shí)驗(yàn)布置圖(Quiroga和Cheung[13])Fig.7 Roughness elements configuration and arrangements(Quiroga and Cheung[13])

數(shù)值計(jì)算模型方面，在動(dòng)量方程中添加阻力項(xiàng)是普遍的做法，一般仍然采用沙質(zhì)海床上常用的底摩阻公式計(jì)算床面阻力，對(duì)于物理模型實(shí)驗(yàn)尺度的計(jì)算通常忽略或采用很小的底摩擦系數(shù)[9-10]，對(duì)于有限的實(shí)際島礁波浪計(jì)算[17]，僅依靠底摩擦系數(shù)不能充分反映島礁海床粗糙度大、可滲透等特征，而且關(guān)于所添加的阻力公式的具體表達(dá)形式以及相關(guān)系數(shù)的選取，眾說紛紜，往往根據(jù)具體的物理實(shí)驗(yàn)結(jié)果去調(diào)整，沒有理論依據(jù)，存在比較大的主觀性。如何在理論上、物理模型實(shí)驗(yàn)和計(jì)算模型中考慮島礁海床的大粗糙度仍是當(dāng)前沒有妥善解決的問題。

4 結(jié)語(yǔ)

島礁周圍地形復(fù)雜，水深變化劇烈，致使島礁附近海域的波浪傳播變形十分復(fù)雜，近島礁波浪表現(xiàn)出強(qiáng)烈的非線性特點(diǎn)，大粗糙度、可滲透的珊瑚礁海床使得復(fù)雜程度進(jìn)一步加劇。研究島礁海域的近岸水動(dòng)力特性主要有三種手段：現(xiàn)場(chǎng)實(shí)驗(yàn)觀測(cè)，物理模型實(shí)驗(yàn)和數(shù)值模型計(jì)算。

開展島礁海域近岸水動(dòng)力特性的現(xiàn)場(chǎng)觀測(cè)實(shí)驗(yàn)，對(duì)于明確該區(qū)域水動(dòng)力的變化過程有積極的意義，國(guó)內(nèi)在這一方面的工作還存在很大的空白。當(dāng)前圍繞島礁的物理模型實(shí)驗(yàn)多在二維簡(jiǎn)單地形上進(jìn)行，且尚無(wú)有效模擬大粗糙度海床的方式，有待進(jìn)一步開展物理實(shí)驗(yàn)，研究波浪在島礁地形上的破碎變形、波高沿程衰減、波浪能量損耗、低頻波浪運(yùn)動(dòng)、礁坪(瀉湖)上波浪共振、災(zāi)害性波浪爬高、波生流及流場(chǎng)等水動(dòng)力現(xiàn)象的運(yùn)動(dòng)特征和產(chǎn)生機(jī)理；研究礁坪粗糙度對(duì)波浪運(yùn)動(dòng)的影響，進(jìn)而概化出拖曳力系數(shù)和慣性力系數(shù)的估算公式；研究島礁海床特征尺度參數(shù)(如礁前斜坡、礁坪或?yàn)a湖水深)對(duì)上述水動(dòng)力的影響機(jī)制。

數(shù)值計(jì)算模型方面，準(zhǔn)確計(jì)算島礁區(qū)域水動(dòng)力要求模型具備以下能力：優(yōu)良的色散性能和非線性、能適應(yīng)海底地形的劇烈變化、能有效處理波浪破碎和水-陸動(dòng)邊界、能捕捉以對(duì)流為主要形態(tài)的流體運(yùn)動(dòng)以及較經(jīng)濟(jì)的計(jì)算效率等。目前海岸工程常用的相位平均模型和相位識(shí)別模型各自都具有局限和不足，由于相位平均模型計(jì)算效率高，適用于海洋環(huán)境業(yè)務(wù)預(yù)報(bào)，圍繞改善相位平均模型的研究也一直在開展，重點(diǎn)是研究如何在相位平均模型中引入波-波非線性相互作用的機(jī)制。相位識(shí)別模型，如Boussinesq類波浪模型，雖在計(jì)算近岸波浪(短波)運(yùn)動(dòng)過程方面具有優(yōu)勢(shì)，經(jīng)過不斷完善后也出現(xiàn)了許多具有優(yōu)良性能的模型，但通常以增加方程復(fù)雜程度和數(shù)值求解難度為代價(jià)，局限于實(shí)驗(yàn)室尺度，很難普及用于相關(guān)工程計(jì)算。對(duì)于近年興起的非靜壓波浪模型，在色散型、非線性、和海底變化劇烈程度方面沒有限制，具備深海島礁水動(dòng)力計(jì)算的潛質(zhì)，但仍需要系統(tǒng)的理論分析支持以及構(gòu)建準(zhǔn)確、高效的數(shù)值通量算法。發(fā)展具備模擬破碎波浪和水-陸邊界能力的非靜壓類模型，并能適用于三維問題，是該類模型的研究方向。另外，與物理模型實(shí)驗(yàn)類似，如何在計(jì)算模型中考慮島礁海床的大粗糙度仍是當(dāng)前沒有妥善解決的問題。