面向?yàn)?zāi)害風(fēng)險(xiǎn)評(píng)估的臺(tái)風(fēng)數(shù)值模擬及可視化

引文格式： 王中山，袁金錦. 面向?yàn)?zāi)害風(fēng)險(xiǎn)評(píng)估的臺(tái)風(fēng)數(shù)值模擬及可視化[J].測(cè)繪通報(bào)，2015(4)：108-110.DOI:10.13474/j.cnki.11-2246.2015.0123

面向?yàn)?zāi)害風(fēng)險(xiǎn)評(píng)估的臺(tái)風(fēng)數(shù)值模擬及可視化

王中山1，袁金錦2

(1. 武漢大學(xué)測(cè)繪學(xué)院，湖北 武漢 430079； 2. 山東省地質(zhì)測(cè)繪院，山東 濟(jì)南 250000)

Typhoon Numerical Simulation and Visualization for Disaster

Risk Assessment

WANG Zhongshan，YUAN Jinjin

摘要：針對(duì)具體臺(tái)風(fēng)對(duì)地區(qū)災(zāi)害風(fēng)險(xiǎn)影響的快速應(yīng)急的特點(diǎn)，運(yùn)用Shapiro臺(tái)風(fēng)風(fēng)場(chǎng)模型，結(jié)合GIS技術(shù)，以201319號(hào)“天兔”臺(tái)風(fēng)為例，實(shí)現(xiàn)了臺(tái)風(fēng)的數(shù)值模擬、該臺(tái)風(fēng)具體時(shí)刻的影響區(qū)域的可視化及對(duì)廣東省的災(zāi)害風(fēng)等級(jí)區(qū)劃，以期為應(yīng)急部門的防災(zāi)減災(zāi)工作提供科學(xué)的參考依據(jù)。

關(guān)鍵詞：臺(tái)風(fēng)；GIS；數(shù)值模擬；可視化

中圖分類號(hào)：P208

收稿日期：2014-09-22

基金項(xiàng)目：城市公共安全脆弱性分析和綜合風(fēng)險(xiǎn)評(píng)估關(guān)鍵技術(shù)研究與示范(2011BAK07B02)

作者簡(jiǎn)介：王中山(1989—)，男，碩士，主要研究方向?yàn)榈乩硇畔⑾到y(tǒng)理論與應(yīng)用。E-mail：wangzhongshan89@163.com

一、引言

臺(tái)風(fēng)是一種破壞力很強(qiáng)的災(zāi)害性天氣，登陸地一般集中在我國(guó)沿海的浙江、福建和廣東等經(jīng)濟(jì)較發(fā)達(dá)的地區(qū)，其伴隨有大風(fēng)、暴雨，以及其引發(fā)的風(fēng)暴潮、滑坡、泥石流和洪水等一系列災(zāi)害，給人們的生命財(cái)產(chǎn)造成巨大的損失。面對(duì)臺(tái)風(fēng)災(zāi)害對(duì)社會(huì)和經(jīng)濟(jì)的嚴(yán)重影響，進(jìn)行臺(tái)風(fēng)風(fēng)場(chǎng)研究，進(jìn)而模擬臺(tái)風(fēng)，促進(jìn)人們深刻認(rèn)識(shí)臺(tái)風(fēng)和制定有效的措施來降低其造成的人員傷亡和財(cái)產(chǎn)損失，具有很強(qiáng)的必要性和現(xiàn)實(shí)意義。

目前，對(duì)臺(tái)風(fēng)風(fēng)場(chǎng)模型的研究，大多數(shù)的研究重點(diǎn)是臺(tái)風(fēng)風(fēng)場(chǎng)模型的數(shù)值模擬及如何提高模擬精度[1-2]。此外，近年來不少學(xué)者也開展了臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)區(qū)劃的相關(guān)研究，研究出災(zāi)害風(fēng)險(xiǎn)的大小及災(zāi)害可能帶來的損失，并且對(duì)等級(jí)劃分，區(qū)分出不同等級(jí)風(fēng)險(xiǎn)的區(qū)域范圍[3-4]。

本文結(jié)合公共安全中應(yīng)急處置快速反應(yīng)的特點(diǎn)，在臺(tái)風(fēng)風(fēng)場(chǎng)數(shù)值模擬模型的基礎(chǔ)上，運(yùn)用GIS空間數(shù)據(jù)管理、GIS空間數(shù)據(jù)疊加分析、GIS空間數(shù)據(jù)插值等技術(shù)，實(shí)現(xiàn)臺(tái)風(fēng)動(dòng)態(tài)時(shí)刻影響區(qū)域及整場(chǎng)臺(tái)風(fēng)對(duì)某個(gè)地區(qū)的風(fēng)險(xiǎn)區(qū)劃，確定風(fēng)險(xiǎn)較高的區(qū)域和較低的區(qū)域，為政府決策部門的資源配置、防災(zāi)減災(zāi)和應(yīng)急決策提供科學(xué)的參考和依據(jù)，從而進(jìn)行必要的應(yīng)急疏散和救援調(diào)度。

二、臺(tái)風(fēng)風(fēng)場(chǎng)模型及其數(shù)值模擬

1. 臺(tái)風(fēng)風(fēng)場(chǎng)模型

本文選取Shapiro風(fēng)場(chǎng)數(shù)值模擬模型來對(duì)臺(tái)風(fēng)進(jìn)行模擬，Shapiro風(fēng)場(chǎng)模型是一組描述大氣運(yùn)動(dòng)的偏微分方程，其理論模型是在不考慮熱帶氣旋豎向?qū)α鞯募僭O(shè)條件下，氣壓場(chǎng)對(duì)稱、定常厚度的邊界層中，通過牛頓第二定律經(jīng)過坐標(biāo)系轉(zhuǎn)換得到的。Shapiro風(fēng)場(chǎng)模型在移動(dòng)柱坐標(biāo)中，動(dòng)量守恒控制方程為[5-6]

(1)

(2)

Shapiro風(fēng)場(chǎng)是一組非線性偏微分方程組，本文采用有限差分法進(jìn)行求解[10]。為了使差分?jǐn)?shù)值具有較好的穩(wěn)定性，所有的空間一階導(dǎo)數(shù)項(xiàng)采用迎風(fēng)格式，空間二階導(dǎo)數(shù)項(xiàng)采用中心差分格式，時(shí)間偏導(dǎo)數(shù)項(xiàng)采用向前差分格式，應(yīng)用時(shí)間推進(jìn)法和牛頓迭代法相結(jié)合的方式求解離散化方程。

迭代初始值由梯度風(fēng)場(chǎng)作為初始值，邊界條件為在臺(tái)風(fēng)中心處，風(fēng)速值為零；最外圍的網(wǎng)格節(jié)點(diǎn)是忽略流體粘性加速度項(xiàng)和阻力加速度項(xiàng)穩(wěn)定狀態(tài)時(shí)的風(fēng)速，也就是科氏力和氣壓梯度力平衡，即

(3)

2. 數(shù)值模擬

基于上述的Shapiro風(fēng)場(chǎng)理論模型，首先根據(jù)臺(tái)風(fēng)的觀測(cè)數(shù)據(jù)得到B、Rmax、ΔP，然后對(duì)臺(tái)風(fēng)進(jìn)行網(wǎng)格劃分[10](如圖1所示)，根據(jù)邊界條件求出初始值(u0,v0)；將其和得到的臺(tái)風(fēng)關(guān)鍵參數(shù)一起代入由式(1)得到的差分控制方程，迭代求解得到臺(tái)風(fēng)風(fēng)場(chǎng)中網(wǎng)格點(diǎn)的臺(tái)風(fēng)風(fēng)速(u,v)。Shapiro風(fēng)場(chǎng)計(jì)算得到的風(fēng)速值為邊界層高度內(nèi)豎向平均值，近似為500m高空中不受地面摩擦影響的梯度風(fēng)速，最后將計(jì)算出的風(fēng)速值采用式(4)轉(zhuǎn)換為近地面10m高度的風(fēng)速值。

圖1　臺(tái)風(fēng)影響范圍內(nèi)的節(jié)點(diǎn)

(4)

式中，v為離地面高度為z的平均風(fēng)速(m/s)；v1為離地面高度為z1的平均風(fēng)速值(m/s)；α為地面粗糙度指數(shù)。

本文選取Shapiro臺(tái)風(fēng)風(fēng)場(chǎng)模型對(duì)201319號(hào)臺(tái)風(fēng)“天兔(Usagi)”進(jìn)行模擬。“天兔”起始于2013年9月15日菲律賓以東；18日5時(shí)加強(qiáng)為強(qiáng)熱帶風(fēng)暴，20時(shí)加強(qiáng)為臺(tái)風(fēng)；19日11時(shí)加強(qiáng)為強(qiáng)臺(tái)風(fēng)，17時(shí)加強(qiáng)為超強(qiáng)臺(tái)風(fēng)；22日19時(shí)40分，“天兔”帶著國(guó)家氣象中心(NMC)的45m/s的評(píng)價(jià)，登上了汕尾市的陸地，汕尾市持續(xù)風(fēng)速達(dá)15級(jí)，陸豐測(cè)得極大風(fēng)速74.7m/s，超過17級(jí)。

為了更好地對(duì)觀測(cè)值與模擬值進(jìn)行比較，選取了臺(tái)風(fēng)路徑附近兩個(gè)觀測(cè)站點(diǎn)，分別是東沙島(21.0213°N, 118.8611°E)、陸豐南塘(22.9093°N, 115.9131°E)。利用Shapiro臺(tái)風(fēng)風(fēng)場(chǎng)模型計(jì)算出這些站點(diǎn)在臺(tái)風(fēng)影響期間的風(fēng)速值，與臺(tái)風(fēng)中心靠近觀測(cè)站點(diǎn)的時(shí)段觀測(cè)到的風(fēng)速值進(jìn)行比較。

如圖2和圖3所示，從中可以看出，東沙島觀測(cè)站只有在1∶00這個(gè)時(shí)刻模擬風(fēng)速值與實(shí)際觀測(cè)值差距較大，達(dá)到11m/s，其他時(shí)刻差異較小，基本吻合，整體變化趨勢(shì)也相同；陸豐南塘站在臺(tái)風(fēng)登陸期間的幾個(gè)時(shí)刻內(nèi)，只有20∶00時(shí)觀測(cè)值與模擬風(fēng)速值差距最大，為12m/s。造成這種個(gè)別點(diǎn)觀測(cè)值與模擬值差距較大的主要原因可能是臺(tái)風(fēng)本身具有的復(fù)雜性及Shapiro風(fēng)場(chǎng)模型在某些地方的參數(shù)的簡(jiǎn)化。總體上，Shapiro風(fēng)場(chǎng)模型能夠滿足應(yīng)急處置快速反應(yīng)的模擬效果。

圖2　東沙島

圖3　陸豐南塘站

三、GIS平臺(tái)的可視化

利用GIS相關(guān)技術(shù)和手段將Shapiro風(fēng)場(chǎng)模型模擬結(jié)果以圖形的形式直觀地顯示出來，可以幫助應(yīng)急指揮決策人員從整體上把握整場(chǎng)臺(tái)風(fēng)對(duì)具體地區(qū)造成的災(zāi)害風(fēng)險(xiǎn)等級(jí)，也可以清楚地了解具體時(shí)刻臺(tái)風(fēng)的影響區(qū)域，以確定科學(xué)合理的防災(zāi)減災(zāi)措施，將臺(tái)風(fēng)災(zāi)害帶來的損失降到最低。本文以“天兔”對(duì)廣東省的影響為例，利用反距離權(quán)重的差值方法，生成“天兔”臺(tái)風(fēng)具體時(shí)刻影響區(qū)域圖，可以直觀地展現(xiàn)不同時(shí)刻“天兔”臺(tái)風(fēng)對(duì)廣東省的影響區(qū)域；以縣為基本單位，生成“天兔”臺(tái)風(fēng)對(duì)廣東省災(zāi)害風(fēng)險(xiǎn)區(qū)劃圖。具體的實(shí)現(xiàn)方法與步驟為：

1) 將Shapiro風(fēng)場(chǎng)模型計(jì)算出的風(fēng)速點(diǎn)值數(shù)據(jù)按照格式生成TXT文件。

2) 將TXT文件中的數(shù)據(jù)導(dǎo)入到ArcMap中，處理成為SHP文件，并對(duì)其進(jìn)行投影變換，使其與系統(tǒng)中的基礎(chǔ)地理信息數(shù)據(jù)具有相同的投影系統(tǒng)。

3) 使用Spatial Analyst 工具，將風(fēng)速值點(diǎn)文件數(shù)據(jù)進(jìn)行插值運(yùn)算，得到風(fēng)場(chǎng)風(fēng)速值的柵格文件，本文采用的差值方法為反距離權(quán)重法。將風(fēng)速值劃分為5個(gè)等級(jí)，對(duì)柵格文件進(jìn)行分類渲染。

4) 將渲染后的風(fēng)速值柵格文件與廣東省的基礎(chǔ)地理數(shù)據(jù)進(jìn)行疊加展示，生成具體某個(gè)時(shí)刻的臺(tái)風(fēng)影響區(qū)域?qū)ｎ}圖。如圖4和圖5分別為2013-09-22 3∶00和2013-09-22 14∶00“天兔”臺(tái)風(fēng)的影響區(qū)域圖。

圖4　2013-09-22 3：00影響區(qū)域圖

圖5　2013-09-22 14:00 影響區(qū)域圖

5) 使用ArctoolBox中的分析工具，將風(fēng)速值點(diǎn)文件數(shù)據(jù)與廣東省基礎(chǔ)地理數(shù)據(jù)進(jìn)行疊加，以縣域?yàn)榛締挝?，獲得該縣域的整場(chǎng)臺(tái)風(fēng)的最大風(fēng)速值，然后以臺(tái)風(fēng)最大風(fēng)速值作為臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)的評(píng)價(jià)指標(biāo)，分級(jí)顯示，生成廣東省災(zāi)害風(fēng)險(xiǎn)區(qū)劃圖，如圖6所示。

圖6　廣東省“天兔”臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)區(qū)劃圖

由圖6可知，“天兔”臺(tái)風(fēng)對(duì)廣東省的影響的災(zāi)害高風(fēng)險(xiǎn)區(qū)主要分布于廣東省東南部(汕尾市、陸豐市、潮陽市、東莞市、惠州市轄區(qū)等地區(qū))，這些地區(qū)位于臺(tái)風(fēng)路徑的兩側(cè)，受臺(tái)風(fēng)影響較大；隨著臺(tái)風(fēng)逐漸向內(nèi)陸移動(dòng)，災(zāi)害風(fēng)險(xiǎn)開始逐漸降低，次高風(fēng)險(xiǎn)區(qū)主要分布于廣東省中部(廣州市、增城市、英德市等地區(qū))；中等風(fēng)險(xiǎn)區(qū)和次低等風(fēng)險(xiǎn)區(qū)主要分布于廣東省中北部、西中部(連州市、陽春市、臺(tái)山市等地區(qū))；風(fēng)險(xiǎn)區(qū)主要分布廣東省西南部(高州市、雷州市、化州市等地區(qū))，這些地區(qū)離臺(tái)風(fēng)路徑較遠(yuǎn)，因此受臺(tái)風(fēng)影響較小。最后將“天兔”對(duì)廣東省影響的可視化成果與廣東省地區(qū)實(shí)際出現(xiàn)的風(fēng)速情況比較，發(fā)現(xiàn)基本吻合。

四、結(jié)束語

本研究基于Shapiro風(fēng)場(chǎng)模型，以201319號(hào)臺(tái)風(fēng)“天兔”為例，借助于GIS在空間數(shù)據(jù)存儲(chǔ)、空間數(shù)據(jù)分析、圖形交互和顯示方面的強(qiáng)大功能，將風(fēng)場(chǎng)模型與GIS結(jié)合，得出了天兔臺(tái)風(fēng)對(duì)廣東省災(zāi)害風(fēng)險(xiǎn)區(qū)劃圖和“天兔”具體時(shí)刻影響區(qū)域圖，為應(yīng)急部門更加直觀地了解相關(guān)情況和決策提供了科學(xué)的參考依據(jù)。

本文只是從臺(tái)風(fēng)最大風(fēng)速值研究了臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)，臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)評(píng)估和區(qū)劃涉及多學(xué)科理論、方法和手段，僅僅一個(gè)因素不能夠較全面地反映實(shí)際情況，同時(shí)還需要考慮承災(zāi)體的空間分布及承災(zāi)能力等其他因素。因此，借助于GIS的強(qiáng)大技術(shù)支持，結(jié)合臺(tái)風(fēng)數(shù)值模擬模型，不斷優(yōu)化和完善臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)評(píng)估的技術(shù)方法，提高臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)評(píng)估和區(qū)劃的精度，能夠更好地為防災(zāi)減災(zāi)事業(yè)服務(wù)，同時(shí)也具有重要的現(xiàn)實(shí)意義。

參考文獻(xiàn)：

[1]江志輝,華鋒,曲平，等.一個(gè)新的熱帶氣旋參數(shù)調(diào)整方案[J].海洋科學(xué)進(jìn)展,2008,26(1):1-7.

[2]林偉,方偉華.西北太平洋臺(tái)風(fēng)風(fēng)場(chǎng)模型中Holland B系數(shù)區(qū)域特征研究[J]. 熱帶地理, 2013(2):124-132.

[3]王勝,田紅,謝五三，等.基于GIS技術(shù)的臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)區(qū)劃研究——以安徽省為例[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào),2012,17(1):161-166.

[4]劉少軍,張京紅,何政偉，等.基于GIS的臺(tái)風(fēng)災(zāi)害損失評(píng)估模型研究[J].災(zāi)害學(xué),2010,25(2):64-67.

[5]李茜,段忠東.Shapiro臺(tái)風(fēng)風(fēng)場(chǎng)模型及其數(shù)值模擬[J].自然科學(xué)學(xué)報(bào),2005,14(1):45-52.

[6]SHAPIRO L J. The Asymmetric Boundary Layer Flow under a Translating Hurricane [J]. Journal of Atmospheric Sciences, 1983, 40(8): 1984-1998.

[7]ROLL H V. Physics of the Marine Atmosphere[M].[S.l.]：Academic Press, 1965:42.

[8]肖玉鳳,肖儀清,段忠東，等.Holland參數(shù)B對(duì)CE風(fēng)場(chǎng)數(shù)值模擬結(jié)果的影響研究[C]∥第十四屆全國(guó)結(jié)構(gòu)風(fēng)工程學(xué)術(shù)會(huì)議論文集.[北京]:[s.n.]，2009：163-168.

[9]肖玉鳳,段忠東,肖儀清，等. 基于數(shù)值模擬的臺(tái)風(fēng)危險(xiǎn)性分析綜述(Ⅰ)——臺(tái)風(fēng)風(fēng)場(chǎng)模型[J].自然災(zāi)害學(xué)報(bào), 2011, 20(2): 82-89.

[10]李茜.我國(guó)臺(tái)風(fēng)區(qū)城市的臺(tái)風(fēng)關(guān)鍵參數(shù)統(tǒng)計(jì)分析與Shapiro風(fēng)場(chǎng)模擬[D].哈爾濱:哈爾濱工業(yè)大學(xué), 2003.

[11]鄭衛(wèi)江,吳煥萍,羅兵，等.GIS技術(shù)在臺(tái)風(fēng)預(yù)報(bào)服務(wù)產(chǎn)品制作系統(tǒng)中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào),2010,21(2):250-255.

[12]吳迪軍,孫海燕,黃全義，等.河道洪水演進(jìn)數(shù)值模擬及其在GIS平臺(tái)上的可視化[J].測(cè)繪科學(xué),2009,34(6):242-243,189.

[13]劉順隆,鄭群.計(jì)算流體力學(xué)[M].哈爾濱：哈爾濱工程大學(xué)出版社,1998.

[14]陸金甫,關(guān)治.偏微分方程的數(shù)值解法[M].北京：清華大學(xué)出版社,2004.

[15]陳香.福建省臺(tái)風(fēng)災(zāi)害風(fēng)險(xiǎn)評(píng)估與區(qū)劃[J].生態(tài)學(xué)雜志,2007,26(6):961-966.

[16]李瑞龍.基于改進(jìn)的臺(tái)風(fēng)關(guān)鍵參數(shù)的臺(tái)風(fēng)極值風(fēng)速預(yù)測(cè)[D].哈爾濱：哈爾濱工業(yè)大學(xué),2007.

[17]王志春,宋麗莉,黃飛，等.沿海地區(qū)風(fēng)速指數(shù)的計(jì)算[J].廣東氣象,2006(4):17-18.

[18]劉仁義,劉南.ArcGIS開發(fā)寶典——從入門到精通[M].北京:科學(xué)出版社,2006.

徠卡測(cè)量新技術(shù)應(yīng)用專欄