玉米糧堆霉變發(fā)熱過程中的溫濕度場變化規(guī)律研究

王小萌，吳文福，尹 君，張忠杰，吳子丹※，姚 渠

?

玉米糧堆霉變發(fā)熱過程中的溫濕度場變化規(guī)律研究

王小萌1,2，吳文福1，尹 君2，張忠杰2，吳子丹1※，姚 渠2

（1. 吉林大學(xué)生物與農(nóng)業(yè)工程學(xué)院，長春 130022；2. 國家糧食和物資儲備局科學(xué)研究院，北京 100037）

為模擬儲糧糧堆局部含水率偏高引起的霉變發(fā)熱現(xiàn)象，進而研究此現(xiàn)象中溫、濕度場的變化規(guī)律，該文在試驗倉內(nèi)濕基含水率14.0%的玉米糧堆中心加入濕基含水率18.2%的玉米，在30 ℃室內(nèi)儲藏40 d。試驗糧堆由于霉變引起自發(fā)熱。試驗過程中，通過計算玉米糧堆中垂面內(nèi)高溫區(qū)和高濕區(qū)的面積變化，從而揭示玉米糧堆霉變發(fā)熱過程中溫、濕度場的變化規(guī)律。試驗結(jié)果表明，糧堆中垂面高濕區(qū)面積緩慢擴大，高溫區(qū)面積開始擴大緩慢，但在與周圍糧溫最高溫差升至3.7 ℃后，面積擴大速率加快，且高溫區(qū)與高濕區(qū)面積的當(dāng)量半徑與溫度差D成正比，此正比關(guān)系經(jīng)過了糧庫淺圓倉的驗證。這為進一步定量分析糧食倉儲過程中的高溫區(qū)和高濕區(qū)擴散提供了依據(jù)。

作物；溫度；濕度傳感器；玉米；溫度場；相對濕度場；霉變發(fā)熱

0 引 言

中國玉米種植范圍廣，產(chǎn)量高，是重要的糧食及經(jīng)濟作物之一。在儲藏過程中玉米糧堆霉變發(fā)熱一直是比較普遍的問題。

對于儲糧糧堆發(fā)熱的起因，國內(nèi)外已經(jīng)做了大量的研究。微生物作用，害蟲大量滋生，糧食呼吸作用，雜質(zhì)聚集等均可以引發(fā)糧堆發(fā)熱[1]。其中，蟲[2]、霉[3]是引起糧堆發(fā)熱的重要因素[4]，但是如果昆蟲密度不高，僅依靠昆蟲自身呼吸作用引發(fā)糧堆發(fā)熱的幾率比高水分糧食霉變引發(fā)糧堆發(fā)熱的幾率低[5-6]。儲糧真菌呼吸作用是糧堆發(fā)熱的主要熱源，且蟲害一般會發(fā)生在霉變發(fā)熱之后[7]。

糧食儲藏過程中影響儲糧的主要霉菌是曲霉和青霉[8]，而影響霉菌生長的重要因素是溫度和水分[9]，因此，儲糧霉變發(fā)熱過程中，溫度、水分與霉變的關(guān)系是需要密切關(guān)注的問題。針對這個問題，國內(nèi)外學(xué)者做了大量的研究。Koehler[10]研究了含水率對玉米霉變的影響。Jian等[11]研究了不同溫度和含水率下蟲害和微生物的演替規(guī)律。Huang等[12]研究了不同溫度和含水率的玉米糧堆CO2產(chǎn)生速率。張航等[13-16]分別研究了中國小麥、玉米、水稻和大豆霉變與溫度和水分的關(guān)系。

在儲糧霉變發(fā)熱過程中，王小萌等[17]對霉變分布與糧堆溫度場和濕度場之間的關(guān)系進行了分析。Wallace等[18]研究了發(fā)熱糧堆的水分遷移，并發(fā)現(xiàn)糧堆發(fā)熱會引起發(fā)熱點上層糧堆含水率增加，下層含水率降低[19]。張中濤[20]利用有限元法研究了小麥呼吸作用對糧堆溫度場、水分場的影響。除此之外，國內(nèi)外學(xué)者也對儲糧發(fā)熱點的溫度變化進行了研究。Zhang等[21]研究了2種高水分小麥的發(fā)熱過程，并建立了發(fā)熱量與時間的指數(shù)模型。Wilson[22]建立了糧堆霉變發(fā)熱過程中熱產(chǎn)生速率和CO2產(chǎn)生速率的方程，但是方程不適用于低水分糧食。Jia等[23-24]分別利用有限元和離散元建立了糧堆內(nèi)存在發(fā)熱點時的溫度分布模型。但是，目前對糧堆霉變發(fā)熱過程中，溫、濕度場在空間上的變化規(guī)律研究甚少。

本文利用儲糧霉變自發(fā)熱引發(fā)糧堆內(nèi)溫度傳遞和水分遷移。通過計算玉米糧堆內(nèi)高溫區(qū)和高濕區(qū)的面積，研究糧堆霉變發(fā)熱過程中溫、濕度場在空間上的變化規(guī)律，以期對糧食倉儲過程中的糧情調(diào)控做出指導(dǎo)。

1 材料與方法

1.1 試驗樣品

本試驗采用的“先玉”玉米樣品初始平均濕基含水率為23%（下文提及的含水率均為濕基含水率），質(zhì)量密度為698 g/L，雜質(zhì)質(zhì)量分?jǐn)?shù)為1.0%，樣品初始真菌孢子數(shù)<104個/g，無儲藏真菌。

1.2 試劑及儀器

多參數(shù)糧情檢測系統(tǒng)[17]：中國深圳市東大恒豐科技有限公司；DHG—9140型電熱恒溫鼓風(fēng)干燥箱：上海精宏實驗設(shè)備有限公司；NIKON E100顯微鏡：日本尼康公司。去離子水作為試驗用水。

1.3 試驗方法

1.3.1 試驗倉裝置

為模擬淺圓倉中糧堆局部含水率偏高引起的局部霉變發(fā)熱現(xiàn)象，依據(jù)相似性設(shè)計原理和淺圓倉倉型，選用圓形的試驗倉[17](圖1)，并進行改進，去除倉底的保溫隔熱材料，將倉底中心位置開1個內(nèi)直徑為0.08 m的通氣孔，目的是有利于試驗倉內(nèi)外氣體進行交換，確保試驗倉內(nèi)氧氣供應(yīng)充足。

1.3.2 樣品處理

將平均含水率為23%的玉米進行除雜、清理，玉米樣品鋪成厚度為0.02 m的薄層，在室外通風(fēng)環(huán)境下，將含水率降低。玉米入倉儲藏安全水分為14%左右[25]。儲糧實踐發(fā)現(xiàn)由于儲糧水分分布不均、雜質(zhì)聚集、水分遷移等因素的影響，玉米糧堆局部含水率可增至16～18%。糧堆內(nèi)的熱傳遞和水分遷移會受到糧堆高度和跨度等的影響[26]，鑒于試驗倉底面直徑和高度較小，為便于觀察糧堆霉變發(fā)熱后的溫濕度遷移現(xiàn)象，因此試驗用糧的高、低含水率可分別設(shè)為18%左右和14%左右。

1.3.3 儲藏模擬設(shè)置

試驗倉內(nèi)玉米糧堆總高度為0.60 m，高水分玉米（18.2%）和低水分玉米（14.0%）的分布如圖1a所示，2種含水率的糧堆體積比為1∶4，將試驗倉儲藏在溫度為30 ℃的室內(nèi)，以期模擬儲糧過程中由于糧堆局部含水率偏高引起的糧堆霉變發(fā)熱現(xiàn)象。2種含水率的玉米的初始糧溫為22 ℃，初始平均相對濕度分別為67%和90%。

1.4 測量指標(biāo)及方法

1.4.1 危害真菌孢子

參考儲糧真菌檢測標(biāo)準(zhǔn)[27]，檢測樣品的真菌孢子數(shù)。

1.4.2 溫濕度

利用溫濕度傳感器實時監(jiān)測玉米糧堆內(nèi)溫濕度變化。為了充分重現(xiàn)溫濕度場在2種含水率玉米中的分布，在借鑒文獻[17]溫濕度傳感器布置的基礎(chǔ)上，增加含水率14.0%玉米中的傳感器個數(shù)。糧堆內(nèi)溫濕度傳感器布置如圖1。

1.4.3 含水率

參考ASAE標(biāo)準(zhǔn)[28]，稱取15.0 g玉米籽粒樣品，放在干燥鋁盒內(nèi)，在103 ℃烘箱內(nèi)干燥72 h，檢測樣品的含水率。

1.5 數(shù)據(jù)處理

1.5.1 溫濕度場的建立

選取糧堆中垂面溫濕度傳感器檢測到的溫度和相對濕度數(shù)據(jù)，再運用溫濕度擬合算法得到糧堆內(nèi)各點的溫度和相對濕度，通過Matlab 模擬軟件重現(xiàn)糧堆的溫度場和濕度場，溫濕度場擬合算法如式（1）～（2）[29]。

由數(shù)值計算方法可知，已知個離散數(shù)據(jù)(x,(x))=1,2,3,…，在[,]上滿足x?[,]，且(x)=(x)，則

注：“●”表示糧堆溫度和相對濕度的監(jiān)測位置，O點為含水率18.2%的玉米糧堆中心。

1.5.2 當(dāng)量半徑的計算

利用Matlab軟件對糧堆中垂面的溫、濕度場云圖進行高溫區(qū)和高濕區(qū)面積求解，主要步驟是：導(dǎo)入云圖圖片，灰度處理，二值化處理，計算高溫區(qū)和高濕區(qū)面積。本文中，將高溫區(qū)和高濕區(qū)的面積按照圓形處理，計算出圓的半徑，作為高溫區(qū)和高濕區(qū)的當(dāng)量半徑。

2 結(jié)果與分析

2.1 溫度和相對濕度變化

圖2和圖3分別是玉米糧堆儲藏0、20和40 d時糧堆中垂面的溫度場和相對濕度場云圖。由圖2可知，儲藏溫度為30 ℃，平均含水率為18.2%的玉米糧堆中由于微生物作用逐漸出現(xiàn)了發(fā)熱點。微生物迅速生長，發(fā)熱點面積不斷擴大，發(fā)熱點最高溫度和周圍糧堆溫度的溫差不斷升高，最大溫差達到8 ℃。研究表明，糧堆內(nèi)的微氣流運動是引起糧堆內(nèi)溫度傳遞和水分遷移的主要因素[30]。從圖3可以看出，糧堆內(nèi)微氣流運動引起糧堆內(nèi)水分向糧堆上層遷移。試驗結(jié)束后，取樣檢測真菌孢子數(shù)，點、點、點、點、點的真菌孢子數(shù)分別為6.3′106，6.9′106，2.1′106，5.4′106，5.7′106個/g，優(yōu)勢真菌為亮白曲霉和黃曲霉。

注：O、A、B,、C、D在溫度場云圖中是共用的,下同。

圖3 玉米糧堆儲藏0、20和40 d時中垂面相對濕度場云圖

2.2 高溫區(qū)與高濕區(qū)面積變化

糧堆內(nèi)發(fā)熱點周圍平均糧溫為30 ℃，因此本文中的高溫區(qū)選取玉米糧堆中垂面溫度高于30 ℃的區(qū)域。亮白曲霉和黃曲霉是引起本文糧堆發(fā)熱的主要霉菌，而相對濕度高于75%時易于亮白曲霉和黃曲霉生長，因此，本文中的高濕區(qū)選取玉米糧堆中垂面相對濕度高于75%的區(qū)域。圖4是試驗過程中，玉米糧堆中垂面、、、、5點的溫度變化圖。圖5是玉米糧堆中垂面的高溫區(qū)和高濕區(qū)面積隨時間的變化圖。由圖5可知，在試驗的40 d內(nèi)，高濕區(qū)面積增加緩慢；高溫區(qū)面積在試驗第25 d后迅速增加，結(jié)合圖4分析，可發(fā)現(xiàn)試驗第25 d后糧溫升高加快，導(dǎo)致高溫區(qū)面積迅速擴大，這主要是因為黃曲霉和亮白曲霉大面積生長，產(chǎn)生大量的熱。

圖4 玉米糧堆中垂面不同點溫度變化

參考當(dāng)量直徑的定義，本文中將高溫區(qū)和高濕區(qū)的面積按照圓形處理，計算出圓的半徑，作為高溫區(qū)和高濕區(qū)的當(dāng)量半徑。溫度差是糧堆內(nèi)濕熱轉(zhuǎn)移的動力源[31-32]。因此，建立玉米糧堆高溫區(qū)和高濕區(qū)當(dāng)量半徑與溫度差D的關(guān)系（圖6）。由圖6可知，高濕區(qū)面積的當(dāng)量半徑與溫度差成正比；當(dāng)溫度差高于3.7 ℃即高溫區(qū)最高溫度高于33.7 ℃后，發(fā)熱區(qū)開始快速向周圍擴散，這與糧堆的導(dǎo)熱系數(shù)小有很大關(guān)系，糧堆導(dǎo)熱系數(shù)小，導(dǎo)致糧堆霉變產(chǎn)生的熱量難以快速傳遞出糧堆，很容易聚集在糧堆內(nèi)；高溫區(qū)面積的當(dāng)量半徑與溫度差成線性正相關(guān)，且D/D=0.058 (2=0.99)。儲糧霉變發(fā)熱過程中，發(fā)熱區(qū)域的擴散速率會受到熱傳導(dǎo)和熱對流的共同作用，由于試驗倉體積小，高度低，熱對流現(xiàn)象微弱，因此糧堆內(nèi)的熱量傳遞和水分遷移主要受到熱傳導(dǎo)的影響。下面用實倉對高溫區(qū)當(dāng)量半徑與溫度差的關(guān)系進行驗證。

圖5 高溫區(qū)和高濕區(qū)面積隨儲藏時間的變化

圖6 高溫區(qū)和高濕區(qū)面積的當(dāng)量半徑與溫度差的關(guān)系

2.3 實倉驗證

以中國某糧庫的淺圓倉糧溫數(shù)據(jù)驗證上述正比關(guān)系。此淺圓倉為鋼屋蓋，倉內(nèi)直徑為25.0 m，裝糧高度為15.6 m，雜質(zhì)質(zhì)量分?jǐn)?shù)為1.5%，糧堆質(zhì)量密度為698 g/L，儲量為5 756.9 t。所儲玉米的平均含水率為13.8%，糧堆局部偏高水分區(qū)域的平均含水量為14.8%。糧堆內(nèi)局部含水率偏高，導(dǎo)致出現(xiàn)發(fā)熱點。通過計算糧堆中垂面高溫區(qū)的面積變化，得到高溫區(qū)當(dāng)量半徑與溫度差的關(guān)系，如圖7。

圖7 淺圓倉中垂面高溫區(qū)的當(dāng)量半徑與溫度差的關(guān)系

在之前研究[17]中，利用外加熱源引起小麥糧堆內(nèi)的溫濕度遷移，試驗過程中沒有出現(xiàn)由于糧堆霉變引發(fā)的發(fā)熱點，以此作者研究了儲糧真菌分布和溫濕度的關(guān)系。此文中，利用高水分玉米自身霉變出現(xiàn)發(fā)熱點，模擬儲糧糧堆局部含水率偏高引起的糧堆霉變發(fā)熱，發(fā)現(xiàn)糧堆高溫區(qū)當(dāng)量半徑與溫度差的正比例關(guān)系，鑒于此文中沒有找到高溫區(qū)當(dāng)量半徑與溫度差的具體定量關(guān)系，還需進行更深入的研究。

3 結(jié) 論

模擬試驗倉中，由18.2%和14.0% 2種含水率的玉米組成的糧堆在30 ℃恒溫室內(nèi)儲藏，4 d后，亮白曲霉和黃曲霉大量生長引發(fā)糧堆出現(xiàn)明顯發(fā)熱點。

玉米糧堆霉變發(fā)熱過程中，主要受熱傳導(dǎo)影響，高溫區(qū)與高濕區(qū)面積不斷擴大。在高溫區(qū)與周圍糧溫最大溫差從3.7升高至8 ℃的過程中，高溫區(qū)和高濕區(qū)面積的當(dāng)量半徑與溫度差呈正相關(guān)。該相關(guān)性也經(jīng)實倉數(shù)據(jù)得到了驗證。

試驗倉和淺圓倉高溫區(qū)與周圍糧溫的溫度差均升高3.7 ℃后，高溫區(qū)面積才開始迅速擴大，但是淺圓倉熱對流強烈，導(dǎo)致淺圓倉高溫區(qū)擴散速率遠大于試驗倉中高溫區(qū)擴散速率。

[1] 房強. 新糧入倉后發(fā)熱的預(yù)防和處理措施探討[J]. 糧食儲藏，2016，45(3)：53－56. Fang Qiang. Discussion on prevention and treatment measures of new grain storage heating[J]. Grain Storage, 2016, 45(3): 53－56. (in Chinese with English abstract)

[2] Sinha R N. Insects and mites associated with hot spots in farm stored grain[J]. Canadian Entomologist, 1961, 93 (8): 609－621.

[3] Wallace H A H, Sinha R N. Fungi associated with hot spots in farm stored grain[J]. Canadian Journal of Plant Science, 1962, 42(1): 130－141.

[4] Sinha R A, Wallace H A H. Ecology of a fungus-induced hotspot in stored grain[J]. Canadian Journal of Plant Science, 1965, 45(45): 48－59.

[5] Cofie-Agblor R, Muir W E, Sinha R N. Comparative heat of respiration of five grain beetles in stored wheat[J]. Postharvest Biology and Technology, 1995, 5 (1/2): 167－175.

[6] Mani S. Computer modelling of insect-induced hot spots in stored wheat[D]. Manitoba: University of Manitoba, 2001.

[7] Francis F L. Integrated management of the risks of stored grain spoilage by seedborne fungi and contamination by storage mould mycotoxins: An update[J]. Journal of Stored Products Research, 2017, 71: 22－40.

[8] Pronyk C, Muir W E, Ndg W，et al. Carbon dioxide production and deterioration of stored canola[J]. Canadian Biosystems Engineering, 2004, 46(3): 25－33.

[9] Mohapatra D, Kumar S, Kotwaliwale N, et al. Critical factors responsible for fungi growth in stored food grains and non-chemical approaches for their control[J]. Industrial Crops and Products, 2017, 108 : 162－182.

[10] Koehler B. Fungus growth in shelled corn as affected by moisture[J]. Journal of Agricultural Research, 1938, 56(4): 291－307.

[11] Jian F, Jayas D S. The ecosystem approach to grain storage[J]. Agricultural Research, 2012, 1 (2): 148－156.

[12] Huang H, Danao M G C, Rausch K D, et al. Diffusion and production of carbon dioxide in bulk corn at various temperatures and moisture contents[J]. Journal of Stored Products Research, 2013, 55 (4) : 21－26.

[13] 張航，黃淑霞，蔡靜平，等. 小麥臨界水分附近微生物活動的特點[J]. 河南工業(yè)大學(xué)學(xué)報(自然科學(xué)版)，2010，31(5)：30－33. Zhang Hang, Huang Shuxia, Cai Jingping, et al. Characteristics of microbial activitity near critical moisture content of wheat[J]. Journal of Henan University of Technology(Natural Science Edition), 2010, 31(5): 30－33. (in Chinese with English abstract)

[14] 唐芳，程樹峰，張海洋，等. 稻谷儲藏真菌危害早期預(yù)測的研究[J]. 糧食儲藏，2015，44(1)：24－27. Tang Fang, Cheng Shufeng, Zhang Haiyang, et al. Study on the forecasting fungal growth in stored paddy[J]. Grain Storage, 2015, 44 (1): 24－27. (in Chinese with English abstract)

[15] 程樹峰，唐芳，歐陽毅，等. 大豆儲藏真菌危害早期預(yù)測的研究[J]. 中國糧油學(xué)報，2015，30(7)：77－80. Cheng Shufeng, Tang Fang, Ouyang Yi, et al. The early prediction of damage of fungus during soybean storage[J]. Journal of the Chinese Cereals and Oils Association, 2015, 30(7): 77－80. (in Chinese with English abstract)

[16] 歐陽毅，祁智慧，李春元，等. 玉米儲藏真菌早期預(yù)測的研究[J]. 糧油食品科技，2017，25(5)：52－55. Ouyang Yi, Qi Zhihui, Li Yuanchun, et al. Early prediction of fungus hazard during corn storage [J]. Science and Technology of Cereals, Oils and Foods, 2017, 25(5): 52－55. (in Chinese with English abstract)

[17] 王小萌，吳文福，尹君，等. 基于溫濕度場云圖的小麥糧堆霉變與溫濕度耦合分析[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(10)：260－266. Wang Xiaomeng, Wu Wenfu, Yin Jun, et al. Analysis of wheat bulk mould and temperature-humidity coupling based on temperature and humidity field cloud map[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(10): 260－266. (in Chinese with English abstract)

[18] Wallace H A H, Sinha R N, Mills T. Fungi associated with small wheat bulks during prolonged storage in Manitoba[J]. Canadian Journal of Botany, 1976, 54(12): 1332－1343.

[19] Sinha R N，Wallace H A H. Ecology of insect-induced hot spots in stored grain in western canada[J]. Researches on Population Ecology, 1966, 8(2): 107－132.

[20] 張中濤. 具有吸濕和呼吸特性的倉儲糧堆內(nèi)熱濕耦合傳遞規(guī)律的研究[D]. 濟南：山東建筑大學(xué)，2015. Zhang Zhongtao. Research of Coupled Heat and Moisture Transfer in Hygroscopic and Respiration Bulk During the Warehouse Storage[D]. Jinan: Shandong Jianzhu University, 2015. (in Chinese with English abstract)

[21] Zhang Q, Muir W E, Sinha R N, et al. Heat production in wet wheat under adiabatic conditions[J]. Canadian Agricultural Engineering, 1992, 34(3): 233－238.

[22] Wilson M T. A model for predicting mould growth and subsequent heat generation in bulk stored grain[J]. Journal of Stored Products Research, 1999, 35(1): 1－13.

[23] Jia Canchun, Sun Dawen, Cao Chongwen. Mathematical simulation of temperature fields in a stored grain bin due to internal heat generation[J]. Journal of Stored Products Research，2000, 43: 227－233.

[24] Rusinek R, Kobylka R. Experimental study and discrete element method modeling of temperature distributions in rapeseed stored in a model bin[J]. Journal of Stored Products Research, 2014, 59: 254－259.

[25] 李興軍，王雙林，付鵬程. 采用水分解吸等溫線分析糧食安全水分的研究[J]. 糧食加工，2017，42(5)：30－36. Li Xingjun, Wang Shuanglin, Fu Pengcheng. A study on the safe storage moisture of cereals and oils with moisture desorption isotherms[J]. Grain Processing, 2017，42(5): 30－36. (in Chinese with English abstract)

[26] 王遠成，潘鈺，尉堯方，等. 倉儲糧堆內(nèi)部自然對流和熱濕傳遞的數(shù)學(xué)分析及驗證[J]. 中國糧油學(xué)報，2017，32(9)：120－125. Wang Yuancheng, Pan Yu, Yu Yaofang, et al. Mathematical analysis and validation for natural convective heat and moisture transfer in grain bulk during sealed storage[J]. Journal of the Chinese Cereals and Oils Association, 2017, 32(9): 120－125. (in Chinese with English abstract)

[27] 國家糧食局. 糧油檢驗-儲糧真菌檢測-孢子計數(shù)法：LS/T 6132-2018 [S]. 北京：中國標(biāo)準(zhǔn)出版社，2018:03.

[28] ASAE. Moisture measurement-unground grain and seeds: S352. 2 -2012 [S]. Michigan: ASABE, 2012:11.

[29] 尹君，吳子丹，張忠杰，等. 不同倉型的糧堆溫度場重現(xiàn)及對比分析[J]. 農(nóng)業(yè)工程學(xué)報, 2015, 31(1): 281－287. Yin Jun, Wu Zidan, Zhang Zhongjie, et al. Comparison and analysis of temperature field reappearance in stored grain of different warehouses[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(1): 281－287. (in Chinese with English abstract)

[30] 尹君，吳子丹，張忠杰，等. 基于多場耦合理論淺析 淺圓倉局部結(jié)露機理[J]. 中國糧油學(xué)報，2015，30(5)：90－95.

Yin Jun, Wu Zidan, Zhang Zhongjie, et al. Analysis of the formation mechanism of partial condensation for the stored grain in steel squat silo based on multi-fields coupling theory[J].Journal of the Chinese Cereals and Oils Association, 2015，30(5)：90－95. (in Chinese with English abstract)

[31] 闞正武. 溫度梯度作用下室內(nèi)天然氣泄漏擴散濃度場研究[D]. 重慶：重慶交通大學(xué)，2013. Kan Zhengwu. Indoor Diffusion of Natural Gas Concentration Field Under Temperature Gradient[D]. Chongqing: Chongqing Jiaotong University, 2013. (in Chinese with English abstract)

[32] 白忠權(quán). 吸濕性倉儲糧堆內(nèi)熱濕耦合傳遞規(guī)律的研究[D]. 濟南：山東建筑大學(xué)，2013. Bai Zhongquan. Research on Coupled Heat and Moisture Transfer in Hygroscopic Grain Bulk During the Warehouse Storage[D]. Jinan: Shandong Jianzhu University, 2013. (in Chinese with English abstract)

Research on temperature and humidity field change during corn bulk microbiological heating

Wang Xiaomeng1,2, Wu Wenfu1, Yin Jun2, Zhang Zhongjie2, Wu Zidan1※, Yao Qu2

(1.130022,; 2.100037,)

Corn is widely planted in China with high yield. It is a major source of foods and materials for industrial processing.Microbiological heating occurs easily during corn bulk storage, which seriously affects grain safety. The problem has also attracted much attention. Stored grain is mostly infected byand. Fungal growth in maize is facilitated by hot and humid conditions. Warm pockets initiated by microorganisms are mainly induced by pockets of wet grain because microorganisms need larger than 0.65 water activity to multiply and develop. It has been reported by many researchers that the relationship between temperature, moisture and fungi growth. Moisture diffusion and migration from a hotspot to its surrounding was recorded by some researchers. They also found that as a hotspot develops, the grain moisture at the top of the hotspot increases while the grain moisture below decreases. Since it is difficult to estimate the production rate of heat without using mathematical models, there are a few models developed to understand the development of hotspot. These models was used to calculate the heat production rate of stored grain. However, there is little research on the quantitative change of temperature and humidity field in space. The objective of this study was to develop a method to measure changes of high temperature and high humidity zones in space during grain microbiological heating. To explore quantitative variation of temperature and relative humidity fields in space in corn heating, different moisture corn (14.0% and 18.2%, w.b.) was stored in a simulated silo in 40 d at 30 ℃ non-airtight. The simulated silo was a cylindrical iron silo with 0.54 m in internal diameter, 0.70 m in height and 0.01 m in thickness, respectively. And its inner wall was provided with insulation layer (0.02 m thickness). Two air pipes (0.08 m internal diameter) on the top and bottom of the silo were applied to exchange the gas inside and outside the silo and ensure adequate oxygen supply. In the experiment, the high moisture corn (cylinder, diameter, 0.30 m; height, 0.30 m; 18.2%, w.b.) in the silo was surrounded by low moisture corn (14.0%, w.b.). After 4 days storage, the growth ofandcaused a hot spot appears in corn bulk. In the paper, high temperature areas were temperature higher than 30 ℃, and high humidity areas were relative humidity higher than 75% due toandincreasing greatly. During the storage, temperature and relative humidity cloud maps of the min-vertical plane were drawn. These cloud maps indicate that areas of high temperature and high humidity expanded under heat conduction. Areas of high temperature zone and high humidity zone were calculated. Then, treated these areas as circles and calculated equivalent radii () of high temperature zone and high humidity zone. Besides temperature difference (D) were equal to the highest temperature in high temperature zone minus 30 ℃. DuringDincreasing from 3.7 to 8 ℃, equivalent radii () had a significant linear correlation with temperature difference (D). However, no noticeable change was observed whenDranged from 0 to 3.7 ℃. The corn temperature data of a squat silo during microbiological heating proved the linear relationship between the equivalent radius of high temperature area and temperature difference. But the diffusion rate of heating area in squat silo was higher than the simulated silo due to height and span of grain bulk. Height and span of grain bulk in squat siloincreased heat convection which was weak in simulated silo. This study lays a foundation for the further quantitative research on the prediction of microbiological heating in grain storage.

crops; temperature; humidity sensors; corn; temperature field; relative humidity field; microbiological heating

王小萌，吳文福，尹 君，張忠杰，吳子丹，姚 渠.玉米糧堆霉變發(fā)熱過程中的溫濕度場變化規(guī)律研究[J]. 農(nóng)業(yè)工程學(xué)報，2019，35(3)：268－273. doi：10.11975/j.issn.1002-6819.2019.03.033 http://www.tcsae.org

Wang Xiaomeng, Wu Wenfu, Yin Jun, Zhang Zhongjie, Wu Zidan, Yao Qu.Research on temperature and humidity field change during corn bulk microbiological heating[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(3): 268－273. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.03.033 http://www.tcsae.org

2018-12-29

2019-01-16

2015年糧食公益性行業(yè)科研專項項目（201513001）

王小萌，博士生，研究方向為糧食信息化與自動化。 Email：wxmhappy99@163.com。

吳子丹，教授，研究方向為糧食儲藏與運輸。Email：wuzidan@263.net。

10.11975/j.issn.1002-6819.2019.03.033

S379

A

1002-6819(2019)-03-0268-06