日光溫室冠層特征溫度時(shí)空變化規(guī)律分析

　　摘要：針對(duì)日光溫室冬季反季節(jié)生產(chǎn)中冠層溫度時(shí)空分布不均的問(wèn)題。以西北型日光溫室為研究對(duì)象,建立40通道PT100溫度場(chǎng)監(jiān)測(cè)系統(tǒng)，實(shí)現(xiàn)了冬季番茄冠層全天候動(dòng)態(tài)溫度數(shù)據(jù)采集。在此基礎(chǔ)上，采用克里金插值算法進(jìn)行冠層溫度場(chǎng)建模，通過(guò)差分進(jìn)化算法獲取溫度極值點(diǎn)，分析了不同天氣條件下冠層特征溫度的時(shí)空變化規(guī)律。結(jié)果表明，晴天、多云、陰(雨)天插值驗(yàn)證的均大于0.9，平均均方根誤差分別為1.34℃、0.95℃、0.40℃，算法更適用于陰(雨)天及夜間低溫冠層溫度場(chǎng)的插值。不同天氣條件下，溫室冠層溫度全天整體呈現(xiàn)西高東低、內(nèi)高外低趨勢(shì)，陰(雨)天溫室內(nèi)溫度整體變化規(guī)律趨于一致，晴天、多云天氣揭被后受外界光輻射等因素影響室內(nèi)溫度分布差異性較大，且晴天的夜間溫度下降程度大于陰(雨)天。進(jìn)一步分析冠層極值點(diǎn)分布結(jié)果表明，在不同天氣條件下，極值特征點(diǎn)在日光溫室中分布區(qū)域基本相同，冠層最高溫點(diǎn)主要出現(xiàn)在溫室中部[22.0m，2.5m]附近，最低溫點(diǎn)主要出現(xiàn)在溫室東部外膜[4.0m，5.48m]附近。日光溫室冠層極值特征點(diǎn)的獲取為溫室栽培、溫度災(zāi)害監(jiān)測(cè)與傳感器部署提供了理論基礎(chǔ)。

　　關(guān)鍵詞：日光溫室;冠層溫度場(chǎng);時(shí)空變化規(guī)律;克里金插值;差分進(jìn)化算法

　　引言

　　冠層溫度是作物生長(zhǎng)發(fā)育過(guò)程中重要環(huán)境因素，直接影響作物的生理狀態(tài)[1]。冬季北方反季節(jié)生產(chǎn)過(guò)程中，溫室溫度的監(jiān)測(cè)和調(diào)控尤為重要，日光溫室通過(guò)后墻對(duì)太陽(yáng)能的蓄放熱特性及保溫覆蓋物的保溫，能有效實(shí)現(xiàn)夜間溫室內(nèi)溫度的提升，在中國(guó)北方地區(qū)廣泛分布[2,3]。但隨著生產(chǎn)需求及建造技術(shù)的提升，日光溫室建造朝著大跨度、大長(zhǎng)度方向發(fā)展，室內(nèi)空間不斷擴(kuò)大[4,5]，溫室內(nèi)溫度不僅受外界環(huán)境影響，同時(shí)受墻體、地溫、作物、濕度等因素影響，導(dǎo)致溫室內(nèi)溫度時(shí)空分布不均的問(wèn)題凸顯[6,7]，特別是受冬季連續(xù)陰(雨)天、雪天等氣候季節(jié)性變化影響時(shí)，極易導(dǎo)致作物冠層局部溫度過(guò)低，造成作物冷害、凍害甚至死亡[8,9]。

　　因此，分析作物冠層溫度分布特征，進(jìn)而尋找溫度最低點(diǎn)，為優(yōu)化日光溫室環(huán)境監(jiān)測(cè)及增溫設(shè)施部署提供理論支撐，對(duì)減少凍害損失具有重要意義。目前針對(duì)日光溫室溫度變化規(guī)律、分布特征及監(jiān)測(cè)預(yù)警已開展大量研究[1011。

　　現(xiàn)有日光溫室小環(huán)境因子建模預(yù)測(cè)方法已具有較高精度，能準(zhǔn)確模擬、預(yù)測(cè)監(jiān)測(cè)點(diǎn)的變化規(guī)律[1215]，然而以單點(diǎn)或少數(shù)監(jiān)測(cè)點(diǎn)數(shù)據(jù)對(duì)溫室內(nèi)整體及冠層溫度進(jìn)行分析存在一定局限性，易忽略全局極值溫度點(diǎn)的分析;近年來(lái)，CFD技術(shù)廣泛應(yīng)用于農(nóng)業(yè)，能動(dòng)態(tài)模擬太陽(yáng)熱負(fù)荷在溫室中的分布與變化，結(jié)合溫室構(gòu)造等因素可實(shí)現(xiàn)對(duì)溫室小氣候整體的變化規(guī)律分析，同時(shí)可進(jìn)行溫室溫度場(chǎng)仿真，為溫室溫度場(chǎng)預(yù)測(cè)與調(diào)控提供新的解決思路[1618]。

　　但由于日光溫室結(jié)構(gòu)復(fù)雜，特別是溫室內(nèi)冠層溫度受作物長(zhǎng)勢(shì)、外界光照、墻體蓄熱、覆蓋材料等因素影響，導(dǎo)致作物冠層溫度的空間分布表現(xiàn)出非線性且變化不定，該方法相對(duì)難以模擬[6,19]，在實(shí)際監(jiān)測(cè)調(diào)控應(yīng)用中存在較大誤差。以上研究雖然在一定程度揭示了日光溫室整體溫度變化，但僅以溫室較少典型位置(溫室中部及兩端)的傳感器部署監(jiān)測(cè)[9,20]，難以全面反映溫室冠層時(shí)空分布不均的實(shí)際情況，導(dǎo)致獲取的特征點(diǎn)存在偏差，無(wú)法為日光溫室生產(chǎn)過(guò)程作出及時(shí)的預(yù)警，防止溫度災(zāi)害的發(fā)生。

　　本文通過(guò)構(gòu)建基于物聯(lián)網(wǎng)技術(shù)的日光溫室溫度場(chǎng)監(jiān)測(cè)系統(tǒng)，采集冠層溫度數(shù)據(jù)，采用插值算法完成番茄冠層溫度場(chǎng)的擬合，并進(jìn)行冠層溫度場(chǎng)特征分析，提取溫度極值點(diǎn)，研究不同天氣條件下溫度極值點(diǎn)分布規(guī)律，以期為日光溫室設(shè)施調(diào)控傳感器及增溫設(shè)施高效部署提供理論依據(jù)。

　　1試驗(yàn)數(shù)據(jù)獲取與處理

　　1.1試驗(yàn)對(duì)象

　　試驗(yàn)采用典型西北下沉式日光溫室為試驗(yàn)對(duì)象，屬于單棟、單坡面、厚土墻式溫室，其坐北朝南，溫室東西長(zhǎng)50m，南北跨度7m，脊高5m，下沉0.5，后墻、東、西兩側(cè)墻體為黏土制成，拱架為鋼架結(jié)構(gòu)，后坡覆以PC板和紡織材料制成的保溫棉被，透光材料為聚氯乙烯薄膜，采用上下通風(fēng)方式建造。試驗(yàn)于陜西省涇陽(yáng)縣西北農(nóng)林科技大學(xué)涇陽(yáng)蔬菜示范基地號(hào)日光溫室進(jìn)行，位于北緯34.4°~34.6°，東經(jīng)108.4°~108.9°，屬于典型的溫帶大陸性季風(fēng)氣候。棚內(nèi)種植作物為番茄，于2019年10月20日定植，溫室內(nèi)壟寬0.8m，壟間距0.8m，每壟種植行番茄，共計(jì)33壟，整個(gè)生育期采用滴灌方式進(jìn)行灌溉，試驗(yàn)期間溫室于8:009:00間揭保溫被，晴天多云天氣10:0010:30間開上通風(fēng)進(jìn)行降溫及外界二氧化碳補(bǔ)充，陰(雨)天一般不通風(fēng)或延遲至12：30后。

　　1.2建模數(shù)據(jù)獲取

　　為采集番茄冠層溫度分布情況數(shù)據(jù)，采用自主開發(fā)的溫度場(chǎng)監(jiān)測(cè)系統(tǒng)，包括溫度采集節(jié)點(diǎn)、485節(jié)點(diǎn)、TU數(shù)據(jù)上傳模塊與物聯(lián)網(wǎng)監(jiān)控平臺(tái)。整套溫度場(chǎng)監(jiān)測(cè)系統(tǒng)分路共計(jì)40個(gè)溫度采集節(jié)點(diǎn)。

　　溫度采集節(jié)點(diǎn)以MSP430FR2433微處理器為核心處理單元，采用德國(guó)賀利氏生產(chǎn)的PT100型薄膜鉑熱電阻進(jìn)行溫度采集，其使用不銹鋼保護(hù)管與四氟高密度鍍銀屏蔽線進(jìn)行封裝，具有良好的抗腐蝕性與抗干擾性能，測(cè)量精度0.15℃，溫度測(cè)量范圍50℃300℃;溫度信號(hào)通過(guò)ADS127424位高精度采樣芯片進(jìn)行轉(zhuǎn)換，并由串口傳送給核心處理單元，每個(gè)處理單元具有單獨(dú)編號(hào)，各個(gè)單元間通過(guò)485總線形式進(jìn)行通信。485節(jié)點(diǎn)采用STM32F103微處理器，負(fù)責(zé)接收處理每個(gè)溫度采集節(jié)點(diǎn)傳輸數(shù)據(jù)。數(shù)據(jù)上傳模塊采用Comway4GDTU全網(wǎng)通透?jìng)髂�塊(模塊)，負(fù)責(zé)將匯集的溫度場(chǎng)數(shù)據(jù)上傳。

　　最終數(shù)據(jù)匯聚至物聯(lián)網(wǎng)監(jiān)控平臺(tái)進(jìn)行管理與分析。試驗(yàn)于2019年11月27日12月日進(jìn)行，番茄處于花期，試驗(yàn)期間植株平均高度4060cm，共連續(xù)采集冠層溫度場(chǎng)數(shù)據(jù)，其中晴天、多云、陰(雨)天天數(shù)分別為 。溫度場(chǎng)監(jiān)測(cè)系統(tǒng)每5min記錄一次溫度場(chǎng)數(shù)據(jù)，通過(guò)4G模塊將數(shù)據(jù)上傳至監(jiān)控平臺(tái)。溫度場(chǎng)監(jiān)測(cè)系統(tǒng)40個(gè)溫度傳感器按圖中所示方案部署，共列，東西距側(cè)墻各4m作為緩沖區(qū)，各列相距6m;每列部署個(gè)傳感器分別距后墻1.5、2.5、3.5、4.5、5.5m，距地面0.5m。本文采用試驗(yàn)期間19數(shù)據(jù)進(jìn)行分析，每30min選取一次數(shù)據(jù)，共獲912組冠層原始溫度數(shù)據(jù)，每組40個(gè)溫度點(diǎn)，總計(jì)36480個(gè)溫度點(diǎn)數(shù)據(jù)用于番茄冠層溫度場(chǎng)插值模型構(gòu)建。

　　1.3數(shù)據(jù)分析

　　以實(shí)測(cè)試驗(yàn)數(shù)據(jù)為基礎(chǔ)，通過(guò)克里金法(Kriging)計(jì)算溫度場(chǎng)數(shù)據(jù)，進(jìn)而采用差分進(jìn)化算法獲取番茄冠層溫度場(chǎng)極值特征點(diǎn)，探尋日光溫室作物冠層溫度特征點(diǎn)分布及變化規(guī)律。

　　1.3.1克里金插值

　　針對(duì)冠層溫度場(chǎng)連續(xù)且分布不均、差異性大、不規(guī)律，以及試驗(yàn)過(guò)程中因試驗(yàn)條件限制導(dǎo)致支撐建模分析數(shù)據(jù)量少的問(wèn)題，本文采用簡(jiǎn)單克里金算法進(jìn)行溫度場(chǎng)建模，克里金法是通過(guò)協(xié)方差函數(shù)對(duì)隨機(jī)過(guò)程或場(chǎng)進(jìn)行空間建模和預(yù)測(cè)的回歸算法，能實(shí)現(xiàn)目標(biāo)位置的無(wú)偏估計(jì)[21,22]，其廣泛應(yīng)用于農(nóng)業(yè)、氣象等空間領(lǐng)域插值，如土壤有機(jī)質(zhì)空間分布、作物需水量變化趨勢(shì)分析、地下地表溫度預(yù)測(cè)、降水時(shí)空分布等[2327]。

　　本文通過(guò)Pycharm軟件與克里金算法庫(kù)進(jìn)行溫度場(chǎng)建模，采用留一法對(duì)19每個(gè)時(shí)刻監(jiān)測(cè)溫度數(shù)據(jù)進(jìn)行交叉驗(yàn)證的數(shù)據(jù)集劃分，即每個(gè)時(shí)刻溫度場(chǎng)的40個(gè)溫度點(diǎn)中每個(gè)點(diǎn)輪流作為測(cè)試集，剩余39個(gè)點(diǎn)為訓(xùn)練集，循環(huán)建模驗(yàn)證40次，每個(gè)時(shí)刻共獲取40個(gè)驗(yàn)證模型及驗(yàn)證數(shù)據(jù)點(diǎn)，19共計(jì)獲取36480個(gè)驗(yàn)證數(shù)據(jù)用于插值模型性能分析。其中克里金算法采用粒子群算法作為激活函數(shù)，由于監(jiān)測(cè)數(shù)據(jù)量有限，預(yù)測(cè)點(diǎn)相鄰要素點(diǎn)數(shù)量設(shè)置為10，即預(yù)測(cè)點(diǎn)數(shù)據(jù)由其最近臨10個(gè)監(jiān)測(cè)點(diǎn)計(jì)算所得。

　　2冠層溫度特征點(diǎn)日變化規(guī)律分析

　　2.1冠層溫度場(chǎng)特征點(diǎn)獲取以克里金插值模型為差分進(jìn)化算法的適應(yīng)度函數(shù)，設(shè)置算法參數(shù)種群規(guī)模為20，迭代次數(shù)1000，縮放因子0.5，交叉因子0.8，對(duì)日光溫室冠層溫度進(jìn)行二維尋優(yōu)，以其中2019年12月03日00時(shí)溫度最大值尋優(yōu)為例其進(jìn)化曲線可知采用DE算法收斂速度較快，在進(jìn)行至11次迭代時(shí)基本平穩(wěn)，能快速獲取溫度極值。本文共獲取連續(xù)19每30min一組極值數(shù)據(jù)，共獲取1824個(gè)溫度極值特征點(diǎn)與位置坐標(biāo)，其中最高溫與最低溫點(diǎn)各912個(gè)，晴天、多云天極值各336組，陰(雨)天極值240組。

　　2.2冠層溫度特征點(diǎn)日變化規(guī)律分析通過(guò)建模尋優(yōu)結(jié)果分析不同天氣條件下溫室內(nèi)冠層溫度特征點(diǎn)的分布及變化規(guī)律，其中溫度變化趨勢(shì)為不同天氣條件下，同一時(shí)刻所有天溫度均值的變化趨勢(shì)。

　　農(nóng)藝師論文投稿刊物：《農(nóng)業(yè)工程學(xué)報(bào)》是由中國(guó)農(nóng)業(yè)工程學(xué)會(huì)主辦的全國(guó)性學(xué)術(shù)期刊，自2005年始為單月刊。主要欄目有：農(nóng)業(yè)水土工程，農(nóng)業(yè)裝備工程與機(jī)械化，農(nóng)業(yè)信息與電氣技術(shù)，農(nóng)業(yè)生物環(huán)境與能源工程，土地整理工程，農(nóng)產(chǎn)品加工工程。

　　3結(jié)論

　　(1)克里金插值能較好獲取符合實(shí)際冠層溫度的插值數(shù)據(jù)，經(jīng)驗(yàn)證其插值后晴天、多云天氣、陰(雨)天的決定系數(shù)分別為0.9553、0.9604、0.9447，均方根誤差分別為1.34℃、0.95℃、0.40℃，由于日光輻射變化及通風(fēng)等農(nóng)事操作影響，克里金插值更適用于陰(雨)天全天及夜間低溫冠層溫度場(chǎng)的數(shù)據(jù)插值。(2)不同天氣條件下，日光溫室內(nèi)冠層溫度變化趨勢(shì)基本一致，揭被后溫室內(nèi)溫度呈上升趨勢(shì)，最低溫主要出現(xiàn)在揭被及日出前后，最高溫出現(xiàn)在中午13:0分左右，夜間溫度下降速率為晴天最快，多云天氣次之，陰(雨)天最慢，最低溫與下降速率規(guī)律相反為陰(雨)天最高;(3)不同天氣條件下溫室內(nèi)冠層最高溫點(diǎn)存在差異，但主要集中出現(xiàn)在溫室中部，即橫向22.0、距后墻2.附近位置，最低溫點(diǎn)主要集中位在溫室橫向.0、距后墻5.48靠近東墻與薄膜附近區(qū)域，分析所得特征點(diǎn)分布位置能為溫室溫度災(zāi)害監(jiān)測(cè)、溫室智能化控制、增溫設(shè)施部署等提供有效參考。

　　參考文獻(xiàn)

　　[1]張宏鳴,王佳佳,韓文霆,等基于熱紅外遙感影像的作物冠層溫度提取[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),201950(4):203210.ZHANGHongming,WANGJiajia,HANWenting,etal.Cropcanopytemperatureextractionbasedonthermalinfraredremotesensingimages[J].TransactionsoftheChineseSocietyforAgriculturalMachinery,201950(4):203210.http://www.jcsam.org/jcsam/ch/reader/create_pdf.aspx?file_no=20190423&flag=1&journal_id=jcsam.DOI:10.6041/.issn.10001298.2019.04.023.(inChinese)

　　[2]張勇,高文波,鄒志榮日光溫室主動(dòng)蓄熱后墻傳熱CFD模擬及性能試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),201531(5):203211.ZHANGYong,GAOWenboandZOUZhirongPerformanceexperimentandCFDsimulationofheatexchangeinsolargreenhousewithactivethermalstoragebackwall[J]TransactionsoftheCSAE,2015,35(5):203211.(inChinese)

　　[3]鮑恩財(cái),曹晏飛,鄒志榮,等不同結(jié)構(gòu)主動(dòng)蓄熱墻體日光溫室傳熱特性[J].農(nóng)業(yè)工程學(xué)報(bào),201935(3):189197.BAOEncai,CAOYanfei,ZOUZhirong,etal.Characteristicofheattransferforactiveheatstoragewallwithdifferentstructuresinhinesesolargreenhouse[J].TransactionsoftheCSAE,201935(3):189197.(inChinese)

　　作者：張軍華1,，沈楷程1,，陳丹艷1,，張明科，張海輝1,，胡瑾1,

轉(zhuǎn)載請(qǐng)注明來(lái)自發(fā)表學(xué)術(shù)論文網(wǎng)：http:///nylw/26938.html