祁曉勇,高玉方
(青海省水文水資源測(cè)報(bào)中心,西寧 810000)
水文觀測(cè)自動(dòng)化是中國(guó)水文現(xiàn)代化發(fā)展的必然要求[1],Bettersize2000W 型激光粒度分布儀等新技術(shù)、新裝備的應(yīng)用為實(shí)現(xiàn)水文現(xiàn)代化奠定了堅(jiān)實(shí)的基礎(chǔ),有助于擴(kuò)大水文資料收集范圍[2],改進(jìn)水文工作向廣度和深度發(fā)展。水體含沙量監(jiān)測(cè)是水文要素觀測(cè)中十分重要的項(xiàng)目[3],含沙量一般是指單位體積的渾水中所含干沙的質(zhì)量,通過(guò)系統(tǒng)地觀測(cè)和收集含沙量資料,研究分析河流內(nèi)泥沙的分布規(guī)律[4]。觀測(cè)站對(duì)單樣含沙量的觀測(cè)主要采用稱重法,該方法簡(jiǎn)單、精度較高[5],但施測(cè)過(guò)程耗時(shí)較長(zhǎng),無(wú)法在短時(shí)間內(nèi)得到水體含沙量數(shù)據(jù)[6]。隨著剛察水文站水文測(cè)驗(yàn)方式、方法的不斷改進(jìn),Bettersize 2000W 型激光粒度分布儀等新儀器、新設(shè)備的推廣應(yīng)用已成為重要課題[7,8],在保證測(cè)驗(yàn)質(zhì)量的基礎(chǔ)上,應(yīng)用先進(jìn)的觀測(cè)設(shè)施、設(shè)備不僅能降低水文測(cè)驗(yàn)的工作強(qiáng)度[9],提高工作效率,而且能節(jié)約一定的人力、財(cái)力,為全面開(kāi)展水文監(jiān)測(cè)改革奠定基礎(chǔ)[10]。
剛察(二)水文站位于青海省剛察縣沙柳河鎮(zhèn),距西寧約196 km,位于東經(jīng)100°08′,北緯37°19′。剛察水文站的測(cè)驗(yàn)河段位于彎道處,為窄深型河槽,左岸陡崖,右岸均為漿砌石護(hù)岸,高約2 m,右岸堤防背側(cè)灘地長(zhǎng)滿雜草及灌木。主流靠右岸,河床由沙礫石組成,測(cè)驗(yàn)斷面有一定的沖淤變化,屬于經(jīng)常性沖淤、整體沖淤,斷面不穩(wěn)定。洪水主要來(lái)源于融雪及降水,以暴雨洪水為主,洪峰過(guò)程多為單峰,暴漲暴落,主峰雨停時(shí)間至洪峰出現(xiàn)時(shí)間約為5 h,洪峰持續(xù)0.5~2.0 h,水沙峰基本相應(yīng),輸沙量年內(nèi)分配不均,年際變化也較大,較大洪水過(guò)程一般持續(xù)1~2 d。
Bettersize2000W 型激光粒度分布儀是由丹東百特儀器有限公司設(shè)計(jì)制造的雙鏡頭激光粒度儀,該系統(tǒng)包括激光粒度儀、自動(dòng)循環(huán)分散系統(tǒng)以及電腦軟件系統(tǒng)等,通過(guò)自動(dòng)循環(huán)分散系統(tǒng)將樣品輸送到激光粒度儀的測(cè)量區(qū)域,激光照射到樣品后將產(chǎn)生散射光,光電探測(cè)器陣列接收這些散射光信號(hào)并轉(zhuǎn)換成電信號(hào),通過(guò)USB 接口將這些電信號(hào)傳輸?shù)诫娔X軟件中,依據(jù)Mie 散射理論對(duì)散射光產(chǎn)生的電信號(hào)進(jìn)行處理,就可以得到所測(cè)樣品的粒度分布及含沙量等結(jié)果。該套自動(dòng)化監(jiān)測(cè)設(shè)備功能齊全、自動(dòng)化程度高,監(jiān)測(cè)數(shù)據(jù)可實(shí)現(xiàn)實(shí)時(shí)測(cè)量及本地存儲(chǔ)。Bettersize2000W 型激光粒度分布儀見(jiàn)圖1。
圖1 Bettersize2000W 型激光粒度分布儀
1)系統(tǒng)誤差。系統(tǒng)誤差是在對(duì)同一被測(cè)量進(jìn)行多次測(cè)量過(guò)程中,出現(xiàn)某種保持恒定或按確定的方法變化的誤差。
式中,P為系統(tǒng)誤差;Ai為某個(gè)量的觀測(cè)值(采用激光粒度分布儀的觀測(cè)值);A為某個(gè)量的真值(采用傳統(tǒng)儀器或方法的觀測(cè)值);n為觀測(cè)次數(shù)。
2)絕對(duì)誤差。絕對(duì)誤差也稱真誤差,是指觀測(cè)值與真值之差,即相對(duì)于真值的誤差稱為絕對(duì)誤差。
式中,Δ 為絕對(duì)誤差。
3)相對(duì)誤差。相對(duì)誤差定義為觀測(cè)值的絕對(duì)誤差與其真值之比。
式中,δ為相對(duì)誤差。
4)標(biāo)準(zhǔn)差。標(biāo)準(zhǔn)差也稱均方差,是各數(shù)據(jù)偏離平均數(shù)的距離的平均數(shù),反映一個(gè)數(shù)據(jù)集的離散程度。
式中,s為標(biāo)準(zhǔn)差,%;Xi為第i個(gè)觀測(cè)數(shù)據(jù);為N次觀測(cè)值的平均值;N為觀測(cè)次數(shù)。
5)合格率。合格率即對(duì)比觀測(cè)時(shí)段內(nèi)自動(dòng)監(jiān)測(cè)數(shù)據(jù)比測(cè)合格的百分比,計(jì)算式如下。
分析資料為2021 年5—9 月收集烘干法和激光粒度分布儀法測(cè)得的115 組含沙量,所有數(shù)據(jù)均經(jīng)過(guò)可靠性、一致性、有效性分析,剔除偶然誤差和人為操作誤差影響的無(wú)效數(shù)據(jù),確保所參與分析的數(shù)據(jù)真實(shí)、準(zhǔn)確、有效。
根據(jù)2021 年5—9 月收集的115 組含沙量對(duì)比數(shù)據(jù)進(jìn)行誤差統(tǒng)計(jì)分析,發(fā)現(xiàn)激光粒度儀所測(cè)數(shù)據(jù)普遍偏大,偏大的數(shù)據(jù)占比為95.7%,絕對(duì)值最小的相對(duì)誤差為5.2%,絕對(duì)值最大的相對(duì)誤差為952.2%,平均相對(duì)誤差為255.8%,相對(duì)誤差絕對(duì)值小于10%(1 個(gè))的占0.87%,相對(duì)誤差絕對(duì)值小于50%(14 個(gè))的占12.17%,相對(duì)誤差絕對(duì)值大于50%(101 個(gè))的數(shù)據(jù)占比為87.83%;絕對(duì)誤差在-0.212~3.412,最小絕對(duì)誤差為-0.212,最大絕對(duì)誤差為3.412,平均絕對(duì)誤差為0.277。原始數(shù)據(jù)及數(shù)據(jù)統(tǒng)計(jì)分析結(jié)果見(jiàn)表1。
表1 激光粒度分布儀與稱重法實(shí)測(cè)含沙量對(duì)比數(shù)據(jù)統(tǒng)計(jì)分析結(jié)果
為探尋二者的相關(guān)性,建立線性方程,以檢驗(yàn)是否滿足測(cè)驗(yàn)精度的要求。對(duì)激光粒度分布儀法測(cè)得的含水量(y)與烘干法測(cè)得的含沙量(x)建立線性相關(guān)關(guān)系,如圖2 所示,發(fā)現(xiàn)二者可以建立一定的線性相關(guān)關(guān)系,并通過(guò)趨勢(shì)線確定最優(yōu)相關(guān)關(guān)系為y=3.014x,由圖2 可見(jiàn),二者線性相關(guān)性一般。
圖2 線性關(guān)系對(duì)比分析
1)根據(jù)分析結(jié)果,激光粒度分布儀測(cè)試含沙量的影響因素主要有以下幾個(gè)方面。一是采樣及操作人員技術(shù)水平不一造成的人為操作誤差。二是各種稱量器具對(duì)水樣體積的稱量誤差。三是若測(cè)試水溫過(guò)低,在激光測(cè)試窗口會(huì)形成凝結(jié)霧狀薄膜對(duì)激光測(cè)試結(jié)果產(chǎn)生影響。
2)激光粒度分布儀測(cè)試分析的過(guò)程中應(yīng)注意以下幾點(diǎn)。一是2 種方法測(cè)含沙量所取的水樣應(yīng)在同一位置及同一時(shí)間進(jìn)行取樣,這樣能排除取樣誤差的影響。二是加樣品量要精確用量杯稱量,并要將量杯中的樣品全部倒入循環(huán)池中,并用清水涮底。三是取樣中的雜質(zhì)問(wèn)題,取樣過(guò)程中,現(xiàn)取現(xiàn)測(cè),防止大顆粒沉淀,出現(xiàn)波動(dòng)較大情況,重新取樣。
從測(cè)試分析過(guò)程看,激光粒度分布儀測(cè)含沙量操作簡(jiǎn)單方便,效率高,含沙量結(jié)果計(jì)算速度快,而人工烘干法測(cè)得的含沙量數(shù)值偏差較大,對(duì)激光粒度分布儀與烘干法測(cè)得含沙量建立線性相關(guān)關(guān)系,二者可以建立一定的線性相關(guān)關(guān)系,并通過(guò)趨勢(shì)線確定最優(yōu)相關(guān)關(guān)系為y=3.014x,此相關(guān)系數(shù)所得成果并不能作為正式資料使用,但在水文應(yīng)急監(jiān)測(cè)等精度要求不是很高的情況下可以參考使用,未來(lái)使用激光粒度分布儀測(cè)試含沙量作為正式資料還需要進(jìn)一步研究和論證。