范雪芬

（遼寧省石佛寺水庫管理局，遼寧 沈陽 110166）

1 流域介紹

石佛寺水庫位于遼河干流的中下游，是遼河干流唯一的控制性水利工程，遼寧省內(nèi)流域面積16 161萬km2，水庫總庫容1.85億m3。石佛寺水庫通過與遼河干流左側(cè)三大支流上的南城子、清河、柴河、榛子嶺4座水庫聯(lián)合調(diào)控，將石佛寺水庫以下遼河干流的防洪標(biāo)準(zhǔn)從30年一遇提高到100年一遇[1]。

根據(jù)石佛寺流域特點和水利工程分布情況，將石佛寺流域劃分為4個子區(qū)間，即松樹子區(qū)間、開原子區(qū)間、鐵嶺子區(qū)間和石佛寺子區(qū)間[2,3]。石佛寺水庫洪水預(yù)報方案中，產(chǎn)匯流預(yù)報采用大伙房模型[4]，洪水演算采用錯時段法。

2 分類洪水預(yù)報模型

2.1 模型建立

采用模糊聚類ISODATA迭代模型對洪水進行聚類，在文獻[5]中詳細介紹了ISODATA模型的推導(dǎo)過程，將歷史洪水劃分為大洪水、中等洪水和小洪水3類，以洪峰和洪量作為模糊聚類的特征指標(biāo)，指標(biāo)權(quán)重向量W=（0.5，0.5）。由于無法提前獲取洪峰和洪量兩個指標(biāo)，使得基于洪峰和洪量指標(biāo)的模糊聚類方法應(yīng)用受到了極大限制。使用模糊推理的方法，來提前判斷洪水所屬類型，以便快速準(zhǔn)確的選取相應(yīng)模型參數(shù)。

基于綜合信息的模糊推理模式是建立在以成因分析、統(tǒng)計分析、經(jīng)驗相關(guān)等方法的基礎(chǔ)上，利用歷史信息歸納生成推理模式語句，再根據(jù)面臨時刻信息與歷史信息相匹配，得到相應(yīng)的結(jié)果作為結(jié)論。模糊推理要求具有代表性好，且足夠多的歷史信息作為輸入因子用來生成規(guī)則。

根據(jù)多規(guī)則、多輸入、單輸出的模型求解思路，多重模糊推理形式為：

式中：A為推理因子；B為推理結(jié)果；m為推理因子個數(shù)；n為規(guī)則個數(shù)。

洪水類型識別模糊推理模型，分為離線階段和在線階段兩部分，洪水分類預(yù)報流程如圖1所示。離線階段是根據(jù)歷史資料生成推理規(guī)則和將歷史洪水進行聚類，在線階段是利用面臨時刻信息進行模糊推理，判斷實時洪水類型。

圖1 基于模糊推理的洪水分類預(yù)報模型流程圖

2.2 模型求解

由于松樹子單元水文資料好，人類活動少，且下墊面條件穩(wěn)定，適合進行模糊推理研究。對松樹子單元1954-2010年，47場歷史洪水作為研究資料進行分析，選擇前期土壤含水量（WM0）、降雨強度（i）、起始流量（Q0）、累計降雨量（P），4個因子作為模型的輸入條件。

具體的計算過程是，將四個因子分級，如表1所示。將歷史洪水的特征指標(biāo)轉(zhuǎn)化為相應(yīng)的級別作為樣本，而獲取實時洪水的特征指標(biāo)后，將特征指標(biāo)轉(zhuǎn)化為相應(yīng)的級別，與樣本洪水中各因子的級別相比較，認為級別最相近的為同一類洪水，采用同一組預(yù)報模型參數(shù)。

表1 模糊推理因子分級

將歷史洪水劃分為率定期和檢驗期，將率定期和檢驗期的分類預(yù)報結(jié)果和不分類預(yù)報結(jié)果列入表2和表3中。經(jīng)過對比分析可以看出，分類預(yù)報方案較原方案相比，模擬精度和合格率都有所提高，預(yù)報誤差減小[6]。用模糊聚類方法能夠根據(jù)特征指標(biāo)準(zhǔn)確的將歷史場次洪水進行聚類，并根據(jù)聚類結(jié)果對模型參數(shù)進行分類優(yōu)選，將不同類型的參數(shù)應(yīng)用于洪水預(yù)報中，模擬精度得到了明顯的提高。

表2 松樹子單元率定期洪水預(yù)報結(jié)果比較

表3 松樹子單元檢驗期洪水預(yù)報結(jié)果比較

3 結(jié)論

為了提高石佛寺流域洪水預(yù)報精度，建立了基于綜合信息的模糊推理模式洪水分類預(yù)報模型，并以松樹子單元為例進行研究。采用洪量和洪峰作為指標(biāo)，應(yīng)用模糊聚類ISODATA模型將歷史場次洪水進行聚類，進行不同類型洪水的模型參數(shù)優(yōu)選。由于實時洪水預(yù)報中無法提前獲取洪峰、洪量指標(biāo)，故選擇前期土壤含水量、降雨強度、起始流量、累計降雨量等，4個因子作為模糊推理的輸入條件，判斷該場洪水所屬的洪水類型。實際結(jié)果表明，洪峰和洪量指標(biāo)與前期土壤含水量、降雨強度、起始流量和累計降雨量等，4個因子存在良好的相關(guān)關(guān)系，并且分類洪水預(yù)報可以有效提高洪水預(yù)報精度。

［1］梁文章，柴曉利，曹煒倫.石佛寺水庫2010720暴雨洪水分析［J］.東北水利水電，2011，624（3）∶49-50.

［2］王猛，彭勇，梁國華.大伙房模型在石佛寺流域洪水預(yù)報方案中的應(yīng)用［J］.南水北調(diào)與水利科技，2012，10（2）∶15-19.

［3］王猛.石佛寺水庫洪水預(yù)報及其實時修正方法研究［D］.大連理工大學(xué)，2012.

［4］大連理工大學(xué)，國家防汛抗旱總指揮部辦公室.水庫防洪預(yù)報方法與應(yīng)用［M］.北京∶中國水利水電出版社.

［5］陳守煜.工程模糊集理論與應(yīng)用［M］.北京∶國防工業(yè)出版社，1998.

［6］中華人民共和國水利部.水文情報預(yù)報規(guī)范（SL-2000）［S］.北京，2000.