目前，常用的蛋白质检测方法有凯氏定氮法、紫外吸收法、电泳法、低场核磁共振法和近红外光谱法等。近红外光谱技术属于间接分析方法，需要通过多元校正方法建立光谱信息与待测成分间的关联，由于仪器性能、外界环境、自身特性等因素影响，使得光谱数据中往往含有共线性、噪声及外界干扰信息等，导致建模结果误差较大甚至无法使用，故此需要优化条件。
由于牛乳介质中存在颗粒不均匀性，使其对近红外光谱有很强的散射作用，多采用正交信号校正、多元散射校正、求导等对光谱进行处理。在已有的建模条件研究中往往只考虑某一方面的优化，而实际问题中各因素完全独立的情况极为少见。因此，河南工业大学粮油食品学院的彭丹、刘亚丽和李林青等人针对各因素间存在相互作用影响的情况，在单一建模条件研究的基础上，运用正交试验设计对建模条件进行优化，建立一个准确度高、重复性好、适用于牛乳蛋白质含量分析的优化检测体系。
牛乳近红外光谱及二维相关分析
如图1所示，脂肪、蛋白质和乳糖等成分近红外吸收分布于整个光谱区域，且相互重叠。在1 460、1 974 nm波长附近有2 个主要吸收峰，分别为牛乳中水、蛋白质、脂肪和乳糖等成分的结构和组成信息。为确定蛋白质的主要吸收位置，以蛋白质含量为外扰，对牛乳近红外光谱进行同步二维相关分析，结果如图2所示。
图2 a中的峰有自相关峰和交叉峰两类，位于对角线上的峰为自相关峰，位于对角线外为交叉峰。自相关谱是由对角线上自相关峰构成的谱图（图2b）。由图2可以看出，在波长978、1 164、1 420、1 524、1 659、1 860 nm和2 238 nm处存在较强的自相关峰，其中2 238 nm处的自相关峰与蛋白质有关；在主对角线以外，在（1 420 nm，1 860 nm）、（1 860 nm，2 238 nm）、（1 420 nm，2 238 nm）、（1 524 nm，2 238 nm）位置处存在明显正交叉峰，表明波长1 420、1 524、1 860 nm和2 238 nm处的吸收峰来源相同，均是由蛋白肽键和氨基中NüH键对光谱吸收形成。可见，1 400～2 300 nm波段的光谱随蛋白质含量变化极显着。短波区域（780～1 100 nm）的光透射性强，但吸收系数较小，与长波相比蛋白质含量信息量相对较少（图2b），增加了数据分析和建模的难度，结合二维相关分析，本研究选择波段1 400～2 300 nm作为蛋白质检测的研究区域。
单因素试验结果
2.1 预处理方法的对比
为降低干扰信息的影响和提升模型预测精度，分别采用矢量归一化、多元散射校正（MSC）、一阶导数（1st）、二阶导数（2nd）、标准正态变量变换（SNV）方法对原始光谱进行预处理后建立PLS模型，结果显示，除2nd方法外，矢量归一化、MSC、1st、SNV方法处理后模型的相关系数和预测准确性均有明显提高，表明这4 种方法能够扣除光谱中与待测成分无关的干扰信号，保留有用信息。其中，经SNV法处理后的蛋白质模型的预测效果最佳，其次为MSC法。与预处理前相比，经SNV法和MSC法处理后模型的RMSEP分别降低了31.1%和27.8%。
通过粒径分析可知，均质后牛乳中脂肪球的粒径分布于0.2～1.0 μm之间，存在明显脂肪球分布的不均匀性，这使其对近红外光有较强的散射作用，由于MSC法和SNV法能够消除颗粒大小及均匀性变化对光谱的影响，所以经MSC法和SNV法处理后可有效消除散射作用引起的干扰。相反，2nd法在去除背景干扰、提高灵敏度的同时放大了噪声，降低了信噪比，使得模型的预测能力下降，导致RMSEP升高了91.9%。可见，预处理方法的选择对蛋白质检测结果有较大影响，其中MSC和SNV是较为理想的预处理方法。
2.2 波段选择的影响
本实验在1 400～2 300 nm波长范围内研究不同波段对蛋白质含量检测的影响。如图3所示，3 个波段1 400～2 300、1 400～1 800 nm和1 800～2 300 nm的RMSEP分别为0.146、0.142和0.125，其中1 800～2 300 nm波段的预测结果明显好于其他波段。虽然1 400～1 800 nm波段也含有蛋白质相关的光谱信息，但是由于存在脂肪较强的散射作用和水的强吸收作用，加之牛乳中水的含量（＞87%）远大于蛋白质，使得此波段存在复杂的背景干扰。因此，1 800～2 300 nm是蛋白质理想的检测波段。
2.3 建模方法比较
采用MLR、PCR、PLS和SVR对同一牛乳样本进行定量分析，结果如图4所示。4 种建模方法中线性校正方法明显好于非线性，且线性建模方法中PLS法和PCR法对蛋白质的预测性能均优于MLR法。而PLS和PCR法能够对光谱数据进行分解和筛选，剔除多重相关信息和无解释意义信息的干扰；SVR法作为非线性定量校正方法，其原理是通过升维方式使原样本空间中的非线性问题转化为特征空间中的线性关系进行建模，但是对于样本空间中线性问题使用SVR法会增加计算的复杂性，甚至引起“维数灾难”，导致模型的预测精度降低。结果显示，PLS法和PCR法建立的蛋白质校正模型的预测结果相差不大，RMSEP分别为0.125和0.132。因此，PLS法和PCR法均适合作为蛋白质检测的建模方法。
正交试验结果
在实际应用中，由于近红外吸收光弱且易受外界干扰，光谱中夹杂许多无用信息（包括噪声、背景等），为提高信噪比、建立稳定的校正模型，需要对光谱数据进行波段、预处理和建模方法的优选组合。本实验选取预处理方法、检测波段和建模方法3 个因素进行正交试验，以F值作为评价指标，F值越大表明模型的性能越好。经方差分析，以上3 个因素对模型的预测精度存在一定交互作用的影响。
结果显示，影响蛋白质模型预测精度的主次因素为建模方法＞检测波段＞预处理方法；根据极值分析结果，得到各因素的较好水平组合为A2B3C1，即建模方法为PCR法、检测波段为1 800～2 300 nm、预处理方法为MSC法，这与单因素试验结果并不完全一致，且A2B3C1没有出现于正交试验设计方案中。为了确定最佳建模条件，分别在A2B3C1、单因素试验以及正交设计试验中最优建模条件下对未知样品进行检测，结果表明A2B3C1条件下的预测效果最好，其R2、RMSEP分别为0.993、0.106。可见，经正交试验设计确定的建模条件能够进一步提高蛋白质模型的准确度，也表明正交设计能够有效优化复杂背景下的近红外建模条件。
结 论
建模条件的选择直接影响近红外光谱分析结果。本实验以复杂背景条件下的蛋白质为研究对象，采用近红外光谱技术对牛乳中蛋白质含量进行检测，通过正交试验设计优化近红外建模条件，即预处理方法、检测波段和建模方法。结果表明，与传统方法对比，该方法不仅避免了各建模条件存在交互作用的影响，覆盖了主要因素的各种组合，而且以较少实验次数得到了预测准确度较高的分析模型，为近红外蛋白质含量建模条件优选提供了一条有效的途径。
