龔 煌，孫大鵬，孫志國

（大連理工大學(xué)海岸和近海工程國家重點(diǎn)實(shí)驗(yàn)室，大連116023）

人口的增長(zhǎng)和社會(huì)的發(fā)展，有限的陸地資源難以滿足人們的用地需求，而廣闊的海洋面積為我們提供了巨大的發(fā)展空間，為此，人工島的建設(shè)成為人們利用海洋資源解決土地需求矛盾的一種有效途徑。早在20 世紀(jì)60 年代，國外就出現(xiàn)了卵石和沙壩式人工島。到20 世紀(jì)80 年代，美國、日本等先后完成了一批鋼筋混凝土人工島和鋼結(jié)構(gòu)人工島的建設(shè)［1］。國內(nèi)外人工島標(biāo)志性工程的完成，如關(guān)西國際機(jī)場(chǎng)、澳門國際機(jī)場(chǎng)以及被稱為“阿拉伯之星”的迪拜帆船酒店等［2-4］，使得人工島的建設(shè)技術(shù)日趨成熟。關(guān)于人工島水動(dòng)力方面，國內(nèi)李孟國、時(shí)鐘［5］等對(duì)如東海域西太陽沙的人工島工程進(jìn)行了潮流模擬分析；秦崇仁［6］、季榮耀［7］等人分別探討了人工島局部沖刷問題和人工島對(duì)水沙動(dòng)力環(huán)境的影響；郭磊［8］等分析了人工島對(duì)周邊海域水動(dòng)力的影響；方偉［9］等人初步分析了灘海工程中的環(huán)境因子與人工島的相互作用；陸敏［10］利用MIKE21 flow model模塊對(duì)嵌入人工島的海灣水體交換進(jìn)行了研究。應(yīng)該說，在人工島對(duì)海灣水環(huán)境影響及評(píng)價(jià)方面的研究還不夠成熟，由于人工島的嵌入，海灣內(nèi)環(huán)流的形成會(huì)極大程度改變?cè)撍虻乃畡?dòng)力特征，污染物擴(kuò)散變得更為復(fù)雜，為避免人工島建設(shè)引發(fā)惡化的環(huán)境代價(jià)，需在前期可行性研究階段對(duì)污染物運(yùn)動(dòng)進(jìn)行模擬和預(yù)測(cè)。

本文基于ECOMSED 三維數(shù)值模型［11］，以大連灣人工島建設(shè)為研究背景，建立全灣的三維潮流模型。采用拉格朗日粒子追蹤技術(shù)，對(duì)建設(shè)人工島前后的大連灣污染物運(yùn)動(dòng)進(jìn)行模擬，并以人工島對(duì)海灣污染物擴(kuò)散影響最小甚至改善其擴(kuò)散狀態(tài)為控制指標(biāo)，從污染物擴(kuò)散的角度推薦可供規(guī)劃參考的人工島建設(shè)方案和污水排放口的布置方式，為人工島的規(guī)劃設(shè)計(jì)提供參考意見。

1 ECOMSED 模型簡(jiǎn)介

ECOMSED 模式是一個(gè)模擬水動(dòng)力、波浪和沉積物輸運(yùn)的三維數(shù)值模型，是Blumberg 等人在海洋數(shù)值模型(Princeton Ocean Model,POM)及其后來的河口、陸架和海洋模式(Estuarine Coastal and Ocean Model,ECOM)基礎(chǔ)上發(fā)展而來。該模型從N-S 方程出發(fā)，采用靜力學(xué)假定和Boussinesq 近似假定下的封閉海洋方程組。模型水平方向采用正交曲線坐標(biāo)，垂向上通過采用σ 坐標(biāo)可以較直角坐標(biāo)更準(zhǔn)確方便地處理淺海的海底變形，并在垂向上通過嵌套2.5 階湍流閉合模型獲得垂向湍流粘滯和擴(kuò)散系數(shù)。為了提高計(jì)算效率，ECOMSED采用內(nèi)外模態(tài)分離技術(shù)，外部模式(垂向積分方程)采用小的時(shí)間步長(zhǎng)，內(nèi)部模式(三維方程)采用大的時(shí)間步長(zhǎng)。

示蹤粒子擴(kuò)散方程在σ 坐標(biāo)下表達(dá)如下

式中：η 為海表面到平均水平面距離；D 為自由表面與底床之間的距離；ξ1、ξ2、σ 分別為水平正交曲線坐標(biāo)和垂向坐標(biāo)；（U1、U2、ω）為（ξ1、ξ2、σ）層上的速度分量；h1、h2為拉梅系數(shù)。KH、AH為垂向和水平方向上的湍流擴(kuò)散系數(shù)。KH由Mellor and Yamada［12］2.5 階湍流閉合模式求解，AH可用Smagornisky［13］參數(shù)化方法求解

其中α 取0.1～0.2 之間的常數(shù)。

圖1 大連灣地形及測(cè)站位置Fig.1 Landform of Dalian bay and station locations

2 數(shù)值模擬

2.1 人工島方案

大連灣是一個(gè)天然的半封閉海灣，位于遼東半島南端，其三面為陸地包圍僅東南面與黃海相通，大連灣地形如圖1 所示，灣內(nèi)自左向右有臭水套、甜水套和紅土堆子灣等3 個(gè)副灣為污染物排放口的規(guī)劃預(yù)設(shè)位置。

人工島方案的確定主要依據(jù)長(zhǎng)度和寬度2 個(gè)指標(biāo)，在大連灣口門外連一條直線，定義直線方向?yàn)槿斯u寬度方向，垂直于該直線方向?yàn)槿斯u長(zhǎng)度方向。沿灣頂陸地突出部分描一條曲線，作為灣頂（圖2）。以L表征灣口到灣頂曲線的距離，以B 表征灣內(nèi)兩岸平行于灣口直線的距離。人工島兩側(cè)岸線作了取直處理，而頂岸線連接側(cè)邊岸線則盡量平滑過渡。本研究以表1所示的3 種人工島方案為比選基礎(chǔ)，對(duì)比分析3 個(gè)副灣為污染物排放口的前提下，人工島建設(shè)對(duì)海灣污染物運(yùn)動(dòng)的影響。

圖2 人工島預(yù)設(shè)方案布置Fig.2 Arranged schemes of artificial island

2.2 計(jì)算范圍

模型計(jì)算范圍為33 km×44 km（圖1），水平方向上采用正交矩形網(wǎng)格，網(wǎng)格長(zhǎng)度為ΔX=250 m，ΔY=270 m；垂向采用σ 分層，等間距分為10 層。模擬時(shí)間從2003 年12 月1 日0：00 至2003 年12 月15 日24：00，共計(jì)算16 d。該模型采用均一的溫鹽度，分別為10℃和31.0 g/L。

為了在三維條件下準(zhǔn)確模擬污染物運(yùn)動(dòng)，先建立整個(gè)渤海和部分黃海海域的大尺度模型，以威海港和小長(zhǎng)山島實(shí)測(cè)潮位資料線性插值給定邊界，驗(yàn)證后提取計(jì)算區(qū)域的邊界資料，驗(yàn)證潮位站資料取自潮汐表。

表1 大連灣人工島的比選方案Tab.1 Scheme of artificial island

2.3 模型驗(yàn)證

模型驗(yàn)證包括潮位站的潮位驗(yàn)證和流速站點(diǎn)0.2、0.6、1.0 倍水深處的速度驗(yàn)證。潮位驗(yàn)證資料取自老虎灘潮位站，流速資料取自D1～D4 站的實(shí)測(cè)值，位置見圖1。這里以老虎灘潮位及D4 大潮速度驗(yàn)證為例（圖3）。

圖3 潮位和流速驗(yàn)證圖Fig.3 Verification of tidal level and flow velocity

驗(yàn)證結(jié)果中，潮位站的潮位過程吻合較好；流場(chǎng)驗(yàn)證中除0.2 h 處流速出現(xiàn)峰值時(shí)計(jì)算值稍偏大外，其余流速和流向均驗(yàn)證良好，這可能跟當(dāng)時(shí)天氣有關(guān)，而該模型忽略了氣候等因素。

3 示蹤粒子運(yùn)動(dòng)模擬

為了研究人工島的嵌入對(duì)污染物運(yùn)動(dòng)規(guī)律的影響，這里僅以人工島方案Ⅰ、甜水套灣內(nèi)設(shè)置污染物粒子點(diǎn)源為例介紹示蹤粒子運(yùn)動(dòng)模擬。粒子釋放時(shí)間為2003 年12 月2 日零時(shí)，共模擬15 d。

3.1 單個(gè)粒子運(yùn)動(dòng)軌跡

為對(duì)比分析嵌入人工島前后單個(gè)粒子的三維運(yùn)動(dòng)規(guī)律，在甜水套灣內(nèi)粒子源的水體表層、中層和底層各放置1 個(gè)粒子進(jìn)行模擬追蹤，粒子位置記錄間隔為1 h，來分析人工島對(duì)海灣污染物運(yùn)動(dòng)的影響。

如圖4 所示，從15 d 的模擬過程來看，無人工島時(shí)，甜水套表層釋放的粒子在灣內(nèi)呈現(xiàn)出沿甜水套海岸線的緩慢運(yùn)動(dòng)，15 d 后仍滯留在副灣內(nèi)；中、底層粒子首先呈現(xiàn)平行岸線的進(jìn)出灣往復(fù)運(yùn)動(dòng)，同時(shí)伴有向東移動(dòng)，繞過岸線運(yùn)動(dòng)至紅土堆子灣灣口水域經(jīng)過短暫徘徊后停止。

人工島嵌入之后，灣內(nèi)水體形成環(huán)流，甜水套在漲落潮流下都成為頂沖點(diǎn)。各層粒子運(yùn)動(dòng)均與無人工島時(shí)明顯不同，在漲落潮流的作用下沿人工島頂端岸線做往復(fù)運(yùn)動(dòng)。經(jīng)過15 d 的模擬，中層粒子經(jīng)過較長(zhǎng)時(shí)間的往復(fù)運(yùn)動(dòng)后，從大連灣左側(cè)運(yùn)動(dòng)至口門外；而表層和底層粒子最終仍滯留在灣內(nèi)。

3.2 粒子群運(yùn)動(dòng)軌跡

鑒于單個(gè)粒子運(yùn)動(dòng)具有較強(qiáng)的隨機(jī)性，故通過模擬污染物粒子群運(yùn)動(dòng)，以更準(zhǔn)確地研究污染物的運(yùn)動(dòng)規(guī)律。在甜水套粒子源表層放置1 000 個(gè)相同粒子，每隔1 h 記錄一次粒子坐標(biāo)位置，每隔5 d 對(duì)粒子群位置輸出一次（圖5）。

由圖5-a 可以看出，在污染物運(yùn)動(dòng)5 d 后，粒子群的分布狀況明顯不同。在原始地形條件下，粒子群主要在甜水套副灣內(nèi)運(yùn)動(dòng)，很少有粒子運(yùn)動(dòng)至副灣外；嵌入人工島之后，除少部分粒子滯留在灣內(nèi)，其余粒子則沿人工島頂端岸線呈對(duì)稱狀分布在人工島兩側(cè)。可見在人工島的挑流作用下，該副灣處漲落潮流速度增大，使得該處污染物擴(kuò)散明顯加快。

在模擬10 d 后，污染物均整體呈現(xiàn)繼續(xù)向外擴(kuò)散的趨勢(shì)。無人工島時(shí)，粒子群在沿岸流的作用下整體呈現(xiàn)紅腿堆子灣方向的緩慢移動(dòng)，其中小部分粒子已運(yùn)動(dòng)至紅土堆子灣；建造人工島之后，一部分粒子繼續(xù)運(yùn)動(dòng)至紅土堆子灣附近，一部分則沿大連灣左側(cè)岸線向?yàn)惩膺\(yùn)動(dòng)，且少量粒子已運(yùn)動(dòng)至大連灣口門外，如圖5-b。

粒子群經(jīng)過15 d 的模擬后，位置布置如圖5-c。無人工島時(shí)，經(jīng)過15 d 的模擬基本沒有粒子運(yùn)動(dòng)至口門外，粒子群整體表現(xiàn)出向紅土堆子灣方向緩慢移動(dòng)的態(tài)勢(shì)，可期出灣路徑為大連灣右側(cè)；有人工島時(shí)，除一部分粒子從大連灣左側(cè)運(yùn)動(dòng)出灣外，較多粒子運(yùn)動(dòng)至大連灣右側(cè)口門靠灣內(nèi)不遠(yuǎn)處，可以預(yù)見其將從右側(cè)出灣。

人工島的嵌入明顯改變污染物運(yùn)動(dòng)軌跡以及污染物的出灣路徑，甜水套的污染物從原來的沿大連灣右側(cè)運(yùn)動(dòng)出灣，變?yōu)閺拇筮B灣右側(cè)和左側(cè)2 條出灣路徑。其次，大連灣在人工島的挑流下形成的環(huán)流，使得灣內(nèi)流速增大，副灣處的污染物擴(kuò)散速度也明顯增加?？傊?，建造人工島后，灣內(nèi)的潮流環(huán)境產(chǎn)生很大差異，污染物的運(yùn)動(dòng)也隨著發(fā)生明顯變化。

圖5 甜水套粒子群位置圖Fig. 5 Location of particles in Tianshuitao

4 人工島方案及排污方式推薦

4.1 人工島方案推薦

海灣的水體交換評(píng)價(jià)方法很多，包括時(shí)間尺度［14］、物質(zhì)交換率［15］、關(guān)聯(lián)矩陣法［16］等。人工島的建造使得海灣內(nèi)的水域面積減少，為使量化評(píng)估的指標(biāo)更加準(zhǔn)確，本研究以前面所介紹的示蹤粒子運(yùn)動(dòng)追蹤為基礎(chǔ)，通過各灣的粒子平均濃度來評(píng)估人工島對(duì)海灣水環(huán)境的影響。在大連灣口門及3 個(gè)副灣灣口處沿水深各取定一個(gè)邊界，設(shè)各邊界與各灣所圍水域的水體體積為V，該區(qū)域水體內(nèi)存留的粒子數(shù)為N，定義平均粒子濃度為C=N/V。

針對(duì)3 種不同人工島方案，通過在3 個(gè)副灣表層、中層、底層污染物粒子排放口各同時(shí)釋放1 000 粒子，對(duì)15 d 的粒子群運(yùn)動(dòng)進(jìn)行追蹤模擬，統(tǒng)計(jì)15 d 后各副灣及大連灣內(nèi)的存留粒子數(shù)并計(jì)算各灣水域的平均粒子濃度。這里以3 個(gè)副灣處的表層粒子群為例，分析各方案下各灣平均粒子濃度的變化。對(duì)于3 個(gè)副灣，副灣粒子群在各副灣內(nèi)的平均粒子濃度如圖6～圖8；對(duì)于全灣水域，15 d 后各粒子群在大連灣水域的平均粒子濃度如圖9。

由圖6 可以看出，臭水套釋放的粒子群經(jīng)15 d 的擴(kuò)散后，無人工島只有在紅土堆子灣內(nèi)的平均粒子濃度比3 種方案條件下小，其余較后者要大。這是由于臭水套污染物經(jīng)過15 d 的擴(kuò)散尚未運(yùn)動(dòng)至紅土堆子灣，在之前的粒子群運(yùn)動(dòng)分析中也討論過。由圖可知方案Ⅲ條件下，臭水套污染物擴(kuò)散最快，方案Ⅱ次之，方案Ⅰ最慢。如圖7 可知，方案Ⅲ可以明顯改善各個(gè)副灣內(nèi)的污染物擴(kuò)散，起到利好的作用；方案Ⅱ在臭水套內(nèi)的平均粒子濃度比無人工島時(shí)稍大，但對(duì)其他2 個(gè)副灣也能起到改善作用；方案Ⅲ則只有在甜水套內(nèi)平均粒子濃度比無人工島小，其余均較其要大。無人工島時(shí)，紅土堆子灣釋放的粒子群只有在紅土堆子灣內(nèi)的平均粒子濃度比方案Ⅲ要大，其余情況下均為最小（圖8）。

而對(duì)于整個(gè)大連灣水域，各個(gè)副灣粒子群經(jīng)過15 d 的模擬后，在全灣水域的平均粒子濃度如圖9。由圖可知，在3 種方案條件下的全灣水域平均粒子濃度均比無人工島時(shí)大，且從方案Ⅲ到方案Ⅰ，其污染物擴(kuò)散速度依次減少?？梢?，人工島的嵌入導(dǎo)致大連灣水域面積減少，對(duì)全灣水域的污染物擴(kuò)散產(chǎn)生不利影響。

表2 不同人工島方案下各粒子源15 d 后在各灣內(nèi)的污染物平均濃度Tab.2 Concentration of the bay with particles in different gulfs after 15 d 個(gè)/萬m3

3 個(gè)副灣處污染物粒子源表層、中層、底層粒子群經(jīng)15 d 運(yùn)動(dòng)擴(kuò)散后，在各副灣及全灣水域的平均粒子濃度如表2 所示。以人工島對(duì)海灣污染物擴(kuò)散影響最小甚至改善其擴(kuò)散狀態(tài)為目標(biāo)，以粒子平均濃度為評(píng)估參數(shù)，經(jīng)綜合分析、比選3 種人工島建設(shè)方案，方案Ⅲ相對(duì)較優(yōu)，宜為人工島建設(shè)的推薦方案。

圖10 臭水套釋放粒子群在各灣內(nèi)的平均粒子濃度Fig.10 Concentration of the bay with particles in Choushuitao

4.2 污染物排放方式推薦

以上推薦的人工島方案Ⅲ為基礎(chǔ)，在3 個(gè)副灣粒子源表層、中層、底層分別預(yù)設(shè)3 個(gè)污染物排放口，每個(gè)排放口放置1 000 個(gè)粒子，通過模擬污染物粒子的運(yùn)動(dòng)，統(tǒng)計(jì)粒子群經(jīng)過15 d 擴(kuò)散后在各副灣及大連灣內(nèi)的平均粒子濃度如圖10～圖12。

圖11 甜水套釋放粒子群在各灣內(nèi)的平均粒子濃度Fig.11 Concentration of the bay with particles in Tianshuitao

圖12 紅土堆子灣釋放粒子群在各灣內(nèi)的平均粒子濃度Fig.12 Concentration of the bay with particles in Hongtuduizi bay

對(duì)于臭水套各層粒子群，在各灣內(nèi)的平均粒子濃度如圖10，可見中層粒子擴(kuò)散速度最快、表層次之、底層最慢，故宜將該副灣處污水排放口布置在水體中層；對(duì)于甜水套（圖11），中層粒子群甜水套灣內(nèi)平均粒子濃度最少，且較快向紅土堆子灣運(yùn)動(dòng)，污染排放口適宜布置在水體中層；而在紅土堆子灣水域，通過其各層粒子群在各灣內(nèi)的粒子濃度（圖12），可知其底層粒子群在該副灣平均粒子濃度最低，宜將污染排放口優(yōu)先布置在水體底層。

5 結(jié) 語

本文借助ECOMSED 建立大連灣三維潮流數(shù)值模型，采用拉格朗日粒子追蹤技術(shù)，對(duì)嵌入人工島前后的海灣污染物運(yùn)動(dòng)進(jìn)行了模擬。針對(duì)3 種人工島比選方案、3 個(gè)污染物預(yù)設(shè)排放口，對(duì)比分析了污染物運(yùn)動(dòng)的變化，以污染物平均濃度為評(píng)估參數(shù)，推薦了可供規(guī)劃參考的人工島建設(shè)方案并對(duì)污染物的排放布置給出建議。本研究從污染物擴(kuò)散的角度，為人工島的規(guī)劃設(shè)計(jì)提供技術(shù)支持，進(jìn)一步完善了人工島對(duì)海灣水環(huán)境影響的評(píng)價(jià)內(nèi)容，可以預(yù)估海灣因建設(shè)人工島對(duì)水環(huán)境的影響，具有一定的工程應(yīng)用價(jià)值。

［1］GERRISTSEN H,SUTHERLAND J,DEIGAARD R,et al. Composite modeling of the interactions between beaches and structure［J］.Journal of Hydraulic Research,2011,47(S1):2-14.

［2］邱毅平.澳門國際機(jī)場(chǎng)海上人工島工程設(shè)計(jì)簡(jiǎn)介［J］.上海水利,1993(3):28-33.QU Y P.The introduction to engineering design of Macau international airport artificial island［J］.Shanghai Water,1993(3):28-33.

［3］謝浩.迪拜奇跡之帆船酒店［J］.綠色建筑,2010(2):3-4.XIE H. Dubai burjAl-Arab［J］.Green building, 2010(2):3-4.

［4］中瀨,明男,呂莉娃.關(guān)西國際機(jī)場(chǎng)人工島建設(shè)［J］.海岸工程,1994(2):197-201.ZHONG L,MING N,LV L W. The construction of Kansai international airport artificial island［J］.Coastal Engineering,1994(2):197-201.

［5］李孟國,時(shí)鐘.江蘇如東海域西太陽沙人工島工程潮流數(shù)學(xué)模型研究［J］.中國海港建設(shè),2006(3):1-4.LI M G, SHI Z.Numerical Modeling Study of Tidal Flow for Construction of An Artificial Island on Xitaiyangsha Shoal at Rudong,Jiangsu Province［J］.China Harbour Engineering,2006(3):1-4.

［6］秦崇仁,肖波,高學(xué)平.波浪水流共同作用下人工島周圍局部沖刷的研究［J］.海洋學(xué)報(bào),1994(3):130-138.QIN C R,XIAO B,GAO X P. The study of scour around the island under the wave flow interaction［J］.Acta Oceanologica Sinica,1994(3):130-138.

［7］季榮耀,徐群,莫思平,等. 港珠澳大橋人工島對(duì)水沙動(dòng)力環(huán)境的影響［J］.水科學(xué)進(jìn)展,2012(6):842-850.JI R Y,XU Q,MO S P, et al. Effects on the hydrodynamics and sediment environment by artificial islands of the Hong Kong-Zhuhai-Macao Bridge［J］.Advances in Water Science,2012(6):842-850.

［8］郭磊,楊樹森.江蘇如東人工島工程對(duì)周邊海域水動(dòng)力的影響［J］.水道港口,2009(5):342-346.GUO L, YANG S S.Influence of artificial island construction on hydrodynamics in waters around Rudong［J］.Journal of Waterway and Harbor,2009(5):342-346.

［9］方偉,李凱雙,張桂榮.灘海工程環(huán)境因子與人工島工程相互作用初探［J］.中國海港建設(shè),2008(3):15-17.FANG W,LI K S,ZHANG G R. Analysis of Interaction of Environmental Factors of Beach Engineering with Artificial Islands［J］.China Harbour Engineering,2008(3):15-17.

［10］陸敏,孫志國,孫大鵬.人工島對(duì)海灣水環(huán)境影響數(shù)值研究［J］.海洋工程,2011(3):88-95.LU M,SUN Z G, SUN D P. Numerical study of the impact of water environment by artificial island in the gulf［J］.The ocean engineering,2011(3):88-95.

［11］Blumberg F A,Mllor L G. A Description of a Three-dimensional Coastal Ocean Circulation Model［C］// In Heaps N5.Threedimensional Coastal Models.Washington DC: American Geophysical Union, 1987.

［12］MELLOR L G, YAMADA T. Development of a turbulence closure model for geophysical fluid problem［J］. Review of Geophysical Space Physics, 1982, 20(4):851-875.

［13］SMAGORINSKY J. General circulation experiments with the primitive equations. I. The basic experiment［J］. Monthly Weather Review, 1963, 91(3): 99-164.

［14］Bolin B, Rodhe H.A note on the concepts of age distribution and transit time in natural reservoirs［J］. Tellus, 1973 (25): 58-62.

［15］Parker D S, Norris D P, Nelson A W. Tidal exchange at golden gate［J］. Journal of Sanit Engineering Division, ASCE, 1972 (2):305-323.

［16］SUN J, TAO J H. Relation matrix of water exchange for sea bays and its application［J］. China Ocean Engineering, 2006(4):529-544.