劉星池，王永學(xué)，陳靜，2

（1.大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧大連116024；2.南昌工程學(xué)院水利與生態(tài)工程學(xué)院，江西南昌330099）

人工島群分階段建設(shè)對(duì)附近水沙環(huán)境影響的數(shù)值研究

劉星池1，王永學(xué)1，陳靜1，2

（1.大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧大連116024；2.南昌工程學(xué)院水利與生態(tài)工程學(xué)院，江西南昌330099）

以龍口灣正在興建的龍口人工島群和招遠(yuǎn)人工島群為例，利用MIKE21軟件建立龍口灣人工島群附近海域水沙數(shù)學(xué)模型。在模型經(jīng)過實(shí)測(cè)資料驗(yàn)證的基礎(chǔ)上，設(shè)計(jì)了兩大人工島群整體建設(shè)和間隔期不同的分期建設(shè)方案，模擬不同方案對(duì)龍口灣潮流變化和海底沖淤演變的影響。結(jié)果表明，間隔期為9 a的分期建設(shè)方案中人工島群水道內(nèi)流速較大，周圍海域泥沙運(yùn)動(dòng)相對(duì)穩(wěn)定，確定為優(yōu)選方案。

龍口灣；數(shù)值模擬；人工島群；分期填海；水沙環(huán)境變化

隨著我國(guó)沿海地區(qū)經(jīng)濟(jì)快速發(fā)展，對(duì)土地需求日益增大，建設(shè)人工島群是開發(fā)利用海洋空間資源的重要方式。通常大規(guī)模人工島群建設(shè)應(yīng)遵循整體規(guī)劃、分期實(shí)施的原則。在海洋動(dòng)力條件作用下，先期建成的人工島群使周邊的海床產(chǎn)生沖淤變化，且這種變化將作為后期建設(shè)的人工島群的初始海床，而整體建設(shè)不存在這一過程，海床的沖淤變化是島群整體共同影響的結(jié)果。也就是說(shuō)，人工島群整體建設(shè)與分階段建設(shè)產(chǎn)生的海床沖淤演變過程是不同的。因此應(yīng)科學(xué)制定填海先后次序與分期建設(shè)方案，以盡量減少對(duì)海洋環(huán)境的影響。

目前，山東半島的龍口灣臨港高端產(chǎn)業(yè)聚集區(qū)正在興建我國(guó)最大的離岸人工島群，該人工島群由龍口人工島群和招遠(yuǎn)人工島群組成。近年來(lái)，國(guó)內(nèi)在人工島（群）附近海域水動(dòng)力條件及泥沙沖淤演變方面開展了一些研究工作，安永寧等（2010）分析了龍口灣沖淤特性對(duì)龍口人工島群建設(shè)的響應(yīng)；周廣鎮(zhèn)等（2014）計(jì)算并對(duì)比了萊州灣東岸近岸海域人工島群建設(shè)前后的最大沖刷深度；英曉明等（2014）研究了海南儋州由五個(gè)人工島組成的?；◢u對(duì)潮流變化及海床沖淤的影響；趙強(qiáng)等（2014）模擬計(jì)算了圍填海工程對(duì)南黃海輻射沙脊群海域潮流和含沙量變化；許婷等（2015）模擬了波流共同作用下韓江河口附近海域汕頭大型離岸人工島群建設(shè)前后的海床沖淤變化；陳靜等（2016）提出了大連灣內(nèi)多個(gè)圍填海工程使得灣內(nèi)納潮量減小。但上述研究大都是基于離岸人工島群整體建成后對(duì)周圍海域水動(dòng)力環(huán)境的改變和產(chǎn)生的泥沙沖淤變化，關(guān)于動(dòng)態(tài)跟蹤建設(shè)過程對(duì)水沙環(huán)境影響的研究較少。

本文將以龍口灣臨港高端產(chǎn)業(yè)聚集區(qū)人工島群建設(shè)為研究區(qū)域，針對(duì)該區(qū)兩組人工島群建設(shè)整體方案和三種不同分期方案，運(yùn)用MIKE21-FM模型數(shù)值模擬了在龍口人工島群建設(shè)后不同間隔時(shí)間建設(shè)招遠(yuǎn)人工島群對(duì)周圍海域流場(chǎng)以及相應(yīng)的泥沙沖淤變化，并進(jìn)行對(duì)比分析，為人工島群分階段建設(shè)施工方案的優(yōu)化提供參考依據(jù)。

1 龍口灣人工島群概況

龍口灣位于山東半島北部，萊州灣之東北隅，是一個(gè)半開敞海灣。灣內(nèi)水深不足10 m，地形平坦，海底比降為1/1 000，沉積物中值粒徑從岸向海逐漸變細(xì)，粒徑值范圍約為4～8Φ。龍口灣的潮汐性質(zhì)為不正規(guī)的混合半日潮，龍口海洋站平均潮差為91 cm；龍口灣為不規(guī)則半日潮流區(qū)，以旋轉(zhuǎn)流為主，漲潮流速略小于落潮流速，灣口流速大于灣內(nèi)流速。海灣內(nèi)常浪向?yàn)镾W，強(qiáng)浪向?yàn)閃SW，全年平均波高為0.4 m，僅為灣外1/3左右（中國(guó)海灣志編纂委員會(huì)，1991）。

龍口灣臨港高端制造業(yè)聚集區(qū)位于龍口灣南部海域，總體規(guī)劃由龍口人工島群和招遠(yuǎn)人工島群組成（圖1）。其中龍口人工島群包括7個(gè)離岸人工島，填海面積35.237 km2，形成人工岸線57.8 km。人工島群之間預(yù)留SW-NE方向兩條寬550～600 m的主水道,以及NW-SE方向兩條寬200～300 m的輔助水道。招遠(yuǎn)人工島群由4座人工島組成，填海面積12.07 km2，形成人工岸線32.8 km。龍口部分人工島已經(jīng)完成整體圍堰工程，正在進(jìn)行海上填筑；招遠(yuǎn)人工島群也已經(jīng)獲得了國(guó)家有關(guān)部門的審批。

圖1 人工島群布置及小模型模擬區(qū)域

2 潮流場(chǎng)數(shù)值模擬方法簡(jiǎn)介

2.1 模擬方法及參數(shù)設(shè)置

MIKE21是丹麥水動(dòng)力研究所開發(fā)的一種平面二維數(shù)值模型，用于模擬河流、湖泊、河口、海灣、海岸及海洋的水流、波浪、泥沙及環(huán)境。

本文采用MIKE21的水動(dòng)力模塊，利用大小模型嵌套的方法，大模型研究區(qū)域?yàn)槿R州灣（37.08°-38.04°N，118.87°-120.89°E），圖2為大模型模擬區(qū)域及實(shí)測(cè)資料的測(cè)點(diǎn)位置圖。小模型研究區(qū)域（圖1）為龍口灣附近海域（37.48°-37.68°N，120.00°-120.62°E）。數(shù)值離散采用非結(jié)構(gòu)化網(wǎng)格，并對(duì)人工島群建設(shè)區(qū)域進(jìn)行加密處理。由M2、S2、K1、O1、N2、K2、P1、Q1 8個(gè)分潮提供大模型開邊界處邊界條件，分潮調(diào)和常數(shù)由渤海、黃海、東海海洋圖集（海洋圖集編委會(huì)，1993）中各分潮同潮圖插值得到，并參考MIKE21自帶全球潮汐模型提供的參數(shù)。小模型邊界從大模型模擬結(jié)果中提取。泥沙運(yùn)動(dòng)數(shù)值模擬以水動(dòng)力模型為基礎(chǔ)，主要模擬人工島群建設(shè)前后，研究區(qū)在正常天氣下年沖淤狀況。

圖2 大模型模擬區(qū)域（萊州灣）及實(shí)測(cè)資料點(diǎn)位置

2.2 大模型潮位潮流驗(yàn)證

數(shù)學(xué)模型潮位驗(yàn)證依據(jù)2010年5月份的潮位資料（國(guó)家海洋信息中心，2010）；數(shù)學(xué)模型潮流驗(yàn)證依據(jù)國(guó)家海洋局第一海洋研究所在萊州灣2007年10月26日10時(shí)-27日11時(shí)測(cè)得的2個(gè)測(cè)站（L1，L2）的大潮實(shí)測(cè)海流數(shù)據(jù)及2007年11月10日12時(shí)-11日13時(shí)測(cè)得的4個(gè)測(cè)站（H1～H4）大潮實(shí)測(cè)海流數(shù)據(jù)（劉建強(qiáng)，2012）。驗(yàn)證結(jié)果表明，潮位過程線和流速流向過程線形態(tài)基本一致，計(jì)算值與實(shí)測(cè)值吻合較好；潮位驗(yàn)證點(diǎn)最高最低潮位值偏差均不超過±10 cm，潮流驗(yàn)證點(diǎn)平均流向均不超過允許偏差，平均流速除H4和L2略超出允許偏差±10%之外，其他各點(diǎn)均滿足精度要求；其中超出允許偏差的測(cè)站（H4，L2）未超出總測(cè)站數(shù)量的20%，滿足相關(guān)規(guī)程（交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所，2010）要求，說(shuō)明模型建立較為合理因篇幅有限，圖3，圖4，圖5，圖6僅列出萊州港和龍口港大潮期的潮位驗(yàn)證曲線和L1、H3站位的流速流向驗(yàn)證曲線。

圖3 萊州港潮位驗(yàn)證曲線

圖4 龍口港潮位驗(yàn)證曲線

圖5 L1站流速流向驗(yàn)證曲線

圖6 H3站流速流向驗(yàn)證曲線

3 龍口人工島群數(shù)值模擬及分析

3.1 小模型合理性驗(yàn)證

在小模型研究區(qū)內(nèi)的龍口部分人工島群周圍選取15個(gè)流速點(diǎn)（圖7），將模擬計(jì)算后的漲潮時(shí)最大流速及對(duì)應(yīng)流向與參考文獻(xiàn)（安永寧，2010）進(jìn)行對(duì)比，結(jié)果吻合較好（圖8），進(jìn)一步說(shuō)明潮流場(chǎng)數(shù)學(xué)模型的參數(shù)設(shè)置較為合理，可進(jìn)行應(yīng)用計(jì)算。

圖7 小模型流速對(duì)比驗(yàn)證點(diǎn)分布圖

圖8 龍口人工島群周圍各點(diǎn)流速及流向驗(yàn)證

3.2 潮流場(chǎng)數(shù)值模擬結(jié)果

圖9和圖10分別是龍口部分人工島群建設(shè)前后附近海域大潮漲急和落急時(shí)刻的潮流場(chǎng)模擬結(jié)果?？梢钥闯?，大潮漲急時(shí)刻潮流繞過屺姆島后向西南流，受到人工島群影響，一部分經(jīng)西側(cè)防波堤繼續(xù)向南流，一部分進(jìn)入人工島群經(jīng)水道向西南流；大潮落急時(shí)刻與漲潮過程大致相反，一部分潮流繞過防波堤后向北流，一部分經(jīng)人工水道后沿西北方向流出。

圖9 工程建設(shè)前附近海域大潮漲急、落急時(shí)刻流場(chǎng)

圖10 龍口人工島群建設(shè)后附近海域大潮漲急、落急時(shí)刻流場(chǎng)

圖11是大潮漲急、落急時(shí)刻龍口部分人工島群建設(shè)后與建設(shè)前流速值之差的等值線圖。由圖11的比較可見，龍口部分人工島群建設(shè)后附近海域的流速將有不同程度的改變。在人工島群建設(shè)后，人工島群西北側(cè)海域流速整體增加，漲急時(shí)刻比落急時(shí)刻增加幅度大；西南側(cè)海域和北側(cè)海域流速均有不同程度的減小，且落急時(shí)刻比漲急時(shí)刻減小幅度大。人工島群附近最大流速增減值為0.12m/s，西南側(cè)減小范圍為0.02～0.08 m/s，北側(cè)減小范圍為0.02～0.12 m/s，龍口部分人工島群水道內(nèi)流速增加幅度在0.02～0.06 m/s。

圖11 大潮漲急、落急時(shí)刻龍口人工島群建設(shè)前后流速差等值線圖

3.3 泥沙沖淤特性模擬結(jié)果

本文在MIKE21-FM中運(yùn)用MT泥沙輸運(yùn)模塊，考慮潮流作用下小模型研究區(qū)在正常天氣下泥沙年沖淤狀況（圖12）。參考中國(guó)海灣志第三冊(cè)（中國(guó)海灣志編纂委員會(huì)，1991）給出的沉積物中值粒徑分布（圖12），采用竇國(guó)仁（1999）的經(jīng)驗(yàn)公式，計(jì)算出研究區(qū)附近海域的沉積物臨界的啟動(dòng)應(yīng)力。

圖12 萊州灣沉積物中值粒徑等值線圖

圖13龍口人工島群建設(shè)前后附近海域年沖淤量分布圖表明：工程建設(shè)前，附近海域近岸地區(qū)大部分處于微淤積狀態(tài)，年淤積量小于0.024 m，離岸較遠(yuǎn)處部分區(qū)域呈微沖刷態(tài)，年沖刷量也小于0.024 m，整個(gè)小模型研究區(qū)基本上沖淤平衡，沖淤趨勢(shì)與陳明波（2012）利用海圖進(jìn)行的20 a水深變化分析以及安永寧（2010）的泥沙模擬計(jì)算結(jié)果基本一致，沖淤厚度量級(jí)一致，說(shuō)明模型參數(shù)基本合理，可用于泥沙沖淤計(jì)算。龍口部分人工島群建成后，島群西北側(cè)海域處于較強(qiáng)的沖刷狀態(tài)，島群西南側(cè)及北側(cè)海域的淤積范圍和深度較未建設(shè)前明顯增大，但屺姆島西南側(cè)的侵蝕狀態(tài)有所減少。

圖13 龍口人工島群建設(shè)前后附近海域年沖淤量分布圖

4 不同實(shí)施方案數(shù)值模擬及分析

鑒于當(dāng)前龍口人工島群已經(jīng)完成整體圍堰施工、招遠(yuǎn)人工島群處于待建狀態(tài)，本文設(shè)計(jì)了龍口人工島群與招遠(yuǎn)人工島群整體建設(shè)，在龍口人工島群建成后以1 a、5 a和9 a為間隔期再建設(shè)招遠(yuǎn)人工島群，共4種實(shí)施方案（簡(jiǎn)稱方案）列于表1。數(shù)值計(jì)算中忽略了招遠(yuǎn)人工島群建設(shè)對(duì)小模型區(qū)域邊界條件的影響，分別在這4種實(shí)施方案下，模擬了龍口灣海域的水動(dòng)力環(huán)境和年沖淤變化，并計(jì)算10年后海底地形變化量，探討人工島群分階段建設(shè)對(duì)海洋水動(dòng)力及海床沖淤演變的影響。

表1 招遠(yuǎn)人工島群不同間隔期實(shí)施方案

4.1 不同實(shí)施方案對(duì)人工島群附近流場(chǎng)變化的影響

在整體人工島群水道內(nèi)取16個(gè)點(diǎn)（A-P），在整體人工島群周邊地形變化較大處取5個(gè)點(diǎn)（QU），點(diǎn)位分布圖見圖14。根據(jù)所處位置，分別提取4種實(shí)施方案漲潮時(shí)各點(diǎn)最大流速繪制成曲線（圖15：其中A、B、C、I、K用于對(duì)比分析龍口人工島群水道內(nèi)的流速值（圖15（a））；M、N、O、P用于對(duì)比分析招遠(yuǎn)人工島群水道內(nèi)的流速值（圖15（b））；Q、R、S、T、U用于對(duì)比分析整體人工島群臨海方向外圍的流速值（圖15（c））；J、F、G、H、L用于對(duì)比分析整體人工島群臨岸水道內(nèi)的流速值（圖15（d））。

圖14 人工島群周圍計(jì)算點(diǎn)分布

圖15中4個(gè)曲線圖比較顯示不同實(shí)施方案對(duì)招遠(yuǎn)人工島群水道和整體人工島群臨岸水道內(nèi)點(diǎn)最大流速值影響最大；對(duì)龍口人工島群水道內(nèi)點(diǎn)最大流速值影響次之；除S點(diǎn)外，對(duì)整體人工島群臨海外圍各點(diǎn)最大流速值影響最小。整體人工島群建成后，值得關(guān)注的是島間水道中水體的交換能力，水體的交換能力與水體流速成正相關(guān)，且這種交換能力直接影響島周圍的水環(huán)境。從水動(dòng)力環(huán)境角度選擇實(shí)施方案時(shí)，需將人工島群內(nèi)水道流速最大化作為目標(biāo)。由圖15（a）龍口人工島群水道內(nèi)各點(diǎn)漲潮最大流速值顯示方案2最佳，其次是方案3；由圖15（b）招遠(yuǎn)人工島群水道內(nèi)各點(diǎn)漲潮最大流速值顯示方案4最佳，其次是方案3；由圖15（c）整體人工島群臨海外圍各點(diǎn)漲潮最大流速值顯示方案2最佳，其次是方案3；由圖15（d）整體人工島群臨岸水道內(nèi)各點(diǎn)漲潮最大流速值顯示方案4最佳，其次為方案3。

若從全局綜合考慮人工島群可選擇次優(yōu)方案3作為實(shí)施方案；若將易淤積區(qū)的招遠(yuǎn)人工島群水道和整體島群近岸水道作為考慮的重點(diǎn)，則應(yīng)選擇方案4。

4.2 不同實(shí)施方案對(duì)整體島群周圍海底沖淤變化的影響

圖15 不同實(shí)施方案龍口灣內(nèi)人工島群附近各點(diǎn)大潮漲潮時(shí)最大流速值

根據(jù)不同實(shí)施方案累計(jì)10 a的龍口灣海底地形變化分布見圖16。其中，陸域河流界河來(lái)沙量較少在此忽略不計(jì)，計(jì)算過程中未考慮研究區(qū)邊界來(lái)沙情況，故龍口灣內(nèi)泥沙沖刷量與淤積量保持不變。圖16中4個(gè)圖均反映整體島群西北側(cè)海域海底沖刷嚴(yán)重，海底標(biāo)高降低大于0.6 m；臨近北側(cè)海域淤積嚴(yán)重，最大量超過0.6 m；在整體島群西南側(cè)海域海底淤積面積較大，但淤積高度較小，約為0.2 m。由于招遠(yuǎn)人工島群所在位置未建設(shè)龍口人工島群之前是泥沙淤積區(qū)，年淤積量在0.024 m/a左右（見圖13（a）），故整體人工島群的建設(shè)對(duì)該區(qū)西南側(cè)未被占用海域海底泥沙淤積影響不大，但對(duì)西北側(cè)兩個(gè)防波堤頭和北側(cè)海域海底沖淤影響相對(duì)嚴(yán)重。

圖16 不同實(shí)施方案附近海域10年海底地形變化量

對(duì)于4種實(shí)施方案，方案2和方案1兩組人工島群的建設(shè)時(shí)間僅相差一年，先期建設(shè)的龍口人工島群對(duì)海底地形的影響較小，圖16（a）（b）顯示方案1和方案2十年后海底地形變化量相差不大；方案3和方案4兩個(gè)人工島群建設(shè)的間隔期較長(zhǎng)，圖16（c）（d）顯示龍口人工島群西側(cè)防波堤頭附近與招遠(yuǎn)人工島群西側(cè)防波堤頭附近的地形變化量與方案1和方案2明顯不同。在龍口人工島群圍堰完成后未建招遠(yuǎn)人工島群的這段時(shí)間間隔期內(nèi)，龍口島群西側(cè)防波堤頭處附近海域處于明顯沖刷狀態(tài)。因此隨著間隔期的變長(zhǎng)，龍口人工島群西側(cè)防波堤頭附近的地形沖刷深度與范圍變大，表現(xiàn)在圖16（d）方案4的沖刷深度與范圍要大于圖16（c）方案3情形。招遠(yuǎn)人工島群建設(shè)后改變了局部的流場(chǎng)，減弱了龍口人工島群西側(cè)防波堤頭附近的地形沖刷深度與范圍，招遠(yuǎn)人工島群西側(cè)防波堤頭附近進(jìn)入沖刷狀態(tài)。

圖17顯示不同實(shí)施方案整體島群建成后龍口島群水道內(nèi)、招遠(yuǎn)島群水道內(nèi)和整體周圍海域各點(diǎn)的海底年沖淤量。其中正值表示淤積量，負(fù)值表示沖刷量，研究點(diǎn)沖淤量值離X軸越近沖淤量越小，海床越穩(wěn)定。（a）（b）（c）（d）4個(gè)圖聯(lián)合起來(lái)綜合考慮，發(fā)現(xiàn)方案4為優(yōu)選方案，其次為方案3，方案1和方案2基本相當(dāng)。但是，由于龍口人工島群東北側(cè)臨岸水道沒有南北貫通留有拐角并且水道較窄，附近J、F點(diǎn)淤積較嚴(yán)重，方案2-4的年淤積量均比方案1大。

圖18顯示不同實(shí)施方案整體人工島群內(nèi)及周圍各點(diǎn)10年后地形變化情況，其中正值表示淤積，負(fù)值表示沖刷。圖18（a）中兩組人工島群建設(shè)間隔時(shí)間越長(zhǎng)龍口人工島群水道內(nèi)沖刷越嚴(yán)重，淤積量越少；圖18（b）中招遠(yuǎn)人工島群水道內(nèi)以淤積為主，兩組人工島群建設(shè)間隔時(shí)間越長(zhǎng)招遠(yuǎn)人工島群水道內(nèi)淤積量越少；在圖18（c）中對(duì)于整體人工島群臨海外圍各點(diǎn)兩組人工島群建設(shè)間隔時(shí)間越短S點(diǎn)沖刷量越小，T點(diǎn)特征與S點(diǎn)相反，Q、R、U點(diǎn)變化不大；圖18（d）整體人工島群臨岸水道中J點(diǎn)淤積量最大，其他各點(diǎn)淤積量較小，與圖17（d）結(jié)果一致，但整體人工島群臨岸水道內(nèi)10年淤積量方案4最小。若從減少人工島群水道內(nèi)泥沙淤積的角度選擇實(shí)施方案時(shí)，可選擇方案4。

圖17 不同實(shí)施方案整體島群建成后附近各點(diǎn)年沖淤量

圖18 不同實(shí)施方案整體人工島群附近各點(diǎn)10年海底地形變化量圖

5 結(jié)論

本文以龍口灣臨港高端制造業(yè)聚集區(qū)的龍口與招遠(yuǎn)兩大人工島群工程為例，運(yùn)用MIKE21軟件計(jì)算并分析了在龍口人工島群圍堰建成后，招遠(yuǎn)部分人工島群相隔不同年份建設(shè)對(duì)工程水域水動(dòng)力與泥沙的影響，可得到如下結(jié)論：

（1）人工島群分期建設(shè)中，先期建設(shè)的人工島群會(huì)引起附近海域的海底地形和海洋動(dòng)力環(huán)境條件較大的改變，因而后期人工島群建成后的附近海域地形沖淤演變不同于兩部分人工島群整體建設(shè)方案。實(shí)際工程中可根據(jù)建設(shè)需求的不同選擇最優(yōu)分階段建設(shè)方案。

（2）在龍口部分人工島群和招遠(yuǎn)部分人工島群分階段建設(shè)過程中，綜合考慮整體島群水道內(nèi)和臨海外圍海域的流場(chǎng)及泥沙沖淤變化情況，分期間隔較長(zhǎng)的方案為優(yōu)選方案。

（3）本文僅考慮潮流場(chǎng)變化下的招遠(yuǎn)部分人工島群相隔不同年份填海建設(shè)對(duì)工程區(qū)域附近海底地形變化的影響。鑒于人工島群工程離岸較近、水深較淺，波浪與水流共同作用下人工島群分階段建設(shè)對(duì)海床沖淤演變的影響需進(jìn)一步的探討。

安永寧，吳建政，朱龍海，等，2010.龍口灣沖淤特性對(duì)人工島群建設(shè)的響應(yīng)，海洋地質(zhì)動(dòng)態(tài)，26（10）：24-30.

周廣鎮(zhèn)，馮秀麗，劉杰，等，2014.萊州灣東岸近岸海域規(guī)劃圍填海后沖淤演變預(yù)測(cè)，海洋科學(xué)，38（1）：15-19.

英曉明，賈后磊，曹玲瓏，等，2014.海南儋州人工島建設(shè)對(duì)海床沖淤影響數(shù)值模擬，海洋通報(bào)，33（5）：571-575.

趙強(qiáng)，何琴燕，楊耀芳，等，2014.人工島工程對(duì)南黃海輻射沙脊群海域潮流泥沙影響研究，海洋通報(bào)，33（4）：397-404.

許婷，劉國(guó)亭，溫春鵬，2015.波流共同作用下大型離岸人工島群工程海床沖淤預(yù)測(cè)，海岸工程，34（4）：11-22.

陳靜，王永學(xué)，2016.岸線變遷對(duì)大連灣海域納潮量的影響，海洋通報(bào)，待發(fā)表.

中國(guó)海灣志編纂委員會(huì)，1991.中國(guó)海灣志第三分冊(cè)（山東半島北部和東部海灣），北京：海洋出版社.

海洋圖集編委會(huì)，1993.渤海、黃海、東海海洋圖集（水文），北京：海洋出版社.

國(guó)家海洋信息中心，2010.潮汐表第一冊(cè)：鴨綠江口至長(zhǎng)江口.北京：海洋出版社.

劉建強(qiáng)，2012.萊州灣海洋工程建設(shè)對(duì)小清河口環(huán)境影響數(shù)值研究.碩士學(xué)位論文.青島：中國(guó)海洋大學(xué).

交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所，2010.JTS/T231-2-2010，海岸與河口潮流泥沙模擬技術(shù)規(guī)程.北京：人民交通出版社.

安永寧，2010.離岸人工島群建設(shè)對(duì)龍口灣沖淤特征的影響研究.碩士學(xué)位論文.青島：中國(guó)海洋大學(xué).

竇國(guó)仁，1999.再論泥沙起動(dòng)流速，泥沙研究，（6）：1-9.

陳明波，2012.萊州淺灘對(duì)萊州灣東部沉積動(dòng)力格局的控制作用研究.碩士學(xué)位論文.青島：中國(guó)海洋大學(xué).

（本文編輯：袁澤軼）

Study on the water-sediment environment of artificial islands constructed in stages by numerical simulation

LIU Xing-chi1,WANG Yong-xue1,CHEN Jing1,2

（1.State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology,Dalian 116024,China; 2.Department of Hydraulic and Eco-engineering,Nanchang Institute of Technology,Nanchang 330099,China）

Hydrodynamic and sediment model for the Longkou bay nearby waters is set up by MIKE21 software in this paper.The overall construction scheme and Installment construction schemes of different intervals are simulated on the basis of the validation by the measured data。The results show that the velocity within the waterways is relatively larger and the sediment movement around the artificial islands is relatively stable in the 9-year interval scheme,so this scheme is preferred.

Longkou Bay;numerical simulation;artificial islands;staging reclamation;water and sediment environment changes

P753

A

1001－6932（2017）03－0302－09

10.11840/j.issn.1001-6392.2017.03.008

2016-01-24；

2016-03-18

國(guó)家自然科學(xué)創(chuàng)新研究群體基金（51221961）。

劉星池（1991－），碩士研究生，主要從事海岸工程與海洋環(huán)境科學(xué)方面研究。電子郵箱：lxc102@mail.dlut.edu.cn。