土地利用類型對濱海城市復合淹沒的影響
——以象山縣為例?

施劭瑜， 楊 波， 江文勝， 于 格， 王有霄

(1. 中國海洋大學海洋環境與生態教育部重點實驗室， 山東 青島 266100； 2. 中國海洋大學環境科學與工程學院， 山東 青島 266100；3. 中國海洋大學海洋與大氣學院， 山東 青島 266100； 4. 中國科學院地理科學與資源研究所， 北京 100101)

風暴潮是最嚴重的海洋災害，世界上許多沿海地區經常受到熱帶氣旋引發的風暴潮的威脅。近年來，受全球氣候變化影響，部分地區強熱帶氣旋活動顯著增加[1]，臺風導致的風暴潮災害不斷增多。臺風導致的風暴潮不僅破壞碼頭、港口和堤防，而且還會大面積的淹沒沿海地區，從而給沿海居民、農田和水產養殖場帶來巨大損失[2]。由于中國東南沿海地區所瀕臨的西北太平洋熱帶風暴頻發，加之該區域人口、財富、資源等要素密集，因此成為世界上受臺風風暴潮災害影響最為嚴重的地區之一。2010—2019年期間中國東南地區共發生71次由臺風導致的風暴潮，造成36人死亡，經濟損失達865.9億元人民幣，占全國海洋災害損失的86.5%[3-5]。

在全球范圍內許多學者對風暴潮的數值模擬進行了研究，這些研究通常集中于提高水動力模型的精度、穩定性和可靠性上[6-11]。在濱海地區的風暴潮漫灘過程中影響因素很多，除了大氣強迫、波浪以及降雨徑流相互作用外[12]，還受該區域地表粗糙度的影響。地表粗糙度是影響地表水體流動的重要參數，僅次于地形[13-14]，在淹沒區里，樹木、草、灌木和建筑物等地面障礙物對水流施加的阻力會耗散水流和波浪的能量，從而影響流速與流向，最終影響風暴潮漫灘模擬的結果[15-16]。在數值模型中，通常將這種地物的阻礙作用參數化，并以底摩擦系數的形式體現[13]，如根據下墊面土地利用類型來選取不同的曼寧系數。

郭洪琳[17]在研究中考慮了土地利用類型所致的底摩擦分布差異，使用ADCIRC(Advanced circulation model for oceanic，coastal and estuarine waters)模型模擬臺風“弗雷德”(9417)事件中溫州地區的風暴潮漫灘過程。張露[18]在寧波市北侖區街區尺度的風暴潮漫灘數值模擬中同樣考慮土地利用類型帶來的底摩擦變化，刻畫了風暴潮漫灘的影響。但濱海區域在遭受臺風和風暴潮襲擊的同時，往往還經歷暴雨過程。在城市區域，上漲的海水阻礙暴雨所致積水通過城市排水系統(地下管網與河流)排出，導致城市出現內澇。因此研究濱海風暴潮漫灘問題時，不僅要考慮底摩擦隨土地利用類型變化的影響，還要考慮暴雨導致的產匯流過程、地下排水系統的水動力過程。

目前已有許多學者利用數值模型研究風暴潮、降雨與排水系統聯合作用下的濱海城市復合淹沒問題[19-20]，但未見同時考慮土地利用類型對淹沒影響的研究，而這對于提高濱海城市風暴潮淹沒模擬的準確性，對指導城市防洪工作有著重要意義。本文將利用Shi等[20]建立的二維風暴潮漫灘一維地下管網耦合模型，研究在考慮城市管網條件下，土地利用類型對風暴潮淹沒的影響。

1 模型簡介及配置

1.1 模型簡介

ADCIRC(Advanced circulation model for oceanic，coastal and estuarine waters)是由美國圣母大學的J.J. Westerink和北卡羅萊納大學的R.A. Luettich于1992年共同開發的有限元數值模型，通常被用來進行風暴潮及漫灘的模擬與風險評估[21-22]。SWMM(Storm water management model)則是由美國環境保護署(EPA)開發的降雨-徑流管理模型，如今被廣泛應用于城市化地區的暴雨洪水模擬和排水系統規劃、分析和設計[20，23]。將ADCIRC模型與SWMM模型通過水位的相互傳遞進行耦合，開發了一種二維風暴潮漫灘一維地下管網的耦合模型，檢驗了水量守恒性，證實了模型模擬漫灘過程的可靠性，并且在臺風“菲特”(1323)期間對象山縣的風暴潮漫灘進行模擬，得到了良好的模擬結果。模型的耦合過程中并未打破原來兩個模型的輸入輸出，所以模型設置均可按照原先模型的說明進行。

模型的配置和檢驗在文獻[20]中有詳細描述，本文僅將模型配置簡要列出，對底摩擦隨土地利用類型變化做詳細敘述。

1.2 模型區域與網格設置

本文選取中國浙江省東部沿海的濱海城市象山縣為研究區域(見圖1)。該縣地勢西北高，微向東南傾斜，三面環山，東南面向大目涂海灣，降雨充沛、易受風暴潮影響，且河流眾多，有著復雜的排水系統。沿岸建有堅固的堤壩將海水阻擋在外，堤壩東西兩端設有兩個可自由開閉的閘門連接著象山縣內部河流與外海。此外，該縣土地利用類型多為城鄉居民用地與耕地，是研究底摩擦變化對復合淹沒影響的理想區域。

整個計算區域的網格采用逐漸加密的非結構化三角網格構建，網格分辨率從開邊界的10 km逐漸增加到象山縣內的10 m，堤壩在ADCIRC模型中被設置為可被漫過的內部障礙邊界(weir類型)。離岸區域水深數據來源于整合了全球陸地地形和海洋水深的ETOPO1[24]，分辨率為1′，近岸水深和海岸線從中華人民共和國海事局下載的海圖中提取(https://www.msa.gov.cn/)。象山縣的地形數據則來源于高分辨率的數字高程模型(DEM，Digital elevation model)(見圖1)。

圖1 象山縣網格與陸地高程示意圖

象山縣內河道與排水井的空間位置數據均根據象山市政府的排水規劃進行設置[25]，共有456個計算節點與515條管線，并以河流和街道將整個區域劃分為312個子匯水區，堤壩上的閘門被設置為排水系統的兩個始終開啟的出口(見圖2)。降水數據來源于象山縣丹城氣象站實測值。計算中的主要水文參數則參照當地水文數據和經驗值進行設置[26]。

圖2 象山縣排水系統示意圖

1.3 模型驅動場

本研究中在開邊界處使用俄勒岡州立大學開發的潮汐預報軟件OTPS提供的M2、S2、K1、O1四個分潮的調和常數來驅動水位變化(http://volkov.oce.orst.edu/tides)。使用日本氣象廳(http://www.data.jma.go.jp)提供的臺風路徑數據集，利用Jelesnianski的圓形臺風模型[27]計算風場，以提供大氣強迫。

1.4 底摩擦方案

本研究在ADCIRC模型中采用二次底摩擦公式，動量方程中的底摩擦項如下：

(1)

(2)

式中：U，V為深度平均的海水水平流速(m/s)；H=ζ+h，為海水水柱高度(m)；Cf為底摩擦系數，可根據以下形式進行調整：

(3)

式中：g為重力加速度；n為曼寧系數，隨下墊面土地利用類型的不同取不同值；H0為截斷水深(m)，本文取1 m；λ為常數，本文取1/3。

根據CGCS2000投影坐標系下2014年的1∶10 000象山縣及其周圍區域土地利用現狀遙感數據獲取象山縣的土地利用類型，主要被分為五類，即海洋、低密度居民區(含學校、政府機構和交通設施)、裸露的巖石與沙地、常綠林與耕地。根據美國NLCD地物數據分類標準[28]，通過對照分別得到每種土地利用類型相應的曼寧系數，代入公式3中可得相應的底摩擦系數。圖3為象山縣內對應土地利用類型的曼寧系數分布示意圖。可以看出，象山縣除三面環繞的海拔較高的山區林地外，其余區域均以城鄉居民用地為主，主要分布在離岸較遠的西北部以及近岸的東南部分，而中部區域主要以耕地為主，包含小部分分布于丘陵的常綠林。灰色區域為被設置成固壁邊界的高層建筑。

圖3 象山縣土地利用類型(2014)及其對應曼寧系數

通常情況下近海風暴潮的模擬中，Cf的取值在1×10-3左右，且地表淹沒區域下墊面與海底底摩擦系數取值相同。本研究中象山縣下墊面通過土地利用類型修正后，Cf取值在1×10-3至0.324之間，平均值為0.077，底摩擦顯著增大，且位于北部與南部沿海的城區底摩擦增長幅度最大，而多為耕地的中部平原則增長幅度略小。

2 敏感性實驗

通常情況下沿海城市都建有復雜的排水系統與海洋連通，Shi等[20]的研究中詳細討論了排水系統對無堤壩情景下復合淹沒過程的影響。結果表明，排水系統增加了陸地與海洋的連通性，對淹沒的影響受制于地形因素，而降雨則會加劇復合淹沒。在此基礎上，為了探究隨土地類型變化的底摩擦對濱海城市淹沒的影響，我們將其作為敏感因子設計不同條件的兩組淹沒實驗，模擬臺風“菲特”(1323)期間象山縣的復合淹沒過程：(1)無堤壩；(2)有堤壩。

2.1 無堤壩情況下隨土地類型變化的底摩擦對淹沒的影響

在極端的臺風事件中可能會出現部分潰壩或者全段潰壩的情況，這意味著發生風暴潮期間大量海水將淹沒陸地，造成最嚴重的損失。在此情況下，風暴潮是造成淹沒的主要因素，而降雨為次要因素。根據模擬中考慮因素的不同，淹沒過程按表1中的6個情景進行模擬。

表1 各情景模擬中考慮的因素

整個模擬過程中各情景的淹沒面積和淹沒體積時間序列如圖4所示；最大淹沒深度與淹沒時長分布以及差值如圖5所示。統計得各情景的最大淹沒面積與體積如圖6所示。

情景a與b中未考慮排水系統隨土地類型變化的底摩擦使每個漲落潮中淹沒面積與體積的峰值減少，谷值增加(見圖4)，這在風暴潮峰值期間這一現象最為明顯，最終導致淹沒的最大面積與體積分別減少72.0%與51.0%，這是由于象山縣沿海與中部被小丘陵阻擋，但丘陵之中存在兩條地勢略低通道(見圖1)。由圖5的淹沒分布可知，在情景a中風暴潮帶來的海水經此通道流入中部低地勢區域造成大面積淹沒；而在情景b中高底摩擦的低密度居民區與常綠林阻礙了海水的流動，淹沒被阻擋在沿岸區域，但這也導致海水滯留帶來更長時間的淹沒。

情景c與d中考慮了排水系統隨土地類型變化的底摩擦使每個漲落潮過程中象山縣的淹沒面積與體積均減少(見圖4)，但減少幅度對比情景a、b有所降低，最終導致淹沒的最大面積與體積分別減少33.4%與26.4%。進一步將情景b與d對比，在同樣都考慮隨土地類型變化的底摩擦的情況下，排水系統的存在會導致淹沒最大面積與體積分別增加119.5%與49.1%。這是由于部分地表淹沒可以經由排水系統流入城區，導致底摩擦對地表淹沒的阻礙作用相對減弱(見圖5)。

圖4 無壩情況下不同情境模擬過程中象山縣的淹沒面積和淹沒體積變化

圖5 無壩情況下考慮隨土地類型變化的底摩擦的各情景在整個模擬過程中的最大淹沒深度與淹沒時長分布以及與相應未考慮底摩擦變化情景的差值

情景e與f則分別在c與d的基礎上增加了發生在第6～7天的強降雨，以探究隨土地類型變化的底摩擦對降雨所致淹沒的影響。結果表明，底摩擦變化使得降雨所致的最大淹沒面積與體積分別增加12%與5%，這是由于底摩擦增大阻礙了降雨匯入排水系統中，從而在地面堆積。圖5最北部城區的淹沒變化展示了這個現象。

2.2 有堤壩情況下隨土地類型變化的底摩擦對淹沒的影響

在堤壩存在的情況下，風暴潮僅能通過閘門(見圖2)處出海口的水位影響內部排水系統水位，內陸的淹沒主要是由于降雨和排水系統溢流導致，而非風暴潮直接淹沒，因此降雨是造成淹沒的主要因素，而風暴潮為次要因素。根據模擬中考慮因素的不同，淹沒過程按表2中的6個情景進行模擬。

圖6 無堤壩淹沒模擬過程中統計的各情景最大淹沒面積與體積

表2 各情景模擬中考慮的因素

整個模擬過程中各情景的淹沒面積和淹沒體積時間序列如圖7所示；淹沒深度與淹沒時長分布以及差值如圖8所示。統計得各情景的最大淹沒面積與體積如圖9所示。

圖7 存在堤壩的情況下不同情境模擬過程中象山縣的淹沒面積和淹沒體積變化

模擬結果表明，在情景a、b中風暴潮單獨作用僅造成沿岸低洼區域被淹沒。風暴潮峰值來臨前每個漲落潮期間的淹沒范圍小，隨土地類型變化的底摩擦顯著降低淹沒面積與體積(見圖7)；風暴潮峰值期間，底摩擦的阻擋作用減弱，象山縣的最大淹沒面積與體積僅分別減少2.8%與18.3%，且淹沒持續時間也有所減少。這是由于風暴潮峰值期間淹沒止步于南部丘陵，底摩擦對淹沒的影響程度降低(見圖8)。

圖8 存在堤壩的情況下考慮隨土地類型變化的底摩擦的各情景在整個模擬過程中的最大淹沒深度與淹沒時長分布以及與相應未考慮底摩擦變化情景的差值

圖9 有壩淹沒模擬過程中統計的各情景最大淹沒面積與體積

在情景c、d中，淹沒由降雨造成，主要集中在象山縣北部城區與部分中部郊區，這里離岸較遠受海水水位影響小；沿岸地區則淹沒較少。由淹沒分布圖8可以看出，降雨峰值來臨前每個漲落潮期間僅沿岸部分地區被淹沒，淹沒范圍小，隨土地類型變化的底摩擦顯著降低淹沒面積與體積。而降雨峰值期間象山縣北部城區與中部郊區被淹沒，底摩擦在這兩個區域對淹沒呈現不同的影響。北部城區下墊面底摩擦較高，阻礙降雨產流匯入排水系統中，因此導致城區淹沒深度、范圍與持續時間增大；反之排水系統往下游輸送的水量就降低，導致下游淹沒程度減輕；同時下游的中部郊區多為耕地，底摩擦較低，對降雨產流的阻礙作用較弱，在這兩種因素的共同作用下，中部郊區的淹沒深度、范圍與持續時間有所減少。最終，隨土地類型變化的底摩擦使得象山縣最大淹沒面積與體積分別減少18.6%與3.6%，且淹沒持續時間也有所減少。

對于情景e、f，降雨與風暴潮聯合作用導致嚴重的復合淹沒，超出二者單獨作用的線性相加。模擬結果表明，風暴潮峰值與降雨峰值來臨前(0～6 d)，淹沒由風暴潮單獨作用導致，故隨土地類型變化的底摩擦的影響與情景a、b類似，即降低了淹沒程度。而風暴潮峰值與降雨峰值期間，淹沒來源于降雨和風暴潮，故可視為情景a、b與情景c、d中影響的疊加，即北部城區淹沒程度增加，中部郊區和沿岸淹沒程度降低，最終隨土地類型變化的底摩擦使得象山縣最大淹沒面積與體積分別減少20.5%與13.1%，且淹沒持續時間也有所減少。

為了進一步探究隨土地類型變化的底摩擦對復合淹沒風險的影響，根據Bilskie和Hagen[29]的研究，通過計算風暴潮與降雨對淹沒的貢獻，淹沒區可被劃分為三個不同的區域：(Ⅰ)潮汐帶、(Ⅱ)水文帶以及(Ⅲ)過渡帶。情景f的劃分結果如圖10(a)所示，其中潮汐帶、水文帶以及過渡區分別占4.5%、25.9%和69.6%，與情景e對比，潮汐帶和過渡帶占比分別減少3.8%和20.0%，水文帶占比增加23.8%。這是由于底摩擦增大阻擋了沿岸區域的風暴潮淹沒，同時也阻礙降雨匯流過程，加劇城區淹沒，因此降低復合淹沒中風暴潮的影響，加劇了降雨的影響，同時削弱它們聯合作用。此外參照Shi等[20]的研究，過渡帶根據風暴潮和降雨的相對貢獻可被進一步劃分為風暴潮主導、降雨主導與共同主導三種區域(見圖10(b))。隨土地類型變化的底摩擦導致過渡帶中風暴潮主導與降雨主導兩區域占比減少，大部分區域為共同主導區，即過渡帶中風暴潮主導區域向潮汐帶轉換，過渡帶中水文主導區向水文帶轉換，故過渡帶中只剩下風暴潮與降雨共同主導的區域，原因同上。

((a)分為水文帶、潮汐帶與過渡帶，(b)將過渡帶劃分為風暴潮主導區、水文主導區以及共同主導區。 (a) Hydrological zone, tidal zone and transition zone; (b) The transition zone is divided into storm surge dominant zone, hydrologic dominant zone and Equal zone.)

3 結語

本文使用一種1D/2D耦合復合淹沒水動力模型，在底摩擦根據土地利用類型進行調整的情況下，模擬了象山縣多種情景下的淹沒過程，以探究底摩擦系數變化對淹沒的影響。總體而言，通過土地利用類型調整后象山縣的底摩擦增大，導致地表水體的流動受阻，對以降雨為來源的地表淹沒而言，流動受阻意味著降雨產流無法及時的匯入排水系統中，最終導致城區的淹沒加劇。而對于以風暴潮為來源的地表淹沒而言，流動受阻意味著海水造成地表淹沒的速度降低，同時地表淹沒往海里排出的速度也降低。而排水系統里的水體流動不受底摩擦的影響。

以象山縣為例，在自然情況下的風暴潮漫灘中，即不考慮堤壩、排水系統與降雨，底摩擦調整后大大減少了淹沒程度。若在此基礎上考慮排水系統，陸地與海洋的連通性增加，海水不僅僅是通過地表還可以通過地下管網或河流流入。因此底摩擦對海水的阻礙作用減弱，導致淹沒減少的幅度大大降低。而由于在此情景下降雨量相較于海水直接淹沒而言為小量，故底摩擦變化對降雨所致的淹沒影響較小，僅能導致未被風暴潮淹沒區域的降雨致地表匯流受阻，加劇淹沒。而在風暴潮期間，現實的沿海城市復合洪水中，堤壩、排水系統與降雨是必須要考慮的因素。在此情景下底摩擦調整減少了主要受風暴潮影響的象山縣中部與南部沿岸的淹沒，而增加了主要受降雨影響的北部城區淹沒，但同時在這兩種作用的疊加下降雨與風暴潮的聯合作用減弱，故最終導致整體上淹沒程度顯著降低。此外淹沒的分區結果顯示，潮汐帶與過渡帶的占比減少，水文帶的占比增加，且過渡帶中多為風暴潮與降雨共同主導區。這進一步說明了底摩擦調整后削弱了風暴潮與降雨的聯合作用，復合淹沒更趨向于它們各自單獨作用的線性疊加。本研究結果有助于指導濱海城市針對不同區域制定對應的防洪策略，從而減少淹沒損失。