孟 迎

(濟(jì)寧市水利事業(yè)發(fā)展中心,山東 濟(jì)寧 272000)

山洪災(zāi)害是一種由降雨或周期性流水引起的特殊自然災(zāi)害,不僅會(huì)導(dǎo)致地表徑流急劇起落,還會(huì)破壞周?chē)霓r(nóng)田、耕地、房屋[1],造成人員傷亡和財(cái)產(chǎn)損失。山洪災(zāi)害往往發(fā)生在山區(qū),而我國(guó)的山區(qū)面積較大[2],在夏季降雨頻發(fā)時(shí)期需要格外注意。為保障人民群眾生命財(cái)產(chǎn)安全,需要設(shè)計(jì)一種有效的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法。

早在2006年,我國(guó)就制定了與山洪災(zāi)害防治相關(guān)的政策[3],在全國(guó)范圍開(kāi)展了大規(guī)模的山洪防治工作,針對(duì)不同區(qū)域的實(shí)際地質(zhì)情況進(jìn)行了分析,初步構(gòu)建了預(yù)警體系。但山洪災(zāi)害的影響因素較多,難以識(shí)別動(dòng)態(tài)變化預(yù)警因素,整體預(yù)警防御能力偏低[4]。盡管目前我國(guó)已經(jīng)有了若干個(gè)監(jiān)測(cè)站點(diǎn),但大多數(shù)監(jiān)測(cè)站點(diǎn)分布稀疏,難以獲取全面的山洪預(yù)警數(shù)據(jù)[5]。相關(guān)研究人員針對(duì)山洪災(zāi)害的特點(diǎn)設(shè)計(jì)了若干種常見(jiàn)的山洪災(zāi)害監(jiān)測(cè)平臺(tái),主要通過(guò)一體化采集硬件采集環(huán)境變化數(shù)據(jù),實(shí)現(xiàn)山洪災(zāi)害預(yù)警[6],受采集硬件敏感性影響,基于現(xiàn)有監(jiān)測(cè)平臺(tái)的動(dòng)態(tài)預(yù)警效果相對(duì)一般。為提高山洪災(zāi)害動(dòng)態(tài)預(yù)警的可靠性,本文考慮降雨數(shù)值變化,設(shè)計(jì)了一種全新的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法。

1 考慮降雨量數(shù)值變化的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法設(shè)計(jì)

1.1 考慮降雨量數(shù)值變化修正山洪災(zāi)害預(yù)警數(shù)據(jù)

降雨量數(shù)值變化與山洪災(zāi)害的動(dòng)態(tài)耦合度正相關(guān),因此,要想提高山洪災(zāi)害動(dòng)態(tài)預(yù)警精度,需要基于降雨量數(shù)值修正山洪災(zāi)害預(yù)警數(shù)據(jù)[7]。首先通過(guò)水文氣象站的觀(guān)測(cè)記錄或相關(guān)機(jī)構(gòu)的網(wǎng)站獲取點(diǎn)降雨量數(shù)據(jù);然后使用TIGGE降雨量數(shù)值變化修正平臺(tái)修正預(yù)警數(shù)據(jù),將其轉(zhuǎn)化為有效的預(yù)警格式。經(jīng)過(guò)山洪災(zāi)害動(dòng)態(tài)預(yù)警數(shù)據(jù)初次采集發(fā)現(xiàn),不同流域、不同采集點(diǎn)采集的數(shù)據(jù)存在一定的重合關(guān)系[8],因此,本文使用降雨量預(yù)報(bào)轉(zhuǎn)化值反映降雨量數(shù)值變化,此時(shí)的降雨量動(dòng)態(tài)計(jì)算公式為

(1)

式中:ω為各個(gè)降雨量采集點(diǎn)的權(quán)重;pi為采集點(diǎn)所在位置的降雨量,mm。此時(shí)利用反距離權(quán)重法統(tǒng)計(jì)信息誤差值。

《降水量等級(jí)》(GB/T 28592—2012)中規(guī)定了降雨量的分級(jí)標(biāo)準(zhǔn),見(jiàn)表1。

表1 不同時(shí)段的降雨量等級(jí)劃分 單位:mm

根據(jù)表1的降雨量分級(jí)標(biāo)準(zhǔn)可以有效地統(tǒng)計(jì)信息誤差值,從而進(jìn)行山洪災(zāi)害預(yù)警數(shù)據(jù)修正,修正計(jì)算式如下:

(2)

式中:TS為預(yù)警數(shù)據(jù)修正誤差;NA為命中個(gè)數(shù);NB為空?qǐng)?bào)個(gè)數(shù);NC為漏報(bào)個(gè)數(shù)。

為降低最終的災(zāi)害預(yù)警數(shù)據(jù)修正誤差,本文設(shè)計(jì)的方法使用分位數(shù)映射法計(jì)算累計(jì)概率密度值,計(jì)算式pdf(x)如下:

(3)

式中:f(a)為降雨量的頻率分布;x為降雨量閾值,mm。使用上述的概率密度計(jì)算式可以確定修正函數(shù),從而有效修正山洪災(zāi)害預(yù)警數(shù)據(jù),提高動(dòng)態(tài)預(yù)警精度。

1.2 確定山洪災(zāi)害動(dòng)態(tài)預(yù)警等級(jí)

待山洪災(zāi)害預(yù)警數(shù)據(jù)修正完畢后,可以確定動(dòng)態(tài)預(yù)警等級(jí),從而提高有效預(yù)警子流域占比。本文設(shè)計(jì)的方法根據(jù)降雨預(yù)報(bào)時(shí)間尺度關(guān)系,構(gòu)建全新的預(yù)警預(yù)報(bào)模式,進(jìn)行合理的產(chǎn)匯流模擬,生成預(yù)報(bào)集合,從而確定山洪災(zāi)害動(dòng)態(tài)預(yù)警等級(jí),進(jìn)行有效的預(yù)警檢驗(yàn)。構(gòu)建的動(dòng)態(tài)預(yù)警模型包含多個(gè)輸入模塊,可以進(jìn)行分區(qū)凈雨計(jì)算、徑流計(jì)算,確定山洪災(zāi)害動(dòng)態(tài)時(shí)間序列關(guān)系,此時(shí)生成的動(dòng)態(tài)預(yù)警模式見(jiàn)圖1。

圖1 研究區(qū)概況

由圖1可知,動(dòng)態(tài)預(yù)警模式可以有效地確定不同時(shí)間段的降雨規(guī)則,使用退水曲線(xiàn)法計(jì)算水位峰值Up,計(jì)算式如下:

(4)

式中:Up為水位峰值,m;C為轉(zhuǎn)換常數(shù);A為洪災(zāi)面積,hm2;Tp為水位處于峰值的時(shí)間,min。

根據(jù)式(4)計(jì)算出的水位峰值,可以對(duì)河道流量進(jìn)行求解轉(zhuǎn)換,確定雨量影響關(guān)系及相對(duì)洪峰偏差,對(duì)現(xiàn)有的預(yù)警結(jié)果進(jìn)行檢驗(yàn)。根據(jù)信息熵原理計(jì)算災(zāi)害動(dòng)態(tài)破壞百分比,轉(zhuǎn)換為山洪災(zāi)害等級(jí),根據(jù)FFAI數(shù)值確定具體的災(zāi)害等級(jí)。為降低信息熵混亂導(dǎo)致的預(yù)警問(wèn)題,本文設(shè)計(jì)的方法使用度量有序原則引入信息論,假設(shè)變量的狀態(tài)完成預(yù)警等級(jí)驗(yàn)證,保證動(dòng)態(tài)預(yù)警的真實(shí)有效性。

1.3 預(yù)警模式實(shí)施建議

為實(shí)現(xiàn)山洪動(dòng)態(tài)預(yù)警,可發(fā)展基于物聯(lián)網(wǎng)(IoT)的氣象觀(guān)測(cè)網(wǎng)絡(luò),建立密集的氣象觀(guān)測(cè)網(wǎng)絡(luò),利用物聯(lián)網(wǎng)技術(shù)將多個(gè)觀(guān)測(cè)站點(diǎn)連接起來(lái),在山區(qū)廣泛布置觀(guān)測(cè)設(shè)備,實(shí)時(shí)收集和傳輸降雨量數(shù)據(jù),提高降雨量的空間分辨率和時(shí)效性,為山洪動(dòng)態(tài)預(yù)警提供更準(zhǔn)確的數(shù)據(jù)支持。還可以強(qiáng)化數(shù)據(jù)集成和分析能力,完善的數(shù)據(jù)集成和分析平臺(tái),能夠?qū)?shí)時(shí)降雨量數(shù)據(jù)與其他相關(guān)數(shù)據(jù)(如地形、土壤水分、河流水位等)進(jìn)行整合和分析,更全面地評(píng)估山洪風(fēng)險(xiǎn),并及時(shí)觸發(fā)預(yù)警機(jī)制。

除了實(shí)時(shí)降雨量預(yù)測(cè),還需要加強(qiáng)信息傳達(dá)和應(yīng)急響應(yīng)能力。建議建立健全的山洪預(yù)警系統(tǒng),包括及時(shí)向當(dāng)?shù)鼐用癜l(fā)布預(yù)警信息、制定有效的應(yīng)急響應(yīng)計(jì)劃等,以最大限度地減少山洪災(zāi)害造成的損失。

2 實(shí)例分析

2.1 概況及準(zhǔn)備

為驗(yàn)證設(shè)計(jì)的考慮降雨量數(shù)值變化的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法的預(yù)警效果,本文選取濟(jì)寧某區(qū)域進(jìn)行了實(shí)例分析。該研究區(qū)域的地形復(fù)雜多樣,北部為沖積平原。山洪災(zāi)害風(fēng)險(xiǎn)區(qū)主要位于濟(jì)寧市東部的曲阜市、泗水縣、鄒城市境內(nèi),由于溝壑縱橫,山體裂隙發(fā)育、斷層較多、地形高差大,巖石裸露多,土層覆蓋薄,植被稀少,遇強(qiáng)降雨極易形成徑流,引發(fā)山洪危及人民群眾的生命財(cái)產(chǎn)安全。據(jù)統(tǒng)計(jì),新中國(guó)成立以來(lái),曲阜市曾在1964年、1965年、1970年、1971年、1973年、2010年發(fā)生山洪災(zāi)害,累計(jì)死亡46人,沖垮、沖毀水庫(kù)、塘壩30座,房屋倒塌241間。泗水縣曾在1957年、1960年、1977年、1991年、2003年等發(fā)生過(guò)較大的山洪災(zāi)害,損失嚴(yán)重。研究區(qū)域的示意見(jiàn)圖2。

圖2 研究區(qū)域

該研究區(qū)域在國(guó)家、山東省大力支持下,開(kāi)始了山洪災(zāi)害防治建設(shè),基本形成了市、縣、鄉(xiāng)一體化的山洪災(zāi)害監(jiān)測(cè)預(yù)警信息管理體系,覆蓋3個(gè)縣(市)34個(gè)鄉(xiāng)鎮(zhèn)713個(gè)村,成為魯水防御濟(jì)寧平臺(tái)內(nèi)容的重要組成部分。全市山洪災(zāi)害風(fēng)險(xiǎn)區(qū)主要涉及曲阜市、泗水縣、鄒城市的38個(gè)鎮(zhèn)725個(gè)村73597人。建有水庫(kù)240座(曲阜市40座,泗水縣87座,鄒城市113座),其中頭頂庫(kù)117座、串聯(lián)庫(kù)71座。1萬(wàn)m3以上的塘壩526座,其中曲阜市108座,泗水縣119座,鄒城市299座。根據(jù)水文氣象站獲取的降雨量數(shù)據(jù),該地區(qū)的氣候特點(diǎn)是四季分明,春冬干旱少雨,年季降雨不均,全市多年平均年降水量700.1mm,極易形成山洪災(zāi)害,特別是在大氣熱力環(huán)流變化、極端天氣頻繁的近些年,防御山洪災(zāi)害的形勢(shì)嚴(yán)峻,任務(wù)更加艱巨。

研究區(qū)域?yàn)?zāi)害防治區(qū)涉及12個(gè)鎮(zhèn)街128個(gè)村,面積約占全市總面積的60%。研究區(qū)域共設(shè)置20個(gè)自動(dòng)雨量站、56個(gè)自動(dòng)水位站、40個(gè)簡(jiǎn)易水位站和50個(gè)簡(jiǎn)易雨量站,以及13個(gè)無(wú)線(xiàn)發(fā)射站和52個(gè)無(wú)線(xiàn)預(yù)警接收站,共同組成縣級(jí)預(yù)警系統(tǒng)平臺(tái)。為防御山洪災(zāi)害,濟(jì)寧市財(cái)政每年列支15萬(wàn)元維修養(yǎng)護(hù)資金,保證系統(tǒng)正常運(yùn)行。濟(jì)寧市水務(wù)局與移動(dòng)公司合作建立山洪災(zāi)害危險(xiǎn)區(qū)短信提醒機(jī)制,在山洪災(zāi)害易發(fā)期,及時(shí)向受山洪災(zāi)害威脅的村莊居民發(fā)送預(yù)警短信;同時(shí),利用大數(shù)據(jù)對(duì)進(jìn)入該區(qū)域的外來(lái)流動(dòng)人群同步發(fā)送短信,提醒做好防御準(zhǔn)備,實(shí)現(xiàn)精準(zhǔn)服務(wù)、全面覆蓋。截至2023年6月,已累計(jì)向群眾發(fā)送短信4萬(wàn)余條,群眾山洪災(zāi)害防范意識(shí)明顯加強(qiáng)。

本文采集的實(shí)例分析數(shù)據(jù)均屬于DEM高程數(shù)據(jù),來(lái)自國(guó)家氣象中心,采用MODIS進(jìn)行歸一化處理,處理后,排除了實(shí)例分析的預(yù)警不確定性,可以構(gòu)建有效的實(shí)例分析模型,計(jì)算式如下:

φ=φkt

(5)

式中:φk為洪水發(fā)生場(chǎng)次;t為降水時(shí)間段,h。

根據(jù)實(shí)例分析模型可以有效模擬研究區(qū)域的山洪過(guò)程。由2022年該區(qū)域降雨量生成模擬預(yù)警徑流量示意圖,見(jiàn)圖3。

圖3 模擬徑流量示意圖

根據(jù)圖3所示的模擬徑流量示意圖可以有效計(jì)算納什效率系數(shù),利用馬爾可夫鏈估計(jì)實(shí)例分析預(yù)警變量。實(shí)例分析參數(shù)確定完畢后,可以劃分山洪動(dòng)態(tài)預(yù)警流域,對(duì)比成功預(yù)警子流域占比,判斷不同動(dòng)態(tài)預(yù)警方法的預(yù)警效果,此時(shí)生成了線(xiàn)性判別函數(shù):

d(x)=w1p+w2s

(6)

式中:w1、w2分別為實(shí)例分析量化參數(shù);p為累計(jì)雨量,mm;s為土壤飽和度,%。

根據(jù)式(6)可以進(jìn)一步生成轉(zhuǎn)移矩陣,確定實(shí)例分析臨界雨量,按照山洪災(zāi)害嚴(yán)重程度將預(yù)警劃分為T(mén)(低等級(jí)藍(lán)色預(yù)警)、2T(中等級(jí)黃色預(yù)警)、3T(高等級(jí)紅色預(yù)警),待全部實(shí)例分析指標(biāo)確定完畢后即可得出最終的分析結(jié)果。

2.2 應(yīng)用結(jié)果與討論

在研究區(qū)域,分別使用本文設(shè)計(jì)的考慮降雨量數(shù)值變化的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法,文獻(xiàn)[6]提出的基于動(dòng)態(tài)臨界雨量的山洪災(zāi)害預(yù)警方法,以及文獻(xiàn)[7]提出的基于Copula函數(shù)風(fēng)險(xiǎn)組合的山洪災(zāi)害預(yù)警方法進(jìn)行預(yù)計(jì),記錄不同預(yù)警方法在不同風(fēng)險(xiǎn)等級(jí)下的有效預(yù)警子流域占比,山洪災(zāi)害動(dòng)態(tài)預(yù)警結(jié)果見(jiàn)表2。

表2 應(yīng)用效果

由表2可知,本文設(shè)計(jì)的考慮降雨量數(shù)值變化的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法在不同預(yù)警等級(jí)下的有效預(yù)警子流域占比較高,文獻(xiàn)[6]提出的基于動(dòng)態(tài)臨界雨量的山洪災(zāi)害預(yù)警方法及文獻(xiàn)[7]提出的基于Copula函數(shù)風(fēng)險(xiǎn)組合的山洪災(zāi)害預(yù)警方法,在不同預(yù)警等級(jí)下的有效預(yù)警子流域占比相對(duì)較低。實(shí) 例分析結(jié)果證明,本文設(shè)計(jì)的考慮降雨量數(shù)值變化的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法的預(yù)警效果較好,具有可靠性,有一定的應(yīng)用價(jià)值。

3 結(jié) 語(yǔ)

綜上所述,我國(guó)山地面積較廣,山洪災(zāi)害牽涉的人群總數(shù)較多,損失較高。常規(guī)的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法難以獲取全部預(yù)警因素,不符合目前的災(zāi)害預(yù)警要求。本文考慮降雨量數(shù)值變化,結(jié)合目前的實(shí)際,設(shè)計(jì)了一種全新的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法。實(shí)例分析結(jié)果表明,設(shè)計(jì)的山洪災(zāi)害動(dòng)態(tài)預(yù)警方法的預(yù)警效果較好,具有可靠性,有一定的應(yīng)用價(jià)值。