李福威，包愛美，疏杏勝，丁 偉

（1.國電電力和禹水電開發(fā)公司，遼寧本溪 117201；2.大連理工大學(xué)水利工程學(xué)院，遼寧大連 116024）

0 引言

長期徑流預(yù)報(bào)對(duì)掌握未來徑流信息，實(shí)現(xiàn)水資源的高效利用、水電站優(yōu)化運(yùn)行、防汛抗旱等具有重要意義。然而，由于長期預(yù)報(bào)的預(yù)見期較長，其受大氣環(huán)流、下墊面情況、人類活動(dòng)等多重因素影響［1］，徑流過程的不確定性大，預(yù)報(bào)難度較大，預(yù)報(bào)精度難以保證。目前，對(duì)于降低較長預(yù)見期徑流預(yù)報(bào)的不確定性、提高預(yù)報(bào)精度的研究，主要考慮兩個(gè)方面，一是預(yù)報(bào)因子篩選方面，探究徑流形成的物理基礎(chǔ)和影響徑流變化的各因素及其間的相互作用，從眾多因子中篩選出與預(yù)報(bào)對(duì)象密切相關(guān)的因子；二是預(yù)報(bào)模型方面，研究預(yù)報(bào)模型的方法原理及其適用性，進(jìn)行多模型的比較分析和綜合利用。

隨著中長期徑流預(yù)報(bào)研究的快速發(fā)展，可利用預(yù)報(bào)方法較多，主要分為兩類，一類是根據(jù)徑流序列自身規(guī)律，僅考慮徑流單要素作為預(yù)報(bào)因子的歷史演變法、時(shí)間序列法等傳統(tǒng)統(tǒng)計(jì)方法［2-4］；一類是挖掘水文和氣象數(shù)據(jù)的潛在規(guī)律，考慮水文、氣象等多要素作為預(yù)報(bào)因子的回歸分析法和機(jī)器學(xué)習(xí)法。后者可選用一切與預(yù)報(bào)對(duì)象存在潛在關(guān)系的影響因素［5］，包括考慮入海平面氣壓、大氣環(huán)流因子等氣象因子用作預(yù)報(bào)因子來實(shí)現(xiàn)對(duì)流域徑流的預(yù)報(bào)，具有較強(qiáng)的物理基礎(chǔ)，獲得較長的預(yù)見期和相對(duì)較高的預(yù)報(bào)精度［6-8］。朱春苗等［9］利用相關(guān)系數(shù)法、主成分分析法和互信息法優(yōu)選松花江流域上游3個(gè)水文代表站的預(yù)報(bào)因子，采用SVR（Support Vector Regression）模型進(jìn)行流域月徑流預(yù)報(bào)，得出適用于該流域月徑流預(yù)報(bào)的最優(yōu)預(yù)報(bào)因子組合模型。雷莉和王超［10］從130項(xiàng)氣候系統(tǒng)指數(shù)中篩選出預(yù)報(bào)因子后，基于篩選的因子構(gòu)建了石羊河流域BP-ANN（Back Propaga‐tion-Artificial Neural Network）、Elman 和PSO-SVR（Particle Swarm Optimization-SVR）三種年徑流預(yù)報(bào)模型并比較各模型預(yù)報(bào)結(jié)果，確定Elman 和PSO-SVR 中長期徑流預(yù)報(bào)模型可為該流域中長期徑流預(yù)報(bào)提供支撐。李伶杰等［11］以龍江水庫入庫徑流預(yù)報(bào)為研究對(duì)象，從環(huán)流指數(shù)、海溫、氣壓和前期月徑流中選取關(guān)鍵預(yù)報(bào)因子，建立隨機(jī)森林與支持向量機(jī)模型，發(fā)現(xiàn)太平洋中北部與西部氣候因子對(duì)徑流預(yù)報(bào)的影響較大。

考慮到長期徑流預(yù)報(bào)受影響因素多，預(yù)報(bào)不確定性較大的問題，本研究以渾江桓仁水庫年入庫徑流和汛期月入庫徑流預(yù)報(bào)為研究對(duì)象，考慮將前期降雨、徑流作為預(yù)報(bào)因子的同時(shí)引入大氣環(huán)流因子，建立包括統(tǒng)計(jì)分析法和機(jī)器學(xué)習(xí)法的多個(gè)長期徑流預(yù)報(bào)模型，比較各模型在桓仁流域的適用性，并分析各模型對(duì)年徑流預(yù)報(bào)及汛期各月徑流預(yù)報(bào)的預(yù)報(bào)水平，給出最優(yōu)預(yù)報(bào)方案，提高桓仁流域長期徑流預(yù)報(bào)的精度。

1 材料與方法

1.1 研究區(qū)概況及數(shù)據(jù)資料

桓仁水庫作為渾江梯級(jí)電站的龍頭水庫，總庫容為36.4 億m3，電站裝機(jī)容量為222.5 MW，兼有防洪、灌溉等綜合作用?；溉仕畮煲陨狭饔蛴兄行∷畮旃?7 座，總庫容合計(jì)2.43 億m3，其中中型水庫6 座，總庫容1.96 億m3，占比80.7%，小型水庫41座，庫容0.48 億m3，占比19.3%。流域多年平均年降水量860 mm，降水年內(nèi)分配不均，70%的雨量集中在6-9 月間，且在7、8月達(dá)到最大，大洪水主要發(fā)生在7 月下旬至8 月中旬；冬季一般從11 月份開始到翌年3 月末或4 月初結(jié)束，期間主要以降雪為主，受積雪影響，冬季徑流一般較枯。因此，桓仁水庫來水從5月份流域發(fā)生春汛開始一直持續(xù)到10月份，汛期作為桓仁水庫年內(nèi)主要來水階段，充分把握汛期各月來水情況，對(duì)于電站做出全局較優(yōu)的發(fā)電調(diào)度過程，提升發(fā)電效益具有重要意義；而年徑流預(yù)報(bào)可為水庫運(yùn)行管理提供更長遠(yuǎn)的參考信息，幫助管理者更早地制定調(diào)度計(jì)劃，從而更高效地利用水力資源。

徑流資料來源于桓仁水庫1967-2017 年月實(shí)測(cè)入庫流量，1967-2017 年74 項(xiàng)大氣環(huán)流因子的逐月數(shù)據(jù)源自國家氣候中心網(wǎng)站（http：∕∕www.ncc-cma.net∕cn∕），桓仁水庫流域1967-2017年逐月實(shí)測(cè)降雨量源自中國氣象數(shù)據(jù)網(wǎng)（http：∕∕data.cma.cn∕）。

1.2 研究方法

1.2.1 預(yù)報(bào)因子篩選

（1）因子初選。采用相關(guān)系數(shù)法對(duì)預(yù)報(bào)因子進(jìn)行初選，按照相關(guān)性大小分別從諸多水文、氣象因子選取預(yù)報(bào)因子。由于預(yù)報(bào)因子對(duì)徑流的影響存在滯后效應(yīng)，預(yù)報(bào)年徑流時(shí)，水文因子取桓仁水庫預(yù)報(bào)年份前3 年的實(shí)測(cè)年降雨量與年徑流量，氣象因子取預(yù)報(bào)時(shí)刻前一年內(nèi)國家氣候中心各月份發(fā)布的74 項(xiàng)大氣環(huán)流因子。預(yù)報(bào)月徑流時(shí)，前12 個(gè)月的水文、氣象因子作為輸入因子。相關(guān)性系數(shù)計(jì)算公式為：

式中：r為相關(guān)系數(shù)；n為資料樣本數(shù)；Xi與分別為輸入因子及其序列均值；Yi與分別為預(yù)報(bào)對(duì)象及其序列均值。r的取值范圍為［-1，1］，其正負(fù)號(hào)表示兩變量之間是正相關(guān)還是負(fù)相關(guān)，絕對(duì)值越接近1表明兩變量的線性相關(guān)程度越高。

（2）因子精選。在使用相關(guān)系數(shù)法初選出因子后，進(jìn)一步結(jié)合預(yù)報(bào)模型分析不同因子組合對(duì)模擬和檢驗(yàn)結(jié)果的影響，以確定最終的預(yù)報(bào)因子。本文采用向前搜索包裹法，以率定期殘差平方和最小為目標(biāo)，通過評(píng)價(jià)預(yù)報(bào)模型在不同因子組合情況下的預(yù)報(bào)結(jié)果，確定最優(yōu)預(yù)報(bào)因子組合，具體步驟如下：

①將初選因子相關(guān)性系數(shù)排序第一的因子x0作為固定因子，分別與其他因子組合輸入到預(yù)報(bào)模型，選擇率定期的殘差平方和f(x0，xi)最小時(shí)加入的因子xi作為選定因子。

②以x0、xi為固定因子，分別與其他因子組合輸入到預(yù)報(bào)模型，確定使f(x0，xi，xj)最小時(shí)的因子xj，若f(x0，xi，xj)小 于f(x0，xi)，則將xj因子加入模型，進(jìn)行步驟③；否則，排除該因子。

③重復(fù)步驟②，不斷增加輸入因子個(gè)數(shù)，直至f不再減小或者f的遞減速度小于一定程度時(shí)，從而確定出最終預(yù)報(bào)因子集。

1.2.2 預(yù)報(bào)模型構(gòu)建

選用多元線性回歸（Multiple Linear Regression，MLR）［12］、多元門限回歸（TR）［13］、逐步回歸（Stepwise Regression，SR）［14］、人工神經(jīng)網(wǎng)絡(luò)（Back Propagation-Artificial Neural Network，BPANN）［15］、支持向量機(jī)（Support Vector Machine，SVM）［16，17］和基于主成分分析的人工神經(jīng)網(wǎng)絡(luò)模型（Artificial Neural Network Model based on Principal Component Analysis，PCA-BP-ANN）六個(gè)模型，以最終篩選的水文、氣象因子為模型輸入，以預(yù)報(bào)徑流的合格率（QR）為評(píng)定指標(biāo)，分別構(gòu)建桓仁水庫流域年徑流預(yù)報(bào)模型和汛期月徑流預(yù)報(bào)模型。據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范GB∕T 22482-2008》，許可誤差限為多年同期變幅的20%，桓仁水庫年徑流多年變幅的20%為14.9 億m3，汛期各月的允許誤差（已轉(zhuǎn)化為流量）見表1。

表1 汛期各月允許預(yù)報(bào)誤差 m3∕sTab.1 Allowable forecast error of each month in the flood season

多元線性回歸（MLR）、多元門限回歸（TR）、逐步回歸（SR）等傳統(tǒng)統(tǒng)計(jì)模型在中長期徑流預(yù)報(bào)中應(yīng)用較早，其結(jié)構(gòu)簡單、易于操作，只需確定各預(yù)報(bào)模型數(shù)學(xué)方程，即可得到較好的預(yù)報(bào)效果。人工神經(jīng)網(wǎng)絡(luò)（BP-ANN）和支持向量機(jī)（SVM）作為在中長期徑流預(yù)報(bào)中應(yīng)用最廣泛且較成熟的機(jī)器學(xué)習(xí)模型，具有較強(qiáng)的非線性映射能力，相對(duì)傳統(tǒng)的統(tǒng)計(jì)方法，能夠更好地刻畫徑流過程非線性及非穩(wěn)態(tài)性的特點(diǎn)。選用具有3層網(wǎng)絡(luò)結(jié)構(gòu)（輸入層、隱含層和輸出層）的人工神經(jīng)網(wǎng)絡(luò)（BP-ANN）模型，以年徑流預(yù)報(bào)合格率為指標(biāo)，采用試錯(cuò)法確定網(wǎng)絡(luò)隱含層，并經(jīng)反復(fù)訓(xùn)練確定每層結(jié)構(gòu)所含神經(jīng)元的個(gè)數(shù)；SVM 模型能夠快速處理小樣本的非線性問題，具有較強(qiáng)的泛化能力，模型需確定的參數(shù)有懲罰函數(shù)C、不敏感損失系數(shù)ε、核函數(shù)及其所含參數(shù)g，本文采用試錯(cuò)法確定核函數(shù)，并采用網(wǎng)格搜索法優(yōu)化確定懲罰系數(shù)C和核函數(shù)參數(shù)g。

主成分分析（Principal Component Analysis，PCA）法是通過解析原來具有一定相關(guān)性的多變量，將原變量進(jìn)行線性組合導(dǎo)出一組新的綜合變量，使這些新的綜合變量能夠盡可能多地反映原始變量的信息，以達(dá)到簡化數(shù)據(jù)和降維的目的［18］。本文構(gòu)建基于主成分分析的人工神經(jīng)網(wǎng)絡(luò)模型（PCA-BP-ANN）時(shí)，采用PCA 法將篩選后的預(yù)報(bào)因子重新組合為新的變量，并根據(jù)計(jì)算的特征值、方差貢獻(xiàn)率和累積貢獻(xiàn)率確定主成分個(gè)數(shù)，以新確定的主成分作為輸入因子輸入到BP-ANN 模型進(jìn)行徑流預(yù)測(cè)。

2 結(jié)果與分析

2.1 預(yù)報(bào)因子篩選

（1）因子初選?？紤]水文和氣象兩類因子對(duì)徑流變化的影響，采用相關(guān)性分析法分別對(duì)流域水文、氣象因子進(jìn)行初選，取相關(guān)性系數(shù)較大的因子作為備選因子。對(duì)桓仁水庫流域月徑流量與前期各月份實(shí)測(cè)降雨、徑流和大氣環(huán)流因子進(jìn)行相關(guān)性分析，得到該流域長期徑流預(yù)報(bào)的備選因子，結(jié)果見表2 和表3。由表可知，對(duì)于月徑流，除了上月月降雨量相關(guān)性系數(shù)達(dá)0.57，其他水文因子與預(yù)報(bào)對(duì)象的相關(guān)性都低于備選的大氣環(huán)流因子；對(duì)于預(yù)見期更長的年徑流，雖然大氣環(huán)流因子與年徑流的相關(guān)性系數(shù)相對(duì)月徑流有所減小，但其與預(yù)報(bào)對(duì)象的相關(guān)性依舊遠(yuǎn)高于水文要素。

表2 相關(guān)性較高的水文因子Tab.2 Hydrological factors with high correlation

（2）因子精選。結(jié)合所選預(yù)報(bào)模型，采用向前搜索包裹法，對(duì)表2 和表3 確定的備選因子進(jìn)一步篩選，分別確定各預(yù)報(bào)模型的最優(yōu)組合因子作為各模型的最終輸入因子，表4 列出了結(jié)合多元線性回歸預(yù)報(bào)模型確定的最終預(yù)報(bào)因子，以此作為多元線性回歸預(yù)報(bào)模型的輸入，同理可確定其他預(yù)報(bào)模型的最終預(yù)報(bào)因子。

表3 相關(guān)性較高的大氣環(huán)流因子Tab.3 Atmospheric circulation factors with high correlation

采用主成分分析（PCA）法將表4 中選定的預(yù)報(bào)因子進(jìn)行線性組合，使其成為一組新的線性無關(guān)的綜合變量，并根據(jù)特征值的累積貢獻(xiàn)率達(dá)90%確定主成分個(gè)數(shù)，各成分的特征值及其方差貢獻(xiàn)率和累積貢獻(xiàn)率結(jié)果見表5。由表5 可知，對(duì)于汛期月徑流預(yù)報(bào)，前3 個(gè)主成分的累積貢獻(xiàn)率為90.5%，而后3 個(gè)主成分的貢獻(xiàn)率較低，對(duì)汛期月徑流預(yù)報(bào)的影響不大，因此選定前3 個(gè)主成分為模型輸入，式（2）～（4）分別給出了3 個(gè)主成分與標(biāo)準(zhǔn)化原始變量的關(guān)系；對(duì)于年徑流預(yù)報(bào)，前5個(gè)主成分的累積貢獻(xiàn)率才能達(dá)90%以上，即需要5 個(gè)主成分方可有效地代表原始變量90%以上的信息。此外，年徑流預(yù)報(bào)確定的主成分個(gè)數(shù)比月徑流預(yù)報(bào)多，說明預(yù)報(bào)預(yù)見期更長的年徑流所受影響因素要比預(yù)報(bào)月徑流更為復(fù)雜，其預(yù)報(bào)不確定性更大。

表4 多元線性回歸預(yù)報(bào)模型最終預(yù)報(bào)因子Tab.4 Determined predictors

表5 應(yīng)用PCA法計(jì)算的各成分特征值、方差貢獻(xiàn)率和累積貢獻(xiàn)率%Tab.5 The eigenvalues，variance contribution rate and cumulative contribution rate of each component calculated by PCA method

式中：Fi（i=1，2，3）表示主成分；Vj（j=1，…，6）表示標(biāo)準(zhǔn)化原始變量。

2.2 模型構(gòu)建及預(yù)報(bào)效果分析

將桓仁水庫流域?qū)崪y(cè)徑流資料以及選中的預(yù)報(bào)因子劃分為率定期和驗(yàn)證期兩部分，率定期為1967-2000 年，驗(yàn)證期為2001-2017 年，以預(yù)報(bào)徑流的合格率（QR）為評(píng)定指標(biāo)，分別構(gòu)建多個(gè)桓仁水庫流域長期徑流預(yù)報(bào)模型并驗(yàn)證。式（5）和式（6）分別給出了流域汛期月徑流和年徑流的多元線性回歸模型（MLR）預(yù)報(bào)方程，而多元門限回歸（TR）和逐步回歸（SR）方法類似；表6 展示了人工神經(jīng)網(wǎng)絡(luò)（BP-ANN）和支持向量機(jī)（SVM）兩機(jī)器學(xué)習(xí)模型的參數(shù)結(jié)果；采用主成分分析（PCA）結(jié)果，將月預(yù)報(bào)確定的3 個(gè)主成分和年預(yù)報(bào)確定的5 個(gè)主成分分別輸入到BP-ANN模型中，構(gòu)建月徑流預(yù)報(bào)PCA-BP-ANN 模型和年徑流預(yù)報(bào)PCA-BP-ANN模型。

表6 BP-ANN和SVM模型參數(shù)Tab.6 The parameters of BP-ANN and SVM models

式（5）和式（6）中xi分別為該模型最終確定的月徑流預(yù)報(bào)輸入因子和年徑流預(yù)報(bào)輸入因子（見表4）。

圖1給出了各月徑流預(yù)報(bào)模型在率定期和驗(yàn)證期的模擬預(yù)報(bào)精度，由圖可知，無論是在率定期還是在驗(yàn)證期，不同模型對(duì)不同月份的模擬預(yù)報(bào)精度均有一定差異。在率定期，TR模型在8、9 月份的模擬預(yù)報(bào)效果最好，合格率達(dá)90%，而在5、6、7 月份的表現(xiàn)卻不如SVM模型；雖然SVM和BP-ANN機(jī)器學(xué)習(xí)模型的合格率在5、6 月份的合格率明顯優(yōu)于傳統(tǒng)回歸模型，但在8、9、10 月份卻不如回歸模型；結(jié)合主成分分析法的PCA-BP-ANN模型在6月份表現(xiàn)最差，但在7月和9月份的合格率要比沒有采用PCA 方法的BP-ANN 模型高10%左右。在驗(yàn)證期，TR 模型在5月和10月份預(yù)報(bào)效果最差，預(yù)報(bào)合格率不足60%，而在7、8月份卻能達(dá)85%以上，其中7 月份合格率達(dá)94%；PCA-BPANN模型雖在6、7月份預(yù)報(bào)合格率不足50%，但在5月份和8月份表現(xiàn)最佳，合格率均在85%以上，且該模型相對(duì)BP-ANN 模型精度提升了10%左右，說明結(jié)合主成分分析法（PCA）進(jìn)行該流域月徑流預(yù)報(bào)可有效提高汛期某些月份的預(yù)報(bào)精度。對(duì)比驗(yàn)證期和率定期，各模型在驗(yàn)證期的預(yù)報(bào)精度相對(duì)率定期均有所降低，而SVM 和BP-ANN 等機(jī)器學(xué)習(xí)模型在汛期各月表現(xiàn)相對(duì)穩(wěn)定。

圖1 汛期月預(yù)報(bào)各模型在率定期和驗(yàn)證期的模擬精度Fig.1 The simulation accuracy of monthly forecast models in flood season for both calibration and verification periods

由此可見，各模型均有“擅長”預(yù)報(bào)的月份，為充分發(fā)揮模型各自價(jià)值，實(shí)現(xiàn)多模型間的相互協(xié)調(diào)和性能互補(bǔ)，綜合各模型的預(yù)報(bào)效果，表7 給出了汛期各月份推薦使用模型的最優(yōu)預(yù)報(bào)方案。在預(yù)報(bào)各月來水時(shí)，采用該月份對(duì)應(yīng)的推薦模型，除6月和10 月份，其他月份合格率均能在70%以上，并且在汛期來水最重要的7、8月份預(yù)報(bào)合格率能達(dá)90%以上。

表7 汛期各月份推薦使用模型Tab.7 Recommended model for each month in the flood season

從各模型在不同月份預(yù)報(bào)的整體效果來看，6 月預(yù)報(bào)效果較差，主要是受流域初始狀態(tài)的影響，6月初流域內(nèi)中小水庫基本呈庫空狀態(tài)，同時(shí)該時(shí)段為流域灌溉期，因而該月徑流受流域蓄水影響較大。7、8月份預(yù)報(bào)精度較高，是因?yàn)樵摿饔蛟?、8月份降水量較大，徑流量也大，從而受其他影響因素（包括流域初始狀態(tài)、水利工程等因素）較小，允許誤差大（見表1），因此預(yù)報(bào)合格率較高，預(yù)報(bào)效果較好。而9 月和10 月份預(yù)報(bào)精度較8月份降低，是因該流域于9月份開始降水逐漸減少，而上游眾多水庫會(huì)在9 月和10 月份蓄水截留導(dǎo)致預(yù)報(bào)水平降低。為進(jìn)一步分析各月份預(yù)報(bào)效果的差異，計(jì)算汛期各月徑流的變差系數(shù)（Cv）分析徑流的不確定性，結(jié)果見圖2，Cv值越大表示該月來水不確定性越大。由圖2 知，5 月份Cv值最低，考慮是5 月徑流大多由流域退水產(chǎn)生，來水相對(duì)穩(wěn)定；而在6 月進(jìn)入汛期后，徑流受天然降水影響較大，不確定性較大，因此5月份預(yù)報(bào)效果優(yōu)于6 月份；而10 月份Cv值最大，反應(yīng)該月受降水和水利工程蓄水等影響導(dǎo)致來水不確定性最大，因而預(yù)報(bào)效果較差。

圖2 汛期各月實(shí)測(cè)徑流不確定性分析Fig.2 Uncertainty analysis of monthly measured runoff in flood season

圖3為年徑流預(yù)報(bào)各模型精度，對(duì)于年徑流預(yù)報(bào)，各模型在率定期的預(yù)報(bào)合格率除了TR 模型外均達(dá)到90%以上，取得了較高的預(yù)報(bào)精度。驗(yàn)證期的預(yù)報(bào)精度有所下降，除了SR 模型外，其他模型預(yù)報(bào)合格率都在80%以上，但BP-ANN 和SVM 等機(jī)器學(xué)習(xí)模型預(yù)報(bào)精度明顯高于傳統(tǒng)的統(tǒng)計(jì)模型，其中SVM 模型的合格率為88%，BP-ANN 和PCA-BP-ANN 模型的合格率在90%以上，比統(tǒng)計(jì)模型高10%左右。此外，由于年徑流預(yù)報(bào)的預(yù)見期較長，不確定性大，為了防止個(gè)別模型的預(yù)報(bào)誤差較大而誤導(dǎo)決策，推薦選用BP-ANN 和PCA-BP-ANN 模型，同時(shí)綜合考慮其他模型的預(yù)報(bào)結(jié)果，最終確定合適的年徑流預(yù)報(bào)結(jié)果。

圖3 年徑流預(yù)報(bào)各模型在率定期和驗(yàn)證期的模擬精度Fig.3 The simulation accuracy of annual runoff forecast models for both calibration and verification periods

3 結(jié)論

耦合相關(guān)性分析方法和向前搜索包裹法確定桓仁流域汛期各月徑流與年徑流預(yù)報(bào)模型輸入，并基于主成分分析重組因子以簡化模型輸入，綜合考慮統(tǒng)計(jì)方法和機(jī)器學(xué)習(xí)方法，分別建立年徑流預(yù)報(bào)模型和汛期月徑流預(yù)報(bào)模型，對(duì)比分析各模型的預(yù)報(bào)水平，給出最佳預(yù)報(bào)方案，以提高桓仁流域長期徑流預(yù)報(bào)的精度。主要結(jié)論如下：

（1）在月徑流預(yù)報(bào)中，篩選出6 個(gè)關(guān)鍵影響因子，且上個(gè)月月降雨量起主要作用，其余6 個(gè)大氣環(huán)流因子與預(yù)報(bào)對(duì)象的相關(guān)性和上月月降雨量水平相當(dāng)，相關(guān)性系數(shù)均在0.5 以上；年徑流預(yù)報(bào)中，篩選出7個(gè)關(guān)鍵影響因素，均為大氣環(huán)流因子。

（2）基于多元線性回歸（MLR）、多元門限回歸（TR）、逐步回歸（SR）、人工神經(jīng)網(wǎng)絡(luò)（BP-ANN）、支持向量機(jī)（SVM）和結(jié)合主成分分析的人工神經(jīng)網(wǎng)絡(luò)（PCA-BP-ANN）六種方法，構(gòu)建了桓仁流域年徑流預(yù)報(bào)模型和汛期月徑流預(yù)報(bào)模型。對(duì)于年徑流預(yù)報(bào)，BP-ANN、SVM 和PCA-BP-ANN 三個(gè)機(jī)器學(xué)習(xí)模型的預(yù)報(bào)效果優(yōu)于傳統(tǒng)統(tǒng)計(jì)模型，其合格率要比統(tǒng)計(jì)模型高10%左右；對(duì)于汛期月徑流預(yù)報(bào)，各模型均有自己“擅長”預(yù)報(bào)的月份，利用最佳預(yù)報(bào)方案預(yù)報(bào)時(shí)，在汛期來水最重要的7、8 月份預(yù)報(bào)合格率均能達(dá)90%以上。

在進(jìn)行流域長期徑流預(yù)報(bào)時(shí)，考慮徑流形成的物理基礎(chǔ)和影響徑流變化的各因素選擇合適的預(yù)報(bào)因子，同時(shí)綜合多種預(yù)報(bào)方法擇優(yōu)選擇最佳預(yù)報(bào)方案，可有效降低流域長期徑流預(yù)報(bào)的不確定性。此外，研究采用預(yù)見期為一個(gè)月，可滿足多模型綜合預(yù)報(bào)方法在該研究區(qū)有效性的論證，而該方法在不同預(yù)見期的有效性可能會(huì)有所區(qū)別，還需結(jié)合實(shí)際調(diào)度需求進(jìn)一步檢驗(yàn)，探索不同預(yù)見期下的最優(yōu)預(yù)報(bào)方案，以更加合理地指導(dǎo)水庫調(diào)度。