王先賀，陳華民，童仁園，葛仁望，陳向陽，李 青*

(1.中國計(jì)量大學(xué)災(zāi)害監(jiān)測技術(shù)與儀器國家地方聯(lián)合工程實(shí)驗(yàn)室，浙江 杭州310018；2.金華市地質(zhì)環(huán)境監(jiān)測站，浙江 金華321000)

我國是一個(gè)多山的國家，是地質(zhì)災(zāi)害頻發(fā)的國家之一。 由于全球氣候異常變化，我國局部地區(qū)短時(shí)間內(nèi)出現(xiàn)暴雨、特大暴雨的情況日漸增多，這使得發(fā)生滑坡和泥石流等地質(zhì)災(zāi)害的隱患也在不斷增加[1]。 短時(shí)間內(nèi)持續(xù)性的強(qiáng)降雨是引發(fā)的滑坡和泥石流等地質(zhì)災(zāi)害的關(guān)鍵因素，因此，如何實(shí)時(shí)、準(zhǔn)確地獲取降雨數(shù)據(jù)，及時(shí)采取防范措施，降低不必要的損失，對于地質(zhì)災(zāi)害的預(yù)警工作顯得尤為重要。

雨量計(jì)是氣象及社會各相關(guān)部門普遍應(yīng)用的測量降雨的儀器，是指能記錄和顯示降雨量的觀測儀器。 目前國內(nèi)應(yīng)用較多的雨量計(jì)有翻斗式、虹吸式、稱重式、壓電式等多種類型[2-3]。 但每種類型的雨量計(jì)在精度或操作上都有一定的局限性。 翻斗式雨量計(jì)受降雨強(qiáng)度大小的限制，雨越大，測量誤差越大[4]；虹吸式雨量計(jì)虹吸過程會漏記降雨量，而且必須定時(shí)到現(xiàn)場去更換記錄紙，操作繁瑣[5-6]；稱重式雨量計(jì)不能自動(dòng)排水或排水不及時(shí)，造成降雨測量數(shù)據(jù)不準(zhǔn)確[7]；壓電式雨量計(jì)受環(huán)境影響大，風(fēng)、灰塵等易引起雨量板振動(dòng)頻率變化，導(dǎo)致數(shù)據(jù)誤差變大[8-9]。 國外上世紀(jì)50年代就開始了降雨自動(dòng)計(jì)量的研究，起步較早，目前已經(jīng)實(shí)現(xiàn)了降雨的連續(xù)、非接觸自動(dòng)測量，國外比較成熟的產(chǎn)品有:芬蘭VAISALA 公司生產(chǎn)的PWD22 天氣現(xiàn)象傳感器[10]，德國OTT 公司的Parsivel 激光降水粒子譜測量系統(tǒng)[11]，美國Borrmann 和Jaenicke 研制的全息液滴和氣溶膠系統(tǒng)。 以上測量系統(tǒng)也有測量不準(zhǔn)確，易受環(huán)境影響等缺點(diǎn)[12]。 國內(nèi)對于雨量計(jì)的研制工作起步較晚，存在一定的差距，因此，研究一款可以精確測量降雨量的裝置意義重大。

針對以上問題，在傳統(tǒng)虹吸式雨量計(jì)的基礎(chǔ)上進(jìn)行了改進(jìn)，設(shè)計(jì)了一種虹吸排水自補(bǔ)償降雨量測量裝置。 提出了該裝置的系統(tǒng)結(jié)構(gòu)及原理技術(shù)，設(shè)計(jì)了硬件電路，實(shí)現(xiàn)了自動(dòng)化監(jiān)測，建立了一套完整的降雨監(jiān)測系統(tǒng)。 針對虹吸排水帶來的計(jì)量誤差問題，建立了一種高精度虹吸補(bǔ)償數(shù)學(xué)模型，有效地提高了測量精度和靈敏度，在降雨監(jiān)測領(lǐng)域具有很好的應(yīng)用前景。

1 系統(tǒng)結(jié)構(gòu)及測量原理

為了解決傳統(tǒng)虹吸式雨量計(jì)虹吸時(shí)產(chǎn)生計(jì)量誤差的問題，實(shí)現(xiàn)對降雨的高精度監(jiān)測，設(shè)計(jì)了虹吸排水自補(bǔ)償降雨量測量裝置，結(jié)構(gòu)示意圖及實(shí)物圖如圖1、圖2 所示。 其工作流程可表述如下，當(dāng)有降雨發(fā)生時(shí)，降雨順著承水漏斗流入盛水器，盛水器中的水位便隨之上升。 當(dāng)盛水器中的水位高度達(dá)到A點(diǎn)時(shí)，盛水器中的雨水便會由于虹吸作用，順著虹吸管自動(dòng)排出到裝置外，直到盛水器中的水位降到O點(diǎn)，虹吸排水才算結(jié)束。 我們將盛水器中的水位從O點(diǎn)開始，上升到A點(diǎn)后，發(fā)生虹吸排水再回落到O點(diǎn)的過程，稱為一次降雨虹吸過程。 當(dāng)發(fā)生虹吸排水時(shí)，測量裝置處于一邊接水一邊放水的狀態(tài)，那么虹吸排水的這段時(shí)間里的降雨也會被隨之排走而漏記，所以需要對虹吸過程中漏記的降雨進(jìn)行補(bǔ)償，本文提供了一種高精度虹吸補(bǔ)償數(shù)學(xué)模型對虹吸排水過程造成的計(jì)量誤差進(jìn)行補(bǔ)償。 虹吸補(bǔ)償算法將在下面章節(jié)詳細(xì)介紹。

圖1 裝置結(jié)構(gòu)示意圖

圖2 裝置內(nèi)部結(jié)構(gòu)實(shí)物圖

本文的測量原理如下，當(dāng)降雨發(fā)生時(shí)，測量裝置收集降雨，位于盛水器底部的稱重傳感器的輸出電壓會隨雨水重量的不同而產(chǎn)生相應(yīng)變化。 然后將稱重傳感器輸出的電信號傳送給數(shù)據(jù)采集系統(tǒng)，由數(shù)據(jù)采集系統(tǒng)中的內(nèi)部電路完成相應(yīng)的差分放大、模數(shù)轉(zhuǎn)換、濾波、數(shù)據(jù)處理、虹吸補(bǔ)償?shù)忍幚砉ぷ鳎玫骄珳?zhǔn)的降雨重量，再把降雨重量轉(zhuǎn)換成降雨量，通過電路搭載的4G 物聯(lián)網(wǎng)模塊，將降雨數(shù)據(jù)上傳到阿里云平臺，實(shí)現(xiàn)遠(yuǎn)程監(jiān)測。

2 測量系統(tǒng)的研究與設(shè)計(jì)

2.1 虹吸管的設(shè)計(jì)

虹吸管作為雨量計(jì)的排水部件，性能的優(yōu)劣直接影響到整套設(shè)備的設(shè)計(jì)要求，因此選取合適的虹吸管非常重要。 虹吸管的工作原理是利用虹吸現(xiàn)象引起排水的自然過程，無需消耗能量，因此非常符合系統(tǒng)低功耗的設(shè)計(jì)要求。 發(fā)生虹吸排水時(shí)，當(dāng)水面升至彎曲段時(shí)，由于水的重力作用，下浸潤邊超前上浸潤邊，當(dāng)水達(dá)到虹吸管頂時(shí)，下浸潤邊可能產(chǎn)生滴水，帶來測量誤差。 通過多次實(shí)驗(yàn)驗(yàn)證后，最終選用玻璃彎管弧度為28°，管內(nèi)徑為4 mm，厚度為1.5 mm 的虹吸管，虹吸管結(jié)構(gòu)示意圖如圖3 所示。 當(dāng)臨近虹吸時(shí)，該虹吸管不會發(fā)生滴水現(xiàn)象；當(dāng)雨水達(dá)到虹吸高度時(shí)，也不會出現(xiàn)虹吸無法發(fā)生的現(xiàn)象。 根據(jù)國家標(biāo)準(zhǔn)《GB/T 21978.3-2008 降雨量觀測儀器》，連續(xù)降雨強(qiáng)度范圍為0.01 mm/min～4.00 mm/min，通過查閱資料了解到，國內(nèi)最大降雨速度的記錄發(fā)生在1971年7月1 日，山西太原地區(qū)古交市梅洞溝，5 min 降雨量達(dá)53.1 mm，平均每分鐘的降雨量為10.62 mm。 該虹吸管的最大排水速度18.3 mm/min，不會發(fā)生阻流現(xiàn)象，滿足設(shè)計(jì)要求。

圖3 虹吸管結(jié)構(gòu)示意圖

2.2 應(yīng)變式稱重傳感器工作原理

傳感模塊采用的是電阻應(yīng)變式稱重傳感器，其工作原理是彈性體在外力作用下產(chǎn)生彈性變形，使粘貼在它表面的電阻應(yīng)變片也隨同產(chǎn)生變形；電阻應(yīng)變片變形后，它的阻值將發(fā)生變化，再經(jīng)相應(yīng)的測量電路把這一電阻變化轉(zhuǎn)換為電信號輸出。 該類型傳感器精度高、頻響特性好、結(jié)構(gòu)簡單、易實(shí)現(xiàn)小型化、耐勞時(shí)間長、能在惡劣條件下工作，符合應(yīng)用條件。 電阻應(yīng)變式稱重傳感器機(jī)械結(jié)構(gòu)和內(nèi)部電路圖如圖4 所示。

圖4 稱重傳感器機(jī)械結(jié)構(gòu)和內(nèi)部電路圖

該傳感器的工作原理主要是利用了應(yīng)變電阻的壓阻效應(yīng)和惠斯通電橋的平衡原理，惠斯通電橋的輸出電壓為:

式中:E為激勵(lì)電壓，R1～R4為應(yīng)變電阻的阻值，ΔR為應(yīng)變電阻因形變產(chǎn)生的變化量。 有R1=R2=R3=R4=R，此時(shí)惠斯通電橋的輸出電壓為:

惠斯通電橋能感應(yīng)微小的電阻變化，使電橋產(chǎn)生一個(gè)與壓力成正比的高度線性、與激勵(lì)電壓也成正比的電壓信號[13]。

2.3 稱重傳感器的電路設(shè)計(jì)

稱重傳感器電路決定著測量系統(tǒng)的工作效果，稱重傳感器的硬件電路主要包括:信號采集電路、放大濾波電路、模數(shù)轉(zhuǎn)換電路、無線傳輸電路、測溫電路、供電電源電路以及微控制器，具體設(shè)計(jì)框圖如圖5所示。

圖5 硬件電路設(shè)計(jì)示意圖

稱重傳感器會將降雨重量轉(zhuǎn)化為微弱電壓信號輸出，經(jīng)過差分放大電路和濾波電路完成信號放大和噪聲濾除，數(shù)模轉(zhuǎn)換電路將采集到的電信號轉(zhuǎn)化為數(shù)字信號以SPI 串行通信的方式傳輸至微控制器，微控制器對采集到的數(shù)字信號進(jìn)行了軟件濾波，本設(shè)計(jì)使用了“滑動(dòng)平均值濾波法”，并且在“滑動(dòng)平均值濾波法”的基礎(chǔ)上先得到一定數(shù)據(jù)寬度的滑動(dòng)窗口數(shù)據(jù)，然后進(jìn)行“中值濾波”，最后求平均值，從而在一定程度上減小干擾信號的影響。

接下來將濾波后的數(shù)字量和砝碼進(jìn)行標(biāo)定，得到砝碼重量與數(shù)字量的關(guān)系。 本設(shè)計(jì)以0.2 g 為一個(gè)單位進(jìn)行標(biāo)定，從0 g 至滿量程1 kg 依次增大砝碼的重量，并進(jìn)行多次重復(fù)測量，記錄砝碼的重量和數(shù)據(jù)采集系統(tǒng)得到的數(shù)字量，并根據(jù)最小二乘法原理，對測量數(shù)據(jù)進(jìn)行曲線擬合，其關(guān)系曲線如圖6所示，橫坐標(biāo)為數(shù)據(jù)采集系統(tǒng)得到的數(shù)字量，縱坐標(biāo)為砝碼重量。

圖6 稱重標(biāo)定曲線圖

其線性表達(dá)式為:

其線性擬合曲線線性相關(guān)系數(shù)為0.999975，非常接近于1，說明數(shù)字量與砝碼重量具有良好的線性關(guān)系。 標(biāo)定完成后，制作完成的電子秤量程為1 kg，最小分度值為0.5 g，精度為1/2000。 電子秤實(shí)物圖如圖7 所示。 通過制作完成的電子秤測量降雨重量，微控制器進(jìn)行數(shù)據(jù)處理和虹吸補(bǔ)償，得到降雨數(shù)據(jù)，并通過物聯(lián)網(wǎng)模塊將降雨數(shù)據(jù)上傳到阿里云平臺實(shí)現(xiàn)遠(yuǎn)程監(jiān)測。

圖7 電子秤實(shí)物圖

根據(jù)國家標(biāo)準(zhǔn)《GB/T 21978.3-2008 降雨量觀測儀器》要求:虹吸式雨量計(jì)集雨桶直徑為200 mm，分辨率要達(dá)到0.1 mm。 該電子秤的分辨率為0.5 g，已知水的密度ρ 為0.001 g/mm3，設(shè)計(jì)完成的雨量計(jì)集雨桶直徑為200 mm，經(jīng)計(jì)算，該雨量計(jì)的分辨率可達(dá)0.02 mm，符合設(shè)計(jì)要求。

考慮到本裝置需要長期工作在野外進(jìn)行降雨量的監(jiān)測，實(shí)際的工況環(huán)境難以實(shí)現(xiàn)有線連接，因此要想保證裝置能夠長期穩(wěn)定工作，裝置所需能量的外部獲取以及低功耗設(shè)計(jì)是非常重要的環(huán)節(jié)。 鑒于該裝置放置在野外環(huán)境中工作，太陽能可作為一種穩(wěn)定可靠的能量來源，所以本設(shè)計(jì)采用太陽能電池板和蓄電池作為裝置的電源。 考慮到功耗問題，本裝置的硬件電路設(shè)計(jì)均采用單電源供電，降低了功耗，簡化了電路的設(shè)計(jì)結(jié)構(gòu)。 另外，降雨過程是一個(gè)階段性非持續(xù)的過程，因此對降雨的監(jiān)測頻率不需要太高，一個(gè)小時(shí)或者幾個(gè)小時(shí)監(jiān)測一次即可。 為了進(jìn)一步降低功耗，設(shè)計(jì)了休眠/喚醒模式，通過STM32 單片機(jī)進(jìn)行控制，使得裝置可以根據(jù)實(shí)際需求工作在休眠和喚醒模式下。 如果沒有發(fā)生降雨，系統(tǒng)便會工作在休眠狀態(tài)，由于遠(yuǎn)程傳輸模塊功耗相對較大，在休眠模式下，單片機(jī)會切斷4G 物聯(lián)網(wǎng)模塊的供電網(wǎng)絡(luò)，僅保留STM32 單片機(jī)、稱重傳感器、儀表運(yùn)放AD8422 和模數(shù)轉(zhuǎn)換芯片AD7606 的供電，休眠模式下，測量裝置不會進(jìn)行監(jiān)測，這樣會大大降低功耗。 當(dāng)數(shù)據(jù)采集系統(tǒng)檢測到降雨發(fā)生時(shí)，單片機(jī)控制系統(tǒng)會打開4G 物聯(lián)網(wǎng)模塊的供電電源，裝置會自動(dòng)對降雨進(jìn)行監(jiān)測。 以上工作狀態(tài)都是由單片機(jī)控制系統(tǒng)自動(dòng)完成的，經(jīng)過試驗(yàn)測試，在太陽能作為直接能量來源的條件下，配合蓄電池在夜晚或者陰天的情況下提供電源，本裝置能夠?qū)崿F(xiàn)長期穩(wěn)定工作。

2.4 系統(tǒng)軟件設(shè)計(jì)

系統(tǒng)軟件設(shè)計(jì)流程如圖8 所示。

圖8 系統(tǒng)軟件設(shè)計(jì)流程圖

系統(tǒng)上電后，初始化各個(gè)模塊，全部初始化完畢后就進(jìn)入待機(jī)狀態(tài)，不停地檢測是否有數(shù)據(jù)傳輸過來。 當(dāng)系統(tǒng)檢測到有降雨發(fā)生時(shí)，便喚醒系統(tǒng)，處理數(shù)據(jù)，并把降雨數(shù)據(jù)發(fā)送至云端。 監(jiān)測系統(tǒng)可以定時(shí)發(fā)送降雨數(shù)據(jù)，本設(shè)計(jì)每隔1 h 發(fā)送一次數(shù)據(jù)，如果前1 h 的降雨超過7 mm，監(jiān)測端自動(dòng)發(fā)出大雨報(bào)警。 如果監(jiān)測系統(tǒng)檢測到?jīng)]有降雨或降雨停止，系統(tǒng)便進(jìn)入待機(jī)休眠狀態(tài)。

3 虹吸補(bǔ)償算法數(shù)學(xué)模型

當(dāng)發(fā)生降雨時(shí)，盛水器內(nèi)的降雨重量逐漸增大，當(dāng)盛水器內(nèi)的水位高度達(dá)到虹吸點(diǎn)高度時(shí)(此時(shí)在虹吸點(diǎn)測得的降雨重量為mN)，由于虹吸作用，盛水器內(nèi)的雨水便會通過虹吸管排出，盛水器內(nèi)的降雨重量隨之減小，一次降雨虹吸過程中盛水器內(nèi)的降雨重量變化曲線如圖9(a)所示。 由于虹吸排水的速度遠(yuǎn)遠(yuǎn)大于降雨的速度，最終盛水器內(nèi)的雨水還是會被排出盛水器，只是有降雨時(shí)發(fā)生虹吸排水的時(shí)間要比沒有降雨時(shí)的虹吸排水的時(shí)間要長一些[14]。

在發(fā)生虹吸排水的過程中，盛水器處在一邊接水一邊放水的狀態(tài)，那么就會致使Δt這段時(shí)間的降雨會被隨之排走而漏記，這里把漏記的降雨重量記為ΔmN，一次降雨虹吸過程中漏記的降雨重量ΔmN分析圖如圖9 所示。 那么如何找到一種補(bǔ)償算法，得到一次虹吸排水過程中漏記的降雨重量ΔmN，成為提高測量精度的關(guān)鍵所在。

根據(jù)上述分析，本文提出了一種在虹吸發(fā)生時(shí)的測量補(bǔ)償方法，以消除誤差[15]。 首先分析一種理想的情況，假設(shè)在虹吸過程降雨速度和虹吸排水的速度都是勻速的，則虹吸排水曲線就是一條固定斜率的曲線，圖9(a)是實(shí)際降雨過程中盛水器中的降雨重量變化曲線圖。 圖9(b)表示無降雨時(shí)虹吸發(fā)生時(shí)盛水器內(nèi)降雨重量變化曲線圖，從虹吸的機(jī)理可知，圖9(b)降雨重量減小過程t1≤t≤t′2的曲線是盛水器和虹吸管制成后的固有曲線，不會隨降雨量的變化而變化。 將有降雨發(fā)生時(shí)的虹吸過程(即t1≤t≤t2)盛水器內(nèi)的降雨重量變化曲線(即圖9(a)的A-B 線段)減去無降雨發(fā)生虹吸時(shí)盛水器內(nèi)的降雨重量變化線段和該線段的延長線段(圖9(b)中t1≤t≤t′2的實(shí)線和t′2≤t≤t2的點(diǎn)劃線)，即可得到在t1≤t≤t2時(shí)間內(nèi)實(shí)際降雨量對應(yīng)的盛水器降雨重量變化曲線，如圖9(c)虹吸過程漏記的降雨重量曲線圖所示，從而得出一次降雨虹吸過程中漏記的降雨重量ΔmN。

圖9 一次降雨虹吸過程中漏記的降雨重量ΔmN 分析圖

圖9 中m(t)是降雨過程實(shí)際的盛水器內(nèi)降雨重量變化的函數(shù)，m1(t)是無降雨發(fā)生虹吸時(shí)盛水器重量變化的直線和延長段直線的函數(shù)，m2(t)是測量過程中t1≤t≤t2時(shí)間內(nèi)虹吸發(fā)生時(shí)的漏記的降雨重量變化的函數(shù)。 單次測量(其中0≤t≤t2)帶補(bǔ)償?shù)臏y量函數(shù)關(guān)系:

式中:m(t)是由圖2 的測量裝置實(shí)際測出的；m1(t)是在虹吸管和盛水器的結(jié)構(gòu)都定下后，專門標(biāo)定實(shí)驗(yàn)測出其實(shí)線部分，再根據(jù)實(shí)線斜率延長，延長的長度由t2時(shí)刻決定，t2則是由降雨速率的大小決定，降雨速率越大則t2越大，降雨速率越小則t2越小，每一個(gè)單次測量循環(huán)的時(shí)間長度0≤t≤t2是不等的；mN是m(t)的最大值，即虹吸點(diǎn)的降雨重量，其大小也是由虹吸管和盛水器的結(jié)構(gòu)決定。 由于實(shí)際降雨是變化的，m(t)的實(shí)際曲線不一定是直線，m2(t)也不一定是m(t)的O-A線段的平移。 根據(jù)上式，任何瞬時(shí)的降雨量測量表達(dá)式為:

式中:ΔmN(j)是由第j＋1 次發(fā)生虹吸現(xiàn)象時(shí)(j∈Z)，在t1≤t≤t2時(shí)間內(nèi)降雨重量的值，即如圖9(c)所示的ΔmN，由于降雨的速率不是恒定量，一般各ΔmN(j)不等，即:ΔmN(j-1)≠ΔmN(j)≠ΔmN(j＋1)，且j＜0 時(shí)ΔmN(j)= 0；n表示已有n次虹吸完成，n∈Z；K是盛水器內(nèi)的降雨重量轉(zhuǎn)換成降雨量的比例系數(shù)；m單(n＋1)(t)是第n＋1 次的m單(t)，參考式(4)，得到m單(n＋1)(t)的表達(dá)式為:

式(6)表述的是降雨過程中任何一次盛水器內(nèi)降雨重量增加和虹吸排水的循環(huán)，當(dāng)n=0 時(shí)表達(dá)是第一次盛水器內(nèi)降雨重量增加和虹吸排水的過程；t′2(n＋1)≤t≤t2(n＋1)是固有的虹吸排水直線的延長段，t2(n＋1)可通過圖9(a)到達(dá)最低點(diǎn)的時(shí)刻確定，這在測量裝置的MCU 可用相應(yīng)的判斷程序?qū)崿F(xiàn)，得到t2(n＋1)后，由MCU 的程序作虹吸排水直線的延長段，m1(t)在t1(n＋1)≤t≤t2(n＋1)內(nèi)完全確定。 將式(6)的結(jié)果代入式(5)，則可得出帶補(bǔ)償關(guān)系的降雨量測量值h降雨量(t)。

4 實(shí)驗(yàn)數(shù)據(jù)處理與分析

4.1 有無降雨虹吸過程模擬實(shí)驗(yàn)

對上述提出的虹吸補(bǔ)償算法模型進(jìn)行試驗(yàn)驗(yàn)證，首先進(jìn)行降雨模擬實(shí)驗(yàn)， 將雨速控制在0.9408 mm/min，根據(jù)實(shí)際降雨過程和虹吸過程得到的降雨重量數(shù)據(jù)，繪制降雨模擬曲線。 一次降雨虹吸過程盛水器內(nèi)的降雨重量變化曲線如圖10 所示，實(shí)際有降雨過程虹吸現(xiàn)象產(chǎn)生時(shí)的起始重量為mN；虹吸過程持續(xù)時(shí)間為t1～t2。

圖10 一次降雨虹吸過程盛水器內(nèi)的降雨重量變化曲線圖

要想對虹吸過程中漏記的降雨量進(jìn)行補(bǔ)償，首先要得到?jīng)]有降雨時(shí)虹吸過程曲線和有降雨時(shí)虹吸過程曲線的實(shí)驗(yàn)數(shù)據(jù)，如圖11 所示，t′2為無降雨虹吸排水過程的時(shí)間長度，t2為有降雨虹吸排水過程的時(shí)間長度。

圖11 有無降雨虹吸過程降雨重量變化對比曲線圖

從虹吸補(bǔ)償算法的理論分析可知，無降雨時(shí)虹吸過程是一條固有曲線，是由虹吸裝置本身決定的。事實(shí)上，在實(shí)驗(yàn)過程中由于虹吸排水過程可能發(fā)生提前虹吸，會導(dǎo)致虹吸點(diǎn)的重量并不完全相同，相應(yīng)的無降雨時(shí)虹吸過程的曲線也就不是唯一的。 為了解決這一問題，就進(jìn)行了多次無降雨時(shí)的虹吸排水實(shí)驗(yàn)，繪制出多條虹吸排水過程重量變化曲線。 然后將無降雨時(shí)虹吸過程中的虹吸點(diǎn)重量劃分了5 個(gè)不同的區(qū)間，根據(jù)不同的虹吸點(diǎn)重量區(qū)間，挑選出了5 條無降雨時(shí)的虹吸過程曲線，有針對性地對有降雨時(shí)虹吸過程中不同的虹吸點(diǎn)重量區(qū)間進(jìn)行虹吸補(bǔ)償。 無降雨時(shí)虹吸過程的曲線在某一虹吸點(diǎn)重量區(qū)間是唯一確定的。

根據(jù)上述分析，虹吸點(diǎn)的重量并不完全相同，降雨總重量的計(jì)算方法是將每次虹吸點(diǎn)重量和每次虹吸過程漏記的降雨重量進(jìn)行疊加，并考慮降雨結(jié)束時(shí)沒有達(dá)到虹吸點(diǎn)的降雨重量。 實(shí)際降雨過程任何瞬時(shí)的降雨量測量表達(dá)式為:

式中:mN(j)是由第j＋1 次發(fā)生虹吸現(xiàn)象時(shí)虹吸點(diǎn)的降雨重量，j∈Z，一般各mN(j)不等，即:mN(j-1)≠mN(j)≠mN(j＋1)，且j＜0 時(shí)mN(j)= 0。

4.2 虹吸補(bǔ)償算法應(yīng)用

由于實(shí)際得到的無降雨時(shí)虹吸過程曲線和有降雨時(shí)虹吸過程曲線是非線性的，而圖9 的虹吸補(bǔ)償算法的理論分析是建立在虹吸過程曲線是線性的基礎(chǔ)上，那么現(xiàn)在需要將虹吸補(bǔ)償算法推廣到單調(diào)非線性的實(shí)際虹吸過程。

首先需要對實(shí)際得到的有無降雨時(shí)虹吸過程的曲線進(jìn)行分析，有降雨時(shí)虹吸過程盛水器內(nèi)降雨重量變化隨雨速的大小而不同，可以根據(jù)實(shí)際降雨數(shù)據(jù)通過程序算法進(jìn)行處理。 有降雨時(shí)虹吸過程的時(shí)間長度為t2，t2可以通過降雨過程降雨重量達(dá)到最低點(diǎn)時(shí)得到，通過最小二乘法對無降雨虹吸過程的曲線進(jìn)行多項(xiàng)式擬合，得到L1(t)，將L1(t)按照自身的函數(shù)關(guān)系延長至t2，得到帶補(bǔ)償關(guān)系的虹吸過程降雨重量變化曲線，令其為L′1(t)。 此虹吸補(bǔ)償方法記為算法A，如圖12 所示。

圖12 算法A 帶補(bǔ)償關(guān)系的虹吸過程降雨重量變化曲線圖

將L(t)和L′1(t)兩者相減的結(jié)果便是一次虹吸過程漏記的降雨重量變化曲線L2(t)，以上算法可以通過軟件編程實(shí)現(xiàn)。 通過算法A 得到的虹吸過程漏記的降雨重量曲線如圖13 所示。

圖13 算法A 求得虹吸過程漏記的降雨重量曲線圖

為了提高測量精度，找到最優(yōu)的補(bǔ)償算法。 經(jīng)分析，由于實(shí)際的虹吸過程是一條曲線，而虹吸補(bǔ)償算法的理論分析是建立在虹吸過程是一條直線的基礎(chǔ)上。 因此，提出了一種新的數(shù)據(jù)處理算法對無降雨時(shí)虹吸過程的曲線L1(t)進(jìn)行處理，首先找到無降雨時(shí)虹吸結(jié)束時(shí)間t′2，過點(diǎn)t′2并做其切線延長至t2，得到帶補(bǔ)償關(guān)系的虹吸過程降雨重量變化曲線，令其為，此虹吸補(bǔ)償方法記為算法B，如圖14 所示。

圖14 算法B 帶補(bǔ)償關(guān)系的虹吸過程降雨重量變化曲線圖

將L(t)和L″1(t)兩者相減的結(jié)果便得到算法B一次虹吸過程漏記的降雨重量變化曲線L′2(t)，以上算法可以通過軟件編程實(shí)現(xiàn)。 通過算法B 得到的虹吸過程漏記的降雨重量曲線如圖15 所示。

圖15 算法B 求得虹吸過程漏記的降雨重量曲線圖

4.3 實(shí)驗(yàn)結(jié)果分析

利用降雨模擬試驗(yàn)場地模擬實(shí)際降雨的過程，調(diào)節(jié)水閥使得模擬降雨速度在0.1 mm/min ～1.0 mm/min的范圍內(nèi)，通過實(shí)驗(yàn)測得各個(gè)降雨速度下降雨虹吸過程的數(shù)據(jù)。 然后，將這些數(shù)據(jù)利用虹吸補(bǔ)償算法A、B 求出虹吸過程中漏記的降雨重量ΔmN，通過MCU 程序計(jì)算得出一次降雨虹吸過程中的降雨量，并與實(shí)際的降雨虹吸過程中的降雨量進(jìn)行對比，比較兩種算法之間哪個(gè)誤差更小。

實(shí)驗(yàn)選取了大雨、暴雨、特大暴雨三個(gè)降雨等級的雨速進(jìn)行降雨模擬實(shí)驗(yàn)，將降雨速度設(shè)定在0.157 mm/min、 0. 314 mm/min、 0. 470 mm/min、0.627 mm/min、0.862 mm/min、0.941 mm/min，對這6 組實(shí)驗(yàn)樣本數(shù)據(jù)采用虹吸補(bǔ)償算法A、B 后，求得在一次降雨虹吸過程中(包括水位上升和虹吸排水過程)的降雨量，降雨量的誤差分析如表1 所示。

表1 一次降雨虹吸過程中降雨量的誤差分析

從實(shí)驗(yàn)結(jié)果可知，在一次降雨虹吸過程中，不作虹吸補(bǔ)償算法的降雨量測量誤差范圍為2.682%～8.147%，經(jīng)過虹吸補(bǔ)償算法A、B 補(bǔ)償后降雨量測量誤差明顯減小，其中算法A 的最大測量誤差為2.127%，算法B 的最大測量誤差為0.704%。 通過對虹吸補(bǔ)償算法A、B 對比可知，算法B 的補(bǔ)償效果較好，其測量誤差范圍在1%以內(nèi)。 根據(jù)以上分析可得，虹吸補(bǔ)償算法有效地消除了降雨隨虹吸排水造成的計(jì)量誤差，提高了降雨量的測量精度。

5 結(jié)論

為了解決傳統(tǒng)虹吸式雨量計(jì)虹吸時(shí)產(chǎn)生的計(jì)量誤差，提高降雨量的測量精度，實(shí)現(xiàn)降雨自動(dòng)化監(jiān)測，本文基于稱重傳感器研制了一種降雨量測量裝置，并融合工業(yè)物聯(lián)網(wǎng)技術(shù)開發(fā)了一套精密在線監(jiān)測系統(tǒng)。 首先選擇內(nèi)部結(jié)構(gòu)為一個(gè)惠斯通電橋的電阻應(yīng)變式稱重傳感器，設(shè)計(jì)出弱小信號差分放大電路及基于單片機(jī)的數(shù)據(jù)采集處理系統(tǒng)，建立了虹吸補(bǔ)償算法數(shù)學(xué)模型，通過對降雨重量信號的實(shí)時(shí)采集及模型計(jì)算，獲取降雨量數(shù)據(jù)。 結(jié)果表明:將虹吸補(bǔ)償算法數(shù)學(xué)模型應(yīng)用于降雨量測量，有效地消除了降雨隨虹吸排水造成的計(jì)量誤差，降雨測量誤差在1%以內(nèi)，比傳統(tǒng)虹吸式雨量計(jì)的測量精度更高。該裝置具有供電方便、通信穩(wěn)定、測量精度高等優(yōu)勢，在降雨量監(jiān)測領(lǐng)域有極大的使用價(jià)值。