離岸式人工島內河水體紊流數值模擬與換水研究

戚 藍，韓 東，汪彭生，張 亞

（天津大學水利工程仿真與安全國家重點實驗室，天津300072）

離岸式人工島內河水體紊流數值模擬與換水研究

戚 藍，韓 東，汪彭生，張 亞

（天津大學水利工程仿真與安全國家重點實驗室，天津300072）

以福建某人工島為例，進行了內河道換水計算研究。為兼顧景觀和環境，人工島內設置內河道，河道兩端設閘門。根據外海潮位和內河保證水位等水文條件，進行人工島內河道換水計算研究。得到開閘時機、換水時間、不同情況下水體交換量和內河道水體流態圖，并推薦了合理的換水方案。

人工島;內河換水;紊流數值模擬;水體交換

隨著現代社會的快速發展，陸域的土地資源漸趨枯竭，而可以用來開發的近海面積非?？捎^，從海洋開辟新的土地資源成為一種途徑。近年來，人們在近海建設人工島已經有了很多成功的經驗，比如迪拜的棕櫚島［1］、澳門國際機場［2］、迪拜帆船酒店［3］、日本大阪關西國際空港［4］等等，為了兼顧景觀和環境希望在人工島內設置內河。以往人們對于人工島的規劃建設多強調其本身的應用價值或者人工島本身的設計及施工等方面的研究［5］，而針對人工島內部的水體交換的研究比較少。國內方偉［6］等人研究了在海灘工程中環境因子與人工島工程的相互作用；陸敏［7］等人研究了人工島對海灣水環境影響；何杰［8-9］等人曾對港池內的水體交換能力作了一些研究。人工島建設在海灣中，區域潮流流速小，人工島內設有細長內河時水動力很弱，造成換水困難。如果人工島建設完成之后，該人工島內河水體無法得到充分交換，島內水體就會停滯、變質、甚至惡化，那么會給人工島島上的生存環境、人工島的具體管理、運營帶來很大的影響，一些已建工程的失敗教訓是前車之鑒。為避免這種情況的發生，在前期規劃階段，就要對人工島島內水體交換進行充分的研究。

本文以福建人工島為例，本人工島工程是經國務院審批通過的國內首個完全用于旅游、居住和生態開發建設的大型離岸式人工島，與傳統的大面積、平推式圍填海方式所造成的生態破壞、海洋污染、價值單一不同，人工島在區域規劃中引入水體循環系統，兩條交相呼應的海豚狀平面形態構成了內河道親水景觀。人工島內河道兩端與外海銜接處設有閘門，在內河道常年保持有清新的內河水，從而構成人工島獨特的河道景觀和環境保障。為了控制內河水位，保證水質，優化開啟河道兩端閘門時機，利用內外水位差進行水體的自然交換，使得內河水體與外來海水水體摻混充分，將污染物盡可能多的帶出島外，增強人工島內河道水體自凈能力。

現利用實際地形數據建立人工島內河道模型的基礎上，結合VOF［10］方法，對人工島內河道進行了數值模擬，并與試驗結果對比［11］。重點研究了人工島內河道的水體交換，通過研究不同方案對內河道水體交換的影響，在定位人工島功能的前提下，基于人工島島內水體充分交換的規劃理念，提供了合理的換水方案。

1 數值模擬方法

1.1 VOF方法概述

VOF法［10］是Hirt和Nichols1981年提出的處理復雜自由表面的有效方法，是目前應用非常廣泛的一種追蹤自由表面的數值方法。VOF方法的基本原理是通過研究網格單元中流體和網格體積比函數F來確定自由面，追蹤流體的變化，而非追蹤自由液面上質點的運動。VOF方法可以處理自由面重組等非線性現象，所需計算時間短、存儲量少，VOF方法根據體積比函數F來構造和追蹤自由面：若F=1，表示該單元被流體完全充滿；若F=0，表示該單元是個空單元；當0＜F＜1時，表示該單元被部分流體充滿。自由表面的位置由當0＜F＜1時所確定的單元決定。

1.2 數學模型

數學模型采用基于結構化矩形網格的FAVOR方法及真實的3步Tru-VOF方法，控制方程中含有體積和面積分數參數。本次模擬應用GMRES方法求解離散方程。

1.3 控制方程

控制方程采用N-S方程，建立內河道三維RNGκ-ε紊流數學模型。

控制方程包括連續性方程、動量方程、能量方程、紊動能κ方程、紊動能耗散率ε方程。

連續方程：

動量方程：

VOF運動學方程：

式中：u，v，w是在x，y，z三個方向上的流速分量，m/s；代表x，y，z三個方向可流動的面積分數，m2；為x，y，z三個方向的重力加速度是x，y，z三個方向的粘滯力；是可流動的體積分數；ρ是流體密度；p是作用在流體微元上的壓力；流體體積函數F=F(x,y,z,t)，表示計算區域內流體的體積占據計算區域的相對比例。

紊動能k方程：

紊動能耗散率ε方程：

式中：k為紊動能，kg·m2/s2；ε為紊動能耗散率為水體動力黏滯系數，為紊動黏性系數為紊動能κ的產生項分別為湍動能和耗散率所對應的Prandtl數是經驗常數=1.42，

1.4 模型網格劃分及邊界條件

按原型（工程圖紙）1：1建立數學模型，模型范圍為整個人工島，規劃總面積為221.67萬m2，半徑為840 m，人工島內河道底高程為-1 m，形狀如圖1為倒S型，內河道水位隨外海潮位變化范圍為2.5～0.0 m。計算區域內劃分連續式區塊，各區塊內采用結構化正交網格來劃分。最小單元尺寸1.0 m×1.0 m×0.5 m(長×寬×高)，網格總數約2.1×106，具體見圖1。

人工島內河道水體交換先開閘排水、再于外海漲潮時進水。圖2為人工島所在區域流態圖，由于人工島是建設在海灣內，所在水域基本為緩流區，水動力較弱，而且內河道進水、排水又有閘門的限制，潮流對內河道水體交換影響甚微，經計算對比認為分析可以忽略潮流對內河道水體交換的影響，只考慮潮位變化的影響，圖3為外海大潮期潮型曲線，以85國家高程為基準面，本文以此潮型為例模擬，排水時段：外海落潮時段排水，內河道與外海水位一致時(2.5m)開閘門，內河水位隨外海潮位降低至一定水位（0m或1m）后關閘。進水時段：外海漲潮時段進水，內河道與外海水位一致時開閘，內河水位隨外海潮位升高至設計常水位（2.5 m）后關閘。（注：以上水位均以85國家高程為基準面）

對于邊界條件作如下設置：相應的進口和出口邊界分別以實測潮位給定變化水位，其他邊界采用固壁邊界設置。

2 數值模擬及分析

2.1 模型驗證

人工島內河道分為北水道和南水道兩部分，兩條水道基本上呈對稱分布，在模型中南、北水道各布置了3個測流點共6個測點，測點布置見圖4。

圖5為進水時不同時刻d1點的實測流速與數值模擬流速的對比圖，從圖中看出，最大數值模擬流速與實測流速基本吻合，其他時刻流速值相差不到0.02m/s；圖6為人工島物理試驗模型，圖7為數值模擬的內河道出水時流態圖，從圖7可看出南水道流速較北水道大，人工島中間表演場地流速較小，數值模擬與模型實測數據吻合較好；表1為內河道各測點最大實測流速［11］與數值模擬流速的對比情況，基本一致；

總體來看，數值模擬中流速、流態驗證情況良好，數值模擬流速、流態過程與實測物理模型的流速、流態過程基本一致。

2.2 換水時間分析

本次計算得到只開南（或北）閘門排水內河道水位由2.5m下降至1.0m時歷時為67min左右；開兩個閘門歷時為63 min左右，具體見表2和3。結合大潮期潮型曲線（圖3），外海潮位由2.5 m下降到1.0m的歷時60 min，可得出排水時，內河水位與外海潮位的下降基本同步（其中開兩個閘門滯后3 min，開一個閘門滯后7 min）。相應的，進水時內河水位與外海潮位的上升也基本同步。因此人工島內河道和島外潮位過程基本一致，沒有明顯滯后現象出現，模擬結果與物理模型試驗結果一致。

2.3 水體交換分析

評價水體交換的指標很多，例如質點示蹤模型［12］、物質交換率［13］等等。本次研究中通過新進入內河道水體占原有內河道水體百分比來評價人工島內河道的水體交換。一次排水、進水為一個換水周期。定義換水率為［14］

式中：C0為引排水前人工島內河道的背景濃度；C為 引排水過程中某時刻內河道的平均濃度值；Ci為引水進入內河道的水體中某溶解物的濃度。設定內河道某溶解物質A初始含量為1，入庫水體含量Ci=0，則換水率

圖8為不同工況下水體交換圖。其中圖8-a為第一個換水周期進水，從中可以看出，由于內河道窄而長，而且內河道水位隨外海潮位穩定上升，所以內河道進水時外海進入的水體與內河道原水體不會有大的摻混，并且會把內河道水體推向南閘門，原內河道水體會持續雍高到設計水位2.5m（t=108min）后進水結束，此周期換水率為40%左右。

第二個換水周期排水情況。（1）當只開南閘門排水時，從圖8-b中可以看出，外海新進的水體由北門逐漸涌向南門，原河道的水體大部分由南閘門排出，當內河道水體由設計常水位2.5 m下降到1.0m（t=69min）時，換水率為85%左右。（2）當只開北閘門排水時，從圖8-c中可以看出，新進入的外海水體大部分又被排出到外海，原河道水體只有部分排出，南閘門附近水體長時間得不到交換，水體可能會惡化、變質，因此開北閘門達不到換水的目的，此時換水率為50%左右。（3）當開南、北兩個閘門排水時，從圖8-d中可以看出，新進入內河道的水體大部分由北閘門排出，原河道的水體大部分由南閘門排出，當內河道水體由設計常水位2.5 m下降到1.0m（t=69 min）時，換水率為70%左右。

依據數值模擬結果和式（8），圖9給出了第二個換水周期人工島內河道換水率的時間過程線。通過對三種工況的比較可見，開兩個閘門比只開北閘門時換水率高，開南閘門的換水效率較之開北閘門和開兩個閘門明顯提高。

2.4 流態分析

圖10為不同工況下具有最大速度時刻排水流態圖。由圖10-a中可以看出，只開南閘門時，南水道流速明顯比北水道大，中間表演場地左側束窄部位局部流速增大，速度達到0.50 m/s左右，整個區域最大流速為1.00m/s。；從圖10-b中可以得出，只開北閘門時，北水道流速明顯比南水道流速大，中間表演場地右側束窄部位局部流速增大，速度達到0.40m/s左右，整個區域最大流速為0.97 m/s。；從圖10-c中可以看出，開南、北兩個閘門時，南、北閘門附近流速明顯比中間表演場地大，中間表演場地左右兩側束窄部位流速較大，為0.20m/s左右，整個區域最大流速為0.69 m/s。

以上分析可以得出只開南閘門和只開北閘門其最大流速基本一致，開兩個閘門比開一個閘門在相同部位的最大流速小0.3 m/s左右，因此開兩個閘門對于內河道的防護有利，但是由于其流速相對比較小，河道內水動力較弱，因此不利于內河道水體交換；單純從河道水體速度流態來說，開南閘門與開北閘門對換水的影響區別不大，但是從2.3節水體交換得出只開北閘門南水道水體不易排出海外，達不到有效換水的目的，因此建議開南閘門排水。

3 結論及建議

本文依據人工島實際地形數據和實測水文資料建立了人工島三維數學模型，且根據實際情況建立了內河道的進、出口邊界條件，采用VOF方法對內河道水體水流自由表面進行追蹤，并采用RNGκ-ε數學模型對內河道水體進行數值模擬。數值模擬和物理模型的對比表明建立的模型和求解方法能夠較好的處理內河道水體的流動。

水體交換模擬結果表明：第一個換水周期，換水量占原內河道的40%左右；雖然只開北閘門的排水時間和最大流速與只開南閘門沒有很大差別，但是由水體交換圖得出只開北閘門換水率只有50%左右，遠小于其他兩種工況的換水率；開南、北兩個閘門的換水時間比只開南閘門的時間少4 min左右，但是其換水率卻比后者小15%，而且其河道內水動力也比后者小，不利于水體的交換。因此從水體交換的角度建議只開南閘門排水，這樣有利于將污染物盡可能多的帶出島外，增強人工島內河道水體自凈能力并保持內河道長期擁有良好的水質，此時第二個周期換水率約為85%，內河道水體由常水位2.5m下降到1.0 m的排水時間為67 min，內河道最大流速約為1.0m/s左右。

由人工島內河道水體流態分析可以看出內河道與閘室連接的部位流速比其他部位相對要大，如有條件可以加強防護，防止沖刷。中間表演場地周邊速度比較小，注意防止淤積現象的發生。本文結論僅限于本工程。

［1］林琛.迪拜棕櫚島:人間奇跡還是人類的遺憾［J］.生態經濟,2007（2）:8-13.LIN Z.Dubai Palm Island:miracle or human regret［J］.Ecological and economic,2007（2）:8-13.

［2］邱毅平.澳門國際機場海上人工島工程設計簡介［J］.上海水利,1999（3）:28-33. QIU Y P.Introduction to engineering design of the Macau International Airport offshore artificial island［J］.Shanghai water conservation，1999（3）:28-33.

［3］謝浩.迪拜奇跡之帆船酒店［J］.綠色建筑,2010（2）:3-4.

［4］中瀨,明男,呂莉娃.關西國際機場人工島建設［J］.海岸工程,1994,13（Z1）:197-201. ZHONG L,MING N,LV LW.The construction of Kansai International Airport artificial island［J］.Coastal Engineering，1994,13（Z1）:197-201.

［5］謝世楞.人工島設計的進展［J］.海岸工程,1995,14（1）:1-7. XIE S L.The progress of artificial island design［J］.Coastal Engineering,1995,14（1）:1-7.

［6］方偉,李凱雙,張桂榮.灘海工程環境因子與人工島工程相互作用初探［J］.中國港灣建設,2008（3）:15-17. FANGW,LIK S,ZHANGG R.Analysis of Interaction of Environmental Factors of Beach Engineering with Artificial Islands［J］. China Harbour Engineering，2008（3）:15-17.

［7］陸敏,孫志國,孫大鵬.人工島對海灣水環境影響數值研究［J］.海洋工程,2011,27（3）:88-95. LU M,SUN ZG,SUN D P.Numerical study of the impact of water environment by artificial island in the gulf［J］.The ocean engineering,2011,27（3）:88-95.

［8］何杰,葉小強,辛文杰.環抱式單口門港池水體交換能力研究［J］.水運工程,2009（2）:87-91. HE J,YE X Q,XINW J.On water exehange aility for encircled basin with one entrance［J］.Port&Waterway Engineering,2009（2）:87-91.

［9］何杰,辛文杰.挖入式港池水體交換能力研究［J］.水利學報,2007（Z1）:330-333. HE J,XINW J.Study on the water exchange ability for the excavated-in harbor basin［J］.Journal of hydraulic engineering，2007（Z1）:330-333.

［10］張健,方杰,范波芹.VOF方法理論與應用綜述［J］.水利水電科技進展,2005,25（2）:67-70. ZHANG J,FANG J,FAN B Q.Advances in research of VOF method［J］.Advances in Science and Technology of Water Resources,2005,25（2）:67-70.

［11］招商局漳州開發區.人工島物理模型試驗研究成果［R］.南京：南京水利科學研究院,2009.

［12］Awaji T,Imasato N,Kunishi H.Tidal exchange through a strait:a numerical experiment using a simple model basin［J］.Journal of Physical Oceanography,1980，10:1 499-1 508.

［13］Parker D S,Norris D P,Nelson AW.Tidal exchange at golden gate［J］.Journal of Sanit Engineering Division,ASCE,1972，2: 305-323.

［14］陳青毅,吳衛,劉樺.河口水庫庫內三維流動數值模擬［J］.水動力學研究與進展:A輯,2011,26（2）:140-149. CHEN Q Y,WUW,LIU Y.Numerical simulation of three dimensional flowsin estuarine reservoir［J］.Journal of hydrodynamics: Ser.A,2011,26（2）:140-149.

Research on turbulent numerical simulation and water exchange in inner river course of offshore artificial island

QI Lan,HAN Dong,WANG Peng-sheng,ZHANG Ya
(State Key Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300072,China)

In this paper,water exchange in the inner river course for an artificial island in Fujian was studied.To build an artificial island,a river course was set taking into account of the landscape and the environment,with two sluices at the end of the river course.According to the hydrological conditions such as the tidal level of open sea and the guaranteed water level in the inner river,water exchange of the inner river course was analyzed.The occasion of sluice opening,the time of water exchange,the exchange capacity and the flow state figures of river course were obtained,and the reasonable scheme for water exchange was proposed.

artificial island;inland river water exchange;turbulent numerical simulation;water exchange

TV 143；O 242.1

A

1005-8443（2013）05-0403-06

2013-01-08；

2013-01-18

戚藍（1955-），女，山東省人，教授，博士生導師，主要從事水工結構工程及巖石力學方面的研究與教學工作。

Biography：QILan（1955-），female，professor.