一种基于x-射线荧光分析的农产品元素定量检测模型构建方法
【技术领域】
[0001] 本发明涉及一种X-射线荧光光谱分析技术，特别是涉及一种X-射线荧光光谱的 峰值法建模。
【背景技术】
[0002] 传统的能量色散型X-射线荧光光谱仪在建立统计模型时大都使用峰面积法。峰 面积法具有以下优点：不需准确知道样品的进样量，操作条件略为变动对试验结果的影响 相对较小，计算方便，适合对多种元素同时进行分析。
[0003] 曾国强等人基于模拟退火算法原理，建立一种新的寻峰模型算法，该算法利用模 拟退火寻找全局最优点的收敛特性，以Metropolis准则作为峰谷判断的基础，引入新的判 定标准和峰谷数组，同时从谱道址两端同时收敛，以收敛到同一最优解为终止条件。同时， 利用该算法与简单比较法、三阶导数法进行对比实验，结果证明该算法在X荧光光谱寻峰 中有较强能力，在实际生产中具有一定价值。刘合凡等人为解决X射线荧光测量过程中的 仪器刻度和校正难题，通过理论研宄并运用蒙特卡洛方法对X射线荧光仪进行建模，得到 了X射线探测器能谱响曲线。结果表明：模拟的X射线探测器能谱响应全谱与实测情况基 本一致，而主成份谱线则高度一致，特征峰半高宽相对误差均在5 %以内。方法为X射线荧 光测量的仪器刻度和校正提供重要的研宄信息，不仅为元素分析工作节省大量人力物力， 也为无标样X射线荧光分析提供了有力的理论支撑。方法也可用于X射线荧光分析过程中 的其他研宄领域。
[0004] 但在采用峰面积建模时，不同的元素在出峰位置上往往会出现重叠现象，且多数 仪器使用的软件在确定峰面积峰界的时候需要根据观察手动输入，因此容易造成较大误 差。对比X射线荧光光谱仪与电耦合等离子体发射光谱仪的数据发现，在测量一些含量较 低的元素，出现了较大的偏差，对比X射线荧光光谱仪与电耦合等离子体发射光谱仪测量 数据可以发现：用峰面积建立的模型的准确度相对较低，对含量值较小的元素产生的误差 较大。除此之外，现有X射线荧光光谱仪大多采用峰面积法建模，建模方法较为单一，选择 性较少。为此提出了用峰值法进行建模分析，以期为X射线荧光建模方法的多样性提供理 论基础。
【发明内容】
[0005] 本发明的目的在于提供一种基于X-射线荧光分析的农产品元素定量检测模型构 建方法，以实现对农产品元素含量进行准确快速地无损检测，且避免使用强酸造成污染。
[0006] 为了解决以上技术问题，本发明采用的具体技术方案如下：
[0007] -种基于X-射线荧光分析的农产品元素定量检测模型构建方法，其特征在于包 括以下步骤：
[0008] 步骤一，收集农产品样品得样品集，对得样品集进行粉碎、过80目筛、压片后，再 测量样品集的X射线荧光光谱谱图；
[0009] 步骤二，用电感耦合等离子体发射光谱仪测定样品集元素，得到样品集元素真值， 并与所述的X射线荧光光谱峰值进行关联；
[0010] 步骤三，对所述X射线荧光光谱谱图进行预处理、划分为校正集样品与验证集样 品并剔除异常样品；所述的预处理包括：微分处理、归一化、多元散射校正、中心化、标准正 态变量交换和光谱平滑，最终得到所述样品集的的X射线荧光光谱值；
[0011] 步骤四，通过主成分分析、偏最小二乘、人工神经网络，利用峰值法建立所述样品 集元素的真值与所述X射线荧光光谱值之间的校正模型；
[0012] 步骤五，通过验证集样品对所述校正模型的精度进行验证；通过所述的预处理及 所述的校正模型建立方式得到不同元素的建模方法及相关性分析，将预测相关系数最高的 模型作为农产品元素定量检测模型；
[0013] 所述元素为六1、卩、5、1(、〇&、]\111、卩6、附、(：11、211、？13中的任一种。
[0014] 所述的微分处理是对所述的X射线荧光光谱值求一阶导数和二阶导数；
[0015] 所述的归一化处理为对X射线荧光光谱值进行归一化，方法如下：设定数据集中 所有的数据向量长度相同，每个向量中所有元素的平方和相似，令：
[0017] 为使数据向量能归一化为N，将向量中每个元素乘以#/#，公式表示为：
[0019] 使向量为单位长度或单位面积，即归一化方法，使用该方法时向量的方差值为1， 即N取为1 ;
[0020] xik表示第i个样本的第k个测量数据，在数据矩阵中xik表示第i行第k列的数 据；
[0021] 所述中心化处理是从数据矩阵中的每个元素中减去该元素所在列的均值的运算， 如式⑶所示
[0023] yik表示第i个样本的第k个测量数据，在数据矩阵中yik表示第i行第k列的数 据；是n个样本的均值。
[0024] 1.根据权利要求1所述的一种基于X-射线荧光分析的农产品元素定量检测模型 构建方法，其特征在于：
[0025] 所述的标准正态变量变换计算公式如式（4)
[0027] Zi为第i个样品光谱的平均值；
[0028] k为波长点数，k= 1，2,…，m;
[0029] i为校正集样品数，i= 1，2,…，n。
[0030] 所述的光谱平滑计算公式如式（4)
[0031] 为，两者具有如下的关系：
[0033] k为被测样品数据点的序号，2m+l为滤波带宽的大小，m为滤波带宽的调节参数；
[0034] y*k为滤波值，yk为原始信号值。
[0035] 所述的人工神经网络中的人工神经元的输入和输出的关系如式（5)
[0037] ai为多个输入，i= 1，2......m；
[0038] Sj为单个输出；
[0039] 0」为阈值；
[0040] Wij为从神经元i到神经元j的连接权重因子；
[0041] f(X)为传递函数，或称激励函数。
[0042] 能量色散型重金属分析仪在建立统计模型时使用的峰面积法。峰面积法具有以下 优点：不需准确知道样品的进样量，操作条件略为变动对实验结果的影响相对较小，计算方 便，适合对多种元素同时进行分析。但由于不同的元素在出峰位置上往往会出现重叠现象， 且多数仪器使用的软件在确定峰面积峰界的时候需要根据观察手动输入，因此容易造成较 大误差。对比X-射线荧光光谱仪与电耦合等离子体发射光谱仪的数据发现，在测量一些含 量较低的元素，出现了较大的偏差。用峰面积建立的模型的准确度相对较低，尤其对含量较 低的元素产生的误差更大。除此之外，现有X-射线荧光光谱仪大多采用峰面积法建模，建 模方法较为单一，选择性较少。为此提出了用峰值法进行建模，以期为X-射线荧光建模方 法的多样性提供理论基础。
[0043] 1谱图预处理
[0044] 为了寻求建立更为准确的定量分析模型，对90个茶叶样品用电感耦合等离子体 发射光谱仪进行测定，得到了 90个茶叶样品的Al、P、S、K、Ca、Mn、Fe、Ni、Cu、Zn、Pb等11 种元素的真值。与提取的X-射线荧光的光谱值进行关联，以期得到预测效果更好的茶叶中 微量元素的X-射线荧光光谱的校正模型。
[0045] 光谱预处理对消除光谱中的噪声，部分消除或降低检测过程中的系统偏差，提高 X-射线荧光光谱信息的有效性有着至关重要的作用。研宄主要采用了微分处理、归一化、多 元散射校正、中心化、标准正态变量交换、光谱平滑等预处理方法对光谱进行了预处理，并 采用主成分分析及偏最小二乘的方法建立校正模型，并对预测集茶叶样品进行了预测及预 测精度分析，寻找出茶叶中各个待测元素的最佳光谱预处理方法。
[0046] 2样品集的划分与异常样品的剔除
[0047] 在用光谱数据进行建模时，异常样品的存在往往会对模型的预测精度产生较大的 影响。在进行样品集的划分和异常样品的剔除时，主要用的数据处理方法是主成分分析 (PCA)及偏最小二乘（PLS)。对茶叶的原始X-射线荧光光谱数据，采用江苏大学近红外319 课题组编写的NIRSAv4. 0数据处理系统进行分析，利用PCA或PLS，结合马氏距离识别出浓 度界外样本，适用于校正集和预测集样品。其中通过PCA方式划分的样本仅仅针对的是光 谱值，没有对化学值进行界外识别，而通过PLS方式划分的样本则是针对光谱值和化学值。
[0048] 3茶叶待测指标的校正模型建立方法分析
[0049] 茶叶建模分析主要采用了两种校正方法：主成分分析+人工神经网络（PCA+ANN)、 偏最小二乘+人工神经网络（PLS+ANN)。
[0050] 当神经元j有多个输入xji= 1，2, "'m)和单个输出y」时，输入和输出的关系可 表示为：
[0052] 其中叫为阈值，为从神经元i到神经元j的连接权重因子，f(x)为传递函数， 或称激励函数。传递函数f(x)可以选择线性函数，但通常选用非线性函数，常见的有：
[0053] (1)阶跃函数：
[0055] (2)Sigmoid(S型）函数：
[0057] 对于茶叶样品的建模，使用不同的预处理方法，经过反复试验，最终确定了相对较 佳的预处理方法，获得了较为准确的样品定标模型，并对预测集样品进行了预测和预测精 度分析。
[0058] 本发明具有有益效果。本发明通过采用峰值法进行建模分析，通过不同的预处理 及建模方式来建立校正模型，从而提高了模型精度，实现对农产品元素含量进行准确快速 地无损检测，且避免使用强酸造成污染。为重金属的快速无损检测提供理论依据。
【附图说明】
[0059] 图1是本发明的茶叶X射线荧光光谱图；
[0060] 图2是本发明的茶叶X射线荧光的一阶微分光谱图；
[0061]图3是本发明的茶叶X射线荧光的二阶微光图谱；
[0062] 图4是本发明的茶叶X射线荧光中心化后的光谱图；