胡尊樂，汪 姍 ，潘 杰，周迎穎 ，朱 凱

（1.江蘇省水文水資源勘測局常州分局，江蘇 常州 213022；2.河南省信陽水文水資源勘測局，河南 信陽 464000）

1 概述

水文測驗是水文工作的基礎(chǔ)，水文測驗質(zhì)量尤其是水文測驗成果質(zhì)量更是關(guān)系到水文測站一年工作的成敗。目前，為了保證水文測驗精度，提高水文測驗成果質(zhì)量，向社會提供更加可靠的基礎(chǔ)資料，同時為了加強(qiáng)水文測驗質(zhì)量管理，規(guī)范水文測驗質(zhì)量檢查評定辦法，水利部水文局編制和頒布了《水文測驗質(zhì)量檢查評定辦法（試行）》，其中對水文測驗成果質(zhì)量檢查評定有了明確的賦分標(biāo)準(zhǔn)（見表1）。

從表1中可以看出，流量（泥沙）項目是水文測驗成果質(zhì)量檢查評定的核心，賦分值占總賦分值的46.7%，主要內(nèi)容涉及流量（泥沙）測次布置、測點(diǎn)控制等?！逗恿髁髁繙y驗規(guī)范》（GB50179-2015）中對流量測驗次數(shù)的布置和測點(diǎn)的控制，有這樣的規(guī)定：“水文站一年中的測流次數(shù)，必須根據(jù)高、中、低各級水文的水流特性、測站控制情況和測驗精度要求，掌握各個時期的水情變化，合理地分布于各級水位和水情變化過程的轉(zhuǎn)折點(diǎn)處”。顯然，良好的流量測驗次數(shù)時空分布，要能準(zhǔn)確反映一年之中的水情變化，要能滿足推算逐日流量和各項特征值的要求。但由于一年之中水情變化（尤其是暴雨洪水變化）的不可預(yù)知性，以及一些客觀條件的限制，水文測站在開展流量測驗時，往往測次分布不夠合理、測點(diǎn)控制不盡到位，有時甚至錯失洪峰，對水位流量關(guān)系曲線的定線精度造成了不可彌補(bǔ)的影響，進(jìn)而影響到水文資料的可靠性。因此，《水文資料整編規(guī)范》（SL 247-2012）和《水文年鑒匯編刊印規(guī)范》（SL 460-2009）要求對水文測站穩(wěn)定的水位流量關(guān)系（臨時曲線法的主要曲線及經(jīng)單值化處理的單一線，下同）進(jìn)行關(guān)系點(diǎn)對關(guān)系線的標(biāo)準(zhǔn)差（Se）和隨機(jī)不確定度（X'Q）計算，且標(biāo)準(zhǔn)差和隨機(jī)不確定度滿足不同精度要求[1]。其中，對于測點(diǎn)在10個以上的上述曲線，均應(yīng)進(jìn)行符號檢驗、適線檢驗和偏離檢驗。另外，為確保水文資料質(zhì)量，對于水文測站因故未能測得洪峰流量或最枯水流量，允許對其當(dāng)年水位流量關(guān)系曲線高水或低水做適當(dāng)延長，以滿足推求全年完整流量過程的要求。其中，高水部分延長不應(yīng)超過當(dāng)年實(shí)測流量所占水位變幅的30%，低水部分延長不應(yīng)超過10%。

表1 水文測驗成果質(zhì)量檢查評定賦分表（水文站，300分）

顯然，《河流流量測驗規(guī)范》對流量測驗次數(shù)的布置和測點(diǎn)的控制有一個定性的規(guī)定，《水文資料整編規(guī)范》和《水文年鑒匯編刊印規(guī)范》對水位流量關(guān)系定線成果有一個定量的規(guī)定，都在一定程度上保證了流量測驗成果的可靠性。但由于各地河流的差異性，對于流量測驗的時空分布還缺乏明確的評價方法。目前，我國各地水文站基于《河道流量測驗規(guī)范》中高、中、低水的有關(guān)規(guī)定，采用一定方法確定某站的高、中、低水位及相應(yīng)的測流要求，如王維志[2]等采用特征值法和漫灘水位指標(biāo)法對河道流量測驗中的水位級劃分方法進(jìn)行了研究，喇承芳[3]采用皮爾遜III型曲線法確定了黃河貴德水文站高、中、低水位，趙建偉[4]采用特征值法確定了沁河孔家坡水文站高、中、低水位及相應(yīng)的流量測驗方案。這些方法對開展流量測驗具有一定的指導(dǎo)意義，但仍沒有對水文測站流量測次的時空分布情況有個準(zhǔn)確合理的評價。本文基于小河新閘等站的流量測驗資料，依次采用累積頻率曲線法、適線檢驗法等方法，并利用概率學(xué)原理，逐步構(gòu)建了流量測驗時空分布合理性的評價方法，進(jìn)而為指導(dǎo)水文測站開展流量測驗工作，建立更為科學(xué)合理的流量測驗成果質(zhì)量評價體系，修訂和完善《河流流量測驗規(guī)范》提供參考。

2 流量測次分布區(qū)間及理想分布

為了便于描述，本文僅對形如：Q=f（z）的水位流量單一曲線（假定其已通過三項檢驗，且高、低水部分延長均滿足規(guī)范要求）進(jìn)行分析。

設(shè)某水文測站當(dāng)年實(shí)測N次流量（從小到大排列：Q1，Q2，…，QN-1，QN）。其中，實(shí)測最大、最小流量分別為Qmax、Qmin，流量變幅為ΔQ=Qmax-Qmin，則流量測次分布區(qū)間為[Qmin，Qmax]，每個子區(qū)間的步長為：

N個子區(qū)間及其集合Φi定義為：

若某一子區(qū)間內(nèi)有流量測次分布，則記為1，即：

顯然，最理想的流量測次分布是上述每個子區(qū)間內(nèi)均有流量測次，即：ΣM（N）=N。

當(dāng) ΣM（N）=N 時，亦有 ΣM（i）=i，則認(rèn)為流量測次分布為理想分布。若流量增幅是均勻的（Q2-Q1=Q3-Q2=…=QN-QN-1），則認(rèn)為流量測次分布為完全理想分布。

3 合格率與不合理

定義水文測站當(dāng)年流量測次分布的合格率p（%）：

可按水文站的精度要求，對流量測次分布的合格率要求如表2規(guī)定：

表2 p（%）指標(biāo)限定表

因此，當(dāng)水文測站當(dāng)年流量測次分布的合格率不滿足表2規(guī)定時，可認(rèn)為流量測次分布是不合理的。另外，令：

當(dāng)水文測站當(dāng)年流量測次分布的大于某一指標(biāo)時，亦可認(rèn)為流量測次分布是不合理的。指標(biāo)同樣按水文站的精度要求，如表3規(guī)定：

表 3 δ（%）指標(biāo)限定表

從上面定義來看，合格率p（%）規(guī)定了集合Φi中有測點(diǎn)分布的子區(qū)間的最少個數(shù)要求，δ（%）規(guī)定了流量測次分布中相鄰測點(diǎn)的最大允許級差要求（最大允許極差占當(dāng)年實(shí)測流量最大幅度的百分比不能超過限定要求）。顯然，合格率p（%）越大，δ（%）越小，流量測次的時空分布越相對合理。但合格率p（%）和δ（%）兩個指標(biāo)只是具體規(guī)定了流量測次分布是否合理的基本要求，仍無法準(zhǔn)確診斷流量測次空間的具體分布情況，如仍無法判斷水位流量關(guān)系點(diǎn)據(jù)集中在高水部分還是低水部分。

圖1為某三類精度要求的水文測站當(dāng)年16個測點(diǎn)所率定的水位流量關(guān)系曲線，已通過三項檢驗，且高、低水延長滿足流量推算要求。

圖1 某站水位流量關(guān)系圖

但從圖1可以看出，水文測站當(dāng)年16個測點(diǎn)所確定的16個子空間中，有7個子空間無測點(diǎn)分布，合格率p（%）=56.2%，δ（%）=25%，均不滿足表2和表3的要求，仍可以認(rèn)為該站流量測次分布是不合理的。

4 累積頻率曲線法

以下僅討論p（%）和δ（%）指標(biāo)滿足表2和表3規(guī)定的情形。

根據(jù)式（2）的規(guī)定，當(dāng)ΣM（N）＜N時，則涉及到有哪些子區(qū)間（記為Ωj，j=1，2，…，N-M 顯然，ΩjΦi，i=1，2，…，N）沒有測次分布（M（j）=0），這些區(qū)間在N個子區(qū)間的分布如何，如圖1所示，第2～5、8、13、15子區(qū)間無測點(diǎn)分布。

累積頻率曲線法是指：對于一組樣本x1，x2…，xn給定某一閾值x0，不大于x0的樣本數(shù)為m，則稱m/n為不大于x0的累積頻率；依次計算就得到累積頻率曲線。因此，建立ΣM（i）～i曲線，近似為累積頻率曲線，且有3種主要分布特性，如圖2所示（橫坐標(biāo)表示i，縱坐標(biāo)表示 ΣM（i））。

圖 2 i～ΣM（i）關(guān)系圖

從圖2中可以看出：

（1）當(dāng)曲線上凸時（I號線），則表示水位流量關(guān)系點(diǎn)據(jù)主要位于中下部；某一段曲線上凸越明顯，則表示點(diǎn)據(jù)越相對集中。

（2）當(dāng)曲線（II號線）下凹時，則表示水位流量關(guān)系點(diǎn)據(jù)主要位于中上部；某一段曲線下凹越明顯，亦表示點(diǎn)據(jù)越相對集中。

（3）當(dāng)曲線（III號線）下凹又上凸（或上凸又下凹）時，則表示水位流量關(guān)系中部點(diǎn)據(jù)相對缺乏。

另外，參照《水文資料整編規(guī)范》（SL247-2012），計算ΣM（i）關(guān)于i的標(biāo)準(zhǔn)差（Se）和隨機(jī)不確定度（X'Q）：顯然，標(biāo)準(zhǔn)差和隨機(jī)不確定度在一定指標(biāo)范圍內(nèi)，可認(rèn)為流量測次分布是相對合理的。標(biāo)準(zhǔn)差和隨機(jī)不確定度的指標(biāo)可按水文站的精度要求如表4。

表4 流量測次分布的標(biāo)準(zhǔn)差和隨機(jī)不確定度指標(biāo)表

從上面分析可以看出，累積頻率曲線法能在一定程度反映流量測次的大致分布情況，尤其通過標(biāo)準(zhǔn)差和隨機(jī)不確定度的計算能夠較好地確定流量測次分布的系統(tǒng)偏差，即：確保點(diǎn)據(jù)不會集中在某一（部分）子區(qū)間。

根據(jù)圖2，計算曲線I、II、III的標(biāo)準(zhǔn)差結(jié)果如表5（按三類精度要求）。

表5 標(biāo)準(zhǔn)差計算成果表

5 適線檢驗法

以下繼續(xù)討論Se（%）指標(biāo)滿足表4規(guī)定的情形。

為進(jìn)一步分析Ωj在Φi中的分布情況，可建立新的集合空間Ψk，Ψk定義如下：

式（7）中，Φ'm表示Φi中與Ωj相鄰的且M（i）=1的子區(qū)間。顯然，Ψk是由0、1兩個元素組成。當(dāng)Ψk中0、1依次排列或排列均勻時，可以認(rèn)為流量測次分布是相對合理的。合理性的判斷方法可借用《水文資料整編規(guī)范》（SL247-2012）中的適線檢驗法。主要步驟如下：

第一步，對 Ψk（=｛Ψ（1），Ψ（2），…，Ψ（k），｝）中已經(jīng)排列好的元素（0、1），從第二個元素開始，統(tǒng)計元素變換，變換符號記1，否則記0。即：

第二步，統(tǒng)計記為“1”的次數(shù)（ΣF（i））。按式（9）計算u值，并與給定的顯著性水平∝及u1-∝（見表6）進(jìn)行比較。當(dāng)u＜u1-∝時，則接受檢驗，流量測次空間分布合理，否則認(rèn)為流量測次空間分布不合理。

表6 臨界值特性表

顯然，通過使用適線檢驗法，在一定程度上構(gòu)建了流量測次分布合理性的評價方法。但對于連續(xù)無測點(diǎn)的子空間的分布要求，還沒有明確的判斷方法。因此，本文繼續(xù)從概率統(tǒng)計角度進(jìn)行分析。主要步驟如下：

第一步，設(shè)定Ψk中連續(xù)為0的元素個數(shù)不超過y。根據(jù)前面流量測次分布，此處y滿足：

第二步，計算Ψk中連續(xù)為0的元素個數(shù)不超過y的可能性py。此處類似：“將J個球（Ψk中為0的元素）放入K個盒子（Ψk中元素的個數(shù)），要求每個盒子至多放y個球”，根據(jù)排列組合與概率學(xué)原理[5]，則有：

第三步，當(dāng)py（%）小于某一給定指標(biāo)pN（%）時，可認(rèn)為流量測次分布是合理。同樣，按水文站的精度要求，pN（%）如表7定義：

表7 pN（%）指標(biāo)限定表

6 實(shí)際應(yīng)用

（1）小河新閘站

小河新閘站為三類精度要求的水文站。2017年2～9月實(shí)測引（長江）水流量16次，10月實(shí)測4次。實(shí)測點(diǎn)據(jù)按流量從小到大排列，如表8：

表8 小河新閘站2017年實(shí)測引水流量成果表

根據(jù)表8中2～9月實(shí)測引水流量點(diǎn)據(jù)，采用一潮推流法率定水位流量關(guān)系為：

三項檢驗成果如表9。

表9 小河新閘站水位流量關(guān)系三項檢驗成果表

（續(xù)表9）

計算該站合格率p（%）=56.2%，δ（%）=25%，均不滿足表2和表3的要求。為此，2017年10月該站又補(bǔ)測了4潮次引水流量（見表8）。經(jīng)計算，該站水位流量關(guān)系沒有發(fā)生改變，也通過三項檢驗（略）。此時，總測次N=20，p（%）=60%，δ（%）=13.2%，均滿足表2和表3的要求，i～ΣM（i）如圖3所示。

圖3 小河新閘站關(guān)系圖

從圖3中可以看出，該站實(shí)測水位流量點(diǎn)據(jù)主要集中在中下部，且標(biāo)準(zhǔn)差Se（%）=3.6%，滿足表4的要求。另外，該站由元素（0、1）組成的集合空間Ψk如表10所示。

表10 小河新閘站Ψk元素表（K=14）

根據(jù)表9和式（9），計算Ψk中元素變換次數(shù)ΣF（h）=10＞7（K/2），通過適線檢驗（亦滿足表 6的要求）。且該站Ψk中連續(xù)為0的元素個數(shù)不超過2，亦滿足式（10）的要求。同樣，根據(jù)式（11）和式（12），計算py=0.02%，滿足表7的要求。

由此可見，經(jīng)過10月有針對性的流量補(bǔ)測，小河新閘水文站流量測次滿足了空間合理分布的要求，全年流量資料整編精度得到了保證。

（2）楓橋站

楓橋站為一類精度要求的水文站。根據(jù)該站2016年實(shí)測流量資料，采用上述評價方法，評價結(jié)果如表11所示。

表11 楓橋站2016年流量測驗時空分布合理性評價結(jié)果表

7 結(jié)語

流量測驗成果質(zhì)量是水文資料可靠性的保證，而流量測驗時空分布（測次的布置和測點(diǎn)的控制）則是流量測驗成果質(zhì)量的基礎(chǔ)?！逗恿髁髁繙y驗規(guī)范》對流量測驗次數(shù)的布置和測點(diǎn)的控制有一個定性的規(guī)定，《水文資料整編規(guī)范》和《水文年鑒匯編刊印規(guī)范》則對水位流量關(guān)系定線成果有一個定量的規(guī)定。但由于各地河流的差異性，目前我國對流量測驗時空分布的合理性尚缺乏明確的評價方法。

本文基于累積頻率曲線法、適線檢驗法等方法，并利用概率學(xué)原理，較為嚴(yán)密地推導(dǎo)和構(gòu)建了流量測次時空分布合理性的評價辦法。根據(jù)小河新閘水文站2017年的流量資料和楓橋水文站2016年的流量資料，驗證效果較為明顯，可為指導(dǎo)水文測站開展流量測驗工作、建立科學(xué)合理的流量測驗成果質(zhì)量評價體系、修訂和完善《河流流量測驗規(guī)范》提供參考。

需要說明一點(diǎn)的是，本文所推導(dǎo)和構(gòu)建的方法僅僅是對流量測驗時空分布合理性評價的一個嘗試。表2、3、4、6中不同指標(biāo)的限定值僅基于十幾個水文站近20年水文資料的驗證結(jié)果，需要用更多的水文資料去驗證和確定。另外，子空間是隨測次變化的，從在小河新閘站上的驗證情況來看，測次增多不一定使流量測次時空分布的合格率增大，對子空間的定義也可能需要更多的水文資料去驗證和確定。