無(wú)實(shí)測(cè)資料斷面設(shè)計(jì)洪水推求方法探討

許繼良 徐 慧

（1. 浙江省水利水電勘測(cè)設(shè)計(jì)院 杭州 310002；2. 常山縣水利局 浙江衢州 324200）

1 問(wèn)題的提出

對(duì)于無(wú)實(shí)測(cè)資料斷面設(shè)計(jì)洪水，小流域洪水一般采用地區(qū)暴雨洪水綜合法，較大流域洪水通常采用面積指數(shù)法（流量法）推求設(shè)計(jì)洪水。實(shí)際運(yùn)用中，由于受水文測(cè)站分布等因素影響，僅采用傳統(tǒng)水文學(xué)方法推求無(wú)實(shí)測(cè)資料斷面設(shè)計(jì)洪水，其參數(shù)的率定、驗(yàn)證通常較為困難，對(duì)設(shè)計(jì)成果精度影響較大。

某流域河長(zhǎng)141km，流域面積2436km2。分大、小兩源，主源集水面積1376km2，支源集水面積 644km2。大、小兩源匯合口下游的沙頭鎮(zhèn)集水面積約2140km2，其斷面設(shè)計(jì)洪水的準(zhǔn)確性關(guān)乎到下游縣城核心區(qū)域防洪安全。

沙頭斷面無(wú)水文站實(shí)測(cè)流量資料，主源設(shè)有石柱水文站，控制流域面積 1273km2，自 1956年有實(shí)測(cè)流量資料，同時(shí)有1853、1905、1920、1922年等歷史洪水資料。支源設(shè)有碧蓮站，控制流域面積 433km2，實(shí)測(cè)流量資料系列為1957～1993年，同時(shí)有1912、1920年等歷史洪水資料。以往規(guī)劃設(shè)計(jì)中均采用面積指數(shù)法推求沙頭斷面設(shè)計(jì)洪水，以石柱水文站作為參證站，面積指數(shù)采用經(jīng)驗(yàn)指數(shù)0.8，20年一遇沙頭斷面設(shè)計(jì)洪峰流量約為 12000m3/s。該設(shè)計(jì)成果被普遍運(yùn)用于下游防洪工程的規(guī)劃設(shè)計(jì)中。從水文設(shè)計(jì)角度而言，該方法符合我省以往水文設(shè)計(jì)經(jīng)驗(yàn)，具有一定的科學(xué)性。但由于斷面無(wú)實(shí)測(cè)資料，且設(shè)計(jì)斷面控制流域面積與參證站控制流域面積相差較大，面積指數(shù)按省經(jīng)驗(yàn)指數(shù)選取，對(duì)設(shè)計(jì)成果準(zhǔn)確性可能帶來(lái)一定的影響。

本文試通過(guò)水文、水力學(xué)結(jié)合分析方法，充分利用石柱水文站、碧蓮水文站實(shí)測(cè)流量資料，以及區(qū)間各降雨站降雨資料，通過(guò)構(gòu)建水文、水動(dòng)力耦合模型，將無(wú)實(shí)測(cè)資料的沙頭斷面水文設(shè)計(jì)模型拓展至有實(shí)測(cè)資料的上塘潮位站水文、水動(dòng)力耦合設(shè)計(jì)模型，以上塘潮位站實(shí)測(cè)水位對(duì)設(shè)計(jì)模型進(jìn)行率定、驗(yàn)證，從而求得沙頭設(shè)計(jì)斷面歷史洪峰過(guò)程，再以傳統(tǒng)的水文統(tǒng)計(jì)方法求得沙頭斷面各設(shè)計(jì)頻率洪峰過(guò)程。

2 流域基本資料

2.1 水文站

流域主源上有石柱水文站，支源上有碧蓮水文站。兩站于 1956年設(shè)立，進(jìn)行水位、流量、降水、泥沙等項(xiàng)目觀(guān)測(cè)，流域面積分別為1273km2和 433km2。

2.2 降水量站

設(shè)計(jì)流域內(nèi)降水量站有上塘、潭頭、里樟等站，各降水量站統(tǒng)轄于浙江省水文勘測(cè)局，觀(guān)測(cè)資料經(jīng)統(tǒng)一整編、校對(duì)，精度可靠，可作為工程設(shè)計(jì)的基本資料。同時(shí)，各降水量站域內(nèi)分布基本合理，基本反映設(shè)計(jì)流域的降水特性，為工程設(shè)計(jì)提供了可靠的依據(jù)。

2.3 水位站

域內(nèi)設(shè)計(jì)斷面沙頭鎮(zhèn)下游設(shè)有上塘潮位站，擁有近 20年水位觀(guān)測(cè)資料，可作為計(jì)算分析工作的驗(yàn)證站。各測(cè)站分布情況見(jiàn)圖1。

圖1 流域水系測(cè)站分布圖

3 分析方法及成果

本文通過(guò)構(gòu)建水文水動(dòng)力耦合計(jì)算模型，以有實(shí)測(cè)資料的碧蓮站及石柱站作為水動(dòng)力模型上邊界，河口潮位站作為計(jì)算下邊界，區(qū)間小流域設(shè)計(jì)洪水采用傳統(tǒng)的水文暴雨推求洪水模型與水動(dòng)力模型進(jìn)行耦合，利用設(shè)計(jì)斷面下游上塘站實(shí)測(cè)水位作為耦合模型率定、驗(yàn)證數(shù)據(jù)，計(jì)算歷史沙頭洪水，再輔以傳統(tǒng)水文統(tǒng)計(jì)方法推求沙頭斷面各頻率設(shè)計(jì)洪水。該方法一定程度上解決了無(wú)實(shí)測(cè)資料斷面推求設(shè)計(jì)洪水參數(shù)選用及模型率定、驗(yàn)證問(wèn)題。

3.1 水文模型參數(shù)率定及驗(yàn)證

本文選用碧蓮有實(shí)測(cè)資料的，距今時(shí)間較近，洪峰流量較大的1987年及1992年兩場(chǎng)較大洪水對(duì)水文模型進(jìn)行率定、驗(yàn)證。驗(yàn)證成果見(jiàn)圖2，圖 3。

3.2 水動(dòng)力計(jì)算模型率定

對(duì)河道水動(dòng)力計(jì)算模型采用 2004年云娜臺(tái)風(fēng)洪水為驗(yàn)證洪水，模擬成果見(jiàn)表1.

由表1結(jié)果可得，各特征點(diǎn)實(shí)測(cè)水位與模型計(jì)算水位差均在 0.1m以?xún)?nèi)，模擬精度較好，計(jì)算模型所采用的特征參數(shù)可用于成果計(jì)算。

3.3 耦合模型計(jì)算驗(yàn)證

采用上述經(jīng)參數(shù)率定后的計(jì)算模型，對(duì)1961～2009年大洪水進(jìn)行模擬，為去除初值計(jì)算帶來(lái)的影響，使得計(jì)算盡量符合實(shí)際發(fā)生情況，選用較大洪水，采用永嘉縣上塘水位站有實(shí)際觀(guān)測(cè)的水位對(duì)計(jì)算成果再次進(jìn)行驗(yàn)證、比較，成果見(jiàn)表2。由表2成果可得，采用水文水動(dòng)力耦合計(jì)算模型模擬歷史各場(chǎng)較大洪水，用上塘水位站實(shí)測(cè)水位進(jìn)行驗(yàn)證對(duì)比，各場(chǎng)洪水計(jì)算值與實(shí)測(cè)值差距均在0.3m以?xún)?nèi)，說(shuō)明模型對(duì)流域產(chǎn)、匯流模擬較好。

圖2 1987年模型率定成果對(duì)比圖

圖3 1992年模型驗(yàn)證成果對(duì)比圖

表1 云娜臺(tái)風(fēng)洪水實(shí)測(cè)、計(jì)算水位對(duì)比表 單位：m

3.4 沙頭設(shè)計(jì)洪水推求

計(jì)算得到的沙頭斷面設(shè)計(jì)歷史年最大洪峰流量成果見(jiàn)表3。對(duì)1960～2009年沙頭段面流量進(jìn)行適線(xiàn)，成果見(jiàn)圖4。各頻率洪峰流量見(jiàn)表4。

表2 歷年上塘水位站水位與計(jì)算水位對(duì)比成果表 單位：m

表3 1960～2009年沙頭斷面計(jì)算最大洪峰流量成果表 單位：m3/s

圖4 沙頭斷面流量適線(xiàn)圖

表4 沙頭斷面各頻率洪峰流量成果表

由于水文水動(dòng)力計(jì)算耦合模型考慮了石柱、碧蓮水文站實(shí)際發(fā)生洪峰的錯(cuò)時(shí)，水文站以下河道對(duì)洪水的坦化作用以及河口潮位的頂托作用，又以能夠掌握的各種實(shí)測(cè)數(shù)據(jù)對(duì)模型進(jìn)行了率定、驗(yàn)證，其所計(jì)算出來(lái)的沙頭流量基本能反映實(shí)際發(fā)生的流量過(guò)程。

經(jīng)分析，沙頭斷面20年一遇流量9230m3/s，50年一遇洪峰流量11200m3/s，100年一遇洪峰流量12700m3/s，300年一遇洪峰流量14900m3/s。相比傳統(tǒng)面積指數(shù)法求得的沙頭斷面洪峰流量，以20年一遇洪峰流量為例，小約20%。

4 結(jié)語(yǔ)

在推求無(wú)實(shí)測(cè)資料設(shè)計(jì)洪水時(shí)，單純以水文學(xué)方法，受流域測(cè)站分布等因素制約，其參數(shù)率定、驗(yàn)證較為困難，給設(shè)計(jì)成果帶來(lái)了更多的不確定性。本文提出構(gòu)建流域水文水動(dòng)力耦合模型的設(shè)計(jì)思路，將流域水文站、河道、參證站等要素通過(guò)模型結(jié)合起來(lái)，推求設(shè)計(jì)斷面歷史洪水，再輔以水文統(tǒng)計(jì)方法推求設(shè)計(jì)洪水的思路，經(jīng)某流域斷面設(shè)計(jì)洪水項(xiàng)目驗(yàn)證，被證明具有相對(duì)較高的精度，可為今后類(lèi)似情況洪水設(shè)計(jì)提供一種設(shè)計(jì)思路。

1 孫映紅. 基于MIKE11HD和NAM耦合模型在河流施工圍堰對(duì)防洪安全影響分析中的應(yīng)用與研究. 浙江水利科技，（162）.

2 劉曾美. 感潮河段洪潮遭遇組合風(fēng)險(xiǎn)研究. 中山大學(xué)學(xué)報(bào)，49（2）.

3 劉連梅. 面積比法計(jì)算設(shè)計(jì)斷面洪水中面積指數(shù)的確定.水利水電技術(shù)，2005，36（4）.

4 陶新紅. 無(wú)實(shí)測(cè)資料小流域推求設(shè)計(jì)洪水探討. 河南水利與南水北調(diào)，2010（8）.

5 李燕紅. 雨洪法計(jì)算設(shè)計(jì)洪水適用面積的探討.人民黃河，26（11）.