【正文】 表。 第三章 預(yù)報(bào)因子挑選 本章從物理成因方面考察各影響因素，選出前述 22 個(gè)因子為初選因子，為了使預(yù)報(bào)因子與預(yù)報(bào)對(duì)象之間具有較好的相關(guān)性，且要求各因子之間能互相獨(dú)立，運(yùn)用統(tǒng)計(jì)學(xué)的方法進(jìn)行第二次篩選，最后選出既有一定的物理意義，又滿足統(tǒng)計(jì)學(xué)原則的因子，作為預(yù)報(bào)因子。 二、單相關(guān)系數(shù) t 檢驗(yàn) t 檢驗(yàn)的計(jì)算公式為： 212 rrnt ??? 檢驗(yàn)時(shí)，可先選一定的信度的 a，再由自由度 f=n2（ n 為樣本數(shù)）查 t分布表，得 ta若按該公式計(jì)算的 tta,可以認(rèn)為在這一信度 下兩者是線性相關(guān)的；若 t≦ ta,則認(rèn)為是不相關(guān)的。 挑選出的預(yù)報(bào)因子及相關(guān)系數(shù)詳見下表 ： 表 31 預(yù)報(bào)因子及相關(guān)系數(shù) 表現(xiàn) 預(yù)報(bào)因子 汛期 月平均 水位 12 月月平均 水位 4 月月平均 水位 時(shí)間 相關(guān)系數(shù) 時(shí)間 相關(guān)系數(shù) 時(shí)間 相關(guān)系數(shù) H500mb(20N、 70、 80、 90E)三點(diǎn)高度和 上一 年 8 月 H500mb(5055N、 7090E+4045N+6585E)巴爾喀什湖區(qū) 1 H500mb(60N、 4050E)兩點(diǎn)高度和 H500mb(120E、 2040N)高度差（沿 120E 線 20N40N） H500mb(25N、 6 7 85E)三點(diǎn)合計(jì) 上一年 12 月 H500mb(4050N120140E)東亞槽區(qū) 500mb 8 點(diǎn)合計(jì) H500mb(2535N、 110130E)長(zhǎng)江中下游區(qū) 7 點(diǎn)合計(jì) H500mb(2030N、 80100E)印緬區(qū) 8 點(diǎn)合計(jì) 上一年 11 月 上一年 11 月 H500mb(120E、 4050N)高度差（沿 129E 線 40N50N） 上一年 1 月 上一年 12 月 上一年 1 月 H500mb(5060N、 100120E)貝加爾湖區(qū) 8 點(diǎn)合計(jì) 上一年 8 月 H500mb(3040N、 8090E)西安高原子 6 點(diǎn)合計(jì) 上一年 3 月 H500mb(1525N、 110130E)南海區(qū) 7 點(diǎn)合計(jì) 上一年 7 月 Q588(105180E)付高強(qiáng)度指數(shù) 上一年 7 月 M588(105180E)付高面和強(qiáng)度指數(shù) 西風(fēng)風(fēng)速 (m/s)(105E、 3540N) 上一年 12 月 西風(fēng)風(fēng)速 (m/s)(105E、 ) 太平洋高壓 Q584 的緯度（ 100120E）的平均位置 烏拉爾 地區(qū)平均高度 H 烏 (6970E、 5060N) 鄂海平均高度 H(135150E、 4560N) 上一年 10 月 上一年 8 月 上一年 1 月 上一年 2 月 西風(fēng)風(fēng)速 (m/s)(105E、 2535N) C102 102 站西風(fēng)指數(shù) (115E、 2530N) 上一年 11 月 上一年 12 月 C836 836 站西風(fēng)指數(shù) (130E、 3035N) 上一年 11 月 上一年 8 月 第四章 多元回歸模型 回歸分析是研究因變量和自變量之間變動(dòng)比例關(guān)系的一種方法，最終結(jié)果一般是建立某種經(jīng)驗(yàn)性的回歸方程。 總殘差平方和為 ： ????????????ntmtmtttntttxbxbxbbyyyQ122211012). .. . ..()?( （ 42） 其依賴于 bi (i=0， 1,2?? m)，要使其最小，則 0???ibQ。即對(duì)式（ 43）進(jìn)行變換。2239。121139。平均 水位 方程成果如下： [汛期 回歸方程 ] y=+03X6+ 式中： y—臨淮關(guān) 水文站 汛期月 平均 水位 X1—上一年 8 月 H500mb(20N、 70、 80、 90E)三點(diǎn)高度和 X2—上一年 12 月 H500mb(25N、 6 7 85E)三點(diǎn)合計(jì) X3—上一年 11 月 H500mb(2030N、 80100E)印緬區(qū) 8 點(diǎn)合計(jì) X4—上一年 1 月 H500mb(120E、 4050N)高度差（沿 129E 線 40N50N） X5—上一年 10 月 鄂海平均高度 H(135150E、 4560N) X6—上一年 11 月 C102 102 站西風(fēng)指數(shù) (115E、 2530N) X7—上一年 11 月 C836 836 站西風(fēng)指數(shù) (130E、 3035N) [12 月 回歸方程 ] y=+++ 式中： y—臨淮關(guān) 水文站 12 月月 平均 水位 X1—上一年 11 月 H500mb(2030N、 80100E)印緬區(qū) 8 點(diǎn)合計(jì) X2—上一年 12月 H500mb(120E、 4050N)高度差（沿 129E 線 40N50N） X3—上一年 3 月 H500mb(3040N、 8090E)西安高原子 6 點(diǎn)合計(jì) X4—上一年 7 月 H500mb(1525N、 110130E)南海區(qū) 7 點(diǎn)合計(jì) X5—上一年 7 月 Q588(105180E)付高強(qiáng)度指數(shù) X6—上一年 8 月鄂海平均高度 H(135150E、 4560N) X7—上一年 8 月 C836 836 站西風(fēng)指數(shù) (130E、 3035N) [4 月 回歸方程 ] y=+++X6 式中： y—臨淮關(guān) 水文站 4 月月 平均 水位 X1—上一年 1 月 H500mb(120E、 4050N)高度差（沿 129E 線 40N50N） X2—上一年 8 月 H500mb(5060N、 100120E)貝加爾湖區(qū) 8 點(diǎn)合計(jì) X3—上一年 12 月西風(fēng)風(fēng)速 (m/s)(105E、 3540N) X4—上一年 1 月鄂海平均高度 H(135150E、 4560N) X5—上一年 2 月鄂海平均高度 H(135150E、 4560N) X6—上一年 12 月 C102 102 站西風(fēng)指數(shù) (115E、 2530N) 四、模型檢驗(yàn) 復(fù)相關(guān)系數(shù) （ R） yyyy SQRSUR ??? 12 或 當(dāng) R≥ Ra，且愈接近于 1 時(shí)，則回歸效果愈好。 檢驗(yàn)結(jié)果及分析 各月回歸方程復(fù)相關(guān)系數(shù) R、剩余標(biāo)準(zhǔn)差 sy、 F 檢驗(yàn)如下： 表 41 臨淮關(guān) 水文站月平均流量多元回歸方程檢驗(yàn)表 通過檢驗(yàn)，發(fā)現(xiàn) 復(fù)相關(guān)系數(shù) R 檢驗(yàn)中，各方程效果都顯著；方程剩余標(biāo)準(zhǔn)差 sy，相對(duì)來說， 汛期 較大，但相應(yīng)月份的 R 檢驗(yàn)卻都是顯著的；在回歸方程的 F 檢驗(yàn)中， 各 方程 效果顯著。由此可見，方程擬合的過程，其實(shí)就是在系列均值左右一個(gè)范圍內(nèi)擬合的過程，所以對(duì)平穩(wěn)系列擬合較好，而對(duì)特大值、特小值等極值頻繁出現(xiàn)的系列，擬合效果是比較差的，只能擬合出變化的趨勢(shì)，往往不能正確的擬合出極值來，具體表現(xiàn)就是相 對(duì)于均值的變化沒有實(shí)測(cè)值大，即振幅較小。之所以形成這種現(xiàn)象，我認(rèn)為主要兩個(gè)方面的原因 ： 一是因?yàn)?水位每天都有觀測(cè)值，相對(duì)于流量實(shí)測(cè)資料來說，水位實(shí)測(cè)資料數(shù)據(jù)多、資料長(zhǎng)，算出的平均水位誤差較小；因此建立起來的回歸方程擬合誤差小，擬合效果較好； 另一方面， 雖然方程的擬合效果較好，但這有可能是 由于 1995— 1999 這 5 年的 實(shí)測(cè)資料正好落在系列均值的范圍內(nèi) ，而由于本站資料年限較短，所以只取了 5 年作預(yù)報(bào)檢驗(yàn)，我們知道，氣象水文因素的變化存在一個(gè)周期，這是因?yàn)樗臍庀笠蛩厥芴柣顒?dòng)影響所 致，其周期在 11 年左右， 而 這個(gè)周期并不是一成不變的，而是一個(gè)長(zhǎng)期大樣本總結(jié)出來的規(guī)律，有可能連續(xù)兩個(gè)豐水年或枯水年組，所以小樣本的分布對(duì)預(yù)報(bào)誤差的影響是巨大的，也是無法估計(jì)的，而這么短的樣本，想要客觀的檢驗(yàn)出模型的好壞，幾乎是不可能的 。 在汛期，特別是主汛期（ 5 月 ~10 月），預(yù)報(bào)的趨勢(shì)接近，只是極值未能成功預(yù)報(bào)出來，從而說明回歸模型對(duì)樣本系列的預(yù)報(bào)，只是在均值左右范圍內(nèi)進(jìn)行擬合，而對(duì)極大值、極小值等極值擬合的效果較差，這也是汛期極值出現(xiàn)頻繁而導(dǎo)致預(yù)報(bào)效果差的一個(gè)重要原因。通過這次畢業(yè)設(shè)計(jì)，運(yùn)用自已學(xué)到的知識(shí)結(jié)合計(jì)算機(jī)能力來完成，提高了自己的業(yè)務(wù)知識(shí)和在計(jì)算機(jī)方面的應(yīng)用能力。 參考文獻(xiàn) [1] 毛紹雄編，《中長(zhǎng)期水文預(yù)報(bào)》，河海大學(xué)函授部印， 2020 [2] 黃振平編著，水文統(tǒng)計(jì)學(xué)，河海大學(xué)出版社， 2020 年 3 月 [3] 于維忠，水文學(xué)原理，河海大學(xué)出版社 [4] 陳其幸，上海市水情長(zhǎng)期預(yù)測(cè)方法比較研究，碩士論文， 202