猶 爽，張寧川

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

大連海上機(jī)場(chǎng)位于渤海遼東灣南部金州灣海域，是我國(guó)內(nèi)地即將修建的第一座離岸式海上機(jī)場(chǎng)，竣工后將成為目前為止世界上最大的海上機(jī)場(chǎng)。機(jī)場(chǎng)人工島在軟土基礎(chǔ)上填筑形成［1］，規(guī)劃人工島總面積為21 km2，平面布置為矩形，寬3.4 km，長(zhǎng)6.4 km。機(jī)場(chǎng)物流區(qū)布置在機(jī)場(chǎng)人工島西側(cè)海岸，機(jī)場(chǎng)施工過程中擬修建從空港物流中心到機(jī)場(chǎng)人工島西南角的全封閉式施工通道。機(jī)場(chǎng)人工島填海完成后，為保障水質(zhì)環(huán)境，施工通道擬全部或部分拆除。

金州灣距離大連市金州區(qū)約2 km，海灣呈橢圓形，灣口朝向西北，水域面積約為342 km2，海岸線長(zhǎng)65.7 km，灣內(nèi)水深自東向西遞增至約10 m。金州灣是金州區(qū)工業(yè)和生活廢水以及北大河、紅旗河等流域的最終歸處。因?yàn)闄C(jī)場(chǎng)人工島位于近灣底水域，加之施工通道，在灣底形成半封閉水域。該水域水質(zhì)保障成為機(jī)場(chǎng)建設(shè)環(huán)境保護(hù)面臨的挑戰(zhàn)性課題。

為保障大連海上機(jī)場(chǎng)建成后金州灣海域水流形態(tài)和污染物輸移順暢，進(jìn)行海上機(jī)場(chǎng)施工通道布置的水動(dòng)力優(yōu)化研究，期望通過優(yōu)化海上機(jī)場(chǎng)施工通道開口寬度，盡可能減小工程建設(shè)對(duì)金州灣海域的水動(dòng)力及環(huán)境影響?；诖四康?，以工程水域物理自凈能力、水體半更換期、局部污染物擴(kuò)散軌跡、局部污染物示蹤粒子第一次到達(dá)灣外的時(shí)長(zhǎng)等參數(shù)為評(píng)價(jià)指標(biāo)，建立非結(jié)構(gòu)網(wǎng)格下的水動(dòng)力及污染物輸移擴(kuò)散數(shù)值模型，對(duì)大連海上機(jī)場(chǎng)施工通道布置方案進(jìn)行水動(dòng)力學(xué)優(yōu)化研究。

1 水動(dòng)力模型的建立及驗(yàn)證

1.1 模型計(jì)算范圍及網(wǎng)格劃分

金州灣位于遼南半島前端黃渤海分界線附近的渤海一側(cè)水域，海灣呈橢圓形，灣口朝向西北開放。漲潮水動(dòng)力主要來自北黃海漲潮流通過遼南半島老鐵山水道的水流動(dòng)力；落潮水動(dòng)力主要來自渤海遼東灣落潮流動(dòng)力。灣內(nèi)潮流受進(jìn)出灣口的水動(dòng)力控制。灣口2個(gè)岬角地理坐標(biāo)分別為39°03′30"N、121°22′25"E和39°11′42"N、121°34′50"E（圖 1）。

人工島工程位于金州灣灣底，距南岸線最短距離僅2 km，距東岸線最短距離約4 km，長(zhǎng)軸走向?yàn)?31°～311°。人工島西南角與西側(cè)空港物流中心由施工通道連接（圖 2）。

針對(duì)上述地理形勢(shì)和工程布置方案，考慮到開邊界對(duì)工程所在海域潮流計(jì)算的影響，將模型計(jì)算區(qū)域開邊界設(shè)置在距離工程區(qū)域西側(cè)約30 km附近的2個(gè)岬角南北垂直連線處。計(jì)算范圍東西跨度約58 km，南北跨度約55 km，計(jì)算范圍的經(jīng)緯度坐標(biāo)為121.266°～121.927°E；38.950°～39.446°N。

為更好地?cái)M合岸線，網(wǎng)格劃分采用三角形網(wǎng)格單元［2-4］，并對(duì)重點(diǎn)關(guān)注區(qū)域進(jìn)行加密，網(wǎng)格單元總數(shù)為19 491，節(jié)點(diǎn)總數(shù)為10 187，最小網(wǎng)格步長(zhǎng)為100 m。

圖1 金州灣地理位置及測(cè)點(diǎn)布置Fig.1 Geographic position of Jinzhou Bay and layout of sites

圖2 機(jī)場(chǎng)人工島及物流中心平面規(guī)劃布置Fig.2 Plane layout of artificial island and logistics center of airport

1.2 計(jì)算模型選取及參數(shù)設(shè)置

計(jì)算模型選取MIKE21中的flow model FM模型。該模型采用淺水連續(xù)性方程和動(dòng)量方程作為控制方程，適用于近岸海域下的水動(dòng)力計(jì)算。本研究不考慮溫度和鹽度影響，模型為正壓模型。

計(jì)算時(shí)，模型渦黏性系數(shù)通過Smagorinsky方程確定，取Smagorinsky系數(shù)Cs為0.28；通常情況下阻力系數(shù)Cf可由曼寧系數(shù)M或謝才系數(shù)C確定，本模型采用曼寧系數(shù)來決定底摩阻，鑒于計(jì)算區(qū)域范圍內(nèi)水深和地質(zhì)類型變化不大，故將整個(gè)計(jì)算范圍內(nèi)的曼寧系數(shù)值設(shè)置為常數(shù)，取值為50；MIKE21采用“凍結(jié)法”處理邊界，通過定義臨界水深來確定干濕單元［5］，本研究中干濕單元的臨界水深取為0.05 m。

模型計(jì)算歷時(shí)長(zhǎng)度取為34 d，啟動(dòng)時(shí)間設(shè)為1 d，模型時(shí)間步長(zhǎng)為10 s。計(jì)算得到的流場(chǎng)如圖3所示。

圖3 金州灣海域流場(chǎng)Fig.3 Flow field of Jinzhou Bay

1.3 模型邊界條件及初始條件

模型計(jì)算區(qū)域驗(yàn)潮站有限，難以提供開邊界資料，本文采用大小區(qū)域嵌套法獲取開邊界潮位。大區(qū)域計(jì)算范圍選取潮位資料豐富的整個(gè)渤海海域，驗(yàn)潮站位包括煙臺(tái)、蓬萊港、龍口港、老虎灘、長(zhǎng)興島等10個(gè)站點(diǎn)，開邊界為威?！￠L(zhǎng)山島連線及其延長(zhǎng)線與岸線相交部分，大范圍模型計(jì)算結(jié)果與實(shí)測(cè)潮位吻合較好，可以為小范圍模型提供可靠的開邊界條件。初始水位及流速均設(shè)置為0。

1.4 模型驗(yàn)證

選取靠近岸邊的棋盤磨潮位驗(yàn)證點(diǎn)（實(shí)測(cè)資料為2009年11月7日～2009年12月11日）和1#～6#潮流驗(yàn)證點(diǎn)分別對(duì)模型的潮位和潮流計(jì)算結(jié)果進(jìn)行驗(yàn)證，各驗(yàn)證點(diǎn)布置如圖1所示。潮流驗(yàn)證分大潮時(shí)段（2009年12月3日～4日）和小潮時(shí)段（2009年11月27日～28日）。

驗(yàn)證結(jié)果如圖4所示，這里潮流驗(yàn)證以2#點(diǎn)和5#點(diǎn)作為示例，從圖4中可以看出，計(jì)算值和實(shí)測(cè)值吻合較好，具有良好的一致性，由于本模型沒有考慮溫度、鹽度及風(fēng)的影響，在個(gè)別時(shí)刻實(shí)測(cè)值和計(jì)算值具有一定差別，但總體符合較好，且能正確地刻畫漲落潮的情況，通過潮位和流場(chǎng)的驗(yàn)證結(jié)果表明，模型是合理的，模擬結(jié)果可以反映計(jì)算海域的水動(dòng)力場(chǎng)特征。

圖4 各站點(diǎn)潮位和潮流驗(yàn)證Fig.4 Tidal level and flow verification of each site

2 施工通道開口寬度對(duì)水環(huán)境的影響

2.1 施工通道開口寬度設(shè)定

工程水域人工島—施工通道—海岸之間形成的半封閉水體與外海的交換，主要取決于施工通道的開口寬度和人工島—海岸之間形成的流口寬度B(圖2)。在人工島寬度和走向確定的條件下，人工島—海岸之間形成的流口寬度B為一確定值。為了使海灣自凈能力、局部水體交換率等指標(biāo)達(dá)到最優(yōu)值，擬將連接人工島與陸地的施工通道從靠近人工島一側(cè)進(jìn)行開挖，并確定最佳開挖寬度。

施工通道總長(zhǎng)度約3 000 m（記為L(zhǎng)，圖2）。為了確定最佳開口寬度，設(shè)定5種開口方案（表1），開口寬度分別為施工通道總長(zhǎng)的1/10、1/3、1/2、2/3以及將施工通道全部開挖，即開口寬度分別為300 m、1 000 m、1 500 m、2 000 m及3 000 m。

表1 施工通道開口方案Tab.1 Schemes of opening-size of construction access

2.2 施工通道開口寬度對(duì)水體交換情況的影響

采用濃度對(duì)流擴(kuò)散模型研究工程水域的水體交換情況。在潮流模型的基礎(chǔ)上，添加對(duì)流擴(kuò)散模塊，求解對(duì)流擴(kuò)散方程得到水體中指標(biāo)物質(zhì)的濃度分布，計(jì)算時(shí)不考慮污染物降解。

在工程水域人工島—施工通道—海岸之間形成的半封閉水域網(wǎng)格點(diǎn)上均布濃度為1 mg/L的保守污染物，邊界污染物初始濃度設(shè)定為0，計(jì)算污染物經(jīng)過30 d輸運(yùn)擴(kuò)散后的濃度分布，并采用物理自凈率和水體半更換期作為水體交換情況的量化指標(biāo)。

（1）物理自凈率。設(shè)經(jīng)過一段時(shí)間i后灣內(nèi)的污染物平均濃度為n%，則物理自凈率定義為

圖5給出了施工通道開口方案四條件下，污染物經(jīng)歷不同時(shí)長(zhǎng)的輸移擴(kuò)散后在全灣水域的數(shù)值計(jì)算結(jié)果示例。其中圖5-a～圖5-f分別為污染物輸移擴(kuò)散5 d、10 d、15 d、20 d、25 d以及30 d后全灣的濃度分布。

圖6匯總給出了不同施工通道開口方案、不同時(shí)間點(diǎn)上的物理自凈率之比較。

圖5 某時(shí)刻污染物濃度分布示例Fig.5 Pollutant distribution in the bay at a certain time

從圖6中可以看出，總體而言，5種施工通道開口寬度的物理自凈率差別不大。具體比較可以看出：施工通道開口寬度與灣內(nèi)污染物擴(kuò)散率并不呈現(xiàn)簡(jiǎn)單的單增或單減關(guān)系，方案一、方案四和方案五的物理自凈率較方案二和方案三略大。這主要是因?yàn)楫?dāng)施工通道開口較小時(shí)，開口處的水流湍急，流速較大，在灣內(nèi)形成環(huán)流，從而有利于灣內(nèi)污染物的輸移；當(dāng)施工通道開口較大時(shí)，盡管此時(shí)開口處的流速相對(duì)較小，但由于開口寬度大，流量也隨之增大，增強(qiáng)了灣內(nèi)水體與灣外水體的交換。

（2）水體半更換期。類比于化學(xué)元素的半衰期，把灣內(nèi)污染物總量降低到起始污染物總量一半所需的時(shí)間定義為水體的半更換期，表明此時(shí)灣內(nèi)水體已被更換一半［6］。根據(jù)此定義，分別統(tǒng)計(jì)各方案的半更換期，分別為 313 h、350 h、361 h、299 h和288 h。其中，方案五的半更換期最短，其次為方案四。

圖6 不同方案之間物理自凈率之比較Fig.6 Comparison of physical self-purification rate among different schemes

圖7 追蹤粒子釋放點(diǎn)Fig.7 Release position of tracer particles

2.3 對(duì)局部污染物擴(kuò)散規(guī)律的影響

灣內(nèi)局部范圍的污染物擴(kuò)散規(guī)律直接影響到部分海域的環(huán)境質(zhì)量，灣內(nèi)水體的交換情況僅反映海灣整體水交換情況，并不能反應(yīng)局部區(qū)域的水體質(zhì)量，因而有必要對(duì)局部污染物的擴(kuò)散規(guī)律進(jìn)行討論。

施工通道掩護(hù)的近岸海域流速較小，若在該海域排放污染物，將對(duì)海域水環(huán)境造成較大不良影響，此處著重關(guān)注污染物近岸排放的情況。為追蹤污染物擴(kuò)散輸移規(guī)律，在施工通道掩護(hù)的水域內(nèi)釋放5個(gè)追蹤粒子，考慮污染物近岸排放，追蹤粒子均在靠近海岸500 m范圍內(nèi)釋放，粒子釋放地點(diǎn)如圖7所示。根據(jù)前人研究結(jié)果［7］，在漲急時(shí)刻投放粒子更為合理，故本模型釋放追蹤粒子的時(shí)刻為2009年11月10日3：00的漲急時(shí)刻。

限于篇幅，圖8中以各方案中3#粒子運(yùn)動(dòng)軌跡作為示例。比較圖8中各圖可以看出：對(duì)應(yīng)開口寬度較小的方案一和方案二，施工通道掩護(hù)的近岸水域流速極小，出現(xiàn)近乎死水區(qū)，該局部污染物輸移非常緩慢，長(zhǎng)時(shí)間不能排出；對(duì)應(yīng)方案三，落潮時(shí)可以將污染物較順利地排出灣外，但漲潮時(shí)將已經(jīng)排出灣外的污染物又重新帶回灣內(nèi)，易造成局部潮流通道水域反復(fù)污染；對(duì)應(yīng)方案四，落潮時(shí)排出到灣外的污染物漲潮時(shí)不會(huì)再帶回灣內(nèi)；對(duì)應(yīng)方案五，由于其開口最大，流速最小，落潮流速幾乎不足以將污染物帶出灣外。

綜上所述，開口寬度為2/3 L時(shí)，更有利于局部污染物的排放，即從施工通道掩護(hù)的近岸局部水域污染物輸移情況看，方案四最佳。

圖8 施工通道掩護(hù)的近岸水域污染物擴(kuò)散軌跡Fig.8 Diffusion trajectory of pollutant in offshore area shielded by construction access

表2匯總給出了各方案中示蹤粒子第一次到達(dá)灣外所經(jīng)歷的時(shí)長(zhǎng)。由表2中可見，除方案三中的4#粒子第一次到達(dá)灣外時(shí)間較長(zhǎng)外，從總體上而言，方案三和方案四中各追蹤粒子第一次到達(dá)灣外的時(shí)間明顯短于其他各方案，說明方案三和方案四施工通道掩護(hù)的近岸水域內(nèi)污染物擴(kuò)散輸移速度較其他方案快。對(duì)于單個(gè)方案來講，這一統(tǒng)計(jì)結(jié)果和追蹤粒子釋放的位置之間存在一定的關(guān)系，粒子釋放的位置不同，統(tǒng)計(jì)結(jié)果也存在一定差異，但從整體規(guī)律上來講，方案一和方案二施工通道較長(zhǎng)，其掩護(hù)的近岸海域中靠近海岸的區(qū)域流速較其他方案小，再加上施工通道開口寬度較小，阻礙了追蹤粒子向?yàn)惩獾倪\(yùn)動(dòng)；方案五中由于施工通道全部開挖，開口寬度較大，其中一些粒子借助落潮流能很快到達(dá)灣外，但工程所在海域漲潮流流速大于落潮流流速，因此方案五中追蹤粒子很快返回灣內(nèi)，且在計(jì)算時(shí)間內(nèi)追蹤粒子長(zhǎng)時(shí)間停留在灣內(nèi)。

2.4 方案推薦

在不同的評(píng)價(jià)指標(biāo)下得到的最優(yōu)方案并不完全一致，若以物理自凈能力為評(píng)價(jià)指標(biāo)，方案一、四、五的計(jì)算結(jié)果較為突出；若以半更換期為評(píng)價(jià)指標(biāo)，方案四和方案五較其他方案好；從局部區(qū)域污染物擴(kuò)散軌跡來看，方案四較其他方案理想；根據(jù)局部區(qū)域追蹤粒子第一次到達(dá)灣外的歷時(shí)來看，方案三和方案四較優(yōu)。綜合以上各評(píng)價(jià)指標(biāo)，方案四是較為理想的方案，因而推薦方案四（即開挖2/3 L）為最終方案。

表2 各方案追蹤粒子第一次到達(dá)灣外經(jīng)歷時(shí)長(zhǎng)Tab.2 The time when tracer particles in local area first get out of the bay of each scheme

3 結(jié)語

本文以大連海上機(jī)場(chǎng)工程為背景，基于MIKE21 flow model FM模塊建立了金州灣二維潮流模型和二維污染物輸移擴(kuò)散模型，計(jì)算結(jié)果與實(shí)測(cè)值的對(duì)比表明模型是可靠的。在此基礎(chǔ)上計(jì)算機(jī)場(chǎng)施工通道不同開口寬度下的物理自凈能力、水體半更換期、局部污染物擴(kuò)散軌跡、局部污染物第一次到達(dá)灣外所經(jīng)歷的時(shí)長(zhǎng)，并以此為指標(biāo)確定施工通道開挖全長(zhǎng)的2/3（即2 000 m）時(shí)為最佳方案。本文計(jì)算結(jié)果可為大連海上機(jī)場(chǎng)施工通道布置的進(jìn)一步改進(jìn)提供可靠的依據(jù)。

［1］王諾,顏華錕,左書華,等.大連海上機(jī)場(chǎng)人工島建設(shè)對(duì)區(qū)域水動(dòng)力及海床沖淤影響分析［J］.水運(yùn)工程,2012（4）：5-11.WANG N,YAN H K,ZUO S H,et al.Impact of Exploitation of Dalian Offshore Airport Artificial Island on Regional Hydrodynamic Condition and Erosion-deposition Variation of Jinzhou bay［J］.Port&Waterway Engineering,2012（4）：5-11.

［2］王志力,陸永軍,耿艷芬.基于非結(jié)構(gòu)網(wǎng)格有限體積法的二維高精度物質(zhì)輸運(yùn)模擬［J］.水科學(xué)進(jìn)展,2008,19（4）：531-536.WANG Z L,LU Y J,GENG Y F.High Resolution Two-dimensional Numerical Model of Scalar Transport Based on Unstructured Grid Finite Volume Method［J］.Advances in Water Science,2008,19（4）：531-536.

［3］賴錫軍,曲卓杰,周杰,等.非結(jié)構(gòu)網(wǎng)格上的三維淺水流動(dòng)數(shù)值模型［J］.水科學(xué)進(jìn)展，2006,17（5）：693-699.LAI X J,QU Z J,ZHOU J,et al.3-D Hydrodynamic Model for Shallow Water on Unstructured Grids［J］.Advances in Water Science,2006,17（5）：693-699.

［4］WANG J H,SHEN Y M,GUO Y K.Seasonal circulation and influence factors of the Bohai Sea：a numerical study based on Lagrangian particle tracking method［J］.Ocean Dynamics,2010（60）：1 581-1 596.

［5］許婷.丹麥 MIKE21 模型概述及應(yīng)用實(shí)例［J］.水利科技與經(jīng)濟(jì),2010,16（8）：867-869.XU T.Calculation Principle and Application Example of a Two-dimensional Flow Model-MIKE 21 HD［J］.Water Conservancy Science and Technology and Economy,2010,16（8）：867-869.

［6］胡建宇,傅子瑯.羅源灣潮流及海水半更換期的數(shù)值研究［J］.廈門大學(xué)學(xué)報(bào)：自然科學(xué)版,1989,28（S1）：24-39.HU J Y,FU Z L.Numerical Studies on Half-exchange Period and Tidal Current in Luoyuan Bay［J］.Journal of Xiamen University：Natural science,1989,28（S1）：24-39.

［7］趙亮,魏皓,趙建中.膠州灣水交換數(shù)值研究［J］.海洋與湖沼,2002,33（3）：23-29.ZHAO L,WEI H,ZHAO J Z.Numerical Study on Water Exchange in Jiaozhou Bay ［J］.Oceanlogia et Limnologia Sinica,2002,33（3）：23-29.