澳門機場填海對九澳灣潮流動力及水體交換的影響*

孫 龍，丁 偉，路川藤，羅小峰，張功瑾

(1.中交第三航務工程勘察設計院有限公司，上海 200032；2.南京水利科學研究院，江蘇 南京 210029)

圍填海工程是沿海地區向外拓展土地、解決用地緊張問題的主要辦法，但各種圍填海工程也會帶來諸如污染、影響通航、改變水動力環境等問題[1-4]。長期以來的圍填工程直接改變沿海岸線形態，從而導致水下地形地貌的變化，影響潮流動力條件。另一方面，地形變化與潮流動力的相互作用又將給地區帶來水環境問題[5-7]，國內外關于圍墾工程對周邊水域水動力及水環境的影響由來已久，大多數學者通過實測、物理模型及數值模擬，就工程附近區域潮汐形式、輸沙環境及海床演變過程進行分析和探討[8-12]。

澳門的海洋水文條件復雜、地理位置特殊，由于多年的圍填海工程，其岸線一直變化，區域的水動力條件因此受到影響。許翔向等[13]、張廣艷[14]根據歷年遙感影像、水下地形等基礎資料，利用遙感影像特征分析等手段，研究澳門附近水域匯流區的水流動力環境、泥沙輸移特征、灘槽格局平面變化等，并分析沖淤原因、預測沖淤趨勢，為該地區綜合治理規劃或其他工程措施提供了重要的基礎數據。黎志均等[15]、陳祖華等[16]和韓保新[17]采用潮流數學模型對澳門海域各階段的潮流、潮汐等水動力條件進行了模擬分析。黃希敏等[18]分析澳門水域十字門水道周邊水域流態和沖淤形式，確定了十字門水道北口治導線的平面布置；方神光[19]通過實測水文資料對澳門水域圍墾最為顯著的澳門機場附近潮流動力環境特性進行了詳細分析和探討。

受惠于澳門旅游業的發展，澳門機場停機坪將接近飽和，澳門機場在保持一條跑道的情況下，跑道可以滿足未來航班需求，但停機坪則會嚴重不足，鑒于機場范圍內及周邊現時沒有適合的土地可供利用，考慮圍填澳門機場跑道西側三角區。為維護九澳灣與外界水域的水體交換能力，圍填區域內將預留一條南北連通的水道，保證九澳灣內水環境、水生態現狀以及水質目標的實現。

本文以九澳灣為研究對象，采用數值模擬的手段，基于經過驗證的CJK3D模型，對澳門機場擴建前后九澳灣內潮動力變化進行模擬分析，并對工程前后的水體交換能力進行研究。

1 研究區域

澳門特別行政區位于廣東省珠海市的東南面，東隔伶仃洋與香港相望，距香港約60 km，南臨南海，西由灣仔水道、十字門水道與珠海市灣仔、橫琴相隔，北連珠海市香洲區，地理位置處于北緯22°06′39″～22°13′06″、東經113°31′33″～113°35′43″。澳門由澳門半島、氹仔島、路環島和路氹城組成。澳門國際機場位于澳門氹仔島和路環島東部的淺海區，于1989年12月動工興建，建成于1995年8月。澳門國際機場由跑道區、機坪和航站樓組成，其中機坪和航站樓位于氹仔島，由開山填海形成陸域，跑道區位于氹仔—路環以東海域，為填海形成的人工島，跑道區和機坪之間通過兩座高樁梁板結構的聯絡橋相連接。澳門機場具體填海范圍見圖1，為機場跑道西側的菱形水域。九澳灣位于澳門路環島東北部，屬伶仃洋西岸淺灘區。九澳灣目前以南口和北側機場三角區與外界水域連通，灣內水體受徑流、潮汐及波浪等多種動力影響。

圖1 研究區域位置及周邊水域水深

2 數學模型的建立

2.1 控制方程

(1)

(2)

(3)

(4)

式中：z為潮位；h為水深；H=h+z為總水深；u和v分別為x和y方向上的流體速度；f=2Ωsinφ為Coriolis系數(其中Ω為地球自轉角速率，φ為當地緯度)；g為重力加速度；C為謝才系數；t為時間；Nx和Ny分別為x和y方向上的水流紊動黏性系數；ψ為水質濃度；px和py分別為x和y方向上的水質擴散系數；k為水質降解系數。

向量形式可表示為：

(5)

式中：U=(H，Hu，Hv)T；E為對流擴散項；Ed為紊動擴散項；M為源項。

采用三角形網格對計算區域進行離散，并將單一的網格單元作為控制單元，水深布置在網格頂點，其他物理變量配置在每個單元的中心。將第i號控制元記為Ωi，在Ωi上對向量式的基本方程組(4)進行積分，并利用Green公式將面積分轉化為線積分，得

(6)

式(6)求解主要分4個部分，一為時變項，二為水平對流項，三為底坡項，四為水平擴散項。對流項數值通量可采用Roe格式的近似Riemann解；水平擴散項含有二次項，采用單元交界的平均值來計算通過該界面擴散項的數值通量；有限體積法底坡項若不加任何處理，則會造成靜水的偽流動現象，本文采用“斜底模型”處理底坡項。

2.2 模型范圍

數學模型包括珠江口主要水系網，北側邊界包含珠江各大入海口門區，西側邊界在黃茅海西側附近，東側邊界至香港附近，南側邊界至-40 m等深線附近，見圖2，模型網格總數為32 516個，最大邊長38 565 m，工程水域網格加密，最小網格約5 m。模型外邊界通過全球潮汐預報模型Tide-Process提供。

圖2 數學模型范圍及局部加密

2.3 模型驗證

數學模型驗證采用2017年10月6—7日枯季大潮水情，包括3個潮位站、6條潮流垂線，驗證內容包括潮位、潮流，驗證點位置見圖3。

圖3 數學模型潮流潮位驗證點布置

圖4為潮位驗證，圖5為潮流驗證。潮位計算偏差基本在0.1 m之內，吻合良好；潮流相位一致，漲落潮過程基本符合水體運動情況，偏差大都在10%之內。總體而言，模型能夠準確反映澳門水域水動力特性。

圖4 潮位驗證

圖5 潮流驗證

3 結果分析

3.1 工程對流場的影響分析

根據2016—2018年澳門水道內潮位站(圖3中潮位站B)3 a洪季潮位資料統計，將3 a累積頻率為10%的潮型作為潮流動力及水體交換的控制潮型，即珠江口洪季典型大潮，見圖6。圖7為澳門機場圍填前后漲落急流態及流場的對比，很明顯，澳門水域整體漲潮流流向為西南向東北，潮流自外海進入澳門水域后，在氹仔島的分流作用下，一股向十字門水道運動、另一股向機場方向運動，由于海域中沒有島嶼或工程的阻擋，因此潮流十分平順；落潮時，水流流向與漲潮相反，由東北向西南，十字門水道的落潮流和機場方向的落潮流匯合后向西南方向運動。工程海域的潮流以往復流為主，外海旋轉流特性相對較強。

圖6 模型控制潮型(珠江口典型大潮)

圖7 工程實施前后漲落急時刻流場流態變化

機場擴建工程實施后，機場圍填工程對整體流場及九澳灣內的流態影響不大，總體來看，漲落急流速變化范圍與趨勢總體一致。漲急時，機場跑道北側延伸段在樁基的阻水作用下，漲急流速呈減小趨勢，機場南側延伸段由于水流的繞流作用，東南側水域流速呈增大趨勢，范圍約為機場東北側0.8 km，大部分水域流速增加幅度在0.05～0.10 m/s，局部最大流速增幅在0.15 m/s左右，現有機場跑道南側東西側水域漲急流速呈減小趨勢，機場東側流速減小范圍長度約為1.6 km，機場西側流速減小范圍至九澳航道轉彎處，最大流速減小幅度約為0.20 m/s。預留連接水道處的北端漲急流速總體呈增大趨勢，主要因為該通道較現狀水深增大，但流速增加幅度相對較小，基本小于0.10 m/s。

根據對九澳灣南進口和機場跑道西水域北進口的潮通量統計可以看出(分析斷面位置見圖8)，由于機場擴建工程的阻流作用，連接九澳灣外水域與澳門水道的水道縮窄，九澳灣南北兩側的進口漲落潮通量大幅度減小，均超過40%，見圖9。

圖8 工程附近潮通量統計斷面位置

圖9 工程實施前后典型斷面漲落潮通量變化

3.2 水體交換影響分析

九澳灣作為連接機場南北水道的海灣，其水體交換能力可以由海水的流動性充分反映。通常，關于水體交換能力大多是通過示蹤劑的濃度變化來判斷水體的交換率和交換時間[20]；本文采用水體半交換周期判斷九澳灣內的水體交換能力。半交換周期是指在潮汐和徑流的共同作用下，灣內水體交換出50%到灣外時的時間，在一定程度上表征了海灣的水交換能力。由九澳灣潮流動力特性可見，受潮汐的影響，九澳灣內的水體交換過程會呈現往復性，灣內各處濃度也會隨之呈現一定的往復性特征，本文以灣內第一次濃度達到初始濃度一半所用的時間作為灣內水體的半交換時間進行研究。

在潮流動力模型基礎上，通過模擬示蹤劑濃度的變化來計算水體交換的速率和半交換周期，進而分析澳門機場填海工程建設對九澳灣的水動力影響。以九澳灣的口門為界，在灣內釋放可溶于水的示蹤劑，不可降解，設置灣內示蹤劑的濃度為1，其他水域為0，模擬在珠江口典型洪季大潮(圖6)潮流作用下濃度場的變化，模擬從高潮漲憩開始，持續24 h。

從水體半交換時間的空間分布(圖10、11)可見，澳門機場擴建工程實施后，九澳灣口門附近水體基本不受工程影響，交換能力最強；受平面形態、地形及距離口門遠近的影響，路環電廠北側內灣的水體交換最弱，在工程前后水體交換能力均較弱；受機場擴建工程影響最顯著的區域為九澳灣內路環電廠往南的凹岸水域，由于機場擴建工程在一定程度上阻礙了九澳灣與澳門水道間的水體交換，加之落潮流的削弱，九澳灣內路環電廠往南的凹岸水域水體交換能力明顯下降。

圖10 工程前九澳灣內濃度分布隨時間變化等值線

圖11 工程后九澳灣內濃度分布隨時間變化等值線

從圖12可以看出，九澳灣內的平均濃度隨時間遞增呈現周期性震蕩。統計結果顯示，工程前，灣內水體的平均半交換周期為2.9 h；珠澳人工島建設后，其水體交換能力有所下降，半交換周期增加約2.2 h，達到5.1 h。

圖12 九澳灣內示蹤劑濃度隨時間變化

4 結論

1)澳門機場圍填海工程前后，海域潮流性質沒有改變，主流仍經機場跑道東側北上南下，小股經九澳灣至澳門水道的潮流受工程束窄作用大幅削弱。

2)澳門機場圍填海工程縮窄了九澳灣與澳門水道的水體流通水域，九澳灣與澳門機場跑道間的水域存在大幅流速降低區域，最大減弱幅度超過0.3 m/s，不利于九澳灣內的水體流動。

3)澳門機場圍填海工程前后，九澳灣內的平均濃度均隨時間遞增呈現周期性震蕩降低的態勢，工程的阻水效應使得九澳灣水體交換能力下降。在珠江口典型大潮水文條件下，灣內水體半交換周期由2.9 h增加至5.1 h。