梅 超，尹明萬，洪 林，李 蒙，李東琴(.武漢大學(xué)水資源與水電工程國家重點(diǎn)實(shí)驗(yàn)室，武漢 4007；.中國水利水電科學(xué)研究院水資源研究所 北京 0008；.華北水利水電大學(xué)， 鄭州 4500)

通過隨機(jī)水文模擬獲得大量徑流序列，再進(jìn)行水文系統(tǒng)分析計(jì)算，以認(rèn)識(shí)系統(tǒng)各組成要素間的復(fù)雜關(guān)系，是規(guī)劃、設(shè)計(jì)、運(yùn)行與管理復(fù)雜水資源系統(tǒng)的重要方式，在水庫群優(yōu)化調(diào)度、水資源系統(tǒng)風(fēng)險(xiǎn)分析和徑流預(yù)報(bào)等水資源規(guī)劃與管理實(shí)踐中具有廣泛應(yīng)用[1-6]。隨機(jī)水文模擬的關(guān)鍵是建立合理可靠的隨機(jī)水文模型，目前常用的隨機(jī)水文模型可以分為回歸類、解集類和具有物理基礎(chǔ)類3種類型[7]。

月徑流等季節(jié)性非平穩(wěn)水文序列隨機(jī)模擬最常用的隨機(jī)水文模型是季節(jié)性自回歸模型SAR(p)，該模型考慮了非平穩(wěn)隨機(jī)序列的明顯季節(jié)變化特征，根據(jù)不同季節(jié)性序列的統(tǒng)計(jì)特征分別建立各自的自回歸模型[8]。季節(jié)性隨機(jī)水文模型能夠較好地保持實(shí)測(cè)樣本季節(jié)序列的各項(xiàng)統(tǒng)計(jì)特征，但由于此類模型將本應(yīng)為有機(jī)統(tǒng)一整體的年內(nèi)分配過程視為相互獨(dú)立的隨機(jī)過程，使得模擬的結(jié)果在保持季節(jié)徑流特性的前提下常不能同時(shí)較好的保持年徑流序列統(tǒng)計(jì)特性，在應(yīng)用中受到限制[9]。年徑流時(shí)間序列可視為平穩(wěn)隨機(jī)過程，對(duì)年徑流的隨機(jī)模擬已有許多精度較高可以應(yīng)用的隨機(jī)水文模型，先隨機(jī)模擬得到年徑流，然后采用一定的方法將年徑流分解得到月徑流，是模擬月徑流的一種重要方法[9]。本文建立了基于年徑流自回歸隨機(jī)模擬和雙層模型分解法、典型解集分解法的月徑流隨機(jī)模型，以平寨水庫為例，對(duì)入庫徑流進(jìn)行隨機(jī)模擬，并比較了兩種分解方法的適用性。

1 基于年徑流分解的月徑流隨機(jī)模擬方法

1.1 年徑流隨機(jī)模擬模型

年徑流隨機(jī)模擬技術(shù)成熟、模擬精度較高，常用于年徑流隨機(jī)模擬的模型有ARMA(p,q)、AR(p)、最鄰近抽樣模型和基于小波分析的組合隨機(jī)模型等[10-12]。本文對(duì)年徑流的隨機(jī)模擬采用自回歸模型AR(p)，并基于AIC準(zhǔn)則確定AR(p)模型的階p，有關(guān)AIC準(zhǔn)則的具體內(nèi)容可以參見文獻(xiàn)[9]。AR(p)模型的基本形式為：

xt=u+φ1(xt-1-u)+φ2(xt-2-u)+…+

φp(xt-p-u)+εt

(1)

式中：xt為原始年徑流隨機(jī)序列；u為xt的均值；φ為自回歸系數(shù)；p為自回歸階數(shù)；εt為均值為0、方差為δ2t的獨(dú)立隨機(jī)變量，該方差與樣本方差有一定聯(lián)系，對(duì)于偏態(tài)序列常用的處理方法有對(duì)數(shù)轉(zhuǎn)換法、獨(dú)立隨機(jī)項(xiàng)變換法和W-H變換法[9]。

1.2 月徑流隨機(jī)模擬模型

由年徑流隨機(jī)水文模型模擬得到年徑流后，在保證年徑流統(tǒng)計(jì)特性不變的前提下，可以采用一定的方法將年徑流分解成月徑流。常采用的年徑流分解方法包括雙層模型法、典型解集法和相關(guān)解集法，不同分解方法有各自的優(yōu)缺點(diǎn)，在實(shí)際應(yīng)用過程中應(yīng)根據(jù)實(shí)際需要進(jìn)行比選后再確定分解方法。

本文采用雙層模型法和典型解集法對(duì)年徑進(jìn)行分解得到月徑流，并對(duì)兩種方法得到的分解結(jié)果進(jìn)行比較分析，模擬流程如圖1所示。

圖1 月徑流隨機(jī)模擬流程

雙層模型法的主要流程是：隨機(jī)模擬年徑流；采用季節(jié)性隨機(jī)模型如SAR(p)等模擬月徑流；基于年徑流對(duì)月徑流進(jìn)行修正，修正公式為:

(2)

式中：Qi為隨機(jī)模擬得到的年徑流；qi,j由季節(jié)性隨機(jī)模型模擬得到的第i年第j(j=1，2,…,12)月的月徑流；q′'i,j為調(diào)整后的第i年第j月徑流量。

典型解集法的主要流程是：隨機(jī)模擬年徑流；計(jì)算實(shí)測(cè)徑流的月徑流分配系數(shù)；在實(shí)測(cè)年徑流尋找與隨機(jī)模擬年徑流最為接近的年份，將年徑流模擬值按該年的月分配系數(shù)解集為月徑流。

2 實(shí)例應(yīng)用

2.1 研究區(qū)概況

平寨水庫位于三岔河干流木底河河段，壩址位于貴州省惠水縣太陽鄉(xiāng)平寨，距離貴陽直線距離85 km，是一座以灌溉、防洪、發(fā)電為主的大型水庫。水庫上游集雨面積3 492 km2,正常水位1 331 m，死水位1 305 m，總庫容10.89 億m3，調(diào)節(jié)庫容4.48 億m3。本文選取平寨水庫40年(1968-2007)年入庫徑流量為基礎(chǔ)數(shù)據(jù)，實(shí)測(cè)徑流數(shù)據(jù)的統(tǒng)計(jì)特性如表1所示，4項(xiàng)統(tǒng)計(jì)特征參數(shù)分別為均值Ut、均方差δt、變差系數(shù)Cv、和偏態(tài)系數(shù)Cs。

表1 實(shí)測(cè)徑流數(shù)據(jù)的統(tǒng)計(jì)特性

2.2 模擬結(jié)果的統(tǒng)計(jì)特性分析

采用AR(p)模型模擬年徑流，由AIC準(zhǔn)則確定p為1，即年徑流采用一階自回歸模型AR(1)模擬，由于年徑流序列呈偏態(tài)分布，采用獨(dú)立隨機(jī)項(xiàng)變換法建立偏態(tài)AR(1)模型。由AR(1)模擬得到1 050組年徑流，為排除初值的影響，舍棄前50組，采用余下1 000組年徑流模擬序列分別采用雙層模型法和典型解集法分解成月徑流。

采用季節(jié)性一階自回歸模型SAR(1)模擬得到1 050組12個(gè)月徑流序列，舍棄前50組，對(duì)余下1 000組模擬值的統(tǒng)計(jì)參數(shù)進(jìn)行分析，表2為采用SAR(1)模擬得到的各項(xiàng)統(tǒng)計(jì)參數(shù)的相對(duì)誤差，該表中年徑流序列通過累加相應(yīng)SAR(1)模擬得到的12個(gè)月月徑流序列取得。

基于SAR(1)模擬得到的月徑流和AR(1)模擬得到年徑流，采用雙層模型法由式(2)修正得到月徑流各項(xiàng)參數(shù)統(tǒng)計(jì)誤差如表3中方法1所示；基于AR(1)模擬得到的年徑流采用典型解集法分解得到月徑流各項(xiàng)參數(shù)統(tǒng)計(jì)誤差如表3中方法2所示。

AR(1)直接模擬得到的年徑流序列各項(xiàng)統(tǒng)計(jì)參數(shù)的相對(duì)誤差分別為-1.36%、-3.59%、-4.81%、-1.26%，四項(xiàng)主要統(tǒng)計(jì)參數(shù)相對(duì)誤差均較小，說明由AR(1)隨機(jī)模擬得到的年徑流序列能夠較好的保持樣本的統(tǒng)計(jì)特性。由SAR(1)模擬各月及由各月徑流相應(yīng)統(tǒng)計(jì)得到年徑流各項(xiàng)統(tǒng)計(jì)參數(shù)相對(duì)誤差如表2所示，其中年徑流的均值、均方差、變差系數(shù)、偏態(tài)系數(shù)的相對(duì)誤差分別為-1.10%、-22.45%、21、59%、-77.35%，除均值模擬效果較好以外，其他3項(xiàng)參數(shù)的誤差均在-20%以上，說明由SAR(1)模擬得到月徑流再統(tǒng)計(jì)得到的年徑流序列不能較好的保持樣本的統(tǒng)計(jì)特性。

表2 SAR(1)模型模擬結(jié)果統(tǒng)計(jì)參數(shù)的相對(duì)誤差 %

對(duì)比表3中雙層模型法和典型解集法模擬得到各月徑流序列各項(xiàng)統(tǒng)計(jì)參數(shù)的相對(duì)誤差，典型解集法所得結(jié)果的相對(duì)誤差普遍較雙層模型法小，且兩種方法對(duì)均值、均方差和變差系數(shù)的模擬效果要好于偏態(tài)系數(shù)的模擬效果，說明在年徑流在解集為月徑流的過程中兩種方法均未能較好的保持月徑流序列的偏態(tài)特性。對(duì)比表2和表3中分解方法1相應(yīng)的各月各項(xiàng)參數(shù)的相對(duì)誤差，可以發(fā)現(xiàn)，表3中相應(yīng)參數(shù)的相對(duì)誤差相對(duì)于表2均有不同程度的增大，說明由SAR(1)模擬得到的月徑流序列經(jīng)雙層模型法用式(2)修正后，在保持年徑流序列統(tǒng)計(jì)特性的同時(shí)，月徑流序列模擬結(jié)果的精度有一定的喪失。其中，雙層模型修正對(duì)均值的統(tǒng)計(jì)特性保持最好，經(jīng)修正后的月徑流序列均值大相對(duì)誤差為-5.79%；而對(duì)其他3項(xiàng)參數(shù)的統(tǒng)計(jì)特性保持最差，各月均方差、變差系數(shù)、偏態(tài)系數(shù)的最大相對(duì)誤差分別達(dá)53.9%、50.22%、910.30%。

表3 兩種分解方法模擬結(jié)果統(tǒng)計(jì)參數(shù)的相對(duì)誤差 %

2.3 模擬結(jié)果年內(nèi)分配過程分析

對(duì)于以月為時(shí)間尺度的水資源系統(tǒng)設(shè)計(jì)、規(guī)劃、調(diào)度管理等水資源分析計(jì)算而言，徑流的年內(nèi)分配同樣具有重要意義，同樣的年徑流量，不同的年內(nèi)分配將導(dǎo)致不同的水資源優(yōu)化配置和調(diào)度決策。

選取10%、25%、50%、75%、95%五種頻率和多年平均條件下雙層模型法和典型解集法兩種方法模擬結(jié)果的年內(nèi)分配過程，將其與實(shí)測(cè)的同頻率下的年內(nèi)分配過程進(jìn)行對(duì)比分析，結(jié)果如圖2所示。

對(duì)比分析圖2中雙層模型法模擬結(jié)果(圖2中過程線1)的年內(nèi)分配過程與實(shí)測(cè)年內(nèi)分配過程，在頻率為10%和25%的偏豐年份模擬結(jié)果與實(shí)測(cè)過程差別較大，而在50%、75%、95%等平水年份和偏枯年份，則基本上能夠保持與實(shí)測(cè)年內(nèi)分配過程大致趨勢(shì)相同；多年平均條件下，模擬結(jié)果的年內(nèi)分配過程與實(shí)測(cè)過程均基本完全一致。

對(duì)比分析圖2中典型解集法模擬結(jié)果(圖2中過程線2)的年內(nèi)分配過程與實(shí)測(cè)年內(nèi)分配過程，10%、25%兩個(gè)偏豐年份下，模擬結(jié)果與實(shí)測(cè)過程不完全一致，但基本能夠保持相同的趨勢(shì)；在50%、75%、95%和多年平均條件下，典型解集法模擬結(jié)果與實(shí)測(cè)過程基本完全一致。

綜合對(duì)比兩種方模擬結(jié)果在5種來水頻率和多年平均條件下年內(nèi)分配過程與實(shí)測(cè)過程的符合性，在偏豐水年份(10%、25%)，兩種方法模擬結(jié)果與實(shí)測(cè)過程的符合性均較差，但典型解集方法模擬結(jié)果比雙層模擬法模擬符合更好；在平水年份(50%)和偏枯年份(75%、95%)，雙層模型法模擬結(jié)果的年內(nèi)分配過程與實(shí)測(cè)過程的符合性則明顯的較典型解集法模擬結(jié)果差；在多年平均條件下兩種方法模擬結(jié)果與實(shí)測(cè)過程的符合性均較好。

進(jìn)一步對(duì)比年內(nèi)各月兩種模擬結(jié)果與實(shí)測(cè)徑流的符合性，可以發(fā)現(xiàn)，模擬結(jié)果與實(shí)測(cè)結(jié)果差別較大的月份主要集中在5-9月，即平寨水庫所處流域的汛期。

表4對(duì)比了兩種模擬結(jié)果汛期(5-9月)徑流量及其占年徑流總量的比例與實(shí)測(cè)過程的符合性，對(duì)比分析可知，雙層模型法在各頻率模擬結(jié)果汛期徑流量占年徑流量的比例均與實(shí)測(cè)過程有一定差別，典型解集法在10%和25%頻率下模擬結(jié)果汛期徑流量占年徑流量的比例與實(shí)測(cè)比例有所差別，但差別較小，且均比雙層模型法的差別小，在50%、75%、95%等頻率下，典型解集法模擬結(jié)果汛期徑流量占年徑流量的比例與實(shí)測(cè)比例完全相符。在多年平均條件下，兩種方法模擬結(jié)果徑流量占年徑流量的比例均能較好符合實(shí)測(cè)比例。

圖2 兩種分解方法模擬結(jié)果不同頻率下的年內(nèi)分配過程

表4 兩種方法模擬結(jié)果汛期徑流量及占年徑流量的比例

綜合分析可知，兩種模擬結(jié)果年內(nèi)分配過程與實(shí)測(cè)過程符合性隨著來水頻率的升高而變好，即在偏豐年符合較差，而在偏枯年符合較好；豐水年份水庫入庫徑流受降雨的影響較大，而降雨的不確定性較大，因此在對(duì)豐水年份月徑流過程進(jìn)行模擬時(shí)不能較好的符合實(shí)測(cè)過程，偏枯年則與之相反；兩種方法在多年平均條件均下的年內(nèi)分配過程與實(shí)測(cè)過程基本完全一致，說明兩種方法在整體上均保持了樣本序列的基本特性。模擬結(jié)果與實(shí)測(cè)過程符合較差的月份主要出現(xiàn)在汛期(5-9月)，由于汛期徑流受降雨影響較大，具有較強(qiáng)的不確定性，因此對(duì)汛期的模擬結(jié)與實(shí)測(cè)過程的符合程度較非汛期低。

3 結(jié) 語

采用AR(1)模型隨機(jī)模擬得到平寨水庫年徑流，分別采用雙層模型法和典型解集法將隨機(jī)模擬的年徑流分解為月徑流，對(duì)年徑流模擬結(jié)果和月徑流分解結(jié)果進(jìn)行了分析。

(1)采用SAR(1)模擬平寨水庫月徑流，并由相應(yīng)月徑流統(tǒng)計(jì)得到年徑流，該結(jié)果能較好地保持樣本的月徑流量的統(tǒng)計(jì)特性，但不能同時(shí)較好的保持年徑流量的統(tǒng)計(jì)特性；由雙層模型法將由AR(1)模型模擬得到的年徑流分解成月徑流，年徑流各項(xiàng)統(tǒng)計(jì)參數(shù)保持AR(1)模型的模擬精度不變，但由雙層模型法修正分解后的月徑流與原SAR(1)模型模擬得到的月徑流相比，各項(xiàng)統(tǒng)計(jì)參數(shù)的精度均有不同程度下降，其中偏態(tài)系數(shù)的精度下降最為顯著。

(2)典型解集法分解得到的各月徑流均值、均方差和變差系數(shù)等三項(xiàng)統(tǒng)計(jì)參數(shù)相對(duì)誤差大部分較雙層模型法模擬結(jié)果相應(yīng)相對(duì)誤差小，其中典型解集法模擬結(jié)果均值、均方差和變差系數(shù)的最大相對(duì)誤差分別為8.30%、16.85%和13.81%；兩種方法模擬結(jié)果偏態(tài)系數(shù)的相對(duì)誤差均較大。

(3)兩種分解方法在偏豐水年的年內(nèi)分配過程與實(shí)測(cè)過程的符合性較差，在偏枯水年符合較好，在多年平均條件下均基本與實(shí)測(cè)過程完全一致，總體上不同頻率下典型解集法模擬結(jié)果的年內(nèi)分配過程與實(shí)測(cè)過程的符合性均優(yōu)于雙層模型法模擬結(jié)果；兩種分解方法分解的月徑流結(jié)果年內(nèi)分配與實(shí)測(cè)過程符合較差的月份均主要集中在汛期(5-9月)。

□

[1] 劉 攀,郭生練,張文選,等. 梯級(jí)水庫群聯(lián)合優(yōu)化調(diào)度函數(shù)研究[J]. 水科學(xué)進(jìn)展,2007,(6):816-822.

[2] 羅啟華,周研來,楊 娜,等. 基于SAR模型和蒙特卡羅法的防洪調(diào)度風(fēng)險(xiǎn)分析[J]. 人民長江,2011,(1):4-8.

[3] 郭旭寧,胡鐵松,曾 祥,等. 基于調(diào)度規(guī)則的水庫群供水能力與風(fēng)險(xiǎn)分析[J]. 水利學(xué)報(bào),2013,(6):664-672.

[4] 曹枝俏,王國利,梁國華,等. 基于隨機(jī)模擬信息的神經(jīng)網(wǎng)絡(luò)洪水預(yù)報(bào)模型[J]. 水力發(fā)電學(xué)報(bào),2010,(4):63-69.

[5] Gaivoronski A, Sechi G M, Zuddas P. Cost/risk balanced management of scarce resources using stochastic programming[J]. European Journal of Operational Research, 2012,216:214-224.

[6] Celeste, A B. & Billib, M. Evaluation of stochastic reservoir operation optimization models[J]. Advances in Water Resources, 2009, 32:1 429-1 443.

[7] 王文圣,金菊良,李躍清. 水文隨機(jī)模擬進(jìn)展[J]. 水科學(xué)進(jìn)展,2007,(5):768-775.

[8] 武蘭婷,何士華,朱劍軍,等. 季節(jié)性自回歸模型適用性分析[J]. 中國農(nóng)村水利水電,2014,(12):17-19,23.

[9] 王文圣，丁 晶，金菊良.隨機(jī)水文學(xué)[M].2版.北京:中國水利水電出版社, 2008:152-156.

[10] 趙麗娜. 年徑流序列隨機(jī)模擬與預(yù)報(bào)研究[D]. 陜西楊凌：西北農(nóng)林科技大學(xué),2010.

[11] 褚永軍,李 明,原 彪,等 永定河上游清水河年徑流序列特性分析及其隨機(jī)模擬[J]. 水文,2011,(2):35-37,45.

[12] 衡 彤,王文圣,李拉丁,等. 基于小波變換的組合隨機(jī)模型及其在徑流隨機(jī)模擬中的應(yīng)用[J]. 水電能源科學(xué),2002,(1):15-17.