郭曉娜,崔春梅,魯素芬
(1.山東省水利勘測(cè)設(shè)計(jì)院有限公司,山東 濟(jì)南 250014;2.山東省水利科學(xué)研究院,山東 濟(jì)南 250014;3.山東省水文中心,山東 濟(jì)南 250002)
根據(jù)全省地質(zhì)巖性情況和一般山丘區(qū)推理公式研究有關(guān)條件要求,經(jīng)綜合分析,選取沂河上游朱家莊流域作為山東省一般山丘區(qū)典型研究流域。朱家莊站符合選用站條件:流域面積較小,產(chǎn)匯流條件比較一致,河道干流比降大且沿程變化不大,流域內(nèi)雨量站數(shù)基本能夠控制流域面降水量。朱家莊小流域?qū)儆谔锴f水庫(kù)流域,位于沂河上游。朱家莊水文站位于高莊河(即田莊河)上游,流域形狀為扇形。流域內(nèi)共有4處雨量站,分別為朱家莊、曹家莊、磨山石和雙廟雨量站。流域內(nèi)無(wú)蒸發(fā)監(jiān)測(cè)資料,移用附近鎮(zhèn)后站蒸發(fā)資料。
基于典型流域產(chǎn)匯流規(guī)律研究獲取的損失參數(shù)μ與主時(shí)段雨強(qiáng)I等新參數(shù)[1],分別采用瞬時(shí)單位線法、推理公式法和地區(qū)經(jīng)驗(yàn)公式法推求小流域相應(yīng)頻率的設(shè)計(jì)洪峰流量。
由洪水水文要素摘錄表分析,朱家莊小流域洪水過(guò)程線多呈陡漲陡落,設(shè)計(jì)暴雨的最長(zhǎng)統(tǒng)計(jì)時(shí)段采用1日。由朱家莊、曹家莊、磨山石和雙廟4雨量站歷年降水量摘錄表,經(jīng)過(guò)雨量資料等時(shí)距插補(bǔ)、站點(diǎn)權(quán)重計(jì)算等前期工作,統(tǒng)計(jì)計(jì)算歷年最大1日面雨量。根據(jù)歷年最大1日面雨量系列,進(jìn)行頻率適線,得朱家莊站200年一遇最大1日面暴雨量379.7 mm,100年一遇最大1日面暴雨量332.9 mm,50年一遇最大1日面暴雨量286.2 mm,20年一遇最大1日面暴雨量225.0 mm,10年一遇最大1日面暴雨量179.2 mm,5年一遇最大1日面暴雨量133.9 mm。
1)設(shè)計(jì)凈雨計(jì)算。查閱山東省水文圖集,確認(rèn)朱家莊流域設(shè)計(jì)前期影響雨量Pa為40 mm。根據(jù)朱家莊產(chǎn)流規(guī)律研究建立的降水徑流關(guān)系線圖[2],分別查算出不同保證率的最大1日面降水對(duì)應(yīng)的設(shè)計(jì)凈雨。
2)設(shè)計(jì)凈雨時(shí)程分配。朱家莊小流域位于山東省沂河流域上游,絕大部分面積為山丘區(qū)。因此,采用水文圖集中沂沭汶泗區(qū)的設(shè)計(jì)雨型進(jìn)行凈雨時(shí)程分配。設(shè)計(jì)凈雨時(shí)段為0.5 h。
根據(jù)朱家莊流域匯流模型研究成果[2],朱家莊流域瞬時(shí)單位線匯流參數(shù)m1=2.3,m2=1.43。因此,求得n=0.70,k=3.29。通過(guò)S曲線轉(zhuǎn)換[3],求得各時(shí)段凈雨的單位線,基于不同保證率的0.5 h時(shí)段設(shè)計(jì)凈雨與相應(yīng)單位線,求得地表徑流洪水過(guò)程,考慮地下徑流過(guò)程的回加,綜合得到朱家莊流域設(shè)計(jì)洪水過(guò)程。朱家莊流域不同頻率設(shè)計(jì)洪水成果見表1。
表1 朱家莊流域?qū)崪y(cè)暴雨法設(shè)計(jì)洪峰流量成果表
1)設(shè)計(jì)暴雨時(shí)程分配?;谏綎|省水文圖集沂沭泗區(qū)概化時(shí)程雨型[2],根據(jù)設(shè)計(jì)面暴雨計(jì)算結(jié)果,對(duì)朱家莊流域不同頻率設(shè)計(jì)暴雨進(jìn)行時(shí)程分配。
2)設(shè)計(jì)凈雨計(jì)算。首先根據(jù)暴雨時(shí)程分配,確定不同頻率降水主時(shí)段雨強(qiáng)。其次由朱家莊流域損失參數(shù)μ與主時(shí)段雨強(qiáng)I關(guān)系[1],推求不同頻率暴雨主時(shí)段雨強(qiáng)對(duì)應(yīng)的損失參數(shù)μ。
根據(jù)設(shè)計(jì)暴雨時(shí)程分配,自雨強(qiáng)最大的雨峰中心開始,向前后相鄰時(shí)段連續(xù)累積,求得不同時(shí)段的最大降雨量Pt,繪制設(shè)計(jì)暴雨P(guān)t~t曲線[4]。在同一坐標(biāo)系中同時(shí)繪制損失曲線μt,平移μt與Pt~t曲線相切,求得切點(diǎn)對(duì)應(yīng)的t即為產(chǎn)流歷時(shí)tc,切線與Pt軸交點(diǎn)即為產(chǎn)流量hR,如圖1所示。
圖1 切線法求解產(chǎn)流量hR示意圖
3)洪峰流量推求。根據(jù)小流域推理公式研究,在參數(shù)m、σ、λ和流域特征值F、L、J已知條件下,采用圖解法[4]推求洪峰流量Qm。在同一坐標(biāo)紙上,Q~τ與Q~t曲線的交點(diǎn),即為所求的Qm值,由此,得到朱家莊流域不同頻率的設(shè)計(jì)洪峰流量,見表2。
表2 朱家莊流域推理公式法設(shè)計(jì)洪峰流量成果表
1)計(jì)算成峰暴雨量。由朱家莊流域最大24 h設(shè)計(jì)暴雨時(shí)程分配,繪制Pt~t曲線(方法同上);采用參數(shù)綜合確定的流域匯流時(shí)間τ[2](3.58 h);在此基礎(chǔ)上,由曲線Pt~t查得相應(yīng)Pt,即為成峰暴雨量。
2)設(shè)計(jì)洪峰流量推求。采用小流域經(jīng)驗(yàn)公式研究成果,分別采用多次相關(guān)圖解公式、SPSS多元線性回歸公式進(jìn)行設(shè)計(jì)洪峰推求。洪峰流量計(jì)算成果見表3。
表3 朱家莊流域地區(qū)經(jīng)驗(yàn)公式設(shè)計(jì)洪峰流量推求m3/s
根據(jù)朱家莊站洪水要素摘錄表,統(tǒng)計(jì)歷年最大洪峰流量。根據(jù)歷年最大流量系列,考慮建站前的兩次大洪水(1922年、1957年)調(diào)查成果(358 m3/s、244 m3/s),進(jìn)行頻率適線。分析得朱家莊水文站200年一遇洪峰流量425.0 m3/s,100年一遇洪峰流量371.2 m3/s,50年一遇洪峰流量317.2 m3/s,20年一遇洪峰流量245.7 m3/s,10年一遇洪峰流量191.4 m3/s,5年一遇洪峰流量136.8 m3/s。
根據(jù)以上研究成果,將朱家莊流域設(shè)計(jì)暴雨推求洪峰成果與流量推求洪峰成果進(jìn)行對(duì)比分析,見表4。
從表4中看出,除極個(gè)別設(shè)計(jì)洪峰誤差較大外,大部分設(shè)計(jì)暴雨推求的洪峰流量與流量推求洪峰流量相對(duì)誤差都控制在±35%以內(nèi)。同時(shí),朱家莊流域設(shè)計(jì)暴雨推求的洪峰流量,三種方法計(jì)算結(jié)果互相對(duì)比,各設(shè)計(jì)頻率成果相差不大。究采用的暴雨洪水系列(1976—2012年),老圖集的所采用的系列(1952—1970年)要偏豐許多,據(jù)此分析綜合的參數(shù)所推求的成果也要偏大很多。
表4 朱家莊站設(shè)計(jì)洪峰流量成果對(duì)比
基于山東省水文圖集已有成果,分別采用瞬時(shí)單位線、推理公式和地區(qū)經(jīng)驗(yàn)公式三種方法推求不同保證率的設(shè)計(jì)洪峰流量,并將設(shè)計(jì)暴雨推求洪峰流量成果與水文圖集成果進(jìn)行對(duì)比分析,見表5。
表5 朱家莊站新、老成果計(jì)算洪峰流量對(duì)比
從表5中可以看出,不管采用哪種計(jì)算方法,設(shè)計(jì)暴雨推求洪峰流量成果與查閱山東省水文圖集老成果相比,設(shè)計(jì)洪峰差值較大。究其原因,這種較大差距是由新老成果采用的不同雨洪系列資料和產(chǎn)匯流參數(shù)決定的。相較于新成果研
基于朱家莊小流域產(chǎn)匯流研究參數(shù)新成果,本研究采用設(shè)計(jì)暴雨推求洪水方法(含瞬時(shí)單位線法、推理公式法、地區(qū)經(jīng)驗(yàn)公式法)分別推求了流域控制斷面保證率為20%、10%、5%、2%、1%和0.5%的設(shè)計(jì)洪峰流量,進(jìn)行了驗(yàn)證。結(jié)果表明,采用新參數(shù)推求的洪峰流量成果,瞬時(shí)單位線法、推理公式法與地區(qū)經(jīng)驗(yàn)公式三種方法相互對(duì)比,各設(shè)計(jì)頻率成果相差不大;與實(shí)測(cè)流量推求洪水成果相比,除極個(gè)別設(shè)計(jì)洪峰誤差較大以外,大部分設(shè)計(jì)暴雨推求的設(shè)計(jì)洪峰流量與流量推求設(shè)計(jì)洪峰流量相對(duì)誤差是相對(duì)合理的;而與查閱水文圖集老成果相比,設(shè)計(jì)洪峰流量普遍偏小很多,這主要是由老水文圖集參數(shù)研究采用的雨洪資料偏豐,而且系列代表性較差造成的。