梁文勝,梅 鹍
(長(zhǎng)沙水文水資源勘測(cè)中心,湖南 長(zhǎng)沙 410014)
長(zhǎng)沙水文站位于長(zhǎng)沙市湘江大道朱張渡下游80 m,東經(jīng)112°58′,北緯28°11′,距河口60 km,集水面積83 020 km2。測(cè)驗(yàn)河段由橘子洲頭分為左分河、右主河。上游4 000 m 以上建有猴子石大橋,下游2 370 m 處建有桔子洲大橋。主河床由卵石夾沙組成。左岸為花崗石護(hù)堤,右岸為鋼筋混凝土防洪墻。基下18 m 處有自記水位井一座,警戒水位36.00 m,保證水位38.37 m,水位達(dá)37.80 m 左右時(shí)桔子洲開(kāi)始漫灘。長(zhǎng)沙水文站為一類(lèi)精度站。歷年最低水位24.63 m(凍結(jié)基面),歷年最低水位24.63 m,歷年最高水位39.51 m(凍結(jié)基面),2019 年最高水位38.35 m。
本次比測(cè)的儀器為Decatur 手持式電波流速儀。長(zhǎng)沙水文站采用Workhouse R5 型ADCP 實(shí)測(cè)流量成果。R5 型平面相控陣變頻智能ADCP 采用頻率合成技術(shù),根據(jù)水深變化自動(dòng)調(diào)節(jié)工作頻率、自動(dòng)調(diào)節(jié)剖面單元大小,以適用大范圍水深變化的測(cè)流斷面。
長(zhǎng)沙水文站為湘江一類(lèi)精度水文站,日常水文測(cè)驗(yàn)項(xiàng)目為:水位、降水,流量項(xiàng)目根據(jù)防汛抗旱要求進(jìn)行施測(cè)。長(zhǎng)沙水文站常規(guī)流量測(cè)驗(yàn)方法為ADCP 測(cè)流,正常情況下可以進(jìn)行湘江流量測(cè)驗(yàn)。為確保高洪時(shí)段ADCP 等常規(guī)流量測(cè)驗(yàn)方法無(wú)法進(jìn)行的情況下,長(zhǎng)沙水文站能及時(shí)施測(cè)流量,為防汛工作提供及時(shí)準(zhǔn)確的水情信息,制定了長(zhǎng)沙水文站高洪流量測(cè)驗(yàn)方案。該方案采用電波流速儀在猴子石大橋上施測(cè)湘江表面流速,計(jì)算斷面虛流量。通過(guò)同步實(shí)測(cè)流量與虛流量建立相關(guān)關(guān)系,率定斷面虛流量系數(shù),由虛流量與虛流量系數(shù)之積,求得湘江長(zhǎng)沙站流量。
由于長(zhǎng)沙水文站基本水尺斷面為復(fù)式斷面,不適合ADCP 測(cè)流。為了提高流量精度,避開(kāi)瀏渭河和撈刀河入流對(duì)湘江流量的影響,ADCP 測(cè)流斷面選擇在湘江長(zhǎng)沙水文站基下7.5 km(瀏渭河入湘江口的上游400 m 處)。
ADCP 過(guò)河方式采用動(dòng)力船拖拽過(guò)河。ACDP 測(cè)流計(jì)算水位采用長(zhǎng)沙水文站水位。
電波流速儀高洪斷面位于長(zhǎng)沙水文站上游4.2 km(猴子石大橋上游50 m 處)。在高洪斷面起點(diǎn)距為150 m、300 m、450 m、600 m、750 m 處布設(shè)測(cè)速垂線(xiàn)。工作人員站在猴子石大橋上,手持Decatur 電波流速儀,面向湘江上游施測(cè)表面流速。猴子石大斷面實(shí)測(cè)成果如表1。
表1 高洪斷面成果表
2019 年7 月9~19 日,水位從31.72 m 上漲到38.35 m。湘江長(zhǎng)沙水文站采用高洪方案測(cè)流,用電波流速儀在猴子石大橋斷面共測(cè)得洪水過(guò)程22 次。
比測(cè)期間繪制長(zhǎng)沙水文站實(shí)測(cè)流量過(guò)程線(xiàn),為了提高長(zhǎng)沙水文站流量過(guò)程資料的合理性和分析精度,參照上游湘潭水文站實(shí)測(cè)流量過(guò)程線(xiàn)趨勢(shì),對(duì)長(zhǎng)沙水文站實(shí)測(cè)流量過(guò)程進(jìn)行修正。
為計(jì)算方便,提高高洪測(cè)流方案可行性、降低方案難度,高洪斷面電波流速儀流量計(jì)算直接采用長(zhǎng)沙水文站水位。雖然長(zhǎng)沙水文站基本水尺斷面和高洪斷面水位存在一定差值,但由于湘江高洪期間,水面落差小且變幅小,其落差影響反映在綜合虛流量系數(shù)中。
根據(jù)電波流速儀測(cè)流時(shí)間,查讀長(zhǎng)沙水文站連實(shí)測(cè)流量過(guò)程線(xiàn)上的流量。通過(guò)建立電波流速儀虛流量和實(shí)測(cè)流量過(guò)程線(xiàn)的同步流量?jī)烧咧g的相關(guān)關(guān)系,從而率定電波流速儀綜合虛流量系數(shù)K電。
7 月9~19 日間,長(zhǎng)沙水文站采用走航式ADCP 在湘江長(zhǎng)沙水文站基下7.5 km(瀏渭河口上游200 處)實(shí)測(cè)洪水過(guò)程16 次。實(shí)測(cè)最大流量23 900 m3/s,相應(yīng)水位38.28 m。實(shí)測(cè)最小流量5 650 m3/s,相應(yīng)水位33.82 m。為了消除測(cè)驗(yàn)誤差,提高長(zhǎng)沙水文站分析流量精度,對(duì)長(zhǎng)沙水文站實(shí)測(cè)流量過(guò)程線(xiàn)中與湘潭水文站的實(shí)測(cè)過(guò)程趨勢(shì)不一致的局部流量過(guò)程,按湘潭水文站流量過(guò)程趨勢(shì)進(jìn)行插補(bǔ),使兩者趨勢(shì)一致。
根據(jù)電波流速儀測(cè)流時(shí)間,查讀長(zhǎng)沙水文站連實(shí)測(cè)流量過(guò)程線(xiàn)上的流量,共得22 組數(shù)據(jù),見(jiàn)表2。
表2 電波流速儀虛流量和查線(xiàn)流量表 m3/s
點(diǎn)繪長(zhǎng)沙水文站高洪方案施測(cè)的電波流速儀虛流量和相應(yīng)時(shí)間的長(zhǎng)沙水文站查線(xiàn)流量數(shù)據(jù)圖發(fā)現(xiàn),二者之間數(shù)據(jù)呈帶狀。兩組數(shù)據(jù)進(jìn)行直線(xiàn)擬合,電波流速儀虛流量系數(shù)為0.610 3,其相關(guān)關(guān)系R2為0.982 1。電波流速儀虛流量與查線(xiàn)流量擬合相關(guān)關(guān)系見(jiàn)圖1。
圖1 長(zhǎng)沙水文站流量~電波流速儀虛流量關(guān)系圖
誤差分析:ADCP 測(cè)流斷面與長(zhǎng)沙水文站高洪方案中電波流速儀測(cè)流斷面區(qū)間匯入的較大支流僅有龍王港。龍王港流域面積154 km2,占長(zhǎng)沙水文站控制流域面積0.19%。本次電波流速儀系數(shù)分析比測(cè)采用2019 年7 月9~17 日特大洪水過(guò)程的實(shí)測(cè)數(shù)據(jù),期間龍王港流域代表站梅溪湖站降雨量29 mm,最大日降雨量13 mm。龍王港流域降水少,受湘江頂托,匯入湘江流量小,故其匯入流量可以忽略不計(jì),可以近似認(rèn)為電波流速儀測(cè)流斷面和長(zhǎng)沙水文站ADCP 測(cè)流斷面流量相等。
根據(jù)公式Q=0.61×Q虛+7 444,計(jì)算電波流速儀流量。將電波流速儀流量與查線(xiàn)流量進(jìn)行分析,進(jìn)行三線(xiàn)檢驗(yàn),結(jié)果見(jiàn)表3。最大誤差為3.5%,最小誤差為0.0%,標(biāo)準(zhǔn)差2.0%,系統(tǒng)誤差0.0%。電波流速儀流量誤差滿(mǎn)足規(guī)范要求。
1)電波流速儀實(shí)測(cè)流量與長(zhǎng)沙水文站查線(xiàn)流量關(guān)系符合直線(xiàn)關(guān)系
最大誤差為3.47%,最小誤差為0.01%,標(biāo)準(zhǔn)差2.0%,系統(tǒng)誤差0.0%。比測(cè)水位變幅占?xì)v年變幅的92%。原始分析資料準(zhǔn)確、可靠,比測(cè)成果合理。
2)長(zhǎng)沙水文站流量計(jì)算公式
Q=0.61×Q虛+7 444
式中 Q——湘江長(zhǎng)沙水文站流量(m3/s);
Q虛——電波流速儀虛流量(m3/s)。
3)測(cè)驗(yàn)條件要求
電波流速儀流量測(cè)驗(yàn)建議避開(kāi)大風(fēng)、大雨期。本方案施測(cè)流量需滿(mǎn)足以下條件:
表3 長(zhǎng)沙水文站電波流速儀流量關(guān)系曲線(xiàn)檢驗(yàn)
①長(zhǎng)沙水文站流量大于14 000 m3/s(適用高于警戒水位36.00 m)。
②電波流速儀測(cè)流斷面在長(zhǎng)沙站基上4.2 km。
③流量測(cè)驗(yàn)方案按5 根測(cè)速垂線(xiàn)均勻分布。
④采用電波流速儀面向湘江上游施測(cè)湘江流速。
⑤斷面、流量計(jì)算均采用長(zhǎng)沙水文站水位(凍結(jié)基面)。