稻谷平衡水分的測定及 E MC/ERH等溫線方程的選擇

李興軍 吳子丹

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

稻谷平衡水分的測定及 E MC/ERH等溫線方程的選擇

李興軍 吳子丹

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

采用靜態(tài)稱重法測定 7個(gè)秈稻和 1個(gè)粳稻品種的平衡水分,并利用常用的 6個(gè) EMC/ERH方程(修正BET、修正 Chung-Pfost、修正 Guggenheim-Anderson-deBoer、修正 Henderson、修正Os win及 Strohman-Yoerger)擬合吸附/解吸等溫線。在 ERH 11.3%～96%范圍內(nèi),修正 Chung-Pfost(MCPE)和 Strohman-Yoerger (STYE)擬合最佳。采用含有 3個(gè)系數(shù)、易于轉(zhuǎn)化為M=f(r.h.,t)或 r.h.=f(M,t)表達(dá)形式的MCPE方程,對 8個(gè)稻谷品種的吸著等溫線進(jìn)行擬合。以M =f(r.h.,t)形式表達(dá)的MCPE方程系數(shù) C1、C2、C3,在擬合吸附等溫線時(shí)分別是 784.894、143.337、0.174,在擬合解吸等溫線時(shí)分別是 588.376、59.026、0.180,吸著平均值分別是 638.444、87.074、0.177。稻谷吸附與解吸等溫線之間存在滯后現(xiàn)象,隨著溫度增加,滯后環(huán)寬度減少。

平衡水分 稻谷 解吸 吸附 數(shù)學(xué)模型 通風(fēng) 干燥

稻谷是世界上重要的谷物,總產(chǎn)量接近小麥和玉米之和。我國稻谷年產(chǎn)量約 1.8億噸[1],居世界首位。新收獲的稻谷含水量高,需干燥后儲存或加工。干燥和儲藏條件均影響稻谷品質(zhì)。調(diào)控這些條件,除了薄層干燥方程(修正BET)外,還需要了解稻谷平衡含水量(EMC)與一定溫度下干燥空氣或通風(fēng)空氣的平衡相對濕度 (ERH)之間的關(guān)系。EMC/ ERH吸著(吸附/解吸)等溫線方程通常用于描述這種關(guān)系。國外已經(jīng)發(fā)展了 200多個(gè)理論的、半理論的及經(jīng)驗(yàn)方程式,用于描述重要糧食種類的吸著等溫線,但是適合稻谷 EMC/ERH的方程報(bào)道有限[2]。Agrawal等[3]采用 11個(gè)模型分析了 4個(gè)稻谷的解吸數(shù)據(jù),發(fā)現(xiàn) Strohman等[4]提出的水分解吸等溫線方程最成功。Pfost等[5]對文獻(xiàn)報(bào)道的玉米吸著數(shù)據(jù)分析了 5個(gè)方程,發(fā)現(xiàn) Strohman-Yoerger方程 (STYE)標(biāo)準(zhǔn)誤差最小,但是含有 4個(gè)系數(shù)的方程,認(rèn)為含有3個(gè)系數(shù)的修正 Henderson(MHE)和修正 Chung-Pfost(MCPE)方程適合分析包含稻谷在內(nèi)糧食吸著行為。Kachru等[6]研究表明稻谷的 EMC是相對濕度的多項(xiàng)式函數(shù),也是溫度的線性函數(shù),而 Putranon等[7]的研究表明稻谷的 ERH是水分的多項(xiàng)式函數(shù),也是溫度的線性函數(shù)。Zuritz等[8]利用 Day-Nelson方程、修正 Henderson及 4個(gè)系數(shù)的 Chung-Pfost方程分析稻谷的 EMC/ERH數(shù)據(jù),依據(jù)均方根誤差、方差檢驗(yàn)分析,發(fā)現(xiàn)這 3個(gè)方程沒有差異,均適合描述EMC數(shù)據(jù)。Chen等[9]比較了 4個(gè) EMC/ERH方程MHE、MCPE、修正 Halsey(MHAE)及修正 Os win (MOE)對 3種來源的稻谷和糙米數(shù)據(jù)的擬合性,除了MHAE通常對高蛋白的油料種子適合外,其他 3個(gè)方程部分或全部地適合稻谷數(shù)據(jù)。陳加忠等[10]指出,55～60℃風(fēng)干的臺灣臺農(nóng) 67號稻谷吸附數(shù)據(jù)以MHE和MCPE為適合模型。就植物來源的農(nóng)產(chǎn)品水分吸著關(guān)系,ASABE標(biāo)準(zhǔn)D245.5推薦適合的方程是MHE、MCPE、MHAE及 MOE[11]。Guggenheim-An2 derson-deBoer(GAB)方程也被 D245.5推薦,但是它不包含吸附/解吸過程中溫度產(chǎn)生的效果。Sun[12]收集整理了 18種資料來源的 763個(gè)稻谷試驗(yàn)數(shù)據(jù),得出 STYE模型最適于描述稻谷的等溫線。胡坤等[13-14]指出,不同品種的秈稻吸附與解吸等溫線之間沒有明顯的差異,而初始含水量、稻谷類型(秈稻、粳稻及糯稻)及加工程度 (原稻、糙米及白米)顯著地影響稻谷的吸附與解吸等溫線,且美國農(nóng)業(yè)工程學(xué)會推薦的MCPE及參數(shù)并不能與中國稻谷的吸附與解吸等溫線數(shù)據(jù)很好的擬合。STYE最適合于擬合秈稻、粳稻的吸附與解吸等溫線及糯稻的吸附等溫線,而MOE最適合擬合糯稻的解吸等溫線。

評價(jià)稻谷吸著行為的MCPE、MHE、MOE及STYE擬合方程,需要大量的稻谷 EMC/ERH數(shù)據(jù)。在小麥平衡水分測定分析基礎(chǔ)上[15],本研究測定分析了我國稻谷主產(chǎn)區(qū)8個(gè)品種的 EMC/ERH數(shù)據(jù),選擇了6個(gè)常用 EMC/ERH方程對適合溫度范圍的數(shù)據(jù)組進(jìn)行擬合,并選擇了最合適的MCPE描述稻谷吸附/解吸數(shù)據(jù)組,以期為我國稻谷儲藏和加工提供基礎(chǔ)性數(shù)據(jù)。

1 材料與方法

1.1 儀器與藥品

AB135-S/FACT分析天平:上海梅特勒 -托利多公司;PRX-350A智能人工氣候箱(溫度偏差 ±0.8℃):寧波海曙賽福實(shí)驗(yàn)儀器廠,經(jīng)常用標(biāo)準(zhǔn)溫度計(jì)(0～50℃范圍)校正;DHG9040A智能干燥箱 (溫度偏差 ±1.0℃):杭州藍(lán)天儀器廠,放置樣品的位置用標(biāo)準(zhǔn)溫度計(jì)校正;干燥器:內(nèi)徑 240 mm;鋁盒:直徑 5 cm、厚度 2 cm;250 mL玻璃廣口瓶;9號橡皮塞;由銅絲網(wǎng)制作的小桶:直徑 2.8 cm、高度 4 cm。醋酸鉀、氯化鎂、碳酸鉀、硝酸鎂、氯化銅、氯化鈉、氯化鉀、硝酸鉀:優(yōu)級純,天津市光復(fù)科技發(fā)展有限公司;氯化鋰:優(yōu)級純,北京國藥集團(tuán)化學(xué)試劑有限公司;五氧化二磷固體:化學(xué)純,北京精益化學(xué)試劑有限公司。

1.2 樣品與試驗(yàn)方法

采用 8個(gè)稻谷品種(表 1)。對于吸附曲線樣品,選取干凈、飽滿、無破裂的稻谷樣品在40.5℃保溫箱內(nèi),將含水率降到 7%～8%,然后在干燥器中用 P2O5固體脫水到 5%干基以下。新收獲的高水分稻谷常德218(含水率 19.5%)、常德 207(含水率 17.9%)直接用于測定解吸曲線樣品,其他稻谷品種由 5.0%干基以下加水調(diào)到 22%,在冰箱平衡 2周后用于測定解吸曲線。

表1 試驗(yàn)用稻谷品種樣品

采用靜態(tài)稱重法測定平衡水分,即利用表 2的飽和鹽溶液在 5種恒定溫度(10、20、25、30及 35℃)下產(chǎn)生恒定的蒸汽壓。將 3.000～5.000 g稻谷籽粒倒入銅絲小桶,懸掛在橡膠塞下方掛鉤上,然后置于盛 60 mL過飽和鹽溶液的廣口瓶,塞緊橡膠塞,保持小桶底部距離鹽溶液上部 2～3 cm。進(jìn)樣后在 25～30 d這一期間每隔 1天稱量1次,兩次連續(xù)稱重不超過 2 mg時(shí),此時(shí)樣品中含水率就是該溫度、ERH下的平衡水分(EMC)。置于飽和硝酸鉀溶液上空的樣品,一旦觀察到發(fā)霉即取出樣品烘干稱重。糧食含水率測定采用整粒烘干法,(103.0±1.0)℃烘干20～28 h。

表 2 九種飽和鹽不同溫度產(chǎn)生的平衡相對濕度(%)[16]

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

2 結(jié)果分析

2.1 自然高水分樣品與調(diào)水樣品解吸行為比較

從圖 1看出,在 30℃、ERH 11.3%～96.0%內(nèi),自然高水分(17.9%)的常徳 207稻谷樣品與調(diào)制樣品(22%)的解吸行為一致。

表3 水分吸著等溫線數(shù)學(xué)模型

圖 1 常德 207稻谷自然高水分樣品與調(diào)制樣品 30℃解吸等溫線

2.2 稻谷吸附與解吸等溫線實(shí)測值分析

表 4和表 5分別是 8個(gè)稻谷品種在 5種溫度(10、20、25、30及 35℃)、9種 ERH下的吸附與解吸平衡含水率均值。標(biāo)準(zhǔn)差分析表明,品種之間平衡水分含量變異不明顯。將表 4和表 5中數(shù)據(jù)作成解吸與吸附等溫線,如圖 2。解吸與吸附等溫線均呈 S型曲線。在一定溫度下,隨著 ERH增大,解吸或吸附數(shù)值增大。在恒定的 ERH下,解吸或吸附數(shù)值均隨著溫度的增加而減少。在 ERH 11%～85%范圍內(nèi),解吸等溫線滯后吸附等溫線 (圖 3),隨著溫度增加,滯后環(huán)跨度變小。

圖 2 8個(gè)稻谷品種平均解吸與吸附等溫線

圖 3 20℃下的 8個(gè)稻谷品種平均解吸與吸附曲線

2.3 稻谷吸附與解吸等溫線數(shù)學(xué)模型分析

從表 6、表 7的 EMC/ERH關(guān)系數(shù)學(xué)模型 r.h.= f(M,t)的參數(shù)估算看出,R2>0.970、MRE<10%,除了MOE方程回歸模型殘差分布對有些數(shù)據(jù)組是模式以外,MCPE、MHE、STYE及MGAB回歸模型殘差分布均是隨機(jī)的,表明所選數(shù)學(xué)模型MCPE、MHE、STYE及MGAB在 RH 11.3%～96%范圍內(nèi)均適合稻谷吸附與解吸等溫線擬合。STYE含有 4個(gè)系數(shù), MCPE、MHE和MGAB含有 3個(gè)系數(shù),可以從這些數(shù)學(xué)模型中探討最合適的方程。

表 4 8個(gè)稻谷品種不同溫度解吸數(shù)據(jù)均值

表 5 8個(gè)稻谷品種不同溫度吸附數(shù)據(jù)均值

表 6 稻谷品種解吸等溫線模型 r.h.=f(M,t)參數(shù)估算

續(xù)表

表 7 稻谷品種吸附等溫線模型 r.h.=f(M,t)參數(shù)估算

續(xù)表

表8 16組稻谷吸著數(shù)據(jù)適合的數(shù)學(xué)模型統(tǒng)計(jì)參數(shù)平均值

表 9 稻谷吸著數(shù)據(jù)最佳擬合方程的參數(shù)及統(tǒng)計(jì)參數(shù)比較

2.4 稻谷吸附與解吸最佳等溫線數(shù)學(xué)模型篩選

為了比較了不同表達(dá)形式的常用 ERH/EMC模型,對 16組稻谷吸著等溫線的擬合結(jié)果,統(tǒng)計(jì)分析結(jié)果總結(jié)如表 8。以 r.h.=f(M,t)形式表達(dá)的 5個(gè)方程,按照平均統(tǒng)計(jì)參數(shù) (R2、RSS、SEE、MRE)對比排序,STYE>MCPE>MHE>MGAB>MOE。以M=f(r.h.,t)形式表達(dá)的 5個(gè)方程中,BET在 RH 11.3%～50%范圍內(nèi)是最佳的吸著等溫線擬合方程,MCPE、MHE、MOE及MGAB均在RH 11.3%～96%范圍內(nèi)適合吸著等溫線。綜合考慮適合吸著數(shù)據(jù)方程的參數(shù)數(shù)目、表達(dá)形式轉(zhuǎn)化等因素,研究結(jié)果表明MCPE方程是最佳的吸著等溫線的擬合方程。

從表 9可以看出,以 r.h.=f(M,t)形式表達(dá)的MCPE和 STYE方程,吸附等溫線與解吸等溫線的擬合方程參數(shù)不同,表明稻谷吸附于解吸之間存在滯后現(xiàn)象。以M=f(r.h.,t)形式表達(dá)的MCPE方程,吸附等溫線與解吸等溫線的擬合方程參數(shù)也不同。

圖4 MCPE方程預(yù)測的10、20及30℃稻谷解吸與吸附等溫線

從圖 4看出,不同溫度下,稻谷吸附曲線與解吸曲線之間存在滯后現(xiàn)象。隨著溫度增加,滯后環(huán)的寬度和跨度變小。從圖5看出,在20℃下,稻谷品種之間吸附或解吸數(shù)據(jù)沒有差異。

圖5 MCPE方程預(yù)測的20℃稻谷解吸與吸附數(shù)值與8個(gè)品種實(shí)測數(shù)值比較

3 討論

Sun[12]分析了包含一系列溫度范圍的 18組稻谷吸附/解吸等溫線的 763個(gè)試驗(yàn)數(shù)據(jù),認(rèn)為 Strohman -Yoerger模型適合描述稻谷的等溫線。胡坤等[14]分析了包括秈稻、粳稻及糯稻在內(nèi)的 5個(gè)品種的吸附與解吸等溫線數(shù)據(jù),認(rèn)為 Strohman-Yoerger模型適合擬合秈稻、粳稻的吸著等溫線,MOE適合擬合糯稻的吸著等溫線。研究表明,Strohman-Yoerger方程以 r.h.=f(M,t)形式、MCPE方程以M =f(r.h., t)形式和 r.h.=f(M,t)形式均適合 7種秈稻和 1種粳稻的吸著等溫線(表 8)。美國農(nóng)業(yè)工程學(xué)會也推薦MCPE方程作為稻谷吸著等溫線擬合方程[11]。MCPE方程的優(yōu)點(diǎn)是含有 3個(gè)系數(shù),易于轉(zhuǎn)化為M =f(r.h.,t)或 r.h.=f(M,t)表達(dá)形式,而 Strohman -Yoerger方程是 4個(gè)系數(shù),不易轉(zhuǎn)化為M =f(r.h., t)的表達(dá)形式。因此,選擇MCPE方程作為稻谷吸著等溫線擬合方程。

通過對比分析 8個(gè)稻谷品種吸附/解吸實(shí)測值與MCPE方程擬合值,實(shí)測值分散地分布在擬合值周圍,品種之間差異不明顯 (圖 5)。因此,測定或收集較多品種的 ERH-EMC數(shù)據(jù),能夠?qū)⑵贩N、測量技術(shù)及測量準(zhǔn)確性的影響最小化,使隨機(jī)誤差和系統(tǒng)誤差平均化,擬合的數(shù)學(xué)模型參數(shù)更加準(zhǔn)確。

解吸等溫線用于稻谷干燥過程。新收獲的稻谷含水量是 16%～28%濕基,干燥是收獲后的必要處理,在濕熱地區(qū)非常重要。干燥處理推遲 1天,目的是讓稻谷籽粒內(nèi)部的水分梯度、單個(gè)籽粒之間的水分梯度能夠達(dá)到平衡。然后在特定條件下緩慢地干燥,將水分降低到 13.5%濕基,以到達(dá)最好的品質(zhì)保持效果。吸附等溫線用于描述加濕過程[17-18]。在現(xiàn)代倉儲管理中,儲藏期間通風(fēng)可能發(fā)生加濕過程。為了降低倉內(nèi)糧食的溫度,在夜間溫度很低時(shí)連續(xù)通風(fēng)。但是夜間相對濕度較高,引起通風(fēng)入口區(qū)域糧食嚴(yán)重吸濕,這種吸濕需要小心地檢測,它是糧食安全儲藏的威脅。因此需要選擇最合適的吸附方程。MCPE方程以M =f(r.h.,t)形式和 r.h.=f(M,t)形式預(yù)測稻谷的吸附和解吸曲線。

7種秈稻和 1種粳稻的吸著等溫線的擬合方程MCPE,以M =f(r.h.,t)形式表達(dá)的方程參數(shù) C1、C2及 C3,在擬合吸附等溫線時(shí)分別是 784.894、143.337、0.174,在擬合解吸等溫線時(shí)分別是 588. 376、59.026、0.180,吸著平均值分別是638.444、87. 074、0.177。

志謝:感謝任錫洪先生對糧食平衡水分測定的指導(dǎo)。

[1]聶振邦.2008中國糧食發(fā)展報(bào)告[M].北京:經(jīng)濟(jì)管理出版社,2008

[2]李興軍,王雙林,王金水.谷物平衡水分研究概況[J].中國糧油學(xué)報(bào),2009,24(11):137-145

[3]Agrawal K,Clary B. Investigations into the theories of de2 sorption isother ms for rough rice and peanuts[J].J Food Sci,1971,36:919-924

[4]Strohman R,Yoerger R.A new equilibrium moisture content equation[J].TransASAE,1967,10:675-677

[5]Pfost H B,Maurer S G,ChungD S,et al.Summarizing and reporting equilibrium moisture data of grains[M].St.Jo2 seph,M I,USA:ASAE,1976.

[6]Kachru R P,Matthes R K.The behaviour of rough rice insorption[J].J Agric Eng Res,1976,21:405-416

[7]Putranon R,BowreyR G,Eccleston J.Sorption isotherms for two cultivars of paddy rice grown in Australia[J]. Food Technology in Australia,1979,31:510-515

[8]Zuritz C,Singh P R,Moini SM,et al.Desorption isotherms of rough rice from 10℃ to 40℃ [J].TransASAE,1979, 22:433-436,440

[9]Chen C C,Morey R V.Comparison of four EMC/ERH equa2 tions[J].TransASAE,1989,32:983-990

[10]陳加忠,賴建洲,曹之祖.稻谷平衡相對濕度性質(zhì)之研究[J].中華農(nóng)業(yè)研究,1990,39:347-366

[11]ASAE Standards.D245.5:Moisture relationships of plantbased agricultural products[S].St.Joseph,M I,USA: ASAE,2002.

[12]Sun D W.Comparison and selection of EMC/ERH isotherm equations for rice[J].J Stored Prod Res,1999,35:249-264

[13]胡坤,張家年.稻谷吸附與解吸等溫線影響因子的研究[J].糧食與飼料工業(yè),2004(10):30-33

[14]胡坤,張家年.稻谷水分吸附與解吸等溫線擬合模型的選擇及其參數(shù)優(yōu)化[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22:153-156

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

[16]JayasDS,Mazza G.Equilibrium moisture characteristics of safflower seeds[J].TransASAE,1991,34:2099-2103

[17]吳子丹.儲糧機(jī)械通風(fēng)的計(jì)算機(jī)控制 [J].糧食儲藏, 1987,16(4):28-31

[18]吳子丹.再談儲糧機(jī)械通風(fēng)的計(jì)算機(jī)控制[J].糧食儲藏,1991,20(3):3-9.

Fiting Parameters of E MC/ERH Model for Chinese Rough Rice

Li Xingjun Wu Zidan
(Academy of the State Administration of Grains,Beijing 100037)

Equilibrium moisture content(EMC)data for rice grain of eight cultivarswere obtained by the gravi2 metric method at 11%to 96%relative humidity(RH)and 10,20,25,30,and 35℃above nine saturated salt so2 lutions.Six commonly used mathematic models,BET,Guggenhei m– Anderson-deBoer(GAB),Modified-Chung-Pfost(MCPE),Modified-Henerson(MHE),Modified-Os win(MOE),and Strohman-Yoerger (STYE),were fitted to the data with evaluating the coefficientof determination,residue sum-of-squares,standard error of estimate,mean relative percent error,and residual plots.Results:The best-fitted equations to the EMC/ ERH data areMCPE and STYE,butMCPE is three-parameter,readily transfor med equation and adopted in this study.The three parameters,C1,C2,and C3of theMCPE in a for m ofM.=f(r.h.,t)are 784.894,143.337,and 0.174 for adsorption isother m of rough rice grain,and 588.376,59.026,and 0.180 for desorption isother m of rough rice grain,respectively.The three parameters,C1,C2,and C3for average values of desorption and adsorption data are 638.444,87.074,and 0.177,respectively.At a constant RH,EMC decreases with increasing temperature. There are hysteresis loops between desorption isotherms and adsorption isother ms.

equilibrium moisture content,rice,adsorption,desorption,mathematicmodel,ventilation,drying

TS210.1 文獻(xiàn)標(biāo)識碼:A 文章編號:1003-0174(2010)06-0001-08

國家教育部留學(xué)歸國科研啟動基金(Z1006)

2009-07-20

李興軍,男,1971年出生,副研究員,糧食生理生化