李 想,江朝暉,陸元洲,潘 煒,余林生
(1.安徽農(nóng)業(yè)大學 農(nóng)業(yè)信息學安徽省重點實驗室,安徽 合肥230036;2.安徽農(nóng)業(yè)大學 蜂業(yè)研究所,安徽 合肥230036)
我國是世界傳統(tǒng)養(yǎng)蜂大國,蜂群擁有量位居世界第一位[1]。溫度是蜜蜂生命活動的重要生態(tài)因子,蜂巢溫度的變化與蜜蜂個體生長發(fā)育、蜜蜂群勢強弱、蜜蜂產(chǎn)卵繁殖等密切相關(guān)[2~4]。人們不斷嘗試新的理論、方法和技術(shù)手段研究蜂巢內(nèi)溫度分布和變化規(guī)律、揭示其恒溫調(diào)控機理,并取得了部分成果[5]。在蜂巢監(jiān)測和分析方面,湖南大學采用ARM9+uITRON 平臺和無線GPRS 模塊,研制了蜜蜂蜂箱專用數(shù)字化監(jiān)測系統(tǒng)[6];云南農(nóng)業(yè)大學研制的蜂巢溫濕度數(shù)據(jù)采集與分析處理系統(tǒng),通過串口接收溫濕度數(shù)據(jù),對接收的數(shù)據(jù)進行檢測與分析[7]。總的來說,這些系統(tǒng)采集點過少很難涵蓋蜂巢整個立體空間,并且數(shù)據(jù)分析僅僅限于最大值、最小值和均值等,不能直觀反映蜂巢內(nèi)不同位置間的聯(lián)系。
本文采用微型溫度傳感器陣列、無線通信和可視化軟件技術(shù),設(shè)計一種多通道蜂巢溫度監(jiān)測與分析系統(tǒng),集采集、顯示和分析等功能于一體,具有便捷、直觀和實時的特點。
系統(tǒng)主要分為微型溫濕度傳感器陣列模塊、多通道無線數(shù)據(jù)采集模塊、數(shù)據(jù)實時顯示模塊、數(shù)據(jù)動態(tài)分析模塊(圖1)。計算機通過無線通信模塊發(fā)送指令給遠端數(shù)據(jù)采集模塊,數(shù)據(jù)采集模塊收到指令后并控制傳感器陣列模塊采集數(shù)據(jù),MCU 控制器處理數(shù)據(jù)后,再由無線模塊發(fā)送給遠端計算機,計算機接收到數(shù)據(jù)后對數(shù)據(jù)進一步的處理并實時顯示和動態(tài)分析蜂巢內(nèi)溫度變化情況。
圖1 系統(tǒng)總體架構(gòu)Fig 1 Total architecture of system
蜂巢空間狹小,能利用的空間非常有限,為減小對蜜蜂正常生活的影響和保證數(shù)據(jù)采集的真實準確,所選取的傳感器必須非常微小。本系統(tǒng)選取的是Sensirion 公司推出SHT21 傳感器,是目前世界上最小的數(shù)字溫濕度傳感器之一。蜂巢是一個空間立體的區(qū)域,單點或簡單多點采集很難反映蜂巢溫度場動態(tài)變化情況。系統(tǒng)均勻選取了蜂巢空間30 個點作為蜂巢溫度場代表點,選取蜂巢外一個點作為環(huán)境溫度代表點。傳感器陣列基本布局是5 行6 列,每行傳感器均勻的分布在矩形平面對角線上(圖2)。整個監(jiān)測裝置,包括箱體、蜂箱蓋、主控板、固定板、巢脾、外部傳感器和內(nèi)部傳感器陣列[8]。主控板貼合固定板上,外部傳感器和內(nèi)部傳感器陣列垂直與主控板相連(圖3),多個巢脾放置于箱體內(nèi),均勻分割傳感器陣列。
圖2 微傳感器陣列布局Fig 2 Layout of micro-sensor array
圖3 微傳感器陣列實物圖Fig 3 Physical map of micro-sensor array
系統(tǒng)選用鼎泰克公司RS—232 轉(zhuǎn)Zig Bee 無線通信模塊,采用RS—232 串行通信方式。微傳感器陣列一共擁有31 只微型溫濕度傳感器,每個傳感器都編有唯一的ID 編號,其中,1~10 編號由MCU1 控制,11~21 編號由MCU2 控制,22~31 編號由MCU3 控制。每個MCU 擁有一個無線節(jié)點,計算機通過無線節(jié)點發(fā)送指令后,蜂巢的三個無線節(jié)點都會收到指令。收到指令后,MCU 會立即采集傳感器溫度濕度數(shù)據(jù)并準備隨時發(fā)送,當MCU1 通過無線節(jié)點一發(fā)送溫度濕度數(shù)據(jù)后,同時MCU2 和MCU3 處于等待之中。當計算機收到數(shù)據(jù)后,向節(jié)點發(fā)送成功指令,MCU2 接收到成功指令后,MCU2 會立即發(fā)送準備好的溫度濕度數(shù)據(jù),以此類推,直至發(fā)送數(shù)據(jù)結(jié)束(圖4)。
圖4 多通道無線采集流程圖Fig 4 Flow chart of multichannel wireless collection
系統(tǒng)軟件采用Matlab GUI 設(shè)計實現(xiàn),系統(tǒng)軟件主要分為數(shù)據(jù)存儲模塊、數(shù)據(jù)實時顯示模塊、數(shù)據(jù)分析模塊。數(shù)據(jù)存儲模塊主要功能是對數(shù)據(jù)識別和存儲,每個傳感器都編有唯一的ID,計算機收到數(shù)據(jù)后對ID 進行識別并加上時間信息對數(shù)據(jù)進行存儲。數(shù)據(jù)實時顯示模塊是蜂巢溫度監(jiān)測系統(tǒng)一個很重要的模塊,當無線傳感器采集的溫濕度更新相應(yīng)的Excel 表后,軟件系統(tǒng)會自動更新相應(yīng)的數(shù)據(jù),對溫濕度實時顯示。系統(tǒng)采用三維立體化的顯示方式,在系統(tǒng)主界面中,可以直觀地看到30 個采集點,并且每個采集點都進行了相應(yīng)編號,與實物相對應(yīng)。
單通道分析主要包括各個采集點的均值、方差以及直方圖,均值主要是反映在額定的時間段內(nèi)溫度場中溫度的分布情況,方差主要是體現(xiàn)額定時間段內(nèi)各采集點溫度變化幅度,直方圖主要是為了了解一組溫度數(shù)據(jù)的分布特征。
蜂巢溫度場是個動態(tài)立體空間區(qū)域,單通道數(shù)據(jù)分析很難反映蜂巢內(nèi)整體溫度分布和相互關(guān)聯(lián)情況。為了對整個蜂巢系統(tǒng)進行立體的分析,系統(tǒng)加入了動態(tài)的多通道數(shù)據(jù)分析。
相關(guān)分析(correlation analysis)是研究數(shù)據(jù)之間是否存在某種依存關(guān)系
偏相關(guān)分析(partial correlation analysis)是在消除其他因素影響的基礎(chǔ)上求相關(guān)
其中,ε1和ε2分別為y=αu+βux2+ε1和y=αu+βux2+ε2的殘部。
選用意大利蜜蜂,在蜂巢中加入足夠的飼料,保證蜜蜂能夠安全地度過低溫環(huán)境。每次實驗從蜂場開始,蜂場數(shù)據(jù)采集完后運往低溫冷庫中采集,冷庫采集完后放回蜂場,完成一個循環(huán)周期,每周期采集10 天左右。經(jīng)過幾個月的采集,系統(tǒng)性能穩(wěn)定,得到過了大量的數(shù)據(jù)。部分數(shù)據(jù)如圖5所示,可以直觀看到蜂巢中心位置和其它位置隨環(huán)境溫度變化情況。
對圖5(a)數(shù)據(jù)進行單通道直方圖分析,結(jié)果如圖6 所示,分別呈現(xiàn)了中心點、第30 點、蜂巢外三個位置的溫度數(shù)據(jù)分布情況。
對圖5(a)數(shù)據(jù)進行多通道相關(guān)分析,可得到互相關(guān)系數(shù)矩陣和偏相關(guān)系數(shù)矩陣。為了更加形象地反映蜂巢內(nèi)各點之間關(guān)聯(lián)程度,將互相關(guān)系數(shù)、偏相關(guān)系數(shù)進行了三維立體顯示(圖6)。
通過定量的單通道和多通道分析,初步得出以下結(jié)論:1)在環(huán)境溫度劇烈變化時,蜂巢中心點溫度波動范圍很小,基本保持在35 ℃左右;2)在環(huán)境溫度驟烈變化時,蜂巢內(nèi)越靠近邊緣的位置,其溫度變化越明顯,但最終會趨于一個穩(wěn)定值;3)蜂巢中心點的周圍、蜂巢邊緣位置間的相關(guān)性較大。
圖5 實驗采集的蜂巢溫度原始數(shù)據(jù)Fig 5 Original data of beehive temperature acquired in experiment
圖6 單通道直方圖分析Fig 6 Histogram analysis of single channel
本文采用SHT21 微型溫濕度傳感器和Zig Bee 無線通信技術(shù),實時采集蜂箱內(nèi)立體均勻分布的30 點溫度和箱外環(huán)境溫度,應(yīng)用Matlab GUI 構(gòu)建多通道數(shù)據(jù)接收、可視化顯示和動態(tài)分析軟件。經(jīng)過長時間實驗和測試,系統(tǒng)性能良好,單通道和多通道分析表明:蜂群為適應(yīng)溫度變化自主地通過抱團聚集來保持中心溫度相對穩(wěn)定,且溫度分布情況與位置具有較大關(guān)聯(lián)性,當環(huán)境溫度劇烈變化時,蜂巢中心溫度為34.8℃左右,而邊緣溫度變化極大。本系統(tǒng)能精確獲取蜂巢內(nèi)溫度分布及其變化規(guī)律,有助于揭示蜂巢恒溫調(diào)控機理,對蜂學研究和蜜蜂養(yǎng)殖具有重要意義。
圖7 空間位置—相關(guān)系數(shù)圖Fig 7 Image of spatial position vs correlation coefficient
[1] 陳瑪琳,趙芝俊,席桂萍.中國蜂產(chǎn)業(yè)發(fā)展現(xiàn)狀及前景分析[J].浙江農(nóng)業(yè)學報,2014(3):825-829.
[2] 汪天澍,劉 芳,余林生,等.蜜蜂蜂群溫濕度調(diào)節(jié)的研究進展[J].生態(tài)學報,2015(10):3172-3179.
[3] Cao Y F,Yu L S,Bi S D,et al.Studying progress of the influence of temperature to honeybee[J].Journal of Bee,2007(4):13-15.
[4] Tautz J,Maier S,Groh C,et al.Behavioral performance in adult honey bees is influenced by the temperature experienced during their larval development[C]∥Proceedings of the National Academy of Sciences of the United States of America,2003:7343-7347.
[5] 江朝暉,張 靜,楊春合,等.越冬蜂群調(diào)溫機理仿真[J].系統(tǒng)仿真學報,2014(6):1301-1307.
[6] 譚慶忠.蜜蜂蜂箱專用數(shù)字化監(jiān)測系統(tǒng)設(shè)計與實現(xiàn)[D].長沙:湖南大學,2012.
[7] 王銳剛,和紹禹,黃 峻,等.蜂巢溫濕度數(shù)據(jù)采集與分析處理系統(tǒng)[J].云南農(nóng)業(yè)大學學報,2007(5):699-705.
[8] 潘 煒,付月生,王遠飛,等.用于研究蜂巢恒溫機制的溫濕度測試裝置:中國,CN203629605U[P].2014—06—04.