張露，傅賜福，董劍希,2，于福江,2

( 1. 國家海洋環(huán)境預(yù)報中心，北京 100081；2. 國家海洋局海洋災(zāi)害預(yù)報技術(shù)研究重點實驗室，北京 100081)

1 引言

我國是世界上少數(shù)幾個海洋災(zāi)害十分嚴重的國家之一，尤其以風(fēng)暴潮災(zāi)害為主，約占總災(zāi)害損失的90%以上，風(fēng)暴潮不僅會對碼頭、港口和堤岸造成破壞，還會在洪水沖毀堤壩后，漫灘淹沒和損壞沿海水產(chǎn)養(yǎng)殖、農(nóng)田和房屋[1]。21世紀以來，每年因風(fēng)暴潮所帶來的經(jīng)濟損失，都高達人民幣數(shù)百億元，我國沿海風(fēng)暴潮災(zāi)害影響出現(xiàn)日益嚴重的趨勢[2–3]。

目前國內(nèi)對風(fēng)暴潮漫灘進行數(shù)值模擬模型的網(wǎng)格分辨率，最高可達到50 m左右。鄭國誕等[4]利用MIKE21FM模塊建立了適用于臺州溫嶺市的高分辨率風(fēng)暴潮漫灘數(shù)值模式，臺州近海及陸上網(wǎng)格邊長約50 m。傅賜福等[5]利用ADCIRC模型建立了一套非結(jié)構(gòu)三角形網(wǎng)格，適用于濱海新區(qū)的高分辨率風(fēng)暴潮漫灘數(shù)值模式，陸地區(qū)域分辨率達到50～80 m，對兩次典型的溫帶風(fēng)暴潮進行模擬。李勇等[6]對渤海灣西部風(fēng)暴潮漫灘進行數(shù)值模擬，計算區(qū)域覆蓋了整個東海海域，對普陀海域進行重點加密，其分辨率為30～200 m。Yang等[7]利用臺風(fēng)“天鴿”（2017年）和臺風(fēng)“山竹”（2018年）臺風(fēng)對澳門的淹沒圖，結(jié)合實測氣象和波浪資料驗證了一套完整的模擬臺風(fēng)期間潮波場耦合的數(shù)值模型，網(wǎng)格最高分辨率為20～100 m。

國內(nèi)尚未有在地形復(fù)雜、街區(qū)密集情況下，基于街區(qū)尺度5～10 m分辨率的風(fēng)暴潮漫灘數(shù)值模擬。但是隨著城市建設(shè)步伐的加快，建筑建設(shè)鱗次櫛比，街巷阡陌縱橫緊密，在近岸和復(fù)雜的海岸地形區(qū)域中使用高分辨率網(wǎng)格，對于精細化的模擬和預(yù)報風(fēng)暴潮漫灘都十分必要。這將是未來風(fēng)暴潮預(yù)報精細化的趨勢，對于制定城市的防災(zāi)減災(zāi)戰(zhàn)略和完善風(fēng)暴潮預(yù)警機制有著重要意義。

本文將利用ADCIRC模型，依據(jù)寧波市北侖區(qū)街區(qū)尺度的高分辨率高程數(shù)據(jù)，建立街區(qū)尺度風(fēng)暴潮漫灘模型，北侖區(qū)陸地分辨率達到5～10 m。利用1211號臺風(fēng)“?？焙?612號臺風(fēng)“溫黛”對街區(qū)尺度風(fēng)暴潮模型進行風(fēng)暴潮模擬驗證。以5612號臺風(fēng)“溫黛”為例，在北侖區(qū)開展基于街區(qū)尺度5～10 m分辨率的風(fēng)暴潮漫灘數(shù)值模擬，并考慮下墊面底摩擦的變化對風(fēng)暴潮漫灘的影響。

2 街區(qū)尺度風(fēng)暴潮漫灘模型的建立

2.1 ADCIRC模型介紹和參數(shù)設(shè)置

本文風(fēng)暴潮模擬采用的是ADCIRC二維模型。其在球坐標系下通過基于垂直平均的原始連續(xù)方程和海水動量方程來求解自由表面起伏、二維流速等3個變量，即 ζ ,u,v。在球坐標系下海水的連續(xù)方程和海水原始動量方程為

式中，λ和φ為經(jīng)度和緯度；ζ為從平均海平面起算的自由表面高度；H為海水水柱的總水深；U和V為深度平均的海水水平流速；R為地球的半徑；η為牛頓引潮勢；g為重力加速度；f為科氏參數(shù)；Ps為海水自由表面的大氣壓強；ρ0為海水密度；Dλ，Dφ為動量方程的水平擴散項；τbλ,τbφ為海底摩擦力分量；τsλ,τsφ海表面應(yīng)力分量。

初始條件為： ζ =u=v=0；海岸邊界條件為：Vn=0這里Vn為 岸邊界的法向深度平均流速。在求解所需物理變量的過程中，空間采用有限元法離散，時間采用有限差分法，將連續(xù)方程和運動方程通過引入了一個時空加權(quán)參數(shù)進行結(jié)合求解[8–9]，提高了計算結(jié)果的穩(wěn)定性。參數(shù)設(shè)置采用冷啟動、球坐標、二維模式、混合底摩擦形式、考慮有限振幅項并采用干濕法、科氏參數(shù)基于β平面近似、考慮潮汐、時間步長為 0.02 s，滿足計算穩(wěn)定和收斂要求。由 K1、O1、P1Q1、N2、M2、S2、K2等 8 個分潮驅(qū)動計算。

2.1.1 氣壓場和風(fēng)場模型

本文選用的臺風(fēng)風(fēng)場模型為Holland（1980）[10]臺風(fēng)氣壓場分布公式為

式中，P(r,θ)為徑向距離r和方位角 θ的函數(shù)，是距臺風(fēng)中心r的海表面氣壓值；Rmax是臺風(fēng)最大風(fēng)速半徑；臺風(fēng)中心氣壓為Pc；臺風(fēng)以外不受干擾的背景氣壓為Pn，設(shè)為 1 012 hPa。

風(fēng)應(yīng)力在ADCIRC中的表達式為

式中，ρa為空氣密度；Cd為拖曳系數(shù)采用Garratt[11]公式

若Cd超過了0.003，則Cd=0.003。式（5）中第1個是代表海面上10 m處風(fēng)速，第2個是海面上10 m處風(fēng)速的大小。

2.1.2 底摩擦設(shè)置

底摩擦的混合形式與二次形式表達式相同，均為

Cf作為底摩擦系數(shù)是變化的?；旌闲问綖?/p>

式中，Cfmin和Hbreak都是常數(shù)；λ和θ也是常數(shù)，λ =1/3，θ=10，因此Cf只隨總水深H的變化而變化。二次形式為

式中，g為重力加速度；n為曼寧系數(shù)；H0為臨界水深，即H0=Hbreak；λ=1/3。

2.1.3 干濕網(wǎng)格處理

風(fēng)暴潮漫灘采用干濕網(wǎng)格法進行計算。此算法在連續(xù)方程之后和動量方程之前進行計算，位于一次時間步長的中間部分。

首先，在模式中設(shè)定一個為較小正數(shù)的最小水深Hmin，計算范圍中總水深大于Hmin的某個格點，模式將該點判斷為濕點。第二，當(dāng)干點變?yōu)闈顸c后，底摩擦力和水位梯度相平衡產(chǎn)生定態(tài)速度U，模式要求設(shè)定的參數(shù)Umin（濕點最小速度）需小于它，為了約束鄰近格點的水位，設(shè)定Umin的值。

2.1.4 風(fēng)暴潮漫灘溢流的計算

水位高度超過堤壩的頂部時，水流將流過堤壩。單位時間內(nèi)的流量表達式為

流量系數(shù)Cd的表達式為

式中，v1為前進流速；g為重力加速度；h1為高于堰頂?shù)乃桓叨取?/p>

2.2 街區(qū)尺度風(fēng)暴潮漫灘模型

寧波市北侖區(qū)地處浙江最東岸，受風(fēng)暴潮和海浪等海洋災(zāi)害影響較大。在歷史上受到臺風(fēng)風(fēng)暴潮災(zāi)害影響次數(shù)較多，臺風(fēng)路徑也有詳細歷史資料，并且在北侖區(qū)有高分辨率的高程數(shù)據(jù)支撐，因此我們選取北侖區(qū)作為實驗區(qū)域。

2.2.1 岸線、水深以及高程數(shù)據(jù)

采用的岸線數(shù)據(jù)為1997–2001年經(jīng)國家測繪局海岸線資料訂正后，比例尺為1∶250 000的海岸線數(shù)據(jù)。寧波市北侖區(qū)以外的水深數(shù)據(jù)是通過范圍5°S～52°N，99°～157°E 之間，分辨率為 2′×2′的水深數(shù)據(jù)，插值得到的。北侖區(qū)近海以及甬江口的水深數(shù)據(jù)采用海圖及實測水深。本文使用的北侖地區(qū)數(shù)字高程模型（DEM）數(shù)據(jù)空間分辨率為5 m，比例尺為1∶10 000，這對于高分辨率的風(fēng)暴潮漫灘模擬提供了重要的基礎(chǔ)數(shù)據(jù)，選取范圍為29.891°～ 29.962°N，121.811°～121.905°E，如圖1所示。

圖1 所選區(qū)域Fig. 1 Selected area

在選取的區(qū)域中，低層建筑多為二層樓房，高程高于5 m的不做處理，低于5 m的高程統(tǒng)一賦值為5 m。高層建筑物所在地高程高于10 m的不做處理，低于10 m的高程統(tǒng)一賦值為10 m。對高層和低層建筑進行劃分如圖2左圖所示，右圖為處理后的效果。

圖2 高低層建筑的劃分示意圖Fig. 2 Division diagram of high and low buildings

2.2.2 網(wǎng)格的建立

本文所選區(qū)域選用SMS v8.1生成街區(qū)尺度非結(jié)構(gòu)三角網(wǎng)格?？紤]到為了較好地模擬大洋潮波、風(fēng)暴潮的傳播過程，減小開邊界對模擬結(jié)果的影響，選取的計算范圍包括東海、黃海和南海部分海域，經(jīng)緯度為 20°～34°N 和 115°～131°E。所選的北侖區(qū)街區(qū)分布密集，城市建筑錯綜復(fù)雜，為達到街區(qū)尺度的分辨率，將北侖區(qū)網(wǎng)格分辨率設(shè)置為5～10 m，沿岸分辨率逐漸向外海增大，能夠滿足對重點地物的刻畫。整套街區(qū)尺度網(wǎng)格包括978 783個三角形網(wǎng)格，共計502 063個節(jié)點，北侖區(qū)陸地占據(jù)了668 162個三角形網(wǎng)格?？紤]到所選臺風(fēng)以及漫灘模擬的要求，不考慮堤壩。

圖3為所選網(wǎng)格區(qū)域、北侖區(qū)附近海域以及北侖區(qū)局部放大的陸地網(wǎng)格示意圖。由于陸地網(wǎng)格分辨率過高，為了滿足CFL條件，時間步長定為0.02 s，但是這樣使計算時間過長，因此將部分陸地高低層建筑的邊界設(shè)置為固體邊界，即圖3中北侖區(qū)陸地網(wǎng)格的空白區(qū)域，減少網(wǎng)格數(shù)量以節(jié)省運算時間。

圖3 大洋所選網(wǎng)格區(qū)域（a）、北侖區(qū)附近海域（b）以及北侖區(qū)局部放大的陸地網(wǎng)格（c）Fig. 3 The selected grid of ocean (a),the sea area near the Beilun District (b),and the land grid magnified locally in the Beilun District (c)

圖4展示了所選區(qū)域、北侖區(qū)附近海域水深分布。圖5為北侖區(qū)高程圖，采用分辨率為5 m的北侖區(qū)高程數(shù)據(jù)進行插值，能進一步細致刻畫北侖區(qū)的城市建筑和街區(qū)分布，清晰的分辨高低層建筑。

圖4 大洋（a）和北侖區(qū)附近海域（b）水深分布Fig. 4 The distribution of water depth in the ocean (a) and nearby seas of Beilun District (b)

圖5 陸地高程插值圖Fig. 5 Interpolation map of land elevation

3 街區(qū)尺度風(fēng)暴潮模擬驗證

3.1 1211號臺風(fēng)“?？憋L(fēng)暴潮過程驗證

1211號強臺風(fēng)“?？庇?012年8月3日在西北太平洋上生成，5日17時加強為強熱帶風(fēng)暴，6日加強為臺風(fēng)，8日凌晨03時20分在浙江省象山縣鶴浦鎮(zhèn)登陸，登陸時中心氣壓為965 hPa，近中心風(fēng)力為14級。8日16時在浙江省杭州市境內(nèi)減弱為強熱帶風(fēng)暴。分別選取鎮(zhèn)海、吳淞、乍浦和定海4個驗潮站。臺風(fēng)風(fēng)場資料來自臺風(fēng)年鑒，臺風(fēng)路徑和驗潮站分布如圖6所示。

圖6 1211號臺風(fēng)“海葵”路徑（a）和驗潮站分布（b）Fig. 6 Track of No. 1211 Typhoon Haikui (a) and distribution of tide gauge stations (b)

采用臺風(fēng)與天文潮的非線性耦合，模擬1211號臺風(fēng)“海葵”風(fēng)暴潮過程。選取時間8月6日0時至8月8日20時，將4個驗潮站風(fēng)暴增水的實測值（藍色散點）與模擬值（紅色曲線）對比（圖7）。表1為4個驗潮站最大增水模擬值與實測值相對誤差分析。

圖7 鎮(zhèn)海、吳淞、乍浦和定海站風(fēng)暴潮實測值和模擬值對比Fig. 7 Comparison of measured and simulated storm surge values at Zhenhai,Wusong,Zhapu and Dinghai stations

表1 4個驗潮站最大增水實測值和模擬值相對誤差分析統(tǒng)計（單位：m）Table 1 Statistical table of relative error analysis of the measured and simulated values of the maximum water increasing at 4 tide gauge stations (unit: m)

由圖7可以看出，4個測站的增水趨勢大體一致，對比最大增水的發(fā)生時刻，鎮(zhèn)海站模擬與實測值的最大增水發(fā)生時刻基本一致，定海站模擬的最大增水推后1 h，乍浦站模擬值的最大增水時刻提前6 h，吳淞站提前7 h。結(jié)合表1來看，鎮(zhèn)海、吳淞、乍浦和定海4個驗潮站的相對誤差都較小，說明模型較好地重現(xiàn)了此次臺風(fēng)風(fēng)暴潮的發(fā)生過程。

3.2 5612號臺風(fēng)“溫黛”風(fēng)暴潮過程驗證

1956年5612號臺風(fēng)“溫黛”生成于菲律賓以東洋面，7月28日02時達到臺風(fēng)強度，30日14時臺風(fēng)強度達到最強（臺風(fēng)中心氣壓905 hPa，近中心最大風(fēng)速90 m/s）。8月1日由石垣島和沖繩之間進入東海，8月2日00時登陸浙江象山，臺風(fēng)中心最大風(fēng)力達到16級，登陸時中心氣壓達到923 hPa，登陸后繼續(xù)向西北方向移動，8月5日消失于陜西境內(nèi)。選取鎮(zhèn)海、吳淞、乍浦和高橋4個驗潮站，驗潮站都位于臺風(fēng)路徑右側(cè)。臺風(fēng)路徑和驗潮站分布如圖8所示。

圖8 5612號臺風(fēng)“溫黛”路徑（a）和驗潮站位置分布（b）Fig. 8 Track of No. 5612 Typhoon Wanda (a) and location distribution of tide gauge stations (b)

采用臺風(fēng)與天文潮的非線性耦合，模擬5612號臺風(fēng)“溫黛”風(fēng)暴潮過程（圖9）。將4個驗潮站的最大增水實測值和模擬值進行相對誤差分析，如表2所示。從圖9和表2可以看出，各驗潮站增水變化趨勢基本一致，增水峰值時間基本同步。最大增水時間對比，鎮(zhèn)海站提前2.5 h，乍浦站提前0.67 h，吳淞站推后1.75 h，高橋站推后1.17 h。結(jié)合臺風(fēng)路徑分析，鎮(zhèn)海、吳淞、乍浦和高橋站位于臺風(fēng)路徑右側(cè)，離臺風(fēng)登陸點位置較近，模擬值偏大一點，相對誤差分別為9.13%、12.96%、5.6%和5.44%，平均誤差為8.28%，模擬效果較好?？傮w來說，街區(qū)尺度風(fēng)暴潮漫灘模型能較準確模擬5612號臺風(fēng)“溫黛”風(fēng)暴潮情況。

圖9 鎮(zhèn)海、吳淞、乍浦和高橋站風(fēng)暴潮實測值與模擬值對比Fig. 9 Comparison of measured and simulated storm surge values at Zhenhai,Wusong,Zhapu and Gaoqiao stations

表2 4個驗潮站最大增水實測值與模擬值相對誤差分析統(tǒng)計表（單位：m）Table 2 Statistical table of relative error analysis between the measured and simulated values of the maximum water increasing at 4 tide gauge stations (unit: m)

4 街區(qū)尺度風(fēng)暴潮漫灘模擬

5612號臺風(fēng)“溫黛”登陸強度大，深入內(nèi)陸深，浙江遭受臺風(fēng)損失最為嚴重。全省共有75個縣（市）和不同程度受災(zāi)作物735萬畝。此次風(fēng)暴潮過程破壞了85萬間房屋，超過4 900人被砸死、淹沒或者觸電身亡，超過15 000人受傷[12]。產(chǎn)生的嚴重風(fēng)暴潮災(zāi)害影響是近幾十年來罕見的，因此選取5612號臺風(fēng)“溫黛”進行漫灘數(shù)值模擬。

4.1 徑流的設(shè)置

甬江位于浙江東部，4 518 km2的流域面積，橫貫寧波市區(qū)，在寧波市區(qū)東北部由奉化江、姚江兩條江匯合后入海，中下游主要為平原河流。由于姚江下游建成姚江大閘，如圖10所示，本文選取了澄浪堰和姚江大閘所處位置為姚江和奉化江集中入流的徑流邊界。

圖10 甬江澄浪堰和姚江大閘徑流邊界位置分布Fig. 10 Runoff boundary location distribution of Chenglang Weir on Yongjiang River and Yaojiang Gate

根據(jù)李文杰和邵學(xué)強[13]關(guān)于甬江徑流流量的資料整理，采用1983?2003年21年間年均尺度徑流流量作為姚江大閘和澄浪堰邊界的徑流輸入數(shù)據(jù)。澄浪堰和姚江大閘的年均流量分別為16.475×108m3和12.321×108m3，本文選用這兩個數(shù)據(jù)作為徑流邊界流量，并把徑流流量當(dāng)作為常數(shù)。

4.2 總潮位數(shù)值模擬

利用已建立的街區(qū)尺度風(fēng)暴潮漫灘模型模擬5612號臺風(fēng)“溫黛”影響下，北侖區(qū)附近驗潮站總潮位的情況，并以此來驗證5612號臺風(fēng)“溫黛”過程中，所建立的街區(qū)尺度風(fēng)暴潮漫灘模型對北侖區(qū)風(fēng)暴潮漫灘模擬的準確性。

將鎮(zhèn)海、吳淞、乍浦和高橋4個驗潮站總潮位實測值與模擬值進行對比（圖11）。由圖11可以看出，對于各站總潮位的模擬，趨勢與實測數(shù)據(jù)大體一致，最高潮位時間稍有偏差。將鎮(zhèn)海、吳淞、乍浦、高橋4個驗潮站模擬的最高潮位和實測最高潮位進行對比，發(fā)現(xiàn)各站絕對誤差為0.03 m、0.34 m、0.79 m、0.36 m，各站相對誤差為1.6%、14.5%、17.97%、15.28%。從模擬結(jié)果以及誤差分析看，鎮(zhèn)海站模擬的誤差較小，吳淞、乍浦和高橋的誤差稍大，是由于吳淞、乍浦和高橋3個站位于臺風(fēng)的右側(cè)，模擬的最高潮位相對于實測值偏大一點。因此從總潮位的趨勢看，建立的街區(qū)尺度模型可以較好地模擬北侖區(qū)的最高潮位。

圖11 鎮(zhèn)海、吳淞、乍浦和高橋站總潮位實測值與模擬值對比Fig. 11 Comparison of measured and simulated total tide level at Zhenhai,Wusong,Zhapu and Gaoqiao stations

4.3 漫灘數(shù)值模擬

在5612號臺風(fēng)“溫黛”影響期間，浙江沿海產(chǎn)生特大海潮，象山縣的最高潮位為4.7 m。據(jù)不完全統(tǒng)計，浙江、上海、安徽和河南省共影響作物6 946萬畝，摧毀220萬戶房屋。但是當(dāng)時沒有對漫灘淹沒范圍和水深進行系統(tǒng)的統(tǒng)計，缺乏當(dāng)年實測數(shù)據(jù)與模擬值對比。所以在風(fēng)暴增水和總潮位模擬較好的效果下，認為此漫灘模擬結(jié)果有一定合理性。

在不考慮海堤的情況下，利用已建立的街區(qū)尺度風(fēng)暴潮漫灘模型對5612號臺風(fēng)“溫黛”作用下北侖區(qū)的漫灘情況進行模擬，其結(jié)果如圖12所示。

圖12 5612號臺風(fēng)“溫黛”漫灘模擬的最大淹沒范圍Fig. 12 Maximum inundation area simulated by No. 5612 Typhoon Wanda

從圖12中可以看出，在5612號臺風(fēng)“溫黛”作用下，北侖區(qū)選定區(qū)域漫灘面積約為628.9 km2，淹沒了近岸大部分區(qū)域，淹沒水深大部分在0.1～0.4 m，最大水深在2.0 m左右。水流繞過高低層建筑，淹沒了建筑周圍的街道，更加清晰地看到淹沒區(qū)域的具體情形。在邊界處淹沒水深較大，是因為此處為陸地邊界，從而導(dǎo)致水流在這附近區(qū)域堆積。在街區(qū)尺度網(wǎng)格下，可以直觀清晰地看到水勢漫延趨勢，以及在高低層建筑物的阻擋下水勢的走向，體現(xiàn)出街區(qū)尺度網(wǎng)格的優(yōu)勢。

為了更直觀地體現(xiàn)5612號臺風(fēng)“溫黛”漫灘過程對北侖區(qū)街區(qū)和建筑的淹沒情況，利用GIS三維可視化工具，在數(shù)字地面模型基礎(chǔ)上疊加高分辨率遙感影像制作了此次漫灘過程的三維立體圖。選取時間是8月1日14時至8月2日10時，時間間隔為4 h。從圖13可以看出此次漫灘過程水深逐漸加深和降低的細致變化，對北侖區(qū)沿海街道和建筑的淹沒比以往的漫灘模擬更加精細?？梢钥闯鏊髟诮謪^(qū)和建筑錯綜復(fù)雜分布時的流動情況，以及不同建筑物附近水深的變化情況，體現(xiàn)了街區(qū)尺度漫灘模擬的顯著特點。

圖13 5612號臺風(fēng)“溫黛”漫灘過程三維圖Fig. 13 Three-dimensional diagram of the inundation process of No. 5612 Typhoon Wanda

5 下墊面底摩擦變化對風(fēng)暴潮漫灘的影響

不同下墊面具有不同的底摩擦系數(shù)，在海水的流動過程中影響著海水的流速和流向，這都直接影響漫灘模擬的結(jié)果[14–16]。

將下墊面的不同地物分類，修改控制方程中底摩擦系數(shù)的表達形式。此次采用二次形式的底摩擦，將下墊面底摩擦的變化考慮在ADCIRC模型中，底摩擦項為

式中，g為重力加速度；n為曼寧系數(shù)；H0為臨界水深，即H0=Hbreak。

利用北侖區(qū)遙感圖像對北侖區(qū)地物進行分類。根據(jù)美國NLCD數(shù)據(jù)的分類標準，我們選取21、22、23、31、43、51、83等 7類地物，不同地物的曼寧系數(shù)

如表3所示。即高、低密度居民區(qū)、商業(yè)區(qū)、裸露的巖石和沙地、四季常青和隨季節(jié)變化的兩種森林的混合、灌木叢和農(nóng)作物稀少地，圖14為所選區(qū)域曼寧系數(shù)n的分布示意圖。

圖14 北侖區(qū)曼寧系數(shù)n示意圖Fig. 14 Schematic diagram of Manning-n in the Beilun District

表3 7類地物的曼寧系數(shù)（n）設(shè)置Table 3 Manning-n settings for seven types of features

在考慮下墊面底摩擦變化的條件下，模擬北侖區(qū)在5612號臺風(fēng)“溫黛”作用下最大水深淹沒范圍，如圖15所示。淹沒面積達到494.37 km2，對比未考慮下墊面底摩擦情況的模擬實驗結(jié)果，淹沒面積減少了21.4%。沿岸大部分區(qū)域，最大淹沒水深在0.2 m左右，部分區(qū)域達到0.5 m，沿岸淹沒較深的局部區(qū)域可以達到1 m以上。

圖15 考慮下墊面底摩擦變化后，北侖區(qū)最大淹沒范圍Fig. 15 The maximum submerged area of Beilun District after considering the change of the underlying surface friction

圖16為對比了考慮下墊面底摩擦變化前后，北侖區(qū)淹沒水深的變化，即考慮下墊面變化后的最大淹沒水深減去未考慮下墊面變化時模擬的最大淹沒水深的差值?？梢钥闯鲈诳紤]下墊面底摩擦變化之后，大部分區(qū)域水深降低0.1～0.2 m。其中大部分水深降低區(qū)域集中在沿岸的低密度居民區(qū)，降低0.1 m左右。在農(nóng)作物、商業(yè)區(qū)局部區(qū)域，水深降低了0.5 m左右。其中水深降低最大的區(qū)域是在北侖區(qū)內(nèi)的兩條河流兩側(cè)。因此可以看出，在考慮下墊面的底摩擦變化之后，其對最大淹沒水深有降低作用，在不同的區(qū)域水深降低程度不同。

圖16 考慮下墊面底摩擦變化前后，北侖區(qū)最大淹沒水深差值Fig. 16 Maximum submerged depth difference in the Beilun District before and after considering the underlying surface friction

圖17為對比考慮底摩擦前后，漫灘面積的變化。藍色為面積減少區(qū)域，約減少138.76 km2，該區(qū)域多為低密度居民區(qū)。紅色為面積增加區(qū)域，約增加15.61 km2，該區(qū)域為巖石區(qū)和低密度居民區(qū)?？傮w而言，淹沒面積減少。極小的面積增加區(qū)域可能是由于在考慮下墊面底摩擦變化之后，街區(qū)細致刻畫使得水流流速變化，致使部分區(qū)域漫灘面積增加。因此在考慮底摩擦情況后，淹沒面積總體減少，下墊面的底摩擦變化對淹沒面積有減少作用。

圖17 考慮下墊面底摩擦變化前后，淹沒區(qū)域變化Fig. 17 The change of the submerged area before and after considering the change of the underlying surface friction

6 結(jié)論

本文采用5～10 m的街區(qū)尺度網(wǎng)格，依據(jù)寧波市北侖區(qū)高分辨率的街區(qū)尺度高程數(shù)據(jù)，細致刻畫街區(qū)道路和建筑分布，建立街區(qū)尺度風(fēng)暴潮漫灘模型，開展基于街區(qū)尺度的風(fēng)暴潮漫灘數(shù)值模擬，并考慮下墊面底摩擦的變化對漫灘的影響。結(jié)論如下：

（1）選取對北侖區(qū)影響嚴重的5612號臺風(fēng)“溫黛”，在考慮甬江徑流和不考慮堤壩的條件下，街區(qū)尺度風(fēng)暴潮漫灘模型可以較好地模擬風(fēng)暴潮和總潮位。依托街區(qū)尺度網(wǎng)格，可以看到水流沿著街道向城市內(nèi)部漫延的情況以及在高低層建筑物的阻擋下水勢的走向，體現(xiàn)出高分辨率街區(qū)尺度網(wǎng)格的優(yōu)勢。

（2）考慮下墊面底摩擦變化對漫灘的影響，漫灘淹沒面積達到494.37 km2，對比未考慮底摩擦情況的模擬結(jié)果，淹沒面積減少了21.4%。在街區(qū)尺度網(wǎng)格下，細致展示了淹沒水深和淹沒面積在不同地物分布情況下的變化。大部分水深降低區(qū)域集中在沿岸的低密度居民區(qū)，降低0.1 m左右。在農(nóng)作物、商業(yè)區(qū)局部區(qū)域，水深降低了0.5 m左右。其中水深降低最大的區(qū)域集中在北侖區(qū)內(nèi)的兩條河流兩側(cè)。淹沒面積減少的區(qū)域多為低密度居民區(qū)。總體而言，考慮下墊面底摩擦變化后，淹沒水深降低，淹沒面積減少。