遼陽市太子河區(qū)山洪災(zāi)害風(fēng)險評價

李煥玉

(遼寧省水利事務(wù)服務(wù)中心,遼寧 沈陽 110003)

山洪災(zāi)害的復(fù)雜性和不可預(yù)測性使其預(yù)測、預(yù)防成為防汛工作中的難點,由于我國地域廣闊,地理特征差異性明顯,加上山丘區(qū)的經(jīng)濟(jì)社會發(fā)展不均衡、降水分布的時空變異性高和面積較大等因素,使得山洪災(zāi)害頻繁發(fā)生,不僅對基礎(chǔ)設(shè)施造成嚴(yán)重破壞,而且對生命和財產(chǎn)安全構(gòu)成嚴(yán)重威脅[1]。近年來,隨著經(jīng)濟(jì)的發(fā)展和人口的增多,洪澇災(zāi)害對山丘區(qū)的威脅也日益加大,并成為制約可持續(xù)發(fā)展的重要問題之一。評估無資料山區(qū)的山洪災(zāi)害風(fēng)險是災(zāi)害防治工作的延續(xù),也是風(fēng)險普查的重點工作[2]。目前,主要有DEM法、水文模型法、參數(shù)移植法和“從有到無”法等水文預(yù)報方法,其中DEM法主要是利用遙感影像提取流域參數(shù),完成相應(yīng)的水文計算,但該方法對遙感影像的精度要求高,在平原區(qū)的適用性較差;水文模型法的關(guān)鍵是提高模型的獨立性,通過建立水文模型降低對數(shù)據(jù)的依賴度;參數(shù)移植法則基于有資料與無資料地區(qū)在地理位置和流域面積上的相似性,通過參數(shù)移植來進(jìn)行水文分析,但該方法無法全面考慮預(yù)報模型參數(shù)受水文特征或下墊面條件的影響,實施應(yīng)用時具有一定局限性;“從有到無”法是常用的方法,將有資料流域的信息通過插值和外推推算到無資料地區(qū)[3]。對于無資料山區(qū),通常選擇內(nèi)外業(yè)調(diào)查相結(jié)合的方法進(jìn)行山洪災(zāi)害風(fēng)險評價,具體而言,在準(zhǔn)備階段可以使用之前的山洪分析成果,從而確定暴雨洪水參數(shù),然后進(jìn)行預(yù)警指標(biāo)分析、防洪現(xiàn)狀評價及設(shè)計暴雨洪水計算等,對以上各環(huán)節(jié)選擇合適的方法是實現(xiàn)山洪風(fēng)險評價的關(guān)鍵。

鑒于此,本文結(jié)合沿河宅基地高程、河道斷面、下墊面條件及成災(zāi)水位等數(shù)據(jù),采用地形圖和MIKE11軟件計算不同頻率洪水淹沒范圍及水位,并進(jìn)一步明確預(yù)警指標(biāo)以及各等級危險區(qū),旨在為當(dāng)?shù)胤篮橹卫砗蜕胶轭A(yù)警提供一定支持。

1 區(qū)域概況

太子河區(qū)位于太子河沖積平原上,總土地面積273.37 km2,地勢東南高,西北低,自東向西傾斜。該區(qū)域?qū)俅箨懶约撅L(fēng)氣候區(qū),年均降水量739.7 mm,年際間變化較大,多集中于6-9月,占全年的73%。境內(nèi)有大小河流8條,其中鐵西街道有太子河段7.4 km、城市總干渠、城市北排水,祁家鎮(zhèn)有南地河、機(jī)場排水、城市南北排水,東寧衛(wèi)鄉(xiāng)有南地河,沙嶺鎮(zhèn)有太子河段14 km、柳壕河5.7 km、兵馬河6 km,王家鎮(zhèn)有太子河段13 km、北沙河7.8 km、渾沙河1 km。根據(jù)《遼陽市城市防洪預(yù)案》,位于太子河右岸的王家鎮(zhèn)輝寧堡村為太子河第二人工分洪點,堤防高約6m,人工分洪長度100 m。滯洪區(qū)為北沙河以西、太子河以北、小小公路以南,小渾河以東地區(qū),集水面積113.0 km2,可蓄水量2.30×108m3,淹沒時間28 h。

太子河上游地區(qū)植樹造林,攔截降水,減少地面徑流,延緩洪峰發(fā)生時間和減小峰量;太子河干流和支流上修建有葠窩、觀音閣、湯河三座大型水庫,蓄水總?cè)萘?6.1×108m3,防洪庫容25.2×108m3,可有效地調(diào)蓄洪水,三庫聯(lián)合調(diào)度削減洪峰,使下游地區(qū)發(fā)生超標(biāo)準(zhǔn)洪水的頻次大大減少。由于特殊的地形和氣候條件,太子河區(qū)容易受到臺風(fēng)和突發(fā)強(qiáng)降雨影響,引發(fā)山洪災(zāi)害,對山區(qū)居民安全構(gòu)成潛在威脅。研究區(qū)主要河流有太子河和北沙河,對防洪能力有重要影響,因此選擇太子河、北沙河進(jìn)行水位流量關(guān)系的評估分析,以更好地掌握太子河區(qū)的山洪災(zāi)害風(fēng)險,并制定相應(yīng)的預(yù)警和防護(hù)措施。

2 研究方法

2.1 洪水計算

(1) 設(shè)計暴雨。根據(jù)當(dāng)?shù)刈儾钕禂?shù)、降雨量均值等值線圖以及定點定面關(guān)系表,計算不同頻率下各流域的設(shè)計面雨量[4]。其中,降雨均值等值線圖可提供不同地點的平均降雨量分布情況,變差系數(shù)等值線圖揭示了降雨的變異性和不確定性程度。

(2) 設(shè)計洪水。針對無資料地區(qū),應(yīng)充分考慮流域面積大小選擇相應(yīng)的計算方法[5]:流域面積<50 km2,采用推理公式推導(dǎo)的關(guān)系式計算設(shè)計洪水;面積>50 km2,使用瞬時單位線法通過考慮不同時段內(nèi)的單位面積洪峰流量,更準(zhǔn)確地估計設(shè)計洪水;面積>300 km2時,為保證計算精度,將流域劃分為多個單元進(jìn)行計算,這還可以降低計算復(fù)雜性。

2.2 防洪現(xiàn)狀分析

通過計算設(shè)計洪水評價區(qū)域防洪能力,結(jié)合評價結(jié)果統(tǒng)計不同危險區(qū)的房屋和人口數(shù),為制定洪災(zāi)防御、人員轉(zhuǎn)移及臨時安置預(yù)案提供數(shù)據(jù)參考。通過分析不同水位下的洪峰流量及其與沿河橋涵、房屋、道路之間的關(guān)系,統(tǒng)計確定設(shè)計洪水下的承載水位、人口和房屋數(shù)全面分析防洪現(xiàn)狀。

洪水位計算。通過率定和驗證水動力學(xué)模型參數(shù)計算不同重現(xiàn)期洪水,結(jié)合分析結(jié)果明確各斷面最高水位,并對比各斷面成災(zāi)水位與洪峰水位,選取防洪最薄弱斷面為控制斷面,按照成災(zāi)水位判定區(qū)內(nèi)防洪能力,參照流量-水位關(guān)系曲線計算成災(zāi)流量。對不同重現(xiàn)期設(shè)計暴雨和洪水進(jìn)行試算,從而確定設(shè)計洪峰量,若設(shè)計洪峰量與某重現(xiàn)期成災(zāi)流量之差處于誤差允許范圍,則該斷面能夠達(dá)到要求的防洪能力。

劃定淹沒危險區(qū)。結(jié)合各頻率洪水位及地形地貌條件,在地圖上劃定極高、高和一般危險區(qū)。

淹沒分析。根據(jù)不同頻率最高洪水位建立流量-水位線,綜合分析DEM數(shù)據(jù)、歷史水災(zāi)、局部地形、堤防高程和沿河居民調(diào)查表等資料,建立淹沒范圍內(nèi)的水位-人口關(guān)系線,將兩者相結(jié)合繪制防洪現(xiàn)狀圖,可以為決策者制定更有效的防洪措施和應(yīng)對策略提供參考。

2.3 預(yù)警分析方法

(1) 預(yù)警指標(biāo)。通常包括臨界水位和雨量,可以分為準(zhǔn)備和立即轉(zhuǎn)移兩個級別。對于無資料地區(qū),臨界水位主要來源于水位監(jiān)測站,但是太子河區(qū)還未設(shè)立監(jiān)測站。所以,本文選擇臨界雨量作為預(yù)警指標(biāo)進(jìn)行研究。

(2) 土壤含水量。流域土壤的濕度一般和較濕兩種情況分別選用Pa=0.75 Wm、Pa=0.90 Wm進(jìn)行表征,其中Wm、Pa為流域最大蓄水量和和流域蓄水量。

(3) 預(yù)警時段。這取決于當(dāng)?shù)氐南聣|面、形狀系數(shù)、平均比降、流域面積及暴雨特性等因素,一般選擇0.5 h、1 h、3 h和6 h作為預(yù)警時段。

(4) 臨界雨量。常用的分析方法有等雨量法、統(tǒng)計歸納法和經(jīng)驗分析法等,由于水文水力學(xué)法充分分析率物理機(jī)制,其預(yù)警精度以及對數(shù)據(jù)資料的要求較高?？紤]到河道特性、下墊面條件和降雨數(shù)據(jù)的充分性,本次選用水文水力學(xué)法來確定臨界雨量,即利用曼寧公式計算預(yù)警水位對應(yīng)的流量,并進(jìn)一步反推出臨界雨量作為預(yù)警指標(biāo)。

結(jié)合以往的轉(zhuǎn)移和預(yù)警情況,短時間轉(zhuǎn)移危險區(qū)全部人員難以實現(xiàn)。所以,將立即轉(zhuǎn)移因子設(shè)定為臨界雨量的90%,并將其向前推0.5 h作為準(zhǔn)備轉(zhuǎn)移的閾值,從而能夠提前預(yù)警,使人們做好轉(zhuǎn)移準(zhǔn)備,降低潛在的災(zāi)害風(fēng)險。

2.4 數(shù)據(jù)來源

3 結(jié)果與分析

3.1 太子河區(qū)洪水分析計算

典型歷時選取10 min、1 h、6 h和24 h,并推算5 a、10 a、20 a、50 a和100 a一遇的暴雨頻率,然后利用遼陽市短歷時暴雨中的數(shù)據(jù)確定太子河區(qū)不同歷時的點暴雨量和變異系數(shù)Cv均值,結(jié)合雨強(qiáng)系數(shù)Kp及點面轉(zhuǎn)換系數(shù)計算設(shè)計面暴雨值如表1。

表1 設(shè)計暴雨量計算值

根據(jù)太子河區(qū)具體情況和小流域山洪特點,在產(chǎn)流部分凈雨量計算過程中選擇蓄滿產(chǎn)流模型,且不考慮蒸發(fā)量損失。研究選用的北沙河、太子河面積均不超過50 km2,故對首個斷面設(shè)計洪水使用推理公式法進(jìn)行計算,結(jié)果如表2。

表2 不同頻率下的設(shè)計洪峰流量

3.2 太子河區(qū)現(xiàn)狀防洪能力

(1) 建立水動力模型。根據(jù)太子河區(qū)所處的河勢、地勢、地形和流域位置等情況,為了達(dá)到洪水計算要求共設(shè)置30個斷面,其中太子河22個和北沙河8個。為了更加客觀地反映太子河區(qū)的河道情況,充分考慮其他支流匯入情況以及堰壩、橋梁等因素影響,采用MIKE11軟件建立一維水動力學(xué)模型并生成該畫圖。概化范圍為下游水位邊界至太子河區(qū)起點斷面,概化斷面為實測斷面。另外,設(shè)定上邊界條件為各頻率的設(shè)計流量,表2為設(shè)計流量的數(shù)據(jù)來源,固定的水位選用洪水的下邊界條件。

(2) 劃定淹沒危險區(qū)。通過洪水計算確定最高水位,通過比較洪峰水位與各斷面的成災(zāi)水位,并考慮太子河區(qū)地形圖可明確各等級下的淹沒范圍。

(3) 淹沒分析。在最薄弱防洪斷面選擇時,將控制斷面選為兩條河流中最先淹沒的斷面,在此基礎(chǔ)上從MIKE系統(tǒng)中提取各斷面最高水位,采用前文所述的計算方法進(jìn)行洪水位的相關(guān)計算,結(jié)果如表3。

表3 各流域控制斷面設(shè)計洪峰水位

根據(jù)洪水計算成果和淹沒范圍內(nèi)的水位-人口關(guān)系,可以利用控制斷面的水位-流量關(guān)系繪制太子河區(qū)的防洪現(xiàn)狀評價圖。

3.3 太子河區(qū)預(yù)警雨量確定

考慮到太子河區(qū)內(nèi)的河流洪水流速較大的特點,使得上、下游斷面之間的洪水傳播時間較短。此外,由于缺乏自動水位站,對上游無法獲得準(zhǔn)確的水位數(shù)據(jù),但該區(qū)域建設(shè)的雨量站較為完善,可以獲得可靠的雨量數(shù)據(jù)。因此,本文選擇預(yù)警指標(biāo)為臨界雨量值,結(jié)合相關(guān)技術(shù)要求確定匯流時間以及0.5 h、1.0 h、3.0 h、6.0 h作為預(yù)警時段。采用前文所述方法反推臨界雨量、立即和準(zhǔn)備轉(zhuǎn)移指標(biāo)(見表4)。

鑒于北沙河預(yù)警指標(biāo)普遍較高的情況,為了保證太子河區(qū)的安全,將北沙河的預(yù)警指標(biāo)作為太子河區(qū)的預(yù)警標(biāo)準(zhǔn),洪水來襲時,太子河區(qū)居民將根據(jù)北沙河的預(yù)警指標(biāo)來采取相應(yīng)的預(yù)防和保護(hù)措施,使其在面對洪水威脅時能夠更加迅速、有效地做出反應(yīng),確保太子河區(qū)整體安全。

3.4 結(jié)果分析

采用上述方法評價太子河區(qū)的防洪現(xiàn)狀,根據(jù)地形圖、沿河宅基地高程以及Arc GIS軟件劃定危險區(qū),并統(tǒng)計各區(qū)域內(nèi)的人口和房屋。通過推算各鄉(xiāng)鎮(zhèn)臨界雨量評定其防洪能力,如表5。

表5 太子河區(qū)防洪現(xiàn)狀

通過危險區(qū)的劃分、人口和房屋分布情況以及臨界雨量的推算,可以有針對性地提出防洪措施和加強(qiáng)應(yīng)急準(zhǔn)備,分析結(jié)果有助于決策者制定更加有效的山洪防災(zāi)方案。

4 結(jié)語

本文結(jié)合歷史山洪災(zāi)害調(diào)查、現(xiàn)場實測以及調(diào)查報告等資料數(shù)據(jù),采用MIKE11軟件和水動力學(xué)模型評價了太子河區(qū)的防洪預(yù)警指標(biāo)及現(xiàn)狀防洪能力,主要結(jié)論如下:

(1) 針對各河流建立水動力模型,經(jīng)參數(shù)率定和驗證分析確定洪水演變過程以及該模型的適用性,結(jié)果顯示該模型適用于無資料山區(qū),其可操作性和適用性良好。

(2) 太子河區(qū)的總體防洪能力為20～50 a一遇,而王家鎮(zhèn)防洪能力較低為8 a一遇,該評價結(jié)果與區(qū)域?qū)嶋H情況基本一致。因此,太子河區(qū)山洪防御建設(shè)的重點為王家鎮(zhèn),可以為區(qū)域山洪防御和預(yù)警提供技術(shù)支持。

(3) 由于涉及的河流缺乏長期水文監(jiān)測數(shù)據(jù),仍需要對本次評價結(jié)果進(jìn)行驗證和率定,使評價成果更加優(yōu)化和完善。此外,由于實際預(yù)警工作受多個因素影響,如防御措施、暴雨時空分布差異以及土壤前期含水量等,應(yīng)考慮實際情況動態(tài)優(yōu)化雨量預(yù)警指標(biāo)。