基于不同數學模型的魔芋水分吸濕等溫線模擬

吳紹鋒，彭桂蘭，邱兵濤，王彩霞（西南大學工程技術學院，重慶400715）

基于不同數學模型的魔芋水分吸濕等溫線模擬

吳紹鋒，彭桂蘭*，邱兵濤，王彩霞
（西南大學工程技術學院，重慶400715）

根據吸附原理，采用靜態調整環境濕度法，測定了魔芋在20、30、40℃3個溫度下，0%～98%水活度范圍內的平衡含水率，繪出魔芋的吸濕等溫線。結果顯示，魔芋的吸濕等溫線屬于II型等溫線；在一定的水活度下隨著溫度的升高魔芋的吸附能力下降。并以平均相對誤差和決定系數為評價指標，用八種數學模型對實驗數據進行擬合，結果表明Peleg模型對魔芋的吸濕等溫線擬合效果最好。

魔芋，平衡含水率，吸濕等溫線

魔芋別名磨芋、鬼芋，屬天南星科魔芋屬。魔芋是我國的一種重要的資源植物，具有重要的藥用、食用及工業價值［1-2］。但是魔芋收獲后，不宜長期儲存［3］，所以除少量鮮用外，大部分魔芋必須在短期內進行干燥，以便于儲藏、包裝和運輸。平衡含水率是研究食物干燥特性和安全貯藏的重要參數，吸濕等溫線對于食品的干燥、包裝和貯藏也有一定的指導意義［4-6］。很多學者對農產品的等溫線進行了研究，Oyelade OJ、Tunde-Akintunde TY等在溫度為27～40℃，水分活度為0.10～0.80的條件下研究了玉米淀粉的吸附和解吸等溫線，并利用五種模型對吸附和解吸等溫線進行了擬合，得出GAB是最適合在系統上描述玉米粉吸附和解吸等溫線的方程［7］；Aviara NA等在水分活度0.07～0.98，溫度40～70℃高溫條件下對大豆的吸附等溫線進行了研究，并用四種等溫線方程對吸附等溫線進行了擬合，結果顯示大豆的吸附等溫線是S型II型等溫線，Oswin模型是預測平衡水分含量的最佳模型［8］；Al-Muhtaseb H，McMinn WAM等研究表明Peleg、GAB和Ferro-Fontan模型對馬鈴薯淀粉、高支鏈淀粉和高直鏈淀粉的水分吸附等溫線都有一個很好的擬合［9］；王辰，劉良忠等取相對濕度10%～90%范圍，測定20、30、40℃條件下芹菜的等溫解吸和吸濕平衡含水率，結果表明，Modified-Hendenson模型對芹菜平衡含水率的擬合效果最好［4］；胡坤、張家年等研究了5種最常用的數學模型對中國不同類型的稻谷吸附與解吸等溫線數據的擬合效果，確定了最佳擬合模型及擬合參數［10］；朱恩龍，楊昭等采用靜態稱重法測定了青豆種子在20～40℃和水分活度在0.112～0.946范圍的等溫線，在所采用的模型中，Halsey模型擬合優度最佳［11］；周洲在11.29%～95.96%相對濕度和4個溫度（10、20、30、40℃）條件下，采用靜態稱量法測定了國產和進口大豆的吸附與解吸平衡含水率，發現MGAB和MOE模型分別是描述國產大豆吸附與解吸等溫線的最佳模型［12］；彭桂蘭測定了玉米淀粉在30、45、60℃溫度不同水活度下的吸濕和解吸等溫線，發現GAB模型、Peleg模型和Henderson模型擬合的好，并且用BP神經網絡建立了一個新的擬合程度優于其他的吸附等溫數學模型［5］；王雙林、郭道林等采用靜態稱重法對6個玉米品種的吸附和解吸等溫線進行測定，并用6個非線性回歸方程描述吸附及解吸等溫線，得出MCPE是最佳的數學擬合模型［13］；此外，李興軍、王雙林等研究了玉米吸濕特性及其等溫線類型和小麥的平衡水分、小麥水分吸著等溫線4-參數GAB方程研究［14-15］；周玉龍、賈富國等研究了貯藏溫濕度對糙米平衡含水率的影響［16］等。但是，目前還沒有學者對魔芋吸濕等溫線進行研究。本研究是在一定的水活度和溫度的條件下，測定魔芋的平衡含水率，分析溫度和水活度對平衡含水率的影響，作出吸濕等溫線，并采用八種數學模型對實驗數據進行擬合，得出最優擬合模型，為魔芋的貯藏和保存提供理論依據。

1　材料與方法

1.1材料與儀器

魔芋片在帶有溫度控制的鼓風干燥箱中干燥成絕干物質作吸濕實驗用；濃硫酸分析純，重慶川東化工有限公司化學試劑廠；干燥器上海精英實驗器材試劑大全；鋁盒河北滄州中建精密儀器有限公司。

DHS-250恒溫恒濕熱實驗箱（控溫精度為0.1℃）上海林頻儀器設備有限公司；AL204電子天平（測量精度為0.0001g） 梅特勒-托利多儀器（上海）有限公司。

1.2實驗方法

1.2.1魔芋吸濕等溫線測定根據實際生產中的貯藏環境條件，選取實驗溫度為：20、30、40℃。將干燥好的魔芋片樣品快速稱取（2.0000±0.01）g放入已經恒重的鋁盒中，然后放入干燥器上部，干燥器下部放入不同濃度的硫酸溶液。硫酸溶液用于保持干燥器中的相對濕度，不同硫酸濃度在不同溫度下的水活度見表1［17］。把干燥器分別放在溫度為20、30、40℃的恒溫箱中。在此期間每天對樣品稱重1次，直到樣品重量的變化不超過0.001g時達到平衡，大約需要6d左右。稱重和重新放置大約需要30s，稱重時在大氣中吸收的水分可以忽略不計，每個實驗重復3次，魔芋的平衡含水率通過靜態稱重法得到［16，18］。

表1　不同濃度和溫度的硫酸溶液的水分活度Table 1　Water activities of sulphuric acid solutions at selected concentrations and temperatures

1.2.2魔芋吸濕等溫線擬合模型的選擇許多模型已被用于預測溫度、平衡水分和水分活度的關系，其中以溫度和水活度為自變量的平衡水分含量的函數在國內外已經進行了大量的測定和分析，常見的函數模型如表2［19］所示。

表2　擬合實驗數據的數學模型Table 2　Mathematical model fit for the experimental data

本實驗選用上述8種數學模型，利用數據分析軟件，對所測實驗數據分別進行擬合，根據評價指標判斷每個模型對吸濕等溫線的擬合效果，并選擇描述魔芋吸濕等溫線的最佳模型。

1.2.3評價指標本實驗采用平均相對誤差E和決定系數R2來衡量模型擬合效果的好壞，平均相對誤差在研究中被廣泛的應用，E值低于10%說明擬合度比較好［20］；R2越大說明模型與等溫線的擬合效果越好［12］。

平均相對誤差：

式中，X：測定值；X’：擬合值；X：平均測定值；n：數據個數。

2　結果與分析

2.1吸濕等溫線的建立

分別在20、30、40℃及不同的水活度條件下，測定魔芋的平衡含水率，得到魔芋的吸濕等溫線，如圖1所示。從圖1中可以看出，在同一溫度下平衡含水率隨著水活度的增加而增加，說明魔芋吸濕等溫線的形狀為S型，屬于II型等溫線［21］。

這些數據表明，在一定的水活度下，隨溫度的升高，平衡含水率值下降，魔芋的吸水性下降。這種趨勢是由于溫度引起物質內部的化學、物理變化而導致水的親合活性點減少的結果［22］。

2.2擬合檢驗

圖1　魔芋吸濕等溫線Fig.1　Experimental adsorption isotherms for konjac

對魔芋在20、30、40℃的實驗數據進行非線性回歸，對八種經典模型進行擬合，計算模型參數、決定系數、平均相對誤差，如表3所示。

表3　數學模型的相關參數及評價指標Table 3　Relevant parameters and evaluation index of the mathematical model

從表3可以看出，除Smith模型外，其余模型的決定系數都在0.98以上，說明這幾種模型都可以用來擬合魔芋的吸濕等溫線，但模型的擬合效果存在差異。由平均相對誤差值可知，Peleg模型在整個水活度范圍內對魔芋的吸濕數據擬合的最好，Peleg模型的相對誤差E值為6.25%～10.01%；其次擬合效果較好是GAB模型，擬合E值在9.95%～14.95%；Ferro-Fontan模型在溫度較低時擬合效果好，較高溫度擬合效果差。Mod-BET模型、Halsey模型、Oswin模型、Henderson模型不能很好的擬合魔芋的吸濕等溫線；Smith模型無法擬合魔芋的吸附等溫線。

3　結論

魔芋吸濕等溫線的形狀為S型，屬于II型等溫線。在同一溫度下，魔芋的平衡含水率隨水活度的增加而增加；在一定的水活度下，魔芋的平衡含水率隨溫度的升高而下降。Peleg模型對魔芋吸濕等溫線的擬合比其他模型好。通過研究魔芋水分吸濕等溫線可以確定魔芋及其制品的安全儲藏的含水率，且可指導產品包裝的方式和材料。

［1］劉佩瑛.魔芋學［M］.北京：中國農業出版社，2004：1-11.

［2］李靜，王志民，張忠，等.魔芋的應用價值與開發前景［J］.西昌學院學報，2006，20（4）：17-19.

［3］董坤.魔芋貯藏期失水率研究［J］.云南農業科技，2011（1）：11-12.

［4］王辰，劉良忠，張金木.芹菜平衡含水率的研究［J］.食品科技，2007（8）：251-254.

［5］彭桂蘭，陳曉光，吳文福，等.玉米淀粉水分吸附等溫線的研究及模型建立［J］.農業工程學報，2006，22（5）：176-179.

［6］邱兵濤，彭桂蘭，賈琳鵬.金銀花水分吸濕等溫線的研究及模型建立［J］.農機化研究，2012（9）：32-34.

［7］Oyelade OJ，Tunde-Akintunde TY，Igbeka JC，et al.Modelling moisture sorption isotherms for maize flour［J］.Journal of Stored Products Research，2008（44）：179-185.

［8］Aviara NA，Ajibola OO，Oni SA，et al.Sorption Equilibrium and Thermodynamic Characteristics of Soya Bean［J］.Biosystems Engineering，2004，87（2）：179-190.

［9］Al-Muhtaseb AH，McMinn WAM，Magee TRA.Water sorption isotherms of starch powders Part 1：mathematical description of experimental data［J］.Journal of Food Engineering，2004（61）：297-307.

［10］胡坤，張家年.稻谷水分吸附與解吸等溫線擬合模型的選擇及其參數優化［J］.農業工程學報，2006，22（1）：153-156.

［11］朱恩龍，楊昭，尹海蛟，等.青豆種子等溫線及平衡含水率模型［J］.天津大學學報，2012，45（5）：401-404.

［12］周洲.大豆水分吸著等溫線及夏季儲藏的安全水分［D］.南京：南京財經大學，2012：1-25.

［13］王雙林，郭道林，李興軍，等.玉米平衡水分測定及等溫線方程確定［J］.中國糧油學報，2010，25（11）：33-36.

［14］李興軍，王雙林，張元娣，等.玉米吸濕特性及其等溫線類型研究［J］.中國糧油學報，2012，27（1）：81-86.

［15］李興軍，張元娣，秦文，等.我國小麥水分吸著等溫線4-參數GAB方程研究［J］.中國糧油學報，2011，26（7）：82-87.

［16］周玉龍，賈富國，張強，等.貯藏溫濕度對糙米平衡含水率的影響［J］.中國糧油學報，2011，26（12）：78-81.

［17］Ruegg M.Calculation of the activity of water in sulphuric acid solutions at various temperature［J］.Lebensmittel-Wissenchaft and Technology，1980（13）：22-24.

［18］李興軍，王雙林，王金水.谷物平衡水分研究概況［J］.中國糧油學報，2009，24（11）：137-143.

［19］彭桂蘭.玉米淀粉濕法加工浸泡機理、工藝及智能專家系統的研究［D］.長春：吉林大學，2006.

［20］Lomaruo CJ，Bakshi AS，Labuza TP.Evaluation of food moisture sorption isotherm equations Part I Fruit，vegetable and meat products［J］.Lebensmittel-Wissenchaft and Technology，1985（18）：111-117.

［21］Ertugay MF，Certel M.Moisturesorptionisotherms of cereals at different temperatures［J］.Nahrung，2000（44）：107-109.

［22］Mazza G，LeMaguer M.Dehydration of onion：Some theoretical practical considerations［J］.Journal of Food Technology，1980（15）：181-194.

Modelling moisture sorption isotherms for konjac based on different mathematical models

WU Shao-feng，PENG Gui-lan*，QIU Bing-tao，WANG Cai-xia
（College of Engineering and Technology，Southwest University，Chongqing 400715，China）

This study was designed for the equilibrium moisture content and moisture sorption isotherms for konjac by means of static-adjusting ambient humidity which based on adsorption theory.It was determined at different temperature 20，30，40℃and different relative humidity ranging from 0%to 98%.Results showed that the isotherms of konjac exhibited Type II and the sorption capacity decreased with the increase of temperature at constant water activity.The experimental data were compared with eight widely recommended models by using the evaluation index of average relative error and coefficient of determination.Results indicated that the Peleg model appeared to be most suitable for fitting the moisture sorption isotherms data for the konjac.

konjac；equilibrium moisture content；sorption isotherms

TS255.1

A

1002-0306（2015）02-0126-04

10.13386/j.issn1002-0306.2015.02.019

2014-04-11

吳紹鋒（1991-），男，碩士研究生，研究方向：農產品加工與品質檢測。

彭桂蘭（1966-），女，博士研究生，教授，研究方向：農產品加工與品質檢測。

西南大學博士基金（SWUB2007021）。