Rudin改良法測定啤酒泡持性條件的優(yōu)化

肖 俠，李 紅，張 珍，*，邵威平，陳光斌，康開萍

(1.甘肅農(nóng)業(yè)大學(xué)食品科學(xué)與工程學(xué)院，甘肅蘭州 730070;2.中國食品發(fā)酵工業(yè)研究院，北京 100027)

啤酒泡沫穩(wěn)定性的測定方法有很多種，如秒表法、Sigma法、Rudin 法、NIBEM 法、瞬時發(fā)泡法等［1］，每種方法各有其利弊。隨著對啤酒泡沫研究的深入，研究啤酒過程產(chǎn)品中的起泡因子也逐漸增多，就會遇到使溶液產(chǎn)生泡沫的問題，只有外源氣體的充入才能解決不含 CO2溶液的起泡問題［2－3］。Rudin法可以排除氣體對酒體起泡物質(zhì)的干擾，其原理:在恒定溫度下，以恒定氣流通入待測樣產(chǎn)生泡沫，由于殘存泡沫相當(dāng)?shù)囊后w體積的對數(shù)與時間呈規(guī)則的線性關(guān)系，所以測定兩特定標(biāo)記刻度間液體形成的時間，便可得知泡沫衰減周期［3］。本實(shí)驗(yàn)應(yīng)用Rudin改良法測定啤酒泡沫的泡持性，以泡持值作為評價指標(biāo)，采用中心旋轉(zhuǎn)組合實(shí)驗(yàn)設(shè)計(jì)分析了CO2流量、氣體總閥出口壓力以及待測試樣的預(yù)處理溫度對測定泡持性的影響，通過響應(yīng)面法對Rudin改良法測定啤酒泡沫穩(wěn)定性的條件進(jìn)行優(yōu)化，以期為實(shí)驗(yàn)室測定啤酒泡沫泡持性及測定不含CO2溶液的起泡問題提供理論依據(jù)。

1 材料與方法

1.1 材料與儀器

原料 某酒廠生產(chǎn)原麥汁濃度8°P的同批次成品酒。

HX－105恒溫水浴箱 上?？戮S儀器有限公司。

1.2 實(shí)驗(yàn)室簡易Rudin改良裝置

如圖1所示。

1.3 實(shí)驗(yàn)方法

1.3.1 Rudin改良法測定的原理 該方法利用一根細(xì)長玻璃管作為發(fā)泡管，底部引入定量啤酒，并恒速緩慢通入CO2盡可能將啤酒全部轉(zhuǎn)變成泡沫，停止充氣，于是泡沫開始塌陷，酒液在細(xì)管中上升。泡沫的塌陷存在這樣一個規(guī)律，即泡沫的體積(可以轉(zhuǎn)化成酒液的體積)的對數(shù)與時間成線性關(guān)系［3－4］，為了得到較為精確的結(jié)果，取泡沫塌陷過程中的數(shù)組時間和體積值，擬合一次線性方程，求出k值。以酒液通過細(xì)管上兩個標(biāo)記之間的時間間隔，或其半衰期t1/2作為泡持值［5］。泡持值的計(jì)算詳見1.3.2。

圖1 Rudin改良法測定啤酒泡沫穩(wěn)定性的簡易裝置Fig.1 The simple apparatus of modified Rudin method to measure the foam stability

1.3.2 泡持值的計(jì)算 以時間為橫坐標(biāo)，泡沫的液體相當(dāng)體積 V(foam)，即 50mL－V(液)的對數(shù)為縱坐標(biāo)，取V(foam)為 10、15、20、25、30、35、40、45mL 八個點(diǎn)作圖，并求出趨勢方程、R2值。由趨勢方程可知直線的斜率，以酒液通過發(fā)泡管上兩個標(biāo)記之間的時間間隔(式1)，或其半衰期t1/2作為泡持值(式2)。

式中:V表示泡沫相當(dāng)于液體的總體積，即通氣前液體的總體積［1］。

1.3.3 單因素實(shí)驗(yàn) 在預(yù)實(shí)驗(yàn)的基礎(chǔ)上，選擇CO2流量、待測試樣的預(yù)處理溫度以及氣體總閥出口壓力三個因素進(jìn)行單因素實(shí)驗(yàn)，分別考察其對測定泡持性的影響，確定各因素的優(yōu)化區(qū)間。量取待測試樣50mL，選擇 CO2流量梯度分別為 220、260、300、340、380、400、500mL/min;溫度梯度為 10、15、20、25、30℃;氣體總閥出口壓力梯度分別為0.1、0.2、0.3、0.4、0.5MPa。每次試樣做3個平行。

1.3.4 響應(yīng)面法優(yōu)化測定條件 根據(jù)單因素實(shí)驗(yàn)，選取CO2流量(mL/min)、待測樣品預(yù)處理溫度(℃)、氣體總閥出口壓力(MPa)三個因素為自變量，以泡持值(s)為響應(yīng)值(Y)，進(jìn)行中心旋轉(zhuǎn)組合實(shí)驗(yàn)設(shè)計(jì)，共17個實(shí)驗(yàn)點(diǎn)，12個析因點(diǎn)和5個零點(diǎn)，由單因素實(shí)驗(yàn)確定各因素變化區(qū)間。每個實(shí)驗(yàn)點(diǎn)均做三個平行，應(yīng)用Design－Expert.8.05b軟件建立方差分析模型，選擇p＜0.05的因素作為主效應(yīng)因素，因素水平編碼見表1。

1.3.5 模型的驗(yàn)證 通過響應(yīng)面分析法優(yōu)化測定的條件，并在優(yōu)化條件下測定試樣的泡持性，通過比較預(yù)測值和實(shí)驗(yàn)值來驗(yàn)證模型的有效性。

式中:V1和V2分別表示兩標(biāo)記處的泡沫的液體相當(dāng)體積。

表1 中心旋轉(zhuǎn)組合設(shè)計(jì)因素水平編碼表Table1 Codes and levels of factors for central composite design

2 結(jié)果與分析

2.1 CO2流量對Rudin改良法測泡沫穩(wěn)定性的影響

將酒樣置于(20±0.5)℃恒溫水浴中0.5h，待測。調(diào)試Rudin改良裝置，確保氣體總閥出口壓力為0.2MPa的情況下，調(diào)節(jié)氣體流量，觀察不同氣體流量對泡持測定的影響，結(jié)果如圖2所示。

圖2 CO2流量對泡持性的影響Fig.2 Effect of CO2flow on foam stability

由圖2所示，CO2流量為340mL/min時，泡持值最大，當(dāng)流量大于或小于340mL/min時，泡持值都會降低，當(dāng)流量太小時，起泡較慢，時間長，泡持值小，當(dāng)流量太大時，雖然耗費(fèi)時間短，但起始泡沫比較粗大，隨時間變化泡沫越來越粗糙，且容易造成泡沫和酒液分界面不清晰，誤差較大［6］。故選擇CO2流量為340mL/min為零點(diǎn)。

2.2 待測試樣預(yù)處理溫度對泡持的影響

CO2流量為340mL/min，氣體總閥出口壓力為0.2MPa 的情況下，將酒樣分別在 10、15、20、25、30℃條件下恒溫0.5h，測其泡持性，觀察不同溫度對泡持測定的影響，結(jié)果如圖3所示。

圖3 待測樣品預(yù)處理溫度對泡持性的影響Fig.3 Effect of temperature of sample pretreatment on foam stability

由圖3可知，泡沫的穩(wěn)定性隨溫度的升高而降低，溫度越高泡沫消失得越快，溫度較低，泡沫消失得越慢，這是因?yàn)闇囟仍礁?，啤酒泡沫中各單個氣泡體積增大，使得表面粘度降低，導(dǎo)致泡沫穩(wěn)定性下降［7］。由于實(shí)際操作中10℃的溫度較難維持，并且國標(biāo)中秒表法測量啤酒泡沫穩(wěn)定性是將酒樣恒溫于20℃水浴中，所以進(jìn)行中心旋轉(zhuǎn)組合設(shè)計(jì)時選擇20℃為零點(diǎn)。

2.3 氣體出口閥壓力對泡持的影響

待測酒樣置于(20±0.5)℃水浴中恒溫0.5h，并將CO2流量調(diào)節(jié)為340mL/min，通過改變氣體出口閥壓力的大小，壓力梯度分別為 0.1、0.2、0.3、0.4、0.5MPa，觀察不同壓力對泡持測定的影響，結(jié)果如圖4所示。

圖4 氣體總閥出口壓力對泡持性的影響Fig.4 Effect of gas outlet pressure of total valve on foam stability

由圖4可知，測定的泡持值呈現(xiàn)先升后降的趨勢，當(dāng)壓力為0.2MPa時，測量的泡持值最大，壓力為0.3MPa時，泡持值下降不明顯，隨后再增加壓力時，泡持值逐漸下降，當(dāng)壓力增大到0.5MPa時，泡持值又呈現(xiàn)小幅度的上升，但壓力過大時，起泡過快，氣泡不均衡，容易破裂［8］。故選擇壓力0.2MPa。

2.4 影響因素的響應(yīng)面分析

以CO2流量、待測樣品預(yù)處理溫度、氣體總閥出口壓力為自變量進(jìn)行響應(yīng)面實(shí)驗(yàn)，實(shí)驗(yàn)設(shè)計(jì)和響應(yīng)值見表2。

表2 中心旋轉(zhuǎn)組合設(shè)計(jì)和響應(yīng)值Table2 Results of response surface design

用Design－Expert.8.05b軟件進(jìn)行二次多項(xiàng)回歸擬合，所得多元二次回歸方程如下:

泡持值的預(yù)測值與實(shí)際實(shí)驗(yàn)值擬合情況見圖5，顯示預(yù)測值和實(shí)驗(yàn)值擬合良好。模型p＜0.0001，決定系數(shù)(R－Squared)為 0.9978，校正系數(shù)(Adj R－Squared)為0.9950，響應(yīng)變量 R2高于0.8000，證明此模型顯著，可充分地反映各變量之間的關(guān)系。

圖5 泡持值的實(shí)際值與預(yù)測值的對應(yīng)關(guān)系Fig.5 The correspondence between the predicted value and the actual value

對模型中的回歸系數(shù)進(jìn)行顯著性檢驗(yàn)，結(jié)果見表3。

表3 中心旋轉(zhuǎn)組合設(shè)計(jì)方差分析表Table3 Variance analysis of regression equation

由表3可以看出，一次項(xiàng)均為極顯著(p＜0.01)即CO2流量、待測樣品預(yù)處理溫度、氣體總閥出口壓力均對泡持值的測定有極顯著的影響;交互項(xiàng)中CO2流量和待測樣品預(yù)處理溫度的交互作用對泡持值的測定有顯著影響(p＜0.05)，CO2流量和氣體總閥出口壓力的交互作用對泡持值的測定有極顯著的影響(p＜0.01);二次項(xiàng)均為極顯著(p＜0.01)。根據(jù)多元二次回歸方程的各個系數(shù)絕對值的大小可判斷三個因素對泡持值測定的影響順序?yàn)闅怏w總閥出口壓力＞CO2流量＞待測樣品預(yù)處理溫度。剔除不顯著項(xiàng)(p＞0.05)，得優(yōu)化回歸模型為:

分別對方程各自變量求極值，得到CO2流量為338.07mL/min，待測樣品預(yù)處理溫度為19.62℃，氣體總閥出口壓力為0.24MPa，在此最優(yōu)條件下，預(yù)測所得的泡持值為111.363s。由于在實(shí)際操作中此三個變量均不能取到小數(shù)點(diǎn)后面的位數(shù)，所以取整進(jìn)行驗(yàn)證性實(shí)驗(yàn)。

圖6～圖8分別給出了三個因素之間兩兩交互作用對泡持值測定的響應(yīng)面曲面圖。

由圖6所示，CO2流量和待測樣品預(yù)處理溫度之間的交互作用顯著，CO2流量控制在328～342mL/min，待測樣品預(yù)處理溫度控制在17～20℃之間，泡持值可達(dá)到最大值。

圖6 CO2流量和待測樣品預(yù)處理溫度對泡持值的影響Fig.6 Response surface diagram of CO2flow and temperature of sample pretreatment

由圖7所示，CO2流量和氣體總閥出口壓力的交互作用極顯著，CO2流量控制在320～352mL/min，氣體總閥出口壓力控制在0.2～0.27MPa內(nèi)，泡持值可達(dá)到最大值。

圖7 CO2流量和氣體總閥出口壓力對泡持值的影響Fig.7 Response surface diagram of CO2flow and gas outlet pressure of total valve

由圖8所示，當(dāng)氣體總閥出口壓力大于0.15MPa，待測樣品預(yù)處理溫度在15～20℃之間時，泡持值達(dá)到最大值，當(dāng)壓力小于0.15MPa，溫度大于20℃時，泡持值下降。

2.5 驗(yàn)證實(shí)驗(yàn)

預(yù)測CO2流量為338.07mL/min，待測樣品預(yù)處理溫度為19.62℃，氣體總閥出口壓力為0.24MPa，實(shí)際操作中此三個變量均不能取到小數(shù)點(diǎn)后面的位數(shù)，所以取整進(jìn)行驗(yàn)證性實(shí)驗(yàn)，即流量取340mL/min，溫度取20℃，壓力取0.3MPa共進(jìn)行5次平行實(shí)驗(yàn)，泡持值的平均值為110.813s，預(yù)測所得的泡持值為111.363s，二者的相對誤差為0.49%。

圖8 待測樣品預(yù)處理溫度和氣體總閥出口壓力對泡持值的影響Fig.8 Response surface diagram of temperature of sample pretreatment and gas outlet pressure of total valve

3 結(jié)論

3.1 三個測定條件對Rudin改良法測定啤酒泡持性影響的排序?yàn)?氣體總閥出口壓力＞CO2流量＞待測樣品預(yù)處理溫度。

3.2 優(yōu)化后的Rudin改良法測定啤酒泡持性條件為:CO2流量為340mL/min，待測樣品預(yù)處理溫度為20℃，氣體總閥出口壓力為0.3MPa。實(shí)際值與預(yù)測值的相對誤差為0.49%。

3.3 驗(yàn)證實(shí)驗(yàn)結(jié)果表明，泡持值的平均值為110.813s，預(yù)測所得的泡持值為111.363s，二者的相對誤差為0.49%。

［1］葉俊華，張峻炎，陸健.啤酒泡沫泡持性測定的研究進(jìn)展［J］.釀酒，2003，30(1):38－39.

［2］A D Rudin.Measurement of the foam stability of beers［J］.Inst Brew，1957，63:506－509.

［3］葉俊華，付兆輝，陸健.Rudin改良法在啤酒泡沫穩(wěn)定性測定中的應(yīng)用［J］.食品工業(yè)科技，2003，24(11):30－32.

［4］Derdelinckx G，Kopmanis J，Haboucha J，et al.Measurement of foaming power by the modified Rudin method［J］.Cerevisia，1984(4):211－214.

［5］管敦儀.啤酒工業(yè)手冊［M］.北京:中國輕工業(yè)出版社，1999:579－582.

［6］Lynch D M，Bamforth C W.Measurement and characterization of bubble nucleation in beer［J］.Food Sci，2002，67:2696－2701.

［7］Bamforth C W.Perc

eptions of beer foam［J］.Inst Brew，2000，106(4):229－238.

［8］董建軍，黃淑霞，王成紅，等.Rudin法測定啤酒泡沫穩(wěn)定性的研究［J］.啤酒科技，2006，12:28－34.

［9］Sharpe F R.Institute ofBrewing AnalysisCommittee Collaborative determination of beer foam stability by Rudin and Nibem［J］.Inst Brew，1997，103(5):277－278.

［10］Charles Bamforth.Beer:Tapping into the Art and Science of Brewing［M］.New York:Plenum Press，1999:179.

［11］Evan Evans D，Surrel Anne.Comparison of foam quality and the influence of hop α－acids and proteins using five foam analysismethods［J］.ASBC，2008，66(1):1－10.

［12］顧國賢.釀造酒工藝學(xué)［M］.北京:中國輕工業(yè)出版社，1996:336－337.

［13］張秀英.啤酒泡沫的測定［J］.釀酒，2001，28(4):21－23.