呂興龍　劉強(qiáng)　孫國(guó)昌　彭海根　毛青鐘　葉芙蓉

摘要：為提高黃酒發(fā)酵過(guò)程常規(guī)指標(biāo)的檢測(cè)效率，提出采用漫反射方式采集黃酒醪液的近紅外光譜，并應(yīng)用近紅外技術(shù)結(jié)合偏最小二乘方法（PLS）建立一種能夠同時(shí)測(cè)定黃酒醪液酒精度、總酸和總糖的定量分析模型。結(jié)果顯示：優(yōu)化后3種指標(biāo)模型的決定系數(shù)（R）分別為0.99、0.91和0.98，交互驗(yàn)證標(biāo)準(zhǔn)偏差（RMSECV）分別為0.42、0.32和1.19;模型進(jìn)行外部驗(yàn)證時(shí)，預(yù)測(cè)標(biāo)準(zhǔn)偏差（RMSEP）分別為0.45、0.35和2.01，預(yù)測(cè)平均相對(duì)偏差分別為3.37%、4.57%和4.66%，滿(mǎn)足黃酒發(fā)酵過(guò)程中對(duì)樣品分析檢測(cè)精度的要求。研究結(jié)果表明該方法可以用于黃酒醪液常規(guī)指標(biāo)的快速分析，對(duì)于解決黃酒自動(dòng)化生產(chǎn)中指標(biāo)檢測(cè)問(wèn)題具有積極意義。

關(guān)鍵詞：近紅外光譜;黃酒;醪液;快速檢測(cè)

中圖分類(lèi)號(hào)：TP274.5

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1674–5124（2019）03–0075–05

Study on rapid detection of physiochemical index of Chinese rice wine fermentation process by near infrared spectroscopy

L? Xinglong1， LIU Qiang2， SUN Guochang1， PENG Haigen3， MAO Qingzhong1， YE Furong1

（1. Kuaijishan Shaoxing Wine Co.， Ltd.， Shaoxing 312030， China; 2. Sichuan Insitute of Product Quality Supervision and Inspection， Chengdu 610100， China; 3. Sichuan Vspec Technology Co.， Ltd.， Chengdu 610041， China）

Abstract： In order to improve the conventional indicators detection efficiency for Chinese rice wine fermentation process， collected the diffuse reflection spectrum and the quantitative analysis models were established that simultaneous determined the alcohol content， total acidity and total sugar of Chinese rice wine by near-infrared （NIR） spectroscopy with PLS method. The results showed that determination coefficient （R） of the models were 0.99， 0.91 and 0.98， root mean square error of cross validation （RMSECV） were 0.42， 0.32 and 1.19. When the models were externally validated， root mean square error of prediction （RMSEP） were 0.45， 0.35 and 2.01， and the average relative deviations were 3.37%， 4.57% and 4.66% respectively， which met the requirement of sample analysis and detection accuracy in Chinese rice wine fermentation. The method can be used for rapid analysis of routine index of Chinese rice wine mash， and has positive significance for solving the problem of physiochemical index detection in automatic production of rice wine.

Keywords： near-infrared （NIR）; Chinese rice wine; mash; rapid detection

0 引言

黃酒具有悠久的歷史傳統(tǒng)和深厚的文化底蘊(yùn)，被譽(yù)為中國(guó)的“國(guó)粹”酒，與啤酒、葡萄酒并稱(chēng)世界三大古酒[1]。近年來(lái)，隨著人們生活水平的提高和消費(fèi)理念的轉(zhuǎn)變，對(duì)酒的需求由“單純嗜好”向“營(yíng)養(yǎng)保健”轉(zhuǎn)變，即高度、烈性的不良飲酒觀念日益為人們所摒棄，黃酒的低度、營(yíng)養(yǎng)、保健的優(yōu)勢(shì)逐漸顯現(xiàn)，慢慢被人們接受。黃酒的產(chǎn)量也有逐年上升的趨勢(shì)，國(guó)家統(tǒng)計(jì)局?jǐn)?shù)據(jù)顯示，黃酒在國(guó)內(nèi)外市場(chǎng)具有廣泛的發(fā)展空間，根據(jù)“十五”計(jì)劃和2015年規(guī)劃，黃酒行業(yè)的總產(chǎn)量從2000年到2015年要從145萬(wàn)噸增長(zhǎng)至250～280萬(wàn)噸[2]。隨著黃酒需求量的不斷增長(zhǎng)，以及計(jì)算機(jī)技術(shù)、自動(dòng)化技術(shù)和信息技術(shù)的發(fā)展，黃酒企業(yè)在生產(chǎn)綜合自動(dòng)化方面進(jìn)行了積極的嘗試，并取得了長(zhǎng)足進(jìn)展，如會(huì)稽山黃酒股份有限公司在紹興柯西建成了一條年產(chǎn)2萬(wàn)千升的黃酒自動(dòng)化生產(chǎn)線(xiàn);紹興女兒紅酒廠與江南大學(xué)合作，建立了可實(shí)現(xiàn)前酵過(guò)程自動(dòng)化的自控系統(tǒng)等[3]。

黃酒產(chǎn)量的提高和生產(chǎn)綜合自動(dòng)化的實(shí)施對(duì)黃酒發(fā)酵過(guò)程醪液理化指標(biāo)檢測(cè)提出全新的要求。黃酒發(fā)酵過(guò)程中，酒精度、總酸和總糖3項(xiàng)指標(biāo)影響黃酒的品質(zhì)、風(fēng)味和口感[4]，需要檢測(cè)上述指標(biāo)含量的變化，標(biāo)準(zhǔn)方法采用GBT13662—2008《黃酒》標(biāo)準(zhǔn)[5]，方法需要多種儀器和化學(xué)試劑，樣品需要預(yù)處理，步驟繁瑣、費(fèi)時(shí)、費(fèi)力，顯然不能滿(mǎn)足現(xiàn)代黃酒生產(chǎn)對(duì)黃酒醪液指標(biāo)快速分析日益增長(zhǎng)的需求。如會(huì)稽山紹興酒股份有限公司實(shí)施自動(dòng)化釀酒之后，發(fā)酵時(shí)間由傳統(tǒng)的60～90天縮短到20～30天，發(fā)酵罐數(shù)量也較之前有所增加，因而每天樣品檢測(cè)量明顯增多。因此，黃酒發(fā)酵過(guò)程中快速分析技術(shù)研究已引起行業(yè)的廣泛關(guān)注。

近紅外技術(shù)具有快速、準(zhǔn)確和幾乎不需要樣品制備等優(yōu)點(diǎn)[6]，非常適合快速分析黃酒發(fā)酵生產(chǎn)過(guò)程中常規(guī)指標(biāo)的含量。蔣詩(shī)泉[7]、胡小邦[8]、周揚(yáng)[9]、謝廣發(fā)[10]，SHENF[11]等分別在黃酒近紅外定性或定量分析進(jìn)行了研究，但是上述研究主要應(yīng)用透射技術(shù)于成品酒酒齡的識(shí)別或常規(guī)指標(biāo)含量的檢測(cè)。應(yīng)用近紅外技術(shù)快速檢測(cè)黃酒發(fā)酵過(guò)程中醪液常規(guī)指標(biāo)的研究尚未見(jiàn)報(bào)道。

本文提出采用積分球結(jié)合漫反射方式測(cè)定黃酒醪液光譜，研究應(yīng)用近紅外技術(shù)同時(shí)快速測(cè)定黃酒醪液酒精度、總酸和總糖常規(guī)指標(biāo)含量，并對(duì)相關(guān)校正模型進(jìn)行綜合評(píng)價(jià)，以期解決目前黃酒發(fā)酵過(guò)程中常規(guī)指標(biāo)檢測(cè)中存在效率低的問(wèn)題，為黃酒發(fā)酵過(guò)程中的黃酒質(zhì)量快速檢測(cè)提供準(zhǔn)確的數(shù)據(jù)支撐，提高產(chǎn)品質(zhì)量，具有重要的實(shí)際意義。

1 材料與方法

1.1 樣品收集

黃酒發(fā)酵醪液樣品由會(huì)稽山紹興酒股份有限公司提供，生產(chǎn)日期為2018年1月-2018年6月，數(shù)量為711個(gè)。樣品涵蓋黃酒前發(fā)酵和后發(fā)酵整個(gè)發(fā)酵過(guò)程。為了保證樣品狀態(tài)的均勻性，取樣時(shí)，使用紗布過(guò)濾后再裝樣，備用。

1.2性質(zhì)測(cè)定常規(guī)方法測(cè)定性質(zhì)前，先將樣品進(jìn)行離心，取上清液，采用GBT13662—2008《黃酒》標(biāo)準(zhǔn)分析方法對(duì)所有樣品的酒精度、總酸和總糖性質(zhì)進(jìn)行測(cè)定。樣品各理化指標(biāo)的參考值統(tǒng)計(jì)結(jié)果見(jiàn)表1，通過(guò)統(tǒng)計(jì)分析，樣品上述性質(zhì)指標(biāo)分布范圍涵蓋黃酒發(fā)酵過(guò)程中性質(zhì)所有指標(biāo)范圍，具有代表性。

1.3 光譜采集

取樣時(shí)，黃酒醪液樣品采用紗布進(jìn)行過(guò)濾，屬于漿狀物，與酸奶形態(tài)類(lèi)似;祝義偉等[12]采用漫反射方法采集得到酸奶近紅外光譜，并建立預(yù)測(cè)效果很好的非脂乳固體近紅外校正模型，因而本文提出采用積分球結(jié)合漫反射方式采集黃酒醪液樣品光譜的方法。光譜采集過(guò)程中，由于黃酒醪液中的酵母菌仍然處于活性狀態(tài)，漫反射方法相較先將樣品離心然后采用透射方式采集光譜的方法，明顯縮短光譜采集時(shí)間，因而可以避免時(shí)間過(guò)長(zhǎng)導(dǎo)致樣品性質(zhì)變化。

在美國(guó)Galaxy公司生產(chǎn)的QuasIR3000傅里葉變換近紅外光譜儀上，采集黃酒發(fā)酵醪液的近紅外光譜。儀器參數(shù)設(shè)置：1）光譜掃描范圍：10000cm–1～4000cm–1;2）分辨率：16cm–1;3）光譜掃描次數(shù)：64次;4）采用積分球漫反射采集光譜，以?xún)x器內(nèi)置背景做參比。將黃酒醪液樣品倒入樣品杯中時(shí)，盡量不要產(chǎn)生氣泡。

圖1是黃酒半成品醪液樣品漫反射吸光度光譜圖，譜圖平滑，吸光度處于正常范圍內(nèi)，且在近紅外譜區(qū)合頻區(qū)域和倍頻區(qū)域都有明顯的吸收峰，進(jìn)一步表明本文采用積分球結(jié)合漫反射采集黃酒醪液樣品的方法可行。

1.4 模型評(píng)價(jià)參數(shù)介紹

在應(yīng)用近紅外技術(shù)建立定量模型過(guò)程中，需要結(jié)合不同參數(shù)評(píng)價(jià)模型，主要包括如下：1）決定系數(shù)R決定系數(shù)主要用來(lái)判斷定量校正模型與待測(cè)組分的線(xiàn)性關(guān)系，計(jì)算方法如下：

式中：yi，actual——第i個(gè)樣品參考方法的測(cè)定值;

y?i，actual——校正集（或驗(yàn)證集）全部樣品參考方法測(cè)定的平均值;

yi，predicted——校正集（或驗(yàn)證集）驗(yàn)證過(guò)程中第i個(gè)樣品的預(yù)測(cè)值;

n——校正集（或驗(yàn)證集）的樣品數(shù)量。

在樣品化學(xué)成分范圍相同的前提下，R越接近1，模型回歸（或預(yù)測(cè)）結(jié)果越好。

2）交互驗(yàn)證標(biāo)準(zhǔn)偏差（RMSECV）和預(yù)測(cè)標(biāo)準(zhǔn)偏差（RMSEP）

建立近紅外光譜定量校正模型時(shí)，樣品分為校正集和驗(yàn)證集，校正集是用來(lái)建立定量校正模型，驗(yàn)證集是用來(lái)對(duì)模型進(jìn)行驗(yàn)證，未參與校正模型的建立過(guò)程，通過(guò)比較參考值與模型預(yù)測(cè)值的差異，來(lái)判斷模型預(yù)測(cè)能力的好壞。

建模過(guò)程中，可以通過(guò)交互驗(yàn)證的方法對(duì)模型進(jìn)行驗(yàn)證。其原理為：模型驗(yàn)證過(guò)程中，每次從校正集中取出一個(gè)或多個(gè)樣品作為臨時(shí)驗(yàn)證樣品，以其余的樣品進(jìn)行建模，然后對(duì)這一個(gè)或多個(gè)樣品進(jìn)行預(yù)測(cè)，如此循壞，則會(huì)分別得到每個(gè)樣品的模型交叉預(yù)測(cè)值，最后，以交叉預(yù)測(cè)值與參考方法值誤差平方和的均方根值為RMSECV，其計(jì)算公式如下：

式中：yi，actual——第i個(gè)樣品參考方法的測(cè)量值;

yi，predicted——校正集交互驗(yàn)證中第i個(gè)樣品的預(yù)測(cè)值;

n——校正集的樣品個(gè)數(shù)。

預(yù)測(cè)標(biāo)準(zhǔn)偏差（RMSEP）是計(jì)算模型得出的預(yù)測(cè)值與參考值之間的誤差平方和的均方根值，其計(jì)算公式如下：

式中：yi，actual——第i個(gè)樣品參考方法的測(cè)量值;

yi，predicted——驗(yàn)證集預(yù)測(cè)過(guò)程中第i個(gè)樣品的模型方法預(yù)測(cè)值;

m——驗(yàn)證集的樣品個(gè)數(shù)。

本文采用交互驗(yàn)證和驗(yàn)證集驗(yàn)證相結(jié)合方法驗(yàn)證模型效果。對(duì)于同一批次樣本，模型RMSECV和RMSEP值越小說(shuō)明模型效果越好。

2 結(jié)果與分析

2.1 建立近紅外定量分析模型

將黃酒醪液樣品的近紅外圖譜和參考值數(shù)據(jù)一一對(duì)應(yīng)，采用VSPEC公司開(kāi)發(fā)VModel近紅外建模軟件分別建立酒精度、總酸、總糖的定量分析模型。

在近紅外定量分析模型建立過(guò)程中，有必要對(duì)光譜數(shù)據(jù)進(jìn)行光譜預(yù)處理。常規(guī)的光譜預(yù)處理方法包括多元散射校正（MSC）、標(biāo)準(zhǔn)正則變換（SNV）、一階導(dǎo)數(shù)、二階導(dǎo)數(shù)和平滑處理等。MSC和SNV具有消除顆粒分布不均勻及顆粒大小產(chǎn)生的散射影響，在固體漫反射和漿狀物透（反）射光譜中應(yīng)用較廣泛[13];導(dǎo)數(shù)處理既可以消除基線(xiàn)偏移，還可以起到一定的放大和分離重疊信息的作用，由于噪聲信號(hào)也被放大，因此通常在微分之前需要對(duì)光譜數(shù)據(jù)做平滑處理。

建模時(shí)，先將光譜全譜段平均分成10等分，采用不同的單個(gè)譜段分別建立不同校正模型，并采用交互驗(yàn)證方式結(jié)合上述評(píng)價(jià)參數(shù)對(duì)各個(gè)模型效果進(jìn)行評(píng)價(jià)，確定最優(yōu)模型的譜段范圍，然后擴(kuò)充譜段數(shù)量和譜段范圍，重新建立校正模型并進(jìn)行模型評(píng)價(jià)，直到找出最優(yōu)的譜段或譜段組合。

偏最小二乘方法（PLS）在主成分分析（PCA）的基礎(chǔ)上，對(duì)光譜矩陣和濃度矩陣同時(shí)進(jìn)行分解，并在分解時(shí)考慮兩者相互之間的關(guān)系，加強(qiáng)對(duì)應(yīng)計(jì)算關(guān)系，從而保證獲得最佳的校正模型。偏最小二乘方法（PLS）是多元線(xiàn)性回歸、典型相關(guān)分析和主成分分析的完美結(jié)合，而且模型效果明顯優(yōu)于其他方法，這也是PLS在光譜多元校正分析中應(yīng)用最為廣泛的原因之一[13]。

綜上所述，首先采用單個(gè)光譜預(yù)處理或預(yù)處理方法組合的方式對(duì)樣品光譜進(jìn)行預(yù)處理，結(jié)合上述建模譜段優(yōu)選方法，采用偏最小二乘法（PLS），分別建立酒精度、總酸和總糖定量校正模型，并使用模型評(píng)價(jià)參數(shù)對(duì)各個(gè)模型進(jìn)行評(píng)價(jià)，從而確定最優(yōu)的預(yù)處理方法。經(jīng)過(guò)優(yōu)化后，黃酒醪液酒精度、總酸和總糖定量模型交互驗(yàn)證結(jié)果和相應(yīng)的建模參數(shù)詳見(jiàn)表2。

表2顯示，模型決定系數(shù)R都大于0.90，即模型相關(guān)性很好，交互驗(yàn)證標(biāo)準(zhǔn)偏差（RMSECV）都很小，模型效果很好，且主因子數(shù)合理。

黃酒醪液酒精度、總酸和總糖定量模型校正集參考值和預(yù)測(cè)值的散點(diǎn)圖如圖2、圖4和圖6所示;主因子數(shù)與RMSECV的變化趨勢(shì)圖如圖3、圖5、圖7所示。

2.2模型外部驗(yàn)證結(jié)果為了進(jìn)一步驗(yàn)證模型效果，隨機(jī)挑選30個(gè)具

有代表性的樣品作為外部驗(yàn)證樣品，并應(yīng)用所建立的黃酒醪液酒精度、總酸和總糖定量模型預(yù)測(cè)樣品上述指標(biāo)的含量值，與標(biāo)準(zhǔn)分析方法值進(jìn)行對(duì)比，結(jié)果如表3所示。

黃酒醪液酒精度、總酸和總糖定量模型驗(yàn)證集參考值和預(yù)測(cè)值的散點(diǎn)圖如圖8、圖9和圖10所示。

通過(guò)對(duì)表3和圖8～圖10分析，本文挑選的驗(yàn)證樣品在各個(gè)指標(biāo)濃度范圍內(nèi)分布均勻，即具有代表性，可以作為模型外部驗(yàn)證樣品;模型預(yù)測(cè)標(biāo)準(zhǔn)偏差（RMSEP）分別為0.45、0.35和2.01，RMSEP和RMSECV相差很小，即模型的穩(wěn)定性很好，且預(yù)測(cè)平均相對(duì)偏差分別為3.37%、4.57%和4.66%，與標(biāo)準(zhǔn)方法重復(fù)性要求5%相近，進(jìn)一步說(shuō)明本文建立的黃酒醪液酒精度、總酸和總糖近紅外定量模型準(zhǔn)確度滿(mǎn)足生產(chǎn)過(guò)程中對(duì)樣品分析檢測(cè)誤差的要求。

3 結(jié)束語(yǔ)

本文提出應(yīng)用漫反射方式采集黃酒醪液光譜，結(jié)合近紅外光譜分析技術(shù)快速測(cè)定黃酒發(fā)酵過(guò)程中酒精度、總酸和總糖指標(biāo)含量的方法可行;黃酒醪液酒精度、總酸和總糖近紅外定量模型經(jīng)過(guò)優(yōu)化后，R分別為0.99、0.91和0.98，RMSECV分別為0.42、0.32和1.19，且模型通過(guò)外部驗(yàn)證，方法精確度滿(mǎn)足黃酒發(fā)酵過(guò)程中常規(guī)指標(biāo)檢測(cè)的要求，可以用于黃酒自動(dòng)化發(fā)酵過(guò)程中醪液常規(guī)指標(biāo)的快速分析，具有積極意義。

參考文獻(xiàn)

[1]駱高遠(yuǎn).紹興黃酒文化與旅游[J].經(jīng)濟(jì)地理，2005（9）：749-752.

[2]楊國(guó)軍.中國(guó)黃酒業(yè)調(diào)研報(bào)告[J].中國(guó)釀造，2005（4）：1-5.

[3]王培璇.中國(guó)黃酒出口競(jìng)爭(zhēng)力研究[D].無(wú)錫：江南大學(xué)，2014.

[4]江超，林峰.陳年紹興黃酒的成分分析與品質(zhì)鑒定[J].食品與發(fā)酵工業(yè)，2009（10）：119-123.

[5]黃酒：GBT13662-2008[S].北京：中國(guó)標(biāo)準(zhǔn)出版社，2008.

[6]陸婉珍主編.現(xiàn)代近紅外光譜分析技術(shù)[M].2版.北京：中國(guó)石化出版社，2006.

[7]蔣詩(shī)泉，周興才，蔣詩(shī)平，等.基于PCA和LS-SVM的傅里葉變換近紅外光譜的黃酒酒齡的鑒別模型[J].光譜實(shí)驗(yàn)室，2012，29（2）：806-811.

[8]胡小邦，呂進(jìn)，劉輝軍，等.基于近紅外透射光譜的黃酒酒精度、酸度檢測(cè)模型研究[J].光電工程，2010，37（10）：122-127.

[9]周揚(yáng)，呂進(jìn)，劉鐵兵，等.小波域高斯混合模型方差估計(jì)近紅外降噪方法[J].光電工程，2011，38（8）：96-100.

[10]謝廣發(fā)，徐榕，樊阿萍，等.近紅外技術(shù)在黃酒理化指標(biāo)快速檢測(cè)中的應(yīng)用[J].中國(guó)釀造，2011，30（11）：182-185.

[11] SHEN F， NIU X Y， YANG D T， et al. Determination of amino acids in Chinese rice wine by fourier transform near-infrared spectroscopy[J]. Agric Food Chem， 2010， 58： 9809-9816.

[12]祝義偉，肖琳，費(fèi)華熙，等.酸奶中非脂乳固體的近紅外定量分析模型建立[J].中國(guó)乳品工業(yè)，2014，42（7）：39-41.

[13]褚小立主編.化學(xué)計(jì)量學(xué)方法和分子光譜分析技術(shù)[M].北京：化學(xué)工業(yè)出版社，2011：4.

（編輯：徐柳）