（山東省煙臺(tái)市水文局 煙臺(tái) 264001）

1 概況

臧格莊水文站位于清洋河流域系大沽夾河西支流，山東省棲霞市臧家莊鎮(zhèn)臧家莊村，東經(jīng)120°59′，北緯37°28′。流域形狀呈扇形，斷面以上流域面積458km2，流域長(zhǎng)度30.5km，流域平均寬度15.0km。干流長(zhǎng)度28.7km，干流平均坡度4.33‰。上游20km 處1961年1月建成庵里中型水庫(kù)一座，基本斷面以下約300m 有山東河匯入，下游約28km有門(mén)樓水庫(kù)水文站。

該流域?qū)贃|亞季風(fēng)區(qū)大陸性氣候，氣候溫和，四季分明。多年平均氣溫11.5℃，七、八月份最高，達(dá)38℃，最低溫度-17.3℃，年無(wú)霜期222 天，歷年最大凍土深0.5m，平均風(fēng)速4.5m/s。多年平均蒸發(fā)量1186.5mm，4—6月達(dá)475.7mm，占全年蒸發(fā)量的40.1%。氣候特點(diǎn)是：春季常刮大風(fēng)，多以西南風(fēng)為主，空氣相對(duì)干燥；夏季麥?zhǔn)涨昂蟀橛斜?zāi)害天氣，汛期濕熱多雨且降水集中，時(shí)有臺(tái)風(fēng)登陸，形成暴雨洪水；秋季天氣涼爽，風(fēng)向較亂，個(gè)別年份甚至有連綿陰雨發(fā)生；冬季干冷、雨雪稀少，多北風(fēng)或西北風(fēng)。

臧格莊水文站流域降水年內(nèi)分配不均勻，多年平均降雨量為667.6mm，年最大降水量為1177.7mm，年最小降水量為344.0mm，極值比3.4。臧格莊水文站的徑流主要由大氣降水補(bǔ)給，上游中型水庫(kù)庵里水庫(kù)流域面積較小，故該站徑流隨時(shí)間的變化過(guò)程十分復(fù)雜，尤其在汛期，徑流變化特點(diǎn)與降水相似。汛期6—9 四個(gè)月多年平均徑流深占全年的80%。從嚴(yán)格意義上來(lái)說(shuō)，臧格莊水文站當(dāng)上游未發(fā)生有效降雨或水庫(kù)放水時(shí)，河道年際分布均處于枯水段。

2 枯季流量測(cè)驗(yàn)發(fā)現(xiàn)的問(wèn)題

取2018年臧格莊水文站枯水期部分實(shí)測(cè)流量資料進(jìn)行分析（見(jiàn)表1），發(fā)現(xiàn)18 號(hào)測(cè)次水位與1 號(hào)測(cè)次水位相同，但實(shí)測(cè)流量較1 號(hào)測(cè)次增加了33%，分析兩次流量測(cè)驗(yàn)過(guò)程可知，兩次流量測(cè)驗(yàn)均在基下250m 處涉水施測(cè)，平均水深及斷面面積等要素未發(fā)生變化，由于處于冬季枯水期，河道無(wú)沖於、水生植物、結(jié)冰、人類(lèi)活動(dòng)等導(dǎo)致河道變化的原因，發(fā)現(xiàn)斷面平均流速增大了33%，得到了流量變大的主要原因。其他要素未發(fā)生變化，而斷面平均流速發(fā)生變化，需進(jìn)一步進(jìn)行分析。

增加14、17 測(cè)次與1、18 測(cè)次進(jìn)行共同比對(duì)分析，繪制水位流量關(guān)系線如圖1所示。根據(jù)臧格莊站實(shí)測(cè)流量成果表，14—16 測(cè)次屬于平水段，17、18 測(cè)次屬于落水段，故17 測(cè)次水位應(yīng)大于18 測(cè)次水位，這從兩次流量數(shù)據(jù)可以佐證。但17 測(cè)次與之前穩(wěn)定流的5.13m 水位時(shí)流量相差較大（20%），故17 測(cè)次水位應(yīng)大于5.12m 但未到達(dá)5.13m，暫記為5.127m；18 測(cè)次水位應(yīng)大于5.12m 但小于5.127m，暫記為5.123m，繪制水位流量關(guān)系線如圖2所示（此處應(yīng)具體分析，臧格莊水文站P3 水尺零點(diǎn)為4.000m，若水尺零點(diǎn)毫米位不為零，應(yīng)按照原水位取舍后再加水尺讀數(shù)得到調(diào)整水位）。

表1 臧格莊水文站枯水期部分實(shí)測(cè)流量特征值表

將17、18 測(cè)次水位變化后，趨勢(shì)線由圖1變?yōu)閳D2，可見(jiàn)線性更加符合實(shí)際情況并且相關(guān)性較高。進(jìn)一步調(diào)整17、18 測(cè)次水位毫米位，得到圖3。

可見(jiàn)，當(dāng)17 測(cè)次水位記為5.126m，18 測(cè)次水位記為5.124m 時(shí)，水位與流量完全相關(guān)，在無(wú)其他條件影響下較為符合實(shí)際情況。由圖1至圖3的變化過(guò)程不難看出，水位的取舍值在枯季測(cè)流過(guò)程中對(duì)流量數(shù)據(jù)合理性檢驗(yàn)起到了一定的作用。

以臧格莊水文站多年定線的經(jīng)驗(yàn)來(lái)年，中高水為單一線，符合定線要求，低水由于水位流量關(guān)系差而采用連實(shí)測(cè)過(guò)程線，流量施測(cè)頻次為3～5 天。通過(guò)上文分析幾個(gè)測(cè)次的實(shí)測(cè)數(shù)據(jù)不難發(fā)現(xiàn)，存在水位精讀不夠而導(dǎo)致的誤差，在進(jìn)行水位精讀的調(diào)整后可以提高水位、流量的相關(guān)性。

3 各要素對(duì)測(cè)驗(yàn)精度的影響

3.1 水位

由于水位讀數(shù)采用的水尺牌為厘米級(jí)，精讀厘米，估讀毫米，所以在實(shí)際水位觀讀時(shí)容易出現(xiàn)讀數(shù)精讀偏差。像調(diào)整后的1.126 實(shí)際讀數(shù)為1.13，1.124 讀數(shù)為1.12，實(shí)際讀數(shù)雖然符合水位觀測(cè)規(guī)范，卻在實(shí)際定線過(guò)程中產(chǎn)生相關(guān)性不高而采用的連實(shí)測(cè)定線方法，使得測(cè)站的工作量加大。

臧格莊水文站基本斷面下游220m有攔河閘一處，攔蓄水量，加之清洋河河道多年無(wú)較大洪水，致使河道內(nèi)蘆葦叢生，斷面死水面積較大。根據(jù)2018年實(shí)測(cè)大斷面資料，水位5.12m 時(shí)，基本斷面面積為102.67m2，水位5.13m 時(shí)為103.88m2（為進(jìn)行數(shù)據(jù)分析，故未按要求進(jìn)行有效位數(shù)的保留），二者相差1.21 m2。這部分基本斷面的面積偏差既可以表現(xiàn)在下游臨時(shí)測(cè)流斷面的水深上，也可以表現(xiàn)在流速上。

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

圖2 水位流量關(guān)系線2 圖

水位是繪制水位流量關(guān)系曲線的主要數(shù)據(jù)，所以測(cè)流時(shí)的水尺觀測(cè)精度對(duì)計(jì)算水位非常重要，應(yīng)盡可能在水尺牌讀數(shù)位于整厘米數(shù)時(shí)實(shí)測(cè)流量。

3.2 水深

枯季臧格莊水文站河道由于攔水閘、水草等因素影響，導(dǎo)致基本斷面死水面積大，不具備測(cè)流條件。在下游涉水施測(cè)流量，同水位觀測(cè)情況相同，測(cè)深桿刻度同樣為1cm，在有水流時(shí)更易出現(xiàn)水深測(cè)量的讀數(shù)精讀偏差。根據(jù)某實(shí)測(cè)流量成果，若當(dāng)次測(cè)深垂線水深測(cè)量值較實(shí)際值均偏大1cm 時(shí)，相對(duì)于斷面原流量數(shù)據(jù)偏大3.3%。

3.3 測(cè)點(diǎn)定位

枯水期河道水深、流速較小，進(jìn)行測(cè)速垂線測(cè)點(diǎn)定位時(shí)應(yīng)更加注意測(cè)點(diǎn)位置準(zhǔn)確，這會(huì)對(duì)流量的測(cè)驗(yàn)精度產(chǎn)生較大影響。

天然河道暢流期常見(jiàn)的垂線測(cè)點(diǎn)流速v 分布曲線示意圖如圖4所示。

圖4中水深為相對(duì)水深，即各點(diǎn)的水深h 歸一化成點(diǎn)位水深h 與總水深H 的比值。

一般情況下，暢流期時(shí)垂線水面流速最大，河底附近流速為零，總體垂線流速分布呈一定形狀的曲線。影響流速曲線形狀的因素有很多，致使垂線流速分布曲線的形狀多種多樣。人們總是努力用曲線函數(shù)近似描述垂線的流速分布，本例采用橢圓流速分布曲線。測(cè)流垂線上的流速分布用橢圓流速分布表示為：

表2 橢圓流速分布公式計(jì)算的測(cè)點(diǎn)流速表

圖3 水位流量關(guān)系線3 圖

圖4 水面無(wú)阻暢流流速沿垂線分布圖

式中：v0—水面流速（η=0），m/s；

P—流速分布參數(shù)，取P=0.6，相當(dāng)于謝才系數(shù)C=40～60；

η—由水面向下起算的相對(duì)水深；

v—相對(duì)于水深為時(shí)的點(diǎn)流速，m/s。

按式1 計(jì)算的幾個(gè)常用相對(duì)水深處的測(cè)點(diǎn)流速值見(jiàn)表2。

由式1 按積分法計(jì)算垂線平均流速：

水文最常采用的計(jì)算垂線平均流速的方法為兩點(diǎn)法與一點(diǎn)法，其計(jì)算結(jié)果vm=0.887（兩點(diǎn)法，相對(duì)水深0.2、0.8）；vm或=0.885（一點(diǎn)法，相對(duì)水深0.6）。由此可見(jiàn)通常采用的計(jì)算垂線平均流速的方法計(jì)算數(shù)值與積分法計(jì)算數(shù)值相近，說(shuō)明流速沿豎向分布對(duì)測(cè)速點(diǎn)位置的選擇有很大實(shí)用價(jià)值。由于測(cè)點(diǎn)位置的差別對(duì)于垂線平均流速的影響也是存在的，且測(cè)點(diǎn)流速定位誤差受水深誤差的綜合影響。

若實(shí)際測(cè)速垂線水深為0.30m 時(shí)，相對(duì)水深0.6m時(shí)為測(cè)點(diǎn)定位在0.18m，若精度不高定位在0.17m 時(shí)，實(shí)際相對(duì)水深為0.57m，則測(cè)點(diǎn)流速差距為0.897，與相對(duì)0.6m 水深誤差達(dá)14%，且測(cè)流斷面水深越小誤差越大。

3.4 起點(diǎn)距

由于臨時(shí)斷面沒(méi)有標(biāo)準(zhǔn)的起點(diǎn)距刻畫(huà)，在實(shí)際工作中大多為人工丈量起點(diǎn)距，該方法受角度誤差與垂線間距離誤差共同作用，若各測(cè)深垂線丈量出現(xiàn)0.05m 誤差，則流量數(shù)據(jù)偏差達(dá)到5.6%，當(dāng)斷面與河流流向不垂直，夾角為85°時(shí)，誤差將增大8.7%。

4 結(jié)語(yǔ)

對(duì)枯水期水位漲落較小的測(cè)站來(lái)說(shuō)，測(cè)時(shí)水位、水深、測(cè)點(diǎn)位置、測(cè)深垂線水深等因素均能導(dǎo)致該次流量偏離值較大，影響到標(biāo)準(zhǔn)差計(jì)算，這種情況就是當(dāng)河道水文站水位變化較小，但是流量卻相差較大的主要原因，這種誤差在水文站是經(jīng)常出現(xiàn)并且會(huì)影響枯水期測(cè)驗(yàn)精度進(jìn)而影響的水位流量關(guān)系曲線定線問(wèn)題，應(yīng)該引起重視■