陳思羽 吳文福 李興軍 吳子丹 張亞秋 張忠杰 韓 峰

(吉林大學(xué)生物與農(nóng)業(yè)工程學(xué)院1，長春 130022)

(佳木斯大學(xué)機(jī)械工程學(xué)院2，佳木斯 154007)

(國家糧食局科學(xué)研究院3, 北京 100037)

谷物濕熱平衡新模型及熱力學(xué)特性的研究

陳思羽1,2吳文福1李興軍3吳子丹3張亞秋1張忠杰3韓 峰1

(吉林大學(xué)生物與農(nóng)業(yè)工程學(xué)院1，長春 130022)

(佳木斯大學(xué)機(jī)械工程學(xué)院2，佳木斯 154007)

(國家糧食局科學(xué)研究院3, 北京 100037)

建立谷物平衡水分與相對(duì)濕度在不同溫度下的關(guān)系模型，并利用小麥、玉米、水稻、大米的靜態(tài)平衡試驗(yàn)數(shù)據(jù)進(jìn)行擬合。結(jié)果表明，該方程適合描述小麥、玉米、水稻、大米的等溫線，可直接求解出谷物的平衡水分。另外，根據(jù)所建立的模型，預(yù)測了不同溫度條件下的安全儲(chǔ)藏水分，并對(duì)谷物的熱力學(xué)性能進(jìn)行了分析：吸著等熱值隨著平衡水分值的減小而增加，解吸等熱值增加更為顯著；解吸和吸附過程中吸著等熱值均高于水汽化潛熱；平衡水分大于25%后，吸著等熱值趨于穩(wěn)定，與汽化潛熱值接近；根據(jù)分析所得數(shù)據(jù)，擬合了小麥、稻谷、玉米解吸與吸附過程的吸著等熱曲線回歸方程，為以能量化的觀念實(shí)施糧倉作業(yè)管理，實(shí)現(xiàn)高效、合理、節(jié)能儲(chǔ)糧提供了有利依據(jù)。

谷物 吸著等溫線 數(shù)學(xué)模型 熱力學(xué)性質(zhì)

谷物在一定溫度、濕度條件下一段時(shí)間后，水分不再發(fā)生變動(dòng)，此時(shí)水分為谷物的平衡水分。Henderson方程[1-2]，Chung-Pfost方程[3], Halsey方程[4],修正Oswin方程[5], 修正GAB方程[6]等經(jīng)典方程都能預(yù)測不同溫度條件下，水分活度與平衡水分關(guān)系。國內(nèi)外學(xué)者通過經(jīng)典方程對(duì)單一種類谷物的等溫線進(jìn)行擬合[7-8]，但單一種類的等溫線方程和參數(shù)僅適合對(duì)應(yīng)谷物收獲后的干燥和儲(chǔ)藏通風(fēng)操作，不同種類谷物的等溫線則差別較大。我國儲(chǔ)糧機(jī)械通風(fēng)技術(shù)規(guī)程(LS/T 1202—2002)[9]采用WU(CAE)模型[10]對(duì)我國濕熱糧堆進(jìn)行調(diào)控。吳子丹等[11]在分析大米和大豆的WU(CAE)模型時(shí)表明，谷物不同種類之間差異較大，在谷物品種之間差異不明顯，因此在通風(fēng)操作中，應(yīng)根據(jù)不同糧食種類分開使用不同的WU(CAE)模型，但不必考慮同一種類糧食的不同品種的影響。該模型給出了糧食平衡相對(duì)濕度、糧食平衡水分及溫度的函數(shù)關(guān)系，奠定了我國糧倉機(jī)械通風(fēng)的技術(shù)基礎(chǔ)。利用WU(CAE)模型可求解糧食平衡相對(duì)濕度，但不能直接求解糧堆內(nèi)部的平衡水分。由于糧堆內(nèi)部導(dǎo)熱性差，易在糧堆內(nèi)形成溫差，使?jié)駳庠诶浼Z處聚集造成水分遷移，糧食水分發(fā)生變化。單靠入倉前檢測的糧食平衡水分不能更好地反映濕熱糧堆內(nèi)部的濕熱交換關(guān)系。

本試驗(yàn)建立了不同種類谷物在不同溫度(10、20、25、30、35 ℃)下糧倉內(nèi)部點(diǎn)糧食平衡水分與溫度和相對(duì)濕度之間的關(guān)系模型，可直接求解糧倉內(nèi)部點(diǎn)的平衡水分，并預(yù)測谷物的安全儲(chǔ)藏水分，分析谷物水分解吸、吸附特性及熱力學(xué)性質(zhì)，為以能量化的觀念實(shí)施糧倉作業(yè)管理，實(shí)現(xiàn)高效、合理、節(jié)能儲(chǔ)糧提供了有利依據(jù)。

1 材料與方法

物料含濕量和空氣濕度之間的平衡關(guān)系，可通過恒定溫度下使物料與空氣經(jīng)過足夠長時(shí)間相接觸的試驗(yàn)來測定。本研究所使用的試驗(yàn)數(shù)據(jù)為文獻(xiàn)[7-8、11-13]中采用靜態(tài)稱重法測定樣品在不同溫、濕度條件下解吸和吸附的平衡水分?jǐn)?shù)據(jù)。

2 谷物EMC/ERH模型建立及參數(shù)擬合

2.1 谷物EMC/ERH模型建立

建立糧食平衡水分、溫度和相對(duì)濕度關(guān)系模型為：

ERH=f(EMC,t)形式

ERH=1/[1+exp(A+B·EMC+Ct)]

(1)

EMC=f(ERH,t)形式

EMC=[ln(1-ERH)-lnERH-A-Ct]/B

(2)

式中：ERH為谷物的相對(duì)濕度(小數(shù)表示)；EMC為平衡水分(小數(shù)表示)；A、B、C為不同糧種解吸和吸附不同狀態(tài)下的擬合系數(shù)；t為糧食溫度/℃。

2.2 谷物EMC/ERH模型擬合結(jié)果分析

表1擬合結(jié)果表明在相對(duì)濕度方程的擬合中，小麥、玉米、水稻、大米的R2均大于0.99，MRE值小于8%，擬合結(jié)果較好。表2擬合結(jié)果表明在平衡水分方程的擬合中，小麥、玉米、水稻、大米的R2均大于97%，MRE值均小于7%，擬合度較好。

表1 ERH=f (EMC ,t)形式模型參數(shù)值

表2 EMC=f (ERH,t)形式模型參數(shù)值

將谷物EMC/ERH模型預(yù)測EMC值并與試驗(yàn)值進(jìn)行對(duì)比，結(jié)果表明谷物EMC/ERH模型對(duì)小麥、玉米、稻谷和大米在不同溫度下解吸和吸附過程中的平衡水分值擬合精度較高，玉米、稻谷和大米在解吸及吸附過程中擬合度R2均大于0.99，解決了WU(CAE)模型不能反求糧食平衡水分的問題。

圖1為谷物EMC/ERH模型預(yù)測20 ℃谷物解吸、吸附等溫線?？梢钥闯觯鞣N谷物預(yù)測與實(shí)際的等溫曲線基本重合。平衡水分EMC隨著相對(duì)濕度ERH的增加而增大，在相對(duì)濕度較高時(shí)平衡水分增長速度更快。

圖1 20 ℃谷物解吸/吸附平衡等溫線預(yù)測值與實(shí)際值

2.3 安全儲(chǔ)藏水分預(yù)測

通過圖1可知，小麥、玉米、稻谷、大豆、大米屬于第Ⅱ種類型等溫線。利用谷物EMC/ERH模型預(yù)測在不同溫度下的安全儲(chǔ)藏水分，如表3。預(yù)測結(jié)果基本與文獻(xiàn)[14-15]中利用經(jīng)典模型預(yù)測的安全水分值相同。

表3 安全儲(chǔ)藏水分

3 谷物的熱力學(xué)特性

積分吸著等熱Qst是指在一定溫度和水分活度條件下，單位質(zhì)量產(chǎn)品從液態(tài)發(fā)生汽化所需的能量。凈吸著等熱Δhst(微分焓)是指在系統(tǒng)溫度下，吸著等熱減去水汽化潛熱。通過經(jīng)典的Clausius-Clayperon方程[16]得到：

(3)

式中：aw為水分活度(小數(shù)表示)；T為溫度/℃；Δhst為凈吸著等熱/kJ·kg-1；R為氣體常數(shù)8.314/J· mol-1·K-1；Tk為溫度/K；C為模型的相關(guān)系數(shù)。假設(shè)Δhst不隨溫度變化，凈吸著等熱可以通過式(4)進(jìn)行計(jì)算[17]。

lnaw=-(Δhst/R)/T+C

(4)

吸著等熱可以通過自由水汽化潛熱L與凈吸著等熱Δhst的和來求解：

Qst=Δhst+L

(5)

式中：Qst為吸著等熱/kJ·kg-1；L為自由水汽化潛熱/kJ·kg-1；Δhst為凈吸著等熱/kJ·kg-1。

L可以通過式(6)求解：

L=2 502.2-2.39t

(6)

式中：t為平均溫度/℃。

吸著微分熵可以通過Gibbs-Helmhotz公式[18]利用ΔG吉布斯自由能進(jìn)行計(jì)算：

ΔS=(Δhst-ΔG)/T

(7)

式中：ΔS為吸著微分熵/kJ·kg-1·K-1；ΔG為自由能/kJ·kg-1·mol-1;可由式(8)求解：

ΔG=RTlnaw

(8)

在自由能上的水分吸著變化影響通常都伴隨著焓和熵的變化，因此將式(8)代入式(7)整理后得到式(9)：

lnaw=Δhst/RT-ΔS/R

(9)

在平衡水分一定時(shí)，利用-lnaw和1/T可作出不同水分的圖形，該圖為一條直線，通過其直線斜率可計(jì)算得到凈吸著等熱Δhst/kJ·mol-1，通過截距可計(jì)算得到ΔS/ kJ·mol·K-1[19-20]。

圖2為小麥、玉米、稻谷凈等量吸著熱和含水量的關(guān)系。在質(zhì)熱交換過程中，吸附熱是測定氣相變液相吸附釋放的熱量，解吸熱是打破水蒸汽分子和吸附劑表面分子間作用力所需的能量。從圖2可以看出，隨著平衡水分逐漸增大，焓Δhst、熵ΔS值迅速下降。解吸、吸附過程中當(dāng)平衡水分增大到20%左右，焓Δhst、熵ΔS值下降趨勢減緩。相同平衡水分的不同谷物凈等量解吸熱均高于凈等量吸附熱，在整個(gè)水分吸著過程中解吸所需的熱量大于吸附釋放的熱量，平衡水分大于25%后解吸與吸附過程的凈等量吸著熱相差較小。

圖2 不同水分谷物在解吸與吸附過程的焓熵值

通過式(6)計(jì)算得到自由水的汽化潛熱為44.28 kJ·mol-1與所求的凈吸著等熱求和可以得到吸著等熱。表4為小麥、稻谷、玉米解吸與吸附過程的吸著等熱曲線回歸方程。

表4 三大谷物吸著等熱回歸方程

注：Qst為吸著等熱/kJ·mol-1；EMC為平衡水分/%。

圖3為不同水分谷物在解吸、吸附過程中的吸著等熱。吸著等熱Qst隨著平衡水分的減小而增加，解吸等熱曲線體現(xiàn)更為顯著。解吸和吸附的吸著等熱值均高于水汽化潛熱，表明水分子和吸附點(diǎn)之間的結(jié)合能高于液相水分子之間的結(jié)合能。另外，不同種類谷物的解吸吸著等熱值均大于吸附吸著等熱值，且隨著平衡水分的增加解吸和吸附過程的吸著等熱的差值逐漸減少，這與Andre′ Luis等[17]和Benado等[21]闡述的一致。平衡水分大于25%后，吸著等熱趨于穩(wěn)定，與汽化潛熱值接近。

圖3 不同水分谷物在解吸與吸附過程的吸著等熱

4 結(jié)論

4.1 谷物EMC/ERH模型對(duì)小麥、玉米、稻谷和大米的平衡水分與相對(duì)濕度擬合精度較高。利用谷物EMC/ERH模型，可以直接求解糧食的平衡水分，解決了WU(CAE)模型不能反求糧食平衡水分的問題。

4.2 隨著平衡水分的逐漸增大，焓Δhst、熵ΔS值迅速下降。解吸、吸附過程中當(dāng)平衡水分增大到20%左右，焓Δhst、熵ΔS值下降趨勢減緩。相同平衡水分的不同谷物凈等量解吸熱值均高于凈等量吸附熱值，在整個(gè)水分吸著過程中解吸所需的熱量大于吸附釋放的熱量。平衡水分大于25%后，解吸與吸附過程中的凈等量吸著熱值相差較小。

4.3 吸著等熱Qst隨著平衡水分的減小而增加，解吸過程的吸著等熱值增加的更為顯著，解吸和吸附過程的吸著等熱值均高于水汽化潛熱。平衡水分大于25%后，吸著等熱值趨于穩(wěn)定，與汽化潛熱值接近。

[1]Henderson S M. A basic concept of equilibrium moisture [J]. Agricultural Engineering，1952, 33：29-32

[2]Thompson T L, Peart R M, Foster G. H. Mathematical simulation of corn drying-a new model [J]. Transaction of the ASAE, 1968, 24(3): 582-586

[3]Chung D S, Pfost H. B. Adsorption and desorption of water vapour by cereal grains and their products [J]. Transactions of the ASAE, 1967, 10: 549-557

[4]Halsey G. Physical adsorption on non-uniform surfaces [J]. Journal of Chemical Physics, 1948, 16: 931-937

[5]Oswin C R. The kinetic of package life III Isotherm [J]. Journal of Chemical Industry, 1946, 65: 419-421

[6]Blahovec J，Yanniotis S. Modified classification of sorption isotherms [J]. Journal of Food Engineering, 2009, 91:72-77

[7]李興軍, 吳子丹. 稻谷平衡水分的測定及EMC/ERH 等溫線方程的選擇[J]. 中國糧油學(xué)報(bào), 2010, 25(6):1-8

[8]李興軍，張?jiān)罚匚? 等. 我國小麥水分吸著等溫線4-參數(shù)GAB 方程研究[J]. 中國糧油學(xué)報(bào), 2011, 26(7): 82-88

[9]LS /T 1202—2002，儲(chǔ)糧機(jī)械通風(fēng)技術(shù)規(guī)程[S]

[10]吳子丹．綠色生態(tài)低碳儲(chǔ)糧新技術(shù)[M]. 北京: 中國科學(xué)技術(shù)出版社, 2011:299-319

[11]吳子丹, 李興軍. 大米和大豆CAE方程及平衡絕對(duì)曲線圖研究[J]. 中國糧油學(xué)報(bào), 2013, 28(1):76-81

[12]李興軍，王雙林，王金水. 小麥的平衡水分與吸著熱研究[J]. 河南工業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版，2009, 30(3):1-6

[13]李興軍，王雙林，張?jiān)? 玉米吸濕特性及其等溫線類型研究[J]. 中國糧油學(xué)報(bào), 2012, 27 (1):79-86

[14]李興軍. 谷物安全水分估算[J]. 糧食加工, 2011,36(3):41-45

[15]李興軍，劉丁，王雙林,等. 采用水分解吸方程估算我國稻谷安全水分[J]. 中國糧油學(xué)報(bào)，2011, 26(4):81-84

[16]Iglesias H, Chirife J. Prediction of the effect of temperature on water sorption isotherms of food material [J]. Journal of Food Technology, 1976, 11: 109-116

[17]Andre’ Luis D G, Paulo C C. Water sorption isotherms and thermodynamic properties of pearl millet grain [J]. International Journal of Food Science and Technology, 2010, 45: 828-838

[18]Rizvi S S H. Thermodynamic properties of foods in dehydration [J]. Engineering Properties of Foods,1995(4):223-309

[19]Tolaba M P, Peltzer M, Enriquez N, et al. Grain sorption equilibria of quinoa grains [J]. Journal of Food Engineering, 2004, 61: 365-371

[20]Mohammad N. Desorption isotherms and thermodynamic properties of fresh and osmotic-ultrasonic dehydrated quinces [J]. Journal of Food Processing and Preservation, 2013, 37: 381-390

[21]Benado A L, Rizvi S S H. Thermodynamic properties of water on rice as calculated from reversible and irreversible isotherms [J]. Journal of Food Science, 1985, 50: 101-105.

Heat and Moisture Balance New Model and Thermodynamic Properties of Grain

Chen Siyu1,2Wu Wenfu1Li Xingjun3Wu Zidan3Zhang Yaqiu1Zhang Zhongjie3Han Feng1

(College of Biological and Agricultural Engineering, Jilin University1, Changchun 130022)(College of Mechanical Engineering, Jiamusi University2, Jiamusi 154007)(Academy of State Administration of Grain3，Beijing 100037)

The equilibrium moisture content and equilibrium relative humidity models of grain are established under different temperatures. The static equilibrium test data of wheat, corn, rough rice and rice are fitted. The fitting results indicate that the model has high precision to describe the isotherm equations of grain and could calculate the equilibrium moisture content directly. In addition, the established model predicates the safe moistures under different temperatures and analyzes the thermodynamic properties of grains, which indicates that the isosteric of sorption heat increases with the decrease of grain equilibrium moisture, especially the increase of isosteric of desorption heat; isosteric of sorption heat is higher than latent heat of vaporization during desorption and sorption heat process; After the equilibrium moisture content reached 25% (w.b.), the isosteric of sorption heat regression equations is to be stable and approaches latent heat of vaporization. Based on the data analysis, sorption isosteric heat curvilinear regression regression equations of wheat, rice and maize are fitted, providing an available basis for the granary energy job management and the efficient reasonable energy-saving ventilation.

grain, moisture sorption isotherm, mathematical modeling, thermodynamic properties

TS210.1

A

1003-0174(2016)03-0110-05

時(shí)間：2016-01-29 14:22:36

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/11.2864.TS.20160129.1706.002.html

國家糧食局公益性行業(yè)科研專項(xiàng)(201313001-07、201313001-03-04)

2015-01-28

陳思羽，女，1980年出生，博士，糧食儲(chǔ)藏

韓峰，男，1978年出生，博士，農(nóng)業(yè)測控理論與技術(shù)、農(nóng)產(chǎn)品加工、智能機(jī)械