藺 瀟

（北京交通大學(xué) 機(jī)械與電子控制工程學(xué)院，北京100044）

雨量的測量在氣象觀測中是一個非常重要的項目， 及時和準(zhǔn)確地獲取雨量數(shù)據(jù)更是對防治洪澇、 干旱以及山體滑坡等自然災(zāi)害具有重要的作用，與人類生活、生產(chǎn)和建設(shè)事業(yè)密切相關(guān)。為了提高測量降雨的時間分辨率，許多自動氣象站都安裝有翻斗式雨量計、虹吸式雨量計等。

雨量計觀測降水具有測量精度高的優(yōu)點(diǎn)，但是對于觀測大面積范圍內(nèi)的降水事件，尤其在人煙稀少、難以廣泛布設(shè)的地區(qū)，其應(yīng)用受到一定限制[1-2]。壓電式雨量傳感器是一種低成本、低功耗的、能夠同步監(jiān)測雨量及降雨微觀特征的雨量傳感器，其輸出信號幅值與雨滴下落的末速度和雨滴的質(zhì)量相關(guān)。但是，由于雨滴下落末速度與雨滴質(zhì)量的非線性關(guān)系以及傳感器結(jié)構(gòu)尺寸的公差，壓電式傳感器的輸出信號幅值和雨滴體積的關(guān)系也是非線性的，因此壓電式雨量傳感器安裝使用前一般需要進(jìn)行準(zhǔn)確標(biāo)定。

1 壓電式雨量傳感器的特點(diǎn)

壓電式雨量傳感器是將雨滴下落的沖擊力信號轉(zhuǎn)換為電信號，間接測量雨滴的粒徑大小。通過采集處理壓電傳感元件的輸出信號，可以計算得到單個雨滴的粒徑， 以及一段時間內(nèi)的雨滴譜分布、降雨強(qiáng)度和累計降雨量，具有實時性好、自動化程度高、成本低廉等優(yōu)點(diǎn)。但是，由于每個壓電式雨量傳感器機(jī)械結(jié)構(gòu)的細(xì)微差別，對輸出信號具有較大的影響， 每個傳感器在安裝使用前必須進(jìn)行標(biāo)定，建立最佳校準(zhǔn)曲線，才能取得較好的測量精度。

現(xiàn)有的壓電式雨量傳感器標(biāo)定方法，按照原理分為直接標(biāo)定法[3]、等動量法[4]、自適應(yīng)標(biāo)定法[5]。直接標(biāo)定法和等動量法采用傳統(tǒng)的單個雨滴的標(biāo)定方法，前者對于試驗條件的需求較為苛刻，后者則依托于在傳感器使用環(huán)境下可以精準(zhǔn)計算水滴動量理論公式，不能保證標(biāo)定結(jié)果在不同地區(qū)、不同氣象環(huán)境下都具有較好的測量精度。

自適應(yīng)標(biāo)定法是采用翻斗雨量計的輸出信號作為標(biāo)準(zhǔn)雨量，對壓電傳感器的輸出信號進(jìn)行自適應(yīng)訓(xùn)練和優(yōu)化，得到試驗條件下最優(yōu)的傳感器響應(yīng)模型，實現(xiàn)不同粒徑雨滴對應(yīng)傳感器輸出電壓的測量[6]。以往的自適應(yīng)標(biāo)定法，在樣本數(shù)據(jù)較少時會產(chǎn)生標(biāo)定效果較差的問題。文中介紹了以翻斗雨量計輸出的差分信號為每一時刻的參考雨量的自適應(yīng)標(biāo)定方法， 增加了同樣試驗條件下的樣本數(shù)據(jù)，在減少了標(biāo)定時間的條件下實現(xiàn)了壓電式雨量傳感器的標(biāo)定，并在實際降雨測量中對標(biāo)定結(jié)果進(jìn)行了試驗驗證。

2 壓電式雨量傳感器的設(shè)計

壓電式雨量傳感器由外殼、 壓電傳感元件、信號采集電路等3 部分組成。傳感器外殼上蓋是直徑120 mm 的圓形金屬蓋，當(dāng)受到雨滴撞擊時，壓電傳感元件會將雨滴撞擊傳感器外殼產(chǎn)生的機(jī)械振動轉(zhuǎn)換為電信號，雨滴撞擊的沖擊力越大，則振動幅值越大，傳感元件輸出電壓的幅值也越大。壓電傳感元件的輸出信號經(jīng)過采集電路處理，然后對照雨滴粒徑和傳感器輸出的標(biāo)定關(guān)系，就可以得到該雨滴的粒徑大小， 然后通過232 總線輸出測量結(jié)果。這種降雨測量儀器， 能夠連續(xù)測量并記錄雨滴粒徑，同時穩(wěn)定獲得降雨強(qiáng)度和某時段的降雨量。

傳感器的外殼由不銹鋼制成，用于接收雨滴的沖擊產(chǎn)生機(jī)械振動。外殼上面的面積就是能采集到雨滴的采樣面積，采樣面積越大，最后的測量結(jié)果受統(tǒng)計誤差的印象就越小， 然而采樣面積越大，同時采集到2 粒以上的雨滴撞擊產(chǎn)生的信號的次數(shù)也就越多，這也會導(dǎo)致測量結(jié)果的不準(zhǔn)確。外殼上表面稍帶有弧度，使得傳感器上的雨滴在重力作用下流走，可以防止在采樣區(qū)域內(nèi)積水較多使雨滴下落的作用力被積水所緩沖。

壓電傳感元件采用圓形PZT 壓電陶瓷片，多個壓電傳感元件并聯(lián)構(gòu)成陣列，使得有效傳感區(qū)域獲得較為均勻的振動響應(yīng)，減小邊界效應(yīng)對測量的影響，減小單個雨滴的測量誤差。傳感器內(nèi)部結(jié)構(gòu)如圖1所示。

圖1 傳感器內(nèi)部結(jié)構(gòu)模型Fig.1 Sensor internal structure model

3 自適應(yīng)標(biāo)定算法

自適應(yīng)標(biāo)定法以雨量計的輸出雨量為參考，將基于傳感器模型的標(biāo)定誤差曲面函數(shù)進(jìn)行自適應(yīng)迭代優(yōu)化，生成傳感器校準(zhǔn)曲線，實現(xiàn)對傳感器輸出的標(biāo)定。自適應(yīng)標(biāo)定法對壓電式雨量傳感器進(jìn)行標(biāo)定，可以避免直接標(biāo)定法、等動量法中所需要的苛刻的試驗條件，以及實際應(yīng)用環(huán)境與標(biāo)定試驗環(huán)境的差別所引起的測量誤差。

根據(jù)經(jīng)驗?zāi)Ｐ?，雨滴體積與傳感器的響應(yīng)關(guān)系不妨設(shè)為

式中：v（n）為第n 個雨滴的體積；x（n）為對應(yīng)的傳感器輸出脈沖幅值；α，β 為自適應(yīng)系數(shù)。設(shè)在采樣時間T 內(nèi)采集到的雨滴數(shù)量為N，所有雨滴的體積和為

采用翻斗雨量計作為標(biāo)定的參考信號源。翻斗雨量計在T 時刻共翻轉(zhuǎn)m 次，tm為第m 次翻轉(zhuǎn)的時刻， 每次翻轉(zhuǎn)信號之間的雨量為近似線性增長，同時刻翻斗雨量計測得的雨滴體積和可以通過時間差分計算，即

壓電傳感器輸出與翻斗雨量計輸出二者做差則可以得到測量誤差曲面函數(shù)，表示不同自適應(yīng)系數(shù)下2 種傳感器的測量誤差為

然后通過最優(yōu)參數(shù)搜索算法就可以計算出E最小時的自適應(yīng)系數(shù)，從而實現(xiàn)傳感器標(biāo)定。在此采用梯度下降法。

通過梯度下降法計算最優(yōu)的自適應(yīng)參數(shù)首先要對自適應(yīng)參數(shù)賦值。該值可以是隨機(jī)的，也可以是一個全零的向量，然后改變其值，使得E 按梯度下降的方向減少。自適應(yīng)系數(shù)向量c 為

步長μ 則通過一維搜索確定， 步長乘以梯度，得到當(dāng)前位置下降的距離， 然后驗證終止條件，判斷當(dāng)前自適應(yīng)系數(shù)是否符合目標(biāo)誤差e， 符合則迭代終止，不符合則更新參數(shù)。重新計算自適應(yīng)系數(shù)向量，進(jìn)入下一次迭代：

4 試驗及數(shù)據(jù)處理過程

水利水電科學(xué)研究院八達(dá)嶺實驗基地人工模擬降雨大廳為傳感器標(biāo)定數(shù)據(jù)的采集試驗提供了所需的模擬降雨環(huán)境。該模擬降雨大廳的有效降雨高度為12 m，降雨強(qiáng)度連續(xù)變化范圍為10～200 mm/h，降雨強(qiáng)度恒定，降雨均勻度＞0.8，總體雨滴直徑變化范圍0.5～4.3 mm，能夠較好地模擬恒定降雨強(qiáng)度的降雨環(huán)境。試驗示意圖如圖2所示。試驗中，使用的翻斗雨量計為精密型雙翻斗式雨量計，分辨力0.2 mm，降雨強(qiáng)度測量范圍0.01～4 mm/min，測量準(zhǔn)確度≤±3%。

圖2 傳感器標(biāo)定試驗示意圖Fig.2 Schematic diagram of sensor calibration experiment

在傳感器標(biāo)定試驗過程中，壓電雨量傳感器和翻斗雨量計對模擬降雨進(jìn)行同步采集，輸出信號并經(jīng)過采集電路后數(shù)據(jù)傳輸?shù)缴衔粰C(jī)，由上位機(jī)計算得到標(biāo)定結(jié)果。壓電傳感器的輸出電壓如圖3所示，圖中信號采樣頻率為1 kHz。

圖3 壓電傳感器輸出的電壓數(shù)據(jù)Fig.3 Voltage data from piezoelectric sensors

計算中，首先提取每個電壓脈沖的峰值，每一個脈沖的峰值對應(yīng)一個被采樣的雨滴，電壓峰值提取效果如圖4所示。將電壓值代入式（1）和式（2），可以得到壓電傳感器測量的累計雨量變化。翻斗雨量計輸出雨量累計如圖5所示，由式（3）可以將階躍信號轉(zhuǎn)換為連續(xù)信號，得到翻斗雨量計測得的每一時刻的累計雨量。再把2 個雨量數(shù)據(jù)代入式（4），就可以得到誤差曲面函數(shù)，進(jìn)而計算出當(dāng)前降雨環(huán)境的自適應(yīng)系數(shù)α，β。

圖4 壓電傳感器輸出的電壓峰值提取結(jié)果（局部）Fig.4 Results of peak voltage extraction from piezoelectric sensor（local）

圖5 翻斗雨量計輸出Fig.5 Output of dump rain gauge

將使用梯度下降法得到的標(biāo)定系數(shù)α，β 代入傳感器模型， 得到壓電式雨量傳感器的標(biāo)定曲線。圖6為一次試驗中經(jīng)過標(biāo)定后的壓電式傳感器輸出數(shù)據(jù)與雨量計輸出數(shù)據(jù)的對比。

圖6 標(biāo)定結(jié)果Fig.6 Calibration results

調(diào)節(jié)降雨大廳控制臺的設(shè)置，在不同的降雨強(qiáng)度下重復(fù)以上試驗步驟，可以得到不同的降雨強(qiáng)度下的標(biāo)定系數(shù)，具體見表1。

表1 不同雨強(qiáng)下的標(biāo)定系數(shù)Tab.1 Calibration coefficient under different rainfall intensities

由表可知參數(shù)的規(guī)律性變化，自適應(yīng)系數(shù)β 穩(wěn)定在0.6 左右，不會隨著雨量變化而大幅變化，但是α 會隨著降雨強(qiáng)度的增大而增大。所以，最終的標(biāo)定結(jié)果若應(yīng)用在不同降雨強(qiáng)度的環(huán)境下，需要結(jié)合已知電壓的累加和來自適應(yīng)系數(shù)α 進(jìn)行修正，才可能得到較為準(zhǔn)確的測量結(jié)果。所以，通過最小二乘法擬合得到電壓的累加和與自適應(yīng)系數(shù)α 的關(guān)系（如圖7所示）， 從而得到不同降雨強(qiáng)度下的標(biāo)定系數(shù)以及最終的標(biāo)定曲線。

圖7 電壓累加和與系數(shù)α 的擬合曲線Fig.7 Curve fitting between accumulative sum of voltage and coefficient α

將最終的標(biāo)定結(jié)果代入不同降雨強(qiáng)度的試驗數(shù)據(jù)，與翻斗雨量計數(shù)據(jù)進(jìn)行比對，可以看到最終標(biāo)定的效果（如圖8所示）。除圖8（a）有較大的誤差，主要源自于翻斗雨量計的沾濕誤差，壓電傳感器的測量結(jié)果與翻斗雨量計能較好地貼合， 其他3組最大誤差均小于6%。

圖8 不同降雨強(qiáng)度下的標(biāo)定結(jié)果Fig.8 Calibration results under different rainfall intensities

5 試驗及其結(jié)果

使用該壓電式雨量傳感器測量， 對2018年7月24日北京市海淀區(qū)某地的一次實際降雨進(jìn)行測量， 并將輸出結(jié)果與翻斗雨量計輸出進(jìn)行對比，由此驗證標(biāo)定曲線。將壓電傳感器和翻斗雨量計同時放置在降雨環(huán)境下，采集二者的輸出信號。在10 min內(nèi)，壓電式傳感器輸出數(shù)據(jù)與雨量計輸出降雨量的對比如圖9所示。

圖9 實際降雨的測量結(jié)果Fig.9 Measured results of actual rainfall

由圖可見，二者具有相似的走勢，而且測量過程中最大的雨量測量誤差＜10%，最終的累積雨量的測量誤差＜6%。試驗中傳感器誤差的主要來源有：

1）標(biāo)定試驗中測量都是在無風(fēng)的試驗環(huán)境下進(jìn)行的，沒有討論風(fēng)速、氣溫、氣壓等其他環(huán)境因素對測量結(jié)果的影響，而實際測量中其他氣象環(huán)境會造成測量誤差；

2）由于傳感器結(jié)構(gòu)和支撐影響，傳感器不同區(qū)域的振動特性不同導(dǎo)致了傳感器輸出特性不均勻，盡管采用多個壓電陶瓷片組成并聯(lián)陣列來減小這種情況所產(chǎn)生的影響， 但是仍然可能產(chǎn)生測量誤差，而且1 個雨滴有可能對傳感器產(chǎn)生2 次沖擊，2 次沖擊的影響會使得雨量測量結(jié)果明顯偏大；

3）由于翻斗雨量計具有滯后性且易受到沾濕誤差的影響，翻斗雨量計的輸出信號具有不確定的延時，無法與壓電傳感器的實時信號相對應(yīng)，使得后續(xù)計算中難以實現(xiàn)翻斗雨量計與壓電式傳感器輸出信號的實時對應(yīng)。

6 結(jié)語

采用基于壓電效應(yīng)原理的降雨量測量傳感器，基于自適應(yīng)信號處理方法，以翻斗雨量計輸出信號為基準(zhǔn)， 對傳感器標(biāo)定方法進(jìn)行了改進(jìn)和試驗驗證，并分析了主要測量誤差的來源。試驗結(jié)果，表明該方法能替代傳統(tǒng)的直接標(biāo)定方法，簡化標(biāo)定試驗的硬件需求，為今后在氣象站及野外環(huán)境的降雨檢測中的應(yīng)用提供了參考價值。