基于單值化分析方法的下巴溝水文站水位流量關(guān)系研究

李曉勝

(甘肅省臨洮水文水資源勘測(cè)局，蘭州 730000)

1 背景資料

1.1 自然地理

洮河是黃河流域較大的支流之一，發(fā)源于青海省海南藏族自治州西傾山東麓，于甘肅省永靖縣匯入黃河干流。洮河全程長(zhǎng)673.1 km，比降0.27%～0.36%，流域面積2.55×104km2。洮河河源與河口落差達(dá)2 630 m以上，蘊(yùn)藏有豐富的水能資源[1]。洮河多年平均徑流量為53×108m3，在黃河眾多支流中，其年均來(lái)水量?jī)H次于渭河，居于第二位。

洮河流域降水量空間分布不均，受海拔、氣象等因素的影響，洮河流域年降水量在300～600 mm范圍內(nèi)波動(dòng)變化。洮河流域年降水量在年內(nèi)分布上的差異也較為明顯，突出表現(xiàn)為5-10月份(夏、秋)氣溫高、降水多，且以7月份為最大[2]，而以1和12月份的降水量為最小。洮河徑流過(guò)程季節(jié)變化線與降水呈相似規(guī)律，但峰值相對(duì)降雨稍顯滯后，以7-10月份為豐水期，11月份開始銳減，12-2月份為枯水期，3-6月份則水量逐漸增多。流域豐水期徑流量占全年徑流總量約為56%。

1.2 水文測(cè)站概況

下巴溝水文站是國(guó)家重要水文測(cè)站，建站于1960年1月，位于E103°01′、N34°42′。該水文站是洮河重要控制站，對(duì)收集洮河流域水雨情信息和流域防災(zāi)減災(zāi)發(fā)揮著不可替代的作用。

1.3 水文測(cè)站斷面情況

下巴溝水文站測(cè)驗(yàn)河段控制較好，上下游河道呈S型狀，斷面較穩(wěn)定。基本水尺斷面上游1 km里處建有安果水電站，最大引水流量96.0 m3/s。下巴溝水文站控制斷面于每年4和9月份分別測(cè)有斷面詳細(xì)高程數(shù)據(jù)。以每年兩次測(cè)得的數(shù)據(jù)均值作為當(dāng)年斷面形狀，得到2010-2017年斷面套圖(圖1)，并據(jù)此可知下巴溝水文站測(cè)驗(yàn)段河道斷面相對(duì)穩(wěn)定。平均河底高程與歷年變化線(圖2)同時(shí)顯示，下巴溝水文站斷面高程受沖淤等影響不大，沖淤交替進(jìn)行，保證了控制斷面良好的一致性。穩(wěn)定的河床斷面有利于該水文站控制斷面水位-流量單值化分析研究的開展。

圖1 下巴溝水文站2010-2017年斷面套圖

圖2 下巴溝水文站2010-2017年斷面平均河底高程與逐年變化線關(guān)系

1.4 單值化分析的必要性與可行性

傳統(tǒng)水文測(cè)驗(yàn)方式必須通過(guò)大量測(cè)次的水文測(cè)驗(yàn)掌握水位流量中的復(fù)式繩套變化，以滿足現(xiàn)行人工整編或電算整編的要求，以至于測(cè)驗(yàn)次數(shù)既難以精簡(jiǎn)次數(shù)，又不能停測(cè)或者改流量站為水位站。傳統(tǒng)水文測(cè)驗(yàn)方式已與現(xiàn)代水文工作的高效要求相違背[3]。

水位-流量單值化分析工作，將大幅減輕水文測(cè)驗(yàn)工作的頻次和強(qiáng)度，顯著提高工作效率。這既可優(yōu)化外業(yè)測(cè)量的實(shí)施，又可以提升測(cè)量成果質(zhì)量，實(shí)現(xiàn)以人為本、科學(xué)發(fā)展的管理模式[4-5]，這將極大推進(jìn)下巴溝水文站的現(xiàn)代化進(jìn)程和信息化程度。

下巴溝水文站建立于1960年1月，已采取駐測(cè)測(cè)驗(yàn)方式不間斷收集58年以來(lái)的逐日水文資料，并在連續(xù)監(jiān)測(cè)的基礎(chǔ)上進(jìn)行數(shù)據(jù)分析，完成了本流域的基本水文特征數(shù)據(jù)收集工作。該站控制斷面經(jīng)多年觀測(cè)和大洪水考驗(yàn)，斷面基本穩(wěn)定；其水位流量關(guān)系亦符合《水文巡測(cè)規(guī)范》[6]等有關(guān)單值化分析資料的相關(guān)要求。因此，洮河下巴溝水文站水位-流量關(guān)系單值化分析是可行的。

2 水位流量關(guān)系單值化分析

2.1 資料選取

此次水位流量關(guān)系單值化分析采用下巴溝水文站2010-2017年共8年的逐日水位流量資料。樣本資料包含了自2010-2017年逐日水位流量資料，數(shù)據(jù)的測(cè)驗(yàn)方法均嚴(yán)格按照《河流流量測(cè)驗(yàn)規(guī)范》(GB 50179-1993)[7]的相關(guān)要求進(jìn)行收集，資料的可靠性、一致性、代表性審查結(jié)果均為良好。選取下巴溝水文站上下游相鄰測(cè)站作為輔助站，兩者距下巴溝水文站分別為13和22 km。

2.2 分析方法選取

現(xiàn)行水位-流量關(guān)系單值化分析方法較為常用的有組合因子法[8]、落差指數(shù)法[9]等?？紤]到目前下巴溝水文站測(cè)驗(yàn)斷面相對(duì)穩(wěn)定，水位-流量關(guān)系主要受洪水漲落和下游回水頂托的影響，利用下巴溝河段現(xiàn)有的水文站網(wǎng)資料，最終采用綜合落差指數(shù)法[10]開展水位-流量關(guān)系單值化研究。

綜合落差指數(shù)法計(jì)算水位-流量關(guān)系的公式為：

(1)

式中：q為校正流量，m3/s；Qm為實(shí)測(cè)流量，m3/s；ΔZm為綜合落差，m；α為落差指數(shù)。

ΔZm為綜合落差，是考慮了上游洪水漲落與下游回水頂托的綜合影響，可將其分解為：

ΔZm=km1·ΔZm1+km2·ΔZm2

(2)

式中：ΔZm1為上游輔助水尺至本斷面的落差，m；ΔZm2為本斷面到下游輔助水尺的落差，m；km1為上游輔助水尺至本斷面的落差系數(shù)；km2為本斷面到下游輔助水尺的落差系數(shù)。

整理得出綜合落差指數(shù)法計(jì)算水位-流量關(guān)系的公式為：

q=Qm/(km1·ΔZm1+km2·ΔZm2)α

(3)

2.3 落差系數(shù)及其確定

落差指數(shù)初選值按照取α=0.5，km1、km2的初選值取測(cè)站上游輔助水尺至下巴溝水文站基本水尺之間的距離占上下游輔助水尺之間總距離的比值，即km1=13/35=0.371，km2=22/35=0.629。采用試錯(cuò)法，通過(guò)MATLAB編程計(jì)算，確定最終km1、km2取值。方法為首先固定α=0.5不變，依據(jù)方差或標(biāo)準(zhǔn)差最小的原則，然后按合適的計(jì)算步長(zhǎng)不斷調(diào)整km1、km2的值，并將其代入式(3)中計(jì)算實(shí)測(cè)流量對(duì)應(yīng)的校正流量，分別繪制水位-校正流量關(guān)系曲線。根據(jù)適線情況反復(fù)調(diào)整參數(shù)，直到各年水位-校正流量關(guān)系曲線有較好的線型分布，且呈單一曲線關(guān)系。各年曲線通過(guò)符號(hào)檢驗(yàn)、適線檢驗(yàn)和偏離數(shù)值檢驗(yàn)，最終取值為km1=0.354，km2=0.646，α=0.63。

由此計(jì)算所得的下巴溝水文站水位-校正流量關(guān)系公式為：

q=Qm/(0.354ΔZm1+0.646ΔZm2)0.63

(4)

3 成果檢驗(yàn)

3.1 歷年水位流量演變

選取2010-2017年下巴溝水文站水位-流量數(shù)據(jù)，并點(diǎn)繪歷年水位-流量關(guān)系圖。由圖3可知，2012-2014年水位-流量曲線呈現(xiàn)逐年左移的趨勢(shì)，原因在于水文站上游附近違法占用河道現(xiàn)象嚴(yán)重，使河道收束。經(jīng)2014年底整治后，恢復(fù)河道原貌。但由于氣候原因，2010-2013年的水位-流量關(guān)系曲線均呈現(xiàn)單一關(guān)系曲線，而2014-2017年均出現(xiàn)2條臨時(shí)低水曲線。

圖3 下巴溝水文站歷年水位-流量關(guān)系圖

3.2 水位-流量關(guān)系曲線檢驗(yàn)

根據(jù)《水文巡測(cè)規(guī)范》(SL 195-1997)的要求，采用 2010-2017年實(shí)測(cè)流量中分別選取2011、2013和2017年作為豐、平、枯3個(gè)代表年份進(jìn)行單值化分析，利用綜合線進(jìn)行流量推定，并推算年徑流總量、汛期總量和一次洪水總量，與原有整編成果進(jìn)行統(tǒng)計(jì)分析，統(tǒng)計(jì)結(jié)果見(jiàn)表1。結(jié)果顯示，按照綜合落差指數(shù)法推算所得的結(jié)果，各時(shí)段徑流量在豐水年結(jié)果較原整編方法稍大，而在枯水年則稍小于原整編成果。但通過(guò)該方法法得出的下巴溝水文站水位-流量單值化分析成果的精度，符合《水文巡測(cè)規(guī)范》(SL 195-1997)的相關(guān)要求。

表1 下巴溝水文站各時(shí)段推算徑流量及其誤差分析表

4 結(jié) 論

利用下巴溝河段現(xiàn)有的站網(wǎng)水位流量資料，采用綜合落差指數(shù)法開展水位流量關(guān)系單值化研究，綜合各方面因素得出下巴溝水文站水位-流量單值關(guān)系公式。并根據(jù)《水文巡測(cè)規(guī)范》(SL 195-1997)對(duì)各年的年總量、汛期總量、一次洪水總量分析顯示，該單值關(guān)系公式的誤差符合相關(guān)規(guī)定，是可靠的。

將該成果應(yīng)用于下巴溝水文站的實(shí)際工作中，加快本站水文現(xiàn)代化進(jìn)程[5]，可將大量的人力資源從傳統(tǒng)的測(cè)驗(yàn)方式中解放出來(lái)，同時(shí)可以將有限的人員調(diào)余濟(jì)缺最大限度優(yōu)化人力資源，亦是將有限的技術(shù)資源集中起來(lái)，做好服務(wù)型水文的工作，為社會(huì)提供更加優(yōu)質(zhì)服務(wù)的重要途徑[11-12]。