匡翠萍， 王彬諭， 顧　杰， 雷　剛

(1.同濟(jì)大學(xué) 土木工程學(xué)院，上海 200092； 2.廣西交通職業(yè)技術(shù)學(xué)院，廣西 南寧 530216；3.上海海洋大學(xué) 海洋科學(xué)學(xué)院，上海 201306； 4.國家海洋局第三海洋研究所，福建 廈門 361005)

?

龍鳳頭海灘修復(fù)工程與波流動(dòng)力響應(yīng)特性

匡翠萍1， 王彬諭2， 顧杰3， 雷剛4

(1.同濟(jì)大學(xué) 土木工程學(xué)院，上海 200092； 2.廣西交通職業(yè)技術(shù)學(xué)院，廣西 南寧 530216；3.上海海洋大學(xué) 海洋科學(xué)學(xué)院，上海 201306； 4.國家海洋局第三海洋研究所，福建 廈門 361005)

基于MIKE 21軟件建立平潭龍鳳頭海灘修復(fù)工程海域雙重嵌套二維潮流和波浪耦合數(shù)學(xué)模型，并將驗(yàn)證后的模型運(yùn)用到工程設(shè)計(jì)方案實(shí)施前后潮流場(chǎng)及波浪場(chǎng)的模擬，比較分析水動(dòng)力的變化.研究結(jié)果表明：平潭海域潮汐類型屬于正規(guī)半日潮，8月大潮潮差約為5 m，潮汐強(qiáng)；海壇灣內(nèi)潮流總體表現(xiàn)為往復(fù)流特征，灣口流速最大，灣內(nèi)流速較小，灣頂流速最?。槐倍朔啦ǖ檀蟠鬁p少了南北水體交換，內(nèi)側(cè)落潮流調(diào)整為沿堤流；堤頭處挑流作用明顯，流速增加，形成明顯的旋轉(zhuǎn)流；南北端防波堤對(duì)其內(nèi)側(cè)一定范圍海域有掩護(hù)作用，波高削減顯著，受掩護(hù)岸線長度略大于防波堤長度；堤頭處受波浪折射作用影響，波能輻聚，波高增大.工程后海灘中部岸段波高及流速削減不明顯，建議增設(shè)離岸堤加強(qiáng)海灘保護(hù).

龍鳳頭海灘； 潮流； 波浪； 防波堤； MIKE 21軟件

福建平潭東部海域有大小海灘215處，多處為具有海濱浴場(chǎng)的優(yōu)質(zhì)海灘[1].上世紀(jì)末，海壇灣(圖1)濱海沙灘中部修建的人工海堤，造成海堤前海灘迅速下蝕.2007年海灘基本不存在干灘，大潮時(shí)灘面出露約300 m，不利于游客全天候利用海灘[2].為了將其打造為國內(nèi)一流的高質(zhì)量旅游、休閑娛樂海灘，亟需開展濱海沙灘修復(fù)工程.

海壇灣位于臺(tái)灣海峽西北側(cè)，是一個(gè)半封閉海灣，灣口朝向正東方向，其東部為強(qiáng)潮強(qiáng)浪區(qū).灣內(nèi)島嶼眾多，水深較淺(20 m以內(nèi))，受冬季大風(fēng)和夏季臺(tái)風(fēng)影響，波浪、潮流動(dòng)力強(qiáng)勁，波流相互作用顯著.Fang等[3]指出臺(tái)灣海峽西北部潮汐振幅明顯大于南部海域，相差可達(dá)一個(gè)數(shù)量級(jí)，最大潮差出現(xiàn)在平潭附近海域.Wang等[4]通過2012年7月22至29日ADV(acoustic doppler velocimetry)和ADCP(acoustic doppler current profiler)測(cè)流數(shù)據(jù)得到，臺(tái)灣海峽西側(cè)主要受M2分潮主導(dǎo).Zhang等[5]通過建立臺(tái)灣海峽潮流和風(fēng)暴潮耦合數(shù)學(xué)模型成功模擬了1999年臺(tái)風(fēng)Dan作用下福建沿岸風(fēng)暴潮的潮位變化過程.Zhang等[6]分析了臺(tái)灣海峽5次臺(tái)風(fēng)作用下的水體輸運(yùn)特性，海峽內(nèi)水體呈明顯的向北輸運(yùn)特征.Zeng等[7]基于COHERENS模型模擬得到臺(tái)灣海峽M2分潮的潮汐動(dòng)態(tài)，潮波系統(tǒng)主要由開爾文潮波和駐波兩個(gè)潮波系統(tǒng)共同控制.田永青等[8]基于統(tǒng)計(jì)方法和潮汐調(diào)和分析方法得出臺(tái)灣海峽西部近海最大潮差從平潭島北部到泉州變化不大，廈門東南到惠來呈線性遞減.Zheng等[9]利用WAVEWATCH-Ⅲ模型計(jì)算得到1988～2011年臺(tái)灣島周圍海域波能密度及其增長速率均較大.劉濤等[10]通過建立波浪數(shù)學(xué)模型計(jì)算得到了福建興化灣深水區(qū)域以及工程區(qū)域的設(shè)計(jì)波浪要素.葉雨穎等[11]對(duì)實(shí)測(cè)1個(gè)月的風(fēng)浪資料進(jìn)行了統(tǒng)計(jì)和頻譜分析，給出了福建省南部特征波要素與頻譜關(guān)系.本文在分析實(shí)測(cè)資料的基礎(chǔ)上，利用MIKE 21軟件建立和驗(yàn)證平潭海域雙重嵌套的潮流和波浪耦合數(shù)學(xué)模型，研究分析潮流和波浪對(duì)平潭海壇灣龍鳳頭海灘修復(fù)工程的響應(yīng)特性.

b 設(shè)計(jì)方案圖1　龍鳳頭海灘地理位置及設(shè)計(jì)方案示意圖Fig.1　Location of Longfengtou Beach and design of beach nourishment project

1　數(shù)學(xué)模型及計(jì)算區(qū)域

MIKE 21模型是丹麥水力研究所開發(fā)的平面二維數(shù)學(xué)模型，可用于模擬河口、海岸、海灣、河流、湖泊及海洋的水流、波浪、泥沙及環(huán)境.MIKE 21 FM(flow model)基于Boussinesq假設(shè)、靜水壓力假設(shè)、淺水條件和邊界條件，通過控制體積法求解由不可壓縮雷諾平均Navier-Stokes方程概化的淺水方程[12-13].在空間上采用有限體積法進(jìn)行離散，在時(shí)間上采用顯性歐拉法進(jìn)行離散.MIKE 21 SW(spectral wave)基于波浪作用平衡方程，可以求解準(zhǔn)定常和非定常波浪作用平衡方程的全譜公式和方向解耦參數(shù)公式，并可與潮流模型耦合[14].應(yīng)力模型是波浪模型和潮流模型實(shí)現(xiàn)耦合的關(guān)鍵部分.本文將潮流模型中計(jì)算得到的水位和流速輸入到波浪模型，比較波流耦合作用下工程前后近岸波浪場(chǎng)分布.

1.1計(jì)算區(qū)域及網(wǎng)格

為精準(zhǔn)模擬工程布置和復(fù)雜地形(島嶼眾多，且岸線極不規(guī)則)，采用大、小雙重嵌套模型進(jìn)行計(jì)算.大、小模型布置如圖2所示.大模型大陸沿岸南起廣東石碑山，北至浙江五溪村，臺(tái)灣島沿岸南起淡水，北至高雄.網(wǎng)格剖分采用三角形網(wǎng)格，節(jié)點(diǎn)數(shù)為15 115，網(wǎng)格單元數(shù)為27 713；兩條開邊界網(wǎng)格較為稀疏，空間步長可達(dá)1 000 m；兩條閉邊界網(wǎng)格空間步長約為200 m.地形資料大部分來源于美國國家海洋與大氣管理局(national geophysical data center，NOAA)官方網(wǎng)站，分辨率為2′.小模型計(jì)算范圍覆蓋平潭東部海域，共劃分為三個(gè)開邊界，計(jì)算區(qū)域共布有8 701個(gè)計(jì)算節(jié)點(diǎn)和16 593個(gè)三角形單元，為研究海壇灣內(nèi)的潮流情況，對(duì)其進(jìn)行了網(wǎng)格加密，網(wǎng)格空間步長約20 m.海壇灣內(nèi)地形采用由國家海洋局第三海洋研究所提供的2011年實(shí)測(cè)地形.

1.2工程方案

2011年7月對(duì)海壇灣龍鳳頭海灘實(shí)施灘肩補(bǔ)沙的海灘修復(fù)工程(一期工程)，于2011年8月底竣工，海灘實(shí)際施工岸線總長度1 359 m，其中最北段龜模嶼附近岸段干灘寬度120 m，中北部100 m為過渡段，中南部及南部岸段干灘寬度60 m.工程后沙灘寬度由0擴(kuò)大到20～30 m，至2011年12月，剖面測(cè)量結(jié)果表明海灘侵蝕嚴(yán)重，岸線平均后退11.1 m，最大后退量達(dá)31.0 m，需要對(duì)沙灘進(jìn)行進(jìn)一步修復(fù).設(shè)計(jì)修復(fù)方案為灘肩補(bǔ)沙和防波堤相結(jié)合[15]，布置如圖1所示.防波堤1離岸距離為300 m，長度為260 m，方向?yàn)闁|偏南60°；防波堤2離岸距離為190 m，長度為110 m，方向?yàn)楸逼珫|20°，防波堤1和防波堤2頂寬均為4 m.

a 大模型網(wǎng)格及測(cè)站位置

b 小模型網(wǎng)格及測(cè)站位置圖2　大、小嵌套模型網(wǎng)格及測(cè)站位置Fig.2　Large and small nested model grids and observation stations

2　模型驗(yàn)證

2.1邊界條件和初始條件

大區(qū)域潮流模型的開邊界條件為MIKE 21軟件包自帶的全球潮汐模型導(dǎo)出的潮位過程，該潮汐模型是在1992～2002年共計(jì)10年的TOPEX/POSEIDON衛(wèi)星高度計(jì)資料基礎(chǔ)上建立的，其分辨率為0.125°，模型考慮的分潮包含M2，S2，K2，N2，S1，K1，O1，P1，Q1和M4，共計(jì)10個(gè)分潮.小區(qū)域潮流模型開邊界條件為大模型提供的潮位過程，模擬時(shí)間為2011年8月14日0：00時(shí)至8月17日18：00時(shí).初始潮位選取計(jì)算開始時(shí)刻的平均潮位，初始流速為0.

入射波浪來源于平潭波浪觀測(cè)站(圖2)實(shí)測(cè)波浪的統(tǒng)計(jì)分析數(shù)據(jù).考慮潮流作用時(shí)，將潮流模型中計(jì)算得到的水位和流速輸入到波浪模型，為了保證潮流和波浪耦合，本文選取8月份的實(shí)測(cè)波浪資料進(jìn)行模型驗(yàn)證，8月平均波周期為4.1 s，有效波高和波向的變化過程如圖3所示，統(tǒng)計(jì)得到平均有效波高為1.04 m，波向全部集中在ENE,E,SE和ESE方向，其中以ENE向最多，因此將北邊界和南邊界設(shè)為側(cè)邊界，設(shè)置從東邊界入射的有效波高Hs=1.04 m，波周期T=4.1 s，ENE向的入射波.

圖3　有效波高和波向變化過程Fig.3　Time history of significant wave height and direction

2.2參數(shù)選取

由于潮灘落潮干出、漲潮被淹，因此采用動(dòng)邊界模擬.干水深hdry、淹沒水深hflood和濕水深hwet分別取0.005,0.05和0.1 m.大區(qū)域潮流模型曼寧數(shù)為25 m1/3·s-1，小區(qū)域潮流模型曼寧數(shù)變化范圍為30～50 m1/3·s-1.大小模型計(jì)算時(shí)間步長均采用變步長，根據(jù)流態(tài)自動(dòng)調(diào)整，下、上限值分別為0.000 1 s和30 s.波浪破碎采用Battjes和Janssen模型，破碎指標(biāo)選擇0.8.根據(jù)龍鳳頭海域表層泥沙中值粒徑，粗糙高度取0.1 mm.

龍鳳頭近岸海域風(fēng)向季節(jié)性變化明顯，一年中秋、冬、春三季以偏北風(fēng)為主，夏季以偏南風(fēng)為主.根據(jù)平潭風(fēng)觀測(cè)站(圖2)統(tǒng)計(jì)得到的歷史各月最多風(fēng)向和平均風(fēng)速統(tǒng)計(jì)可知，8月份主要風(fēng)向?yàn)镾SW向，平均風(fēng)速為5.7 m·s-1[2].圖4為對(duì)應(yīng)潮流模擬時(shí)間內(nèi)的風(fēng)向和風(fēng)速變化過程，數(shù)據(jù)來源于NOAA平潭東部海域東經(jīng)120°，北緯25.71°測(cè)點(diǎn)的6 h一次的風(fēng)觀測(cè)數(shù)據(jù)，統(tǒng)計(jì)8月14日0：00時(shí)至17日18：00時(shí)風(fēng)速和風(fēng)向平均值可知與歷史各月最多風(fēng)向和平均風(fēng)速相一致，因此潮流模型和波浪模型中風(fēng)驅(qū)動(dòng)力均采用定常風(fēng)(5.7 m·s-1，SSW向).

2.3潮流模型驗(yàn)證

平潭海域潮汐屬于正規(guī)半日潮，平均漲、落潮歷時(shí)相當(dāng)，8月大潮潮差約為5 m，潮汐強(qiáng).大模型潮位驗(yàn)證資料采用2011年潮汐表上平潭、福清灣、三江口和廈門4個(gè)測(cè)點(diǎn)(圖2)8月16日10：00時(shí)至8月17日12：00時(shí)的潮位過程.平潭測(cè)點(diǎn)的計(jì)算潮位值與潮汐預(yù)報(bào)值的比較如圖5所示.潮位驗(yàn)證結(jié)果表明：臺(tái)灣海峽大模型的潮位計(jì)算值與潮汐預(yù)報(bào)值吻合較好，能較為有效地反映臺(tái)灣海峽的潮位變化過程.

圖4　風(fēng)速和風(fēng)向過程Fig.4　Time history of wind speed and direction

a 大模型(平潭)水位驗(yàn)證

b 小模型(T1)水深驗(yàn)證圖5　大模型(平潭)水位及小模型(T1)水深驗(yàn)證

Fig.5Verification of water level in large scale model and water depth in small scale model

平潭東部海域小模型基于T1和T2測(cè)流站(圖2)一次天文大潮(2011年8月16日10：00時(shí)至8月17日12：00時(shí))的實(shí)測(cè)垂向平均流速、流向及水深數(shù)據(jù)進(jìn)行驗(yàn)證.T1測(cè)流站的水深驗(yàn)證如圖5所示，流速與流向的驗(yàn)證如圖6所示.本文選擇Wilmott[16]提出的統(tǒng)計(jì)學(xué)方法來評(píng)價(jià)數(shù)學(xué)模型，其計(jì)算方法為

(1)

a 流速驗(yàn)證

b 流向驗(yàn)證圖6　小模型(T1)流速和流向驗(yàn)證

Fig.6Verification of flow velocity and direction using small scale model

圖7中矢量表示海壇灣漲、落急時(shí)刻的潮流場(chǎng)，等值線為有效波高分布.潮流總體表現(xiàn)為往復(fù)流特征，灣口流速最大，灣內(nèi)流速較小，灣頂流速最小.漲潮時(shí)水流由灣口北部進(jìn)入，以灣內(nèi)中部蛇嶼為界，分為南北兩支，北支呈順時(shí)針旋轉(zhuǎn)流進(jìn)入海壇灣北部，南支呈逆時(shí)針旋轉(zhuǎn)流從東南角流出；落潮時(shí)灣內(nèi)落潮流和灣外由南向北的落潮流在灣口北側(cè)匯合流出.

2.4波浪模型驗(yàn)證

根據(jù)統(tǒng)計(jì)風(fēng)場(chǎng)資料，對(duì)波流耦合作用下海壇灣內(nèi)波浪場(chǎng)進(jìn)行了模擬，潮流歷時(shí)對(duì)應(yīng)于2011年8月16日10：00時(shí)至8月17日12：00時(shí)的驗(yàn)證潮流過程，其中漲急時(shí)刻為2011年8月16日23：00時(shí)，落急時(shí)刻為2011年8月17日5：00時(shí).模擬結(jié)果與觀測(cè)波浪數(shù)據(jù)統(tǒng)計(jì)結(jié)果均顯示主波向?yàn)镋NE向，觀測(cè)站實(shí)測(cè)波高和波周期分別為1.04 m和4.1 s，計(jì)算波高和波周期分別為1 m和4.2 s，二者一致.圖7中有效波高分布顯示波浪在近岸區(qū)發(fā)生折射，逐漸垂直于岸線方向.由于波浪繞射作用，島嶼后側(cè)波高削減顯著，且落急時(shí)刻波掩區(qū)范圍明顯大于漲急時(shí)刻，波高削減更為明顯.由于潮流對(duì)波浪折射、波浪自身在海壇灣南、北岬頭折射和水深不均勻引起的折射等綜合作用，海壇灣灣口中北部及其內(nèi)側(cè)海域漲急波高明顯弱于落急波高，而南部岬頭及其內(nèi)側(cè)海域漲急波高強(qiáng)于落急波高.

a 漲急時(shí)刻

b 落急時(shí)刻圖7　海壇灣漲、落急時(shí)刻潮流場(chǎng)和波高分布

Fig.7Tidal current fields and distribution of wave height at maximum flood and maximum ebb in Haitan Bay

3　波流動(dòng)力對(duì)工程的響應(yīng)模擬

3.1潮流對(duì)工程的響應(yīng)特性分析

通過對(duì)工程區(qū)及近岸海域的潮流場(chǎng)進(jìn)行模擬，得到漲、落急時(shí)刻工程實(shí)施前后流場(chǎng)分布，分別如圖8和圖9所示.圖10和圖11分別為漲、落急時(shí)刻工程前后流速等值線分布.龍鳳頭海域漲潮流受海壇灣南支漲潮流控制，該由北向南的水流進(jìn)入龍鳳頭海域，受北部龜模嶼及其東側(cè)島嶼和南部龍王頭岬頭影響，形成一個(gè)大的順時(shí)針環(huán)流.工程后龍鳳頭海灘北部及中部岸段近岸流速略有減少，流速等值線向岸推進(jìn).北端防波堤1的阻水作用，使其近岸側(cè)的流速減小，減幅約為0.02 m·s-1，其東側(cè)區(qū)域存在一個(gè)逆時(shí)針環(huán)流.位于南端的防波堤2阻擋了一部分西北向來流，由于堤頭挑流作用，改變了水流結(jié)構(gòu)，內(nèi)側(cè)形成較明顯的逆時(shí)針旋轉(zhuǎn)流，堤頭流速增加至0.12 m·s-1，進(jìn)而可能導(dǎo)致該處沖刷.近岸海域落潮流向主要為東北向.由于龍鳳頭海灘北端存在一個(gè)水深極淺區(qū)域，落潮時(shí)灘地出露，大量落潮流從灘地與龜模嶼之間的深槽流出，流速可達(dá)0.28 m·s-1，工程后防波堤1東側(cè)海域最大流速約為0.18 m·s-1，相較于工程前有所減小.北端防波堤1大大減少了南北水體交換，改變了水流方向，內(nèi)側(cè)落潮流多為沿堤流，堤頭處流速增加至0.14 m·s-1.整個(gè)離岸區(qū)和近岸區(qū)的流速都有所減小，最大減幅小于0.02 m·s-1.

a 工程前

b 工程后圖8　龍鳳頭海域漲急時(shí)刻工程前后流場(chǎng)

Fig.8Maximum flood current fields before and after the project at Longfengtou Beach

a 工程前

b 工程后圖9　龍鳳頭海域落急時(shí)刻工程前后流場(chǎng)Fig.9　Maximum ebb current fields before and after the project at Longfengtou Beach

a 工程前

b 工程后圖10　龍鳳頭海域漲急時(shí)刻工程前后流速等值線分布Fig.10　Isolines of maximum flood current velocity magnitude before and after the project at Longfengtou Beach

3.2波浪對(duì)工程的響應(yīng)特性分析

通過波流耦合作用下波浪模擬，得到漲、落急時(shí)刻工程實(shí)施前后波高分布如圖12和圖13所示.漲急時(shí)刻工程區(qū)及近岸海域波高明顯大于落急時(shí)刻，由于落潮時(shí)灘地出露，龍鳳頭海灘近岸區(qū)域波高極小，約為0.04 m.從圖12可知，漲急時(shí)刻北端防波堤1僅對(duì)其內(nèi)側(cè)海域有掩護(hù)作用，沿岸距北端約350 m范圍內(nèi)波高有所削減，南端防波堤2也僅對(duì)其內(nèi)側(cè)海域有掩護(hù)作用，沿岸距南端約130 m范圍內(nèi)波高有所削減.通過對(duì)工程區(qū)及近岸海域13個(gè)監(jiān)測(cè)點(diǎn)(圖1中J1～J8和N1～N5)工程實(shí)施前后的波高及波向進(jìn)行對(duì)比得到：工程后堤頭位置的測(cè)點(diǎn)波高略有增大，漲落急時(shí)刻J2和J6測(cè)點(diǎn)波高增加在0.1 m以內(nèi)，是由于堤頭處波浪發(fā)生折射，導(dǎo)致波能輻聚，波高增大.由于堤頭繞射作用，波能輻散，防波堤內(nèi)側(cè)為波掩區(qū)，波高下降顯著，因此掩護(hù)區(qū)內(nèi)測(cè)點(diǎn)的波高得到有效削減，J3，J7和J8測(cè)點(diǎn)波高減幅均接近或超過50%，受掩護(hù)岸線長度略大于防波堤長度.中部岸段近岸各測(cè)點(diǎn)波高均只有小幅減小，漲急時(shí)刻N(yùn)2～N4測(cè)點(diǎn)波高減幅均小于1.5%，且越靠近岬灣頂部，測(cè)點(diǎn)波高削減幅度越小，可以看出南北端防波堤對(duì)灣頂近岸海域波浪削減效果不明顯，建議增設(shè)離岸堤加強(qiáng)海灘保護(hù).

a 工程前

b 工程后圖11　龍鳳頭海域落急時(shí)刻工程前后流速等值線分布Fig.11　Isolines of maximum ebb current velocity magnitude before and after the project at Longfengtou Beach

a 工程前

b 工程后圖12　龍鳳頭海域漲急時(shí)刻工程前后波高分布Fig.12　Distribution of wave height at maximum flood before and after the project at Longfengtou Beach

a 工程前

b 工程后圖13　龍鳳頭海域落急時(shí)刻工程前后波高分布Fig.13　Distribution of wave height at maximum ebb before and after the project at Longfengtou Beach

4　結(jié)論

基于MIKE 21軟件建立了平潭龍鳳頭海灘及其鄰近海域的雙重嵌套二維潮流和波浪耦合的水動(dòng)力數(shù)學(xué)模型，通過實(shí)測(cè)資料對(duì)數(shù)學(xué)模型進(jìn)行了驗(yàn)證，并將驗(yàn)證后的數(shù)學(xué)模型運(yùn)用到工程設(shè)計(jì)方案下潮流場(chǎng)和波浪場(chǎng)的預(yù)測(cè)，比較分析波流動(dòng)力對(duì)海灘修復(fù)工程的響應(yīng)關(guān)系，得到以下主要結(jié)論：

(1) 平潭海域潮汐屬于正規(guī)半日潮，平均漲、落潮歷時(shí)相當(dāng)，8月大潮潮差約為5 m，潮汐強(qiáng).海壇灣內(nèi)潮流總體表現(xiàn)為往復(fù)流特征，灣口流速最大，灣內(nèi)流速較小，灣頂流速最小.

(2) 近岸流速受南北防波堤影響較小，離岸區(qū)域流速分布幾乎不受影響；北端防波堤大大減少了南北水體交換，改變了水流方向，內(nèi)側(cè)落潮流多為沿堤流；堤頭處挑流作用明顯，流速增加，形成較明顯的旋轉(zhuǎn)流，進(jìn)而可能導(dǎo)致堤頭沖刷.

(3) 南北側(cè)防波堤只對(duì)其內(nèi)側(cè)一定范圍海域有掩護(hù)作用，受掩護(hù)岸線長度略大于防波堤長度，掩護(hù)區(qū)域波高削減幅度接近50%；堤頭處受波浪折射作用影響，波能輻聚，波高增大.

(4) 工程后中部岸段波高及流速削減不明顯，建議增設(shè)離岸堤加強(qiáng)海灘保護(hù).

[1]林軍, 施文耀. 福建平潭島旅游地質(zhì)資源特征淺析[J]. 福建地質(zhì), 2000, 19(2): 111.

LIN Jun, SHI Wenyao. Discussion on the characteristics of geological tourist resources in Pingtan Island, Fujian Province[J]. Geology of Fujian, 2000, 19(2): 111.

[2]國家海洋局第三海洋研究所. 平潭龍鳳頭海灘養(yǎng)護(hù)實(shí)驗(yàn)工程演變監(jiān)測(cè)報(bào)告[R].北京:國家海洋局第三海洋研究所,2012.

Third Institute of Oceanography, State Oceanic Administration. Report of the bed evolution after the nourishment project in Longfengtou Beach, Pingtan County[R]. Beijing: Third Institute of Oceanography, State Oceanic Administration, 2012.

[3]FANG Guohong, YANG Jingfei, Thao Yucai. A two-dimensional numerical model for tidal motion in the Taiwan Strait[J]. Marine Geophysical Researches, 1984, 7(1/2): 267.

[4]WANG Aijun, YE Xiang, DU Xiaoqin,etal. Observations of cohesive sediment behaviors in the muddy area of the northern Taiwan Strait, China[J]. Continental Shelf Research, 2014, 90: 60.

[5]ZHANG Wenzhou, HONG Huasheng, SHANG Shaoping,etal. A two-way nested coupled tide-surge model for the Taiwan Strait[J]. Continental Shelf Research, 2007, 27:1548.

[6]ZHANG Wenzhou, HONG Huasheng, YAN Xiaohai. Typhoons enhancing northward transport through the Taiwan Strait[J]. Continental Shelf Research, 2013, 56:13.

[7]ZENG Ganning, HU Jianyu, HONG Huasheng,etal. Numerical study on M2tidal system in the Taiwan Strait[J]. Procedia Environmental Science, 2012, 12: 702.

[8]田永青, 潘愛軍. 臺(tái)灣海峽西部近海潮汐特征[J]. 臺(tái)灣海峽, 2011, 30(4): 483.

TIAN Yongqing, PAN Aijun. Tidal characteristics in offshore sea of western Taiwan Strait[J]. Journal of Oceanography in Taiwan Strait, 2011, 30(4): 483.

[9]ZHENG Chongwei, LI Chongyin. Variation of the wave energy and significant wave height in the China Sea and adjacent waters[J]. Renewable and Sustainable Energy Reviews, 2015, 43: 381.

[10]劉濤, 陳漢寶, 劉海源. 福建興化灣波浪條件數(shù)值模擬研究[J]. 水運(yùn)工程, 2007(3): 400.

LIU Tao, CHEN Hanbao, LIU Haiyuan. Numerical simulation research on wave condition of Xinghuawan Bay, Fujian[J]. Port & Waterway Engineering, 2007(3): 400.

[11]葉雨穎, 潘偉然, 張國榮, 等. 福建東山灣海浪現(xiàn)場(chǎng)觀測(cè)的統(tǒng)計(jì)特征[J]. 廈門大學(xué)學(xué)報(bào):自然科學(xué)報(bào), 2007, 46(3): 386.

YE Yuying, PAN Weiran, ZHANG Guorong,etal. The field observations and statistical characteristics of waves in Dongshan Bay[J]. Journal of Xiamen University: Natural Science, 2007, 46(3): 386.

[12]馮靜. MIKE21數(shù)值模型在海洋工程環(huán)境影響評(píng)價(jià)中的應(yīng)用研究[D].青島:中國海洋大學(xué), 2011.

FENG Jing. Application of MIKE21FM numerical model in environmental impact assessment of ocean engineering[D]. Qingdao: Ocean University of China, 2011.

[13]DHI Hydraulics. MIKE21 & MIKE3 flow model hydrodynamic and transport module science documentation[M]. Denmark: DHI Water & Environment, 2009.

[14]DHI Water & Environment & Health. MIKE21 spectral wave module scientific documentation[M]. Denmark: DHI, 2009.

[15]匡翠萍, 單云馳, 王彬諭, 等. 海壇灣龍鳳頭海灘養(yǎng)護(hù)工程方案研究[J]. 同濟(jì)大學(xué)學(xué)報(bào):自然科學(xué)版, 2015, 43(5): 714.

KUANG Cuiping, SHAN Yunchi, WANG Binyu,etal. Study on beach nourishment in Longfengtou Beach of Haitan Bay[J]. Journal of Tongji University: Natural Science, 2015, 43(5): 714.

[16]Wilmott C J. On the validation of models[J]. Physical Geography, 1981, 2(2): 184.

Study of Wave-Current Coupled Hydrodynamic Responses to a Beach Nourishment Project at Longfengtou Beach

KUANG Cuiping1, WANG Binyu2, GU Jie3, LEI Gang4

(1. College of Civil Engineering, Tongji University, Shanghai 200092, China; 2. Guangxi Vocational and Technical College of Communications, Nanning 530216, China; 3. College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China; 4.Third Institute of Oceanography, State Oceanic Administration, Xiamen 361005, China)

A 2D wave-current coupled model for a beach nourishment project sea area at Longfengtou Beach was set up based on MIKE 21 by using a two-level nested grid. The calibrated model was then adopted to simulate the flow fields and wave fields before and after the project, and the hydrodynamic responses to the nourishment project were analyzed. The simulation results show that the tide in Pingtan sea area is regular semidiurnal and the average tidal range at spring tide in August is about 5 m; the tidal current in Haitan Bay is of typical reversing current with a large flow velocity in the mouth of the bay, and a small and the smallest flow velocity in the middle and top of the bay respectively; the northern breakwater, which greatly reduces the north-south water exchange, adjusts the ebb current along the breakwater; the flow velocity at the tip of the breakwaters increases due to the deflecting flow, and an anticlockwise circulation generates obviously; the northern and southern breakwaters only have a limited effect on the area behind the breakwaters where the wave height significantly decreases, and the influenced length of shoreline is slightly greater than that of breakwater; the wave height at the tip of breakwaters increases due to the energy convergence caused by wave refraction. Overall, the reduction of wave height and flow velocity in the middle part of the beach is negligible, thus the construction of offshore breakwater is suggested to strengthen the protection of Longfengtou Beach.

Longfengtou Beach; tidal current; wave; breakwater; MIKE 21

2015-02-07

匡翠萍(1966—), 女, 教授, 工學(xué)博士, 主要研究方向?yàn)楹０豆こ? E-mail: cpkuang@#edu.cn

顧杰(1961—), 男, 教授, 工學(xué)博士, 主要研究方向?yàn)楹Ｑ蟓h(huán)境. E-mail：jgu@shou.edu.cn

TV148+.5

A