唐保山,李 坤, 張 弘, 張?chǎng)?關(guān)慶芳, 史正軍
(1.西南林業(yè)大學(xué) 林學(xué)院,云南 昆明 650224; 2.中國(guó)林業(yè)科學(xué)研究院資源昆蟲(chóng)研究所國(guó)家林業(yè)和草原局特色森林資源工程技術(shù)研究中心,云南 昆明 650233;3.安寧戴科精細(xì)化工有限公司,云南 昆明 650301)
紫膠又名蟲(chóng)膠、赤膠、紫草茸,是紫膠蟲(chóng)吸食寄主植物汁液后分泌的一種純天然樹(shù)脂混合物[1-4]。紫膠是我國(guó)重要的工業(yè)原料,因其具有優(yōu)良的抗酸、防銹、耐油、電絕緣性和熱塑性能,廣泛用于國(guó)防、涂料、機(jī)械、醫(yī)藥、印刷和食品行業(yè)[5-7]。依據(jù)紫膠的不同加工階段和產(chǎn)品形態(tài),紫膠又可分為紫膠原膠、顆粒紫膠、片膠、漂白紫膠等。紫膠原膠是直接從寄主樹(shù)上采收下來(lái)的成熟膠梗,也是紫膠其他產(chǎn)品加工的原料,主要由紫膠樹(shù)脂、蠟質(zhì)和色素組成[8-10]。長(zhǎng)期以來(lái),紫膠原膠的質(zhì)量評(píng)價(jià)多依據(jù)國(guó)標(biāo)法,盡管?chē)?guó)標(biāo)法具有可靠性、普適性強(qiáng)的特點(diǎn),但操作復(fù)雜、對(duì)人員和熟練程度要求高、檢測(cè)耗時(shí)等缺點(diǎn)已在很大程度上妨礙了國(guó)標(biāo)法在行業(yè),尤其在流通領(lǐng)域及企業(yè)中的普及和應(yīng)用。在工業(yè)快速發(fā)展的今天,如何提高紫膠檢測(cè)方法的準(zhǔn)確性、時(shí)效性和便攜性已成為提高工作效率的重要環(huán)節(jié)[11-12]。近紅外光譜法是近年發(fā)展起來(lái)的一種非破壞性分析技術(shù),具有檢測(cè)成本低、速度快、樣品預(yù)處理簡(jiǎn)單等優(yōu)點(diǎn)[13-15]。近紅外光譜分析法是一種間接分析技術(shù),是用統(tǒng)計(jì)學(xué)方法在樣品待測(cè)屬性值與近紅外光譜數(shù)據(jù)之間建立一個(gè)關(guān)聯(lián)模型,之后通過(guò)驗(yàn)證完成模型的建立,利用近紅外光譜儀測(cè)得的樣品光譜數(shù)據(jù)即可通過(guò)軟件模型預(yù)測(cè)被測(cè)組分含量[16-17]。目前,近紅外光譜法已在藥品、紡織物、酒類(lèi)、飼料、糧油、煙草等領(lǐng)域得到廣泛的應(yīng)用[18-19]。其中飼料中水分、粗蛋白質(zhì)、粗纖維等含量的測(cè)定[20],玉米中淀粉含量的測(cè)定[21]已成為國(guó)標(biāo)中采用的標(biāo)準(zhǔn)檢測(cè)方法。本研究依據(jù)化學(xué)法測(cè)定紫膠原膠中樹(shù)脂、蠟質(zhì)和色素含量的真實(shí)值,探索基于近紅外光譜分析模型的構(gòu)建方法,以期將近紅外光譜法用于紫膠的檢測(cè),實(shí)現(xiàn)快速定量檢測(cè)紫膠原膠中主組分含量的目的,為建立科學(xué)化、自動(dòng)化的紫膠現(xiàn)代質(zhì)量標(biāo)準(zhǔn)評(píng)價(jià)體系提供理論依據(jù)。
紫膠原膠樣本(中國(guó)21、泰國(guó)2、印度2、緬甸2和越南2),由墨江縣洪森蟲(chóng)膠有限公司、墨江滇力林化廠、墨江森源科技有限公司、綠春興龍紫膠有限公司、云縣凝鑫紫膠產(chǎn)業(yè)發(fā)展有限公司和安寧戴科精細(xì)化工有限公司提供。為保證測(cè)試樣品的均勻性,把收集到的29個(gè)紫膠原膠樣本分別用高速粉碎機(jī)粉碎,取粒徑≤0.38 mm的樣品備用。無(wú)水碳酸鈉、鹽酸、鄰苯二甲酸氫鉀、四氯化碳、95%乙醇、無(wú)水乙醇和氫氧化鉀均為分析純。
MPA傅里葉變換近紅外光譜儀,德國(guó)Brucker公司;Agilent Cary Seties紫外(UV-vis)可見(jiàn)分光光度計(jì),Agilent Technologies公司;HH- 601水浴磁力攪拌器;Testo-206便攜式pH計(jì);SHZ-D(Ⅲ)循環(huán)水式多用真空泵。
1.2.1紫膠樹(shù)脂 紫膠樹(shù)脂的測(cè)定參照國(guó)標(biāo)GB/T8143—2008[22]的方法進(jìn)行。精確稱(chēng)取樣品5 g放入用95%乙醇萃取過(guò)的定性濾紙中,用棉線捆緊后放入事先在100 ℃質(zhì)量恒定的萃取瓶中,置入索氏提取器,用95%乙醇在80 ℃提取4 h,將萃取瓶中的溶劑減壓蒸發(fā),在100 ℃下干燥至質(zhì)量恒定,最后扣除蠟質(zhì)質(zhì)量,即得紫膠樹(shù)脂質(zhì)量,平行測(cè)定3次。樹(shù)脂質(zhì)量分?jǐn)?shù)按式(1)計(jì)算。
1.2.2紫膠蠟 紫膠蠟的測(cè)定參照國(guó)標(biāo)GB/T8143—2008[22]中的方法進(jìn)行。精確稱(chēng)取樣品10 g放入燒杯中,之后加入150 mL溶有2.5 g碳酸鈉的熱水,放在沸水浴中加熱攪拌。待試樣溶解后再加熱2~3 h(不要攪拌),然后將燒杯從水浴中取出,冷卻至室溫,蠟質(zhì)結(jié)晶析出,用四氯化碳萃取過(guò)的濾紙過(guò)濾,用水洗滌至濾液無(wú)色,將濾紙取出放在60 ℃的烘箱中烘去水分,再用一張濾紙包好,用棉線捆緊,放入事先在100 ℃質(zhì)量恒定的萃取瓶中,置入索氏提取器,用四氯化碳在80 ℃提取4 h,將萃取瓶中的四氯化碳減壓蒸發(fā),放入干燥箱中,100 ℃下干燥至質(zhì)量恒定,平行測(cè)定3次。蠟質(zhì)量分?jǐn)?shù)按式(2)計(jì)算。
ω1=(m2-m1-m0)/m×100%
(1)
ω2=(m3-m1)/m×100%
(2)
式中:ω1—紫膠中樹(shù)脂質(zhì)量分?jǐn)?shù),%;m—樣品質(zhì)量,g;m0—蠟質(zhì)質(zhì)量,g;m1—萃取瓶質(zhì)量,g;m2—萃取瓶與樹(shù)脂質(zhì)量和,g;ω2—紫膠中蠟的質(zhì)量分?jǐn)?shù),%;m3—萃取瓶和蠟質(zhì)質(zhì)量和,g。
1.2.3紫膠色素 參照文獻(xiàn)[23]測(cè)定水溶物。精確稱(chēng)取樣品10 g放入250 mL燒杯中,加入150 mL去離子水浸泡過(guò)夜,之后抽濾,用水洗至無(wú)色,最后濾液倒入燒杯中在80 ℃鼓風(fēng)干燥箱烘至質(zhì)量恒定,平行測(cè)定3次。水溶物質(zhì)量=燒杯與水溶物質(zhì)量之和-燒杯質(zhì)量。
參照國(guó)標(biāo)GB1886.17—2015[24]方法進(jìn)行色價(jià)測(cè)定。稱(chēng)取干燥水溶物0.1 g置于50 mL燒杯中,加入10 g/L碳酸鈉溶液10 mL攪勻、溶解,傾入100 mL容量瓶中,用少量水洗滌燒杯,洗滌液并入100 mL容量瓶中,用水稀釋至刻度、搖勻。準(zhǔn)確吸取5 mL置于100 mL容量瓶中,用0.1 mol/L鹽酸溶液調(diào)pH值至3.0左右,用pH值3.0緩沖溶液稀釋至刻度、搖勻,即為試樣溶液。取出試樣溶液置于1 cm比色皿中,于紫外分光光度計(jì)為490 nm波長(zhǎng)處測(cè)量吸光度,平行測(cè)定3次。色價(jià)按式(3)計(jì)算。
參照文獻(xiàn)[25]方法測(cè)定色素含量。已知色素標(biāo)準(zhǔn)品的含量,測(cè)出色素標(biāo)準(zhǔn)品的色價(jià),樣品通過(guò)測(cè)定其色價(jià),按照色素標(biāo)準(zhǔn)品折合試樣色素含量,最后計(jì)算出紫膠原膠中的色素,平行測(cè)定3次。色素按式(4)計(jì)算。
(3)
ω3=m4/m×ω4
(4)
1.3.2模型建立的方法和步驟 一個(gè)完整的近紅外預(yù)測(cè)模型包括6個(gè)部分,具體包括:樣本數(shù)據(jù)采集及樣品集劃分、近紅外光譜測(cè)定條件設(shè)定及光譜圖采集、確定數(shù)據(jù)樣本預(yù)處理方法、樣品近紅外光譜區(qū)域的劃分及最優(yōu)波段的選取、主成分?jǐn)?shù)的選擇和模型質(zhì)量評(píng)價(jià)。
1.3.2.1樣本數(shù)據(jù)采集及樣品集劃分 所收集的樣本見(jiàn)表1,囊括了中國(guó)、越南、泰國(guó)、緬甸和印度等紫膠主產(chǎn)國(guó),為基于近紅外法的紫膠原膠主組分快速預(yù)測(cè)模型構(gòu)建提供了數(shù)據(jù)支撐。
表1 不同來(lái)源紫膠樣本中樹(shù)脂、蠟質(zhì)及色素含量1)
紫膠樣本按4.5 ∶1隨機(jī)分選為校正集和驗(yàn)證集,通過(guò)表2可知,校正集和驗(yàn)證集中樹(shù)脂、蠟質(zhì)和色素的含量范圍、平均值和極差表明,所選樣品中各組分含量范圍較寬,代表性較強(qiáng),且校正集樣品范圍包含驗(yàn)證集所有樣品的組成范圍,符合建立模型樣品選擇的原則。
表2 紫膠原膠中樹(shù)脂、蠟質(zhì)和色素樣品集的劃分
1.3.2.2測(cè)定條件設(shè)定及光譜圖采集 為保證儀器的穩(wěn)定性,開(kāi)機(jī)后預(yù)熱30 min;設(shè)置掃描時(shí)儀器分辨率為8 cm-1,掃描次數(shù)為64次,光譜范圍為波數(shù)12500~4000 cm-1;以空IN312-SH型樣品杯(Φ=97 mm)作為參比樣品進(jìn)行光譜采集,上樣量以樣品體積的2/3為宜,為確保樣品采集的均勻性,采用旋轉(zhuǎn)模式,每個(gè)紫膠原膠樣品重復(fù)采集3次[27],取其平均光譜作為樣品代表性光譜。
1.3.2.3確定數(shù)據(jù)樣本預(yù)處理方法 近紅外光譜易受到樣品的均勻性、噪音信號(hào)以及雜散光等影響,為保證譜圖基線平穩(wěn)和所建立校正模型的準(zhǔn)確性,必須對(duì)原始光譜進(jìn)行預(yù)處理。常用的預(yù)處理方法有無(wú)光譜處理、多元散射校正(MSC) 、消除常量偏移量(ECO)、矢量歸一化(VN)、最小-最大歸一化(M-MN)、一階導(dǎo)數(shù)+矢量歸一化 (FD+VN)[28]。
1.3.2.4區(qū)域劃分及最優(yōu)波段選取 偏最小二乘法(PLS) 雖然可以處理全譜信息,但全譜中往往包含了過(guò)多的冗余信息,這既影響了建模的準(zhǔn)確性,過(guò)大的數(shù)據(jù)量也降低了模型預(yù)測(cè)的計(jì)算速率[29]。因此,需要對(duì)光譜波段進(jìn)行提取,改善模型性能。依據(jù)紫膠原膠主組分近紅外譜圖的吸收特征,將原始光譜分為5個(gè)波段范圍進(jìn)行交互建模,并考察各個(gè)波段的模型與組分含量間的相關(guān)性,探索最優(yōu)波段,5個(gè) 波段分別為9400~7500、7500~6100、6100~5450、5450~4600和4600~4250 cm-1。
1.3.2.5主成分?jǐn)?shù)的選擇 在建模過(guò)程中,采用不同的主成分?jǐn)?shù),模型的預(yù)測(cè)能力也會(huì)有較大的差異。因此合理確定參加建模的主成分?jǐn)?shù)是充分利用光譜信息和濾除噪音的有效方法之一。通常數(shù)值過(guò)大,所建模型包含太多的噪音,會(huì)出現(xiàn)過(guò)擬合現(xiàn)象;主成分太小,則導(dǎo)致建模信息不全,模型預(yù)測(cè)能力相對(duì)較差。當(dāng)RMSECV最小時(shí),對(duì)應(yīng)主成分?jǐn)?shù)最佳。
(5)
(6)
(7)
(8)
利用近紅外儀器自帶的OPUS 7.5分析軟件采集的紫膠樣品的近紅外光譜,原始光譜如圖1(a)所示。從圖1(a)可以看出,紫膠主組分在峰位、峰強(qiáng)和峰形等方面重疊嚴(yán)重,很難直觀分析,就需要對(duì)近紅外光譜信號(hào)進(jìn)行分析處理,不同的預(yù)處理方法得到的譜圖見(jiàn)圖1(b)、1(c)和1(d)所示。
通過(guò)對(duì)不同波段的篩選,由表3可知,樹(shù)脂、蠟質(zhì)和色素的最優(yōu)建模波段分別為9403.7~5446.3、4601.6~4246.7 cm-1,9403.7~7498.3、6102~5446.7 cm-1和7502.1~6098.1、5450.1~4246.7cm-1,說(shuō)明在該波數(shù)范圍內(nèi),組分與模型的相關(guān)性最佳。
為更好的驗(yàn)證上述最優(yōu)波段選擇的可靠性,選擇了脫蠟脫色片膠、紫膠蠟和紫膠紅色素作為紫膠樹(shù)脂、蠟質(zhì)和色素的純凈物分別進(jìn)行全譜掃描。
a.紫膠原始圖original spectra of sticklac; b.矢量歸一化處理vector normalization (VN);c.最小-最大歸一化處理
表3 樹(shù)脂、蠟質(zhì)和色素最佳預(yù)處理方法和最優(yōu)波段的選擇
續(xù)表3
a.樹(shù)脂resin; b.蠟質(zhì)wax; c.色素pigment圖2 不同樣品的近紅外光譜圖Fig.2 Near infrared spectra of different samples
由圖2可知,純樹(shù)脂的近紅外光譜圖與紫膠原膠的光譜圖幾乎一致,是因?yàn)樽夏z原膠中樹(shù)脂是主要組分,樹(shù)脂吸收強(qiáng)度最強(qiáng),在譜圖疊加中占據(jù)主導(dǎo),故從純樹(shù)脂的光譜圖上很難找到最優(yōu)波段。由表3得到蠟質(zhì)的最優(yōu)建模波數(shù)為9403.7~7498.3 cm-1、6102.0~5446.7 cm-1,對(duì)比圖2蠟質(zhì)圖譜可知,此波數(shù)范圍正好在蠟質(zhì)特征吸收范圍內(nèi)。由色素近紅外光譜可知,色素主要吸收波段為7378.0~6258.0 cm-1和5411.0~4511.0 cm-1,而通過(guò)不同光譜段處理得到最優(yōu)波段為7502.1~6098.1、5450.1~4246.7 cm-1,與色素主要吸收波段吻合。
樹(shù)脂、蠟質(zhì)和色素的主成分?jǐn)?shù)與RMSECV的關(guān)系如圖3所示,主成分?jǐn)?shù)的大小影響模型的質(zhì)量,一般認(rèn)為RMSECV最小時(shí),對(duì)應(yīng)主成分?jǐn)?shù)最優(yōu)。由圖3可知,樹(shù)脂、蠟質(zhì)和色素主成分?jǐn)?shù)分別為8、10、6時(shí),對(duì)應(yīng)的RMSECV最小,說(shuō)明在該主成分?jǐn)?shù)下,該模型的誤差最低。
a.樹(shù)脂resin; b.蠟質(zhì)wax; c.色素pigment
利用以上建模條件,對(duì)樹(shù)脂、蠟質(zhì)和色素校正集進(jìn)行內(nèi)部交叉驗(yàn)證,建立校正模型的樣品數(shù)為24個(gè),分別剔除0、0、2個(gè)光譜異常的樣品,實(shí)際采用的樣品只有24、24、22個(gè),利用建立的樹(shù)脂、蠟質(zhì)和色素定量模型預(yù)測(cè)驗(yàn)證集為5個(gè)樣品。樹(shù)脂、蠟質(zhì)和色素的校正集和驗(yàn)證集的預(yù)測(cè)值和真實(shí)值的關(guān)系如圖4所示。
△校正集calibration sets; ○驗(yàn)證集validation sets
所建模型的各項(xiàng)數(shù)理統(tǒng)計(jì)指標(biāo)都達(dá)到預(yù)期值,所有樣品的近紅外光譜測(cè)定結(jié)果和標(biāo)準(zhǔn)方法之間的相關(guān)性以及誤差基本都在要求的再現(xiàn)性范圍以?xún)?nèi),近紅外光譜所建紫膠原膠中3種主組分模型預(yù)測(cè)樣品的預(yù)測(cè)值與國(guó)標(biāo)方法測(cè)得的化學(xué)值基本一致,說(shuō)明所建模型的預(yù)測(cè)效果良好。
由表4可知,對(duì)樹(shù)脂、蠟質(zhì)和色素的預(yù)測(cè)值和真實(shí)值做了比較,每個(gè)樣品的5次預(yù)測(cè)結(jié)果平均差異都小于0.10%;并且對(duì)每個(gè)樣品的5次預(yù)測(cè)值做了顯著性分析,結(jié)果P>0.05,均無(wú)顯著性差異,可見(jiàn)所建模型精密度和穩(wěn)定性良好。
表4 模型預(yù)測(cè)樣品的精密度和穩(wěn)定性檢驗(yàn)(n=5)