薛夢楠，王麗學(xué)，張 靜，王明亮，王永東

(1. 沈陽農(nóng)業(yè)大學(xué)水利學(xué)院，沈陽 110866；2. 遼寧省水文水資源勘測局，沈陽 110003； 3. 遼寧省撫順水文局，遼寧 撫順 113015)

0 引 言

流量是反映水資源和江河湖庫水量變化的基本信息。進行流量測驗時應(yīng)在保證資料完整與精度的條件下，針對具體的測驗條件選用不同的測流方法。對于水文站而言，測取到完整的洪水過程及相應(yīng)的水文信息是非常重要的。河道流量演算是分析研究洪水波在河道中的運動規(guī)律，以圣維南(Saint-Venant)方程組為理論基礎(chǔ)，基于河段流量的蓄泄情況和水量平衡原理，通過定量計算河段上游的流量過程進而得到河段下游的流量過程的方法[1]。

本次所研究的姚家山水電站下游10 km處建有占貝水文站，其距姚家山水電站較近且流量主要受姚家山水電站壩下流量影響，因此，本文針對姚家山電站至占貝水文站之間河段，選取姚家山壩下泄流資料和占貝測驗斷面的測流數(shù)據(jù)，在改進的馬斯京根法的基礎(chǔ)上對流量演進系數(shù)進行推求，求得最優(yōu)值，進而通過建立流量演進模型對研究河段的流量演進進行模擬，利用上游姚家山水電站泄流計算流量推算占貝水文站測流斷面流量數(shù)據(jù)，為水文站提供備用資料，從而達到減少水文站觀測，降低水文站運行成本及資料整編難度的目的。目前還未發(fā)現(xiàn)有相關(guān)人員研究相近電站與水文站河段的流量演進規(guī)律，進行流量預(yù)測的研究工作。這一過程的實現(xiàn)，可以對水文站流量整編資料進行校核，并為類似水文站提供借鑒，因此本研究具有重要的現(xiàn)實意義。

1 研究方法

河道流量演算是基于圣維南方程組(Saint - Venant)的基本理論，將河段上斷面流量過程進行演算從而得到下斷面的流量過程。而流量從上斷面至下斷面?zhèn)鞑サ臅r間確定則是進行流量預(yù)報、洪水演進等的基礎(chǔ)。在建立預(yù)報模型進行流量模擬的過程中，首要任務(wù)是流量傳播時間的確定，且傳播時間的準確性直接關(guān)系到預(yù)報模型的模擬結(jié)果及預(yù)報的精度。在研究河流洪水預(yù)報的過程中，流量的傳播時間大多是在合成流量法、相應(yīng)流量(水位)法等基礎(chǔ)上確定的[2]。本文采用灰色系統(tǒng)理論中的灰色關(guān)聯(lián)分析法[3]，計算河段上、下斷面同一時段、向前一個或幾個時段的流量時間序列間的灰色關(guān)聯(lián)度值，對比計算的關(guān)聯(lián)度值的大小，估算河段的流量傳播時間。同時運用馬斯京根改進的建模方法推求最優(yōu)演算參數(shù)，建立河段的流量演算模型。

1.1 灰色關(guān)聯(lián)分析法

灰色關(guān)聯(lián)度可以用來描述兩個或多個序列之間關(guān)系的強弱、大小以及相應(yīng)的次序，是衡量各個事物、因素之間關(guān)系的一個指標，同時通過比較它們序列曲線之間的相似程度來判斷關(guān)聯(lián)程度的好壞。如果要比較的兩個序列曲線的變化幅度、變化方向及速度等基本一致，則說明兩者的關(guān)聯(lián)程度好；否則，兩者關(guān)聯(lián)程度差。

絕對值關(guān)聯(lián)度和速率關(guān)聯(lián)度是灰色關(guān)聯(lián)度的兩種主要的計算方法。在絕對值關(guān)聯(lián)度的計算過程中會涉及數(shù)據(jù)資料序列中的極值，因而計算結(jié)果容易受到數(shù)據(jù)資料序列中最大值和最小值的影響，如果數(shù)據(jù)的序列中存在極值，則極值處的關(guān)聯(lián)度就會發(fā)生變化，數(shù)據(jù)序列的關(guān)聯(lián)程度就會產(chǎn)生偏差，因此有時計算的絕對值關(guān)聯(lián)度不能準確地反映數(shù)據(jù)序列之間關(guān)聯(lián)程度的大小[4]。而速率關(guān)聯(lián)度主要是用來衡量兩個事物、因素在發(fā)展變化的過程中相對變化速率的關(guān)聯(lián)程度的大小。如果兩個事物、因素在發(fā)展變化過程中的相對變化速率基本一致，則說明兩者的關(guān)聯(lián)程度很好；如果兩個事物、因素在發(fā)展過程中的相對變化速率差距較大，則說明兩者的關(guān)聯(lián)程度很差。因此，速率關(guān)聯(lián)度在計算的過程中考慮了兩者變化的一致性，彌補了絕對值關(guān)聯(lián)度的易受極值影響的缺點。鑒于以上原因，本文采用灰色速率關(guān)聯(lián)度[5]對研究河段的流量傳播時間進行計算。

灰色速率關(guān)聯(lián)度的計算，假定參考的時間序列為：Y0= [Y0(1)，Y0(2)，…，Y0(n)]；比較的時間序列為：Xi=[Xi(1)，Xi(2)，… ，Xi(n)]，i= 1，2，…，N。則關(guān)聯(lián)函數(shù)的表達式為：

(1)

其中，ΔXi(t)=Xi(t+1)-Xi(t)；ΔY(t)=Y0(t+1)-Y0(t)。

式中：ΔXi(t)為第i個斷面t時的比較流量系列的相對變化；ΔY(t)為同一斷面t時參考流量系列的相對變化；Xi(t)為第i個斷面t時比較的流量系列；Y(t)為同一斷面t時參考的流量系列；Δt為計算的時段長度。

灰色速率關(guān)聯(lián)度的表達式為：

(2)

1.2 基于馬斯京根改進的建模方法

改進的建模方法判斷最優(yōu)流量演進參數(shù)的標準為：下斷面流量的演算值和實測值的相對偏差平方和最小，即目標函數(shù)為：

(3)

式中：ΔQi為i時段流量演算值跟實測值之間的誤差；Qi為i時的實測流量；Q為i時的模擬流量；C0、C1、C2皆為馬斯京根法流量演算系數(shù)；Ii為i時河段的上斷面的實際入流量。

依據(jù)最小二乘法、馬斯京根的基本原理C0+C1+C2=1，可以得到C0、C1的最優(yōu)估計值[6]分別為：

(5)

1.3 模型的檢驗方法

(1)采用相關(guān)系數(shù)R、徑流總量相對誤差E、合格率QR及一致性指標d對模型模擬的精度進行評定。

相關(guān)系數(shù)R。R越接近1，表明模擬的流量越符合實際情況[7]。R按照下式計算：

(6)

徑流總量的相對誤差E。E< 20%為合格[8]。E按照下式計算：

(7)

合格率QR。當模擬的誤差小于允許的誤差時，為合格[9]。合格模擬的次數(shù)與總次數(shù)之比的百分數(shù)為合格率。QR≥85.0%為甲等，85.0% >QR≥ 70.0%為乙等，70.0% >QR≥ 60.0%為丙等。QR按照下式計算：

(8)

一致性指標d。d越接近1，表明模擬的流量越符合實際情況[10]。d按照下式計算：

(9)

(2)采用流量相關(guān)關(guān)系[11]的數(shù)學(xué)表達式，將其模擬的結(jié)果與模型模擬的結(jié)果進行比較，對模型模擬的優(yōu)良性進行檢驗。

2 實例應(yīng)用

2.1 研究區(qū)基本概況

姚家山水利樞紐工程位于蘇子河干流上，遼寧省撫順市的新賓縣上夾河鎮(zhèn)南溝村姚家山河北屯下游的1.2 km處，是蘇子河流域水能開發(fā)規(guī)劃中穆家水庫(大伙房水庫輸水工程出口)以下河段4級水電開發(fā)的第二級。姚家山水利樞紐工程是一座以發(fā)電為主，兼有防洪作用的中型水利水電樞紐工程。自左岸至右岸的主要建筑物依次為擋水壩段、電站壩段及電站廠房、底孔壩段、溢流壩段。姚家山水利樞紐工程的校核洪水位為189.17 m，總庫容為5 880 萬m3。正常蓄水位183.0 m，相應(yīng)庫容為3 008 萬m3。溢流壩為開敞式，堰頂高程為183.00 m，當上游水位高于183.00 m，溢流壩開始泄洪。底孔每年排沙一次。

占貝水文站始建于2006年1月，位于遼寧省撫順市撫順縣湯圖鄉(xiāng)，東經(jīng)124°24′，北緯41°50′。本站處于渾河一級支流蘇子河上，集水面積為1 902 km2，距離河口28.0 km。按水文站網(wǎng)規(guī)劃原則，確定為國家基本站、省級輔助報汛站。測站配有纜道自動測流控制系統(tǒng)、水情自動化遙測報訊系統(tǒng)、智能語音報訊電話等先進測報設(shè)施設(shè)備。本站的凍結(jié)基面為假定基面。流量測驗精度為二類，泥沙測驗精度為二類。蘇子河屬山溪性河流，屬東部山區(qū)，植被較好。上游有紅升中型水庫一座，控制面積為78.5 km2，穆家電站一座，控制面積為1 711 km2，用于發(fā)電為當年調(diào)節(jié)，對該站削峰約10%左右。沿岸耕地為黏土和沙壤土，河道坡度較大。本站測驗河段河道順直約800 m，無支流匯入，左岸為峭壁，右岸為沙壤土。斷面呈矩形，均為單式河槽，河床質(zhì)為粗沙及卵石組成，河床多年基本穩(wěn)定。高水時河床有沖淤變化，年際間水位流量關(guān)系多有變化，資料整編時多用臨時曲線法。

姚家山水電站上游建有穆家水電站及大伙房引水工程，建設(shè)及運行中人為對河道流量的調(diào)整，使得水文站的測驗工作面臨著流量變化迅速而難以精準把握的難題，同時水利工程建設(shè)改變了水文測驗河段天然河流的變化規(guī)律及測驗條件，對水文站的測流工作影響也越來越大。

2.2 模型建立

以姚家山占貝河段為本次研究對象。此河段長度約為10 km，區(qū)間內(nèi)有一支流南嘉河但觀測數(shù)據(jù)顯示其流量極小，在研究中可忽略不計。綜合考慮姚家山電站的泄流時間和發(fā)電泄流資料的連續(xù)性、完整性后選取2015年8月、9月、10月份的姚家山電站壩下逐日平均下泄流量資料和占貝水文站對應(yīng)的流量資料作為模型的演算時段，進行流量演進參數(shù)的率定。

2.2.1 傳播時間的推求

鑒于洪水波的運動原理，在運動過程中同一位相的流量在上斷面和下斷面出現(xiàn)的時間會有延遲，這一時間即為流量的傳播時間。根據(jù)灰色關(guān)聯(lián)度的概念可以得出，上、下斷面對應(yīng)的流量過程的相關(guān)程度是比較大的。因此，在進行計算時選取河段下斷面的流量過程作為參考的時間序列，上斷面相應(yīng)時段、向前一個、兩個、三個時段…的流量過程作為進行比較的時間序列[12]。同時，計算參考時間序列和各個比較時間序列之間相應(yīng)的灰色關(guān)聯(lián)度值(分別用rt，rt-1，rt-2，…進行表示)，計算得到的關(guān)聯(lián)度的值越大，說明兩個對應(yīng)時間序列之間的關(guān)聯(lián)程度越大，計算得到關(guān)聯(lián)度的最大值對應(yīng)的比較時間序列的向前時段數(shù)即為流量從河段上斷面至下斷面的傳播時間。因此，需要比較這些時間系列相應(yīng)的灰色關(guān)聯(lián)度的計算值，最大值對應(yīng)的比較時間序列的向前時段數(shù)就是流量從上斷面至下斷面的傳播時間。

選取姚家山占貝河段下游的占貝水文站8、9、10月份測驗斷面的逐日平均流量的時間序列作為參考的時間序列，選取上游姚家山電站相應(yīng)時期的對應(yīng)日、向前一日、兩日、三日、…、的壩下逐日平均泄流量的時間序列作為比較的時間序列。利用灰色速率關(guān)聯(lián)度的計算公式，分別計算姚家山占貝河段8、9、10月份的參考時間序列與對應(yīng)的比較時間序列的灰色關(guān)聯(lián)度值，進而對計算結(jié)果進行比較，得出相應(yīng)月份流量從姚家山壩下傳播至占貝測驗斷面的時間。傳播時間的計算結(jié)果見表1。

表1 姚家山占貝河段的灰色速率關(guān)聯(lián)度值Tab.1 Correlation value of gray rate between Yaojiashan and Zhanbei

根據(jù)表1，比較計算的5組灰色速率關(guān)聯(lián)度值，值越大說明關(guān)聯(lián)程度越大，最大值對應(yīng)的時段的關(guān)聯(lián)程度最大，其對應(yīng)的時段為相應(yīng)月份的流量傳播時間。由表1可以看出，姚家山占貝河段的流量傳播時間為對應(yīng)日，即本文確定的姚家山占貝河段的傳播時間為對應(yīng)日。

同時應(yīng)該注意，采用灰色速率關(guān)聯(lián)度法計算得到的是流量從上斷面至下斷面?zhèn)鞑r間的平均值，雖然和特定的某一個時間或時段的傳播時間可能不完全一致，但是確定的傳播時間很大程度上是可以反映實際傳播時間的平均值。

2.2.2 流量演進參數(shù)的推求

根據(jù)灰色關(guān)聯(lián)分析法得到姚家山占貝河段的流量傳播時間為對應(yīng)日，且區(qū)間南嘉河流量經(jīng)觀測數(shù)據(jù)顯示流量值極小即忽略區(qū)間入流，因此，按照無支流的單一河道對該河段進行流量演算，同時采用前文介紹的馬斯京根改進的建模方法對流量演進參數(shù)進行率定，求解最優(yōu)值。推求的流量演算參數(shù)見表2。

表2 河段的流量演算參數(shù)值Tab.2 Calculation of the parameter values

2.3 結(jié)果分析

選擇2015年7-10月的資料進行分析，繪制姚家山電站的壩下逐日平均泄流量過程和占貝水文站的逐日平均流量過程線來分析該河段的流量傳播時間，對應(yīng)傳播時間見圖1。雖然河段沿程存在的一些損耗以及小支流的匯入等不確定因素會對流量的傳播有一定的影響，但由圖1可以看出峰谷的對應(yīng)時間還是比較明顯的，據(jù)此可以得到流量的傳播時間基本都為對應(yīng)日，與由灰色關(guān)聯(lián)分析法計算的流量傳播時間一致。因此，本文運用灰色速率關(guān)聯(lián)度法推求的流量傳播時間是合理、可靠的。

圖1 2015年7-10月姚家山電站和占貝站的日平均流量過程Fig.1 The average daily flow of Yaojiashan Hydropower Station and Zhanbei Station from July to October in 2015

2.4 模型檢驗

對于模型精度的檢驗，選取2015年7-10月的實測平均流量資料進行參數(shù)的率定和模型的驗證，同時運用模型對2016年3月、4月占貝站的流量進行模擬。模型驗證期與模擬期的評定結(jié)果見表3。

對于模型優(yōu)良性的檢驗，同樣選取2015年7-10月的實測平均流量點繪流量關(guān)系圖，確定函數(shù)表達式為：y=0.976 6x+ 1.564 9，見圖2，同時模擬2016年3月、4月的占貝日平均流量，驗證期與模擬期的精度評定見表4。

圖2 姚家山電站流量與占貝水文站流量相關(guān)示意Fig.2 The correlation of flow between Yaojiashan and Zhanbei station

由表3可以得到，本研究建立的模型在驗證期和模擬期的實測流量與模擬流量的相關(guān)系數(shù)分別為0.90、0.83，都在0.83以上；徑流總量的相對誤差分別為-0.04%、6.90%，在容許的范圍內(nèi)；合格率分別為87%、77%，都在70%以上，達到乙級預(yù)報精度；一致性指標分別為0.93、0.89，趨近于1。因此可以得出模型模擬的效果較好。將表3、表4對比可得，不論在驗證期還是模擬期模型模擬的精度較傳統(tǒng)的數(shù)學(xué)表達式模擬的精度相比，前者精度較高且結(jié)果也更為可靠。同時，由模型的精度評定結(jié)果對照水文預(yù)報中流量預(yù)報的精度評定等級可知，模擬的結(jié)果符合進行作業(yè)預(yù)報的精度要求，進而檢驗了用模型進行流量預(yù)報的有效性和可行性。因此，本論文的研究結(jié)果可以應(yīng)用到實際的作業(yè)預(yù)報中。

表3 模型精度評定結(jié)果Tab.3 Model accuracy evaluation results

表4 函數(shù)表達式精度評定結(jié)果Tab.4 Function expression accuracy evaluation results

3 結(jié) 語

(1)利用灰色關(guān)聯(lián)分析法對河段的流量傳播時間進行推算，避免了復(fù)雜相關(guān)曲線的繪制和對影響因素的分析；采用改進的馬斯京根的建模方法，減少了復(fù)雜煩瑣的參數(shù)試算，使得流量演算系數(shù)推求更加簡單易行，同時減小了模擬的誤差，提高了模擬的精度。

(2)目前應(yīng)用水電站泄流量推求水文站測驗斷面流量的方法研究還停留在利用相關(guān)流量建立數(shù)學(xué)函數(shù)關(guān)系的傳統(tǒng)方法中，本文利用流量演進模型推求水文站流量，較傳統(tǒng)的數(shù)學(xué)函數(shù)關(guān)系更為科學(xué)、可靠，且精度較高。

(3)近年來對于流量演算的研究大多是基于各個相近水文站之間河段的演進，本研究將流量演進的河段擴展到水電站與水文站之間，將流量演進模型應(yīng)用到電站泄量的演進過程中，不僅為占貝水文站的測驗斷面流量過程提供參考和備用資料，也為有相似情況的電站和水文站測驗資料的獲取提供思路和依據(jù)。

(4)對于下游水文站測驗斷面流量的推求是基于上游電站的發(fā)電泄量，即為受單一泄量影響的水文站且有連續(xù)發(fā)電泄量下的流量演進模型，有一定的局限性。

□

[1] 包為民. 水文預(yù)報[M]. 北京: 中國水利水電出版社, 2009.

[2] 王燕生. 工程水文學(xué)[M]. 北京: 水利電力出版社, 1992.

[3] 賀北方, 劉正才. 灰色系統(tǒng)理論方法與應(yīng)用[M]. 北京: 氣象出版社, 1995.

[4] 張紹良, 張國良. 灰色關(guān)聯(lián)度計算方法比較及其存在問題分析[J]. 系統(tǒng)工程, 1996,14(3):45-49.

[5] 左其亭, 崔永華. 運用灰色關(guān)聯(lián)分析法推求河流流量傳播時間[J]. 水電能源科學(xué), 2005,23(2):7-10.

[6] 衛(wèi)小蓉, 李 彬, 劉廷璽, 等. 黃河內(nèi)蒙段干流流量演進模型[J]. 水電能源科學(xué), 2015,33(9):22-25.

[7] 王力磊, 袁晶瑄, 徐 煒, 等. 桓仁流域汛期旬徑流預(yù)報方法研究[J]. 水文, 2012,32(6):52-55,27.

[8] 孫美云, 劉 俊, 左 君, 等.分段馬斯京根法在黃河龍潼河段洪水預(yù)報中應(yīng)用[J]. 水電能源科學(xué), 2012,30(8):44-46.

[9] 胡宇豐, 張 峰, 李 匡, 等. 珠江梧州水文站枯季月徑流預(yù)報[J]. 南水北調(diào)與水利科技, 2012,10(1):40-44.

[10] 鄒銘盛. 利用水電站泄流量作為水文站流量測驗備用方案的分析[J]. 廣東水利水電, 2014,(7):79-82.

[11] Willmott C J. Some Comments on the Evaluation of Model Performance[J]. Bulletin of the American Meteorological Society, 1982,63(11):1 309-1 313.