CN106198447A - 基于近红外光谱技术的复混肥主要成分无损定量检测方法 - Google Patents

基于近红外光谱技术的复混肥主要成分无损定量检测方法 Download PDF

Info

Publication number
CN106198447A
CN106198447A CN201610551981.8A CN201610551981A CN106198447A CN 106198447 A CN106198447 A CN 106198447A CN 201610551981 A CN201610551981 A CN 201610551981A CN 106198447 A CN106198447 A CN 106198447A
Authority
CN
China
Prior art keywords
mixed fertilizer
sample
chemical mixed
spectrum
content
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610551981.8A
Other languages
English (en)
Inventor
汪玉冰
林志丹
王儒敬
宋良图
张正勇
汪六三
鲁翠萍
刘洋
朱利凯
郭红燕
高钧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN201610551981.8A priority Critical patent/CN106198447A/zh
Publication of CN106198447A publication Critical patent/CN106198447A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/28Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/286Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/28Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/286Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
    • G01N2001/2866Grinding or homogeneising

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供一种基于近红外光谱技术的复混肥主要成分无损定量检测方法,首先对复混肥样本集进行采集,然后将每个样本分为两份,一份利用标准化学方法测定总氮含量、总磷含量和总钾含量,另一份利用近红外光谱仪测量其反射光谱,再对光谱进行预处理,选取校正样本集,利用该校正样本集的光谱数据和标准化学方法所测得的含量数据,通过多元回归算法建立校正模型或者优选出特征波长直接建立多元校正模型,通过该校正模型对复混肥样品的主要成分含量进行测定。本发明大大缩短了分析测量时间,过程中无需大量反应试剂,不仅节省了大量的人力物力,同时不会造成环境污染,能有效解决传统复混肥生产中传统离线或抽样检测耗时长、效率低的问题。

Description

基于近红外光谱技术的复混肥主要成分无损定量检测方法
技术领域
本发明涉及化肥生产质量控制技术领域,具体涉及一种基于近红外光谱技术的成分含量测定方法,用于快速、无损、定量测量复混肥的多种主要成分。
背景技术
化肥是现代农业中必不可少的重要物资,施用化肥可有效增加粮食产量。研究表明,多种肥料混合施用的效果最佳,因此目前大部分化肥生产厂家都是生产含有两种以上作物生长必需营养成分的复混肥。由于生产厂家良莠不齐以及工艺上的问题,导致复混肥中营养成分并不是完全固定的,而是在某一范围内波动,因此对复混肥主要有效成分快速、准确地测量,对化肥质量把控及后续精准施肥具有重要意义。
复混肥的主要成分包括总氮、总磷和总钾等,传统的检测方法多在实验室中进行,其中总氮含量测定主要使用凯氏定氮法、总磷含量测定主要采用磷钼酸喹啉重量法,总钾含量测定主要采用四苯硼酸钾重量法等,这些化学方法涉及称重、溶解、消化、蒸馏、滴定以及计算等一系列复杂的操作步骤,不仅费时费力、成本昂贵,还会对实验操作人员造成一定的危险性并造成环境污染。
近红外光谱位于可见光和中红外波段之间,波长范围为780~2526nm。近红外光谱属于分子振动光谱,产生于共价化学键非谐能级振动,是非谐振动的倍频和组合频。原则上能产生近红外光谱的物质都有用近红外光谱分析的可能,并且近红外光谱分析方法展现出一系列优点,诸如快速、简便、低成本、非破坏性和多组分同时测定等。随着化学计量学的快速发展,近红外光谱分析技术在石油化工、食品、药物等领域得到了广泛的应用,这些都为我们将该技术应用于复混肥的主要成分检测提供了可行性。
发明内容
本发明提供一种可基于近红外光谱技术的复混肥主要成分无损定量检测方法,以有效满足复混肥生产过程中多种主要成分的快速测定,有效监管复混肥质量。
为解决上述技术问题,本发明采用如下技术方案:
一种基于近红外光谱技术的复混肥主要成分无损定量检测方法,包括如下步骤:
1)复混肥样本集的采集:采集不同批次的复混肥样本,标明采样时间和生产批次,再研磨成粉末状;
2)复混肥样本主要成分含量测定和光谱采集:将每个样本分为两份,一份利用标准化学方法测定总氮含量、总磷含量和总钾含量,其中总氮含量的测定采用凯氏定氮法,总磷含量测定采用磷钼酸喹啉重量法,总钾含量测定使用四苯硼酸钾重量法;另一份利用近红外光谱仪测量其反射光谱,多次重复测量(3~4次),取平均光谱近似作为该样本的标准光谱;
3)对光谱进行预处理,即近红外光谱仪所采集的光谱信息除了有用信息之外,还包括仪器本身的噪声、背景光和杂散光引起的噪声、样品颗粒大小引起的光程误差等,因此,在后续利用化学计量学建立模型时,必须使用合适的预处理方法;优选地,所述光谱预处理方法包括平滑、导数、标准正态变量转换、多元散射校正和正交信号校正;
4)校正模型的建立:选取校正样本集,利用该校正样本集的光谱数据和标准化学方法所测得的含量数据,通过多元回归算法建立校正模型或者优选出特征波长直接建立多元校正模型;
优选地,所述校正样本集的选取方法包括Kennard-Stone方法,即基于样品光谱的欧氏距离选择校正样本集;或者先将光谱进行主成分分析,根据主成分得分作为特征变量进行Kennard-Stone方法选择校正样本集;
优选地,所述多元回归算法包括多元线性回归方法和多元非线性回归方法,其中多元线性回归方法包括主成分回归和偏最小二乘法,多元非线性回归方法包括人工神经网络和支持向量回归;
优选地,所述特征波长选择方法包括相关系数法、逐步回归法、无信息变量消除法以及遗传算法;
5)复混肥样品的主要成分含量的测定:首先扫描复混肥样品的近红外光谱,将所采集到的复混肥样品的近红外光谱带入校正模型,计算即可得到复混肥样品的主要成分含量。
本专利所涉及的分析对象是复杂的样品体系,不能通过人工配制获得校正样品,必须收集实际样本,在生产线上所采集到的复混肥样本含有大量的重复样本,因此有必要挑选出具有代表性的样品建立校正模型;通过所选取的代表性样品进行建模不仅可以减少模型库的存储空间、提高建模速度,更能通过较少的样品提高模型的适用范围,便于模型更新和维护;因此在步骤4)和步骤5)之间还需要设置校正模型的验证,具体方法如下:
选取已知成分含量的多个预测集样本,将预测集样本的近红外光谱进行预处理后代入所建校正模型中,分别计算预测集样本的总氮、总磷和总钾含量,分别使用相关系数(R)、预测均方根误差(RMSEP)和相对分析误差(RPD)作为标准对所建立的回归模型预测性能进行评价。
由以上技术方案可知,本发明具有简单、快速、准确、高效、无污染的特点,该方法大大缩短了分析测量时间,过程中无需大量反应试剂,不仅节省了大量的人力物力,同时不会造成环境污染;测量结果较为准确,误差小,能有效解决传统复混肥生产中传统离线或抽样检测耗时长、效率低的问题,为复混肥生产过程中在线质量监控提供技术支撑。
附图说明
图1为本发明的原理框图。
具体实施方式
本发明是通过以下具体实施方式实现的,但本发明并不限于此。
实施例1
1.复混肥样本的收集与预处理。
采集不同批次的174个复混肥样本,每份样品控制在0.5kg左右,取样后放入透明自封袋,标明采样时间和生产批次,再研磨成粉末状。将每个样品分为两份,一份送往化学实验室进行测量,另一份样品利用近红外光谱仪测量其反射光谱。
2.复混肥样本主要成分含量的测定。
分别使用凯氏定氮法、磷钼酸喹啉重量法和四苯硼酸钾重量法测定复混肥中的总氮含量、总磷含量和总钾含量。
3.复混肥样本光谱的采集。
使用海洋光学公司的USB4000可见光谱仪和滨松公司的C9914GB近红外光谱仪以及卤钨灯光源采集复混肥样本的漫反射光谱。整个光谱范围为342~2221nm,每个样本重复测量3~4次,取平均光谱作为该样品的标准光谱。
4.校正样本集的选择。
首先我们剔除异常样本,采用主成分分析与马氏距离相结合的方法(PCA-MD),取阈值为马氏距离均值加上三倍的标准偏差,则编号为8和136的样本判断为异常样本,将其删除,剩余172个样本。再采用KS方法挑选124个样品作为校正样本集,剩下的43个样本作为预测集。
5.校正模型的建立。
对光谱进行不同预处理,再使用主成分回归方法(PCR),针对不同的预测量分别采用留一交叉验证方法,根据预测残差平方和选取最佳主成分数。
6.校正模型的验证。
选取已知成分含量的43个预测集样本,将预测集样本的近红外光谱进行相同的预处理(多元散射校正和平滑)后代入所建校正模型中,分别计算预测集样本的总氮、总磷和总钾含量。分别使用相关系数(R)、预测均方根误差(RMSEP)和相对分析误差(RPD)作为标准对所建立的回归模型预测性能进行评价。
表1为采用PCR方法所建复混肥主要成分含量模型的分析结果,其中对复混肥中总氮、总磷和总钾含量进行了计算。
实施例2
1.复混肥样本的收集与预处理同实施例1。
2.复混肥样本主要成分含量的测定同实施例1。
3.复混肥样本光谱的采集同实施例1。
4.校正样本集的选择同实施例1。
5.校正模型的建立。
对光谱进行不同预处理,再使用遗传算法(GA)优选特征波长,根据所优选出的特征波长建立多元线性回归模型。
6.校正模型的验证同实施例1。
表2是基于GA方法所建立复混肥主要成分含量模型的分析结果。
实施例3
1.复混肥样本的收集与预处理同实施例1。
2.复混肥样本主要成分含量的测定同实施例1。
3.复混肥样本光谱的采集同实施例1。
4.校正样本集的选择同实施例1。
5.校正模型的建立。
对光谱进行不同预处理,再使用连续投影算法(SPA)优选特征波长,根据所优选出的特征波长建立多元线性回归模型。
6.校正模型的验证同实施例1。
表3是基于SPA方法所建立复混肥主要成分含量模型的分析结果。
综上,上述方法所建预测模型效果较好,可用于复混肥主要成分的在线检测。
以上所述实施方式仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案作出的各种变形和改进,均应落入本发明的权利要求书确定的保护范围内。

Claims (7)

1.一种基于近红外光谱技术的复混肥主要成分无损定量检测方法,其特征在于,包括如下步骤:
1)复混肥样本集的采集:采集不同批次的复混肥样本,标明采样时间和生产批次,再研磨成粉末状;
2)复混肥样本主要成分含量测定和光谱采集:将每个样本分为两份,一份利用标准化学方法测定总氮含量、总磷含量和总钾含量;另一份利用近红外光谱仪测量其反射光谱,多次重复测量,取平均光谱近似作为该样本的标准光谱;
3)对光谱进行预处理;
4)校正模型的建立:选取校正样本集,利用该校正样本集的光谱数据和标准化学方法所测得的含量数据,通过多元回归算法建立校正模型或者优选出特征波长直接建立多元校正模型;
5)复混肥样品的主要成分含量的测定:首先扫描复混肥样品的近红外光谱,将所采集到的复混肥样品的近红外光谱带入校正模型,计算即可得到复混肥样品的主要成分含量。
2.根据权利要求1所述的复混肥主要成分无损定量检测方法,其特征在于,所述校正样本集测量主成分含量的标准化学方法为:总氮含量测定采用凯氏定氮法,总磷含量测定采用磷钼酸喹啉重量法,总钾含量测定使用四苯硼酸钾重量法。
3.根据权利要求1所述的复混肥主要成分无损定量检测方法,其特征在于,所述光谱预处理方法包括平滑、导数、标准正态变量转换、多元散射校正和正交信号校正。
4.根据权利要求1所述的复混肥主要成分无损定量检测方法,其特征在于,所述校正样本集的选取方法包括Kennard-Stone方法,即基于样品光谱的欧氏距离选择校正样本集;或者先将光谱进行主成分分析,根据主成分得分作为特征变量进行Kennard-Stone方法选择校正样本集。
5.根据权利要求1所述的复混肥主要成分无损定量检测方法,其特征在于,所述多元回归算法包括多元线性回归方法和多元非线性回归方法,其中多元线性回归方法包括主成分回归和偏最小二乘法,多元非线性回归方法包括人工神经网络和支持向量回归。
6.根据权利要求1所述的复混肥主要成分无损定量检测方法,其特征在于,所述特征波长选择方法包括相关系数法、逐步回归法、无信息变量消除法以及遗传算法。
7.根据权利要求1所述的复混肥主要成分无损定量检测方法,其特征在于,步骤4)和步骤5)之间还包括有校正模型的验证,具体方法如下:
选取已知成分含量的多个预测集样本,将预测集样本的近红外光谱进行预处理后代入所建校正模型中,分别计算预测集样本的总氮、总磷和总钾含量,分别使用相关系数(R)、预测均方根误差(RMSEP)和相对分析误差(RPD)作为标准对所建立的回归模型预测性能进行评价。
CN201610551981.8A 2016-07-13 2016-07-13 基于近红外光谱技术的复混肥主要成分无损定量检测方法 Pending CN106198447A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610551981.8A CN106198447A (zh) 2016-07-13 2016-07-13 基于近红外光谱技术的复混肥主要成分无损定量检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610551981.8A CN106198447A (zh) 2016-07-13 2016-07-13 基于近红外光谱技术的复混肥主要成分无损定量检测方法

Publications (1)

Publication Number Publication Date
CN106198447A true CN106198447A (zh) 2016-12-07

Family

ID=57477822

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610551981.8A Pending CN106198447A (zh) 2016-07-13 2016-07-13 基于近红外光谱技术的复混肥主要成分无损定量检测方法

Country Status (1)

Country Link
CN (1) CN106198447A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108051466A (zh) * 2017-12-17 2018-05-18 中国科学院合肥物质科学研究院 基于x射线荧光光谱分析的化肥成分无损定量检测方法
CN109374556A (zh) * 2018-12-14 2019-02-22 中国科学院合肥物质科学研究院 基于可见-近红外光谱的复合肥中水分含量快速检测方法
CN109696407A (zh) * 2019-01-22 2019-04-30 中国农业大学 一种基于特征波长的椰糠基质有效氮光谱检测方法
CN109708973A (zh) * 2018-12-19 2019-05-03 武汉大学 一种光聚合反应中材料化学-力学参量实时追踪测试***、方法及力学参量优化***、方法
CN112086137A (zh) * 2020-08-18 2020-12-15 山东金璋隆祥智能科技有限责任公司 一种发酵液中山梨糖含量定量分析的方法
CN112945899A (zh) * 2021-01-29 2021-06-11 燕山大学 一种聚谷氨酸复合肥鉴别方法
CN114354534A (zh) * 2021-12-30 2022-04-15 中国航空油料有限责任公司 利用二元线性分类器建立航空煤油性质预测模型的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1712937A (zh) * 2004-06-21 2005-12-28 深圳市无公害农产品质量监督检验站 非破坏快速检测复合肥中多种元素组分含量的方法
CN101762569A (zh) * 2010-01-08 2010-06-30 中国农业大学 一种畜禽粪便工厂化堆肥发酵过程的无损监测方法
CN101769867A (zh) * 2010-01-08 2010-07-07 中国农业大学 一种堆肥产品品质的无损检测方法
CN102353648A (zh) * 2011-10-24 2012-02-15 深圳市芭田生态工程股份有限公司 化肥成分检测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1712937A (zh) * 2004-06-21 2005-12-28 深圳市无公害农产品质量监督检验站 非破坏快速检测复合肥中多种元素组分含量的方法
CN101762569A (zh) * 2010-01-08 2010-06-30 中国农业大学 一种畜禽粪便工厂化堆肥发酵过程的无损监测方法
CN101769867A (zh) * 2010-01-08 2010-07-07 中国农业大学 一种堆肥产品品质的无损检测方法
CN102353648A (zh) * 2011-10-24 2012-02-15 深圳市芭田生态工程股份有限公司 化肥成分检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
任惠芳: "复混肥中氮磷钾养分的近红外定量分析和特征波长的选择", 《中国优秀硕士学位论文全文数据库 农业科技辑》 *
周新奇 等: "近红外光谱结合偏最小二乘法测定复合肥中氮磷钾含量", 《科学仪器服务民生学术大学论文集》 *
王献忠: "近红外光谱法快速测定复混肥总氮含量", 《萍乡高等专科学校学报》 *
陈晓翔 等: "近红外光谱分析技术在肥料成分分析中的应用及前景", 《光谱实验室》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108051466A (zh) * 2017-12-17 2018-05-18 中国科学院合肥物质科学研究院 基于x射线荧光光谱分析的化肥成分无损定量检测方法
CN109374556A (zh) * 2018-12-14 2019-02-22 中国科学院合肥物质科学研究院 基于可见-近红外光谱的复合肥中水分含量快速检测方法
CN109708973A (zh) * 2018-12-19 2019-05-03 武汉大学 一种光聚合反应中材料化学-力学参量实时追踪测试***、方法及力学参量优化***、方法
CN109696407A (zh) * 2019-01-22 2019-04-30 中国农业大学 一种基于特征波长的椰糠基质有效氮光谱检测方法
CN109696407B (zh) * 2019-01-22 2020-11-03 中国农业大学 一种基于特征波长的椰糠基质有效氮光谱检测方法
CN112086137A (zh) * 2020-08-18 2020-12-15 山东金璋隆祥智能科技有限责任公司 一种发酵液中山梨糖含量定量分析的方法
CN112945899A (zh) * 2021-01-29 2021-06-11 燕山大学 一种聚谷氨酸复合肥鉴别方法
CN112945899B (zh) * 2021-01-29 2022-10-14 燕山大学 一种聚谷氨酸复合肥鉴别方法
CN114354534A (zh) * 2021-12-30 2022-04-15 中国航空油料有限责任公司 利用二元线性分类器建立航空煤油性质预测模型的方法

Similar Documents

Publication Publication Date Title
CN106198447A (zh) 基于近红外光谱技术的复混肥主要成分无损定量检测方法
CN101382488B (zh) 利用可见-近红外漫反射光谱技术检测茶鲜叶氮含量的方法
CN102590129B (zh) 近红外检测花生中氨基酸含量的方法
Tan et al. Analysis of different hyperspectral variables for diagnosing leaf nitrogen accumulation in wheat
Huang Chemometric methods in analytical spectroscopy technology
CN101915744A (zh) 物质成分含量的近红外光谱无损检测方法及装置
CN106706553A (zh) 一种快速无损测定玉米单籽粒直链淀粉含量的方法
CN104142311B (zh) 一种用近红外光谱技术预测火炬松松脂产量的方法
CN101762569A (zh) 一种畜禽粪便工厂化堆肥发酵过程的无损监测方法
CN107421911A (zh) 一种基于便携式近红外光谱仪的土壤氮素检测的预处理方法
Lazaar et al. The application of proximal visible and near-infrared spectroscopy to estimate soil organic matter on the Triffa Plain of Morocco
Srivastava et al. Development of hyperspectral model for rapid monitoring of soil organic carbon under precision farming in the Indo-Gangetic Plains of Punjab, India
CN106951720A (zh) 基于典型相关性分析及线性插值的土壤养分模型转移方法
Lu et al. A feasibility quantitative analysis of NIR spectroscopy coupled Si-PLS to predict coco-peat available nitrogen from rapid measurements
CN109738380A (zh) 一种土壤盐渍化程度的高光谱遥感判断方法
Wang et al. Rapid detection of chlorophyll content and distribution in citrus orchards based on low-altitude remote sensing and bio-sensors
CN113049499A (zh) 水体总氮浓度间接遥感反演方法、存储介质和终端设备
CN101614663A (zh) 可见和近红外光谱的蜂花粉品种鉴别方法
CN101769867A (zh) 一种堆肥产品品质的无损检测方法
Cuq et al. Assessing macro-element content in vine leaves and grape berries of vitis vinifera by using near-infrared spectroscopy and chemometrics
CN103760130B (zh) 近红外光谱测定复方麝香注射液中吐温-80含量的方法
Wang et al. Predicting organic matter content, total nitrogen and ph value of lime concretion black soil based on visible and near infrared spectroscopy
CN104502307A (zh) 一种快速检测长牡蛎糖原和蛋白质含量的方法
CN107271389B (zh) 一种基于指标极值的光谱特征变量快速匹配方法
CN105092522A (zh) 一种棉籽粉中铬元素含量的测定方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20161207

RJ01 Rejection of invention patent application after publication