段江勇 潘暢 陳大玄

摘要：相應(yīng)流量報汛是長江水利委員會水文局水情分中心（以下簡稱“水情分中心”）一項極其重要的日常工作，要做到精準(zhǔn)發(fā)布相應(yīng)流量信息卻是一項技術(shù)難題。通過對不同類型水文測站水位流量關(guān)系進(jìn)行分析研究，提出了在極其復(fù)雜水情條件下的水文測站相應(yīng)流量推算方法與精度控制措施，可供各地水情分中心同行借鑒與參考。

關(guān)鍵詞：相應(yīng)流量;報汛精度;流量推算;控制措施

中圖法分類號：TV131.4

文獻(xiàn)標(biāo)志碼：A

1研究背景

長江水利委員會水文局2005年率先在全國實行自動報汛之后，水文信息從內(nèi)部掌控走向網(wǎng)絡(luò)信息平臺，成為社會共享資源而進(jìn)入了公眾視野，得到了各級政府、各相關(guān)部門和普通百姓的廣泛關(guān)注。確保水文信息發(fā)布的及時性和準(zhǔn)確性無疑是對水文人的嚴(yán)峻考驗。水位、降水實時數(shù)據(jù)已全部實現(xiàn)自動記錄、存儲、傳輸和報汛，而實時相應(yīng)流量報汛更多受斷面等因素所限。目前，水文儀器設(shè)備暫時還無法達(dá)到自動測報條件的流量測站，更多地依賴于建立水位一流量關(guān)系，以水位為參數(shù)推求出相應(yīng)流量，將相應(yīng)流量作為實時流量連同實時水位信息一同發(fā)布。

2相應(yīng)流量連時序法推求方法

在天然河道中，水位實測流量關(guān)系曲線呈單一線的測站數(shù)量極少，在河流中下游地區(qū)，尤其是洞庭湖區(qū)的水文測站特性影響因素極其復(fù)雜，大多水位一流量關(guān)系曲線呈非單一線。在年度整編推算時段流量時，通過實測流量測點按施測時間順序在水位一流量關(guān)系圖上連接成光滑曲線，以時間、水位為參數(shù)在水位一流量關(guān)系連時序曲線上查讀流量，完成年度水文資料整編，行業(yè)稱這種推算流量方法為連時序法或俗稱繩套線法。

相應(yīng)流量報汛與時段流量整編具有相同之處，都是以時間、水位為參數(shù)推算流量，差異在于：年度整編時推求時段流量的水位一流量關(guān)系連時序曲線，是由已知實測流量測點組成的光滑曲線;而某時刻相應(yīng)流量推求的臨時預(yù)報線只掌握前面的實測流量測點，其后的流量實測點的準(zhǔn)確位置未知，需要水情報汛人員預(yù)測。因此，為報汛需要而推求的相應(yīng)流量僅為預(yù)報值。相應(yīng)流量的成果精度取決于水情報汛人員所作預(yù)報線的準(zhǔn)確度，就目前而言，大多數(shù)水情分中心還是依靠工作人員憑經(jīng)驗作線，相應(yīng)流量精度很大程度上依賴于工作人員的責(zé)任心與實踐經(jīng)驗。然而，僅憑這些無法滿足相應(yīng)流量的報汛精度，應(yīng)從專業(yè)角度和技術(shù)層面來解決相應(yīng)流量推求難題，下面介紹兩種在極其復(fù)雜水情條件下連時序法推算相應(yīng)流量的方法。

2.1水位落差連時序關(guān)系線相關(guān)法

對于水位一流量關(guān)系非單一線的測站，在繪制臨時預(yù)報線時一般根據(jù)曲線趨勢延伸，它僅適應(yīng)于水位漲落變化平穩(wěn)，上游來水單一，下游渲泄順暢，支流影響小的洪峰漲落水過程中的中間線段，對于峰頂、峰谷轉(zhuǎn)折處若按曲線趨勢延伸往往會得到相反的結(jié)果，方法上存在局限性，采用水位落差連時序關(guān)系線相關(guān)法可彌補(bǔ)上述方法的不足。

根據(jù)流量變化與水位落差成正比關(guān)系即Q=f（Zm），通過點繪上下游站連時序法水位落差曲線圖，參照其落差變化趨勢來確定同期水位一流量關(guān)系線（即預(yù)報線）。例如，南咀站2016年7月3日10：00水位達(dá)34.62m并持續(xù)上漲，按趨勢延伸曲線查得7月4日08：00水位34.83m，相應(yīng)流量為10400m2/s。如果參照落差變化，7月4日08：00落差為0.82m，比7月3日08：00落差1.03m偏小0.21m，落差減少而流量增大明顯不合理。通過南咀站2016年6月9日至8月15日的落差圖與其同期的水位一流量關(guān)系曲線圖比較（見圖1），不難發(fā)現(xiàn)，水位一流量關(guān)系曲線與落差曲線走向基本一致，落差減小時其相應(yīng)流量減小，反之亦然。再用南咀站2017年6—7月的落差圖與其同期的水位一流量關(guān)系曲線圖進(jìn)行比較，同樣得到上述結(jié)論（見圖2）。利用測站連時序水位落差曲線與水位一流量關(guān)系曲線的相關(guān)性，根據(jù)落差變化趨勢適時修訂預(yù)報線，可適時糾正相應(yīng)流量的偏差，大幅提高相應(yīng)流量報汛精度，是保證精準(zhǔn)預(yù)報水位一流量關(guān)系曲線的最有效的輔助方法。

2.2樣本曲線模型法

在極其復(fù)雜水情條件下推求時段流量時，一般采用連時序法，盡管連時序法形成的繩套曲線極其復(fù)雜，但只要認(rèn)真分析，不難發(fā)現(xiàn)每個測站的連時序法水位一流量關(guān)系繩套曲線都有其獨(dú)特性。若將多年的水位一流量關(guān)系繩套曲線進(jìn)行比較，許多洪峰過程的線型非常類似，有大部分不同年份的曲線可以重合在一起。根據(jù)這一規(guī)律，若把歷史資料中具有代表性的洪峰過程曲線提取出來，配上其影響參數(shù)（如時間、不同河流來水、上下游落差、漲落水趨勢、水位變化區(qū)間和漲落率等），設(shè)計編制專用軟件，將這些繩套曲線作為模型錄入計算機(jī)數(shù)據(jù)庫，水情分中心人員可根據(jù)當(dāng)前測站的影響參數(shù)輸入專用計算機(jī)軟件中進(jìn)行搜索，把類似的曲線模型一一顯示出來，通過比對、篩選、優(yōu)化，選取符合本站當(dāng)前水情的模型，稍作修改或不予修改作為臨時預(yù)報線，再不斷通過實測流量的校正、總結(jié)經(jīng)驗和修改參數(shù)等，相應(yīng)流量報汛不再是難事。如比較南咀站2016年與2017年的高水期6～7月水位一流量關(guān)系曲線（見圖1～2）可知，2017年6月28日至7月5日與2016年7月2-12日兩條曲線非常相似，均是受到下游湘水與資水的嚴(yán)重頂托，水位一流量關(guān)系曲線呈現(xiàn)非常態(tài)的反線套曲線。若2017年7月3—8日預(yù)報線通過上下游落差及區(qū)域水情分析，不難發(fā)現(xiàn)與2016年同期水位落差圖也非常相似，見圖1～2。參照2016年同期水位一流量關(guān)系線來預(yù)測2017年的關(guān)系線，如果把南咀近lOa水位一流量關(guān)系曲線檢索出來，可以找出許多非常相似或完全相似的曲線，再加以分析歸類，提煉出一些比較典型的曲線模型，為水情分中心制作預(yù)報線提供參考，可大大提高相應(yīng)流量報汛精度。

3相應(yīng)流量單值法數(shù)學(xué)模型研究方法

前面介紹的方法能在很大程度上對推求相應(yīng)流量起到定性作用，要精準(zhǔn)推求相應(yīng)流量最為有效的方法還是水位一流量關(guān)系單值化數(shù)學(xué)模型，其基本公式為式中，Qm為實測流量;Qx為校正流量因素;Zm為落差;B為落差指數(shù)。通過公式（1）就可以將Z～Q非單值相關(guān)關(guān)系轉(zhuǎn)化成單值相關(guān)關(guān)系，但僅限于單一河道的水文測站。對于非單一河道，同時用n個落差組成綜合落差，可歸納為一種通用標(biāo)準(zhǔn)數(shù)學(xué)模型為

計算校正流量：

推算時段流量：式中，Zm為綜合落差;Zml、Zm2、…、Zmn為各站輔助落差;Km1、Km1、…、Kmm為各站輔助落差權(quán)重系數(shù);B為綜合落差改正值;B'為落差指數(shù);Qm為實測流量;Qx為校正流量;K2為綜合落差改正值;ZO為基本站水位;Z1、Z2、…、Zn為各輔助站水位;bi、b2、…、bn為落差站凍黃差;對于順逆流站引入K1參數(shù)的概念，K1為順逆流改正系數(shù);Q為時段流量（或稱相應(yīng)流量、或稱實時流量）;Qc為水位校正流量關(guān)系線上流量值。在上述數(shù)學(xué)模型中，共有3+n個待定參數(shù)Ki、K2、B、km1、Km2、…、kmm，通過歷史資料，用水位一流量關(guān)系數(shù)學(xué)模型處理系統(tǒng)分析計算可求得這些參數(shù)。

3.1參數(shù)的確定

在數(shù)學(xué)模型中，有3+n個參數(shù)Ki、K2、Km1、Km2、Km3、…、Kmm、B。其中，Km1、Km2、Km3、…、Kmm為落差水尺的權(quán)重系數(shù)，當(dāng)單一落差時，Km2、Km3、…、Kmn為0，Km1一般取1或某一常數(shù)或某變量Km1=f（Zo或Zm）。當(dāng)輔助水尺采用2～3組或n組時，Km1、Km2、Km3、…、Kmm一般采用距離加權(quán)求得，絕大多數(shù)情況采用試算，Km1至Kmm可以是常數(shù)或變量函數(shù)Kmm=f（Zn或Zm）。數(shù)學(xué)模型中的權(quán)重系數(shù)和落差指數(shù)均采用變動系數(shù)，是基本水位的函數(shù)。也就是說，n組落差值的權(quán)重系數(shù)和落差指數(shù)值隨基本站水位變化而變化。一旦確定輔助水尺，某測站水位一流量關(guān)系數(shù)學(xué)模型的分析就是對n+3組參數(shù)的設(shè)置、計算與優(yōu)化處理過程。

3.1.1 Km1至Kmm/B的試算方法

數(shù)學(xué)模型分析的具體步驟是：①用單站落差作水位流量相關(guān)分析，擇優(yōu)選取代表性強(qiáng)的輔助站作輔助水尺，確定輔助水尺組數(shù)n。②按距離加權(quán)初設(shè)落差權(quán)重系數(shù)，建立水位一流量關(guān)系數(shù)學(xué)模型，擬合水位流量關(guān)系曲線，通過調(diào)整落差權(quán)重系數(shù)，使擬合曲線與水位實測流量關(guān)系曲線的圖型比例與趨勢保持一致。③設(shè)置每個參數(shù)的取值區(qū)域與增減變幅，對Km1至Kmm個參數(shù)進(jìn)行遞增或遞減作3+n層大循環(huán)試算，數(shù)學(xué)模型處理系統(tǒng)自動完成曲線擬合、3種檢驗、誤差統(tǒng)計和方案評估，大循環(huán)試算結(jié)束后從若干個方案中挑選出最佳方案。④再將最佳方案中關(guān)系曲線和Km2至Km4固定反算Km1，建立Km1與基本站水位（或輔助站水位）關(guān)系曲線，得到擬合曲線公式Km1=f（ZO或Zn），再將新的Km1=f（Zo或Zm）替代原有的Km1，反復(fù)進(jìn)行第四項，直到Km1最優(yōu)為止。如此類推，優(yōu)化Km2至Kmm。

B值是基本水位的函數(shù)，可經(jīng)過試算求得，采用變動值對提高數(shù)學(xué)模型成果精度的作用不明顯，但有改正水位校正流量關(guān)系線線型的作用，可將“S”型水位校正流量關(guān)系線改正為二次拋物線，使水位校正流量關(guān)系線趨勢合理。

3.1.2 K1和K2的確定方法

k2為綜合落差改正系數(shù)，當(dāng)測站有順逆時，通過設(shè)置k2來解決正負(fù)流量問題。當(dāng)Zm為負(fù)值時，取k2為負(fù)值，將綜合落差改為正值，便于指數(shù)運(yùn)算;當(dāng)Zm為正值時，取k2為正值，k1為負(fù)流量的改正系數(shù)，實測流量測點大致對稱于正負(fù)流量坐標(biāo)兩側(cè)。對負(fù)流量測點乘以k1時，可將負(fù)流量改為正值與正流量合并定為一條正負(fù)流綜合線，由于正負(fù)流兩線不完全對稱而K1并非常數(shù)，而是隨基本水位變化的變動系數(shù)，通過曲線擬合，當(dāng)負(fù)流時，K1=f（Z0）;當(dāng)正流時K1=1（常數(shù)），順逆站按此方法處理可獲得比較適用的數(shù)學(xué)模型。

下面通過兩個實例介紹說明測站水位一流量關(guān)系數(shù)學(xué)模型的分析過程。

（1）實例一。小河咀站位于西洞庭湖湖口南端（見圖3），水情極其復(fù)雜，上受沅水、澧水、長江來水漲落影響，下受湘、資水的回水頂托影響，同時還受西、南洞庭湖湖泊調(diào)蓄作用，上游選擇二組輔助水尺，即沅水洪道周文廟站Z1，西洞庭湖湖泊代表站沙灣站Z2，下游選南洞庭湖湖泊代表站沅江站Z3，下泄主洪道東南湖站Z4，其單值化方案為

計算校正流量：

（2）實例二。南咀站位于西洞庭湖湖口北端，其測站水情特性極其復(fù)雜，上游受沅、澧水、長江來水的影響，下游受湘、資水的回水頂托。在水位一流量關(guān)系單值化模型中，上游選定了蒿子港、白蚌口站、西洞庭湖泊代表站沙灣站，下游選擇了東南湖與黃茅洲站，共選擇了5組落差輔助水尺，通過單值化數(shù)學(xué)模型處理后，在中高水位其精度可滿足相應(yīng)流量報汛的要求。要達(dá)到年度整編推求時段流量的精度要求，還有待對單值化數(shù)學(xué)模型作進(jìn)一步的分析論證。南咀站2016年與2017年的6～7月均受資水、湘水嚴(yán)重頂托，水位一流量曲線極其反常，給相應(yīng)流量報汛帶來較大難度，造成相應(yīng)流量誤差偏大，但通過單值化處理后，完全滿足相應(yīng)流量報汛要求（見圖5）。

水位一流量關(guān)系單值化數(shù)學(xué)模型分析有專用軟件處理系統(tǒng)，其設(shè)計有多個功能模塊，通過這些模塊可快速尋找水位一流量關(guān)系規(guī)律，準(zhǔn)確判斷3+n，個待定參數(shù)與基本水位、綜合落差、輔助站水位、不同來水時期、漲落率等因素是否存在相關(guān)關(guān)系。通過對基本站水位、輔助站水位、落差、時間的排序，判斷流量相對誤差大小變化，從中尋找3+n個參數(shù)的變化規(guī)律，分析出某一測站的水位一流量關(guān)系單值化數(shù)學(xué)模型。

4相應(yīng)流量推算軟件設(shè)計創(chuàng)新

相應(yīng)流量的精確度是水情分中心最難掌控的一項技術(shù)指標(biāo)，也是防汛部門最為關(guān)注的重點，盡管上級主管部門加大了對水情分中心管理的獎懲措施，對提高成果質(zhì)量只能起到一定督促作用，要真正提高相應(yīng)流量的報汛精度，必須從技術(shù)角度對相應(yīng)流量推求軟件進(jìn)行改進(jìn)與創(chuàng)新。

連時序法相應(yīng)流量推求軟件可設(shè)計為綜合模塊，連時序法曲線顯示框采用米格形式，縱軸為水位，橫軸為流量兼落差，曲線框可無限放大與縮小，某時刻以前已上報的相應(yīng)流量以連時序法（繩套線）形式顯示在圖上，某時刻以后的相應(yīng)流量是未知數(shù)據(jù)，是預(yù)報值。繪制該段曲線可通過點繪上下游站連時序法水位落差關(guān)系線作為其參照值，在模塊中設(shè)計選擇窗口，只要選擇上下游站，就可用不同顏色將連時序法水位落差關(guān)系線與連時序法水位流量點繪在一起，通過對過去一段時間的兩曲線進(jìn)行比較分析就能確定連時序法預(yù)報線的趨勢，最終確定預(yù)報線。經(jīng)過不斷調(diào)整實測流量修訂預(yù)報線，不斷提高相應(yīng)的流量精度。

在模塊中設(shè)計測站連時序法歷史曲線模型選擇窗口，通過對測站名、時間、水位等參數(shù)的選擇，結(jié)合多種方法進(jìn)行綜合分析最終確定預(yù)報線。按此思路設(shè)計的軟件，可改盲目作線為直觀繪圖;可改靠經(jīng)驗主觀作線為以多種參數(shù)為旁證綜合分析定線，這對提高相應(yīng)流量預(yù)報線成果精度是一種非常有效的技術(shù)手段和措施。

5相應(yīng)流量報汛精度控制措施

5.1過程控制措施

5.1.1 建立水情分中心管理制度

建立科學(xué)合理、切合實際的水情分中心管理制度，強(qiáng)化分中心人員工作責(zé)任心，以此來保障報汛產(chǎn)品的成果質(zhì)量，其基本內(nèi)容措施是：①規(guī)定分中心人員在崗工作時間、地點和內(nèi)容。②規(guī)定當(dāng)班人員對設(shè)備定時段、定項目巡視檢查登記。③規(guī)定對硬件故障處理的要求時限。④規(guī)定對報汛軟件掌握的熟練程度。⑤依據(jù)條款規(guī)定對值班人員過錯進(jìn)行責(zé)任追究處理。

5.1.2制定應(yīng)急預(yù)案

制定各項應(yīng)急預(yù)案是保障水情分中心報汛質(zhì)量的重要措施，應(yīng)急預(yù)案分為硬件、軟件及測驗方式3類。硬件包括數(shù)據(jù)接收及處理設(shè)備、服務(wù)器、網(wǎng)絡(luò)傳輸設(shè)備、電源設(shè)備等，平常應(yīng)加強(qiáng)對上述儀器設(shè)備的維護(hù)保養(yǎng)工作，發(fā)生故障時能及時處理，必要時更換備品備件。軟件包括水情數(shù)據(jù)接收系統(tǒng)、水情數(shù)據(jù)處理系統(tǒng)、水情信息交換系統(tǒng)和MicrosoftSQL Server數(shù)據(jù)庫等。這些軟件雖已成熟，但必須預(yù)先在備用機(jī)上裝好，關(guān)鍵時刻能夠及時切換，平時應(yīng)認(rèn)真做好各軟、硬件的測試和應(yīng)急演練工作。測驗方法應(yīng)急預(yù)案是指當(dāng)出現(xiàn)特殊水情或測區(qū)出現(xiàn)潰堤決垸等極端情況時，原有水位一流量關(guān)系被改變，單值化數(shù)學(xué)模型失效，如何改變測驗方式，相應(yīng)流量推求方法應(yīng)有相應(yīng)的應(yīng)急預(yù)案可用。

5.1.3 內(nèi)部培訓(xùn)交流機(jī)制

加強(qiáng)水情分中心報汛人員培訓(xùn)是一項長期的重要工作內(nèi)容，相應(yīng)流量報汛是極其重要的工作，因此應(yīng)定期對報汛相關(guān)人員進(jìn)行培訓(xùn)和開展內(nèi)部交流，水情報汛人員必須熟悉了解本測區(qū)各測站水情特性，熟練掌握水情分中心軟、硬件的操作技能，尤其是新職工必須經(jīng)過嚴(yán)格的培訓(xùn)并考核合格之后方可上崗值班。

5.1.4痕跡材料留存?zhèn)洳闄C(jī)制

水情分中心應(yīng)該制定一套完整齊備的工作手冊，設(shè)計規(guī)范化、標(biāo)準(zhǔn)化的操作工序，工作人員交換班時，應(yīng)規(guī)定下班時應(yīng)交接的事項、接班人員應(yīng)完善工作程序，每道程序的操作過程必須用文字記錄留存痕跡材料。認(rèn)真做好相應(yīng)流量報汛過程控制痕跡材料保存工作，尤其是水位一流量關(guān)系進(jìn)行率定和流量預(yù)報線修定及審查情況。對于實測流量必須在規(guī)定時效內(nèi)接收、處理并發(fā)送，同時做好記錄。保證所有流量報汛痕跡材料規(guī)范完備，便于主管部門例行檢查與事后追責(zé)備查。

5.2責(zé)任追究與考評機(jī)制

5.2.1 大錯成因分析與追責(zé)

水情分中心出現(xiàn)數(shù)據(jù)遲報、漏報或錯報時，都應(yīng)查清問題出現(xiàn)的原因，可能是人為因素、自然事故或儀器突發(fā)故障，查明原因后尋求最佳處置方法，必要時應(yīng)依照上級主管部門制定的有關(guān)條款予以追責(zé)，以杜絕類似事故發(fā)生。

5.2.2相應(yīng)流量計算方式方法缺陷修復(fù)機(jī)制

相應(yīng)流量報汛精度是目前水情分中心最難掌握的項目。不管測站水情多么復(fù)雜，只要認(rèn)真做輔助性分析，均可找到一些自然規(guī)律，認(rèn)真分析后就能夠大幅度提高相應(yīng)流量報汛精度。如某測站相應(yīng)流量誤差長期超限，應(yīng)及時改變方法方式，及時修正存在的缺陷，創(chuàng)新出適合測站特性的方法。

5.2.3 成果考評機(jī)制

每年報汛成果質(zhì)量考評應(yīng)由主管部門組織實施，主要有精度考核、時效考核和痕跡材料考核。精度考核主要針對相應(yīng)流量而言，相應(yīng)流量成果精度應(yīng)是水情分中心重點需要攻克的難題。時效考核主要是實測流量在規(guī)定的時間之內(nèi)是否報出的問題。痕跡材料考核主要是檢驗分中心人員是否按照上級主管部門要求對各項檢查記錄材料存檔。

6結(jié)語

在當(dāng)前的防汛測報工作中，精準(zhǔn)推算極其復(fù)雜水情條件下水文測站的相應(yīng)流量是水情分中心的工作重點和技術(shù)難題。通過不斷摸索與經(jīng)驗積累，提出的水位落差連時序關(guān)系線相關(guān)法、樣本曲線模型法、相應(yīng)流量單值法數(shù)學(xué)模型研究方法及相應(yīng)流量報汛精度控制措施，可為水文同行提供技術(shù)支撐與參考，不同水文測站可選擇適合測站特性的相應(yīng)流量推求方式方法，為各級防汛部門提供準(zhǔn)確的水文數(shù)據(jù)。

參考文獻(xiàn)：

[1]李厚永，張潮，吳瓊.水位流量關(guān)系單值化分析綜合模型研究及應(yīng)用[J]冰文，2011，31（S1）：152-153，157.