微真空條件下基于色澤參數(shù)變化的西蘭花葉綠素質(zhì)量分?jǐn)?shù)預(yù)測(cè)模型

王麗嬌， 張圣杰， 李文香， 王士奎， 孫樹(shù)杰

（青島農(nóng)業(yè)大學(xué) 食品科學(xué)與工程學(xué)院，山東 青島266109）

西蘭花，俗稱青花菜、莖椰菜等，屬十字花科蕓薹屬甘藍(lán)變種。其營(yíng)養(yǎng)物質(zhì)豐富，同時(shí)，西蘭花富含VC，可延緩衰老、提高免疫力、增強(qiáng)肝臟的解毒功能；食用富含高纖維的西蘭花還可有效降低腸胃對(duì)葡萄糖的吸收，是糖尿病患者的天然藥膳[1]，而且還含有抗癌作用的活性成分，如硫代葡萄糖苷[2]。西蘭花價(jià)格隨季節(jié)變化幅度大，出口額大，是速凍蔬菜的主要代表[3]，但其呼吸旺盛，常溫下葉綠素迅速降解，花球極易黃化變質(zhì)。評(píng)價(jià)西蘭花品質(zhì)的指標(biāo)主要包括新鮮度、色澤及營(yíng)養(yǎng)成分等，其衰老的原發(fā)過(guò)程和真正標(biāo)志就是葉綠素的降解。而葉綠素是葉綠酸的酯，在采后貯藏的過(guò)程中，葉綠素在葉綠素酶和脫鎂葉綠素酶的作用下[4]，極易脫去分子中的植醇和卟啉環(huán)中的鎂離子，生成黃褐色的脫鎂脫植醇葉綠素和脫鎂葉綠素，導(dǎo)致西蘭花褪綠黃化。對(duì)于葉綠素質(zhì)量分?jǐn)?shù)的測(cè)定，現(xiàn)階段多采用分光光度計(jì)法，除用于葉綠素提取的丙酮試劑有一定毒性外，操作也比較麻煩，且測(cè)量過(guò)程樣品的采集會(huì)對(duì)西蘭花造成一定的損失，為此有學(xué)者通過(guò)探討光譜特性來(lái)研究果蔬表面色澤的變化[5]。

在果蔬采后貯藏加工過(guò)程中，CIE-Lab系統(tǒng)是一種有效可信賴的評(píng)價(jià)色澤的方法[6-7]。它利用色差計(jì)測(cè)定色澤參數(shù)來(lái)研究果蔬表面的色澤變化，不僅操作簡(jiǎn)便沒(méi)有破壞性，而且可利用動(dòng)力學(xué)進(jìn)行研究和分析[8-10]。由于貯藏期間西蘭花葉綠素不斷降解的過(guò)程中a值（紅度，正值代表紅色，負(fù)值代表綠色）和b值（黃度，正值代表黃色，負(fù)值代表藍(lán)色）[11]均有明顯的變化，所以可通過(guò)a值和b值的變化了解采后西蘭花褪綠黃化的情況。

果蔬采后貯藏加工過(guò)程中動(dòng)力學(xué)的研究大多是從動(dòng)力學(xué)變化的角度探討果蔬品質(zhì)的損失。其品質(zhì)因子Q（如色澤、硬度及固形物等）的損失率可以表示為：dQ/dt=-k（Q）n，式中 n 為反應(yīng)級(jí)數(shù)；k 為反應(yīng)速率常數(shù)；t為貯藏時(shí)間。刁恩杰等（2010）[12]建立了脫水菠菜在貯藏過(guò)程中顏色變化的動(dòng)力學(xué)模型；劉春泉等（2011）[13]建立了毛豆仁低溫凍藏期間葉綠素降解和色澤變化的動(dòng)力學(xué)模型；任珂等 （2005）[14]研究發(fā)現(xiàn)表面色澤參數(shù)a、b值的變化可衡量包裝及不同溫度對(duì)青花菜采后色澤變化的影響。大多數(shù)都集中在不同采后處理對(duì)葉綠素降解動(dòng)力學(xué)模型的研究，而微真空貯藏條件下西蘭花采后色澤變化與葉綠素含量的相關(guān)性及其預(yù)測(cè)模型的研究尚屬空白?！拔⒄婵召A藏”是在專利（CN 1530290A）[15]技術(shù)基礎(chǔ)上開(kāi)發(fā)的一種輕度減壓貯藏技術(shù)，周莎莎等（2011）[16]采用微真空設(shè)施對(duì)萊陽(yáng)梨進(jìn)行了貯藏試驗(yàn)，并取得了良好的貯藏效果。作者在此基礎(chǔ)上，通過(guò)對(duì)微真空條件下西蘭花葉綠素質(zhì)量分?jǐn)?shù)與色澤變化的相關(guān)性研究，探討葉綠素的降解規(guī)律，建立葉綠素降解與色澤變化的預(yù)測(cè)模型，以期為采后西蘭花貯藏過(guò)程中的無(wú)損檢測(cè)提供理論支持。

1 材料與方法

1.1 材料

西蘭花：于2012年5月10日訂購(gòu)于青島市城陽(yáng)區(qū)蔬菜批發(fā)市場(chǎng)的經(jīng)銷商，選擇新鮮潔凈，花球緊密，花蕾未開(kāi)，成熟度一致，大小適中，無(wú)病蟲(chóng)害及其它傷害的個(gè)體。

1.2 儀器

微真空貯藏設(shè)備：青島農(nóng)業(yè)大學(xué)食品科學(xué)與工程學(xué)院研制；電子分析天平：上海奧豪斯國(guó)際貿(mào)易有限公司；754型紫外可見(jiàn)分光光度計(jì)：上海光譜儀器有限公司；CR-400色彩色差計(jì)：日本柯尼卡美能達(dá)有限公司。

1.3 材料的處理

1.3.1 試驗(yàn)組的處理 將西蘭花貯于壓力分別為53.33～66.66、66.66～79.99、79.99～93.32 kPa 及常壓，溫度為 （3±0.5）℃的微真空設(shè)備中，分別在第0，7，14，21，28，35，42，49 天進(jìn)行相關(guān)指標(biāo)的測(cè)定。

1.3.2 驗(yàn)證組的處理 將西蘭花貯于壓力分別為53.33～66.66、66.66～79.99、79.99～93.32 kPa 及常壓，溫度為 （3±0.5）℃的微真空設(shè)備中，分別在第1，8，15，22，29，36，43，50 天進(jìn)行相關(guān)指標(biāo)的測(cè)定。

1.4 葉綠素質(zhì)量分?jǐn)?shù)的測(cè)定

采用分光光度計(jì)法，參照曹建康等（2007）[17]的方法。

1.5 色澤的測(cè)定

西蘭花表面色澤采用色彩色差計(jì)測(cè)定，用標(biāo)準(zhǔn)陶瓷板標(biāo)定X=76.74，Y=81.07，Z=80.89作為工作標(biāo)準(zhǔn)[18]。將保鮮膜覆于西蘭花上，選擇其表面較為平整的點(diǎn)作為測(cè)量點(diǎn)，分別測(cè)a值和b值，計(jì)算色澤參數(shù)-a/b值，每個(gè)樣品取5個(gè)不同的點(diǎn)，取平均值。

1.6 色澤動(dòng)力學(xué)模型建立的方法

采用EXCEL軟件，利用回歸方程建立色澤動(dòng)力學(xué)模型。

2 結(jié)果與分析

2.1 葉綠素質(zhì)量分?jǐn)?shù)與色澤參數(shù)-a/b值的變化及其動(dòng)力學(xué)模型的建立

葉綠素是植物進(jìn)行光合作用的重要色素物質(zhì)。新鮮的西蘭花中含有大量葉綠素，但隨著貯藏時(shí)間的延長(zhǎng)，葉綠素質(zhì)量分?jǐn)?shù)不斷下降，導(dǎo)致西蘭花褪綠黃化，是影響西蘭花外觀品質(zhì)的關(guān)鍵[19]。由圖1可以看出，不同壓力條件下西蘭花貯藏期間葉綠素質(zhì)量分?jǐn)?shù)均呈不斷下降的趨勢(shì)，貯藏前期下降速度較慢 ， 貯 藏 至 28 d 時(shí) ，53.33～66.66、66.66～79.99、79.99～93.32 kPa及常壓下葉綠素的降解率分別為22.39%、32.81%、38.71%、40.30%，后期下降速度變快，至貯藏結(jié)束時(shí)，葉綠素的降解率分別為64.18%、70.31%、80.65%、85.07%，說(shuō)明微真空條件下采后西蘭花葉綠素的降解呈先慢后快的變化趨勢(shì)，且在53.33～66.66 kPa條件下的葉綠素降解速率最小。

圖1 不同壓力條件下葉綠素質(zhì)量分?jǐn)?shù)隨貯藏時(shí)間的變化Fig.1 Changes of chlorophyll content with storage time under different pressure conditions

經(jīng)回歸分析得，不同壓力條件下葉綠素質(zhì)量分?jǐn)?shù)隨時(shí)間的變化均符合零級(jí)動(dòng)力學(xué)反應(yīng)，53.33～66.66、66.66～79.99、79.99～93.32 kPa 及常壓下葉綠素質(zhì)量分?jǐn)?shù)的變化速率分別為：dQ1/dt=-0.008 5，dQ2/dt=-0.070 6，dQ3/dt=-0.071 4，dQ4/dt=-0.081 2，葉綠素降解的動(dòng)力學(xué)模型分別為=-0.008 5t+0.945 4，=-0.070 6t+0.756 4，=-0.071 4t+0.726 4，=-0.081 2t+0.782 9。

綠色色澤是西蘭花最為重要的外觀品質(zhì)之一，在西蘭花的貯藏過(guò)程中，色澤發(fā)生變化時(shí)，a值和b值也會(huì)隨之發(fā)生改變，其中，-a值主要是葉綠素呈現(xiàn)的顏色，葉綠素的大量損失造成-a值的下降[20]，同時(shí)b值也不斷改變，因此，可利用-a/b值的變化表示其色澤的變化。從圖2可以看出，貯藏前期在不同壓力條件下，-a/b值下降緩慢，貯藏至28 d時(shí)，53.33～66.66、66.66～79.99、79.99～93.32 kPa 及常壓下-a/b值分別下降了 21.37%、24.03%、35.71%、41.98%，說(shuō)明貯藏前期西蘭花的黃化速度較慢；后期下降較快，至貯藏結(jié)束時(shí)，-a/b值分別下降了54.2%、61.24%、69.05%、75.57%。表明微真空條件下，采后西蘭花在貯藏過(guò)程中綠色色澤不斷消退而逐漸黃化。

圖2 不同壓力條件下顏色參數(shù)-a/b值隨貯藏時(shí)間的變化Fig.2 Changes of color parameter-a/b value with storage time under different pressure conditions

經(jīng)回歸分析得，不同壓力條件下-a/b值隨時(shí)間的變化均符合零級(jí)動(dòng)力學(xué)反應(yīng)，53.33～66.66、66.66～79.99、79.99～93.32 kPa及常壓下其色澤的變化速率分別為：dQ1/dt=-0.014 3，dQ2/dt=-0.106 8，dQ3/dt=-0.127 1，dQ4/dt=-0.146，-a/b值變化的動(dòng)力學(xué)模型分別為：-a1/b1=-0.014 3t+1.370 8，-a2/b2=-0.010 68t+1.456 8，-a3/b3=-0.127 1t+1.452 1，-a4/b4=-0.146t+1.519 3。

2.2 基于色澤參數(shù)-a/b值變化的西蘭花葉綠素質(zhì)量分?jǐn)?shù)預(yù)測(cè)模型的建立

色澤參數(shù)-a/b值與葉綠素質(zhì)量分?jǐn)?shù)之間的關(guān)系見(jiàn)圖3。由圖3可以看出，不同壓力條件下，西蘭花-a/b值與葉綠素質(zhì)量分?jǐn)?shù)的變化均呈線性正相關(guān)； 經(jīng)回歸分析，53.33～66.66、66.66～79.99、79.99～93.32 kPa及常壓條件下建立基于西蘭花-a/b值變化的葉綠素質(zhì)量分?jǐn)?shù)的預(yù)測(cè)模型分別為：Cchl1=0.598 5（-a1/b1）-0.108 7（R2=0.997 6），Cchl2=0.608 7（-a2/b2）-0.149 3（R2=0.981 7），Cchl3=0.559 3（-a3/b3）-0.083 7（R2=0.997 7），Cchl4=0.531 9 （-a4/b4）-0.035 1 （R2=0.986 9）。

2.3 模型驗(yàn)證

根據(jù)微真空條件下-a/b值的動(dòng)力學(xué)模型，可得西蘭花采后在 1、8、15、22、29、36、43、50 d 時(shí)-a/b 的值，即預(yù)測(cè)值。與實(shí)測(cè)值進(jìn)行比較，見(jiàn)圖4。經(jīng)相關(guān)性分析得，53.33～66.66、66.66～79.99、79.99～93.32 kPa及常壓下回歸方程分別為：y1=0.975 1x1+0.034 8，y2=0.926 3x2+0.084 5，y3=0.940 9x3+0.067，y4=0.954 6x4+0.075 5， 復(fù)相關(guān)系數(shù)分別為：R12=0.973 6，R22=0.949 3，R32=0.982 1，R42=0.993 5。證明應(yīng)用西蘭花色澤參數(shù)動(dòng)力學(xué)模型能較好的預(yù)測(cè)-a/b值的變化。

圖3 西蘭花-a/b值與葉綠素質(zhì)量分?jǐn)?shù)的關(guān)系圖Fig.3 Diagram between broccoli-a/b value and chlorophyll content

圖4 －a/b值預(yù)測(cè)模型預(yù)測(cè)值與測(cè)量值的相關(guān)性Fig.4 Correlation between actual and predicted values of-a/b value

將-a/b預(yù)測(cè)值代入微真空條件下基于西蘭花-a/b值變化的葉綠素質(zhì)量分?jǐn)?shù)預(yù)測(cè)模型中，得到西蘭花葉綠素質(zhì)量分?jǐn)?shù)的預(yù)測(cè)值，將其與實(shí)測(cè)值進(jìn)行比較，見(jiàn)圖5。

圖5 葉綠素模型預(yù)測(cè)值與測(cè)量值的相關(guān)性Fig.5 Correlation between actual and predicted values of chlorophyll content in broccoli

經(jīng)相關(guān)性分析得，53.33～66.66、66.66～79.99、79.99～93.32 kPa及常壓下回歸方程分別為：y1=0.958 6x+0.027 6，y2=0.937 8x+0.035 1，y3=0.961 3x+0.024 6，y4=0.955 4x+0.029 8，復(fù)相關(guān)系數(shù)分別為：R12=0.968 7，R22=0.98，R32=0.989 57，R42=0.989 6。 證明不同微真空壓力及常壓條件下，基于-a/b值變化的西蘭花葉綠素質(zhì)量分?jǐn)?shù)預(yù)測(cè)模型均能很好地預(yù)測(cè)葉綠素質(zhì)量分?jǐn)?shù)的變化。

3 結(jié)語(yǔ)

從葉綠素的降解過(guò)程來(lái)看，葉綠素質(zhì)量分?jǐn)?shù)呈一直下降的趨勢(shì)，且貯藏前期西蘭花葉綠素降解速率較慢，后期降解速率加快。其原因可能是在西蘭花采摘的過(guò)程中，由于產(chǎn)生了較多的機(jī)械傷，果蔬自身的修復(fù)機(jī)制啟動(dòng)，這種機(jī)制通常是通過(guò)自身旺盛的代謝功能促進(jìn)傷口的愈合，在此過(guò)程中會(huì)引起內(nèi)源抗氧化酶活性的升高，而酶活的升高降低了引起葉綠素氧化漂白的活性氧代謝，從而抑制了葉綠素的降解。故西蘭花在貯藏前期其葉綠素降解速率較慢。

不同壓力條件下，西蘭花在采后貯藏過(guò)程中，隨著貯藏時(shí)間的延長(zhǎng)，葉綠素質(zhì)量分?jǐn)?shù)及-a/b值均呈逐漸下降的趨勢(shì)，由試驗(yàn)數(shù)據(jù)的分析結(jié)果可得：53.33～66.66、66.66～79.99、79.99～93.32 kPa 及常壓條件下葉綠素、-a/b值與貯藏時(shí)間之間的動(dòng)力學(xué)模型分別為=-0.008 5t+0.945 4，=-0.070 6t+0.756 4，Cchl3=-0.071 4t+0.726 4，Cchl4=-0.081 2t0.782 9和-a1/b1=-0.014 3t+1.370 8，-a2/b2=-0.010 68t+1.456 8，-a3/b3=-0.127 1t+1.452 1，-a4/b4=-0.146t+1.519 3，均符合零級(jí)動(dòng)力學(xué)模型；通過(guò)線性回歸分析，建立基于色澤參數(shù)-a/b值變化的葉綠素質(zhì)量分?jǐn)?shù)預(yù)測(cè)模型為：=0.598 5（-a/b）-0.108 7（R2=0.997 6），Cchl2=0.608 7（-a/b）-0.149 3（R2=0.981 7），Cchl3=0.559 3（-a/b）-0.083 7 （R2=0.997 7），Cchl4=0.531 9 （-a/b）-0.035 1（R2=0.986 9）；進(jìn)一步通過(guò)模型的驗(yàn)證，證明該預(yù)測(cè)模型能較好的反映貯藏過(guò)程中葉綠素質(zhì)量分?jǐn)?shù)的變化，為進(jìn)一步應(yīng)用無(wú)損檢測(cè)技術(shù)預(yù)測(cè)微真空條件下西蘭花的品質(zhì)變化提供了理論依據(jù)。

[1]馬海燕，張慜，孫金才.臭氧水和超聲波協(xié)同作用在速凍西蘭花中的應(yīng)用研究[J].食品與生物技術(shù)學(xué)報(bào)，2010，29（4）：538-542.MA Haiyan，ZHANG Min，SUN Jincai.The application study of ultrasonic-ozone combined treatments in quick-frozen broccoli[J].Journal of Food Science and Biotechnology，2010，29（4）：538-542.（in Chinese）

[2]沈蓮清，蘇光耀，王奎武.西蘭花種子中硫苷酶解產(chǎn)物蘿卜硫素的提純與抗腫瘤的體外試驗(yàn)研究[J].中國(guó)食品學(xué)報(bào)，2008，8（5）：16-21.SHEN Lianqing，SU Guangyao，WANG Kuiwu.Studies on purification and antitumor in vitro experiment of sulforaphane enzymolyzed from glucosinolate in broccoli seed[J].Journal of Chinese Institute of Food Science and Technology，2008，8（5）：16-21.（in Chinese）

[3]張素文，張慜，孫金才.水分含量對(duì)西蘭花玻璃化轉(zhuǎn)變溫度的影響[J].食品與生物技術(shù)學(xué)報(bào)，2008，27（3）：28-32.ZHANG Suwen，ZHANG Min，SUN Jincai.Effect of water content on the glass temperature of broccoli[J].Journal of Food Science and Biotechnology，2008，27（3）：28-32.（in Chinese）

[4]高持都，胡小松，廖小軍，等.高壓二氧化碳技術(shù)對(duì)葉綠素酶活性的影響研究[J].中國(guó)食品學(xué)報(bào)，2010，10（1）：55-60.GAO Chidu，HU XiaoSong，LIAO Xiaojun，et al.Effect of high pressure dioxide on chlorophyllase activity[J].Journal of Chinese Institute of Food Science and Technology，2010，10（1）：55-60.（in Chinese）

[5]Tijskens L M M，Schijvens E P H M，Biekman E S A.Modelling the change in colour of broccoli and green beans during blanching[J].Innovative Food Science and Emerging Technologies，2001（2）：303-313.

[6]Tian M S，Davies L，Downs C G，et al.Effects of floret maturity，cytokinin and ethylene on broccoli yellowing after harvest[J].Postharvest Biology and Technology，1995，6（2）：29-40.

[7]Tian M S，Woolf A B，Bowen J H，et al.Changes in color and chlorophyll fluorescence of broccoli florets following hot water treatment[J].Journal of the American Society for Horticultural Science，1996，121（2）：310-313.

[8]Jhas N，Chopra S，Kingsly A R P.Modeling of color values for nondestructive evaluation of maturity of mango[J].Journal of Food Engineering，2007，78（1）：22-26.

[9]Godoy C，Monterubbianesi G，Tognetti J.Analysis of highbush blueberry （Vaccinium corymbosum L.） fruit growth with exponential mixed models[J].Scientia Horticulturae，2008，115（4）：368-376.

[10]Goncalves E M，Pinheiro J，Abreu M，et al.Modelling the kinetics of peroxidase inactivation，colour and texture changes of pumpkin（Cucurbita maxima L.） during blanching[J].Journal of Food Engineering，2007，81（4）：693-701.

[11]徐鳳，楊震峰，裴嬌艷.基于顏色參數(shù)變化的青花菜葉綠素含量預(yù)測(cè)模型[J].食品科學(xué)，2011，32（13）：54-57.XU Feng，YANG Zhenfeng，PEI Jiaoyan.A prediction model for chlorophyll content in post-harvest broccoli based on color parameter changes[J].Food Science，2011，32（13）：54-57.（in Chinese）

[12]刁恩杰，李向陽(yáng)，丁曉雯.脫水菠菜貯藏過(guò)程中顏色變化動(dòng)力學(xué)[J].農(nóng)業(yè)工程學(xué)報(bào)，2010，26（8）：350-354.DIAO Enjie，LI Xiangyang，DING Xiaowen.Kinetics of colour changes of dehydrated spinaches during storage[J].Transactions of the Chinese Society of Agricultural Engineering，2010，26（8）：350-354.（in Chinese）

[13]劉春泉，卓成龍，宋江峰.毛豆仁低溫凍藏期葉綠素和色澤降解動(dòng)力學(xué)[J].江蘇農(nóng)業(yè)學(xué)報(bào)，2011，27（4）：858-862.LIU Chunquan，ZHUO Chenglong，SONG Jiangfeng.Degradation kinetics of color and chlorophyll in frozen edamames[J].Jiangsu Journal of Agricultural Sciences，2011，27（4）：858-862.（in Chinese）

[14]任珂，屠康，潘磊慶.青花菜貯藏期間顏色變化動(dòng)力學(xué)模型的建立[J].農(nóng)業(yè)工程學(xué)報(bào)，2005，21（8）：146-150.REN Ke，TU Kang，PAN Leiqing.Modeling of the kinetics of color change of broccoli during storage[J].Transactions of The Chinese Society of Agricultural Engineering，2005，21（8）：146-150.（in Chinese）

[15]王世清.微真空貯藏設(shè)施[P].中國(guó)專利：CN 1530290A，2004-09-22.

[16]周莎莎，吳愛(ài)現(xiàn)，李文香.微真空貯藏條件對(duì)萊陽(yáng)茌梨保鮮效果的影響[J].食品科學(xué)，2011，32（12）：320-323.ZHOU Shasha，WU Aixian，LI Wenxiang.Effect of micro vacuum storage on the fresh-keeping of laiyang pear[J].Food Science，2011，32（12）：320-323.（in Chinese）

[17]曹建康，姜微波，趙玉梅.果蔬采后生理生化實(shí)驗(yàn)指導(dǎo)[M].中國(guó)輕工業(yè)出版社，2007.CAO Jiankang，JIANG Weibo，ZHAO Yumei.Experiment guidance of postharvest physiology and biochemistry of fruits and vegetables[M].China Light Industry Press，2007.（in Chinese）

[18]裴嬌艷，楊震峰，徐鳳.基于綠熟番茄果實(shí)L*值的番茄紅素含量預(yù)測(cè)模型建立[J].食品科學(xué)，2010，31（10）：296-300.PEI Jiaoyan，YANG Zhenfeng，XU Feng.Prediction model，based on L*value，of lycopene content of mature green tomato[J].Food Science，2010，31（10）：296-300.（in Chinese）

[19]Ku V V，Wills R B H.Effect of 1-methylcyclopropene on the storage life of broccoli[J].Posthary Biol Technol，1999，17（2）：127-132.

[20]王冬梅，江連洲，趙曉燕.干燥過(guò)程中甘藍(lán)葉綠素降解動(dòng)力學(xué)的研究[J].食品科學(xué)，2012，33（17）：64-67.WANG Dongmei，JIANG Lianzhou，ZHAO Xiaoyan.Kinetic model of chlorophyll degradation in cabbage during drying process[J].Food Science，2012，33（17）：64-67.（in Chinese）