戚衛(wèi)偉

(遼寧省營(yíng)口水文局，遼寧 營(yíng)口 115000)

1 概 況

丹東市位于遼寧省東南部，與朝鮮隔江相望，屬于遼東山地丘陵的一部分，地勢(shì)總體北高南低，總體地形以山地和丘陵為主，局部還存在階地、盆地、臺(tái)地等小型地貌單元，地形結(jié)構(gòu)非常復(fù)雜。丹東市屬于溫帶亞濕潤(rùn)季風(fēng)型氣候，年平均降雨量可達(dá)800～1200 mm，且2/3集中在夏季。特殊的地形條件加上較多的降水量，導(dǎo)致丹東地區(qū)特別容易發(fā)生汛期洪水災(zāi)害，因此建立洪水監(jiān)測(cè)信息系統(tǒng)很有必要。

2 研究方法

2.1 降雨監(jiān)測(cè)方法

(1)本項(xiàng)目相關(guān)降水?dāng)?shù)據(jù)來(lái)源于丹東市氣象局，共提取2000—2019年這20 a間的降水?dāng)?shù)據(jù)，以各水文站點(diǎn)平均測(cè)得數(shù)據(jù)為基準(zhǔn)；

(2)本項(xiàng)目相關(guān)地形數(shù)據(jù)來(lái)源于地理空間數(shù)據(jù)云，分辨率為25 m，以這些數(shù)據(jù)為基礎(chǔ)，利用GIS技術(shù)建立丹東市DEM模型；

(3)最后利用GIS系統(tǒng)中的空間分析模塊功能來(lái)提取各區(qū)域的坡度、坡向數(shù)據(jù)，河道分布統(tǒng)計(jì)為全市范圍內(nèi)所有5級(jí)以上河流。根據(jù)GIS雷達(dá)回波情況估計(jì)各區(qū)域降水量，對(duì)可能遭受洪水災(zāi)害區(qū)域提前發(fā)出預(yù)警信息[1]。

2.2 影響因子的Person和RDA分析方法

Person相關(guān)性分析法是用于度量?jī)蓚€(gè)變量X和Y之間線性相關(guān)程度的方法，其值介于-1～1之間，其中“1”表示變量完全正相關(guān)，“0”表示無(wú)關(guān)，“-1”表示完全負(fù)相關(guān)[2]。

RDA分析(Redundancy analysis)即冗余分析，該方法本質(zhì)是約束化的主成分分析，目的是尋找新的變量來(lái)代替原來(lái)變量。由于汛期洪水主要來(lái)源為降水量，所以本項(xiàng)目以降水量為目標(biāo)變量，以各類(lèi)環(huán)境因子為解釋變量，定量分析不同的環(huán)境因子對(duì)降水量的影響。

通過(guò)翻閱相關(guān)文獻(xiàn)可知：影響丹東市不同區(qū)域降水量的三個(gè)因子為：地形(海拔、坡度、坡向等)、下墊面(植被、河道分布)、海陸位置(距離、經(jīng)度、緯度)，其中坡向指數(shù)ASP按照(1)式計(jì)算[3]，該值越大，表明向陽(yáng)面積越大。

(1)

式中：a為坡向?qū)崪y(cè)值，(°)。

3 丹東市降水分布特征及影響因子冗余分析

3.1 汛期丹東市GIS雷達(dá)波回波圖分析

GIS雷達(dá)預(yù)警系統(tǒng)可以充分結(jié)合區(qū)域地理信息特征將回波信息準(zhǔn)確覆蓋到具體鄉(xiāng)鎮(zhèn)水平，而且可實(shí)現(xiàn)對(duì)暴雨、臺(tái)風(fēng)等災(zāi)害天氣預(yù)測(cè)預(yù)警能力。

3.1.1 水體識(shí)別原理

本項(xiàng)目應(yīng)用MODIS數(shù)據(jù)衛(wèi)星接收區(qū)域影像，之后通過(guò)遙感模塊處理來(lái)獲取并傳輸至用戶GIS系統(tǒng)，因此遙感數(shù)據(jù)首先需要能夠準(zhǔn)確識(shí)別水體。在MODIS數(shù)據(jù)衛(wèi)星中，第1波段紅光區(qū)(0.62～0.67 μm)下，水體比地表反射率明顯偏高；在第2波段紅外區(qū)(0.84～0.87 μm)，水體比地表反射率明顯偏低。根據(jù)兩個(gè)波段下地表和水體的反射率不同，利用歸一化差異植被指數(shù)NDVI來(lái)增強(qiáng)兩者差別，計(jì)算見(jiàn)式(2)[4]。

(2)

式中：CH1為第1波段下地表反射率；CH2為第2波段下地表反射率。

3.1.2 GIS雷達(dá)回波與降雨強(qiáng)度關(guān)系

遙感影像僅僅分別出地表和水體區(qū)域還遠(yuǎn)遠(yuǎn)不夠，各區(qū)域降雨強(qiáng)度才是預(yù)警汛期洪水的核心?？萍既藛T通過(guò)對(duì)比大量GIS雷達(dá)回波強(qiáng)度Z和降雨強(qiáng)度I關(guān)系，得出兩者之間基本符合如(3)所示關(guān)系，利用誤差分析法來(lái)確定A、b的值[5]。

(3)

式中：A、b為相關(guān)參數(shù)，與降雨類(lèi)型、季節(jié)、地形相關(guān)。

圖1是2019年9月8日拍攝的遼寧地區(qū)GIS雷達(dá)降雨回波圖，由圖可知：丹東地區(qū)的西部、中部、南部均被大到暴雨覆蓋。之后通過(guò)進(jìn)一步提高地圖精度，對(duì)可能存在洪水危險(xiǎn)的居民區(qū)提前12 h發(fā)出預(yù)警，并對(duì)可能被洪水淹沒(méi)地區(qū)的居民分批有序撤離至安置點(diǎn)。

圖1 2019年9月8日遼寧地區(qū)GIS雷達(dá)回波圖(顏色越深，降雨強(qiáng)度越大)

3.1.3 丹東市平均降水分布特征分析

通過(guò)對(duì)2010—2019年間丹東地區(qū)平均降水量分布圖分析，可得出以下結(jié)論：(1)總體來(lái)看，以東北—西南條狀方向?yàn)橛^察線，降水量從東部到西部依次減少。(2)降水量最大的地區(qū)為丹東市區(qū)，其次為東港市一小部分，中部相對(duì)于東、西部分降雨量要大。(3)通過(guò)對(duì)比平均降雨量分布圖和GIS雷達(dá)回波圖，確認(rèn)了兩者一致性達(dá)到80%，也直接驗(yàn)證了GIS雷達(dá)回波圖的準(zhǔn)確性[6]。

3.2 影響丹東地區(qū)降雨量的因子計(jì)算分析3.2.1 Person相關(guān)性分析

在此首先應(yīng)用Person相關(guān)性分析法來(lái)研究地形(海拔、坡度、坡向等)、下墊面(植被、河道分布)、海陸位置(距離、經(jīng)度、緯度等)與降雨量大小之間的關(guān)系，部分分析結(jié)果如下：

(1)降雨量與坡向關(guān)系

為簡(jiǎn)化計(jì)算，本項(xiàng)目用“4方向法”將坡向分為：陽(yáng)坡、半陽(yáng)坡、半陰坡、陰坡(如圖2所示)，由于丹東市屬于溫帶季風(fēng)氣候，所以陽(yáng)坡均為迎風(fēng)破，陰坡均為背風(fēng)坡。由表1數(shù)據(jù)可知：隨著陽(yáng)坡向陰坡轉(zhuǎn)化，年平均降雨量總體呈減少趨勢(shì)[7]。

表1 丹東地區(qū)降雨量與坡向數(shù)據(jù)(部分)

圖2 “4方向法”劃分坡向示意圖

本項(xiàng)目經(jīng)過(guò)Person相關(guān)性分析法來(lái)分析地形因素，得出降雨量與海拔、坡向、坡度的相關(guān)性系數(shù)分別為：0.62、0.51、0.22。

(2)降雨量與海洋間距關(guān)系

由丹東市地圖可知：丹東市市區(qū)南部緊鄰黃海，分析降水分布圖可知整個(gè)丹東地區(qū)由南向北，降雨量是逐漸減少趨勢(shì)。

本項(xiàng)目經(jīng)過(guò)Person相關(guān)性分析法來(lái)分析海陸位置因素，得出降雨量與海洋間距、緯度、經(jīng)度的相關(guān)性系數(shù)分別為：0.78、0.21、0.01。

3.2.2 RDA冗余分析

經(jīng)過(guò)對(duì)影響因子進(jìn)行RDA排序，得出三大主要影響因素獨(dú)立解釋信息具體數(shù)值，如圖3所示。其中海陸位置解釋了0.224的信息，在所有值中是最大值，因此可知該因素是影響丹東地區(qū)降水量的主因，其次為地形因素、下墊因素。

圖3 RDA冗余分析影響因子結(jié)果

根據(jù)RDA分析結(jié)果，在無(wú)法改變海陸位置、地形因素下，只能通過(guò)人為改變下墊因素來(lái)調(diào)節(jié)降雨量。包括：封山育林、打造海綿城鎮(zhèn)，一方面通過(guò)森林蒸騰作用阻滯雨水，另一方面加強(qiáng)居民區(qū)吸收積水能力。此外，針對(duì)丹東市區(qū)雨量較大問(wèn)題，只能通過(guò)改善城市排水系統(tǒng)應(yīng)對(duì)。

4 結(jié) 語(yǔ)

汛期洪水預(yù)警信息化建設(shè)是當(dāng)前信息化和現(xiàn)代化背景的要求，同時(shí)也是利用新型計(jì)算機(jī)手段維護(hù)居民生命財(cái)產(chǎn)安全的重要途徑。丹東地區(qū)通過(guò)這一系列操作，大大提高了對(duì)洪水預(yù)警水平，與之前相比，可使洪水經(jīng)濟(jì)損失數(shù)額減少50%以上，取得了很好的社會(huì)及經(jīng)濟(jì)效益。并且對(duì)以后城市建設(shè)、突發(fā)性事件管理積累了寶貴經(jīng)驗(yàn)。