張繼東,康青勇,唐 斌
(四川省雅安水文水資源勘測中心,四川 雅安 625000)
隨著水文站網(wǎng)不斷建設(shè)完善,水文自動(dòng)化測驗(yàn)在減少人為錯(cuò)誤和節(jié)省人力資源方面發(fā)揮著重要作用。近年來隨著時(shí)差法測流技術(shù)的逐步應(yīng)用和推廣,時(shí)差法已廣泛應(yīng)用于大洪水及惡劣環(huán)境流量測驗(yàn),有效提高了測量精度及工作效率[1-2]。陳武等[3]通過分析時(shí)差法測流系統(tǒng)在劉老澗泵站的應(yīng)用,表明該監(jiān)測系統(tǒng)具有自動(dòng)化和信息化優(yōu)勢,具備較高的應(yīng)用與推廣價(jià)值;段炎沖等[4]研究表明在含沙條件下,超聲時(shí)差法的流量測量值總體偏大,且偏差量隨含沙量增大呈增大趨勢;李正最等[5]采用二層測速方式構(gòu)建模型進(jìn)行流量計(jì)算參數(shù)的率定,取得的流量成果精度能夠滿足規(guī)范要求。本文對蘆山水文站時(shí)差法測流系統(tǒng)進(jìn)行數(shù)據(jù)比測,探討系統(tǒng)在該站的試驗(yàn)成果是否合理可靠,促進(jìn)系統(tǒng)投產(chǎn)應(yīng)用,針對性解決蘆山水文站水文監(jiān)測工作手段單一的問題。
蘆山河發(fā)源于蘆山縣與大邑縣交界處的羊子嶺大雪塘峰。上源稱黑水河,納大尖峰及大雪峰所出眾山溪,南流至中嘴,右納黃水河,以下稱玉溪河。入蘆山縣境又叫蘆山河。河長113 km,流域面積1 397 km2。
流域內(nèi)以丘狀高原為主,海拔高程在3 000~4 000 m,地表丘狀起伏,丘頂渾圓,谷地寬闊,低洼處多為沼澤、草甸。丘谷間相對高差不大,一般在50~100 m。流域內(nèi)土壤以深棕壤為主,間以褐土、腐殖土等。自然植被組成群落以草本植物為主,覆蓋較淺。
蘆山水文站地處四川省雅安市蘆山縣蘆陽鎮(zhèn)黎明村,坐標(biāo)為東經(jīng)102°55′50.9″,北緯30°07′38.1″,于2002年5月1日建站。該站是長江流域大渡河青衣江水系青衣江上源寶興河左岸支流蘆山河控制站,集水面積1 381 km2,至河口4 km,屬省級重要水文站、二類精度流量站,測驗(yàn)項(xiàng)目為水位、流量、降水。
蘆山河流域徑流主要由降雨所形成,由于該流域植被較好,使下墊面具有較強(qiáng)的滯留能力,枯期、汛前過渡期有少量融雪水和地下水補(bǔ)給,因此,地下徑流較豐沛,徑流年際變化小,年內(nèi)分配均勻。其年內(nèi)分配主要情況:汛期(5-10月)徑流占全年徑流的85.8%,主汛期(6-9月)徑流占全年徑流的71.5%。
根據(jù)蘆山水文站多年水文資料統(tǒng)計(jì),蘆山河多年年平均流量均值21.5 m3/s,最大年平均流量37.2 m3/s(2020年),最小年平均流量12.4 m3/s(2009年),最大與最小相比為3.0,年際間變化較大,月年平均徑流量統(tǒng)計(jì)見表1。
表1 蘆山水文站月年平均徑流量及年內(nèi)分配
蘆山河低水較穩(wěn)定,中高水流速偏大。多受上游來水、本地暴雨影響,且具有峰高、量大的特點(diǎn)。洪水期,坡陡、流急、浪大,又因洪水夾帶大量泥沙水草和漂木等。洪水漲落變化迅速,浪大水急。根據(jù)洪水特性屬于暴漲暴落洪水,洪峰歷時(shí)短,約10 min至1h,洪水過程一般為1~2d,并且一般出現(xiàn)洪峰時(shí)間在凌晨4-9時(shí)。蘆山河洪水季節(jié)的變化與暴雨時(shí)程分布基本一致,洪水發(fā)生在6-9月,年最大流量一般在7-9月。最大洪峰流量各月出現(xiàn)次數(shù)統(tǒng)計(jì)見表2。
表2 蘆山水文站年最大洪峰流量各月出現(xiàn)次數(shù)統(tǒng)計(jì)
超聲波時(shí)差法的測量原理[6]基于兩個(gè)超聲波傳感器之間傳播的聲信號持續(xù)時(shí)間的直接測量。這種方法敘述如下:至少安裝兩個(gè)超聲波傳感器,從而使它們之間傳播的聲信號與水流成對角線。對著水流持續(xù)的聲信號其傳播時(shí)間長于沿著水流持續(xù)的聲信號。傳播時(shí)間差與流速成正比,并且如果知道橫截面,則與流量成正比。
可通過頻率法和脈沖法兩種方法測量傳播時(shí)間。如果使用頻率法,將通過河道在發(fā)送和接收傳感器之間發(fā)送一系列既定頻率信號,這些信號的傳播時(shí)間將被視為測量對象;如果使用脈沖法,將在既定頻率下測量短聲脈沖的傳播時(shí)間。使用脈沖法時(shí)的傳播時(shí)間原理具有可以連續(xù)測量的優(yōu)點(diǎn)。正因?yàn)檫@個(gè)優(yōu)點(diǎn),這一傳播時(shí)間原理已被成功用于河流、渠道和管道的超聲波流量測量。這是Quantum使用的方法,如今,已被水文學(xué)家認(rèn)可為最具有價(jià)值的方法。
為了確定傳播時(shí)間,超聲波傳感器接收方波信號(正弦波)形式的聲信號;然后,傳感器將通過河道發(fā)送的這些信號轉(zhuǎn)變?yōu)槁暡ㄈ?最后,對該波群的傳播時(shí)間進(jìn)行測量。傳播時(shí)間t1-2和t2-1是通過直接傳播時(shí)間法確定的,而制造商以前憑借模擬或數(shù)字信號識別器確定的零交叉是計(jì)算得出的。使用傳播時(shí)間法時(shí),時(shí)間測量必須以毫微秒為單位,從而在所需精度下獲得最小的水流速度。
蘆山水文站超聲波時(shí)差法設(shè)備由德國的Quantum Hydrometrie GmbH公司生產(chǎn),為單路徑系統(tǒng),有線主機(jī),安裝于蘆山水文站基本水尺斷面處,在基本水尺斷面上下游兩岸分別安裝一個(gè)換能器,詳細(xì)參數(shù)見表3。
表3 主要技術(shù)參數(shù)
蘆山水文站超聲波時(shí)差法設(shè)備2020年6月10日安裝調(diào)試完成,收集數(shù)據(jù)到2020年8月10日止,共收集11 000多條數(shù)據(jù)。2020年8月11日被洪水沖毀,由于廠家配件原因,還未恢復(fù)。本次分析采用安裝調(diào)試完成后2020年6月10日至8月10日的數(shù)據(jù),施測涵蓋水位級797.43~799.54 m,共2.11 m的變幅,流量級3.31~340 m3/s。
根據(jù)收集的數(shù)據(jù),超聲波時(shí)差法每5 min一個(gè)數(shù)據(jù),對實(shí)測數(shù)據(jù)按照不同水位進(jìn)行提取,得到有代表性的實(shí)測數(shù)據(jù)系列,超聲波時(shí)差法為兩組系列,分別是6月10日至6月27日和6月27日至8月10日。用轉(zhuǎn)子流速儀施測的流量繪制一條水位流量關(guān)系線,2020年4月3日至9月26日蘆山水文站水位流量關(guān)系線為同一系列,為單一線。將3條水位流量關(guān)系曲線繪制在同一張圖上,見圖1所示。
圖1 蘆山站水位流量關(guān)系線
從圖1可以看出,6月27日前一個(gè)系列,6月27日后一個(gè)系列。6月27日前與轉(zhuǎn)子流速儀施測流速基本一致,但6月27日后較轉(zhuǎn)子流速儀數(shù)據(jù)偏大。
根據(jù)收集的資料,經(jīng)計(jì)算機(jī)線性函數(shù)回歸分析點(diǎn)繪斷面實(shí)測流量與超聲波時(shí)差法相關(guān)關(guān)系圖,見圖2和圖3。從圖2和圖3,得出相關(guān)方程:6月27日前Y=1.07×X,相關(guān)系數(shù)R2=0.986,為顯著相關(guān),相關(guān)關(guān)系較好;6月27日后Y=0.82×X,相關(guān)系數(shù)R2=0.992,為顯著相關(guān),相關(guān)關(guān)系較好。其中Y為流速儀實(shí)測流量,即斷面流量;X為超聲波時(shí)差法流量。
圖2 6月27日前實(shí)測流量與超聲波時(shí)差法流量相關(guān)關(guān)系
圖3 6月27日后實(shí)測流量與超聲波時(shí)差法流量相關(guān)關(guān)系
超聲波時(shí)差法實(shí)測流量點(diǎn)水位流量與實(shí)測水位流量關(guān)系曲線水位查流量比值,6月27日前比值平均值為1.10;6月27日后為比值平均值為0.82,見表4和表5。
表4 6月27日前水位流量結(jié)點(diǎn)比值
表5 6月27日后水位流量結(jié)點(diǎn)比值
6月27日前和6月27日后與線性函數(shù)相關(guān)關(guān)系1.07、0.82接近,故成果可靠。
(1)通過比測分析,本次超聲波時(shí)差法雖為兩個(gè)系列,但6月10日安裝完成到6月27日之間無漲水過程,且時(shí)間較短,6月27日后該系列較穩(wěn)定代表性好,故系數(shù)率定成果采用0.82。
(2)根據(jù)原始數(shù)據(jù)分析,6月25日23時(shí)開始起漲到6月27日0時(shí)段無流速,6月27日0時(shí)(該時(shí)段為退坡中間)后流速數(shù)據(jù)又恢復(fù),蘆山站漲坡含沙量較大,說明超聲波時(shí)差法受水體含沙量影響較大。
(3)蘆山站屬山區(qū)性河流,洪水暴漲暴落,含沙量大,本次采用的超聲波時(shí)差法工作頻率為200 kHz/28 kHz ,該工作頻率對于天然河道流量測驗(yàn)工作頻率相對較低,不足以使兩換能器波束穿透含沙量較大時(shí)的水體。
(4)由于超聲波時(shí)差法自身不能對水深進(jìn)行測量,需結(jié)合水位和斷面數(shù)據(jù)對流量進(jìn)行計(jì)算,兩岸安裝的換能器還需入水一定深度,因此,其對斷面變化要求較高,測驗(yàn)河段盡量順直、沖淤變化小、有一定的水深。
(5)超聲波時(shí)差法目前正在試用階段,設(shè)備的安裝對系數(shù)影響很大,特別是其換能器安裝于水下,有可能受漂浮物等的撞擊,故系數(shù)有待在以后的使用中不斷修正。