• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 20
    • 下载费用:30 金币  

    重庆时时彩后二缩水器: 一种检测塑化剂含量的方法.pdf

    关 键 词:
    一种 检测 塑化剂 含量 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201110162558.6

    申请日:

    2011.06.16

    公开号:

    CN102230888A

    公开日:

    2011.11.02

    当前法律状态:

    撤回

    有效性:

    无权

    法律详情: 发明专利申请公布后的视为撤回IPC(主分类):G01N 21/31申请公布日:20111102|||实质审查的生效IPC(主分类):G01N 21/31申请日:20110616|||公开
    IPC分类号: G01N21/31; G01N21/35; G01N21/33 主分类号: G01N21/31
    申请人: 浙江大学
    发明人: 王靖岱; 陈美娟; 施小妹; 蒋斌波; 黄正梁; 廖祖维; 阳永荣
    地址: 310027 浙江省杭州市西湖区浙大路38号
    优先权:
    专利代理机构: 杭州天勤知识产权代理有限公司 33224 代理人: 胡红娟
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201110162558.6

    授权公告号:

    ||||||

    法律状态公告日:

    2014.06.25|||2011.12.14|||2011.11.02

    法律状态类型:

    发明专利申请公布后的视为撤回|||实质审查的生效|||公开

    摘要

    本发明公开了一种检测塑化剂含量的方法,包括如下步骤:(1)利用光谱仪获得多个已知塑化剂含量的校正样品的光谱,对校正样品的光谱进行预处理,得到预处理后的校正样品的光谱;(2)将步骤(1)中预处理后的校正样品的光谱通过数据处理方法进行处理,建立塑化剂含量的预测模型;(3)在与获得校正样品的光谱相同的色谱条件下利用光谱仪检测待测样品的光谱,对待测样品的光谱进行预处理,将预处理后的待测样品的光谱输入预测模型得到待测样品中塑化剂的含量;所述的塑化剂为邻苯二甲酸酯类塑化剂。该方法具有简便、快速准确、安全环保等特点,且无须对样品进行前处理,适用于工厂在线或离线质量控制和检验机构样品检验。

    权利要求书

    1.一种检测塑化剂含量的方法,其特征在于,包括如下步骤:
    (1)利用光谱仪获得多个已知塑化剂含量的校正样品的光谱,对校
    正样品的光谱进行预处理,得到预处理后的校正样品的光谱;
    (2)将步骤(1)中预处理后的校正样品的光谱通过数据处理方法进
    行处理,建立塑化剂含量的预测模型;
    (3)在与获得校正样品的光谱相同的色谱条件下利用光谱仪检测待
    测样品的光谱,对待测样品的光谱进行预处理,将预处理后的待测样品的
    光谱输入预测模型得到待测样品中塑化剂的含量;
    所述的塑化剂为邻苯二甲酸酯类塑化剂。
    2.如权利要求1所述的检测塑化剂含量的方法,其特征在于,所述
    的光谱仪为拉曼光谱仪或近红外光谱仪。
    3.如权利要求2所述的检测塑化剂含量的方法,其特征在于,所述
    的拉曼光谱仪的光源为可见光、近红外光或紫外光。
    4.如权利要求2所述的检测塑化剂含量的方法,其特征在于,所述
    的拉曼光谱仪的光谱范围为0~3600cm-1;
    或者,所述的近红外光谱仪的光谱范围为4000~14300cm-1。
    5.如权利要求1所述的检测塑化剂含量的方法,其特征在于,所述
    的光谱仪采用的分辨率为0.5cm-1~10cm-1。
    6.如权利要求1所述的检测塑化剂含量的方法,其特征在于,所述
    的预处理的方法为平滑、扣减、微分、归一化、标准化、标准正态变量变
    换、多元散射校正、分段多元散射校正、正交信号校正、去趋势、傅里叶
    变换、小波变换、净分析信号中的一种或两种以上。
    7.如权利要求1所述的检测塑化剂含量的方法,其特征在于,所述
    的数据处理方法为多元线性回归法、主成分回归法、偏最小二乘法、人工
    神经网络法、支持向量机法中的一种或两种以上。
    8.如权利要求1所述的检测塑化剂含量的方法,其特征在于,所述
    的邻苯二甲酸酯类塑化剂为邻苯二甲酸二(2-乙基己基)酯、邻苯二甲酸二
    丁酯、邻苯二甲酸甲苯基丁酯、邻苯二甲酸二异壬酯、邻苯二甲酸二异癸
    酯、邻苯二甲酸二正辛酯、邻苯二甲酸二异戊酯中的一种或两种以上。
    9.如权利要求1所述的检测塑化剂含量的方法,其特征在于,所述
    的校正样品和待测样品均为食品、包装材料、塑料制品中的一种。
    10.如权利要求9所述的检测塑化剂含量的方法,其特征在于,所述
    的食品为液体加工食品、粘稠胶体加工食品、固体加工食品或天然食品;
    所述的包装材料为聚乙烯、聚丙烯、聚氯乙烯或聚苯乙烯;
    所述的塑料制品为玩具、日常用品或医疗用品。

    说明书

    一种检测塑化剂含量的方法

    技术领域

    本发明涉及拉曼和近红外光谱检测领域,尤其涉及一种利用拉曼或近
    红外光谱检测邻苯二甲酸酯类塑化剂含量的方法。

    背景技术

    邻苯二甲酸酯(Phthalic?Acid?Esters,简称PAEs,别名钛酸酯)一般
    为无色透明的油状液体,难溶于水,易溶于甲醇、乙醇等。PAEs主要用
    作塑料袋及塑料制品中的塑化剂,增大产品可塑性和提高产品强度,也可
    起到软化的作用,被普遍应用于玩具、食品包装材料、医用血袋和胶管、
    乙烯地板和壁纸、清洁剂、润滑油、个人护理用品(如指甲油、头发喷雾
    剂、香皂和洗发液)等数百种产品中。但是该类化合物已经被证明对人的
    生殖系统有损害,具有生殖毒性、胚胎毒性和遗传毒性。欧盟早在
    2005/84/EC指令中即列出了六种邻苯二甲酸酯类物质邻苯二甲酸二(2-乙
    基己基)酯(DEHP)、邻苯二甲酸二丁酯(DBP)、邻苯二甲酸甲苯基丁酯
    (BBP)、邻苯二甲酸二异戊酯(DIPP)、邻苯二甲酸二异壬酯(DINP)、
    邻苯二甲酸二异癸酯(DIDP)、邻苯二甲酸二正辛酯(DNOP)进行限制,
    其中前三种DEHP、DBP、BBP不得用于儿童玩具和用品中,每种在塑料
    中的重量百分含量每种不得超过0.1%;后三种DINP、DIDP、DNOP不得
    用于能入口的儿童玩具及儿童类物品中,每种重量百分含量不得超过
    0.1%。

    近日,台湾塑化剂风波酿成重大食品安全?;?,卫生部专门为此公布
    第六批食品中可能违法添加的非食用物质和易滥用的食品添加剂名单,包
    括邻苯二甲酸酯类物质等共17种。食品中塑化剂的来源除了非法食品添
    加剂外,还包括环境污染和食品包装材料渗入等。喷涂涂料、塑料垃圾的
    焚烧、农用薄膜中塑化剂的挥发等使得PAEs进入自然环境,在土壤、水
    体和大气等环境介质中迁移、转化、积累,并且可以进入生物体内,在家
    禽海鲜蔬菜瓜果等天然食品内蓄积,从而造成极大的危害。在食品包装材
    料方面,PAEs与塑料基质之间没有形成化学共价键,因而在接触到包装
    食品中所含有的油脂等时便会溶出,污染食品。为了防止及避免塑化剂对
    人体造成危害,研究开发一种快速检测塑化剂含量的方法具有十分重大的
    意义。

    PAEs在环境样品中含量极低,并且基质复杂,所以分析测定困难。
    目前,PAEs分析检测方法主要有气相色谱法(GC)、液相色谱法(LC)、
    红外光谱法(IR)和薄层色谱法(TLC)和气相色谱-质谱联用法(GC-MS)。
    此外,随着高效液相色谱仪和液-质联用仪(LC-MS)的推广使用,对HPLC
    和LC-MS的研究也越来越多。GB/T?21911-2008邻苯二甲酸酯的测定方
    法采用气相色谱-质谱联用法(GC-MS),适用于食品中16种邻苯二甲酸
    酯含量的测定,对于含油脂样品的PEAs检测限为1.5mg/kg,不含油脂样
    品中的检测限为0.05mg/kg;主要原理是将不含油脂样品混匀后,加入适
    量水,用正己烷提纯,静置分层取上清液进行GC-MS分析,含油脂样品
    则需先用凝胶色谱柱进行净化后再按不含油脂样品的前处理方法进行前
    处理。上述提到的方法不仅需要对样品进行前处理,而且凝胶净化系统和
    分析系统中很多阀体及管路有大量的塑料或胶体部件,而这些部件在不同
    程度上都添加了各种塑化剂成分,从而造成了在塑化剂样品处理过程中带
    来污染,造成本底较高,影响分析的检出限和数据的可靠性。

    拉曼散射,又称拉曼效应,是光通过介质时由于入射光与分子运动相
    互作用而引起的频率发生变化的散射。目前,拉曼光谱仪主要可以分为两
    类:傅里叶拉曼光谱仪和激光拉曼光谱仪。傅里叶拉曼光谱借助Nd:YAG
    激光器(λ=1064nm)来激发,具有近红外线系统的干涉仪用于检测拉曼
    辐射,借助于陷波滤波器来抑制非波长偏移的瑞利辐射。激光拉曼光谱仪,
    又可称为分散型拉曼光谱仪,可应用不同类型的激光器激发拉曼辐射,比
    较常见的类型包括He:Ne激光器(λ=632nm)和半导体激光器(λ=785nm),
    通过使用格栅和CCD检测器来进行光谱衰变和检测,瑞利散射辐射是借
    助于陷波滤波器来阻断。拉曼光谱通过检测物的拉曼散射得到其相应的结
    构和组成信息,塑化剂含量的变化导致样品拉曼散射的变化,这为拉曼光
    谱检测塑化剂含量提供了理论基础。近红外光谱区与有机分子中含氢基团
    (OH、NH、CH)振动的合频和各级倍频的吸收区一致,通过扫描样品的
    近红外光谱,同样可以得到样品中有机分子含氢基团的特征信息。

    因此,发展能够快速、准确测定塑化剂含量的拉曼或近红外光谱检测
    方法,能够实现食品、包装材料及塑料制品生产过程中实时监控及质检部
    门的快速检测,对加工生产及人类健康有着重大的意义。

    发明内容

    本发明提供了一种利用拉曼或近红外光谱与数据处理方法相结合检
    测邻苯二甲酸酯类塑化剂含量的方法,该方法快速、准确,有效地解决了
    现有技术存在的局限,对实现生产过程中实时监控及质检部门的快速检测
    有着积极的意义。

    一种检测塑化剂含量的方法,包括如下步骤:

    (1)利用光谱仪获得多个已知塑化剂含量的校正样品的光谱,对校
    正样品的光谱进行预处理,得到预处理后的校正样品的光谱;

    (2)将步骤(1)中预处理后的校正样品的光谱通过数据处理方法进
    行处理,建立塑化剂含量的预测模型(即回归模型);

    (3)在与获得校正样品的光谱相同的色谱条件下利用光谱仪检测待
    测样品的光谱,对待测样品的光谱进行预处理,将预处理后的待测样品的
    光谱输入预测模型得到待测样品中塑化剂的含量;

    所述的塑化剂为邻苯二甲酸酯类(PAEs)塑化剂。所述的邻苯二甲酸
    酯类塑化剂可选用邻苯二甲酸二(2-乙基己基)酯(DEHP)、邻苯二甲酸二
    丁酯(DBP)、邻苯二甲酸甲苯基丁酯(BBP)、邻苯二甲酸二异壬酯(DINP)、
    邻苯二甲酸二异癸酯(DIDP)、邻苯二甲酸二正辛酯(DNOP)、邻苯二甲
    酸二异戊酯(DIPP)等中的一种或两种以上。

    所述的校正样品和待测样品均为食品、包装材料、塑料制品中的一种。

    所述的食品包括:茶饮料、果汁饮料、运动饮料等液体类加工食品;
    果酱、果浆、果冻等粘稠胶体类加工食品;胶囊锭状粉状等固体类加工食
    品;家禽、海鲜、蔬菜、瓜果等天然食品等中的一种。

    所述的包装材料包括聚乙烯(PE)、聚丙烯(PP)、聚氯乙烯(PVC)、
    聚苯乙烯(PS)等聚合物中的一种。

    所述的塑料制品包括玩具、日常用品(如指甲油、头发喷雾剂、香皂
    或洗发液)、医疗用品等中的一种。

    所述的光谱仪优选采用拉曼光谱仪或近红外光谱仪,可选用傅里叶拉
    曼光谱仪(FT-Raman光谱仪)、激光拉曼光谱仪、紫外激光拉曼光谱仪、
    傅里叶近红外光谱仪等现有色谱仪中的一种。通过拉曼光谱仪或近红外光
    谱仪获取液体类样品的光谱时,可采用浸入式测量或非浸入式测量。其中
    浸入式测量是指将探头等设备插入到液体样品中,获得光谱,浸入式测量
    无需专门采样,可以实现在线检测。非浸入式测量是指液体样品与光谱仪
    或探头不直接接触,而是通过透射获得光谱,例如通过光谱仪检测样品池
    中的液体就属于非浸入式测量的一种。本发明的检测方法,液体类样品的
    光谱可以通过浸入式测量也可以通过非浸入式测量获得,粘稠胶体类、固
    体类及天然食品及包装材料和塑料制品等非液体类样品的光谱通过非浸
    入式测量获得。所述的光谱仪使用不同的激发光源对于同一个基团而言,
    产生的拉曼位移位置不会变,只是强度不同而已。因此,光谱仪激发光源
    的选择依据不会损伤(烧掉、降解)样品及能够得到拉曼信号的标准。对
    于所述的塑化剂涉及到的基团,各种波长均能得到较好的结果。

    所述的拉曼光谱仪的光源为可见光、近红外光或紫外光。

    通过拉曼光谱仪或近红外光谱仪获取样品光谱时,设置光谱仪的分辨
    率为0.5cm-1~10cm-1为宜。当设置最高分辨率为0.5cm-1时,可以提高检
    测结果的准确性;当设置最高分辨率为10cm-1时,检测结果的准确性有
    所下降,但检测时间却大大缩短;当设置最高分辨率高于10cm-1时,光谱
    中一些有用信息会被掩盖,对检测结果产生一定影响。

    一般而言,扫描和分析检测物质(校正样品和待测样品)拉曼光谱时,
    选用的拉曼光谱范围为0~3600cm-1;扫描和分析检测物质近红外光谱时,
    选用的光谱范围为0~10000cm-1。PEAs的特征拉曼频率如表1所示。

    表1塑化剂特征拉曼频率表


    拉曼光谱仪或近红外光谱仪所采集的光谱除样品的自身信息外,还包
    括了其他无关信息和噪音,如电噪音、样品背景和杂散光等。因此,消除
    光谱数据无关信息和噪音的预处理方法变得十分关键和必要。本发明的检
    测方法,所述的预处理的方法可采用平滑、扣减、微分、归一化、标准化、
    标准正态变量变换、多元散射校正、分段多元散射校正、正交信号校正、
    去趋势、傅里叶变换、小波变换、净分析信号等现有分析方法中的一种或
    两种以上。

    所述的平滑可以提高分析信号的信噪比,最常用的方法是移动式平均
    平滑法和Savizky-Golay多项式平滑。

    所述的扣减是选择一吸收波长,将其吸收值对所有波长处的吸收值进
    行扣除,一般用于消除背景影响。

    所述的微分可以消除基线漂移、强化谱带特征、克服谱带重叠,是常
    用的光谱预处理方法,一阶微分可以去除同波长无关的漂移,二阶微分可
    以取出同波长线性相关的漂移。

    所述的标准正态变量变换(Standard?Normal?Variate,SNV)主要用于
    消除固体颗粒大小、表面散射以及光程变化对光谱的影响,其算法如下:

    X i , SNV = X i , k - X i Σ k = 1 m ( X i , k - X i ) 2 m - 1 ]]>

    式中,Xi,k为原始光谱,为第i样品光谱的平均值,k=1,2,...,m,
    m为波长点数;i=1,2,...,n,n为校正样品数,Xi,SNV为经SNV预处理
    后第i个样品的拉曼光谱。

    所述的多元散射校正(Multiplicative?Scatter?Correction,MSC)通过数
    学方法将光谱中的散射光信号与化学吸收信息进行分离,它可以去除光谱
    中样品的镜面反射及不均匀性造成的噪声,消除基线的不重复性,其具体
    算法如下:

    1)计算校正集样品光谱的平均光谱

    2)将Xi与进行线性回归,求取ai和bi;其中,Xi为
    第i个样品的光谱矩阵,为平均光谱矩阵,ai,bi为拉曼光谱矩阵的多元
    散射校正拟合系数;

    3)Xi,MSC=(Xi-bi)/ai,其中Xi,MSC为经MSC校正后的第i个样品的拉
    曼光谱矩阵;

    所述的分段多元散射校正是对多元散射校正的进一步优化,它能消除
    全波长范围内单体光谱和平均光谱的非线性散射。

    所述的正交信号校正是在建立定量校正模型前,将光谱阵用浓度阵正
    交,滤除光谱与浓度阵无关的信号,再进行多元校正,达到简化模型及提
    高模型预测能力的目的。

    所述的去趋势算法通常用于标准正态变量变换处理后的光谱,用来消
    除漫反射光谱的基线漂移。

    所述的傅里叶变换能够实现谱域函数与时域函数之间的转换,其实质
    是把原光谱分解成许多不同频率的正弦波的叠加和,它可以用来对光谱进
    行平滑去噪、数据压缩以及信息的提取。

    所述的小波变换能将化学信号根据频率的不同分解成多种尺度成分,
    并对大小不同的尺度成分采取相应粗细的取样步长,从而能够聚焦于信号
    中的任何部分。

    所述的净分析信号算法的基本思想与正交信号校正基本相同,都是通
    过正交投影除去光谱阵中与待测组分无关的信息。

    所述的数据处理方法采用的多元数据处理方法,可采用多元线性数据
    处理方法和/或多元非线性数据处理方法,如多元线性回归法、主成分回归
    法、偏最小二乘法、人工神经网络法、支持向量机法等现有数据处理方法
    中的一种或两种以上。其中多元线性回归法、主成分回归法、偏最小二乘
    法属于线性算法,人工神经网络法和支持向量机法属于非线性算法。

    所述的多元线性回归法是最早期定量分析常用的校正方法,但是存在
    诸多局限性,在运算过程中可能会出现病态矩阵,无法求逆;此外由于直
    接在波长集合上建模,获取有用信息时可能引入噪声,导致模型过拟合,
    降低了模型预测能力。

    所述的主成分回归法结合了对光谱的主成分分析和多元线性回归校
    正技术获得对复杂样品的定量模型,克服了多元线性回归的局限性,但是
    主成分分析只代表了校正集样本光谱自身的最大差异性,而这种差异性不
    一定完全是由于待测组分的不同引起的;另外主因子数是一个重要参数,
    如果选取的主因子数过多,有可能造成模型过拟合;反之,如果主因子数
    过少则会丢失有效光谱信息,造成模型拟合误差大。偏最小二乘法由Wold
    于1975年提出,现已成为光谱定量分析中应用最为广泛的一种方法。

    所述的偏最小二乘法(Partial?Least?Square,PLS)在从光谱数据矩阵中
    提取信息的同时,还对样品性质矩阵进行分解,这样既保留了较多的方差,
    又保留了较多的与因变量的相关性,从而在消去光谱矩阵复共线性的同
    时,使建立的回归模型仍能充分地反映出自变量与因变量之间的相关关
    系。通常偏最小二乘法采用非线性迭代偏最小二乘法(NIPALS)进行计
    算,它在分解自变量数据矩阵的同时,也在分解因变量数据矩阵,并设法
    使自变量数据矩阵中提取的成分尽可能靠近因变量数据矩阵中的成分,亦
    使它们的相关性尽量大。本发明方法中光谱矩阵X为自变量数据矩阵,浓
    度矩阵(即含量矩阵)Y为因变量数据矩阵,偏最小二乘法首先对光谱矩
    阵X和浓度矩阵Y进行分解,其模型为:

    X=SP+E

    Y=UQ+F

    式中S和U分别为光谱矩阵X和浓度矩阵Y的得分矩阵;P和Q分
    别为光谱矩阵X和浓度矩阵Y的载荷矩阵;E和F分别为光谱矩阵X和浓
    度矩阵Y的偏最小二乘法拟合残差矩阵。

    接下来,偏最小二乘法将S和U作线性回归:

    U=SB,B表示系数矩阵;

    B=(STS)-1STY,ST表示得分矩阵S的转置矩阵;

    在预测时,首先根据P求出预处理后的待测样品的光谱X未知得分S未
    知,然后由下式得到浓度预测值:

    Y未知=S未知BQ

    在实际的偏最小二乘法回归中,矩阵分解与回归是并为一步实现的,
    即光谱矩阵X和浓度矩阵Y矩阵的分解同时进行,并且将浓度矩阵Y的信
    息引入到光谱矩阵X的矩阵分解过程中,在每计算一个新主成分前,将光
    谱矩阵X的得分矩阵S与浓度矩阵Y的得分矩阵U进行交换,使得光谱
    矩阵X的主成分直接与浓度矩阵Y关联。这就克服了主成分回归法只对光
    谱矩阵X进行分解的缺点。同时,区别于主成分回归法,偏最小二乘法还
    可以根据不同的待测组分选取最佳的主成分数。一般偏最小二乘法选用
    1~5个偏最小二乘组分。

    所述的人工神经网络法是近些年发展起来的数学处理方法,它模仿人
    脑神经的活动力图建立脑神经活动的数学模型,即把对信息的存储和计算
    同时储存在神经单元中,所以在一定程度上神经网络可以模拟动物神经系
    统的活动过程,具有自学习、自组织、自适应能力、很强的容错能力、分
    布储存与并行处理信息的功能及高度非线性表达能力,这是其他传统方法
    所不具备的。支持向量机法是在20世纪90年代以后在Vapnik等建立的
    统计学系理论的基础上发展起来的一种新的统计学习方法,可以被应用于
    定性和定量分析。

    所述的支持向量机法可以有效克服人工神经网格训练速度慢、易于过
    拟合、预测能力差等缺点,在小样本、非线性和高维数据空间条件下的样
    本分类和回归问题中,表现出了传统方法不具有的优势。

    本发明与现有技术相比具有的优点如下:

    1)操作简便、能够快速地检测各种体系中塑化剂含量,且检测结果
    精度较高;

    2)无需对样品进行前处理,安全环保,适用于工厂在线或离线质量
    控制和检验机构样品检验;

    3)拉曼或近红外光谱检测装置灵敏度高,且检测范围宽;

    4)同时适用于浸入式测量或非浸入式测量;

    5)拉曼或近红外光谱检测可以实时在线;

    6)拉曼或近红外光谱检测装置的输出信号可直接与计算机联用,便
    于塑化剂含量的监控。

    附图说明

    图1是本发明的检测流程图;

    图2是三层BP神经网络结构示意图;

    图3是不同训练次数下奶粉中DNOP含量预测结果的平均相对误差示
    意图;

    图4是PVC中BBP含量预测结果示意图;

    图5是洗发液中DINP含量预测结果示意图;

    图6是炼乳中DEHP含量预测结果示意图;

    图7是钙片中DNOP含量预测结果示意图;

    图8是青菜中DBP含量预测结果示意图;

    图9是PE中DIPP含量预测结果示意图。

    具体实施方式

    实施例1

    采用德国Bruker公司生产的MultiRAM型傅里叶拉曼光谱仪检测不同
    DEHP含量的绿茶饮料样品的拉曼光谱,预测DEHP含量。该光谱仪选用
    激光光源波长为1064nm,设定光谱仪光谱扫描范围0~3600cm-1,采用浸
    入式测量。一系列不同DEHP含量的样品(校正样品和待测样品)是通过
    在市场销售的某品牌绿茶中(不含PEAs化合物)加入一定量DEHP,充
    分混合配置得到。样品中DEHP含量的具体值分别为0.520mg/kg、0.617
    mg/kg、0.713mg/kg、0.810mg/kg、0.907mg/kg、1.003mg/kg、1.100mg/kg、
    1.197mg/kg和1.294mg/kg。具体步骤如下:

    (1)利用FT拉曼光谱仪获得校正样品(DEHP含量:0.520mg/kg、
    0.617mg/kg、0.810mg/kg、1.003mg/kg、1.197mg/kg和1.294mg/kg)绿
    茶饮料的拉曼光谱,通过多元散射校正进行校正样品光谱预处理,其具体
    过程如下:

    1)计算所需校正样品绿茶饮料拉曼光谱集的平均光谱:

    X = Σ i = 1 8 X i 8 ]]>

    式中,Xi为第i个样品的拉曼光谱矩阵,为平均光谱矩阵。

    2)对Xi与进行线性回归:

    X i = a i X + b i ]]>

    式中,ai,bi为拉曼光谱矩阵的多元散射校正拟合系数。

    3)对每一条光谱作多元散射校正:

    X i = ( X i - b i ) a i ]]>

    式中,Xi’为经多元散射校正后的第i个样品的拉曼光谱矩阵。

    (2)根据方差解释,选择两个PLS成分进行偏最小二乘法分析,得
    到绿茶饮料样品中DEHP含量的回归模型,其过程如下:

    1)程序开始,计数器h置为1:h=1;

    2)将校正样品的含量矩阵Y送入列向量ss中,作为其初值:ss=Y;

    3)将校正样品的光谱矩阵X投影于列向量ss上:uh=XTss/(ssTss),
    uh、XT、ssT分别表示载荷向量、X和ss的转置向量;

    4)将向量uh归一化:uh=uh/||uh||;

    5)将校正样品的光谱矩阵X投影于行向量上,得到X的第h个
    PLS得分向量th:

    6)计算校正样品的光谱矩阵X的载荷向量ch:thT
    表示th的转置向量;

    7)将列向量ss对得分向量th进行回归:bh、sh分别
    表示载荷向量和得分向量;

    8)从校正样品的光谱矩阵X中除去第h个PLS成分:
    chT表示ch的转置向量;

    9)从校正样品的含量矩阵Y中除去回归项:Y=bhthss;

    10)计数器h加1,回到步骤3),再次开始计算;

    最终得到绿茶饮料样品中DEHP含量的预测模型:

    wt(mg/kg)=-6.03+1191.97×t1+7542.37×t2

    (3)对3个待测绿茶饮料样品中DEHP含量进行预测,其步骤如下:

    1)对待测样品的拉曼谱图进行多元散射校正预处理,得到预处理后
    的待测样品的拉曼光谱;

    2)将预处理后的待测样品的拉曼光谱代入绿茶饮料样品中DEHP含
    量的预测模型便可得出待测样品中的DEHP含量分别为0.695mg/kg、
    0.872mg/kg和1.080mg/kg,与3个待测绿茶饮料样品中DEHP含量真实
    值(真实值根据GB/T?21911-2008测得)的平均偏差为2.98%。

    采用不同谱图预处理方法对最终检测结果的影响如表2所示:

    表2不同预处理方法对最终检测结果的影响比较


    实施例2

    采用法国JDbin-yvon公司生产的LabRamHRUV型紫外激光拉曼光谱
    仪检测不同DINP含量的草莓果酱样品的拉曼光谱,预测DINP含量。该
    光谱仪选用激光光源波长为325nm,设定光谱仪光谱扫描范围0~3600
    cm-1,采用非浸入式测量。一系列不同DINP含量的样品(校正样品和待
    测样品)是通过在市场销售的某品牌草莓果酱中(不含PEAs化合物)加
    入一定量DINP,充分混合配置得到的。样品中DINP含量的具体值分别
    为0.005mg/kg、0.053mg/kg、0.102mg/kg、0.150mg/kg、0.199mg/kg、
    0.247mg/kg、0.295mg/kg、0.344mg/kg和0.392mg/kg。具体步骤如下:

    (1)利用紫外激光拉曼光谱仪获得校正样品(DINP含量:0.005
    mg/kg、0.053mg/kg、0.150mg/kg、0.247mg/kg、0.344mg/kg和0.392mg/kg)
    草莓果酱的拉曼光谱,通过标准正态变量变换对校正样品原始光谱进行预
    处理,其具体算法如下:

    X i , SNV = X i , k - X i Σ k = 1 m ( X i , k - X i ) 2 m - 1 ]]>

    式中,Xi,k为原始拉曼光谱,为第i样品光谱的平均值。k=1,2,...,
    m,m为波长点数;i=1,2,...,n,n为校正集样品数。

    (2)利用主成分回归法拟合得到草莓果酱中DINP含量的预测模型。
    主成分回归法主要包括主成分分析和多元线性回归两个阶段,实现的具体
    步骤为:

    1)数据预处理(对X、Y按列标准化等);

    2)求相关系数矩阵R;

    3)求R的特征根(λi)和特征向量P;

    4)根据特征根值从大到小重排特征根和特征向量P;

    5)计算主成分贡献率;

    6)计算累计贡献率,当主成分累计贡献率>85%,记录主成分个数k;

    7)计算主成分矩阵:T=XRk;

    8)计算回归系数。

    通过商业软件Matlab中提供的主成分分析程序实现多元数据处理,
    其程序文件名为princomp.m,对预处理后的校正样品的拉曼光谱X进行主
    成分分析:

    [Coeff?Score]=princomp(X)

    得到主成分载荷(Coeff)和主成分得分(Score)。选取校正样品的拉曼谱
    图的主成分得分前两项与DINP含量进行关联,通过主成分回归拟合得到
    草莓果酱中DINP含量的预测模型;

    wt(mg/kg)=0.89+1491.57×Score?1-464.94×Score?2

    (3)对3个待测样品草莓果酱中DINP含量进行预测,其步骤如下:

    1)对待测样品拉曼谱图进行标准正态变量变换预处理,其具体算法
    同校正样品中预处理的算法;

    2)经上述预处理后的待测样品的拉曼光谱与校正样品主成分载荷前
    两项计算得到新的主成分得分;

    3)将待测样品主成分得分代入草莓果酱中DINP含量的回归模型便
    可得出待测样品中NIDP含量分别为0.098mg/kg、0.193mg/kg、0.286
    mg/kg,与3个待测样品草莓果酱中DINP含量的真实值(真实值根据GB/T
    21911-2008测得)的相对误差分别为3.50%、2.61%和3.10%。

    实施例3

    采用美国Thermo?Fisher公司生产的DXR智能拉曼光谱仪检测不同
    DNOP含量的奶粉样品的拉曼光谱,预测DNOP含量。该光谱仪选用激光
    光源波长为780nm,设定光谱仪光谱扫描范围0~3600cm-1,分辨率为0.5
    cm-1,采用非浸入式法测量。一系列不同DNOP含量的奶粉样品(校正样
    品和待测样品)是通过在市场销售的某品牌奶粉中(不含PEAs化合物)
    加入一定量DNOP,充分混合配置得到。样品中DNOP含量的具体值分别
    为0.008mg/kg、0.057mg/kg、0.107mg/kg、0.156mg/kg、0.205mg/kg、
    0.255mg/kg、0.304mg/kg、0.353mg/kg和0.402mg/kg。具体步骤如下:

    (1)利用DXR智能拉曼光谱仪获得校正样品(DNOP含量为:0.008
    mg/kg、0.057mg/kg、0.156mg/kg、0.255mg/kg、0.353mg/kg和0.402mg/kg)
    奶粉的拉曼光谱,对校正样品的拉曼谱图进行微分法预处理,得到预处理
    后的校正样品拉曼光谱矩阵X,其中微分法的过程如下:

    X=d(XO)/d(g)

    式中XO为原始拉曼光谱,X为预处理后的拉曼光谱,g为微分窗口宽
    度。

    (2)通过人工神经网络法拟合得到奶粉中DNOP含量的预测模型,
    其具体过程如下:

    1)根据三层BP神经网络结构,如图2所示,首先确定输入层。若将
    拉曼光谱直接作为网络的输入,则网络规?;岱浅E哟?。因此,选择主成
    分分析(PCA)得到的主成分得分(其过程参见实施例2)作为输入。根
    据校正样品的预测效果(如表3),确定最佳主成分得分项为3;

    表3主成分得分项对最终检测结果的影响

    ??1
    ??2
    ??3
    ??4
    ??5
    ??最大偏差
    ??4.72%
    ??3.14%
    ??2.11%
    ??2.23%
    ??2.08%
    ??平均偏差
    ??3.58%
    ??2.13%
    ??1.64%
    ??1.70%
    ??1.68%

    2)BP神经网络是通过输入与输出之间不断循环迭代建立起的一种基
    于连接权重的非线性关系。在取3个主成分的情况下,比较训练次数对标
    定样品预测结果的影响(如图3),发现训练次数为2000次时预测效果最
    佳;

    3)由于其他网络参数的确定缺乏理论依据,因此经多次实验反复验
    证,最终确定隐含层节点数为4,控制误差为0.0001,学习速率为0.05,
    态势因子取0.1,输出节点数为1;

    4)将主成分分析的3个主成分得分作为输入,奶粉中DNOP含量作
    为输出,通过BP神经网络进行反复训练建立PCA-BP校正集模型,即得
    到预测模型。

    (3)对3个待测样品奶粉中DNOP含量进行预测,其步骤如下:

    1)对待测样品的拉曼谱图进行微分变换预处理,其微分法的过程同
    校正样品中微分法的过程的过程;

    2)将上述预处理后的待测样品的拉曼光谱代入奶粉中DNOP含量的
    预测模型便可得出待测样品奶粉中DNOP含量分别为0.105mg/kg,0.201
    mg/kg和0.298mg/kg,与3个待测样品奶粉中DINP含量的真实值(真实
    值根据GB/T?21911-2008测得)的相对误差分别为1.87%、1.95%和1.97%。

    实施例4

    采用美国Thermo?Fisher公司生产的Antaris?II傅里叶近红外光谱仪检
    测不同DBP含量的猪肉样品的近红外光谱,预测DBP含量。该光谱仪设
    定光谱仪光谱扫描范围0~10000cm-1,分辨率为8cm-1,扫描次数16次,
    采用非浸入式法测量。一系列不同DBP含量的猪肉样品(校正样品和待
    测样品)通过在铰碎了的市场销售的猪肉中(不含PEAs化合物)加入一
    定量DBP,充分混合配置得到。样品中DBP含量的具体值分别为0.002
    mg/kg、0.054mg/kg、0.106mg/kg、0.159mg/kg、0.211mg/kg、0.263mg/kg、
    0.315mg/kg、0.367mg/kg和0.420mg/kg。具体步骤如下:

    (1)利近红外光谱仪获得校正样品(DBP含量:0.002mg/kg、0.054
    mg/kg、0.159mg/kg、0.263mg/kg、0.367mg/kg和0.420mg/kg)猪肉的近
    红外光谱,对校正样品的谱图进行傅里叶变换预处理,得到处理后的校正
    样品的光谱矩阵X。傅里叶变换预处理可利用Matlab中的fftshift.m程序
    实现:

    X=fftshift(XO)

    式中XO为原始光谱,X为预处理后的光谱。

    (2)通过最小二乘法和支持向量机法进行多元数据处理,得到猪肉
    中DBP含量的预测模型。具体步骤为:

    1)将函数拟合问题描述为以下最优化问题:

    min w , e J ( w , e ) = 1 2 w T w + 1 2 γ Σ k = 1 N e k 2 ]]>


    其中xk∈Rm为输入数据,yk∈R为输出数据,是将输入
    数据映射到高位特征空间的函数,权值向量为误差变量和偏置值
    分别为ek∈R及b∈R;γ>0为惩罚系数,可以调整误差所起的作用,WT表
    示W的转置向量;

    2)最优化问题可以变换到对偶空间加以解决,得到Lagrange函数,


    其中Lagrange乘子αk∈R,称为支持值;J表示1)中的最优化问题;

    3)对各变量求偏导,得到以下的条件等式:


    L b = 0 Σ k = 1 N α k = 0 , ]]>

    L e k = 0 α k = γ e k , k = 1,2 , . . . , N ]]>


    4)消去变量w和e后,可得线性方程组:

    0 1 T 1 Ω + γ - 1 1 b α = 0 y ]]>

    其中y=[y1,...,yN]T,1=[1,...,1]T,α=[α1,...,αN]T,Ω={Ωkl|k,l=1,...,N},
    K(xk,xl)为RBF核函数。

    5)得到最小二乘支持向量机法拟合模型为:

    y ( x ) = Σ k = 1 N α k K ( x k , x l ) + b ]]>

    其中α,b为4)中线性方程组的解。

    (3)对3个待测样品猪肉中DBP含量进行预测,其步骤如下:

    1)对待测样品的近红外谱图进行傅里叶变换预处理,其傅里叶变换
    预处理的过程同校正样品的傅里叶变换预处理的过程;

    2)将上述处理后的待测样品的光谱代入猪肉中DBP含量的预测模型
    便可得出待测样品猪肉中DBP含量分别为0.105mg/kg,0.209mg/kg和
    0.312mg/kg,与3个待测样品猪肉中DBP含量的真实值(真实值根据GB/T
    21911-2008测得)的相对误差分别为1.21%、0.96%和1.16%。

    实施例5

    采用美国Thermo?Fisher公司生产的DXR智能拉曼光谱仪检测一系列
    不同BBP含量的聚氯乙烯(PVC)样品的拉曼光谱,预测BBP含量。该
    光谱仪选用激光光源波长为532nm,设定光谱仪光谱扫描范围0~3600
    cm-1,分辨率为0.5cm-1,采用非浸入式法测量。

    采用傅里叶变换法对校正样品光谱进行预处理,方法与实施例4一致;
    通过偏最小二乘法进行多元数据处理,方法与实施例1一致,选择2个主
    成分建立预测模型;对待测样品拉曼光谱进行傅里叶变换预处理,代入预
    测模型。预测结果如图4所示。从图4中可看出,针对同一样品,本发明
    检测方法测得的PVC样品中BBP的预测含量与根据GB/T?21911-2008
    测得的PVC样品中BBP的真实含量基本一致,误差在允许的范围之内。

    实施例6

    采用法国JDbin-yvon公司生产的LabRamHRUV型紫外激光拉曼光谱
    仪检测不同DINP含量的洗发液样品的拉曼光谱,预测DINP含量。该光
    谱仪选用激光光源波长为514nm,设定光谱仪光谱扫描范围0~3600cm-1,
    采用非浸入式测量。实施方法与实施例1一致。预测结果如图5所示。从
    图5中可看出,针对同一样品,本发明检测方法测得的洗发液样品中DINP
    的预测含量与根据GB/T?21911-2008测得的洗发液样品中DINP的真实
    含量基本一致,误差在允许的范围之内。

    实施例7

    采用德国Bruker公司生产的MultiRAM型傅里叶拉曼光谱仪检测一系
    列不同DEHP含量的炼乳样品的拉曼光谱,预测DEHP含量。该光谱仪选
    用激光光源波长为1064nm,设定光谱仪光谱扫描范围0~3600cm-1,采用
    浸入式测量。实施方法与实施例2一致。预测结果如图6所示。从图6中
    可看出,针对同一样品,本发明检测方法测得的炼乳样品中DEHP的预测
    含量与根据GB/T?21911-2008测得的炼乳样品中DEHP的真实含量基本
    一致,误差在允许的范围之内。

    实施例8

    采用美国Thermo?Fisher公司生产的DXR智能拉曼光谱仪检测不同
    DNOP含量的口服钙片样品的拉曼光谱,预测DNOP含量该光谱仪选用激
    光光源波长为780nm,设定光谱仪光谱扫描范围0~3600cm-1,分辨率为
    0.5cm-1,采用非浸入式法测量。实施方法与实施例3一致。预测结果如图
    7所示。从图7中可看出,针对同一样品,本发明检测方法测得的口服钙
    片样品中DNOP的预测含量与根据GB/T?21911-2008测得的口服钙片样
    品中DNOP的真实含量基本一致,误差在允许的范围之内。

    实施例9

    采用美国Thermo?Fisher公司生产的Antaris?II傅里叶近红外光谱仪检
    测不同DBP含量的青菜样品的近红外光谱,预测DBP含量。该光谱仪设
    定光谱仪光谱扫描范围0~10000cm-1,分辨率为8cm-1,扫描次数16次,
    采用非浸入式法测量。实施方法与实施例4一致。预测结果如图8所示。
    从图8中可看出,针对同一样品,本发明检测方法测得的青菜样品中DBP
    的预测含量与根据GB/T?21911-2008测得的青菜样品中DBP的真实含量
    基本一致,误差在允许的范围之内。

    实施例10

    采用美国Thermo?Fisher公司生产的DXR智能拉曼光谱仪检测一系列
    不同DIPP含量的聚乙烯(PE)样品的拉曼光谱,预测DIPP含量。该光
    谱仪选用激光光源波长为780nm,设定光谱仪光谱扫描范围0~3600cm-1,
    分辨率为0.5cm-1,采用非浸入式法测量。实施方法与实施例5一致。预
    测结果如图9所示。从图9中可看出,针对同一样品,本发明检测方法测
    得的PE样品中DIPP的预测含量与根据GB/T?21911-2008测得的PE样
    品中DIPP的真实含量基本一致,误差在允许的范围之内。

       内容来自专利网重庆时时彩单双窍门 www.4mum.com.cn转载请标明出处

    关于本文
    本文标题:一种检测塑化剂含量的方法.pdf
    链接地址://www.4mum.com.cn/p-5870782.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 欠发达地区如何赚钱 广东快乐10分万能码 001期香港六合彩 闲来陕西麻将俱乐部 双色球蓝球2017142 拉菲彩票首页 江苏11选5走势图 58彩票安卓 江西时时彩开奖 3d试机号 赛马会官方网 淘宝做墙纸赚钱吗 吉林11选5开奖结果公告 U宝娱乐首页 生肖时时彩计划 电子游戏网站大全注册