吳志勇,徐 梁,唐運憶,何 海
(1.河海大學水文水資源學院,江蘇 南京 210098; 2.江蘇省水資源服務中心,江蘇 南京 210029)
流量是單位時間內(nèi)通過江河某一橫斷面的水體體積,是反映江河等水體水量變化的基本數(shù)據(jù),也是河流最重要的水文特征值[1-2]。傳統(tǒng)的流量測驗方法是將測驗斷面分成多個部分,并在每部分上布設(shè)測速垂線,利用流速儀等設(shè)備測量各垂線上測點的流速,得到垂線平均流速并計算該部分流量,進而累加得到整個斷面的流量。這種測驗方法雖能在精度上滿足流量測驗規(guī)范的要求,但操作復雜、歷時長、工作強度大,效率低下,已跟不上我國水文站點快速增加和功能擴展的需求。隨著社會、經(jīng)濟和科技的快速發(fā)展,水文監(jiān)測現(xiàn)代化的重要性日漸凸顯,其中河流流量監(jiān)測自動化是當今社會發(fā)展對水文工作的必然要求之一,而水文測驗方式和方法的自動化又是整個水文信息化的前提和保障。為此,改進傳統(tǒng)流量測驗方式,廣泛推進流量在線監(jiān)測,提高水文自動化程度和效率,具有十分重要的意義。
流量在線監(jiān)測是指結(jié)合水位及大斷面數(shù)據(jù),根據(jù)采集到的實時流速信息,計算斷面實時流量,并將這些數(shù)據(jù)保存到數(shù)據(jù)庫中的過程。在線監(jiān)測系統(tǒng)主要由水位和流速數(shù)據(jù)的自動采集、數(shù)據(jù)傳輸和分析處理3部分組成,涉及現(xiàn)代傳感技術(shù)、自動測量技術(shù)、自動控制技術(shù)、計算機應用技術(shù)以及相關(guān)的專用分析軟件和通信網(wǎng)絡(luò)[3-4]。其中,現(xiàn)代信息技術(shù)的蓬勃發(fā)展推動了整個水文測驗行業(yè)的進步,新測流技術(shù)的出現(xiàn)不僅提高了測流的時效性,采集到的數(shù)據(jù)量也大大增加。與此同時,計算機處理技術(shù)的發(fā)展讓我們可以充分地挖掘數(shù)據(jù)中的信息,對斷面水流特性有更好的把握。在線監(jiān)測技術(shù)雖經(jīng)幾十年的發(fā)展,逐漸走向生產(chǎn)應用,但其核心問題還未完全解決,如何通過局部流速計算斷面平均流速依然是個難點。
在線流量監(jiān)測的不斷推進將極大提高測流效率,提升我國水文測驗行業(yè)的能力和水平。水利部水文司資料顯示,2018年全國水文系統(tǒng)已裝備在線測流系統(tǒng)1 616套、聲學多普勒流速儀(acoustic doppler current profiler, ADCP)2 665臺,先進儀器設(shè)備較2016年顯著提升。然而,目前我國流量在線監(jiān)測設(shè)備主要還依賴進口,由于在線監(jiān)測技術(shù)本身成熟度不高,加之測站工作人員水平限制、專業(yè)運維團隊匱乏,在線監(jiān)測設(shè)備及系統(tǒng)在許多水文站點應用效果并不好[5-9],流量的精度、測流設(shè)備的穩(wěn)定性及應用的問題依然突出。本文在總結(jié)現(xiàn)有幾種在線監(jiān)測方法優(yōu)點及局限性的基礎(chǔ)上,依據(jù)流量在線監(jiān)測中流量計算的核心技術(shù),對在線監(jiān)測方法中的流速面積法進行分類,并針對精度、穩(wěn)定性和應用等問題對流速面積法及水力學法進行深入探討。
從測驗原理來看,流量在線監(jiān)測可通過流速面積法和水力學法來實現(xiàn)。流速面積法是目前應用最廣泛的方法,該方法通過實測斷面流速與斷面面積來推求流量;在流速面積法中,電磁流速儀、多普勒流速儀等儀器被廣泛應用于流量在線監(jiān)測工作中。水力學法利用水位流量關(guān)系,由測得的水位通過水力學公式計算得到流量,分為量水建筑物法和水工建筑物法[10];水力學方法多用于較小的河流,大江大河流量在線監(jiān)測中應用較少。
現(xiàn)階段基于流速面積法的流量在線監(jiān)測方法還沒有一個明確的分類,通常采用測流儀器進行命名,例如水平式ADCP(H-ADCP)法、坐底式ADCP(V-ADCP)法、雷達法等。這種分類方法雖能直觀地顯示各儀器測流方式之間的差異,但卻無法反映流量計算方法及適用性之間的差異,不利于在線監(jiān)測的推廣應用?;谝陨显颍疚囊罁?jù)斷面平均流速計算原理提出分類方法,將現(xiàn)有在線監(jiān)測方法分為3類:指標流速法、斷面流速分布模型法和表面流速法。該分類方法能夠突出在線流量監(jiān)測的原理和關(guān)鍵技術(shù),有利于在線流量監(jiān)測技術(shù)的完善和推廣應用。
1.1.1指標流速法
目前,國內(nèi)外使用最多的在線監(jiān)測方法為指標流速法[11],該方法專指要通過建立指標流速與斷面平均流速相關(guān)關(guān)系計算斷面流量的方法。這里的指標流速通常指難以明確具體對象或物理含義的實測局部流速,例如H-ADCP采集的斷面某一層的局部流速,該流速一般為某一深度下H-ADCP量程范圍內(nèi)一個錐形水體對象的平均流速。
現(xiàn)階段,結(jié)合H-ADCP開發(fā)流量在線監(jiān)測系統(tǒng)的研究較多,越來越多測流方案的提出,也使解決大型感潮河段以及大江大河測流難點問題成為可能[12-16]。通常情況下,指標流速與斷面平均流速直接構(gòu)建線性相關(guān)關(guān)系[15]。研究發(fā)現(xiàn),受水位變化的影響,不同水位級下建立的回歸方程通常有所差異。為解決此問題,也有將水位作為變量考慮,構(gòu)建基于指標流速、水位的回歸方程;杜耀東等[17]則建議將儀器安放在周期性轉(zhuǎn)動的平臺上,使測量范圍由線變成面獲取指標流速。對于特別復雜的河流,也可同時采用多個指標流速,例如南京水文試驗站采用H-ADCP與定點ADCP結(jié)合的方式獲取指標流速,構(gòu)建多變量相關(guān)方程實現(xiàn)在線監(jiān)測[18-19]??傮w而言,國內(nèi)對H-ADCP的應用研究多停留在單站試驗,尚未有系統(tǒng)的測驗誤差控制及參數(shù)優(yōu)化方法。
1.1.2斷面流速分布模型法
斷面流速分布模型法是指通過采集到的水層或垂線流速,結(jié)合流速分布模型,進行流量實時在線監(jiān)測的工作。該方法又可分為水平流速法和垂線流速法,前者采用水層平均流速通過逐層流速積分或流速模型計算斷面流量,后者則直接采用垂線流速模型計算斷面流量。
a. 水平流速法。該方法通過獲取測流斷面上的水層流速,采用流速分布模型計算得到斷面平均流速,進而得到流量。目前水層流速的獲取方法主要采用超聲波時差法,該法利用超聲波在河段上、下游兩個固定點之間順水和逆水傳播的時間差,計算得到水層平均流速[1,20],該法有單層測流法和多層測流法兩類[10]。在計算流量時,若斷面規(guī)則,可直接通過逐層流速積分計算;反之,則要通過斷面水深-面積模型,根據(jù)水層流速和水深自動計算。水平流速分布模型法適用于較寬的河流,能夠?qū)崿F(xiàn)實時流量監(jiān)測,適合無人值守的測站[21],但該方法土建成本高,運行維護不便,水草較多對其精度也有一定的影響。日本、英國、美國、德國、加拿大、瑞士、法國、韓國、荷蘭等國均采用該法測流[8,10],國內(nèi)對其應用也有研究,但都集中在人工渠道上,天然河道上的成功應用鮮有報道[22-24]。
b. 垂線流速法。該方法主要通過實測一條或多條垂線平均流速,通過流速模型反演其他垂線流速,或建立與斷面平均流速相關(guān)關(guān)系,結(jié)合斷面面積,計算斷面流量。盡管這類方法中也有與指標流速法類似構(gòu)建的相關(guān)關(guān)系,但用于建立相關(guān)的變量有著本質(zhì)的區(qū)別,即垂線平均流速具有明確的物理意義。垂線流速法通常配合V-ADCP使用[25],可有效解決H-ADCP法流量在線監(jiān)測時存在的相關(guān)線建立時間周期長、受航運影響大、低水位時探頭外露和高水外延等問題[26-27]。典型的垂線流速法有兩種:代表線法和二線能坡法。代表線法通過歷史資料,選定至少一條與斷面平均流速相關(guān)性較好的測速垂線作為代表線,并建立相關(guān)關(guān)系,結(jié)合斷面面積,計算斷面流量。國內(nèi)對該方法有較多應用,例如徐六涇水文站用該方法解決潮流量測驗問題[28]、前垾村水文站利用該方法為尋找典型點流速提供支撐[29]、龍門水文站利用該方法降低高洪作業(yè)風險[30]。但國內(nèi)大部分應用還存在原始資料系列較短、斷面穩(wěn)定性及流量精度評價指標不完整、流速相關(guān)關(guān)系使用條件不明確等問題。二線能坡法以曼寧公式為基礎(chǔ),通過虛擬矩形、三角形斷面中垂線平均流速與斷面平均流速之間關(guān)系,建立垂線流速模型。實際計算中,以實測的兩條垂線流速為已知條件,代入垂線流速模型,反求得到能坡,作為均勻流條件下曼寧公式中的水面比降。最后通過垂線流速計算模型,結(jié)合水位面積關(guān)系計算斷面流量,國內(nèi)已有少量應用成功的例子[31-32]。
1.1.3表面流速法
表面流速法指通過獲取測流斷面上水面流速,并通過轉(zhuǎn)換得到特征垂線平均流速或斷面平均流速進行流量計算的方法。表面流速的獲取可以采用纜道雷達、側(cè)掃雷達、粒子圖像法等非接觸式測流儀器。非接觸式儀器不受污水及泥沙影響、不影響水流狀態(tài),尤其在陡漲險落的山區(qū)性河流的監(jiān)測中具有明顯優(yōu)勢[33-35]?,F(xiàn)階段國內(nèi)應用較廣的表面流速法測流有雷達法和粒子圖像法。
a. 雷達法。該法又稱電波法,是基于微波多普勒效應測量水面流速的一種方法,主要通過雷達波傳感器自動發(fā)射和接收電磁波,利用多普勒原理測定水面流速,通過流速系數(shù)或回歸分析得到斷面平均流速,從而得到斷面流量。雷達法測流不受含沙量及漂浮物的影響,可用于巡測,其測驗精度高于浮標法,可以較好地應用于在線監(jiān)測工作中。但雷達法對測量條件敏感,易受大風大雨及測量角度影響,低流速測量誤差大[33]。非接觸式雷達測速儀的應用可有效提高測流效率,但目前大部分測驗均存在比測時間短、測次少、比測結(jié)果不具有代表性等缺陷,而且這些站多是中小河流水文站,在大江大河上適用性還有待驗證[36-38]。國產(chǎn)側(cè)掃雷達的研究也取得一定成果,并已在黃河蘭州站、花園口站、長江南京站、廣西南寧站等站點得到良好應用[39]。該設(shè)備環(huán)境適應性強,受氣象條件影響很小,適用于惡劣天氣、高洪多漂浮物、漫灘等惡劣環(huán)境。國產(chǎn)側(cè)掃雷達性能優(yōu)良,表面流速測驗精度較高,但由于自然河道流體運動的復雜性,絕大部分情況相關(guān)性較差,流量精度還難以滿足測驗要求。
b. 粒子圖像法。粒子圖像法利用拍攝到的河流水面視頻,估計示蹤物在圖像序列中運動矢量的大小,分析測流區(qū)域內(nèi)局部流體的運動規(guī)律,通過流速系數(shù)將水面流速轉(zhuǎn)換為斷面平均流速,從而得到斷面流量[40]。該方法在快速獲取水面流場及紊動特性、高洪及極端條件流量監(jiān)測有明顯優(yōu)勢。目前,日本、韓國等國家正在大量開展該方法的試驗研究工作[35]。粒子圖像法又包括極坐標攝影浮標法、衛(wèi)星遙感圖像法、大尺度粒子圖像測速法和時空圖像法。極坐標攝影浮標法采用攝影照片代替人工目測,雖然精度和自動化程度較高,但由于部署工作量較大,實際應用較少。衛(wèi)星遙感圖像法對地面信息和歷史數(shù)據(jù)的依賴過大,同時測量誤差較大,還無法實用化。大尺度粒子圖像測速法具備時空分辨率高、測量范圍廣、成本低廉等優(yōu)點[8],不僅可用于常規(guī)條件下測流工作,更具有極端條件下河道水流監(jiān)測的應用潛力[33,41],韓國[8]以及我國吉安市坳下坪水文站[42]已開展相關(guān)研究。時空圖像法是一種高空間分辨率的一維時均運動矢量估計方法,目前也處于試驗階段[43-45]。
水力學法不受人為因素及河流特性的影響,無須建設(shè)測流纜道,容易實現(xiàn)遙測,是一類重要的方法并且具有推廣價值。該方法可進一步分為量水建筑物法和水工建筑物法。量水建筑物法通常采用量水堰、量水槽、量水池等量水建筑物測流;英國有一半水文站采用量水建筑物法,該方法對建筑物的建設(shè)要求比較嚴格,適用于小型河流,需在實驗室做模擬試驗計算流量系數(shù),或者在典型水文站上校測流量系數(shù)[46]。水工建筑物法采用堰、水電站、泵站等水工建筑物測流[4],該法根據(jù)能量轉(zhuǎn)換和守恒原理,通過實測水頭差、閘門開啟高度等水力因素,確定流量系數(shù)或效率系數(shù),用水力學公式計算流量。水工建筑物法精度較高,程序簡便,但同時存在泄水建筑物造價昂貴的缺點,韓國大量利用該方法開展流量測驗和在線監(jiān)測工作[9],我國三河閘[47]、柳河[48]、琚灣[49]、飛來峽[50]等水文站均有相關(guān)研究;近年來,許多省份有在適宜的河流上推廣該方法的意向,在引黃灌區(qū)、多地供水設(shè)施還采用不同形式的堰閘在線測流設(shè)施設(shè)備[51]。
沒有一種設(shè)備是可以滿足所有斷面和水流情況的,具體應用時應結(jié)合測站任務及河道自然屬性,選擇合適的測流方法,從而確定所需的設(shè)備儀器。流量測驗首要保證測流的精度,但現(xiàn)階段在線測流技術(shù)在該方面還存在不足。例如,中泓最大水面流速低于0.2 m/s時,暫時還沒有一個精準的測流方法;如何進行點流速比測也是一個難點問題;側(cè)掃雷達使用時可能受回水及船只影響,造成流速測量誤差較大,還需不斷優(yōu)化流量計算模型;同時,當?shù)厮臈l件、測站工作人員水平等問題也是影響測驗精度的重要因素。
無論采用哪種測流方法,如何通過實測的局部流速獲得斷面平均流速一直是個難點。要提高流量測驗精度,關(guān)鍵還是要根據(jù)測驗河段的水力學特性,充分應用現(xiàn)代化測流手段和數(shù)值計算技術(shù),確定特定測站的流量計算方法和模型,率定其所需要的參數(shù)。隨著聲學多普勒測流技術(shù)及三維水動力模型的逐漸成熟,可將固定式ADCP與水力學模型結(jié)合,建立三維流速模型來推求斷面上每一處水流影響因子和程度[52-54]。結(jié)合ADCP采集到的流速,進行流量在線監(jiān)測,這將是解決在線監(jiān)測精度問題的突破方向。對于非接觸測流技術(shù)尤其需要注意風速的影響,未來還需建立風場與水面流場的相關(guān)關(guān)系,以提高設(shè)備精度。
在線監(jiān)測儀器的運行穩(wěn)定性受流量、流速、風速等外界條件影響。雷達法、二線能坡法等在水量較小時測量誤差較大,并且水量較小也會加大風速的影響。側(cè)掃雷達還可能受到周邊相同頻率的干擾源的影響,大大增加測流誤差。聲學多普勒測流技術(shù)已相當成熟,但H-ADCP只能測得儀器安裝處水層的局部流速分布,而且是近岸一段,而水位不斷變化勢必會造成水層流速代表性的變化;若多次移動儀器,不僅增加了測流工作量,還可能影響測流精度。因此,可基于歷史水文資料構(gòu)建水位與H-ADCP安裝位置及傾角之間的相關(guān)關(guān)系,采用可自動調(diào)整傾角的底座,大大減輕工作量。大江大河的寬度遠大于H-ADCP的測流范圍,除采用組合式監(jiān)測系統(tǒng)外,還可將采集到的數(shù)據(jù)與歷史水文數(shù)據(jù)結(jié)合,補全斷面流速,再計算得到斷面流量[55]。
在線監(jiān)測技術(shù)已經(jīng)歷幾十年的發(fā)展,但缺乏專業(yè)技術(shù)人員進行比測率定、定期維護等工作,造成測驗設(shè)備在應用過程中問題層出[9,33,56]。同時,感潮河段水文現(xiàn)象復雜,河口區(qū)河面寬闊、航運繁忙,國內(nèi)雖有相關(guān)研究,但河口流量測驗還需尋找具有快速全場測量的能力并能適應復雜的測流環(huán)境的測流方案。我國非接觸測流設(shè)備主要依賴于進口,許多國產(chǎn)ADCP性價比有待提高,高端儀器產(chǎn)品基本為空白[7],缺少更加系統(tǒng)全面的技術(shù)依據(jù)及相關(guān)技術(shù)指標,這些問題限制了該設(shè)備的廣泛應用。側(cè)掃雷達、國產(chǎn)ADCP等技術(shù)還未成熟,配置技術(shù)參數(shù)還不夠明晰,比測時均存在水位變幅較小的問題。
衛(wèi)星遙感圖像法以及低空遙感無人機測流技術(shù)[57]雖有相關(guān)研究,但測流精度依然較低,僅在試驗河段及站點測流結(jié)果較好,無法滿足實際測流需要,還需一段時間的研究發(fā)展,短期內(nèi)無法為在線監(jiān)測服務,但這些方法對于洪水、堰塞湖、泥石流等條件下的應急監(jiān)測具有重要意義。
在自動采集水文信息后通過現(xiàn)代通訊技術(shù)進行數(shù)據(jù)傳輸,利用計算機自動接受、處理、存儲分發(fā)水文信息,是水文信息化的發(fā)展方向[58]。為了加快水文信息化的發(fā)展,要走在線監(jiān)測、遠程控制、巡駐結(jié)合的道路,應用多源傳感器信息融合技術(shù)[28],加快物聯(lián)網(wǎng)和5G數(shù)據(jù)傳輸技術(shù)的研究,加快水文儀器國產(chǎn)化進程,真正實現(xiàn)“有人看管,無人值守”,實現(xiàn)互聯(lián)網(wǎng)+水文深度融合。
本文總結(jié)了目前大江大河流量在線監(jiān)測中廣泛采用的幾種方法,依據(jù)流速數(shù)據(jù)采集及流量計算原理對監(jiān)測方法進行了分類。目前在線監(jiān)測推廣應用工作正在如火如荼地展開,但缺乏專業(yè)的技術(shù)團隊對監(jiān)測儀器誤差控制、精度分析等方面進行深入研究,這種失衡是導致在線監(jiān)測推廣應用停滯不前的一大原因。本文針對在線監(jiān)測精度、穩(wěn)定性和應用方面的問題進行詳細論述,并指出固定式ADCP與水力學模型結(jié)合將是解決我國在線監(jiān)測問題的突破方向;同時,在具備條件的地區(qū),使用測流建筑物的水力學法測流也是一條可行的途徑。