錐栗脆球微波干燥動(dòng)力學(xué)模型研究

孫輝 毛志幸 陳宗道

摘? 要：為推動(dòng)閩北錐栗精深加工技術(shù)的發(fā)展，提高錐栗脆球的干燥效率和品質(zhì)，以干基含水率、干基失水速率和感官質(zhì)量評(píng)分為指標(biāo)，探索微波強(qiáng)度、裝載量及初始微波功率密度（微波功率與裝載量的比值，Mp0）對(duì)蒸煮熟化錐栗球脫水規(guī)律的影響。結(jié)果表明，微波功率和裝載量均對(duì)干燥時(shí)間有顯著影響，Mp0越大，干燥速率越快，其中恒速階段隨著Mp0增大而縮短，而且Mp0對(duì)熟化錐栗質(zhì)地結(jié)構(gòu)的影響更顯著，Mp0在5.5 W/g時(shí)錐栗脆球感官評(píng)分最高（80），膨化率最大（2.92），脆度（3.51 N）最佳，孔隙率為61.4%，水分活度低于0.3。然后通過Matlab建立錐栗脆球水分比與時(shí)間的動(dòng)力學(xué)模型，發(fā)現(xiàn)Page模型擬合度很高（0.99），可以闡明錐栗脆球微波脫水各階段對(duì)其干燥速率的影響;干燥后樣品的顆粒含水量隨機(jī)分布符合六次多項(xiàng)式分布，與平均含水量一致。說明根據(jù)錐栗脆球微波干燥的動(dòng)力學(xué)規(guī)律，可以通過分段變功率的方式進(jìn)一步優(yōu)化錐栗脆球干燥制品的干燥參數(shù)。

關(guān)鍵詞：錐栗脆球;微波功率密度;轉(zhuǎn)換含水率;動(dòng)力學(xué)模型

中圖分類號(hào)：TS255.6????? 文獻(xiàn)標(biāo)識(shí)碼：A

Microwave Drying Dynamic Model of Castanea Henryi Crisp Ball

SUN Hui1， MAO Zhixing2， CHEN Zongdao3*

1. College of Tea and Food Science， Wuyi University， Wuyi， Fujian 354300， China; 2. Chongqing City Agricultural Machinery Appraisal Station， Chongqing 402160， China;3. Southwest University， Chongqing 400715， China

Abstract： In order to obtain Castanea henryi crisp ball products with a fast drying rate and high quality by microwave， a single factor experiment was used to explore the effects on dry basis moisture content， dry basis water loss rate and sensory quality of different microwave power， loading and microwave power per unit load， ie microwave power density （Mp0）. The results showed that both the microwave power and loading had a significant influence on the drying time， higher the microwave power would lead to higher vacuum degree and faster the drying. In the test range， the Mp0 was the decisive factor affecting the moisture removing properties of cooked chestnut. When the Mp0 was 5.5 W/g， the optimal sensory quality score of the ripened chestnut was 80 points， high expansion， brittleness rate and porosity was 2.92， 3.51 N and 61.4% respectively， water activity was less than 0.3. Based on the relationship between moisture ratio and time of microwave moisture removing of wild chestnut， Matlab software was used to determine the drying properties accorded with Page equation， and the model fitting degree was very high （0.99）， indicating that the equation could clarify the stage of microwave moisture removing of chestnut crisp ball. The results would provide theoretical references for the drying production of chestnuts. The random distribution of the particle moisture content of the sample after drying conformed to the sixth-order polynomial distribution， which was consistent with the average moisture content， indicating that the model be used to improve the quality stability of dried C. henryi crisp balls by phased varying power microwave

Keywords： Castanea henryi crisp ball; microwave power density; moisture content after pre-drying; dynamic model

DOI： 10.3969/j.issn.1000-2561.2021.07.034

錐栗（Castanea henryi）屬殼斗科栗屬堅(jiān)果類植物，素有“干果之王”的美稱，其果實(shí)營(yíng)養(yǎng)豐富，蛋白質(zhì)含量較高，富含多種限制性氨基酸。因?yàn)樯r錐栗水分含量高，淀粉水解快，極易出現(xiàn)發(fā)芽、霉?fàn)€、蟲蛀等情況，因而貯藏保鮮難度大[1]。另外，市面上的錐栗仍多以糖炒板栗為主，不過目前其深加工制品也涌現(xiàn)了錐栗酥、錐栗飲料等系列產(chǎn)品，但仍存在安全貯藏期和品質(zhì)貨架期短的問題，亟需引進(jìn)新型干燥技術(shù)提升錐栗類深加工食品的生產(chǎn)效率和品質(zhì)。

目前，傳統(tǒng)薄層恒溫?zé)犸L(fēng)干燥法是錐栗及錐栗深加工制品干制貯藏的主要手段，而熱風(fēng)干燥屬于熱傳導(dǎo)式加熱方式。由于溫升速度慢、熱導(dǎo)效率低、介質(zhì)表面過熱，易造成產(chǎn)品品質(zhì)劣變[2-3]。所以，為探索錐栗干燥加工過程中的脫水特性，提高干制品品質(zhì)，利用新型的干燥技術(shù)建立高效簡(jiǎn)便的干制加工動(dòng)力學(xué)模型，對(duì)于指導(dǎo)和推動(dòng)錐栗精深加工產(chǎn)業(yè)化發(fā)展是十分必要的。

微波干燥作為新型干燥技術(shù)，其突出的優(yōu)勢(shì)包括脫水速度快、能量轉(zhuǎn)化率高、操作控制簡(jiǎn)單和選擇性加熱等。所以，不同農(nóng)產(chǎn)品及加工食品多采用微波技術(shù)進(jìn)行干燥保藏[4]。農(nóng)產(chǎn)品干燥的動(dòng)力學(xué)模型符合多孔物料的薄層干燥過程。毛志幸等[5]以微波干燥方式對(duì)白果進(jìn)行干燥，發(fā)現(xiàn)白果微波干燥的失水規(guī)律可以通過Page方程來描述。曹小紅等[6]對(duì)板栗片的微波干燥過程看出其干燥規(guī)律與水分?jǐn)U散系數(shù)有關(guān)。周潔等[7]使用微波干燥工藝干燥板栗原料，發(fā)現(xiàn)起初物料溫度過高會(huì)造成產(chǎn)品質(zhì)量嚴(yán)重降低，因而采用微波-熱風(fēng)聯(lián)合干燥提高其干后品質(zhì)。物料的傳熱、傳質(zhì)特性是決定其脫水干燥特性的重要影響因素[8]。栗仁是以水分、淀粉、蛋白質(zhì)和脂肪等多種營(yíng)養(yǎng)素為主的多孔物料，不同干制過程會(huì)明顯影響其感官及營(yíng)養(yǎng)質(zhì)量[9]。因此，為避免微波干燥的邊角化效應(yīng)，結(jié)合蒸煮預(yù)處理的方式將錐栗熟化碾碎后，加入適量膨化劑搓成球形，然后摸索微波強(qiáng)度和物料裝載量這2個(gè)關(guān)鍵因素對(duì)控制物料干燥特性及干后品質(zhì)的規(guī)律，但是熟化錐栗球的微波干燥機(jī)理尚不清楚。所以，充分利用微波干燥技術(shù)研究熟化錐栗不同脫水進(jìn)程的動(dòng)力學(xué)規(guī)律，能夠?yàn)閮?yōu)化干燥工藝、提高農(nóng)副產(chǎn)品質(zhì)量等方面提供重要的理論支持和技術(shù)依據(jù)。

1? 材料與方法

1.1? 材料

錐栗，購自武夷山市農(nóng)貿(mào)市場(chǎng)。

EM7KCGW3-NR美的微波爐，廣東美的電器股份有限公司;JA2003電子天平，上海舜宇恒平科學(xué)儀器有限公司;HD-E804-AT電熱恒溫鼓風(fēng)干燥箱，上海精宏儀器有限公司。

1.2? 方法

1.2.1? 錐栗球的制作及干燥? 新鮮錐栗用剝殼器開口后，去內(nèi)皮，于高壓滅菌鍋中121 ℃處理20 min，加入5%白砂糖0.5%食鹽和0.2%的NaHCO3混勻粉碎成泥[10]，冷卻搓成直徑2 cm的球。按GB 5009.3—2016直接烘干法測(cè)定熟化錐栗球的初始濕基含水率為（52.34±1.00）%，相當(dāng)于干基含水率為（109.82±1.00）%。

稱取一定裝載量的錐栗球，手動(dòng)設(shè)定不同微波強(qiáng)度條件后，將樣品單層均勻平鋪在容器中置于微波爐爐盤上進(jìn)行干燥，以每次間隔2 min時(shí)間快速取出樣品盤稱量后記錄，樣品含水率≤8.00%即為干燥終點(diǎn)。

1.2.2? 單因素試驗(yàn)? 固定錐栗裝載量為50.0 g，微波強(qiáng)度依次設(shè)定為231、385、539 W，考察微波強(qiáng)度對(duì)熟化錐栗球微波干燥特性的影響。

固定微波強(qiáng)度385 W，裝載量依次設(shè)定為50.0、70.0、90.0 g，考察裝載量對(duì)錐栗脆球微波干燥特性的影響。

1.2.3? 相關(guān)指標(biāo)的計(jì)算? 干基含水率計(jì)算公式：

Mg=（m-mg/mg）×100%? （1）

式中：Mg為干基含水率，%;mg為物料絕干質(zhì)量，g;m為物料最初質(zhì)量，g。

干基失水速率計(jì)算公式：

DR=Mt/t???? （2）

式中：DR為干基失水速率，%/min;Mt為相鄰兩次干基含水率的變化量，%;t為相鄰兩次稱重的間隔時(shí)間，min。

初始微波功率密度：

Mp0=P/m?? （3）

式中：Mp0為微波功率密度，W/g;P為微波強(qiáng)度，W;m為裝載量，g。

1.2.4? 感官品質(zhì)評(píng)定? 由經(jīng)過專業(yè)培訓(xùn)的評(píng)價(jià)員8人對(duì)干燥后的熟化錐栗球按照表1所示評(píng)分標(biāo)準(zhǔn)進(jìn)行感官品質(zhì)評(píng)價(jià)。

1.2.5? 膨化率? 采用游標(biāo)卡尺測(cè)定干燥后錐栗球的直徑，計(jì)算體積。膨化率[10]的計(jì)算公式如下：

（4）

式中：ER為膨化率;V0為膨化前體積，cm3;Vt為膨化后體積，cm3。

1.2.6? 脆度? 簡(jiǎn)易脆度測(cè)定裝置：精度0.5 g托盤天平，于天平的一端放燒杯，測(cè)定過程中以一定的流速添加水于水杯中，天平另一端放樣品，并在樣品上方固定1個(gè)直徑為2 mm的長(zhǎng)釘。樣品脆度以樣品斷裂時(shí)所受壓力（N）來計(jì)算，每1組實(shí)驗(yàn)取6個(gè)樣品測(cè)量[11]。

1.2.7? 孔隙率? 錐栗球體積用油菜籽置換量筒里錐栗球間空隙的方法[12]測(cè)定。

D0= m0 /V0????? （5）

Dt = mt /Vt （6）

N=（1–Dt/D0） ×100%???? （7）

式中：D0為初始密度，g/mL;m0為初始質(zhì)量，g;V0為初始體積，mL;Dt為表觀密度，g/mL;mt為干燥后質(zhì)量，g;Vt為干燥后體積，mL;N為孔隙率，%。

1.2.8? 水分活度測(cè)定? 根據(jù)康衛(wèi)氏皿擴(kuò)散法[13]測(cè)定。

1.2.9? 常用干燥動(dòng)力學(xué)模型? 常見的薄層干燥模型包括：

單項(xiàng)擴(kuò)展模型： MR=aexp（–kt）?? （8）

指數(shù)模型：? MR=exp（–kt）???? （9）

Page模型： MR=exp（–ktn）?? （10）

水分比：???? MR=（Mt–Me）/（M0–Me）?? （11）

式中：MR為水分比;Mt為錐栗球在t時(shí)刻瞬時(shí)含水率;Me為錐栗球在t時(shí)刻的平衡含水率，干燥過程中Me近似為0;M0為0時(shí)刻錐栗球的起始干基含水率;a、k、n為待定系數(shù)。

所以水分比的簡(jiǎn)化公式為：

MR=Mt/M0? （12）

將公式（8）～（10）取對(duì)數(shù)：

單項(xiàng)擴(kuò)展模型： ln（MR）=lna–kt?? （13）

指數(shù)模型：? ln（MR）= –kt????? （14）

Page模型： ln[–ln（MR）]=lnk+nlnt????? （15）

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

每組樣品指標(biāo)重復(fù)測(cè)定3次。采用Excel 2016軟件進(jìn)行數(shù)據(jù)處理和制圖，采用SPSS 19.0軟件進(jìn)行差異顯著性分析。利用Matlab軟件進(jìn)行干燥動(dòng)力學(xué)建模、回歸擬合和準(zhǔn)確性校驗(yàn)。

2? 結(jié)果與分析

2.1? 微波強(qiáng)度對(duì)熟化錐栗微波干燥特性的影響

微波強(qiáng)度對(duì)錐栗球脫水特性的影響見干基含水率與時(shí)間關(guān)系曲線（圖1）和干基失水速率與時(shí)間的關(guān)系曲線（圖2）。

從圖1可以看出，在微波強(qiáng)度231、385、539 W時(shí)錐栗球達(dá)到干燥終點(diǎn)的時(shí)間等比例縮短，分別是32、16、8 min。當(dāng)微波強(qiáng)度逐漸減小時(shí)，干基含水率隨干燥時(shí)間延長(zhǎng)下降速度越來越慢，由此可知微波強(qiáng)度與物料脫水效率正相關(guān)。

由圖2可知，385 W和539 W微波條件下，錐栗球脫水過程可分為升速脫水階段和降速脫水階段;而231 W微波條件下還包含恒速脫水階段。一方面微波強(qiáng)度與升速脫水階段的周期正相關(guān)，隨著微波強(qiáng)度大于385 W時(shí)恒速脫水過程甚至逐漸消失，從而出現(xiàn)轉(zhuǎn)入降速脫水過程的時(shí)間更短。另一方面，微波強(qiáng)度與升速階段脫水加速度正相關(guān)，干燥周期隨著微波強(qiáng)度增大等比縮短。在231、385、539 W微波強(qiáng)度條件下，物料對(duì)應(yīng)的干基失水速率最大值依次是4.26、8.70、10.90 %/min，呈遞增的趨勢(shì)。

2.2? 裝載量對(duì)熟化錐栗微波干燥特性的影響

錐栗球裝載量分別是50.0、70.0、90.0 g時(shí)，對(duì)干燥曲線和干燥速率曲線的影響如圖3和圖4所示。

由圖3可知，裝載量明顯影響錐栗球微波干燥過程中的干基含水率。隨著裝載量增大，干基含水率隨時(shí)間延長(zhǎng)下降速度逐漸變小，干燥時(shí)間增加。與微波強(qiáng)度對(duì)物料干燥干基含水率的影響趨勢(shì)相比，裝載量對(duì)物料干基失水速率的加速度的影響更小。這可能與單位質(zhì)量微波能量對(duì)物料溫升速度的影響有關(guān)，單位質(zhì)量錐栗球吸收微波能量后發(fā)生梯度升溫則可以有效控制水分蒸發(fā)速度，進(jìn)一步影響物料的干燥規(guī)律和干燥效率。

如圖4所示，錐栗球的裝載量同樣影響其分界干燥過程。當(dāng)裝載量50.0 g時(shí)，呈現(xiàn)升速和降速過程，而裝載量為70.0、90.0 g的干燥過程還包含了恒速階段，尤其是裝載量為70.0 g與90.0 g時(shí)，錐栗球的脫水特性趨勢(shì)相似，說明超過70.0 g后裝載量對(duì)物料干燥特性的影響不明顯。

由圖1～圖4可推算出Mp0對(duì)錐栗球脫水規(guī)律的影響，如圖5所示，錐栗球在不同Mp0區(qū)間的干燥過程存在明顯差異。Mp0越大，水分子吸收能量越多，分子在單位時(shí)間內(nèi)方向改變次數(shù)越多，錐栗球干燥速率曲線越陡。當(dāng)Mp0區(qū)間為4.3～7.7 W/g時(shí)，錐栗球干燥過程分為升速、恒速和降速3個(gè)階段;而Mp0在10.8 W/g時(shí)，則僅包含升速和降速2個(gè)階段。當(dāng)Mp0為4.52 W/g時(shí)，錐栗球干燥過程由升速向恒速階段轉(zhuǎn)換的含水率為79.76%;Mp0為5.61 W/g時(shí)，錐栗球干燥過程由升速到恒速階段轉(zhuǎn)變的含水率為85.71%。當(dāng)Mp0 7.7 W/g以下時(shí)，物料存在一個(gè)相對(duì)緩慢的升速階段，其進(jìn)入減速階段的含水率范圍是20.07%～ 24.45%;而Mp0為7.7 W/g及以上時(shí)，物料進(jìn)入減速過程的含水率在41.38%～54.45%之間。可以看出，無論是升速到恒速階段，還是由恒速到降速階段，轉(zhuǎn)換含水率都是隨Mp0的增大而增大。因此，由于升速階段，物料質(zhì)量減少，Mp0增加過速，導(dǎo)致恒速階段缺失會(huì)造成干燥產(chǎn)品質(zhì)量下降;而降速階段，隨著Mp0不斷增加，微波極易對(duì)錐栗球的非水成分產(chǎn)生較大影響。后期Mp0超過10.0 W/g時(shí)，錐栗球減速過程明顯縮短。因此，Mp0為4.3～10.8 W/g時(shí)，錐栗球微波干燥抵達(dá)恒速過程的含水率為79.76%～85.71%;轉(zhuǎn)換含水率在20.07%～54.45%時(shí)抵達(dá)干燥降速階段。把握合適的含水率進(jìn)行Mp0的調(diào)整，是保證錐栗球干燥品質(zhì)的有效手段。

2.3? 微波功率密度對(duì)錐栗球水分活度的影響

由不同初始微波功率密度下錐栗球在不同干基含水率時(shí)的水分活度，確定錐栗球的內(nèi)部水勢(shì)（圖6）。從圖6可知，相同干基含水率下，隨著微波功率密度的升高，錐栗脆球的水分活度呈逐漸減小的趨勢(shì)，但變化并不明顯，表明錐栗脆球在不同Mp0下具有較好的穩(wěn)定性;相同Mp0下，錐栗脆球的水分活度隨干基含水率的增大而增大，呈J型曲線。

2.4? 干燥后的錐栗球感官品質(zhì)和物性評(píng)價(jià)

不同的微波強(qiáng)度和裝載量條件下，錐栗球微波干燥制品的感官品質(zhì)得分如圖7、圖8所示。

由圖7和圖8可知，錐栗球的感官品質(zhì)各屬性評(píng)分均在10分以上，總得分在60～80之間，符合生產(chǎn)要求。隨著微波強(qiáng)度的增大，錐栗球的香味逐漸加強(qiáng)，但是外觀完整性下降，色澤逐漸由棕黃色向焦褐色轉(zhuǎn)變;隨著裝載量增大，錐栗球水分汽化受熱不均勻，色澤隨之加深，且有雜色斑點(diǎn)出現(xiàn)。綜合分析可知，錐栗球處于微波強(qiáng)度385 W、裝載量50.0 g（即Mp0 7.7 W/g）時(shí)的干燥產(chǎn)品感官評(píng)分最高（80分），此時(shí)干燥樣品質(zhì)量最佳。從感官特性分析結(jié)果可知，中高火（539 W）、裝載量50.0 g條件下，錐栗球發(fā)出焦香，碎裂嚴(yán)重，外形不完整，色澤呈焦褐色且炭化爆點(diǎn)較多，錐栗球的感官品質(zhì)下降嚴(yán)重。當(dāng)Mp0超過14.00 W/g（即高火700 W、裝載量50.0 g）則因劇烈脫水使樣品變成焦黑色，感官質(zhì)量驟然下降，因此數(shù)據(jù)中未體現(xiàn)此組數(shù)據(jù)。而前期Mp0太低造成物料的恒速脫水過程過長(zhǎng)，進(jìn)而干燥效率降低，且錐栗脆球香氣不足，可能是微波能量不足導(dǎo)致干燥時(shí)間延長(zhǎng)，致使呈香基團(tuán)減少造成的。

由圖9可知，隨著Mp0增大，錐栗脆球的感官品質(zhì)發(fā)生明顯變化，尤其是形態(tài)完整性出現(xiàn)先上升后下降的趨勢(shì)。膨化率能夠反映干燥制品的質(zhì)地疏松程度;脆度可作為評(píng)判干燥制品口感改良程度的一個(gè)重要指標(biāo);孔隙率是指固體材料中孔隙體積與總體積的百分比，與物料的密實(shí)度呈反比，是用于描述多孔介質(zhì)滲流阻力的一個(gè)重要參數(shù)。在干燥過程中，錐栗脆球硬度變大往往與失水速率正相關(guān);而物料的膨化性能越弱，其脆度則越低，孔隙率越小，則水分遷移阻力越大。根據(jù)膨化率、脆度和孔隙率來考察不同Mp0對(duì)錐栗脆球質(zhì)地結(jié)構(gòu)的改良程度，結(jié)果見表2。從表2可以看出，當(dāng)Mp0為4.6～5.5 W/g，錐栗脆球孔隙率比較大，而膨化率相對(duì)較小，其中當(dāng)功率是5.5 W/g時(shí)，膨化率最大，孔隙率最大，且脆度較好，水分活度低于0.3，產(chǎn)品形態(tài)感官評(píng)分達(dá)到最高（24.6）且綜合評(píng)分為80.9。

2.5? 錐栗脆球微波干燥動(dòng)力學(xué)模型

2.5.1? 干燥模型的確定? 微波干燥的動(dòng)力學(xué)模型能夠解析各干燥參數(shù)之間的變化規(guī)律，可以在生產(chǎn)、加工中發(fā)揮切實(shí)有效的指導(dǎo)作用。在不同微波強(qiáng)度及裝載量條件下，按照實(shí)驗(yàn)數(shù)據(jù)繪制以–ln（MR）為縱坐標(biāo)與以時(shí)間t為橫坐標(biāo)的關(guān)系曲線（圖10A和圖10B）、以ln[–ln（MR）]為縱坐標(biāo)與以lnt為橫坐標(biāo)的關(guān)系曲線（圖11A和圖11B）。

由圖10可知，–ln（MR）與t的關(guān)系曲線表現(xiàn)為冪增長(zhǎng)趨勢(shì);從圖11中可以看出，ln[–ln（MR）]與lnt存在明顯的線性相關(guān)關(guān)系。由此可知，錐栗脆球的微波干燥動(dòng)力學(xué)規(guī)律可以用MR=exp（–ktn）的Page方程來預(yù)測(cè)，然后，使用Matlab確定其模型擬合度。

2.5.2? 建立Page模型和統(tǒng)計(jì)分析? 利用Matlab軟件對(duì)lnt與ln[–ln（MR）]之間的變量關(guān)系進(jìn)行線性回歸擬合分析，代入公式ln[–ln（MR）]=lnk+nlnt中，計(jì)算得出其中的2個(gè)模型常數(shù)lnk和n（表3）。

由表3可知，統(tǒng)計(jì)量F檢驗(yàn)依次為2578.5910、2051.9730、3148.9770、1459.4490、2800.4870，每個(gè)對(duì)應(yīng)P值都是0，說明自變量與因變量之間在回歸模型中存在極顯著線性相關(guān)關(guān)系。因此，錐栗脆球微波干燥的動(dòng)力學(xué)數(shù)學(xué)模型可根據(jù)各擬 合方程來確立。

根據(jù)表3把lnk、n當(dāng)作因變量，Mp當(dāng)作自變量，利用Matlab軟件對(duì)lnk、n與Mp之間的變量關(guān)系進(jìn)行非線性擬合，計(jì)算求出錐栗脆球微波干燥動(dòng)力學(xué)模型回歸擬合結(jié)果（表4）。

從表4可以看出，回歸模型的擬合程度很高（R2>0.97），說明該回歸模型適合描述錐栗脆球微波干燥的脫水過程。由表3可求得待定系數(shù)，代入Page方程：ln[–ln（MR）]=lnk+nlnt中，可知：

（12）

（13）

2.5.3? 干燥模型的驗(yàn)證? 錐栗脆球干燥動(dòng)力學(xué)模型準(zhǔn)確性的驗(yàn)證可以任選一組試驗(yàn)數(shù)據(jù)來完成。本研究選擇微波強(qiáng)度385 W、裝載量50.0 g的參數(shù)進(jìn)行校驗(yàn)，模型方程為ln[–ln（MR）]= –5.7699+ 2.3512lnt，根據(jù)水分比MR的與預(yù)測(cè)值和實(shí)驗(yàn)值的線性回歸關(guān)系分析模型預(yù)測(cè)的準(zhǔn)確性（圖12）。從圖12可以看出實(shí)驗(yàn)值和預(yù)測(cè)值基本處于同一直線上，表明錐栗脆球的微波干燥過程中水分的變化規(guī)律可以通過干燥動(dòng)力學(xué)模型Page方程來準(zhǔn)確描述。

圖13是5.5 W/g條件下微波干燥樣品平均含水率為9.00%的每個(gè)錐栗球含水率隨機(jī)分布情況。從圖13中可以看出，干燥后的錐栗球含水率分布較為集中，平均含水率為9%的同批次錐栗脆球的水分分布從5.00%～13.00%，近似服從于多項(xiàng)式分布，且其決定系數(shù)達(dá)到0.993，說明此水分分布符合六次多項(xiàng)式隨機(jī)分布。

3? 討論

錐栗脆球的干燥過程包括升速、恒速和減速3個(gè)階段。其中，錐栗脆球升速脫水階段是干基含水率從0升至最大失水速率的過程，物料在脫水升速階段的熱能積累來源于微波能量聚集轉(zhuǎn)化，由于熟化錐栗初始水分含量很高，物料內(nèi)部水勢(shì)較大，內(nèi)部水分極易受熱變?yōu)樗魵馍l(fā)，因此干燥初期，限制內(nèi)部水分遷移速度是確保錐栗初始升速階段產(chǎn)品疏松多孔質(zhì)地的關(guān)鍵。緊接著由于物料水分部分散失后，錐栗球水分含量下降，雖然熟化后組織破潰對(duì)于水分的截留能力下降，但是水分活度降低增大了內(nèi)部水分穩(wěn)定性，造成物料在恒速階段干基失水速率增大到最大值。最后，物料的大量水分散失后碳酸氫鈉發(fā)生熱解使得錐栗球質(zhì)地結(jié)構(gòu)疏松，水分活度越來越低，使得內(nèi)外水勢(shì)差越來越小，造成物料失水速率逐漸減小。雖然水分穩(wěn)定性的增加與高微波功率直接相關(guān)，但強(qiáng)迫性脫水會(huì)造成脫水過度，進(jìn)而導(dǎo)致物料干燥品質(zhì)下降，曹小紅等[6] 利用微波干燥板栗片的試驗(yàn)發(fā)現(xiàn)其脫水過程以降速階段為主，熟化錐栗微波干燥過程則存在明顯的升速過程。說明微波干燥除了與物料成分和微波條件相關(guān)外，也會(huì)受到物料外形和薄層厚度的明顯影響。另外，多孔介質(zhì)干燥層的孔隙率對(duì)干燥速率有顯著影響[14-15]。錐栗球在微波干燥的過程中，樣品的內(nèi)部水勢(shì)高于環(huán)境的外部水勢(shì)，水分由樣品遷移到環(huán)境中去，當(dāng)內(nèi)部的孔隙率較大時(shí)，水分遷移的物理阻隔較小，樣品的內(nèi)部水勢(shì)較高，形成較大的內(nèi)外水勢(shì)差，水分遷移速率增加，可以有效分散極性分子水的能量聚集效應(yīng)，進(jìn)而控制水分散失速度，故微波干燥時(shí)，生栗裝載量越大，雖然吸收能量越大，但過大的裝載量會(huì)導(dǎo)致單位質(zhì)量物料能量不足，增加干燥時(shí)間，干燥效率相繼下降。由此可知，錐栗球的水分穩(wěn)定性和物性特征是決定其微波干燥過程中水分散失速度的關(guān)鍵因素。

根據(jù)微波功率與裝載量的比值（即微波功率密度）Mp0對(duì)錐栗脆球品質(zhì)的影響規(guī)律，錐栗脆球的分段微波干燥過程和干燥品質(zhì)與Mp0和轉(zhuǎn)換含水率密切相關(guān)。當(dāng)Mp0大于7.7 W/g，恒速干燥階段變短甚至消失。Mp為4.3～10.8 W/g時(shí)，錐栗脆球含水率79.76%～85.71%時(shí)抵達(dá)恒速階段;其含水率在20.07%～54.45%時(shí)抵達(dá)降速階段。在Mp0低于5.5 W/g時(shí)，干燥時(shí)間過長(zhǎng)，產(chǎn)品中非水組分發(fā)生變化，斷裂力增大，脆度下降;且失水過程中，產(chǎn)品中蒸發(fā)掉的水分空間還在，因而錐栗脆球仍具有多孔結(jié)構(gòu)，質(zhì)地疏松性較好;微波干燥功率過高（10.8 W/g）時(shí)，由于物料急劇失水，呈香物質(zhì)流失[15]，使纖維結(jié)構(gòu)明顯并導(dǎo)致焦糊現(xiàn)象發(fā)生，雖然膨化率較大，但產(chǎn)品易碎，因此，錐栗脆球感官品質(zhì)也更差[16]。因此，控制錐栗脆球Mp0在5.5 W/g結(jié)合干燥動(dòng)力學(xué)特征，進(jìn)而控制干燥失水速率在掠奪式脫水速率（10.9%/min）以內(nèi)，對(duì)于調(diào)節(jié)物料的溫升速度，平衡熱敏性物料的能量聚集規(guī)律，提升其質(zhì)地結(jié)構(gòu)和干燥品質(zhì)至關(guān)重要。通過Matlab軟件考察微波干燥條件對(duì)錐栗脆球脫水規(guī)律的薄層干燥動(dòng)力學(xué)模型發(fā)現(xiàn)，Page方程適合闡述錐栗脆球微波干燥動(dòng)力學(xué)規(guī)律，相關(guān)系數(shù)0.99，最終建立的動(dòng)力學(xué)方程： 且實(shí)驗(yàn)值與預(yù)測(cè)值線性擬合度較高，錐栗脆球在微波干燥過程中水分變化的規(guī)律的描述可以采用此方程進(jìn)行預(yù)測(cè)。結(jié)合感官品質(zhì)、膨化率、脆度和孔隙率評(píng)價(jià)結(jié)果發(fā)現(xiàn)，Mp0為5.5 W/g時(shí)，錐栗脆球干制品品質(zhì)最佳。當(dāng)前大部分干燥動(dòng)力學(xué)研究提及到的物料含水率均為樣品的平均含水率，但在實(shí)際過程中各個(gè)薄層物料的含水率之間都存在著一定的差距。由于隨著干燥時(shí)間延長(zhǎng)，物料溫度的升高，導(dǎo)致食品中水分子的動(dòng)能增加，失水速率增加，導(dǎo)致錐栗脆球?qū)λ肿拥奈侥芰ο陆?，?dǎo)致干基含水率降低[17-20]。因此本研究建立干燥過程中物料的水分隨機(jī)分布圖可以有效反映多孔物料顆粒內(nèi)部的水分分布情況。研究結(jié)果證明，同一批次的物料顆粒含水率與平均含水率一致。根據(jù)熟化錐栗脆球的干燥動(dòng)力學(xué)規(guī)律，后期可通過分段變功率控制物料失水速率和顆粒內(nèi)部孔隙率的方式提高錐栗球微波干燥效率，縮短干燥時(shí)間及提升干燥品質(zhì)。

參考文獻(xiàn)

[1]? 魏宗烽， 邵? 穎， 魏明奎. 板栗深加工產(chǎn)品的研究現(xiàn)狀及前景分析[J]. 農(nóng)產(chǎn)品加工（學(xué)刊）， 2014， 15（8）： 69-71.

[2]? Silva A C D， Sarturi H J， DallOglio E L， et al. Microwave drying and disinfestation of brazilnut seeds[J]. Food Control， 2016， 70（12）： 119-129.

[3]? Cuccurullo G， Giordano L， Albanese D， et al. Infrared thermography assisted control for apples microwave drying[J]. Journal of Food Engineering， 2012， 112（4）： 319-325.

[4]? Vandeweyer D， Lenaerts S， Callens A， et al. Effect of blanching followed by refrigerated storage or industrial microwave drying on the microbial load of yellow mealworm larvae （Tenebrio molitor）[J]. Food Control， 2017， 71（1）： 311-314.

[5]? 毛志幸， 孫? 輝， 陳宗道. 白果微波干燥特性及干燥動(dòng)力學(xué)模型研究[J]. 食品工業(yè)科技， 2017， 38（22）： 18-23， 28.

[6]? 曹小紅， 常學(xué)東. 板栗的微波干燥特性及其對(duì)干后品質(zhì)的影響[J]. 食品工業(yè)科技， 2005， 26（1）： 63-65.

[7] 周? 潔， 孔曉玲， 高? 楠. 微波與熱風(fēng)組合干燥對(duì)板栗干燥特性的影響[J]. 包裝與食品機(jī)械， 2009， 27（2）： 11-13.

[8]? 應(yīng)巧玲， 勵(lì)建榮， 傅玉穎， 等. 食品薄層干燥技術(shù)的研究進(jìn)展[J]. 中國糧油學(xué)報(bào)， 2010， 25（5）： 115-119， 128.

[9]? 張? 樂， 趙守渙， 王趙改， 等. 板栗微波真空干燥特性及干燥工藝研究[J]. 食品與機(jī)械， 2018， 34（4）： 206-210.

[10]????? 常學(xué)東， 朱京濤， 高海生， 等. 微波膨化板栗脆片的工藝研究[J]. 食品與發(fā)酵工業(yè)， 2006， 32（9）： 78-81.

[11]????? 蔣予箭， 蔣家新， 吳彩珍. 微波膨化米餅生產(chǎn)工藝條件的研究[J]. 糧食與飼料工業(yè)， 2001（5）： 44-45.

[12]????? 田? 方， 鄭志堅(jiān)， 張仁杰， 等. 預(yù)處理對(duì)微波膨化芋頭片斷裂力和表觀密度的影響[J]. 食品工業(yè)， 2012， 33（12）： 50-51.

[13]????? 趙世民. 康衛(wèi)氏皿擴(kuò)散法測(cè)定蕎麥粉的水分活度[J]. 安徽農(nóng)業(yè)科學(xué)，2012， 40（5）： 2901-2902， 2905.

[14]????? 姚? 荷， 譚亦成， 譚興和， 等. 微波-熱風(fēng)聯(lián)用制取筍干工藝條件優(yōu)化[J]. 食品科學(xué)， 2019， 40（12）： 260-266.

[15]????? 宋瑞凱， 張付杰， 楊? 薇， 等. 馬鈴薯微波干燥動(dòng)力學(xué)建模與仿真[J]. 湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2018， 44（2）： 204-209.

[16]????? Prommas P， Rattandecho P， Cholaseuk D. Energy and exergy analyses in drying process of porous media using hot air[J]. International Communications in Heat and Mass Transfer， 2010， 37（4）： 372-378.

[17]????? Roca E，Guillar D V， Gulber S， et al. Moisture migration in a cereal composite food at high water activity： Effects of initial porosity and fat content[J]. Journal of Cereal Science， 2006， 43（2）： 144 -151.

[18]????? 張鐘元， 李麗娟， 李大婧， 等. 蓮藕片薄層真空微波干燥特性及動(dòng)力學(xué)模型[J]. 江蘇農(nóng)業(yè)科學(xué)， 2019， 47（2）： 201-205.

[19]????? 周鳴謙， 劉春泉， 李大婧. 不同干燥方式對(duì)蓮子品質(zhì)的影響[J]. 食品科學(xué)， 2016， 37（9）： 98-104.

[20]????? Figen K E， Gedik A. Sorption isotherms and isosteric heat of sorption for grapes， apricots， apples and potatoes[J]. Food Science and Technology， 2004， 37（4）： 429-438.

責(zé)任編輯：崔麗虹