磨盤山水庫(kù)中長(zhǎng)期入庫(kù)徑流模糊預(yù)測(cè)模型研究

李 鵬

(哈爾濱市供水工程有限責(zé)任公司，哈爾濱150010)

1 徑流中長(zhǎng)期預(yù)測(cè)方法選擇

徑流預(yù)測(cè)是水庫(kù)興利調(diào)度中非常重要的部分，現(xiàn)階段基本分大氣與非大氣因子分析法、數(shù)理統(tǒng)計(jì)法和模糊數(shù)學(xué)法3大類。目前國(guó)內(nèi)研究比較多的方法和模型主要有模糊推理法、灰色系統(tǒng)預(yù)測(cè)方法、周期疊加外推法、平穩(wěn)時(shí)間序列外推法、人工神經(jīng)網(wǎng)絡(luò)模型和遺傳算法等。論文分析選擇單要素模糊模式推理法和平穩(wěn)時(shí)間序列法2種方法對(duì)水庫(kù)壩址以上流域入庫(kù)徑流過(guò)程進(jìn)行預(yù)測(cè)。

單要素模糊模式推理法方法簡(jiǎn)單，具有資料收集容易，預(yù)見(jiàn)期長(zhǎng)，尤其是能分級(jí)預(yù)報(bào)的優(yōu)點(diǎn)，符合徑流預(yù)測(cè)作業(yè)的實(shí)際情況;并且，隨著資料的積累，推理模式的子式可逐漸增加，克服了數(shù)理統(tǒng)計(jì)法中因補(bǔ)充新的資料需重新分析計(jì)算的缺點(diǎn)。關(guān)于平穩(wěn)時(shí)間序列外推法，該方法假定水文要素的時(shí)序演變過(guò)程是一個(gè)隨機(jī)過(guò)程現(xiàn)實(shí)，并且屬于平穩(wěn)隨機(jī)過(guò)程。一般水文要素通常屬于非平穩(wěn)時(shí)間序列，但是考慮到:①要與單要素模糊模式法的預(yù)報(bào)結(jié)果進(jìn)行對(duì)照;②磨盤山水庫(kù)流域特性較好，徑流年際變化不是太大。因此，選擇了該方法。

2 單要素模糊模式推理法

2.1 單要素模糊模式推理方法原理

單要素模糊模式推理法利用水文或氣象要素自身前m年資料作預(yù)報(bào)因子(輸入)，通過(guò)模糊模式推理模式(模糊控制器)，以特征展開(kāi)近似推理，求出待預(yù)報(bào)的第m+1年要素值(輸出)。本方法的基本原理是:

1)統(tǒng)計(jì)待預(yù)報(bào)的水文或氣象要素x(如年徑流量，汛期流量等)的歷史資料，其最大、最小值(或通過(guò)相似流域水文比擬法適當(dāng)外延)作為論域的域界，U［xmin，xmax］。

2)將論域上U上的值離散化，并分成l級(jí)(如l=5級(jí)，則分別表示豐、偏豐、正常、偏枯、枯)，每級(jí)對(duì)應(yīng)一個(gè)模糊子集，{ Yk}?U(x)，k=1，2，…，l。模糊子集的隸屬函數(shù)分布亦可取簡(jiǎn)單的梯形分布，以Yk(x)公式表示。

如果有n年歷史演變資料，則每個(gè)值都對(duì)應(yīng)論域U上的一個(gè)模糊子集，即

3)設(shè)想前期連續(xù)m年流量與第m+1年年徑流量存在模糊關(guān)系。利用前期連續(xù)m年模糊信息(亦稱預(yù)報(bào)因子)作輸入，亦模糊推理模式作控制器，通過(guò)變換輸出第m+1年模糊量(亦稱預(yù)報(bào)量)，便可根據(jù)n年資料歸納出(n-m)條m重復(fù)合模糊條件語(yǔ)句，即

上式就是具有m重(n-m)次復(fù)合模糊蘊(yùn)涵命題特性的“單要素預(yù)報(bào)的模糊推理模式”，即模糊控制器。它描述了前m個(gè)要素(亦可比擬為預(yù)報(bào)因子)模糊集與第m+1個(gè)要素(待預(yù)報(bào)量)模糊集之間的模糊關(guān)系，可表示為:

式中:代表第 i個(gè)模式，第 j個(gè)因子為第 k級(jí)模糊子集，?U(x)，i=1，2，…，n-m，j=1，2，…，m+1;k=1，2，…1。根據(jù)預(yù)留年份前m年實(shí)測(cè)年徑流總量來(lái)預(yù)報(bào)第m+1年年徑流總量，然后，將預(yù)報(bào)值與實(shí)際值進(jìn)行比較，分析誤差，如果誤差較大，調(diào)整m，返回第(3)步，重新建立模糊控制器，最確定精度最高的m值。

5)根據(jù)確定的m值，將所有年份徑流資料全部統(tǒng)計(jì)在內(nèi)建立模糊控制器，將未來(lái)年份前m年年徑流量資料做輸入，預(yù)報(bào)未來(lái)年年徑流量。

2.2 磨盤山水庫(kù)單要素模糊模式推理模型

1)確定流域有1952—2003年共52 a資料(2001—2003年徑流資料由電子文檔資料得到)，見(jiàn)表1。由歷史資料可以得到年最大、最小年徑流量為:902.99×106m3和281×106m3，對(duì)年徑流量最大最小值的適當(dāng)外延作為論域的域界U［200，936］(論域單位:106m3)。

2)級(jí)別劃分，將論域分為枯、偏枯、平、偏豐、豐5種等級(jí)，當(dāng)年徑流總量＜400×106m3時(shí)視為枯水年;當(dāng)年徑流總量在400×106m3與540×106m3時(shí)視為偏枯年;當(dāng)年徑流總量在540×106m3和760×106m3之間時(shí)為平水年;在760×106m3和800×106m3之間時(shí)為偏豐年，年徑流總量＞800×106m3時(shí)為豐水年。

3)確定系列周期，系列周期的確定影響模型的精度，有經(jīng)驗(yàn)的預(yù)報(bào)人員可以根據(jù)以往的預(yù)報(bào)經(jīng)驗(yàn)來(lái)取定，本處采用3—11周期分別試預(yù)報(bào)，選擇預(yù)報(bào)精度較高的5 a周期作為周期，記為m(也可根據(jù)系列資料繪制歷年演變過(guò)程線，或使用方差分析、方差普密度圖等方法分析周期)。

4)建立模型，根據(jù)1952—2003年52 a歷史系列資料和五級(jí)分級(jí)標(biāo)準(zhǔn)歸納出推理模式(大前提，記為R)，周期為3的部分推理模式語(yǔ)句見(jiàn)表2。

表1 磨盤山水庫(kù)歷年年徑流量表 106m3

2.3 推理預(yù)報(bào)

以周期為3 h預(yù)報(bào)1994年徑流量為例，可由1994年前3 a的年徑流量資料確定小前提，將小前提輸入到大前提R中，推理過(guò)程如下:

2.3.1 趨勢(shì)一致選擇

以D1表示趨勢(shì)值，D1越大表示小前提與該子句變化趨勢(shì)越接近。用 C i，t(i=1，2，3，…，n-mk;t=1，2，3，…，mk)表示小前提與大前提第i子句中級(jí)別差，則 Ci，t，特別的當(dāng)級(jí)別相同，即=0時(shí)ci，t=1;;用 di，t(i=1，2，3，…，n-mk;t=1，2，3，…，mk)表示小前提中與大前提第i子句中級(jí)別差系數(shù)，級(jí)別相同 di，t=1 ，差一級(jí) di，t=0.5 ，差兩級(jí)或兩級(jí)以上di，t=0，則趨勢(shì)一致性Di可以表示為:

如果Di＞0，說(shuō)明大前提中第i子句與待預(yù)報(bào)年份前mk年水文要數(shù)變化趨勢(shì)基本相似，第i子句的推理結(jié)果可能用來(lái)預(yù)報(bào)。

推理模式中與小前提具有趨勢(shì)一致性子句見(jiàn)表2。

表2 推理模式表

2.3.2 相似性選擇

相似性選擇是在表2中列出的Di＞0子句中選出最大Di，將其輸出信息加權(quán)疊加得出預(yù)報(bào)結(jié)果=［1 × 3+3 × 1+5 ×1］/［1+1+1］ =3，即

即得出1994年的預(yù)報(bào)級(jí)別屬于三級(jí)。根據(jù)前面的分級(jí)標(biāo)準(zhǔn)，1994年來(lái)水量在540×106～760×106m3，為平水年，取區(qū)間的中點(diǎn)值作為預(yù)報(bào)值，得到1994年徑流量650×106m3。該方法可行。

3 結(jié)語(yǔ)

由于中長(zhǎng)期徑流方法有很多種，以上我僅分析2種預(yù)測(cè)方法，建議廣大水利工作者在實(shí)際預(yù)報(bào)過(guò)程中采取多種方法預(yù)報(bào)，針對(duì)自己的水庫(kù)積累預(yù)報(bào)經(jīng)驗(yàn)，選擇適合流域的預(yù)報(bào)模型，以提高徑流預(yù)報(bào)精度。