最新某流域水文預(yù)報(bào)方案編制

1、水文預(yù)報(bào)大作業(yè)某流域水文預(yù)報(bào)方案編制本文在9場(chǎng)洪水資料的基礎(chǔ)上，完成對(duì)某流域的水文預(yù)報(bào)方案的編制。運(yùn)用由一層土壤蒸發(fā)模式、蓄滿產(chǎn)流模型、時(shí)段單位線模型以及地下線性水庫模型組成的二水源新安江模型進(jìn)行參數(shù)率定及檢驗(yàn)工作。通過雨洪關(guān)系的合理性分析，取得其中8場(chǎng)合格的洪水，運(yùn)用遺傳算法對(duì)其中6場(chǎng)洪水進(jìn)行參數(shù)率定，并預(yù)留兩場(chǎng)洪水進(jìn)行參數(shù)檢驗(yàn)。經(jīng)檢驗(yàn)，產(chǎn)流參數(shù)，匯流參數(shù)對(duì)兩場(chǎng)預(yù)留洪水的產(chǎn)流模擬、 匯流模擬均取得了 100% 合格的效果。但由于率定及檢驗(yàn)的洪水資料較少，率定的參數(shù)在預(yù)報(bào)適用的時(shí)間上并不 具有代表性，且預(yù)報(bào)模型具有一定的缺陷。因此，本預(yù)報(bào)方案僅可用來參考。關(guān)鍵詞：水文預(yù)報(bào)方案；遺傳算法；二水源

2、新安江模型摘 要弓I 言01流域及水文概況11.1 流域水系概況11.2 下墊面概況11.3 洪水特性情況12水文資料處理22.1 降雨資料處理22.2 流量資料處理32.3 雨洪關(guān)系的合理性分析63預(yù)報(bào)模型參數(shù)率定83.1 模型的選擇83.2 參數(shù)率定及步驟103.2.1 率定方法選擇 103.2.2 參數(shù)及范圍 103.2.3 參數(shù)率定步驟 114預(yù)報(bào)方案134.1 產(chǎn)流預(yù)報(bào)方案134.2 匯流預(yù)報(bào)方案144.3 預(yù)報(bào)方案分析及評(píng)定16水文預(yù)報(bào)是根據(jù)已知的信息對(duì)未來一定時(shí)期內(nèi)的水文狀態(tài)作出定性或定量的預(yù)測(cè)。 水文預(yù)報(bào)方法以水文基本規(guī)律、水文模型研究為基礎(chǔ)，結(jié)合生產(chǎn)實(shí)際問題的需要，構(gòu)成 具體

3、的預(yù)報(bào)方法或預(yù)報(bào)方案，服務(wù)于生產(chǎn)實(shí)際。一般水文預(yù)報(bào)的研究的重點(diǎn)和關(guān)鍵有兩 部分：(1)共性規(guī)律研究，即具有一定普遍性的水文基本規(guī)律模擬方法和流域水文模型 研究；(2)個(gè)性問題研究，對(duì)反映問題的特征、方法進(jìn)行了解，構(gòu)成具有解決各種具體 實(shí)際問題的、具有較高預(yù)報(bào)精度的預(yù)報(bào)方案。其在防汛?？购怠⑺Y源開發(fā)利用、國(guó)民經(jīng)濟(jì)建設(shè)和國(guó)防等領(lǐng)域都有廣泛的應(yīng)用，經(jīng)濟(jì)效益巨大，應(yīng)用單位眾多。水文預(yù)報(bào)在生產(chǎn)上的應(yīng)用領(lǐng)域十分廣泛，主要有：流域或區(qū)域性洪水與旱情預(yù)測(cè)， 河道、水庫、湖泊等水體的封凍。開動(dòng)狀況及冰凌冰清預(yù)測(cè)，積雪冰川的融雪徑流預(yù)報(bào)， 水利工程施工期的施工預(yù)報(bào)，水庫運(yùn)行管理要求的入庫水流過程預(yù)報(bào)，河道航運(yùn)要

4、求的 沿程水位變化預(yù)報(bào)等。本文針對(duì)某流域，通過對(duì)其9場(chǎng)洪水資料的審查處理，以及水文預(yù)報(bào)模型參數(shù)的率 定，簡(jiǎn)要編制某流域水文預(yù)報(bào)方案。由于資料數(shù)據(jù)較少，無法進(jìn)行系統(tǒng)和全面的預(yù)報(bào)方 案分析，進(jìn)而結(jié)果也無法得到有效保證，因此本文的目的就是在 9場(chǎng)洪水的基礎(chǔ)上，熟 悉編制水文預(yù)報(bào)方案的具體工作及流程，達(dá)到基本掌握預(yù)報(bào)方案編制目的。1 流域及水文概況1.1 流域水系概況某地處大連市東北部，發(fā)源于鞍山市龍?zhí)多l(xiāng)老北溝，從北向南穿過莊河市中部地區(qū)， 流經(jīng)莊河市和岫山市的十一個(gè)鄉(xiāng)鎮(zhèn)，最后流入黃海。河流全長(zhǎng)94km。某一級(jí)支流6條,依次為沙河、滾子河、漢屯河、二道河、溝門子河、史家卜河。主要水系分東西兩大支 流，

5、東支流為大沙河，匯流點(diǎn)以上積水面積 324km2,西支流為主河流，匯流點(diǎn)以上面 積273km2,與莊河三架山鄉(xiāng)沙河口處會(huì)合。某流域共 7個(gè)雨量站，一處水文站，一處水 庫站，各類遙測(cè)點(diǎn)8處。1.2 下墊面概況某流域?yàn)榈蜕角鹆陞^(qū)。流域的地貌成因和形態(tài)主要受構(gòu)造和巖性控制，進(jìn)而形成了 北高南低的地貌形態(tài)，山脈大體呈北東南西向展布。按其形態(tài)成因類型主要有石英砂巖、 花崗巖、片麻巖及混合巖所組成構(gòu)造侵蝕低山、呈北向或近東西向展布的山間谷地。1.3 洪水特性情況某流域大多數(shù)暴雨集中大多數(shù)暴雨主要集中在 7、8月份，一般7月下旬到8月上 旬，持續(xù)時(shí)間12天，更長(zhǎng)達(dá)3天。由于暴雨強(qiáng)度大、降雨時(shí)間歷時(shí)較短、以及

6、流域面 積不大的特點(diǎn)，使得暴雨造成的洪水具有峰高量大且峰現(xiàn)時(shí)間段的特點(diǎn)。2.2流量資料處理本方案采用9場(chǎng)洪水資料進(jìn)行預(yù)報(bào)方案的編制，9場(chǎng)洪水的降雨資料均由三道溝雨 量站、廟東溝雨量站、大堡子雨量站、大房子雨量站、前梨樹雨量站、姜家堡子雨量站、 某雨量站測(cè)得；流量資料由某水庫日入庫流量資料和時(shí)段入庫流量資料獲得，若無水庫 入庫流量資料，則依據(jù)水庫水位一流量關(guān)系曲線和水量平衡反推入庫流量過程。2.1 降雨資料處理某流域有7個(gè)雨量測(cè)站來記錄降雨資料，這7個(gè)雨量站所控制的面積以及相應(yīng)的權(quán) 重如表2.1。各測(cè)站記錄的各時(shí)刻雨量值乘以相應(yīng)測(cè)站的權(quán)重，再將各測(cè)站該值相加， 即可得到該時(shí)刻的雨量值。表2.1雨

7、量控制面積及控制權(quán)重三道溝廟東溝大堡子大房子前梨樹姜家堡子面積(km2118881077810814053權(quán)重0.170.130.150.110.160.20.08對(duì)9場(chǎng)降雨資料的初步審查中，發(fā)現(xiàn)有些雨量站發(fā)生故障，以及數(shù)據(jù)記錄延遲，造 成降雨量的記錄不是很準(zhǔn)確，因此對(duì)降雨資料的處理，針對(duì)“無問題雨量站降雨處理” 與“有問題雨量站降雨處理”，分兩部分進(jìn)行。(1)無問題雨量站降雨處理20120710號(hào)、20120803號(hào)、20130730號(hào)、20130816號(hào)洪水過程中，7個(gè)雨量站工 作正常。因此，對(duì)各時(shí)段各站的雨量按照表2.1采用面積加權(quán)平均法計(jì)算各時(shí)段的雨量。(2)有問題雨量站降雨處理201

8、10626號(hào)廟東溝站、20130702號(hào)前梨樹站，姜家堡子站、 20130716號(hào)姜家堡子 站發(fā)生故障，故采取舍棄測(cè)站的記錄值，將剩余幾站的面積重新分配權(quán)重計(jì)算各時(shí)段雨 量；20120424號(hào)、20120728號(hào)廟東溝雨量站發(fā)生雨量傳遞延遲，根據(jù)流域降雨同時(shí)刻 其他正常測(cè)站的峰值特征，進(jìn)行該站雨量值移動(dòng)，兩場(chǎng)洪水分別前移2個(gè)時(shí)段、4個(gè)時(shí)段。對(duì)各場(chǎng)洪水，首先進(jìn)行流量過程線平滑處理，由于個(gè)別洪水的流量過程是由水庫水 位一流量關(guān)系曲線和水量平衡反推入庫流量過程求得的，因此存在鋸齒狀的不穩(wěn)定現(xiàn) 象，需要在水量平衡的前提下，對(duì)流量過程線進(jìn)行平滑處理；其次，進(jìn)行次洪徑流分割, 由于各場(chǎng)洪水流量過程中即存在

9、前期洪水的退水過程部分，也存在基流部分，因此本文 采用平均退水曲線求出本次洪水過程，以及運(yùn)用平割法求次洪徑流深。(1) 流量過程線修勻?qū)?場(chǎng)洪水修勻前與修勻后徑流深的對(duì)比見表 2.2,徑流深變化在誤差允許范圍內(nèi)； 并以20110626號(hào)洪水的流量過程線為例進(jìn)行修勻，見圖 2.1。250.0|pi111150.01D0.050.Q黑 200,0-25口3D.D35020110626修勻-40/045.D-50.0表2.2修習(xí)前后徑流深變化表洪號(hào)修習(xí)前修勻后備注2011062629.229.82012042483.182.32012071074.271.6201207285655.11.單位：mm

10、20120803280.2287.92.沒有扣除2013070299.298.1基流20130716155.9155.72013073081.576.12013081670.471.30 10 20 3。 40 50 60 70 a0 90 100 110 120 130 140 150 1&0 170 130 190圖2.1 20110626 修勻流量圖(2)分割本次洪水洪水分割采用退水曲線法。退水曲線是流域蓄水量的消退過程線。對(duì)某一流域而言, 地下徑流退水過程比較穩(wěn)定，故去 9場(chǎng)洪水的洪水過程退水部分，使他們尾部重合，做 光滑的下包線，即可求得流域的地下水退水曲線。運(yùn)用該退水曲線，

11、在本次洪水的流量過程圖中，于本次洪水起漲點(diǎn)處描繪出上一場(chǎng) 洪水的尾水退水部分，然后將該部分扣除，即為本次洪水的流量過程(未扣除基流)。 以20110626號(hào)洪水為例，如圖2.2黑色虛線部分。圖2.2本次洪水流量過程及前期退水20110626降雨修為(3)平割法求次洪徑流深2.3中在各場(chǎng)洪水的起漲點(diǎn)處，做一平行線，于本次洪水退水部分有一交點(diǎn)，則圖ABCDEA的面積即為本次洪水徑流深(未扣除基流)。設(shè)基流為 7m3/s,則扣除基流后 的各場(chǎng)洪水徑流深見表2.3。圖2.3平割法求次洪徑流深圖表2.3各場(chǎng)洪水徑流深(已扣除基流)徑流深2011062624.262012042477.4720120710

12、66.812012072852.891.單位：mm20120803281.812.已扣除基流2013070290.043.20130716只取條-個(gè)洪峰2013071681.02013073066.172013081667.40(4)二水源劃分根據(jù)退水曲線，找出每一場(chǎng)洪水的地面徑流(直接徑流)終止點(diǎn)，與起漲點(diǎn)相連， 即可將實(shí)測(cè)流量過程線進(jìn)行地表徑流 Qs與地下徑流Qg的分割。以20110626號(hào)洪水為 例,如圖2.4。翔slMBInmQ 0 口 o S 115 020.025.030.035.040.045.050.0O1 10 20 30 &0 50 50 7Q 80 90 100

13、110 120 130 HO 150 160 170 180 190圖2.4水源分割圖一水源分割鏤2.3 雨洪關(guān)系的合理性分析對(duì)于一場(chǎng)合格的洪水，必須要滿足兩點(diǎn)要求：(1)符合水量平衡對(duì)于一場(chǎng)洪水，必須符合水量平衡方程式：P=E+R +Wm- Wo(2.1)其中：P:降雨，E:蒸發(fā)，R次洪徑流深，Wm：流域平均蓄水容量，W0:流域初 始土壤含水量。由于在對(duì)洪水資料進(jìn)行模型參數(shù)率定之前，Wm與W0可能未知，因此，對(duì)于一場(chǎng)洪水起碼要符合(2.2)的關(guān)系式：P> E+R(2.2)(2)徑流系數(shù)小于1合格的洪水徑流系數(shù)必須小于1,但不排除各別洪水徑流系數(shù)大與 1的情況，如過 徑流系數(shù)大與1不多也

14、是可以接受的。根據(jù)以上兩個(gè)指標(biāo)對(duì)9洪水進(jìn)行雨洪關(guān)系合理性分析，發(fā)現(xiàn)由3場(chǎng)洪水徑流系數(shù)大 與1。分析原因?yàn)榻涤炅坑涗浿灯』蛘呤橇髁恐灯?。?jīng)過對(duì)某流域水工建筑物及當(dāng) 時(shí)降雨現(xiàn)實(shí)情況的分析，上游沒有攔蓄工程，不會(huì)發(fā)生防水造成流量增大的現(xiàn)象，因此, 只有降雨量記錄上發(fā)生了問題。經(jīng)分析，降雨是由翻斗式雨量計(jì)記錄的，當(dāng)雨強(qiáng)較大時(shí)， 反斗會(huì)因?yàn)榉D(zhuǎn)效率問題造成記錄的雨量比實(shí)際的雨量小。因此，對(duì)問題洪水降雨強(qiáng)度 較大的降雨適當(dāng)增大。處理后發(fā)現(xiàn) 20120803號(hào)洪水徑流系數(shù)為1.1,仍大與1,因此， 將該洪水舍棄，不參與參數(shù)率定工作。結(jié)果如表 2.4。表2.4洪水徑流系數(shù)原始降雨徑流深(修勻后)徑流系數(shù)是

15、否止常20110626106.7024.260.23正常20120424115.8577.470.67正常20120710115.3666.810.58正常2012072867.9552.890.78正常20120803257.33281.811.10否20130702115.6490.040.78正常2013071686.2780.960.94正常2013073065.3666.171.01正常20130816102.2167.400.66正常3預(yù)報(bào)模型參數(shù)率定3.1 模型的選擇某流域?yàn)榘霛駶?rùn)地區(qū)，產(chǎn)流具有蓄滿產(chǎn)流和超滲產(chǎn)流兩種特性。當(dāng)前期土壤較為濕 潤(rùn)時(shí)，容易發(fā)生蓄滿產(chǎn)流，當(dāng)前期比較干旱時(shí)

16、，容易發(fā)生超滲產(chǎn)流。根據(jù)參與參數(shù)率定 的8場(chǎng)洪水的前期降雨及土壤含水量情況，發(fā)現(xiàn)較為濕潤(rùn)，因此適合蓄滿產(chǎn)流模式，故 決定采用二水源新安江模型進(jìn)行產(chǎn)匯流參數(shù)的率定。本文采用的二水源新安江模型包括 產(chǎn)流模型、二水源劃分模型、匯流模型三部分組成（圖 3.1）。（1）產(chǎn)流模型產(chǎn)流模型為蓄滿產(chǎn)流模型，其中土壤蒸發(fā)采用最簡(jiǎn)單的一層蒸發(fā)模式，因此需要率 定三個(gè)參數(shù)：流域平均蓄水容量 Wm；土壤含水量分布曲線指數(shù) b;蒸發(fā)能力折 算系數(shù)k。（2）水源劃分模型將率定出的蓄滿產(chǎn)流參數(shù)帶入產(chǎn)流模型中， 算出各時(shí)刻的徑流量Ri,再運(yùn)用試算法 求穩(wěn)定下滲率fc,進(jìn)而將總凈流量R劃分為地表徑流Rs和地下徑流Rg。（3）匯

17、流模型匯流模型分為地表匯流和地下匯流。地表匯流根據(jù)地表徑流Rs,通過時(shí)段單位線模型可求得地表流量Qs；地下匯流根據(jù)地下徑流 Rg,通過地下線性水庫模型可求得地 下徑流Qg。將Qs與Qg相加即總流量Q,與實(shí)測(cè)流量進(jìn)行對(duì)比，調(diào)參率定。由于時(shí)段 單位線可由瞬時(shí)單位線推求，而瞬時(shí)單位線只需要n （反映流域調(diào)蓄能力的參數(shù)），和K （線性水庫蓄泄系數(shù)），即可確定。并且地下線性水庫只需要Kg （地下水蓄水常數(shù)）即可確定，因此回流模型需要率定 3個(gè)參數(shù)：n；K;Kg。匯 流 模 型n, KKg地表單位線模型時(shí)段單位線UH地表流量Qs地下線性水庫地下徑流演算地下流量Qg輸出總出流Q圖3.1二水源新安江模型計(jì)算流

18、程3.2參數(shù)率定及步驟3.2.1 率定方法選擇本模型采用遺傳算法及人工調(diào)參相結(jié)合的方法進(jìn)行參數(shù)率定。遺傳算法的過程就是從一個(gè)初始種群出發(fā)，經(jīng)過不斷地選擇、交叉、變異，使種群 越來越接近某個(gè)目標(biāo)。如果把種群看作是超空間中的一組點(diǎn)，則選擇、交叉、變異的過 程便是在超空間中點(diǎn)集之間的變換過程，通過這些信息的變換使種群不斷地進(jìn)化，直至 達(dá)到目標(biāo)。標(biāo)準(zhǔn)遺傳算法的基本運(yùn)算過程是一個(gè)典型的迭代過程，其必須完成的工作和基本步 驟如下：選擇編碼策略：根據(jù)實(shí)際問題定義適應(yīng)度函數(shù)；確定遺傳策略，包括種群 大小、選擇、交叉、變異方法，以及確定交叉概率、變異概率等遺傳參數(shù)；隨機(jī)初始 化生成初始群體；計(jì)算當(dāng)前種群中個(gè)體編

19、碼用解碼后的適應(yīng)度；（按照遺傳策略，運(yùn)用選擇、交叉和變異算子作用于群體，形成下一代種群；）判斷種群性能是否滿足某一指標(biāo)，或者己完成預(yù)定迭代次數(shù)滿足則輸出最佳個(gè)體，退出；不滿足則返回。遺傳算法本身具有一定的隨機(jī)性，結(jié)果也具有一定的不確定性，并且，有些參數(shù)采 用試算法率定，因此，需要人為的對(duì)率定出的參數(shù)進(jìn)行調(diào)試，最終得出合理的參數(shù)。3.2.2 參數(shù)及范圍本模型需要的參數(shù)有10個(gè)，其中有7個(gè)需要率定。其取值范圍及性質(zhì)如表 3.1。表3.1模型各參數(shù)匯集表名稱符號(hào)取值范圍敏感性備注土壤初始含水量P W00<常量計(jì)算時(shí)段t1h流域面積F2692km2流域平均蓄水容量Wmw0<Wm<15

20、0mm不敏感算法率定土壤含水量分布曲線指數(shù)b0Vb<0.5不敏感蒸發(fā)能力折算系數(shù)k0.1<k<1敏感穩(wěn)定卜滲舉fc0<fc敏感試算率定地下水蓄水常數(shù)kg40<kg<100敏感算法率定反映流域調(diào)蓄能力的參數(shù)N1<N<2敏感線性水庫蓄泄系數(shù)K0<K<15敏感3.2.3 參數(shù)率定步驟整個(gè)模型分三個(gè)部分進(jìn)行參數(shù)率定， 各參數(shù)率定方式及范圍如表3.1,率定流程（圖3.2） 如下：（1）蓄滿產(chǎn)流模型的 Wm, b, k蓄滿產(chǎn)流模型由參數(shù)率定模型及結(jié)果輸出模型組成。 輸入資料有降雨資料P,蒸發(fā)資料E, 土壤初始含水量 W0 （W0可由前期 影響雨量

21、公式計(jì)算及日降雨資料得到）。經(jīng)參數(shù)率定模型中的遺傳算法循環(huán)率定， 輸出相應(yīng)參數(shù)值。 將率定的結(jié)果帶入結(jié)果輸出模型求出 PE、R時(shí)間序列值。 2）水源劃分模型fc水源劃分模型分為試算fc與求Rs、Rg兩步構(gòu)成水源劃分模型的fc可用手動(dòng)試算法求得。通過輸入上一步求得的PE、R,及Qg（Qg 由洪水分割線分割求得），求得凈流深 R ,與分割出的徑流深R （由Qg求得）進(jìn)行對(duì) 比，通過不斷調(diào)整fc,直到兩者接近獲相等為止，即可獲得 fc。將試算的fc作為參數(shù)帶入到接下來的匯流模型中，即可輸出Rs、Rg。（3）匯流模型kg、N、K匯流模型同樣由參數(shù)率定模型及結(jié)果輸出模型組成。在參數(shù)率定模型輸入實(shí)測(cè)流量過

22、程線 Q以及求出的PE、R、fc,經(jīng)過地表時(shí)段 單位線模型及地下線性水庫模型的計(jì)算及循環(huán)率定后，即可求出相應(yīng)的參數(shù)。 將求出的參數(shù)帶入結(jié)果輸出模型，同樣經(jīng)過地表時(shí)段單位線模型及地下線性水 庫模型的計(jì)算，即可輸出 Rs、Rg、UH、Qs、Qg的時(shí)間序列。PE, R時(shí)間序列蓄 滿 產(chǎn) 流 模 型水 源 劃 分 求 fcI圖3.2參數(shù)率定流程圖匯 流 模 型4預(yù)報(bào)方案本文在 8 場(chǎng)洪水中取 2011062& 20120424 20120710 20120718 20130702、20130716 號(hào)6場(chǎng)洪水進(jìn)行產(chǎn)匯流參數(shù)率定(其中，20120803號(hào)洪水的雨洪關(guān)系不正常，故舍棄不 用。)預(yù)留

23、2013073R 20130816號(hào)2場(chǎng)洪水進(jìn)行參數(shù)檢驗(yàn)。4.1 產(chǎn)流預(yù)報(bào)方案運(yùn)用遺傳算法率定蓄滿產(chǎn)流模型參數(shù)，得出6場(chǎng)洪水的產(chǎn)流參數(shù)如表4.1。表4.1蓄滿產(chǎn)流模型產(chǎn)流模擬參數(shù)洪號(hào)W0 (mm)Wm (mm)bk實(shí)測(cè)徑流深(mm)率定徑流深(mm)(%)20110626431320.210.8124.2623.70-2.302012042437740.40.677.4772.50-6.4220120710781150.230.8566.8171.306.71201207281001140.290.852.8952.80-0.18201307021151390.320.8990.0488.40

24、-1.82201307161251300.20.8280.9670.70-12.68表示率定值小于實(shí)測(cè)值預(yù)留的20130730和20130816號(hào)洪水發(fā)生在7、8月份，因此主要選擇6場(chǎng)洪水中 處于7、8月份的洪水率定的參數(shù)，取參數(shù)的平均值(表 4.2),并置入產(chǎn)流結(jié)果輸出模 型中。將預(yù)留洪水的資料輸入模型，模擬結(jié)果如表 4.3。表4.2輸入產(chǎn)流結(jié)果輸出模型的參數(shù)平均值Wmbk20120710115.00.230.8520120728114.00.290.8020130702139.00.320.8920130716130.00.200.82平均值124.50.260.84表4.3蓄滿產(chǎn)流模型模

25、擬2場(chǎng)洪水的結(jié)果洪號(hào)降雨量(mm)實(shí)測(cè)徑流深(mm)徑流系數(shù)預(yù)報(bào)徑流深(mm)絕對(duì)誤差(mrm相對(duì)誤差(%)是否合格2013073065.466.21.064.2-2.0-2.98是20130816102.267.40.772.95.58.16是根據(jù)水文情報(bào)預(yù)報(bào)規(guī)范 GB/T22482-2008-T，產(chǎn)流誤差范圍需滿足以下其中1項(xiàng)誤差評(píng)定要求即可合格：徑流深相對(duì)誤差20%且徑流深絕對(duì)誤差小于20mm；徑流深絕對(duì)誤差小于3mm。故模擬的兩場(chǎng)洪水都合格，合格率 100%,精度等級(jí)為甲。 因此當(dāng)暴雨時(shí)間為7、8月份時(shí)，可采用表4.4的產(chǎn)流參數(shù)。由于缺少其他月份的洪水資料，故無法驗(yàn)證產(chǎn)流參數(shù)在其他月份

26、的有效性。因此參 數(shù)僅適用于7、8月份。表4.4 7、8月份蓄滿產(chǎn)流參數(shù)Wmbk數(shù)值124.50.260.844.2 匯流預(yù)報(bào)方案運(yùn)用試算法計(jì)算出每一場(chǎng)洪水的穩(wěn)定下滲率，并用遺傳算法對(duì)匯流模型的參數(shù)進(jìn)行率定，得出6場(chǎng)洪水的匯流參數(shù)，如表4.5。根據(jù)水文情報(bào)預(yù)報(bào)規(guī)范GB/T22482-2008-T, 匯流誤差范圍需都滿足3項(xiàng)誤差評(píng)定要求才算合格：降雨徑流預(yù)報(bào)以實(shí)測(cè)洪峰量值的 20%作為許可誤差；峰現(xiàn)時(shí)間以預(yù)報(bào)根據(jù)時(shí)間至實(shí)測(cè)洪峰出現(xiàn)時(shí)間之間時(shí)距的30%作為許可誤差，若許可誤差小于一個(gè)時(shí)段或3h,取一個(gè)時(shí)段或3h;確定性系數(shù)主要根據(jù)精度的等級(jí)要求而定。(本方案設(shè)定為丙級(jí))故率定的匯流參數(shù)合格。表4.

27、5洪水匯流參數(shù)率定結(jié)果洪號(hào)fckgnK實(shí)測(cè)洪峰(m3/s )模擬洪峰 (m3/s)洪峰絕對(duì)誤差 (m3/s)洪峰相對(duì)誤差 (溝峰現(xiàn) 誤差 (h)DC是否合格201106260.3338.291.108.00217216-1-0.4600.74是201204241.1061.551.017.001000984-16-1.6000.89一是201207101.3058.421.2015.0048148761.2500.94是201207287.0051.441.107.50661643-18-2.72-20.81是201307024.60100.001.044.3514001342-58-4.14-20.7是2013071612.00100.001.207.2361662371.1400.91是表示率定值小于實(shí)測(cè)值由某流域的洪水特性情況得知，暴雨造成的洪水具有峰高量大且峰現(xiàn)時(shí)間短的特點(diǎn)。對(duì)6場(chǎng)率定參數(shù)的洪水的流量過程線進(jìn)行觀察后發(fā)現(xiàn)，其洪峰及峰型的確都符合上述特點(diǎn)，并且峰現(xiàn)時(shí)間基本上相似。故可推斷，預(yù)留模擬的20130730和20130816號(hào)洪水同樣具有上述特征。因此，在匯流階段，取 6場(chǎng)洪水率定參數(shù)各項(xiàng)的平均值作為模擬參數(shù)(表4.6)置入?yún)R流結(jié)果輸出模型中。