王新華
(山東省棗莊市水文局 棗莊 277100)
隨機水文模擬技術(shù)在水文預(yù)報中的應(yīng)用
王新華
(山東省棗莊市水文局 棗莊 277100)
水文系統(tǒng)受氣候和人類活動影響,呈現(xiàn)出非常復(fù)雜的行為特征,在現(xiàn)有技術(shù)條件下,還不能用準確的數(shù)理方程描述并求解,要了解水文系統(tǒng)各組成間的相互關(guān)系,分析水文要素的發(fā)展趨勢,當前可行的一類方法就是隨機水文模擬。本文利用隨機水文模擬技術(shù)模擬長系列水文數(shù)據(jù),并將其應(yīng)用于水文預(yù)報,可克服水文預(yù)報中實測系列較短的困難,提高水文預(yù)報精度。
隨機 模擬 預(yù)報 展望 周期
水文預(yù)報是水文系統(tǒng)水情服務(wù)的一項重要內(nèi)容,包括洪水報、徑流預(yù)報、降水量預(yù)報等,降水量展望是降水量預(yù)報的重要組成部分,在此以降水量展望為例談隨機模擬長系列的應(yīng)用。降水量展望就是根據(jù)已知多年降水量系列,點繪降水量—時間關(guān)系線,根據(jù)關(guān)系線反映的降水量發(fā)展趨勢預(yù)測今后幾年甚至幾十年降水量豐、枯狀況。降水量系列一般具有周期性,了解其周期變化,對預(yù)測精度有很大幫助,但由于實測系列資料較短,很難反映其周期變化,因此利用水文隨機模擬技術(shù)模擬水文長系列數(shù)據(jù)是非常必要的。
1.1隨機水文現(xiàn)象
現(xiàn)實世界存在兩點現(xiàn)象:確定現(xiàn)象和隨機現(xiàn)象。水文過程顯示著水文現(xiàn)象隨時間而變化的特性,由于水文現(xiàn)象受眾多因素的影響,所以水文過程常呈現(xiàn)出隨機變化性。在全面研究隨機水文變化的基礎(chǔ)上建立起反映水文現(xiàn)象主要變化特性的隨機水文模型,就可模擬與實際水文數(shù)據(jù)統(tǒng)計特性相近的大量水文序列,用于統(tǒng)計、預(yù)測,滿足水利工程設(shè)計、水文預(yù)報對長系列資料的需求,能夠更加合理、可靠地表征未來水文現(xiàn)象可能出現(xiàn)的各種情況。
1.2均勻隨機數(shù)模擬
理論上,只要有了任何一個連續(xù)分布隨機數(shù),就可模擬出其他任意分布隨機數(shù)。由于(0,1)區(qū)間上的均勻分布是最簡單、最基本的連續(xù)分布,所以在此使用(0,1)區(qū)間上的均勻分布隨機數(shù),其產(chǎn)生方法是利用計算機標準函數(shù),見式(1)。
式中:X—(0,1)區(qū)間隨機數(shù);RND()—標準函數(shù)。
1.3服從P-Ⅲ分布純隨機序列模擬
1.3.1服從正態(tài)分布純隨機序列模擬
對均勻隨機數(shù)作以下變換:
式中,ξ1,ξ2相互獨立的標準正態(tài)分布〔N(0,1)〕變量,這樣由均勻隨機數(shù)Xt即可模擬出標準正態(tài)變量ξt。
1.3.2服從P-Ⅲ分布純隨機序列模擬
利用標準化隨機數(shù)ξt和標準化P-Ⅲ分布隨機數(shù)φt之間的關(guān)系(4式)、實測系列的Cs便可模擬出標準P-Ⅲ分布隨機數(shù)φt,再由已知系列參數(shù)平均值、Cv和(5式)模擬出P-Ⅲ分布隨機系列。
圖1 隨機模擬年降水量系列滑動均值過程線圖
圖2 隨機模擬降水量系列滑動均值累積距平值過程線圖
表1 實測年降水量系列與隨機模擬年降水量系列統(tǒng)計參數(shù)對比表
2.1年降水量周期性分布
根據(jù)隨機模型和實測年降水量系列統(tǒng)計參數(shù)平均值、Cv、CS,可以模擬300年年降水量系列,該系列具有周期性變化特征,可將其看作一種波動,是多種振動合成的結(jié)果,設(shè)想一種濾波器,年降水量系列輸入后,有的波被過濾掉,有的波被保留,或者有的波被削弱,有的波更加明顯,可以采用這種濾掉短波,保留長波和超長波的濾波器來分析降水量周期。
2.1.1平滑濾波器周期分析對于上述系列X(t),平滑濾波器模型為:
式中:k—步長,k=1,2,……,n;
i—序號,i=1,2,……,n-k+1。選定合適步長k,依次求出滑動均值,就可得到最合適滑動均值過程線,使周期反映最清晰。取步長5年,繪制出模擬年降水量系列滑動均值過程線(見圖1),可以看出模擬降水量系列城豐枯交替變化,豐枯變化周期為15年左右居多,最大不超過20年,最小為8年。經(jīng)滑動濾波器濾波后,其周期變化可以作為短期洪水預(yù)報參考,提高短期預(yù)報精度。為更好地分析系列中長期變化周期及規(guī)律,采用累積濾波器對系列進行第二次濾波。
2.1.2累積濾波器周期分析
對于以上滑動均值系列X(t),累積濾波器模型為:
式中:n—序列中數(shù)據(jù)總數(shù)。
滑動均值系列過濾方法:求出系列均值,然后求系列中每一項與均值之差即“距平值”,再計算累計距平值,就可繪制出年降水量系列累積距平值過程線(見圖2),如曲線上兩點連線坡度為正,表示這一時期為豐水期;反之,為枯水期;豐、枯水組合基本上就經(jīng)歷了一個周期。從過程線上可以比較直觀的看出,模擬降水量系列有共有三個大的周期變化,年降水量系列周期80年左右。
2.2模擬系列代表性分析
模擬系列畢竟不是實測數(shù)據(jù),所以有必要分析一下它對實測數(shù)據(jù)的代表程度,是否能反映年降水量豐、枯水組合情況,以及系列與實測系列統(tǒng)計參數(shù)符合程度。選用嶧城水文站1956~2015年實測年降水量系列與隨機模擬年降水量系列進行代表性分析,統(tǒng)計參數(shù)見表1。
由表1可得出模擬年降水量系列統(tǒng)計參數(shù)與實測系列較符合,具有較好的代表性,模擬系列可應(yīng)用于降水中、長期展望。
隨機水文模擬技術(shù)是在水文計算基礎(chǔ)上發(fā)展起來的先進技術(shù),利用隨機模擬技術(shù)模擬水文長系列,用于中長期水文預(yù)報,能夠提高預(yù)報精度。同時,模擬長系列解決了現(xiàn)行水文計算中實測系列較短的難題,減少計算中的主觀性,能夠有效提高設(shè)計暴雨、設(shè)計洪水等計算精度。因此,隨機水文模擬技術(shù)可在水文預(yù)報和水文計算中得到較好的應(yīng)用,具有很好的研究價值和發(fā)展前景■