蘭芝梅，朱 詠，陳學(xué)林

（1.甘肅省張掖水文站,甘肅 張掖 734000；2.甘肅省水文站,甘肅 蘭州 730000）

1 前言

水文現(xiàn)代化建設(shè)主要是將先進(jìn)的自動化監(jiān)測設(shè)備應(yīng)用到水文測驗(yàn)中,實(shí)現(xiàn)基礎(chǔ)水文數(shù)據(jù)的自動化監(jiān)測、傳輸、分析等功能。應(yīng)用歷年綜合水位～流量關(guān)系曲線進(jìn)行推流計(jì)算,是滿足流量自動化監(jiān)測的關(guān)鍵環(huán)節(jié)。傳統(tǒng)的綜合關(guān)系曲線率定方法主要是采用整編規(guī)范中的t 檢驗(yàn)對相鄰年份實(shí)測流量數(shù)據(jù)（2 個(gè)獨(dú)立樣本）的同分布顯著性檢驗(yàn)分析,將通過t 檢驗(yàn)的3 年～4 年的流量數(shù)據(jù)作為繪制綜合線的依據(jù)[1-3]。為進(jìn)一步保證綜合線的科學(xué)性、合理性,本研究提出采用spss 軟件中單樣本k-s 非參數(shù)檢驗(yàn)和多樣本k-w 非參數(shù)檢驗(yàn)對實(shí)測流量數(shù)據(jù)進(jìn)行分析研究。非參數(shù)檢驗(yàn)在總體分布未知時(shí)有很大的優(yōu)越性,總是比傳統(tǒng)的參數(shù)檢驗(yàn)安全[4],通過非參數(shù)檢驗(yàn)的歷史實(shí)測流量樣本數(shù)據(jù),表明該樣本數(shù)據(jù)總體無顯著性差異,每年的水位流量關(guān)系變化較穩(wěn)定。對于穩(wěn)定的水位流量關(guān)系,可以對其進(jìn)行數(shù)學(xué)公式的選配,以一個(gè)數(shù)學(xué)方程式來表達(dá)這種水位流量關(guān)系[5]。最小二乘法曲線擬合方法是使實(shí)測值與擬合值差的平方和最小,可較準(zhǔn)確地?cái)M合水位流量關(guān)系曲線[6]。

2 實(shí)測流量成果的非參數(shù)檢驗(yàn)分析

正義峽水文站位于張掖市高臺縣羅城鄉(xiāng)天城村,是黑河中游和下游的分界點(diǎn),屬國家重要水文站、省界測驗(yàn)斷面,測驗(yàn)斷面無明顯沖淤變化；高崖水文站位于臨澤縣板橋鄉(xiāng)土橋村以上1.5 km 處,是黑河中游的省級重要站、中央報(bào)汛站,屬縣界測驗(yàn)斷面,測驗(yàn)河段基本順直,河床及兩岸系礫石細(xì)沙組成,屬流沙河床,沖淤變化較大[7]。

2.1 實(shí)測流量單樣本的總體分布分析

k-s 非參數(shù)檢驗(yàn)是基于連續(xù)型計(jì)量資料樣本分布未知的前提下,用樣本本身的信息來檢驗(yàn)樣本來自同一個(gè)總體假設(shè)的一種統(tǒng)計(jì)檢驗(yàn)方法,其原理主要是運(yùn)用某隨機(jī)變量x 的順序樣本來構(gòu)造樣本分布函數(shù),使得能以一定的概率保證x 的理論分布函數(shù)Fn（x）落在某個(gè)范圍內(nèi)[8]。k-s 非參數(shù)檢驗(yàn)統(tǒng)計(jì)量Dn表示樣本所有變量xi,在經(jīng)驗(yàn)分布與理論分布對應(yīng)值之差的最大絕對值[9],表達(dá)式為：

根據(jù)文獻(xiàn)[6]提出當(dāng)樣本容量n<100 時(shí),顯著性水平=0.05,按照由Dn的精確分布制定的D（n,α）Lilliefors 修正臨界值檢驗(yàn)表查驗(yàn),當(dāng)樣本實(shí)際統(tǒng)計(jì)量Dn

表1 單樣本k-s 非參數(shù)檢驗(yàn)結(jié)果

由表1 可得正義峽和高崖水文站單年實(shí)測流量樣本數(shù)據(jù)的實(shí)際統(tǒng)計(jì)量Dn均大于查表值D（n,α）,即單樣本的總體均不服從正態(tài)分布,故傳統(tǒng)方法采用t 檢驗(yàn)分析率定歷年綜合線理論上不嚴(yán)密。

2.2 實(shí)測流量多樣本同分布分析

多獨(dú)立樣本k-w 檢驗(yàn)用于在樣本來自的總體分布未知情況下由多個(gè)獨(dú)立樣本數(shù)據(jù)檢驗(yàn)多個(gè)總體的分布是否存在顯著差異,其原理是將多組樣本數(shù)據(jù)混合后按照升序排序,求出各變量的秩,檢驗(yàn)各組秩的平均值是否存在顯著差異,如果各組秩的均值不存在顯著差異,則推斷多個(gè)總體的分布無顯著差異；如果各組秩的均值存在顯著差異,則推斷多個(gè)總體的分布存在顯著差異,至少有一個(gè)樣本不同于其他樣本[10]。k-w 秩和檢驗(yàn)統(tǒng)計(jì)量為：

式中：N 為總樣本量；k 為樣本組數(shù)；ni為第i 組的樣本量；Ri為第i 組樣本中的秩總和。一般數(shù)組中會存在相等秩值的情況,需按照校正系數(shù)C 進(jìn)行調(diào)整：

當(dāng)每個(gè)獨(dú)立樣本中實(shí)測數(shù)據(jù)大于等于5 時(shí),KWC接近自由度為k-1 的卡方分布,用卡方分布進(jìn)行統(tǒng)計(jì)量檢驗(yàn)[11]。取顯著性水平=0.05,當(dāng)實(shí)際檢驗(yàn)顯著性值>0.05 時(shí),接受原假設(shè),即多個(gè)樣本的總體分布無顯著性差異?；谝陨戏治龇椒?應(yīng)用spss 軟件中k-w 檢驗(yàn)對正義峽水文站和高崖水文站7 個(gè)獨(dú)立樣本進(jìn)行檢驗(yàn)分析,結(jié)果見表2。

表2 多獨(dú)立樣本k-w 非參數(shù)檢驗(yàn)結(jié)果

由以上檢驗(yàn)結(jié)果可知,正義峽水文站基本順直河道的測驗(yàn)斷面相對穩(wěn)定,歷年實(shí)測流量數(shù)據(jù)的總體分布無顯著性差異,即在水力要素變化不大的前提下,每年的水位流量關(guān)系具有穩(wěn)定性；高崖水文站基本順直河道的測驗(yàn)斷面不穩(wěn)定,歷年實(shí)測流量數(shù)據(jù)的總體分布存在顯著性差異,即每年的水位流量關(guān)系會存在一定的差異,不具備率定歷年綜合關(guān)系曲線的條件。故可以采用正義峽水文站2012 年～2018 年實(shí)測流量數(shù)據(jù),率定歷年綜合水位-流量關(guān)系曲線。

3 歷年綜合水位流量關(guān)系曲線擬合

3.1 擬合數(shù)據(jù)分析

正義峽水文站2012 年～2018 年暢流期實(shí)測流量數(shù)據(jù)共239 測次,經(jīng)分析舍棄在水位～流量關(guān)系上不合理的測次11 個(gè),實(shí)際參與分析的共228 測次。其中實(shí)測最低水位2.17 m,最高水位5.70 m,水位變幅3.53 m；實(shí)測最小流量0.123 m3/s,最大流量570 m3/s。根據(jù)文獻(xiàn)[12]中黑河非汛期生態(tài)需水量最小達(dá)到多年平均流量的10%,在水位低于2.44 m 時(shí)需獨(dú)立率定綜合線。綜合以上分析,采用水位高于2.44 m（共202 測次）的實(shí)測流量數(shù)據(jù)擬合暢流期的常規(guī)綜合線；水位低于2.44 m（共26 測次）的實(shí)測流量數(shù)據(jù)擬合暢流期的生態(tài)流量綜合線。

曲線擬合是指采用某種方法求出某條光滑的曲線最佳擬合數(shù)據(jù),最小二乘法的多項(xiàng)式擬合原理是從整體上考慮函數(shù)q（x）=anxn+an-1xn-1+an-2xn-2……+a0同實(shí)測數(shù)據(jù)點(diǎn)（xi,yi）的誤差平方和（SSE）最小[13],即：

從幾何意義上講就是尋求與給定點(diǎn)（xi,yi）的距離平方和最小的函數(shù)曲線y=q（x）,具體由多元函數(shù)求解極值的方法對式（5）中多項(xiàng)式系數(shù)進(jìn)行求解,最終得到最佳多項(xiàng)式擬合函數(shù)。本文衡量擬合程度的指標(biāo)為均方根（RMSE）、確定系數(shù)R-square,其表達(dá)式為：

式中：SSR 為預(yù)測數(shù)據(jù)與原始數(shù)據(jù)均值之差的平方和；SST為原始數(shù)據(jù)與原始數(shù)據(jù)均值之差的平方和,確定系數(shù)越接近1,表明方程的自變量對函數(shù)的解釋能力越強(qiáng)。

3.2 暢流期常規(guī)歷年綜合關(guān)系曲線

參與綜合線擬合的數(shù)據(jù)點(diǎn)共202 個(gè),設(shè)自變量為水位z,因變量為流量q,經(jīng)過擬合分析可得最優(yōu)擬合函數(shù)：

綜合線的RMSE=1.2978,R-square=0.9999。為滿足高水推流的要求,將擬合函數(shù)高水延伸到歷年最大實(shí)測流量值850 m3/s。擬合結(jié)果見圖1。

圖1 暢流期常規(guī)歷年綜合水位-流量關(guān)系曲線

3.3 暢流期生態(tài)流量歷年綜合關(guān)系曲線

參與綜合線擬合的數(shù)據(jù)點(diǎn)共26 個(gè),由于在生態(tài)流量期對推流精度要求高,因此根據(jù)數(shù)據(jù)點(diǎn)的散點(diǎn)圖分布情況,取其中15 個(gè)率定綜合線1,剩余11 個(gè)率定綜合線2。設(shè)自變量為水位z,因變量為流量q,經(jīng)過擬合分析可得綜合線1 的最優(yōu)擬合函數(shù)：

綜合線2 的最優(yōu)擬合函數(shù)：

綜合線1 的RMSE=0.1178,R-square=0.9877；綜合線2的RMSE=0.0772,R-square=0.9938。擬合結(jié)果見圖2。

圖2 暢流期生態(tài)流量歷年綜合水位-流量關(guān)系曲線

4 誤差分析

根據(jù)《水文資料整編規(guī)范》中定線精度要求,穩(wěn)定的水位流量關(guān)系曲線應(yīng)計(jì)算實(shí)測點(diǎn)對關(guān)系線的隨機(jī)不確定度和系統(tǒng)誤差,且一類精度水文站的系統(tǒng)誤差絕對值小于1%；隨機(jī)不確定度小于8%。經(jīng)對前文得到的綜合水位-流量關(guān)系曲線誤差計(jì)算,可得表3。

表3 水位流量綜合關(guān)系線型定線精度

正義峽水文站2019 年1 月～2020 年4 月暢流期實(shí)測流量共5 測次,應(yīng)用本文暢流期常規(guī)歷年綜合線進(jìn)行推流,可得數(shù)據(jù)對比表4。

表4 實(shí)測流量與綜合線推流對比

由表3 可得,本文擬合的綜合線定線精度符合規(guī)范要求；由表4 可得,相對誤差均小于規(guī)范要求。

5 結(jié)論與討論

5.1 結(jié)論

本文采用多樣本k-w 非參數(shù)檢驗(yàn)對河床基本穩(wěn)定的正義峽和測驗(yàn)斷面處河床不穩(wěn)定的高崖水文站的參與綜合線擬合的多樣本數(shù)據(jù)進(jìn)行對比分析,得出正義峽水文站歷年實(shí)測流量多樣本數(shù)據(jù)的總體無顯著性差異,而高崖水文站歷年實(shí)測流量多樣本數(shù)據(jù)的總體有顯著性差異,故可以采用正義峽水文站歷年實(shí)測流量率定歷年綜合水位-流量關(guān)系曲線；按照q=f（z）函數(shù)對實(shí)測數(shù)據(jù)進(jìn)行擬合分析,率定暢流期水位高于2.44 m 綜合線一條；水位低于2.44 m 的生態(tài)流量期綜合線2條,擬合指標(biāo)均較好；經(jīng)過相關(guān)檢驗(yàn)誤差分析,均符合規(guī)范要求。

5.2 討論

根據(jù)本文分析結(jié)果,可得出歷年綜合線在低水期并非單一線；參與率定綜合線的樣本數(shù)據(jù)來自一個(gè)總體,在水力要素不發(fā)生大的變化前提下,關(guān)系線穩(wěn)定不變。該方法適用于測驗(yàn)斷面處河床基本穩(wěn)定的天然河道。本研究擬合的歷年綜合水位-流量關(guān)系曲線為單一函數(shù)曲線,可作為流量自動化監(jiān)測的技術(shù)支撐。

不足之處：本研究只是針對暢流期進(jìn)行歷年綜合關(guān)系線定線分析,沒有對非暢流期的關(guān)系線分析；測驗(yàn)斷面基本穩(wěn)定,而在水位較低時(shí)會出現(xiàn)不同的歷年綜合關(guān)系線,本文并沒有對其產(chǎn)生的深層原因做分析說明。