糯米粉的水分等溫解吸及孔特性

張 華，趙學偉，白艷紅，段瑞謙

（1.鄭州輕工業學院食品與生物工程學院，河南 鄭州 450002；2.河南省食品生產與安全協同創新中心，河南 鄭州 450002）

糯米粉的水分等溫解吸及孔特性

張 華1,2，趙學偉1,2，白艷紅1,2，段瑞謙1

（1.鄭州輕工業學院食品與生物工程學院，河南 鄭州 450002；2.河南省食品生產與安全協同創新中心，河南 鄭州 450002）

糯米粉的等溫解吸特性一方面對分析其與周圍環境之間的水分傳遞十分必要，另一方面，還可用于計算糯米粉內孔的特性，借此加深對水分吸附機理的了解。本實驗采用靜態稱質量法在10、20、30 ℃條件下測定糯米粉在10 個水分活度解吸后的平衡水分含量，然后采用4 個等溫吸附模型對實驗結果進行擬合，并分析糯米粉的孔特性。研究結果表明：糯米粉中水分的解吸屬于典型的Ⅱ型吸附；Lewicki模型最能描述同一溫度條件下的等溫解吸過程，而GDW模型能夠同時描述溫度和水分活度對平衡水分含量的影響。糯米粉中同時存在微孔和介孔；孔徑分布為單態分布，僅在微孔區域出現一單峰，溫度變化僅對該峰的峰值產生影響；微孔體積隨溫度降低而增大；糯米粉的吸附面具有分形特性，隨溫度降低吸附表面的多孔性增強，單位吸附面積增大，最終導致吸附能力增強。

糯米粉；解吸；孔徑；孔徑分布；單位表面積；微孔體積；分形維數

食品及其組分的水分吸附特性（包括吸濕和解吸）對涉及水分傳遞的食品加工、保藏過程十分重要。如對于最為常見的干燥過程，水分吸附特性將決定干燥時間、成品的貨架期以及不同水分活度的組分混合后的平衡水分[1]。水分吸附特性對于控制不期望的化學反應、酶反應[2]以及多域（multi-domain）食品內的水分轉移[3]十分必要。水分等溫吸附數據可用來計算吸附過程中的熱力學特性（包括微分焓和微分熵、積分焓和積分熵、焓-熵補償、鋪展壓等），從而更為深刻地理解等溫吸附現象和吸附機理[4-5]。

等溫吸附數據還可用于計算多孔性固體和粉體的孔特性，包括單位吸附面積、孔徑及孔徑分布、孔體積、吸附表面的分型維數等[6]。很多物理特性（如：吸附、溶解、分散等）的平衡過程和速率與表面結構有很大關系[7]。近年來有研究顯示增大微孔部分的體積，可以保證食品雖有較高水分含量但卻不發生品質劣變[8-9]。而分形維數則能提供吸附質與吸附劑之間相互作用的物理機理（是范德華作用或是毛細管凝結）[10]。

氣體吸附法是測定多孔性材料孔特性的簡單而有效的方法，氣態N2是最為常見的吸附質[6]。近年來有研究者根據食品對水蒸氣的吸附結果，分析計算食品的吸附面積、孔徑分布和微孔體積[8-9,11-12]。

糯米幾乎不含直鏈淀粉（0%～2%）[13]。糯米粉是由糯米用水浸泡后碾磨而成。許多以糯米粉為主料的食品加工主要是利用蒸煮后其柔軟且黏度高的特性。糯米粉在我國食品工業中的一個重要應用是制作湯圓的皮層。速凍是實現湯圓這一傳統食品工業化的重要技術手段。糯米粉的水分解吸特性決定湯圓速凍及保存過程中的失水量，從而影響其最終品質如凍裂率，也將決定湯圓皮與湯圓芯之間水分分配比。因此，增進對糯米粉水分解吸特性的了解是對湯圓速凍進行過程分析和品質控制的基礎。

前人研究建立了米淀粉、米粉的水分吸附等溫方程并計算出凈吸附熱[14-16]，而關于糯米粉的等溫吸附特性鮮見報道。本研究的目為：1）根據糯米粉在不同溫度、不同水分活度下的平衡水分測定結果，確定描述糯米粉等溫解吸的最適模型；2）計算糯米粉內孔的特性，包括吸附面積、孔徑及其分布、微孔體積和分形維數，以增進對糯米粉的孔結構特性和水分吸附機理的了解。

1 材料與方法

1.1 材料與試劑

糯米粉購自鄭州市華聯超市，經測定其組分含量為：蛋白質8.08%（以干基計，下同）、粗脂肪1.11%、粗纖維0.92%、灰分0.90%、碳水化合物88.99%（差減法）。將糯米粉與適量的水混合，使最終的水分含量為45%（以濕基計），調質3 h后用于解吸實驗。

1.2 水分解吸實驗設計

采用靜態稱質量法在10 個水分活度（water activity，aw）條件下測定糯米粉的平衡水分含量（Me）。將約5 g的糯米粉放入盛有不同飽和鹽溶液的干燥器中。所用的鹽溶液包括：LiBr、LiCl、CH3COOK、MgCl2、K2CO3、Mg(NO3)2、KI、NaCl、KCl和K2SO4，不同溫度條件下所對應的水分活度參見文獻[17]。在aw＞0.75時，干燥器內放一盛有麝香草粉的小試管，以防止實驗期間微生物在糯米粉上生長。然后將干燥器放入恒溫培養箱中，溫度分別控制在10、20、30 ℃（波動范圍±1 ℃），平衡水分一直到樣品質量不發生明顯變化 （±0.001 g）。平衡后的樣品采用105 ℃烘箱法測定水分含量（kg/kg（以干基計））。所有實驗進行3 次，取其平均值。

2 結果與分析

2.1 糯米粉的水分等溫解吸

圖1 糯米粉的解吸平衡水分及Lewicki模型的擬合曲線Fig. 1 Equilibrium moisture contents of glutinous rice fl our and curves fi tted with the Lewicki model

糯米粉在10、20、30 ℃不同水分活度條件下的平衡水分見圖1，其中的實線是根據最適模型給出的預測值?？梢钥闯?，糯米粉的等溫解吸曲線呈Sigma形，根據Brunauer等[18]提出的分類方法，屬Ⅱ型吸附等溫線，這是以有機大分子為主要成分的食品的典型類型?？梢愿鶕胶馑蛛S水分活度變化的特點將平衡曲線分為3 個區域[19]：在aw較低（＜0.07）時，水分與固體物質之間的結合十分緊密，這部分水包括結構水和單分子層水分；在aw為中等水平（約0.07～0.80）時，水分子與固形物的結合作用相對較弱，這部分水分通常存在于細小的毛細管內；在高aw區域（＞0.8），水分與固體的結合十分微弱，其特性接近自由水，主要存在于固體內的孔洞、縫隙和大的毛細管內。由圖1還可明顯看出，在相同水分活度條件下，平衡水分含量隨溫度降低而增大。通常認為是高溫引起的物理或化學變化導致了可吸附位點數量降低所致，降低的幅度取決于食品組分的自身特性[20]。

描述水分活度與其對應的平衡水分之間數量關系的方程稱為吸附等溫式。本研究采用Chung-Pfrost模型（式（1））、Oswin模型（式（2））、Lewicki模型（式（3））和GDW模型（式（4））模擬同一溫度條件下的平衡水分含量與水分活度之間的關系。其中Chung-Pfrost模型和Oswin模型為美國農業與生物工程師協會所推薦[21]。Lewicki模型是Lewicki于1998年提出的，該模型的一個特點是aw增大至1時平衡水分含量趨向無窮大[22]。GDW模型是Furmaniak等于2005年提出的，并將其成功用于描述食品的水分吸附特性[23]。

式中：Me為平衡水分含量/（kg/kg）；aw為水分活度；a、b、c、f為模型參數。

本研究采用決定系數（R2）、標準估計誤差（SEE）、平均相對偏差（mean relative deviation，MRD）為統計分析指標評價上述4 個模型對糯米粉解吸平衡水分的預測性能。各模型的擬合參數及統計分析指標值見表1?？梢钥闯?，就糯米粉而言，在本研究所涉及的溫度和水分活度范圍內，Lewicki模型的預測性能最好，因其R2最大而SEE、MRD值最??；與Chung-Pfrost模型相比，GDW模型的MRD值較大但SEE值較??；Oswin模型的SEE、MRD值比GDW模型的大。

表1 4 種模型擬合所得的模型參數及統計分析指標Table 1 Statistical indicators and best-f i t parameters of four models fi tted to isotherm data of glutinous rice fl our at 10, 20 and 30 ℃

一個好的模型不僅要具有好的統計學指標，其參數還應有合理的物理意義。上述4 個模型中，GDW模型為理論模型，該模型的理論假設為：存在初級吸附位點，其吸附符合朗繆爾吸附；部分吸附在初級位點的水分子可以作為二級吸附位點；二級位點上的吸附遵循Dubinin-Serpinsky假設[23]。GDW模型中，a為吸附在初級吸附位點上的最大吸附量，f為可提供二級吸附位點的數量與初級位點總量的比；b、c為初級、二級吸附的動力學常數，其與溫度的關系為：

式中：b0、c0分別為b、c的指前因子；ΔEb、ΔEc為一級、二級吸附焓；R為氣體常數（8.314×10-3kJ/（mol·K））；T為溫度/K。

很多食品加工過程中的水分傳遞為非穩態過程，要對該過程進行分析計算，需要吸附等溫方程必須能體現溫度的影響。式（4）與式（5）、（6）聯合可實現這一點。式（5）、（6）代入式（4）后擬合實驗數據，得b0=3.23×10-5、c0=0.91，ΔEb=28.92 kJ/mol、ΔEc=0.06 kJ/mol。ΔEc接近于零，說明解析二級位點上的水分幾乎不需要吸收熱量。模型預測結果及測定結果的比較見圖2，圖中實線為y=x?？梢钥闯鍪剑?）聯合式（5）、（6）可以很好預測一定溫度、水分活度條件下糯米粉的平衡水分含量。

圖2 糯米粉平衡水分含量測定值與預測值Fig. 2 Measured equilibrium moisture contents of glutinous rice fl our at three temperatures and their predicted values

2.2 孔徑

Kelvin方程（式（7））是由水分活度計算孔徑的經典方法[11]。

式中：rc為臨界半徑/m，為水分在圓柱形孔內發生凝結（對吸濕）或蒸發（對解吸）時孔的半徑；γ為吸附質的表面張力（γ水=75.64×10-3N/m）；VM為液態吸附質的摩爾體積（對于水VM=1.8×10-5m3/mol）。

Kelvin方程沒有考慮在水分凝結之前多孔表面上已吸附的水分子層厚度。食品研究中，多數文獻采用式計算該吸附層的厚度t（×10-10m）[25-27]，該式實際上是Halsey方程針對氮氣作為吸附質的結果。Halsey方程的一般式為[6]：

式中：m為單分子吸附層的厚度（×10-10m），等于單分子層排列時的摩爾體積與面積之比[28]。以水為吸附質時，取截面積=10.6×10-20m2[29]，計算出m=2.82；A、r為模型系數。由于該方程是用于描述多層吸附的，采用式（8）擬合3 個溫度條件下的5 個較高水分活度的實驗結果，得A =-4.87，r =0.299 5。則對于水分作為吸附質時，應采用式（9）計算多分子層厚度（t）。

孔徑（rp）為臨界半徑（rc）與多層厚度（t）之和：

圖3 糯米粉內孔徑與平衡水分含量的關系Fig. 3 Radius of the pores in glutinous rice fl our as a function of equilibrium water content

aw在0.05～0.95范圍內，以0.05為間隔取各aw點，采用式（10）計算各水分活度時糯米粉內的孔半徑，采用Lewick模型計算相應水分活度條件下的平衡水分含量。平衡水分含量與孔徑之間的關系見圖3。在本研究所涉及的溫度和水分含量范圍內，糯米粉的平均孔徑在0.8～23.8 nm之間。國際理論與應用化學聯合會將孔分為3 大類：1）直徑大于50 nm的稱為大孔；2）直徑在2～50 nm之間的稱為介孔；3）直徑小于2 nm的稱為微孔[27]。根據該分類標準，糯米粉中的孔為微孔和介孔。

相同溫度條件下，隨解吸進行，平均孔徑越來越??；在水分低于某臨界值（與溫度有關）時，孔徑減小的速率降低（圖3）；進一步觀察發現該臨界水分含量對應的aw約為0.65。對于解吸過程，曲線出現折點可能意味著毛細管蒸發作用已經結束[25]。比較圖中的rp曲線和t曲線（注意t曲線的縱坐標刻度）可以看出，t對rp的貢獻較小，尤其在水分含量很低時。rp曲線的走勢主要是由rc決定的。

圖3還顯示，隨溫度升高孔徑會變大。在發生毛細管蒸發作用的區域內，溫度變化對孔徑的影響較大，這可能是固體材料受熱膨脹的微觀結果。Singh等[25-26]關于雞肉香腸和羊肉的研究也有類似發現。

2.3 孔徑分布

孔徑分布（PSD）有多種表示方法，最為常見的是孔體積V相對于孔徑rp尺寸的分布，由式（11）定義[6]。

由式（10）推導出以rp為自變量的aw的表達式，將其代入式（3），同時根據水的密度，求出以rp為自變量的水分體積V的表達式。求導即得dV/drp與rp的關系。作圖得孔徑分布曲線，見圖4。

圖4 糯米粉在3個溫度水平下的體積孔徑分布Fig. 4 Pore size distribution of glutinous rice fl our at three different temperatures

圖4 顯示3 個溫度條件下糯米粉的孔徑分布均為單態分布，僅在約0.6～2.0 nm范圍內有一個峰，而在更高rp水平時PSD值一直維持在較低水平。峰值對應的孔徑幾乎不隨溫度而變化。不同溫度條件下的峰寬也基本不變，說明溫度變化對糯米粉內孔結構的均勻性基本沒有影響。溫度越低峰值越高（圖4），說明吸水能力越強，換言之，糯米粉在較低溫度時有更強的水分吸附能力。而這種吸附能力的不同主要是微孔部分不同造成的，因為不同溫度條件下孔徑分布的不同主要體現在微孔部分。

Yano等[7]以N2為吸附氣體測定馬鈴薯淀粉的孔徑分布，發現孔徑分布隨孔徑增大而快速降低，在孔徑大于10 nm時變化很小，這種變化趨勢與本研究結果有一定的相似性。不同溫度條件下小扁豆的孔分布形態、峰寬保持不變，而馬鈴薯的峰寬隨溫度升高而增大[11]。Karathanos等[30]采用水銀孔隙儀對支鏈淀粉的測定結果表明，淀粉中存在3 類大小的孔，并指出孔徑為6～8 μm的孔可能對應于淀粉粒聚集體之間的孔隙；第2個峰位于l.0～3.5 μm處，可能對應于淀粉粒之間的孔隙，第3個峰位于約3 nm處，大多數存在于淀粉粒內部。本研究沒有在更高的孔徑范圍（rp＞2 nm）內發現有峰，說明對糯米粉來說，孔徑分布在該范圍內是均勻的。

2.4 微孔體積

通常采用Dubinin-Radushkevich方程（D-R方程，式（12））描述微孔吸附過程中被吸附物質的量與水分活度之間的關系[8-9]：

式中：n為被吸附的水分的質量（g /100 g（以干基計，下同））；n0為與微孔體積相對應的被吸附水分的質量（g/100 g）；B為一個與吸附劑微孔結構有關的常數。根據該方程，在微孔吸附范圍內lgn與lg2（1/aw）應當呈線性關系。

圖5 根據Dubinin-Radushkevich方程分析糯米粉的解析過程Fig. 5 Dubinin-Radushkevich graphs for water desorption from glutinous rice fl our

圖5 給出了以D-R方程描述糯米粉等溫解吸過程的情況?？梢钥闯?，對于糯米粉來說，在水分活度較低時，lgn與lg2（1/aw）之間呈線性關系，D-R方程成立；在水分活度較高時，曲線開始向上彎曲，說明可以利用的微孔體積已經被完全填充[31]，D-R關系式不再適用。為此采用水分活度較低的3 個點，進行線性回歸，由回歸直線的截距求出微孔內水的質量，分別是11.70、10.19、9.03 g/100 g。再根據水的密度ρ（10、20、30 ℃時分別為999.77、998.49、995.71 kg/m3[32]）計算出微孔體積。

在10、20、30 ℃時糯米粉內的微孔體積分別是11.70、10.21、9.05 mL/100 g。微孔內的水分含量比根據BET方程（式（13））計算出的單分子層水分含量（分別為7.60、6.85、6.28 g/100 g）要高，Azuara[8]和Viganó[33]等也有類似發現。實際上，由于微孔的直徑比水分子的直徑大得多，水分子在微孔內的吸附除第一層的單層吸附外，還會進行多層吸附。盡管微孔體積對應的水分含量與單分子層水分含量是兩個不同的概念，仍然可以推測，溫度變化對微孔體積的影響與對單分子層水分含量的影響應歸于相同的某個原因。

式中：M0為BET單分子層水分含量/（g/g（以干基計））；C為BET常數。

2.5 單位吸附面積

最常用的估算單位吸附面積的方法是BET法。先采用BET方程計算出單分子吸附層水分，然后由式（14）計算[6]，這里用SBET表示。

式中：MW為水的摩爾分子質量（18 g/mol）；N0為阿伏伽德羅常數（6.02×1023/mol）；AW為水分子的截面積（10.6×10-20m2）。

除BET方法外，在食品類文獻中還經常采用另外兩個方法計算單位吸附面積。其中一個是基于Caurie在1981年提出的等溫吸附模型（式（15）），由式（16）計算[34]，這里以SC1981表示。

式中：CC為Caurie常數，在數值上等于單分子吸附層內被吸附的水的密度；MC0為根據式（16）求得的單分子層水分；d為水分子的直徑，取值3.673×10－10m。以Sl表示ln（1/Me）對ln（（1－aw）/aw）所得直線的斜率，則：式（16）轉變為簡潔的式（17）。該式更為常用[25-26]。

另一個方法是Caurie于2012年[12]根據其2005年[35]提出的等溫吸附方程（式（18））推導出來的。這里以SC2012表示。以水為吸附劑時，由式（19）計算。

式（19）中：n為方程（18）中的參數，可由實驗測得的Me-aw數據組，通過模型擬合求得；CC’為模型參數。Caurie針對多種材料的研究顯示：分別以水、N2為吸附劑時采用BET方法計算出的吸附面積相差很大，不在一個數量級；而采用該作者提出的方法計算時，兩個吸附面積在一個數量級上[35]。

本研究也采用式（17）、（19）計算糯米粉的吸附面積，結果見圖6。

圖6 采用3 種方法的糯米粉內的單位吸附表面積Fig. 6 Specif i c surface area of glutinous fl our at three different temperatures calculated using three different methods

糯米粉的SEBT（222.58～269.29 m2/g）低于馬鈴薯淀粉的（287～320 m2/g）[36]以及玉米淀粉的（283.60 m2/g）[35]。圖6顯示，隨溫度升高，不論采用哪種計算方法，糯米粉內可供水分吸附的面積都減小。式（18）中的n=M0/CC’，而CC’可認為表示單分子吸附層的密度[35]，同時結合式（15）、（17）可以看出，單位吸附面積隨溫度升高而降低主要是由于單分子吸附層水分含量降低。SBET計算的是被吸附分子所占據的實際面積之和，其值大小取決于有效吸附位點的數量?？梢哉J為，溫度升高引起的某種理化變化導致一級吸附位點失活。也可從熱力學的角度理解，前面的結果顯示一級吸附的活化能為正值，這說明解吸是一個吸熱過程；因此，低溫條件下不利于解吸，使得單分子吸附層水分含量較高。

2.6 吸附面的分形維數

可由多種方法根據等溫吸附數據計算吸附面的分形維數，其中FHH模型因適用于多種情況且簡單方便而最為常用。FHH理論是由Frenkel、Halsey和Hill提出的，用于描述氣體分子在分形表面多層吸附的一個模型。根據該模型，被吸附水的質量與水分活度之間存在如下關系[37]：

式中：AM0為一個與單分子層吸附水分含量有關的常數；1/m為一個與分形維數有關的指數。根據式（20）中1/m的值可區分吸附質與吸附劑之間的作用是范德華引力或是毛細管凝結作用[38]。如果1/m＜1/3，是范德華引力，分形維數為：如果1/m＞1/3，則是毛細管凝結作用，分形維數為：

需要注意的是，根據上述方法計算出的分形維數為吸附后形成的吸附膜表面的分形維數，而不是固體表面的分形維[6]。

圖7 糯米粉的等溫吸附過程的分形分析Fig. 7 Fractal analysis of water adsorption isotherm of glutinous rice fl our at three different temperatures

根據糯米粉的等溫解吸實驗結果，lnMe對ln（－lnaw）作圖，結果見圖7?？梢钥闯?，aw值較高時兩者之間呈線性關系。Neimark[37]指出，FHH方程適用于毛細管凝結作用占主導的微孔吸附或解吸，并建議在aw＞0.7范圍內使用。在此aw范圍內進行線性回歸（圖7），根據回歸直線的斜率求出糯米粉在10、20、30 ℃時的分形維數分別為2.19、2.11、2.05。分形維數界于2～3之間, 說明糯米粉的吸附面確有分形特性。

糯米粉吸附面的分形維數與Yano等[7]關于小麥粉、大米粉的分形維數測定結果（2.34、2.29）較為接近。需要說明的是，Yano等是根據顆粒表面積與粒徑之間的關系估算分形維數的。Aguerre等[39]采用與BET理論相同的假設前提，推導出一個分形表面的多層吸附等溫模型，據此估算分形維數，得出多數淀粉的分形維數接近3，小麥粉的為2.6；Resio等[40]根據該方法估算出莧菜籽淀粉在25、35、45 ℃時的分形維數分別為2.80、2.75、2.78，均高于本研究的結果。隨溫度升高糯米粉的分形維數略有降低，這一點與Aguerre等[39]關于淀粉的研究結果一致。

對于糯米粉，由于1/m＞1/3，說明糯米粉中水分的解吸過程主要是由毛細管凝結作用控制的；而馬鈴薯淀粉的水分吸附是由范德華力和毛細管凝結共同決定的[36]。

表面分形維數反映表面結構的不均勻性和自相似性，值越大說明其結構越接近三維多孔性結構。因此，分形維數的值越高，吸附表面越粗糙，吸附能力也就越強[41]，對比本研究關于不同溫度條件下的分形維數和平衡水分含量，也印證了這一點。

3 討 論

由式（17）、（19）計算的結果均低于采用BET方法計算的結果（圖6）；式（19）給出的結果更接近BET方法結算的結果。盡管Caurie等[35]關于馬鈴薯淀粉和玉米淀粉的結果也顯示出相同的趨勢，至于為何不同計算方法所得結果有很大差異，該作者并未關注。下面對此做深入討論。

3.1 BET法與Caurie1981法的比較

不同的單分子層水分含量：由前面的計算方程知，Caurie1981法與BET法都是根據單分子層水分含量計算單位吸附面積的，但BET方法中的單分子層水分是采用BET方程得到的，而Caurie1981法中是采用方程（16）得到的。Caurie1981法計算出的糯米粉的單分子層水分含量在10、20、30 ℃時分別為9.56、8.37、7.70 g/100 g，大于BET方程的計算結果（見2.4節）。僅就20 ℃條件下糯米粉的單分子層水分含量而言，則有M0=0.818MC0。

不同的單分子層密度：在由單分子層水分質量計算其對應體積的過程中，Caurie1981方法中采用的是被吸附單分子層的容積密度（式（16）中的CC）。針對本研究，該值在1.382～1.393 kg/m3之間，可取其平均值1.388 kg/m3。而BET方法在計算單位吸附面積時是以單個水分子為計算單位的，可以認為采用的是單個水分子的密度。關于被吸附的單個水分子的體積，可以參考與氨基酸結合的水分子的體積，25.8×10－30m3[42]。由水分子的質量，可計算出水分子的密度為1.159 kg/m3。則相同質量時，VBET=1.198VC1981。

不同的面積計算策略：BET法計算的實際上是單分子吸附層中所有水分子的截面積（實際上為水分子在吸附面上的正投影面積）之和，并不包括投影圓之間的間隙。Caurie1981法則是根據單分子層的總體積除以層高（認為等于水分子的直徑，關于這種假設的合理性，后面有詳細論述）。根據Caurie1981法，假設水分子球的半徑為rw，單位吸附面積為SC1981，則體積為V=2rwSC1981。由于兩種方法采用相同的分子直徑取值（BET法中面積取值為1.06×10-19m2，該面積的圓所對應的直徑為3.673×10-10m，正是Caurie1981法的分子直徑），則BET方法中體積為V的吸附層中，半徑為rw的球體的數量為單位吸附面積為，將上述V=2rwSC1981代入該式得SBET=1.5SC1981。

針對本研究中20 ℃時的情況，上述3 個因素疊加的結果則是：SBET=0.818×1.198×1.5SC1981=1.470SC1981。而根據BET法和Caurie1981法計算出的單位吸附面積分別為242.74、164.35 m2/g，兩者之比為1.477，十分接近1.470。這說明上述關于BET法和Caurie1981法產生差異原因的分析是合理的。

3.2 Cauire1981法與Caurie2012法的比較

由式（16）推導出式（17）的過程中假設水分子的直徑為3.673×10-10m，而推導出式（20）的過程中認為水分的直徑為2.75×10-10m。如果這兩個方法中采用相同的分子直徑取值，計算出的單位吸附面積幾乎相等（結果略）。

3.3 水分子直徑的取值

水分子直徑的取值為3.673×10-10m還是2.75×0-10m更合理？Caurie1981法中水分子的直徑3.673×10-10m是根據水分子的截面積10.6×10-20m2計算出來的，并以此作為單分子層的厚度[34]。該文獻并未給出該截面積數值的出處。而文獻[29]中給出的吸附在固體表面的水分子的截面積正是這個值。問題的關鍵是吸附在固體表面的水分子未必呈球形，僅根據其截面積推算其在另一維上的長度是不合適的，因而，Caurie1981法中采用3.673×10-10m作為單分子吸附層的厚度顯得不妥。Caurie等[35]的研究中取水分子的直徑為2.75×10-10m，這個值十分接近水分子的范德華半徑（2.82×10-10m[43]），但顯然不應當將此理解為水分子的直徑，因為如此的話，水分子的體積為10.88×10-30m3，比文獻[42]報道的值低得多。將2.75×10-10m理解為被吸附的單分子層的厚度更合理。如將水分子簡單看作截面為10.6×10-20m2、高為2.75×10-10m的柱體，其體積為29.15×10-30m3。該值比上述的與氨基酸結合的水分子的體積明顯要大，而僅略小于大量水體中水分子的體積29.70×10-30m3[42]。實際上，處于被吸附狀態的水分子的截面積會因吸附質的不同而有所變化。有報道顯示，吸附在纖維素上的水分子的投影面積為14.8×10-20m2[44]。

最后需要說明的是，盡管BET法與Caurie法之間存在前述的3 個差異，但兩者之間有著共同點：都認為吸附面積是被吸附的分子所占據的實際面積，并不包括每個分子所占面積之間的“空白”面積。當然，在Caurie法中，認為分子之間不存在“空白”面積。所以，單位吸附面積并不等同于吸附表面的實際面積。實際面積=吸附面積+“空白”面積。

4 結 論

糯米粉中水分解吸屬于Ⅱ型吸附。溫度越高，其平衡水分含量越低，Lewicki模型能更好低描述其解吸過程中平衡水分含量隨水分活度的變化，而GDW模型可以同時預測平衡水分含量隨溫度和水分活度的變化。

糯米粉中的孔屬于微孔和介孔，孔徑分布為單態分布，不同溫度下孔徑分布的不同主要體現在微孔部分。微孔體積、SBET分別在11.70～9.05 mL/100 g、223～269 m2/g之間，隨溫度升高而減小。吸附面具有分形特性，分形維數在2.05～2.19之間，隨溫度降低而增大。

導致BET法與Cauire法計算單位吸附面積結果顯著不同的根本原因在于：單分子層水分含量、單分子層密度、面積計算策略3 個方面的差異。

[1] CHIRIFE J, IGLESIAS H A. Equations for fitting water sorption isotherms of food: part I-a review[J]. Journal of Food Technology, 1978, 13(3): 159-174. DOI:10.1111/j.1365-2621.1978.tb00792.x.

[2] ROCKLAND L B, BEUCHAT L R. Water activity: theory and applications to food[M]//Practical Applications of Water Activity and Moisture Relations in Foods. New York: Marcel Dekker Inc., 1987: 369-395.

[3] LABUZA T P, HYMAN C R. Moisture migration and control in multi-domain foods[J]. Trend in Food Science and Technology, 1998, 9(2): 47-55. DOI:10.1016/S0924-2244(98)00005-3.

[4] BERISTAIN C I, AZUARA E, VERNON-CARTER E J. Effect of water activity on the stability to oxidation of spray-dried encapsulated orange peel oil using mesquite gum (Prosopis julif l ora) as wall material[J]. Journal of Food Science, 2002, 67(1): 206-211. DOI:10.1111/j.1365-2621.2002.tb11385.x.

[5] 邱光應, 彭桂蘭, 吳紹鋒, 等. 花椒吸附等溫線及熱力學性質[J]. 食品科學, 2015, 36(21): 1-5. DOI:10.7506/spkx1002-6630-201521001.

[6] LOWELL S, SHIELDS J E, THOMAS M A. Characterization of porous solids and powders: surface area, pore size and density[M]. New York: Springer Science+Business Media, 2004: 94-100. DOI:10.1007/978-1-4020-2303-3.

[7] YANO T, NAGAI T. Fractal surface of starchy materials transformed with hydrophilic alcohols[J]. Journal of Food Engineering, 1989, 10(2): 123-133. DOI:10.1016/0260-8774(89)90031-9.

[8] AZUARA E, BERISTAIN C I. Enthalpic and entropic mechanisms related to water sorption of yogurt[J]. Drying Technology, 2006, 24(11): 1-7. DOI:10.1080/07373930600961173.

[9] PASCUAL-PINEDA L A, FLORES-ANDRADE E, ALAMILLABELTRAN L, et al. Micropores and their relationship with carotenoids stability: a new tool to study preservation of solid foods[J]. Food Bioprocess Technology, 2014, 7(4): 1160-1170. DOI:10.1007/s11947-013-1162-0.

[10] JARZEBSKI A B, LORENC J, PAJAK L. Surface fractal characteristics of silica aerogels[J]. Langmuir, 1997, 13(5): 1280-1285. DOI:10.1021/la960011z.

[11] AOUAINI F, KNANI S, YAHIA M B, et al. Statistical physics studies of multilayer adsorption isotherm in food materials and pore size distribution[J]. Physica A: Statistical Mechanics and Its Applications, 2015, 432: 373-390. DOI:10.1016/j.physa.2015.03.052.

[12] CAURIE M. A method to correct the wide discrepancy between Brunauer, Emmett and Teller water and N2surface areas of adsorbents[J]. International Journal of Food Science and Technology, 2012, 47(11): 2366-2371. DOI:10.1111/j.1365-2621.2012.03111.x.

[13] LI J R, HSIEH P Y H. Traditional Chinese food technology and cuisine[J]. Asia Pacific Journal of Clinical Nutrition, 2004, 13(2): 147-155.

[14] DURAKOVA A G, MENKOV N D. Moisture sorption characteristics of rice flour[J]. Die Nahrung, 2004, 48(2): 137-140. DOI:10.1002/ food.200300387.

[15] BRETT B, FFGUEROA M, SANDOVAL A J, et al. Moisture sorption characteristics of starchy products: oat flour and rice flour[J]. Food Biophysics, 2009, 4(3): 151-157. DOI:10.1007/s11483-009-9112-0.

[16] 劉成梅, 周國輝, 萬婕, 等. 大米淀粉解吸等溫線與吸附等溫線的擬合模型研究[J]. 食品工業科技, 2014, 35(10): 198-201. DOI:10.13386/j.issn1002-0306.2014.10.036.

[17] GREENSPAN L. Humidity fi xed points of binary saturated aqueous solutions[J]. Journal of Research of the National Bureau of Standards, Section A: Physics and Chemistry, 1977, 81(1): 89-96.

[18] BRUNAURE S, DEMING L S, TELLER E. On a theory of van der Waals adsorption of gases[J]. Journal of the American Chemical Society, 1940, 62(7): 1723-1732. DOI:10.1021/ja01864a025.

[19] DEVAHASTIN S. Mujumdar’s practical guide to industrial drying[M]// Fundamental Principles of Drying. Montreal: Exergex Corp, 2000: 1-22.

[20] RAO M A, RIZVI S S H. Engineering properties of foods[M]// Thermodynamics of Foods in Dehydration. New York: Marcel Dekker, 1995: 133-214.

[21] American Society of Agricultural and Biological Engineers. ASAE D245.6 Moisture relationship of plant-based agricultural products[S]. St. Joseph: American Society of Agricultural and Biological Engineers, 2007: 496-512.

[22] LEWICKI P P. A three parameter equation for food moisture sorption isotherms[J]. Journal of Food Process Engineering, 1998, 21(2): 127-144. DOI:10.1111/j.1745-4530.1998.tb00444.x.

[23] FURMANIAK S, TERZYK A P, GOLEMBIEWSKI R, et al. Searching the most optimal model of water sorption on foodstuffs in the whole range of relative humidity[J]. Food Research International, 2009, 42(8): 1203-1214. DOI:10.1016/j.foodres.2009.06.004.

[24] GARCIA-PEREZ J V, CARCEL J A, CLEMENTE G, et al. Water sorption isotherms for lemon peel at different temperatures and isosteric heats[J]. LWT-Food Science and Technology, 2008, 41(1): 18-25. DOI:10.1016/j.lwt.2007.02.010.

[25] SINGH R R B, RAO K H, ANJANEYULU A S R, et al. Moisture sorption properties of smoked chicken sausages from spent hen meat[J]. Food Research International, 2001, 34(2/3): 143-148. DOI:10.1016/S0963-9969(00)00145-9.

[26] SINGH R R B, RAO K H, ANJANEYULU A S R, et al. Water desorption characteristics of raw goat meat: effect of temperature[J]. Journal of Food Engineering, 2006, 75(2): 228-236. DOI:10.1016/ j.jfoodeng.2005.04.013.

[27] ROSA G S, MORAES M A, PINTO L A A. Moisture sorption properties of chitosan[J]. LWT-Food Science and Technology, 2010, 43(3): 415-420. DOI:10.1016/j.lwt.2009.09.003.

[28] de ALBANI M I D, ARCIPRETE C P. A study of pore size distribution and mean pore size on unsupported gamma-alumina membranes prepared by modifications introduced in the alkoxide hydrolysis step[J]. Journal of Membrane Science, 1992, 69(1/2): 21-28. DOI:10.1016/0376-7388(92)80164-F.

[29] LIVINGSTON H K. Cross-sectional areas of molecules adsorbed on solid surfaces[J]. Journal of the American Chemical Society, 1944, 66(4): 569-573. DOI:10.1021/ja01232a021.

[30] KARATHANOS V T, SARAVACOS G D. Porosity and pore size distribution of starch materials[J]. Journal of Food Engineering, 1993, 18(3): 259-280. DOI:10.1016/0260-8774(93)90090-7.

[31] FLETCHER A J, THOMAS K M. Compensation effect for the kinetics of adsorption/desorption of gases/vapors on microporous carbon materials[J]. Langmuir, 2000, 16(15): 6253-6266. DOI:10.1021/ la9916528.

[32] SPIEWECK F, BETTIN H. Review: solid and liquid density determination[J]. Technisches Messen, 1992, 59(7/8): 285-292. DOI:10.1524/teme.1992.59.78.285.

[33] VIGANó J, AZUARA E, TELIS V R N, et al. Role of enthalpy and entropy in moisture sorption behavior of pineapple pulp powder produced by different drying methods[J]. Thermochimica Acta, 2012, 528: 63-71. DOI:10.1016/j.tca.2011.11.011.

[34] CAURIE M. Derivation of full range moisture isotherms[M]// ROCKLAND L B, STEWART G F. Water Activity: Influences on Food Quality, New York: Academic Press, 1981: 63-87.

[35] CAURIE M. The unimolecular character of the classical Brunauer, Emmett and Teller adsorption equation and moisture adsorption[J]. International Journal of Food Science and Technology, 2005, 40(3): 283-293. DOI:10.1111/j.1365-2621.2004.00892.x.

[36] CZEPIRSKI L, KOMOROWSKA-CZEPIRSKA E, SZYMONSKA J. Fitting of different models for water vapor sorption on potato starch granules[J]. Applied Surface Science, 2002, 196(1/2/3/4): 150-153. DOI:10.1016/S0169-4332(02)00050-8.

[37] NEIMARK A V. Determination of the surface fractal dimensionality from the results of an adsorption experiment[J]. Russian Journal of Physical Chemistry, 1990, 64(10): 1397-1403.

[38] SOKOLOWSKA Z, BOROWKO M, RESZKO-ZYGMUNT J, et al. Adsorption of nitrogen and water vapor by alluvial soils[J]. Geoderma, 2002, 107(1/2): 33-54. DOI:10.1016/S0016-7061(01)00137-9.

[39] AGUERRE R J, VIOLLAZ P E, SUAREZ C. A fractal isotherm for multilayer adsorption in foods[J]. Journal of Food Engineering, 1996, 30(1/2): 227-238. DOI:10.1016/S0260-8774(97)89377-6.

[40] RESIO A C, AGUERRE R J, SUAREZ C. Analysis of the sorptional characteristics of amaranth starch[J]. Journal of Food Engineering, 1999, 42(1): 51-57. DOI:10.1016/S0260-8774(99)00103-X.

[41] XIONG J, LIU X, LIANG L X. An investigation of fractal characteristics of marine shales in the southern china from nitrogen adsorption data[J]. Journal of Chemistry, 2015: 1-12.

[42] PERKINS S J. Protein volumes and hydration effects. The calculations of partial specif i c volumes, neutron scattering matchpoints and 280-nm absorption coefficients for proteins and glycoproteins from amino acid sequences[J]. European Journal of Biochemistry, 1986, 157(1): 169-180. DOI:10.1111/j.1432-1033.1986.tb09653.x.

[43] FRANKS F. Water: a matrix of life[M]. 2nd ed. Cambridge: Royal Society of Chemistry, 2000: 9-14. DOI:10.1021/ed078p593.1.

[44] TOPGAARD D, SODEMAN O. Diffusion of water absorbed in cellulose fibers studied with1H-NMR[J]. Langmuir, 2001, 17(9): 2694-2702. DOI:10.1021/la000982l.

Water Desorption Isotherm and Pore Characteristics of Glutinous Rice Flour

ZHANG Hua1,2, ZHAO Xuewei1,2, BAI Yanhong1,2, DUAN Ruiqian1
(1. School of Food and Bioengineering, Zhengzhou University of Light Industry, Zhengzhou 450002, China；2. Henan Province Collaborative Innovation Center for Food Production and Safety, Zhengzhou 450002, China)

Water desorption isotherm of glutinous rice flour is necessary for understanding water transfer from it to environment. On the other hand, the isotherm data can be used to explore pore characteristics a solid matrix, which will allow a deeper understanding of water adsorption mechanism. Equilibrium water contents of glutinous rice flour were determined after conditioning to 10 different water activity levels at 10, 20 and 30 ℃ using a static gravimetric method. The obtained desorption data were fi tted with four selected isotherm models, and were further mined to estimate pore properties of glutinous rice fl our. Glutinous rice fl our exhibited a typical type-Ⅱ adsorption isotherm. The Lewicki model performed best in predicting water desorption of glutinous rice fl our at a constant temperature, while the GDW model could predict the moisture content as a function of water activity and temperature. The pores in glutinous rice fl our included micropore and mesopore. Pore size distribution posed a single peak located in micropore zone whose value increased with the decrease in temperature. Consequently, micropore volume also increased as temperature decreased. The pore surface in glutinous rice fl our was fractal in nature. The more porous surface structure and the larger specif i c surface area at lower temperature contributed to the increased adsorption capacity.

glutinous rice flour; desorption; pore size; pore size distribution; specific surface area; micropore volume; fractal dimension

10.7506/spkx1002-6630-201707015

TS213.2

A

1002-6630（2017）07-0088-08

張華, 趙學偉, 白艷紅, 等. 糯米粉的水分等溫解吸及孔特性[J]. 食品科學, 2017, 38(7): 88-95. DOI:10.7506/spkx1002-6630-201707015. http://www.spkx.net.cn

ZHANG Hua, ZHAO Xuewei, BAI Yanhong, et al. Water desorption isotherm and pore characteristics of glutinous rice flour[J]. Food Science, 2017, 38(7): 88-95. (in Chinese with English abstract)

10.7506/spkx1002-6630-201707015. http://www.spkx.net.cn

2016-03-12

“十二五”農村領域國家科技計劃項目（2012BAD37B06-05）；2014年河南省重大科技專項（141100110400）；鄭州市重大科技專項（152PZDZX026）

張華（1975—），男，副教授，博士，研究方向為米面食品加工。E-mail：zhh7510@126.com