采用近紅外光譜快速測(cè)定法對(duì)生物柴油的成分(脂肪酸甲酯、單甘酯、二甘酯、三甘酯和甘油)進(jìn)行了研究。采用氣相色譜方法獲得其成分的基礎(chǔ)數(shù)據(jù)，通過(guò)偏最小二乘方法與近紅外光譜數(shù)據(jù)進(jìn)行回歸運(yùn)算，分別建立以文冠果油生物柴油為例的單原料油校正模型及多種原料油生物柴油的混合校正模型，并以花椒油生物柴油為例考察了校正模型的適用性。結(jié)果表明：通過(guò)偏最小二乘方法可以建立適合多種原料油生物柴油的通用校正模型。對(duì)于新型生物柴油，向校正集中添加10個(gè)以上樣本，擴(kuò)充校正模型后，便可較為準(zhǔn)確地測(cè)定這類新生物柴油樣本的成分含量。此方法分析速度快、成本低、操作便捷、重復(fù)性好，適合于生物柴油生產(chǎn)過(guò)程的中間控制分析。

　　【關(guān)鍵詞】近紅外光譜;生物柴油;化學(xué)計(jì)量學(xué);甲酯

　　1引言

　　動(dòng)植物油脂主要是各種脂肪酸甘油三酯的混合物，在與甲醇的反應(yīng)過(guò)程中會(huì)生成脂肪酸單甘酯、二甘酯等中間產(chǎn)物和副產(chǎn)物甘油，影響生物柴油的產(chǎn)率及品質(zhì)。因此，在生產(chǎn)過(guò)程中，需要及時(shí)準(zhǔn)確地測(cè)定脂肪酸甲酯、單甘酯、二甘酯、三甘酯和甘油的含量，以便精確控制反應(yīng)過(guò)程及產(chǎn)品質(zhì)量。目前，對(duì)這些組分的測(cè)定均采用氣相色譜方法[3,4]，分析前需對(duì)樣品進(jìn)行衍生化(如硅烷化)處理，定性和定量過(guò)程較為復(fù)雜，操作難度大，分析時(shí)間長(zhǎng)，不適用于生物柴油生產(chǎn)過(guò)程的中間控制分析。

　　近紅外光譜法(NIR)不需要對(duì)樣品進(jìn)行預(yù)處理，操作簡(jiǎn)單,分析快速，非常適合油品的定量和定性分析[5]，石油化工業(yè)中得到了較為廣泛的應(yīng)用[6～10]。已有將NIR光譜用于測(cè)定生物柴油中甲酯、甲醇和甘油等成分的報(bào)道[11～14]，均得到了滿意結(jié)果。但這些研究所用樣本的原料相對(duì)單一，建模樣本較少，所測(cè)定的成分也不完全，且大都是實(shí)驗(yàn)室合成的樣品，不具有實(shí)際的工業(yè)應(yīng)用價(jià)值。

　　本研究從我國(guó)多個(gè)中試和工業(yè)裝置中收集了200余個(gè)生物柴油樣本，涉及6種原料油及多種加工工藝，其中包括來(lái)源復(fù)雜的地溝油和酸化油，樣本的覆蓋范圍廣，所得的結(jié)果更具參考性，建立的校正模型具有實(shí)際應(yīng)用價(jià)值。

　　2實(shí)驗(yàn)部分

　　2.1生物柴油樣品和儀器

　　收集了224個(gè)生物柴油樣品，從反應(yīng)器取樣后取其甲酯相，并蒸餾除去(120℃)其中的甲醇及少量的水，置于20mL密封小瓶中保存。其原料油來(lái)源分別為花椒油(47個(gè))、苦山杏油(16個(gè))、文冠果油(44個(gè))、地溝油(40個(gè))和酸化油(石煉廠43個(gè)、中糧集團(tuán)34個(gè))。

　　氣相色譜儀，程序升溫冷柱頭柱上進(jìn)樣系統(tǒng)，HT1不銹鋼毛細(xì)柱。光譜儀。

　　2.2基礎(chǔ)數(shù)據(jù)測(cè)定

　　采用氣相色譜方法[4]測(cè)定生物柴油成分的基礎(chǔ)數(shù)據(jù)，色譜條件：分別采用油酸甲酯、單油酸甘油酯、1,3�捕�油酸甘油酯、三油酸甘油酯作為定性和定量標(biāo)準(zhǔn)。采用雙內(nèi)標(biāo)法進(jìn)行定量分析，以正十五烷為內(nèi)標(biāo)測(cè)定甲酯含量，以三癸酸甘油酯為內(nèi)標(biāo)測(cè)定甘油、單甘酯、二甘酯和三甘酯的含量。

　　以空白樣品池為參比，采集生物柴油樣品的NIR譜圖，光譜范圍4000～10000cm-1，分辨率8cm-1，掃描次數(shù)128。

　　2.3數(shù)據(jù)分析

　　采用石油化工科學(xué)研究院編制的“化學(xué)計(jì)量學(xué)軟件3.0版”在PC計(jì)算機(jī)上處理光譜數(shù)據(jù)。用K�睸方法將生物柴油樣品分為校正集和驗(yàn)證集，其中校正集用于建立分析模型，驗(yàn)證集用于檢測(cè)模型的準(zhǔn)確性。采用馬氏距離、光譜殘差和最鄰近距離3個(gè)指標(biāo)判斷模型對(duì)未知樣本的適用性。

　　校正集的光譜經(jīng)一階微分處理后分別與氣相色譜測(cè)定的成分?jǐn)?shù)據(jù)(甘油、甲酯、單甘酯、二甘酯以及三甘酯的質(zhì)量分?jǐn)?shù))通過(guò)偏最小二乘方法(PLS1)進(jìn)行回歸運(yùn)算，建立校正模型，PLS所用的最佳主因子數(shù)由留一法交互驗(yàn)證所得的預(yù)測(cè)殘差平方和(PRESS值)確定。為考察校正模型的適用性和準(zhǔn)確性。本研究對(duì)不同種類的生物柴油樣品分別進(jìn)行建模與驗(yàn)證研究，包括：以文冠果油生物柴油為例的單原料油的校正模型與驗(yàn)證;多種原料油的通用校正模型的建立與驗(yàn)證;以花椒油生物柴油為例，研究了NIR光譜校正模型的適用性。

　　3結(jié)果與討論

　　3.1基礎(chǔ)數(shù)據(jù)獲取

　　采用2.2節(jié)的方法測(cè)定生物柴油成分的基礎(chǔ)數(shù)據(jù)。典型的生物柴油樣品色譜圖見(jiàn)圖1;224個(gè)生物柴油的成分分布見(jiàn)表1�？梢�(jiàn)樣本的覆蓋范圍較寬，具有較強(qiáng)的代表性。與標(biāo)準(zhǔn)方法[15]相比，此方法的準(zhǔn)確性和重復(fù)性較好。表1樣品集性質(zhì)統(tǒng)計(jì)(略)

　　3.2近紅外光譜解析

　　圖2是所收集的6類原料油生物柴油的NIR光譜圖。由圖2可見(jiàn)，光譜中主要包含CH倍頻和合頻的吸收信息，如8609和8316cm-1處的吸收是甲基CH的二級(jí)倍頻峰，7216和7108cm-1處的吸收是亞甲基CH的第一組合頻峰，5828和5766cm-1處的吸收分別是甲基和亞甲基CH的一級(jí)倍頻峰，6005cm-1處為甲酯中的CH吸收峰。

　　盡管6類原料油生物柴油的NIR光譜非常相近，但其成分上卻存在一定差異。例如，花椒油生物柴油主要是亞麻酸、亞油酸和油酸甲酯;苦山杏油生物柴油主要是油酸和亞油酸甲酯;酸化油和地溝油生物柴油則主要是多種脂肪酸甲酯的混合，具體視原料油種類而定。這種差異可以在其NIR光譜的主成分分析(PCA)得分圖上表現(xiàn)出來(lái)。圖3為6類52個(gè)有代表性生物柴油樣品的前3個(gè)主成分得分分布圖。從圖3可見(jiàn)，不同原料生產(chǎn)的生物柴油有聚類分布的傾向。由于酸化油原料來(lái)源復(fù)雜，所以其分布較寬。其余4種生物柴油的分布則相對(duì)集中。但從整體來(lái)看，盡管不同種類生物柴油的NIR光譜之間存在差異，但這些差異并非十分顯著，說(shuō)明脂肪酸基團(tuán)對(duì)光譜會(huì)產(chǎn)生一定的影響，但NIR光譜主要反映的是甲酯基團(tuán)的信息。因此，有可能建立適合多種原料油生物柴油的通用校正模型。

　　3.3校正模型的建立

　　3.3.1文冠果油單原料油模型與驗(yàn)證

　　將44個(gè)文冠果油生物柴油樣品分為校正集(36個(gè)樣品)和驗(yàn)證集(8個(gè)樣品)，考察了不同光譜區(qū)間對(duì)模型參數(shù)的影響，以及單種原料油所建模型的預(yù)測(cè)準(zhǔn)確性。

　　生物柴油在4000～10000cm-1范圍的光譜包含了OH鍵的一級(jí)倍頻、CH鍵的一級(jí)倍頻、二級(jí)倍頻及其合頻信息。從圖2可以看出，CH一級(jí)倍頻區(qū)(6000～6550cm-1)和部分CH第一合頻區(qū)(4550～4000cm-1)的吸光度值超過(guò)了1.5AU，這些區(qū)間的光譜吸收強(qiáng)度可能與成分濃度之間存在較強(qiáng)的非線性關(guān)系。另外，10000～9000cm-1區(qū)間的光譜存在較大的噪聲，且沒(méi)有顯著的特征吸收。若這些光譜區(qū)間參與模型建立，將會(huì)使校正模型預(yù)測(cè)能力下降。表2給出了利用全譜和去除這些光譜區(qū)間(即選取4550～5500cm-1和6550～9000cm-1區(qū)間)所建模型及其預(yù)測(cè)的對(duì)比結(jié)果�？梢钥闯觯�選用4550～5500cm-1和6550～9000cm-1光譜區(qū)間建立模型的預(yù)測(cè)能力明顯優(yōu)于全譜的結(jié)果。表2文冠果油生物柴油不同光譜區(qū)間的校正和預(yù)測(cè)結(jié)果比較(略)注(Note)：R2為校正集交互驗(yàn)證得到的相關(guān)系數(shù)(Correlationcoefficient);SECV為校正集交互驗(yàn)證得到的預(yù)測(cè)標(biāo)準(zhǔn)偏差(Standarderrorofcrossvalidation);SEP為驗(yàn)證集預(yù)測(cè)標(biāo)準(zhǔn)偏差(Standarderrorofprediction)。

　　從單種原料油的建模和預(yù)測(cè)結(jié)果可以看出，對(duì)測(cè)定生物柴油中的甘油、甲酯、單甘酯、二甘酯和三甘酯的含量，NIR光譜方法可以給出比較準(zhǔn)確的結(jié)果。

　　3.3.2多種原料油模型與驗(yàn)證

　　若能建立適合所有原料油來(lái)源的生物柴油全局校正模型(Globalmodels)，則對(duì)該方法的推廣和使用都將帶來(lái)很大的便利。將6種原料油來(lái)源的224個(gè)生物柴油樣品分為校正集(167個(gè)樣品)和驗(yàn)證集(57個(gè)樣品)，以考察建立全局模型的可行性。表3多種原料油生物柴油的校正和預(yù)測(cè)結(jié)果(略)

　　采用與文冠果油生物柴油建模相同的條件建立了多種原料油生物柴油的校正模型。表3給出了模型建立過(guò)程的交互驗(yàn)證結(jié)果和驗(yàn)證集的預(yù)測(cè)結(jié)果;圖4中A,B和C分別為甲酯、單甘酯和二甘酯的NIR光譜交互驗(yàn)證預(yù)測(cè)值與色譜方法的相關(guān)圖。與單原料的文冠果油結(jié)果相比，由于成分變得相對(duì)復(fù)雜，光譜之間的干擾更為顯著，建立模型所用的主因子數(shù)明顯增加，同時(shí)校正和預(yù)測(cè)結(jié)果都有不同程度的降低。盡管如此，仍得到了滿意的結(jié)果，完全滿足工業(yè)生產(chǎn)對(duì)過(guò)程分析誤差的要求。

　　上述結(jié)果表明，建立多種原料油來(lái)源的生物柴油全局校正模型是完全可行的，但隨著模型中原料油種類的增多，成分的復(fù)雜性將引起光譜與濃度間的非線性變動(dòng)關(guān)系，PLS方法所建模型的預(yù)測(cè)準(zhǔn)確性將下降。為得到準(zhǔn)確性更高的預(yù)測(cè)結(jié)果，需要分別建立單種原料油的PLS校正模型，或采用更為復(fù)雜的算法，如局部權(quán)重回歸、神經(jīng)網(wǎng)絡(luò)或支持向量機(jī)等建立非線性校正模型。

　　3.3.3模型適用性的研究

　　以花椒油生物柴油為例，研究了多種原料油生物柴油的NIR校正模型對(duì)新原料油生物柴油的適用性。若校正集中不含待測(cè)類型的樣本，其預(yù)測(cè)結(jié)果將會(huì)存在一定的系統(tǒng)誤差，且模型的適應(yīng)性判據(jù)也會(huì)給出提示。為考察向校正集添加少量樣本對(duì)擴(kuò)充模型適應(yīng)性的效果，從47個(gè)花椒油生物柴油中任意選取10個(gè)樣品，用于模型的擴(kuò)充，剩余37個(gè)作為驗(yàn)證集樣本。

　　將花椒油生物柴油樣本從多種原料油校正集中剔除，建立剩余5種原料油生物柴油的校正模型。然后，再將任意選取的10個(gè)花椒油生物柴油樣本添加到校正集中，建立含有10個(gè)花椒油生物柴油的多原料油生物柴油的校正模型。用以上建立的兩類模型分別對(duì)37個(gè)驗(yàn)證集樣本進(jìn)行預(yù)測(cè)分析(表4)。由表4可見(jiàn)，未含花椒油生物柴油的模型仍具有較好的預(yù)測(cè)趨勢(shì)，但存在一定的系統(tǒng)誤差，預(yù)測(cè)準(zhǔn)確性也有明顯下降。當(dāng)模型中添加了10個(gè)花椒油生物柴油樣本后，消除了系統(tǒng)誤差，且預(yù)測(cè)準(zhǔn)確性有了顯著改善，與3.3.2節(jié)建立的模型基本相當(dāng)。

　　上述結(jié)果說(shuō)明，對(duì)于模型中未含的新類型原料油生物柴油樣本，只需向多種原料油校正集中添加10個(gè)以上樣本便可較為準(zhǔn)確地測(cè)定新類型生物柴油的成分。表4不同校正模型對(duì)37個(gè)花椒油生物柴油的預(yù)測(cè)結(jié)果(略)

　　3.4近紅外分析方法的重復(fù)性

　　由于光譜測(cè)量的穩(wěn)定性，NIR分析方法具有良好的重復(fù)性。隨機(jī)選取一個(gè)樣本重復(fù)掃描7次的NIR光譜，并分別調(diào)用建立的多種原料油的生物柴油模型，其質(zhì)量分?jǐn)?shù)的相對(duì)標(biāo)準(zhǔn)偏差分別為：甘油0.74%，甲酯0.26%，單甘酯2.86%，二甘酯0.70%，三甘酯1.41%。說(shuō)明采用NIR分析方法測(cè)定混合油生物柴油的成分具有較高的精密度。

　　3.5小結(jié)

　　通過(guò)以上研究，可以得到以下結(jié)論：(1)通過(guò)偏最小二乘方法可以建立適合多種原料油生物柴油的通用近紅外校正模型，其結(jié)果完全滿足過(guò)程控制分析的要求。但若想得到更準(zhǔn)確的預(yù)測(cè)結(jié)果，需要建立單種原料油的PLS校正模型，或者采用更為復(fù)雜的校正算法;(2)對(duì)于新型生物柴油，只需向校正集添加少量樣本，擴(kuò)充模型覆蓋范圍后，便可較為準(zhǔn)確地對(duì)該類樣本進(jìn)行預(yù)測(cè)分析。這樣，基于本研究建立的這些基礎(chǔ)校正模型，針對(duì)新類型原料油或加工工藝，只需添加少量樣本，便可直接應(yīng)用于生物柴油的實(shí)際工業(yè)生產(chǎn)中;(3)近紅外光譜方法有望較好解決生物柴油傳統(tǒng)分析方法的諸多弊端，可以快速準(zhǔn)確地測(cè)定生物柴油的主要成分含量，適合于生物柴油生產(chǎn)過(guò)程的中間控制分析。