廖一丁 輝

(1.南京信息工程大學大氣科學學院，江蘇 南京 210044；2.南京信息工程大學應用水文氣象研究院，江蘇 南京 210044)

近幾十年來，我國由暴雨造成的洪澇災害事件頻繁發(fā)生，且呈顯著加劇的趨勢，對人民的生命財產(chǎn)安全造成了嚴重的威脅[1]。廣東省是中國大陸最南部的沿海省份，東西北三面連接內(nèi)陸，南面瀕臨南海，鄰近西太平洋，受西風帶天氣系統(tǒng)和熱帶天氣系統(tǒng)影響頻繁，加上南嶺及其他山脈的地形作用，降水量十分充沛，暴雨洪澇成為廣東省最主要的氣象災害之一。尤其是短歷時暴雨所引發(fā)的洪澇災害突發(fā)性強，危害性更高。由于廣東省地形復雜，不同區(qū)域受天氣系統(tǒng)的影響程度各有差異，使得降水的時空分布不均勻[2]，給暴雨洪澇災害的預警和防范造成一定難度。因此，對廣東省的暴雨進行頻率分析并研究其空間分布特征，有助于防災減災工作的開展，減少暴雨洪澇災害帶來的損失。

過去我國常用的水文頻率分析方法是“單站、單時段、單一線型、基于常規(guī)矩的適線法”，線型通常選擇皮爾遜Ⅲ型曲線[3]。然而單站分析需要較長的序列資料，序列長度的不足和觀測站點的缺乏很大程度上影響了降雨頻率估計值的可靠性[4]，而且不能反映空間分布規(guī)律；常規(guī)矩估計的統(tǒng)計參數(shù)具有很大的偏態(tài)性，因而得到的降雨頻率估計值極不穩(wěn)定[3]。地區(qū)線性矩法是Hosking等[5]在線性矩的基礎上提出的一種地區(qū)頻率分析方法。美國國家海洋大氣管理總署(NOAA)自1991年起在美國開展地區(qū)線性矩法在防洪設計標準的應用研究，從1997年起分區(qū)分批對全國的暴雨頻率圖集進行更新，于2006年提出了一套基于次序統(tǒng)計量的線性矩法結(jié)合基于水文氣象一致區(qū)的地區(qū)分析法進行暴雨頻率分析的完整系統(tǒng)[3,6-7]。近年來，國內(nèi)有不少學者利用降雨資料進行了地區(qū)線性矩法的應用研究[8-16]，證明了該方法的優(yōu)越性：線性矩相比常規(guī)矩具有良好的估計參數(shù)的不偏性和對特大值的穩(wěn)健性；地區(qū)分析法能夠充分利用鄰近站點的信息，提高降雨頻率估計值的準確性，并且能更準確地獲得其空間分布。

暴雨高風險區(qū)劃是基于水文氣象地區(qū)線性矩頻率分析法得到的一個地區(qū)內(nèi)某一定歷時、一定頻率(重現(xiàn)期)降雨事件中最大雨強的空間分布[12,15]。實際上，暴雨高風險區(qū)劃的概念為定量估算“強降雨雨強、降雨落區(qū)、概率”三者的關(guān)系奠定了理論基礎，填補了設計暴雨理論中強降雨雨強空間分布的空白；而現(xiàn)有的暴雨空間分布僅探討暴雨的點—面關(guān)系。暴雨高風險區(qū)劃為中小流域洪水尤其是山洪防治、預警提供了理論基礎。

關(guān)于廣東省降雨頻率分析的研究，Yang等[17]以珠江流域為研究區(qū)，采用地區(qū)線性矩法對1 d、3 d、5 d和7 d年最大降水量進行頻率分析，并聯(lián)系氣候背景及地形特征，探析流域極端降雨的時空分布規(guī)律；黃強等[18]利用廣東省年最大日降水量資料，應用地區(qū)線性矩法，得到不同流域分區(qū)100 a或500 a以下重現(xiàn)期可靠性較高的降雨頻率估計值。這些都是對1 d及以上長歷時極值降雨的頻率分析，而對小時歷時降雨的頻率分析比較少。

本文應用水文氣象地區(qū)線性矩頻率分析法以及在此基礎上發(fā)展的暴雨高風險區(qū)劃概念和技術(shù)，根據(jù)1 h、6 h、12 h和24 h年最大降水量資料對廣東省進行暴雨高風險區(qū)劃研究，繪制暴雨高風險區(qū)劃圖并分析其空間分布特征，以期為廣東省的防洪規(guī)劃設計和洪澇災害早期預警提供參考。

1 資料與方法

1.1 研究資料

廣東省陸地面積約17.97萬km2，為提高研究區(qū)邊界處降雨頻率估計值的可靠性，本文將研究區(qū)廣東省向內(nèi)陸擴大約30～50 km的范圍作為緩沖區(qū)。在研究區(qū)和緩沖區(qū)收集水文雨量站1 h、6 h、12 h 和24 h的歷史年最大降水量資料。經(jīng)質(zhì)量控制，去掉明顯錯誤的數(shù)據(jù)值，舍去不一致的資料序列，最終篩選出資料序列長度在20 a及以上的站點用于分析，其中具備1 h年最大降水量資料的有202個站點(其中廣東省內(nèi)有92個站點)，具備6 h、12 h和24 h年最大降水量資料的有291個站點(其中廣東省內(nèi)有157個站點)。研究區(qū)及周圍的地形和所有站點的分布如圖1，可見廣東省內(nèi)站點的空間分布比較均勻，緩沖區(qū)站點密度略大于廣東省內(nèi)站點密度。站點降水序列范圍從1939—2015年，長度不等，平均長度約34 a。

圖1 廣東省及周圍緩沖區(qū)地形和站點

1.2 研究方法

Hosking[19]將線性矩定義為次序統(tǒng)計量線性組合的期望值，并定義了類似于常規(guī)矩的統(tǒng)計特征參數(shù)：線性矩離差系數(shù)CvL、線性矩偏態(tài)系數(shù)CsL和線性矩峰度系數(shù)CkL。基于水文氣象地區(qū)線性矩法的暴雨高風險區(qū)劃研究，以線性矩作為參數(shù)估計方法，在水文氣象一致區(qū)劃分的基礎上，選擇各一致區(qū)的最優(yōu)分布線型，再根據(jù)地區(qū)分析法推求各站點的降雨頻率估計值，并進行時空一致性調(diào)整，最終獲得不同重現(xiàn)期下降雨的空間分布。

1.2.1水文氣象一致區(qū)劃分

地區(qū)分析法的前提是劃分水文氣象一致區(qū)。由于不同時段分析的站點不一定完全相同，降雨的分布特征也不相同，一致區(qū)的劃分應分不同時段進行。水文氣象一致區(qū)的劃分和判定主要從以下3個方面進行：

a.氣象相似性。一致區(qū)內(nèi)滿足降雨的水汽入流和氣象成因一致。

b.水文相似性。一致區(qū)內(nèi)所有站點的統(tǒng)計參數(shù)CvL和CsL在一定的容忍度內(nèi)一致。采用基于CvL計算值的異質(zhì)性檢驗指標H1來判斷所劃分的子區(qū)是否為一致性區(qū)域[5]。當H1<1時，表示該子區(qū)為可以接受的一致區(qū)。同時，利用各站點的CsL對子區(qū)進行判別和調(diào)整，選出CsL特大值和特小值對應的站點，分別考察刪除該站點前后，整個子區(qū)站點是否存在100 a重現(xiàn)期降雨頻率估計值小于實測資料序列中的最大降水量值的不合理情況。若刪除該站點后，這種不合理情況有較大改善，則需將該雨量站點移至相鄰子區(qū)分析，否則可以保留。

c.不和諧性檢驗。檢驗所劃分的子區(qū)內(nèi)是否存在不和諧的站點，若站點的不和諧性指標Di超過一定的臨界值，則可認為該站點是不和諧站點。Di具體計算方法和子區(qū)內(nèi)不同站點數(shù)對應的臨界值可參考Hosking等[5]的研究。對不和諧站點進一步檢查原數(shù)據(jù)的可靠性，并考慮將其調(diào)整至相鄰子區(qū)或單獨分區(qū)。如果該站點較大的Di值是由局部極端氣象事件所引起的，也可以保留在當前區(qū)域。

1.2.2一致區(qū)線型選擇

選擇5種三參數(shù)分布線型作為一致區(qū)的候選分布線型，根據(jù)其尾端形態(tài)由厚至薄依次為：廣義邏輯分布(GLO)、廣義極值分布(GEV)、廣義正態(tài)分布(GNO)、廣義帕累托分布(GPA)和皮爾遜Ⅲ型分布(P-Ⅲ)。采用以下3種擬合優(yōu)度檢驗方法來選擇各一致區(qū)的最優(yōu)分布線型：

a.蒙特卡洛模擬檢驗。該檢驗通過比較一致區(qū)內(nèi)區(qū)域平均的線性矩峰度系數(shù)與分布函數(shù)的線性矩峰度系數(shù)之間的差異來考察分布函數(shù)擬合的質(zhì)量。檢驗的統(tǒng)計量ZDIST的具體計算方法可參考Hosking和Wallis的研究[5]。若統(tǒng)計量滿足|ZDIST|≤1.64，認為該分布函數(shù)擬合結(jié)果是合理可接受的，并且|ZDIST|越接近于0，認為擬合效果越好。

b.樣本線性矩的均方根誤差檢驗。該檢驗利用一致區(qū)內(nèi)各站點的線性矩峰度系數(shù)與分布函數(shù)的線性矩峰度系數(shù)之間的差值，再根據(jù)站點資料序列長度進行加權(quán)平均得到的均方根誤差(RMSE)來比較分布函數(shù)的擬合效果[7]，具有最小RMSE的分布函數(shù)擬合效果最好。

c.實測數(shù)據(jù)檢驗。該檢驗首先分別計算一致區(qū)內(nèi)各站點不同重現(xiàn)期下的經(jīng)驗頻率與候選分布函數(shù)理論頻率之間的相對誤差(RE)，取所有站點的RE平均值[7]；再將5種候選分布函數(shù)在2 a、5 a、10 a、25 a和50 a重現(xiàn)期的RE平均值從大到小排列，各重現(xiàn)期下的排列序號求和作為相應分布函數(shù)的RE分數(shù)SRE，SRE越高、RE越小，表明該分布函數(shù)擬合效果越好。

最后，綜合3種檢驗方法的結(jié)果，確定各一致區(qū)擬合效果最佳的分布函數(shù)。

1.2.3降雨頻率估計值計算及時空一致性調(diào)整

b.時段間一致性調(diào)整。由于不同時段分布函數(shù)的適線是獨立進行的，獲得的不同時段的降雨頻率估計值曲線可能會出現(xiàn)交叉的情況，即從交叉點之后，時段較短的降雨頻率估計值比時段較長的降雨頻率估計值大，這與實際情況不符。本文采用“誤差分攤”的方法[7]調(diào)整時段間不一致的降雨頻率估計值，其思路為：計算不一致起點前一頻率下相鄰的較長時段與較短時段降雨頻率估計值的比值，將比值大于1的誤差部分，按頻率步長權(quán)重分配到不一致起點之后的各不一致頻率點上，加上1作為相應頻率下新的比值，乘以原較短時段降雨頻率估計值，即得到較長時段調(diào)整后的降雨頻率估計值，而其他正常的降雨頻率估計值部分不改變。

c.空間一致性調(diào)整。同理，由于雨量站點資料有限、站點分布不均勻，以及各一致區(qū)的數(shù)據(jù)是獨立進行分布函數(shù)選擇和參數(shù)化的，可能會造成計算得到的降雨頻率估計值在相鄰一致區(qū)邊界處出現(xiàn)不連續(xù)、梯度較大的現(xiàn)象。本文采用“往返兩次”空間平差法[20]來調(diào)整降雨頻率估計值的空間不一致性，其步驟如下：①構(gòu)造一個與站點分辨率大致相同的空間網(wǎng)格，采用克里金插值方法，將不規(guī)則的站點上的降雨頻率估計值插值到規(guī)則網(wǎng)格點上；②采用反距離加權(quán)插值法，利用規(guī)則網(wǎng)格點上的降雨頻率估計值反向插值回各站點，即得到站點空間平差校正后的降雨頻率估計值。

2 結(jié)果與分析

2.1 水文氣象一致區(qū)劃分

首先，根據(jù)廣東省的地形和氣候初步劃分水文氣象一致區(qū)。從廣東省地形(圖1)來看，整體地勢從北部山地向南部沿海呈逐步降低趨勢。北部群山是南嶺的組成部分，東部山地由三列東北—西南走向的山脈構(gòu)成，分別為九連山、羅浮山和蓮花山[1]。廣東省暴雨的主要成因是鋒面類暴雨和臺風類暴雨，水汽主要來源于南面的南海[21]。通常南嶺南側(cè)、蓮花山東南坡等迎風面降水量較大，而背風面的谷底和內(nèi)陸盆地降水量較少。另外珠江三角洲平原是一個尺度很大的南開喇叭口地形，對氣流有輻合抬升作用，使降水量和降雨強度加大[1]。綜合以上分析，考慮將北部、東部和西南部劃分不同的一致區(qū)，并且喇叭口地形、山脈迎風面和背風面之間要進行區(qū)分。

在初步分區(qū)的基礎上，分不同時段、利用各站點的CvL和CsL進一步細分和調(diào)整。圖2和圖3給出了站點CvL和CsL的空間分布，盡可能將CvL和CsL接近的站點劃分為一個子區(qū)，對子區(qū)內(nèi)CsL特大值和特小值的站點進行考察，并用異質(zhì)性檢驗指標H1檢驗子區(qū)是否滿足一致區(qū)標準。經(jīng)過多次反復調(diào)整，最終確定了各時段水文氣象一致區(qū)的劃分方案，針對 1 h、6 h、12 h和24 h年最大降水量序列分別劃分了13、15、16和16個一致區(qū)(圖2中以藍色實線劃分的區(qū)域表示)。不同時段的分區(qū)總體形態(tài)上相似，其中存在的差異可能是由于長歷時站點數(shù)比短歷時站點數(shù)多，且CvL和CsL的梯度較大，因此在一些區(qū)域需要更細的劃分。表1列出了各時段各子區(qū)的站點數(shù)、異質(zhì)性檢驗指標H1值和不和諧指標Di超過臨界值的站點數(shù)。所有子區(qū)H1<1，表明都可以認為是一致區(qū)。1 h、6 h和12 h年最大降水量序列的部分子區(qū)存在少數(shù)不和諧站點，其中1 h時段的第三個一致區(qū)(1 h-3區(qū))的大仚站Di值最高，達4.41(相應的臨界值為3)。通過對原數(shù)據(jù)序列分析發(fā)現(xiàn)，該站2003年的最大1 h降水量為94.5 mm，同時也是該子區(qū)所有站點數(shù)據(jù)序列中的最大值，造成了該站的CvL和CsL值較周圍站點偏大，所以Di值偏高。但該站序列長度僅為21年，不考慮單獨分區(qū)，因此還是將其保留在1 h-3區(qū)。其他子區(qū)的不和諧站點的Di值只是略高于臨界值，且分區(qū)過程中已經(jīng)將這些站點劃入不同的子區(qū)進行比較，從中選擇了站點不和諧度最小的方案。

(a) 1 h

表1 廣東省各水文氣象一致區(qū)的站點數(shù)及異質(zhì)性檢驗和不和諧性檢驗結(jié)果

(a) 1 h

2.2 一致區(qū)內(nèi)最優(yōu)頻率分布線型的選擇

對所劃分的一致區(qū)，采用蒙特卡洛模擬檢驗、均方根誤差檢驗和實測數(shù)據(jù)檢驗，考察5種候選分布線型GLO、GEV、GNO、GPA和P-Ⅲ的擬合效果。以24 h水文氣象一致區(qū)為例，各一致區(qū)3種檢驗的結(jié)果見表2。如24 h時段的第一個一致區(qū)(24 h-1區(qū))，3種檢驗中都是GEV分布表現(xiàn)最好，因此該一致區(qū)的最優(yōu)分布為GEV。又如24 h-5區(qū)，GEV分布在蒙特卡洛模擬檢驗和實測數(shù)據(jù)檢驗中表現(xiàn)最好，而RMSE最小和次小的分布分別為GNO和GEV，但兩者RMSE相差不大，綜合3種檢驗結(jié)果，選取GEV為該一致區(qū)的最優(yōu)分布。對各時段、各一致區(qū)進行類似的分析，最終可以確定所有一致區(qū)的最優(yōu)分布(表3)。由表3可以看出，選擇作為最優(yōu)分布最多的是GEV(占60%)，其次是GNO(占38%)，只有24 h-12區(qū)選擇了GPA。

表2 廣東省24 h水文氣象一致區(qū)擬合優(yōu)度檢驗結(jié)果

表3 廣東省水文氣象一致區(qū)的最優(yōu)分布線型

2.3 降雨頻率估計值計算及時空不一致性調(diào)整

由上一步選擇的各時段各一致區(qū)的最優(yōu)分布函數(shù)推求的地區(qū)增長因子，與相應時段各站點的年最大降水量平均值相乘，即得到各站點在各重現(xiàn)期下的降雨頻率估計值(1 h、6 h、12 h和24 h的降雨頻率估計值分別表示為Q1 h、Q6 h、Q12 h和Q24 h)。圖4給出了紫洞站和犁市站各時段降雨頻率估計值隨重現(xiàn)期變化的曲線。可以看出，同一時段，降雨頻率估計值隨重現(xiàn)期的增加而增加，相同重現(xiàn)期下，降雨頻率估計值也隨時段的增加而增加，這是符合統(tǒng)計規(guī)律的。

(a) 紫洞站

根據(jù)前文的分析，此時計算得到的降雨頻率估計值會存在時間或空間不一致的問題，因此需要進行時空一致性檢驗及調(diào)整。從較短時段到較長時段，依次檢查各站點兩相鄰時段降雨頻率估計值的一致性，對站點存在的較長時段比較短時段降雨頻率估計值小的異常部分，采用“誤差分攤”的方法進行調(diào)整。最終對5站的Q12 h、34站的Q24 h進行了調(diào)整。這里以百候站為例分析時段間一致性調(diào)整的過程。圖5給出了該站調(diào)整前后降雨頻率估計值曲線的對比，該站缺少1 h年最大降水量數(shù)據(jù)，因此只分析Q6 h、Q12 h和Q24 h，曲線數(shù)據(jù)點上的數(shù)值為對應重現(xiàn)期下Q24 h和Q12 h的比值。調(diào)整前，Q24 h和Q12 h曲線在重現(xiàn)期100 a到200 a之間出現(xiàn)了交叉，從200 a到10 000 a，Q24 h小于Q12 h，即兩者的比值小于1。于是將100 a重現(xiàn)期處Q24 h和Q12 h的比值誤差分配到200 a到10 000 a之間各重現(xiàn)期上，得到新的比值及Q24 h。調(diào)整后，Q24 h和Q12 h曲線不再交叉，Q24 h和Q12 h的比值都大于1。經(jīng)過時段間一致性調(diào)整后的降雨頻率估計值更符合降雨的統(tǒng)計特性，增強了不同時段間降雨頻率估計值的可比性。

(a) 調(diào)整前

在時段間調(diào)整的基礎上，采用“往返兩次”空間平差法調(diào)整空間的不一致性。根據(jù)各時段的站點密度，第一步插值的網(wǎng)格分辨率1 h數(shù)據(jù)取0.33°×0.33°(約33 km×33 km)，6 h、12 h和24 h數(shù)據(jù)取0.25°×0.25°(約25 km×25 km)。這里以24 h數(shù)據(jù)的50 a重現(xiàn)期為例考察降雨頻率估計值空間調(diào)整的效果(圖6)，調(diào)整前降雨頻率估計值的等值線在部分一致區(qū)邊緣處會存在梯度突然增大的現(xiàn)象，而調(diào)整后，一致區(qū)邊緣處的梯度有所緩和，但整體上主要的分布形勢和大值中心不變。

(a) 調(diào)整前

2.4 暴雨高風險區(qū)劃圖分析

將時空一致性調(diào)整后各站點在不同時段、不同重現(xiàn)期下的降雨頻率估計值，通過克里金法空間插值，獲得相應的降雨頻率估計值的空間分布，即暴雨高風險區(qū)劃圖，其中最大雨強區(qū)域為暴雨高風險區(qū)。圖7～10給出了在1 h、6 h、12 h和24 h時段25 a、50 a 和100 a重現(xiàn)期下的暴雨高風險區(qū)劃圖。總體上看，同一時段不同重現(xiàn)期下，降雨頻率估計值的空間分布態(tài)勢基本一致；不同時段相同重現(xiàn)期進行比較，6 h、12 h和24 h時段下的降雨頻率估計值的空間分布態(tài)勢較相似，而與短歷時1 h降雨頻率估計值的空間分布存在一定差異。所有時段最明顯的暴雨高風險區(qū)位于廣東西南沿海陽江附近(圖7～10中的C1)，隨著重現(xiàn)期增大，其范圍逐漸延伸至西南沿海陽江—江門一帶。6 h、12 h和24 h降雨頻率估計值的空間分布還呈現(xiàn)出兩個較明顯的暴雨高風險區(qū)，分別位于珠江三角洲廣州到北部山區(qū)南側(cè)清遠附近區(qū)域(圖8～10中的C2)和東部沿海陸豐附近(圖8～10中的C3)，且重現(xiàn)期越大，兩個區(qū)域的范圍越大、中心值越大。不同時段不同重現(xiàn)期下暴雨高風險區(qū)中心的降雨頻率估計值見表4，可見4個時段C1處中心值都最大，為最主要的暴雨高風險區(qū)，6 h、12 h和24 h時段下位于C2和C3處的中心值相當，C2處中心值在重現(xiàn)期較小時小于C3處，而隨著重現(xiàn)期增大變?yōu)榇笥贑3處中心值。6 h、12 h和24 h時段下的降雨頻率估計值在3個暴雨高風險區(qū)的南側(cè)都出現(xiàn)了低值區(qū)，分別位于西部的云浮—肇慶附近、北部的樂昌附近和東部的五華附近。這3個暴雨高風險區(qū)和3個低值區(qū)的位置與廣東省的3個多雨中心和3個少雨中心位置相一致，形成的原因與不同地區(qū)主要受影響的天氣系統(tǒng)不同，以及地形的阻擋、輻合抬升作用有關(guān)[22-23]，表明此暴雨高風險區(qū)劃圖符合廣東省降雨的空間分布特征。6 h、12 h和24 h時段暴雨高風險區(qū)C2所在的珠江三角洲區(qū)域，在1 h時段沒有呈現(xiàn)明顯的高值中心，說明該區(qū)域主要受較長歷時暴雨的影響。相反，1 h時段雷州半島在不同重現(xiàn)期的降雨頻率估計值都比較高，說明雷州半島主要受短歷時暴雨影響。另外，廣東東部由于1 h序列所用資料的站點比較稀疏，一些局地分布的變化可能顯示不出，無法確定是否存在高風險區(qū)。由此可以說明，不同區(qū)域主要的成災暴雨歷時是不同的[24]，實際應用時，應充分調(diào)查不同區(qū)域歷史主要暴雨、洪澇災害事件，得到成災暴雨的平均歷時和平均降水量，在此基礎上選擇相應設計時段和重現(xiàn)期下的暴雨高風險區(qū)劃圖進行重點分析。

(a) 25 a

(a) 25 a

(a) 25 a

(a) 25 a

表4 廣東省暴雨高風險區(qū)中心降雨頻率估計值

3 結(jié) 論

a.依據(jù)水文氣象一致區(qū)劃分準則，基于1 h、6 h、12 h和24 h年最大降水量序列，分別將廣東省及其周圍緩沖區(qū)劃分為13、15、16和16個水文氣象一致區(qū)；采用3種擬合優(yōu)度檢驗方法，確定了各水文氣象一致區(qū)的最優(yōu)分布線型，其中GEV和GNO是選擇最多的兩種分布。

b.根據(jù)地區(qū)分析法，計算各站點不同時段、不同重現(xiàn)期下的降雨頻率估計值，并采用“誤差分攤”法和“往返兩次”空間平差法進行時段間和空間的一致性調(diào)整，最終獲得時間上更符合統(tǒng)計特性、空間分布更合理的降雨頻率估計值成果。

c.相同時段降雨頻率估計值的空間分布態(tài)勢在不同重現(xiàn)期下基本一致，而1 h時段與6 h、12 h和24 h時段的降雨頻率估計值空間分布存在一定差異。廣東省多小時時段最主要的暴雨高風險區(qū)位于西南沿海陽江—江門區(qū)域；6 h、12 h和24 h時段的第二、第三個暴雨高風險區(qū)分別位于珠江三角洲廣州到北部山區(qū)南側(cè)清遠附近和東部沿海陸豐附近區(qū)域。該暴雨高風險區(qū)劃圖符合廣東省降雨的空間分布特征，可為防洪規(guī)劃設計以及洪澇災害風險早期預警提供參考。