鐘 華，嵇海祥，王旭瀅，商華嶺，4

(1.南京水利科學(xué)研究院，江蘇 南京 210029；2.水利部南京水利水文自動(dòng)化研究所，江蘇 南京 210029；3.上?？睖y(cè)設(shè)計(jì)研究院有限公司，上海 200335；4.河海大學(xué)水文水資源學(xué)院，江蘇 南京 210024)

1 研究背景

山區(qū)流域性洪水每年在我國(guó)發(fā)生的各類災(zāi)害中都占有很大比例，對(duì)地區(qū)的社會(huì)經(jīng)濟(jì)健康發(fā)展和群眾生命安全構(gòu)成了極大的威脅[1]。與平原地區(qū)相比，山區(qū)地形環(huán)境更為復(fù)雜，導(dǎo)致洪水流速變化快，水位暴漲暴落，其帶來(lái)的洪水災(zāi)害因此具有突發(fā)性強(qiáng)、破壞力大等特點(diǎn)[2]。借助洪水?dāng)?shù)值模擬模型對(duì)流域不同量級(jí)洪水進(jìn)行模擬分析和災(zāi)害風(fēng)險(xiǎn)評(píng)估，能夠?yàn)樯絽^(qū)洪水防御及避險(xiǎn)轉(zhuǎn)移提供重要的技術(shù)支撐，提高山區(qū)防災(zāi)減災(zāi)的管理能力，從而保障流域安全[3-4]。

近些年來(lái)，我國(guó)對(duì)山區(qū)流域洪水災(zāi)害已開展了大量的研究工作，對(duì)洪水的成因、損失評(píng)估、流域分析評(píng)價(jià)、風(fēng)險(xiǎn)區(qū)劃管理等都取得了較大的進(jìn)展[5-8]。一般認(rèn)為，高強(qiáng)度、歷時(shí)長(zhǎng)的降雨一般是引起流域洪水災(zāi)害的主要原因；因此，我國(guó)東部季風(fēng)區(qū)更易受到洪水災(zāi)害的侵?jǐn)_[9]。而東部沿海地區(qū)人口密集，社會(huì)經(jīng)濟(jì)活動(dòng)廣泛，洪水災(zāi)害發(fā)生后其潛在的影響范圍也會(huì)更廣。對(duì)東部山區(qū)流域進(jìn)行有效的洪水分析和風(fēng)險(xiǎn)評(píng)估，有助于進(jìn)一步加快補(bǔ)齊水旱災(zāi)害防御體系建設(shè)的短板，可為防洪規(guī)劃、防洪工程措施、非工程措施等提供科學(xué)依據(jù)。

本文以浙江省甌江流域青田段作為研究對(duì)象，通過(guò)構(gòu)建一、二維水文水動(dòng)力模型，對(duì)設(shè)計(jì)洪水情景進(jìn)行分析模擬，并依據(jù)社會(huì)經(jīng)濟(jì)和人口等資料，對(duì)其風(fēng)險(xiǎn)影響進(jìn)行定量評(píng)估。

2 研究區(qū)概況

青田縣位于浙江省東南部，地處甌江中下游，全縣總面積2 493.8 km2。其中，丘陵低山占89.7%，河溪、塘、庫(kù)占5%，其余為平地，有“九山半水半分田”之稱。

甌江流域?qū)賮啛釒Ъ撅L(fēng)氣候區(qū)，流域多年平均降水量一般在1 500～2 100 mm之間。梅雨和臺(tái)風(fēng)暴雨是甌江流域洪水的主要成因，青田段主要以8月～9月份的臺(tái)風(fēng)暴雨造成的洪澇為主，尤其是正面襲擊的臺(tái)風(fēng)，伴隨而來(lái)的風(fēng)暴潮，使甌江河道水位暴漲，極易造成洪澇災(zāi)害。

甌江穿過(guò)青田縣境內(nèi)，主要有大溪(甌江干流)、甌江及支流小溪。本文研究范圍涉及大溪、甌江及支流小溪等，涉及河長(zhǎng)77.6 km，覆蓋河道兩岸面積65.2 km2(見(jiàn)圖1)。沿江青田縣城現(xiàn)狀防洪能力約為20年一遇；沿江鄉(xiāng)鎮(zhèn)現(xiàn)狀防洪能力約為5～20年一遇。

圖1 研究范圍示意

3 模型構(gòu)建

洪水來(lái)源為上游洪水(大溪、小溪)和區(qū)間暴雨洪水，為了能夠模擬河道洪水，建立一維二維水動(dòng)力學(xué)模型。河道洪水采用一維水動(dòng)力學(xué)模型、河道外洪水淹沒(méi)采用二維水動(dòng)力學(xué)模型，通過(guò)一二維耦合計(jì)算，模擬分析河道洪水演進(jìn)及淹沒(méi)。

3.1 一維河道建模

河道一維模型搭建范圍上游為大溪開潭水電站、下游為甌江溫州潮位站，河道建模范圍包括大溪、小溪、四都港、甌江干流等。模型構(gòu)建中河道斷面間距不超過(guò)500 m。

模型的水文邊界的準(zhǔn)確控制是提高洪水風(fēng)險(xiǎn)分析精度的必要條件。為確定合理的邊界條件，必須充分考慮河道的連續(xù)性與整體性(見(jiàn)圖2)。一維模型主要水文邊界條件如下：大溪上邊界采用開潭水電站壩址斷面，設(shè)計(jì)洪水移用五里亭水文站設(shè)計(jì)洪水流量過(guò)程；小溪上邊界采用灘坑水庫(kù)壩址斷面，設(shè)計(jì)洪水采用灘坑水庫(kù)設(shè)計(jì)下泄流量過(guò)程疊加灘坑壩址至三溪口斷面的區(qū)間設(shè)計(jì)洪水流量過(guò)程；支流四都港上邊界為秋蘆斷面，設(shè)計(jì)洪水采用秋蘆站設(shè)計(jì)洪水流量過(guò)程；支流船寮港上邊界為紅光斷面，設(shè)計(jì)洪水采用紅光斷面以上集雨面積設(shè)計(jì)暴雨推求的設(shè)計(jì)洪水流量過(guò)程；官莊源、臘溪源、禎埠港、官坑源、?？谠?、雄溪源、高市源、芝溪源、石蓋源、大路源、石溪源等11個(gè)匯入大溪的支流，湖邊源、水碓坑、石郭源、港頭坑、貴岙源、石洞源等6個(gè)匯入甌江的支流均作為支流邊界，采用洪水流量過(guò)程線；甌江下邊界采用甌江溫州站，采用溫州站設(shè)計(jì)潮位過(guò)程(見(jiàn)圖3)。按照同頻原則，即干流發(fā)生某一設(shè)計(jì)頻率洪水、支流發(fā)生相應(yīng)頻率設(shè)計(jì)洪水進(jìn)行組合。其余支流均為流量邊界，采用洪水流量過(guò)程線，由區(qū)間設(shè)計(jì)暴雨推求設(shè)計(jì)洪水得到，共設(shè)置19個(gè)。

圖2 洪水計(jì)算邊界條件設(shè)置示意

圖3 一維河道建模范圍

3.2 二維模型概化

二維模型以DEM數(shù)字高程模型為依據(jù)，充分考慮線性工程的阻水(公路、鐵路、堤防等線狀物)的阻水及導(dǎo)水影響?？紤]到甌江洪水、小溪洪水及流域內(nèi)暴雨、灘坑水庫(kù)泄洪等，二維模型范圍以大溪、小溪沿岸按照75 m等高線、甌江沿岸按照45 m等高線為依據(jù)確定(見(jiàn)圖4)。

圖4 二維模型網(wǎng)格劃分

在對(duì)河道兩側(cè)堤防外部的保護(hù)區(qū)進(jìn)行網(wǎng)格劃分時(shí)，根據(jù)GIS地圖使用2D模擬多邊形對(duì)象繪出兩岸洪水可能淹沒(méi)范圍的外形輪廓，自動(dòng)計(jì)算出面積大小，并根據(jù)實(shí)際情況修訂面積值及形狀。根據(jù)面積大小，確定最大三角形面積、最小網(wǎng)格面積和最小角度；在經(jīng)濟(jì)與人口分布密集區(qū)段，進(jìn)行適當(dāng)網(wǎng)格加密(見(jiàn)圖5)。

圖5 二維模型網(wǎng)格高程插值

3.3 一二維耦合設(shè)置

在創(chuàng)建了斷面、2D多邊形、溢流單元這些單一對(duì)象之后，使用特定對(duì)象將這些單一對(duì)象連接起來(lái)：離散的斷面使用連接相連；斷面和溢流單元之間使用堰流公式連接相連；干支流交匯處使用交叉點(diǎn)進(jìn)行連接。同時(shí)，利用計(jì)算模型自帶的模型工具對(duì)斷面的方向、角度以及左右岸標(biāo)記進(jìn)行修正，完成洪水模擬計(jì)算網(wǎng)絡(luò)的構(gòu)建。

4 模型率定與驗(yàn)證

考慮地形及水文資料的實(shí)際條件，選取“20140820”洪水，根據(jù)模型范圍內(nèi)的鶴城水文站、禎埠水文站的實(shí)測(cè)水位、流量過(guò)程，對(duì)模型進(jìn)行率定；同時(shí)，沿河選取多個(gè)具有實(shí)測(cè)最高水位的斷面，進(jìn)行水位最高值的比對(duì)分析。

“20140820”洪水是指從2014年8月18日開始，連日暴雨引發(fā)的洪水。其對(duì)麗水市的影響至20日8時(shí)，麗水主要干支流重要水位站大部分已超警戒水位；甌江干流持續(xù)高水位并不斷增高，麗水市主城區(qū)大水門站水位51.91 m，流量已達(dá)50年一遇洪水標(biāo)準(zhǔn)；20日21時(shí)青田縣禎埠站水位34.88 m，鶴城站水位12.06 m，鶴城站流量10 140 m3/s。

選擇沿江水文站斷面，取實(shí)測(cè)的水位、流量過(guò)程與模型計(jì)算的水位、流量過(guò)程進(jìn)行比較分析。本次我們選取模型范圍內(nèi)的鶴城水文站、禎埠水文站的實(shí)測(cè)水位、流量過(guò)程進(jìn)行率定。

模型模擬成果見(jiàn)圖6。各站水位及流量計(jì)算值與實(shí)測(cè)值吻合良好，基本滿足模型驗(yàn)證要求。同時(shí)，計(jì)算的長(zhǎng)歷時(shí)潮位及流量相關(guān)關(guān)系與實(shí)測(cè)同步性也很好，再現(xiàn)了各站水位流量時(shí)間上的變化特征。

圖6 鶴城站實(shí)測(cè)結(jié)果與計(jì)算結(jié)果過(guò)程對(duì)比

從水文站實(shí)測(cè)水位、流量過(guò)程看，模型模擬結(jié)果基本與實(shí)測(cè)過(guò)程吻合；從沿江12個(gè)斷面實(shí)測(cè)最高水位值與計(jì)算值比較(見(jiàn)表1)，其中有10個(gè)斷面的計(jì)算最高水位值與實(shí)際洪水的最高水位誤差≤30 cm。

表1 模型計(jì)算結(jié)果比較

選用2019年7月9日至10日的洪水對(duì)模型進(jìn)行驗(yàn)證。2019年7月9日至10日，甌江流域發(fā)生暴雨洪水，青田縣境內(nèi)大溪沿岸部分鄉(xiāng)鎮(zhèn)受淹嚴(yán)重。上邊界開潭水庫(kù)由于暫無(wú)流量監(jiān)測(cè)，其入庫(kù)洪水流量的估算根據(jù)上游小白巖和秋潭水文站的流量進(jìn)行線性疊加，并考慮洪水傳播時(shí)間。根據(jù)水文站至開潭水庫(kù)的距離，洪水傳播時(shí)間約為1 h。下邊界溫州潮位站選用實(shí)測(cè)潮位數(shù)據(jù)作為邊界控制條件。此外，考慮支流小溪上游灘坑水庫(kù)的滯蓄作用，水庫(kù)出流按照200 m3/s設(shè)置。

模型計(jì)算結(jié)果和實(shí)測(cè)數(shù)據(jù)比較結(jié)果顯示，模型能夠較好地模擬甌江流域在洪水及下游潮位頂托共同作用下的水動(dòng)力分布特性(見(jiàn)圖7)。綜上分析，所建立的一二維水動(dòng)力模型基本符合甌江青田段的洪水分析要求。

圖7 控制斷面實(shí)測(cè)水位與模擬過(guò)程對(duì)比

5 洪水情景模擬和風(fēng)險(xiǎn)評(píng)估

5.1 洪水情景模擬

5.1.1 設(shè)計(jì)洪水

設(shè)計(jì)洪水以流量資料推求(流量法)為主，區(qū)間支流設(shè)計(jì)洪水用暴雨資料推求(雨量法)。本文選取甌江青田段100年一遇設(shè)計(jì)洪水進(jìn)行情景模擬和洪水風(fēng)險(xiǎn)評(píng)估(見(jiàn)圖8)。

圖8 設(shè)計(jì)洪水流量過(guò)程線

河道一維模型上邊界為大溪開潭水庫(kù)下泄流量，由于開潭水庫(kù)集雨面積為8 544 km2，而五里亭水文站的集雨面積為8 870 km2，面積比為0.96。因此，移用五里亭水電站相關(guān)頻率下的設(shè)計(jì)流量作為大溪開潭水庫(kù)的設(shè)計(jì)洪水，采用面積比進(jìn)行修正。典型洪水選擇1992年8月31日歷史洪水過(guò)程。計(jì)算支流小溪的設(shè)計(jì)洪水時(shí)，由于小溪上建有灘坑水庫(kù)，因此分為兩部分計(jì)算：灘坑水庫(kù)設(shè)計(jì)下泄流量以及灘坑至三溪口區(qū)間設(shè)計(jì)洪水。按照最不利原則，將二者進(jìn)行線性疊加至灘坑壩址作為小溪的計(jì)算邊界。

灘坑水庫(kù)的下泄流量

Q泄=Q安全-DQ1

(1)

式中，Q泄為灘坑水庫(kù)下泄流量，包含機(jī)組發(fā)電流量，調(diào)洪計(jì)算時(shí)按400 m3/s計(jì)；Q安全為下游圩仁站斷面處的安全泄量；DQ1為考慮洪水傳播時(shí)間后的灘坑—圩仁站區(qū)間流量(包括大溪干流洪水流量)。此外，計(jì)算考慮灘坑—圩仁站區(qū)間洪水流量的預(yù)報(bào)誤差約為6%。

灘坑至三溪口區(qū)間及其他19個(gè)支流設(shè)計(jì)洪水用暴雨資料推求。暴雨取樣為年最大值法，取最大24 h雨量，面雨量計(jì)算采用面積加權(quán)平均法。暴雨系列選用1953年～2011年共59年。暴雨頻率分析采用P-Ⅲ型理論曲線擬合適線，求得流域設(shè)計(jì)暴雨成果，計(jì)算使用24 h概化雨型。流域產(chǎn)流計(jì)算采用蓄滿產(chǎn)流原理的簡(jiǎn)易扣損法，匯流計(jì)算主要采用浙江省瞬時(shí)單位線法。

5.1.2 設(shè)計(jì)潮位

甌江河口屬于山溪性強(qiáng)潮河口，潮汐為非正規(guī)淺海半日潮。干流上中游水位主要受徑流控制，下游河口潮流段，主要受潮汐控制。資料統(tǒng)計(jì)表明，溫州站具有較好的代表性，選用測(cè)站年最高水(潮)位進(jìn)行頻率分析(見(jiàn)圖9)，經(jīng)P-Ⅲ型理論曲線適線，求得年最高設(shè)計(jì)潮位成果。設(shè)計(jì)潮型采用典型年法；選取偏不利(1992年8月30日～9月2日)實(shí)測(cè)潮型進(jìn)行同倍比放大。

圖9 溫州站設(shè)計(jì)潮位過(guò)程線

5.1.3 洪潮組合

甌江流域下游洪水宣泄速度與洪潮遭遇情況相關(guān)密切。若遇上大潮，洪水受大潮頂托，江道水位會(huì)明顯抬高。研究區(qū)域內(nèi)年最大洪水與年最高潮位相遇的年份有：1971年、1987年、1990年、1992年、1996年等。經(jīng)水文(潮位)站資料分析，甌江青田段洪水有可能受大潮的頂托。因此，在選取洪潮組合時(shí)，選取5年一遇高潮位進(jìn)行模擬計(jì)算。

5.2 洪水模擬分析

模擬甌江100年一遇設(shè)計(jì)洪水遭遇溫洲站5年一遇高潮位的情景，甌江青田段防洪體系中可能出現(xiàn)的洪水風(fēng)險(xiǎn)，為評(píng)估青田現(xiàn)狀防洪能力提供參考。下圖為根據(jù)模擬成果繪制的洪水風(fēng)險(xiǎn)圖，包括淹沒(méi)水深和淹沒(méi)歷時(shí)(見(jiàn)圖10、11)。

圖10 甌江流域100年一遇設(shè)計(jì)洪水淹沒(méi)水深

模擬結(jié)果分析顯示，洪水淹沒(méi)主要發(fā)生在大溪、甌江兩岸的鄉(xiāng)鎮(zhèn)、村莊；這是因?yàn)樾∠嫌斡袨┛铀畮?kù)的防洪調(diào)蓄，因而兩岸洪水淹沒(méi)較小。大溪、甌江兩岸的低洼區(qū)域，如禎埠、?？凇⒋嫉妊亟l(xiāng)鎮(zhèn)，在100年一遇設(shè)計(jì)洪水情景下發(fā)生了不同程度的洪水淹沒(méi)。

5.3 損失評(píng)估

洪水風(fēng)險(xiǎn)評(píng)估的基礎(chǔ)資料主要包括基礎(chǔ)地理信息(行政區(qū)劃、居民點(diǎn)、高程、道路交通、流域水系等)、水文資料(降雨、水位流量等)、構(gòu)筑物及工程調(diào)度規(guī)則、洪澇災(zāi)害資料(歷史洪水、淹沒(méi)范圍、水深)和社會(huì)經(jīng)濟(jì)數(shù)據(jù)等。

結(jié)合洪水分析的結(jié)果，對(duì)受洪水淹沒(méi)影響的人口、資產(chǎn)、設(shè)施及洪水災(zāi)害造成的經(jīng)濟(jì)損失進(jìn)行評(píng)估。甌江青田段江沿岸情況較特殊，主要沿江行政村聚集區(qū)作為防洪關(guān)注重點(diǎn)，分析時(shí)重點(diǎn)考慮甌江干流段(青田縣城區(qū)域)的沿江行政村，進(jìn)行經(jīng)濟(jì)、人口等的損失評(píng)估。

圖11 甌江流域100年一遇設(shè)計(jì)洪水淹沒(méi)歷時(shí)

在100年一遇設(shè)計(jì)洪水情景下，青田縣城兩岸研究范圍內(nèi)如圖12所示。受淹村鎮(zhèn)集中在支流小溪和四都港匯入干流的區(qū)域，洪水位已高過(guò)絕大多數(shù)堤防標(biāo)準(zhǔn)，發(fā)生淹沒(méi)的區(qū)域涉及4個(gè)街道(三溪口、甌南、鶴城、油竹)共14個(gè)社區(qū)(村)，共計(jì)影響1.35萬(wàn)人，淹沒(méi)區(qū)GDP約328 720.76萬(wàn)元，由于洪水造成的損失約131 488.3萬(wàn)元。

圖12 青田縣100年一遇設(shè)計(jì)洪水時(shí)受洪水影響區(qū)域示意

在此情景下，青田縣各級(jí)防控指揮單位應(yīng)結(jié)合實(shí)際組織人員巡堤查驗(yàn)，做好應(yīng)急物資儲(chǔ)備、搶險(xiǎn)隊(duì)伍準(zhǔn)備，保障青田縣城區(qū)等最重要保護(hù)對(duì)象的防洪安全，提前組織淹沒(méi)區(qū)域的人員及時(shí)轉(zhuǎn)移避險(xiǎn)。

6 結(jié) 論

本文以浙江省甌江青田段為研究區(qū)，通過(guò)構(gòu)建一、二維水文水動(dòng)力模型，對(duì)研究范圍內(nèi)的洪水風(fēng)險(xiǎn)進(jìn)行模擬分析，評(píng)估山區(qū)流域性洪水對(duì)地區(qū)的社會(huì)經(jīng)濟(jì)影響。情景模擬分析充分考慮了臺(tái)風(fēng)型洪水影響的實(shí)際情況，通過(guò)洪潮組合對(duì)水文邊界條件進(jìn)行設(shè)置，并借助高精度DEM數(shù)據(jù)構(gòu)建二維水動(dòng)力學(xué)模型。經(jīng)歷史洪水的率定驗(yàn)證，本研究所建模型能夠較好地模擬研究流域內(nèi)的洪水淹沒(méi)特征，這為未來(lái)青田縣的流域洪水防御工作提供了技術(shù)支撐。

根據(jù)模型模擬的100年一遇設(shè)計(jì)洪水淹沒(méi)成果，大溪、甌江等兩岸沿江低洼區(qū)域仍然是流域防洪重點(diǎn)關(guān)注區(qū)域，是加強(qiáng)防洪短板和應(yīng)急搶險(xiǎn)工作的優(yōu)先對(duì)象。