树龄约50年的落叶松[Larix gmelinii (Rupr.) Kuzen.]气干材板材锯制成规格为70 mm(纤维方向)× 30 mm(弦向)× 5 mm(径向)的片状试样。 选择无树节、 早晚材分布均匀的试样共105件, 平均分为五组, 其中一组为对照, 其余四组试材在弦切面进行紫外光照射人工老化试验。

紫外老化测试采用美国Q-Lab公司生产的型号为QUV/SPRAY的人工加速紫外老化箱以模拟户外自然老化进程, 参照标准《ASTG G154-06》进行老化试验。 每个老化周期由8 h紫外辐射和4 h冷凝过程。 8 h的紫外辐射波长340 nm, 照射强度为0.89 W· m^-2, 老化箱黑板温度为60 ℃; 冷凝阶段紫外灯自动关闭, 老化箱内部温度为50 ℃。 本实验中分别在老化时间180, 540, 900和1 080 h后取出试样用于后续的色度学等检测与分析。

通过色差仪(CM-2300d, Konica Minolta, INC, Japan)测量光老化木材各项色度学参数。 每个试样上取9个点测量色度学参数, 结果取平均值。 色度学参数Δ L^* , Δ a^* , Δ b^* 参照CIE1976L^* a^* b^* 系统确定, 总色差(Δ E^* )表示光老化试材与对照组相应差值, 计算公式如式(1)

ΔE* =(ΔL* 2+Δa* 2+Δb* 2)1/2(1)

式(1)中: Δ L^* , Δ a^* 和Δ b^* 分别表示试材老化前后的明度值差、 红绿值差和黄蓝值差。

采用MPA型近红外光谱仪(德国 Bruker)采集近红外光谱, 光谱采用直径为20 mm的积分球漫反射方式, 波谱范围8 000~4 000 cm^-1, 分辨率为8 cm^-1, 每个扫描点面积约7 mm², 采用漫反射模式, 每个位点自动扫描次数为32次, 获得一张图谱。 落叶松试材弦切面上采集9个点扫描获得光谱信息, 结果取平均值。

将105块试材采集到的色度学参数(L^* , a^* , b^* )与其对映的近红外谱图联立, 参数模型的样品数量比为校正模型:验证模型=2:1。 使用OPUS 7.6 (Bruker Optics, Germany)定量分析模块基于PLS建立落叶松颜色预测模型, 使用留一法交叉验证(leave-one-out cross validation)检验模型。 模型建立采用多种预处理方法: 矢量归一、 多元散射校正、 一阶导数。 选取决定系数(R²)、 校正均方根误差(RMSEC)、 交叉验证均方根误差(RMSECV)及相对分析误差(RPD)等指标对模型进行评价。

木材经180 h紫外光照射后, 颜色从黄褐色变成棕褐色, 继续经900和1 080 h照射后, 最终变为深褐色、 黑褐色。 不同老化时间后试材颜色变化的色度学指数(L^* , a^* , b^* )变化见表1。

表1

Table 1

表1(Table 1)

表1 未处理和光老化落叶松试材的色度学参数 Table 1 Chromaticity parameters of larch wood with and without ultraviolet light irradiation 紫外老化时间 L* a* b* Δ E* 未处理 72.01 11.29 23.74 - 180 h 53.61 19.57 33.96 22.62 540 h 49.04 21.07 33.97 26.98 900 h 42.78 20.88 27.73 31.02 1 080 h 39.19 19.66 23.00 33.88

表1 未处理和光老化落叶松试材的色度学参数 Table 1 Chromaticity parameters of larch wood with and without ultraviolet light irradiation

由表1可以看出, 随着照射时间的延长, 明度值L^* 持续下降, 与对照材的72.01相比, 经1 080 h紫外光照射后, L^* 值降至39.19, 表明落叶松表面的颜色加深约一倍。 与未处理材相比较, L^* 值的变化分为两个阶段: (1)紫外照射时间180 h时, 变化较明显, 180 h→ 540 h变化不大; (2)照射时间达到1 080 h时, 再次发生明显的变化。 紫外光照射使木材表面a^* 与b^* 明显增加, 与对照材相比, 木材表面变红、 变黄, 照射540 h时达到最大值, 其后红色、 黄色略微变浅。 此外, 色差值Δ E^* 与照射时间呈正相关, 随着照射时间的延长, a^* 和b^* 增加的趋势有所减弱。 上述结果表明, 紫外光照射下, 木材材色变红、 变黄、 材色加深。

图1为谱带7 100~5 450 cm^-1 (a)和谱带5 000~4 000 cm^-1 (b)范围内未处理、 紫外光照射至1 080 h试材的NIRs二阶导数谱。 在NIRs区, 7 200~6 000 cm^-1范围主要表征纤维素上羟基伸缩振动一次倍频, 6 000~5 400 cm^-1范围主要表征脂肪族和芳香族C—H基团伸缩振动一次倍频与O—H基团的合频, 而5 000~4 000 cm^-1谱带范围主要表征羟基伸缩振动与变形振动的合频、 以及与次甲基各种振动模式的合频^[8]。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 未处理和光老化落叶松试材二阶导数谱图 (a): 7 100~5 450 cm-1; (b): 5 000~4 000 cm-1Fig.1 The NIRs second derivative spectra of larch wood with and without ultraviolet light irradiation (a): 7 100~5 450 cm-1; (b): 5 000~4 000 cm-1

随着照射时间的延长, 照射时间180和540 h的试材两者吸光度之间差别变化不大, > 900 h时吸光度则发生明显的变化, 该结果与2.1中色度学的测量结果相一致。 图1(a)中, 从高波数往低波数方向, 观察到5个明显的吸收峰, 图1(b)中, 观察到4个明显的吸收峰(吸收峰归属如表2所示)。 在图1(a, b)中, 首先, 除d外, 所有的吸收峰随照射时间延长, 相应基团相对数量的增加。 其次, 谱带d为木质素芳香环骨架上C—H基团伸缩振动的一次倍频, 该吸收峰随着紫外光老化时间的延长而降低, 表明木质素相对含量减少, 木质素发生降解。 吸收带a, b, c, f和i归属于纤维素, 随光老化时间的延长, 吸光度增加, 表明试材中纤维素含量均呈上升趋势。 其中b, c和f谱带中由于结晶区和半结晶区O—H的伸缩振动, 纤维素含量上升的趋势相对不明显。

表2

Table 2

表2(Table 2)

表2 落叶松典型吸收带归属 Table 2 Assignment of representative absorption bands in Larch wood No. 波数/cm-1 吸收带归属 文献 a 6 996 1st OT. O—H str./纤维素非结晶区 [8] b 6 773 1st OT. O—H str./纤维素半结晶区 [9] c 6 287 1st OT. O—H str./纤维素结晶区 [8] d 5 986 1st OT. C—H str.木质素芳香环骨架 [9] e 5 800 1st OT. C—H str./半纤维素的呋喃和吡喃 [8] f 4 852~4667 O—H str.+C—H def./纤维素结晶区或半结晶区 [9] g 4 404 O—H str.+C—H2 def./纤维素和半纤维素 [8] h 4 281 C—H str.+C—H def.纤维素和半纤维素、 木聚糖 [8] i 4019 C—H str.+C-C str./纤维素 [8] 注: OT.(overtone): 倍频; str.(stretching vibration): 伸缩振动; +: 组合频; def.(deformation vibration): 变形振动

表2 落叶松典型吸收带归属 Table 2 Assignment of representative absorption bands in Larch wood

通过对照材二阶导数谱带减去1 080 h试材二阶导数谱带, 获得NIRs二阶导数差谱, 见图2。 差谱体现木材化学组分相对含量的增加或减少, 正值表示相对含量的增加, 负值表示相对含量的减少。 图2(a, b)分别表示光谱范围7 100~5 450和5 000~4 000 cm^-1的差谱(虚线)和二阶导数光谱(实线), 字母编号见表2。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 未处理和紫外光照射1 080 h试材差谱(虚线)及其近红外二阶导图(实线) (a): 7 100~5 450 cm-1; (b): 5 000~4 000 cm-1Fig.2 The NIRs difference spectrum (dotted line) and Second-derivative spectra (solid lines) of without and with ultraviolet light irradiation of 1 080 h (a): 7 100~5 450 cm-1; (b): 5 000~4 000 cm-1

图2(a)中, 二阶导数谱带出现的5个吸收峰谱带, 分别位于6 996, 6 773, 6 287, 5 986和5 800 cm^-1处。 其中, 6 996 cm^-1为纤维素的非结晶区自由羟基以及弱氢键结合的羟基伸缩振动一次倍频, 谱带6 773和6 287 cm^-1分别表示纤维素半结晶区和结晶区羟基伸缩振动的一次倍频, 而5 800 cm^-1则表示半纤维素中呋喃环上C—H伸缩振动的一次倍频。 上述表示木材纤维素以及半纤维素的各个吸收谱带特征峰差谱均为正值, 表明紫外光辐射使木材表面纤维素和半纤维素相对含量增加。 谱带5 986 cm^-1处为木质素芳香环骨架上C—H伸缩振动一次倍频, 该处差谱值为负, 表明木质素相对含量降低。 图2(b)中, 二阶导数谱带出现4个吸收谱带, 分别位于4 852~4 667, 4 404, 4 281和4 019 cm^-1处, 均归属于纤维素或半纤维素中相应基团的伸缩振动或者变形振动。 由差谱图(虚线)可知, 以上4个吸收谱带均为正值, 表明经紫外光照射后, 落叶松表面纤维素和半纤维素相对含量均有增加。 上述NIRs分析结果表明, 紫外光照射木材表面木质素发生降解, 而纤维素和半纤维素的相对含量则增加。 有研究表明, 木材表面经紫外光老化后, 以木质素为主, 纤维素和半纤维素将出现不同程度的降解^[10]。 本研究中, 通过差谱观察到木质素降解、 相对含量减少外, 纤维素和半纤维等多糖类组分的解聚没有体现在NIRs所反映的信息中, 考虑该结果归结于木质素降解的速度较快, 而多糖类组分发生不同程度的解聚时, 由于速度较慢, 则随着木质素的降解, 多糖类组分相对含量增加。

PLS结合适宜预处理方法能减少基线漂移、 降噪等影响, 使建立的模型更稳健。 根据1.5节建模及评价方法, 选择低的RMSECV, 低维数, 高R²参数作为预测模型评价参数。 表3列出了落叶松材色预测模型的各项参数, 参考Zornoza等的研究总结^[11]。

表3

Table 3

表3(Table 3)

表3 落叶松试材颜色参数预测模型的各项参数 Table 3 Detailed parameters of color prediction models of larch wood 参数 预处理方法 光谱范围/cm-1 维数 RC2 RV2 RMSEC RMSECV RPDC RPDV L* 矢量归一+ 一阶导数+ 多元散射校正 6 102~4 598 5 450~4 598 4 602~4 247 10 0.970 8 0.948 8 2.14 2.68 5.85 4.42 a* 一阶导数 5 450~4 598 6 102~5 446 9 0.954 2 0.928 0 0.878 1.04 4.67 3.73 b* 一阶导数 5 450~4 247 6 102~5 446 9 0.9448 0.831 0 1.31 2.18 4.26 2.43 注: C: 校正模型; V: 交叉验证模型

表3 落叶松试材颜色参数预测模型的各项参数 Table 3 Detailed parameters of color prediction models of larch wood

在本研究中, 色度学参数L^* , a^* , b^* 在适宜的光谱范围及预处理方法内所建预测模型的 RV2和RPD_V分别为0.948 8和4.42; 0.928 0, 3.73和0.831 0, 2.43。 L^* 和a^* 预测模型均为“ 优秀” ; b^* 预测模型介于“ 好” 及“ 还可以” 之间。

对颜色参数测量-预测模型进行拟合得到线性关系图(图3)。 图中L^* 和a^* 预测值与测量值散点靠近斜率为45° 直线, 表明所构建的模型具有准确回归结果和良好的预测精度, b^* 相对分散。 本研究对老化落叶松表面材色建立的NIRs定量预测模型, 满足木材颜色测定的要求。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 色度学参数a* , b* , L* (a, b, c)预测值和测量值相关关系图Fig.3 Correlation between colorimetric parameters values a* , b* , L* (a, b, c) measured value and predicted value