海平面上升對珠江河口東四口門鹽水入侵的影響

劉忠輝，宏波

(華南理工大學，廣東廣州510641)

鹽水入侵是由于河口鄰近海域的高鹽水團隨潮汐漲潮流沿著潮汐通道向上游推進，咸淡水混合使河口上游區域水體變咸的現象，會對沿途區域的水環境、農業灌溉和生活供水等產生重要影響。珠江三角洲城市群擁有城鎮人口5 500萬，已成為全世界最大的城市片區[1]。該區域約占中國國土面積0.5%的土地，但其生產總值卻占中國總量的近20%[2]，是中國最重要的經濟區域之一。在全球海平面上升的背景下，珠江三角洲地區已成為全球受鹽水入侵嚴重影響的地區之一[3]。從2000年開始，珠江三角洲沿海城市幾乎每年都會受到由于鹽水入侵導致的淡水供應不足問題[4]。2005年12月，珠海、中山和廣州等地均受鹽水入侵影響，取水口水閘氯度嚴重超標，最高達7 500 mg/L[5]。

眾多研究表明，海平面上升會強化鹽水入侵。Rice等[6]利用HEM-3D模型對美國東海岸切薩皮克灣的鹽水入侵長度進行模擬預測，研究發現海平面上升1 m將使得10 ppt等鹽線向上游推進10 km，這將改變現有的水質梯度，威脅當地飲用水供應。Hong等[7]利用EFDC模型研究美國東海岸切薩皮克灣的河口鹽度與物質輸運過程對海平面上升的響應，模擬結果顯示，海平面上升將使水體分層加強，同時引起鹽度上升，鹽水入侵長度增加，并且影響程度與淡水流量有明顯關系。孫志林等[8]利用FVCOM模型對錢塘江河口進行模擬，以2012年為基準年，模擬結果表明海平面上升1 m的情況下，大潮期0.45 ppt等鹽線向上游推進了6 km，小潮期推進了9.8 km。陳維等[9]預測海平面上升1 m時，1月長江河口的南支受北支倒灌過來的咸水的影響將明顯增大，南支水域平均鹽度均超過0.45 ppt。

海平面上升對珠江河口鹽水入侵的影響也非常顯著。何用等[10]利用磨刀門咸潮整體物理模型試驗，對海平面上升引起的咸潮運動變化特點進行分析，發現海平面上升使磨刀門河道沿程鹽度均增大；以2012年的海平面為基準，每升高10 cm最大咸界的上移距離大致為6.8 km。關帥等[11]利用FVCOM模型預測得海平面平均每上升1 cm，磨刀門水道鹽水入侵長度將增加0.85 km。

前人對珠江河口磨刀門鹽水入侵的研究報告較多[10]，而對東四口門的研究相對較少，且東四口門的鹽水入侵問題將影響到廣州、東莞、佛山、中山等重要城市超過3 000萬居民的生產生活。鑒于此，本文采用三維水動力模型EFDC(Environmental Fluid Dynamic Code)模擬珠江河口區域的水動力過程，以2007年為基準年，通過設定多種海平面上升幅度的模型計算模式，對珠江河口的鹽水入侵情況進行模擬預測，分析珠江河口東四口門(虎門、蕉門、洪奇瀝、橫門)鹽水入侵受海平面上升的響應，為相關部門應對海平面上升背景下該區域淡水資源管理、利用提供有效參考。

1 研究方法

1.1 研究區域

珠江河口是位于中國南海北部的一個大型河口，擁有范圍廣闊且復雜的河網，共有8個入?？陂T，分別是：虎門、蕉門、洪奇瀝、橫門、磨刀門、雞啼門、虎跳門和崖門，形成“三江匯流，八口入?！钡乃堤卣鱗11]，而東四口門是指八大口門中向東注入伶仃洋的4個口門[12]。珠江河口位于亞熱帶季風氣候區，夏季高溫多雨是洪季，冬季溫和少雨是枯季。其承接西江、北江、東江三江的徑流量，多年平均總徑流量為281.1億m3[13]。

為了得到真實的水動力過程，本文將模擬區域由珠江河口擴展至其鄰近海域，模擬網格見圖1。

a)模型網格

b)珠江河口圖1 珠江河口及鄰近海域的模擬網格

珠江河口的東四口門區域見圖2，該區域河網復雜，有多條支流，其中R1—R12是東四口門的12條主要支流，圖上還標注了東四口門的4個自來水廠的位置，即新塘水廠、石溪水廠、沙灣水廠和東涌水廠。

圖2 東四口門及其主要河網

1.2 研究方法

本文采用的三維水動力模型EFDC是由美國威廉瑪麗大學維吉尼亞海洋科學研究所John Hamrick等人開發的三維地表水水質數學模型[14]，可實現河流、湖泊、河口和海洋等水體的水動力和水質模擬，是一個多參數有限差分模型，已經在很多河口海岸地區有過成功應用[15-17]。

模型利用曲線正交網格模擬珠江河口及其鄰近海域，垂向上有10個σ層。模型以珠江河口2007年的自然條件為基準，包括地形、風場、徑流量等。潮汐諧波數據由美國俄勒岡州立大學潮汐預測軟件(OTPS)算得，并用于模型開邊界的驅動。開邊界的溫度和鹽度數據取自世界海洋數據庫，風力數據來源于香港機場觀測站(圖1b)采集的2007年風場數據。河流徑流量數據采用珠江三條主要支流的上游各站點2007年的徑流量觀測數據，即石角站(北江)、高要站(西江)和博羅站(東江)。各個淡水源的流量分配按照珠江流域水資源保護局公布的比例進行，見表1。

近幾十年來，全球海平面一直呈現上升趨勢，并且海平面上升的速度不斷增加。根據美國Climate Change Science Program(CCSP) 2009年工作報告指

表1 珠江河口各口門徑流量分配比例

出，到2100年海平面可能升高1.0～1.1 m[18]。Intergovernmental Panel on Climate Change(IPCC) 預測到21世紀末海平面將上升0.82 m[19]。從規避風險的角度來說，對未來海平面上升估計的結果應該更偏重于預測上升值的99分位[20]，即區間內至少有99%的概率包含真值?；诘湫蜐舛嚷窂絉CP8.5，到2100年全球平均海平面上升的99分位是0.4～1.8 m[21]。故本文以2007年為基準，設定5種模式來研究珠江河口全年對不同海平面上升幅度的響應，模式命名分別為Base、P30、P50、P100、P150，分別代表基準情況，海平面上升30、50、100、150 cm的情況。各模式之間除了海平面高度不同以外，其余設置均相同，各模式從河道上游輸入的徑流量都按照2007年珠江河口的實際徑流量進行。

1.3 模型驗證

將珠江河口2007年的模型結果與實際觀測數據進行比較，驗證模型的可靠性。

1.3.1水位驗證

取珠江河口區域的8個站點(三灶、香港、MO1、九洲港、內伶仃、桂山島、燈籠山、黃埔，圖1b)作為水位驗證點，所使用的水位驗證實測數據均為2007年在各站點的實際觀測數據。本文利用Wilmott[22]提出的Skill值來反映模型模擬數據與實測數據之間的匹配程度。Skill按下式計算：

若模型模擬結果與實際觀測吻合越好，Skill值越接近于1，反之趨近于0。本文所選的8個驗證點的Skill值見表2。結果表明模型的Skill值與1相近，模型模擬結果與實際觀測吻合程度較好。

表2 所選驗證點的水位驗證數據

各站點2007年12月的水位實測值和模擬的水位值隨時間的變化過程曲線見圖3。各個站點的水位實測值和模擬值都非常接近，模型模擬的水位與實際水位能較好吻合。

a)三灶 b)香港

c)MO1 d)九洲港

e)桂山島 f)燈籠山

g)黃埔 h)內伶仃

圖3 各個站點水位驗證曲線

1.3.2潮流驗證

以圖1b中的MO1站點作為潮流驗證點，把驗證點處的2007年7月第1周的模擬數據與2007年該站點實測潮流數據進行對比。驗證點水深8.3、9.3 m處的潮流矢量變化情況見圖4。圖中可以看出，8.3、9.3 m處的實測潮流與模擬潮流的流速與流向都吻合良好。

a) 水深為8.3 m

b)水深為9.3 m

圖4 MO1站點的潮流驗證

1.3.3鹽度驗證

以圖1b中的MO23站點作為潮流驗證點，驗證該處在珠江典型枯水期模型模擬的鹽度變化情況與實際觀測結果是否相符。具體驗證的時間為2007年第359.5～361.5天，驗證點表層、中層和底層鹽度變化見圖5，驗證結果反映模型模擬的鹽度在表層、中層和底層都與實際觀測的鹽度非常吻合，模型模擬效果良好，滿足使用要求。

a) 表層圖5 MO23站點的鹽度驗證

b) 中層

c) 底層

續圖5MO23站點的鹽度驗證

2 結果

2.1 鹽度變化

珠江河口各個口門的上游就是珠江三角洲城市群，若出現鹽水入侵將對周邊區域居民的生產生活造成嚴重影響。模型模擬的結果包含整個模擬區域全年的4D鹽度數據。通過對鹽度進行篩選，對各口門1年內鹽度超過1 ppt的天數進行統計，得到了各個口門在不同海平面上升幅度情況下，鹽度超過1 ppt的總天數，見圖6。海平面上升對珠江河口各個口門的鹽度有明顯的影響，海平面上升幅度越大，鹽度超過1 ppt的天數越多。海平面上升1 m時，各個口門的高鹽天數相比基準情況增多的天數統計見表3。

圖6 天數統計

口門增加天數/d虎門63蕉門103洪奇瀝60橫門39磨刀門110雞啼門84虎跳門155崖門153

圖6和表3中的結果表明，在相同海平面上升幅度的影響下，不同口門的響應情況是不同的，其中虎跳門和崖門對海平面上升最敏感，海平面上升1 m時高鹽天數相比基準情況分別增加155、153 d。通過研究各口門鹽度超過1ppt情況的出現時間，發現在海平面上升的情況下，各口門高鹽現象仍是多發生于春、冬兩季，也是珠江的枯水期，這與基準情況相似，但高鹽現象的持續時間會隨海平面上升幅度的增加而延長。

當自來水廠取水口處水體的鹽度超過自來水取用標準(0.5 ppt)時，自來水廠將可能采取停止自來水供應的措施。珠江上游區域4個自來水廠，水廠位置見圖2，全年最高鹽度和超標天數的統計見表4。將海平面上升1 m的模擬結果與基準情況進行對比可以看出，各個水廠鹽度受海平面上升影響非常明顯，鹽度超標現象開始出現，尤其是東涌水廠的最高鹽度將增加7.84 ppt，同時鹽度超標天數增加了114 d。

表4 天數統計

2.2 鹽水入侵長度變化

珠江河口各口門上游區域河流分叉較多，單以某一條支流的鹽水入侵長度來代表一個口門的鹽水入侵并不充分，故將各個口門的主要支流(圖2)都作為研究對象進行分析。模擬的珠江河口在不同海平面上升幅度情況下，2007年1月各主要支流平均的0.5 ppt等鹽線位置見圖7。海平面上升幅度越大，0.5 ppt等鹽線向上游推進越多，鹽水會入侵到越上游區域，加劇鹽水入侵。

圖7 平均0.5 ppt 等鹽線所在位置

海平面上升對珠江河口的鹽水入侵有明顯的強化作用?；鶞是闆r下，洪季各口門都不出現鹽水入侵，但當海平面上升150 cm時，即使在洪季，各口門也會出現不同程度的鹽水入侵問題。把圖2中標定的各口門位置作為計算鹽水入侵長度的起算點，以各支流0.5 ppt等鹽線沿河道至對應口門的距離作為鹽水入侵長度。海平面上升100、150 cm時，東四口門各支流鹽水入侵長度相比基準情況向上游的最大推進距離見表5。由于蕉門(R7、R8)、洪奇瀝(R9、R10)和橫門(R11、R12)各處鹽水都未入侵到各自的河流分叉點(圖2中A、B、C為分叉點)，故這3個口門各自的支流的入侵情況相同。表 5的結果可以看出海平面上升1 m時，各支流的0.5 ppt等鹽線相較于基準情況都向上游推進了較大距離，其中R2的推進距離最大，達21.371 km；蕉門最大推進9.644 km；洪奇瀝最大推進9.753 km，橫門最大推進4.819 km。

表5 0.5 ppt等鹽線相比基準情況最大向上游推進的距離 km

3 討論

受Yang等[23]和Chen等[24]將海平面上升幅度與鹽度變化通過數據回歸的方式得到擬合度較好的直線或曲線的啟發，并考慮到本文模型設定的幾個計算模式之間，僅有海平面上升幅度一個變量，便將各條支流鹽水入侵長度進行月平均，與海平面上升幅度進行線性回歸，得到海平面上升幅度與各支流鹽水入侵長度的擬合直線。選擇入侵情況最嚴重的1月份作為線性回歸分析對象，得到的回歸直線(圖8)可以用于預測未來海平面上升某一特定幅度時，珠江河口各口門上游區域鹽水入侵的距離。每條回歸直線的誤差范數見表6。誤差范數都很小，則回歸直線可靠，說明珠江河口各支流由于海平面上升幅度變化引起的鹽水入侵長度變化與海平面上升幅度之間存在著較好的線性關系。得到的線性回歸關系可用于預測特定海平面上升幅度下某支流的鹽水入侵長度，指導工程實踐，如自來水廠遷址、選址等。

表6 R1-R12線性回歸直線誤差范數

a)虎門R1 b)虎門R2 c)虎門R3

d)虎門R4 e)虎門R5 f)虎門R6

g)蕉門R7-R8 h)洪奇瀝R9-R10 i)橫門R11-R12

圖8 海平面上升幅度與鹽水入侵長度的回歸直線

4 結論

海平面上升會加劇珠江河口東四口門的鹽水入侵問題。在海平面上升1 m的情況下，東四口門高鹽天數(鹽度≥1 ppt)都有明顯增加，其中蕉門增加最多，達3個多月。而建在東四口門上游的自來水廠也會出現由于海平面上升引起的取水口處鹽度超標的問題，其中新塘水廠和東涌水廠受到的影響最大，海平面上升1 m情況下，全年鹽度超標天數將增加3個多月，所以在未來各個水廠必須要考慮遷址。

模型模擬結果還表明海平面上升會使東四口門的鹽水入侵長度增加，其中虎門各支流鹽水入侵長度增加最多，海平面上升150 cm情況下，0.5 ppt等鹽線最多將向上游推進27 km。通過對珠江河口東四口門主要支流在一月份的平均鹽水入侵長度與海平面上升幅度進行線性回歸，得到了可靠的回歸直線。借助回歸直線，可以預測未來海平面上升特定幅度時，珠江河口各支流的鹽水入侵長度，可以為未來水廠遷址提供依據。