CN111398243A - 用于钨浮选精矿品位在线检测的光纤拉曼系统及方法 - Google Patents

用于钨浮选精矿品位在线检测的光纤拉曼系统及方法 Download PDF

Info

Publication number
CN111398243A
CN111398243A CN202010255653.XA CN202010255653A CN111398243A CN 111398243 A CN111398243 A CN 111398243A CN 202010255653 A CN202010255653 A CN 202010255653A CN 111398243 A CN111398243 A CN 111398243A
Authority
CN
China
Prior art keywords
spectrum
tungsten
raman
flotation
raman spectrum
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
CN202010255653.XA
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.)
Central South University
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN202010255653.XA priority Critical patent/CN111398243A/zh
Priority to PCT/CN2020/097212 priority patent/WO2021196416A1/zh
Publication of CN111398243A publication Critical patent/CN111398243A/zh
Pending legal-status Critical Current

Links

Images

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/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/65Raman scattering
    • 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/01Arrangements or apparatus for facilitating the optical investigation

Landscapes

  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
  • Spectrometry And Color Measurement (AREA)

Abstract

本发明公开了一种用于钨浮选精矿品位在线检测的光纤拉曼系统及方法。本发明基于拉曼光谱技术,设计了钨浮选中浮选槽泡沫层拉曼光谱在线测量装置,构建了泡沫层矿物的拉曼光谱数据在线处理系统,从而提取实时采集的钨浮选泡沫层拉曼光谱相关特征并构建精矿品位在线定量计算模型,最终在线获取钨浮选精矿品位结果。本发明所设计的光纤拉曼系统基于不同钨精矿品位下的拉曼光谱特征向量建立钨精矿定量分析模型,在浮选过程中实时获取泡沫层中的钨精矿品位值,能够在线监测钨浮选中重要生产和考核指标,为钨浮选全流程控制提供了重要参数。

Description

用于钨浮选精矿品位在线检测的光纤拉曼系统及方法
技术领域
本发明涉及在线光谱检测技术领域,具体涉及一种用于钨浮选精矿品位在线检测的光纤拉曼系统及方法。
背景技术
矿物浮选是钨原矿应用最广泛的一种选矿方法,通过极其复杂的物理化学反应过程进行矿化起泡。应用矿物浮选技术是为了提高钨原矿品位,达到钨矿还原冶炼的要求。钨浮选中精矿品位是精选工序中有用物质的含量,是钨浮选中重要的工艺指标。在线实时获取钨浮选精矿品位能够指导浮选过程中药剂的合理添加和工况参数的设置,使得浮选过程处于最优生产状态。
钨浮选中产生形态和大小不等的各种气泡来携带矿粒,是一种包含有气、液、固的三相复杂体,其流动性强且容易破裂。钨精矿品位与矿物泡沫状态直接相关,同时也是浮选过程中加药和液位控制的重要指标。现有的钨精矿品位分析方法只能离线获取检测值,即人工舀取少量泡沫并晒干,将矿物粉末应用化学方法化验得到精矿品位数据。这种方法受人为干扰较大、化验过程复杂、成本较高,往往一天只能化验一份样本,难以实时指导钨浮选中的加药和工况参数调整,直接影响最终钨精矿的回收率。少数浮选企业应用X射线荧光光谱技术直接测量浮选槽泡沫的矿物品位,但对待检测泡沫样品的制备要求较高,且由于仪器会产生X射线,从而检测过程中具有较强放射性影响,对检测环境的密封性和使用者的安全操作要求较高,提高了单次检测成本和复杂度。拉曼光谱是一种对样本制备要求较低的非接触式光谱检测技术,获取光谱速度快,且不受水的影响,从而能够直接检测精选过程中泡沫层钨矿品位值,具有实现钨浮选精矿品位在线检测的潜力。
发明内容
本发明的目的在于解决钨浮选中精矿品位难以在线检测的问题,设计了一种用于钨浮选精矿品位在线检测的光纤拉曼光谱系统,同时针对系统给出了其在线检测钨精矿品位的方法。
本发明的技术方案是,
一种用于钨浮选精矿品位在线检测的光纤拉曼系统,包括激光发生装置、光谱装置和工业计算机,所述的激光发生装置产生激光并竖直的照射至钨浮选槽内的浮选泡沫层上以激发泡沫层矿物产生拉曼光谱,所述的光谱装置采集浮选泡沫层被照射后产生的拉曼光谱,所述的工业计算机通信连接至光谱装置并接收采集的拉曼光谱信息。
所述的系统,所述的激光发生装置和光谱装置通过光学装置来发射激光和收集光谱,所述的光学装置包括复合型拉曼光纤探头和光学准直透镜,所述的复合型拉曼光纤探头设置于浮选槽上方并竖直朝向泡沫层,所述的光学准直透镜固定于复合型拉曼光纤探头下,所述的激光发生装置通过复合型拉曼光纤探头朝向浮选泡沫层发射激光,所述的光谱装置通过光学准直透镜和复合型拉曼光纤探头采集拉曼光谱。
一种用于钨浮选精矿品位在线检测的方法,包括以下步骤:
步骤S1:设置采集钨浮选泡沫层拉曼光谱的系统;
步骤S2:在预设时间区间内,设定固定的光谱积分时间,从此时间区间内获得泡沫层样本的多个拉曼光谱采样结果,从而构成拉曼光谱集合;
步骤S3:持续采集一段时间内泡沫层样本的拉曼光谱,基于步骤S2构建多个拉曼光谱集合,同时人工采集与拉曼光谱集合相对应的钨精选泡沫层样本,基于ICP-MS方法获取其精矿品位值;
步骤S4:对各个集合中所包含的拉曼光谱进行预处理,滤除拉曼光谱的随机噪声及扣除其荧光散射背景信号;
步骤S5:对各个集合中预处理后的拉曼光谱求取其平均光谱作为代表性光谱,消除泡沫层矿物分布不均带来的检测误差;
步骤S6:基于含钨矿物的分子结构特征及其拉曼光谱表征方法获取泡沫层中含钨矿物的特征峰峰位信息;
步骤S7:结合步骤S6中所获取的特征峰峰位,基于高斯谱峰拟合方法获取含钨矿物特征谱峰的峰高、峰面积和谱峰区域信息,从而获取含钨矿物的特征谱段;
步骤S8:利用所述泡沫层样本的代表性拉曼光谱作为训练样本集合,提取其每个样本所对应光谱的特征谱段信息,基于偏最小二乘回归方法建立泡沫层拉曼光谱与钨浮选精矿品位间的回归模型;
步骤S9:对钨浮选中在线采集的泡沫层拉曼光谱基于所述步骤S3至所述步骤S7进行预处理,得到滤除随机噪声、扣除荧光散射背景和冗余光谱波长影响的含钨矿物拉曼光谱的特征谱段信息,其后将该特征谱段信息作为所述步骤S8中所建立回归模型的输入,从而获取钨浮选中的实时精矿品位。
所述的方法,所述步骤S3中人工采样的泡沫层样本的精矿品位值基于ICP-MS方法离线测定。
所述的方法,所述步骤S4中对泡沫层矿物拉曼光谱的预处理包括平滑滤波、荧光背景扣除和光谱归一化。
所述的方法,所述步骤S6中含钨矿物的分子结构特征及其拉曼光谱表征方法是通过匹配含钨矿物分子键与拉曼光谱特征峰获取。
所述的方法,所述步骤S7中高斯谱峰拟合方法是基于高斯谱峰的数学模型实现,数学模型公式为:
Figure BDA0002437211100000041
其中,R是拉曼位移,I是一定拉曼位移下的拉曼散射强度,I0为基线强度,当原始拉曼光谱经过基线校正后,I0的值为0,S为单个高斯谱峰的峰面积,w为高斯谱峰的半高宽,Rc为最大光谱散射强度所对应的谱峰处的拉曼位移,π为圆周率,因此,单个高斯峰由S、w和Rc三个参数决定。
所述的方法,所述步骤S7中获取含钨矿物的特征谱段,是提取光谱的特征谱段基于高斯分峰拟合结果中特征谱峰的峰面积和半高宽,来确定其特征谱段范围,即多个特征谱峰所包含的拉曼光谱波长点的范围。
所述的方法,所述步骤S8中基于偏最小二乘回归方法建立泡沫层拉曼光谱与钨浮选精矿品位间的回归模型,是基于偏最小二乘回归方法对钨浮选泡沫层样本的代表性拉曼光谱筛选特征谱段后的数据集合进行分析处理从而建立定量模型,包含以下三个步骤:
第一步,随机选取总样本中的70%作为训练集,而剩余30%样本作为测试集,将其划分成训练样本集和测试样本集两个光谱数据集;
第二步,基于训练样本集,采用偏最小二乘回归法建立钨精矿品位和拉曼光谱的对应回归模型;
第三步,通过测试样本对已建立的定量模型进行验证:将测试样本的拉曼光谱数据代入定量模型中进行预测和拟合分析,通过预测值与实际值之间的误差来判断模型性能,若两者之间的误差超过实际浮选过程中精矿品位检测所要求的误差,则通过增大建模样本数、调整偏最小二乘模型中的主成分个数的方式调整模型参数,直至误差符合要求,从而优化回归模型。
所述的方法,所述步骤S8中第二步,采用偏最小二乘回归法建立钨精矿品位和拉曼光谱的对应回归模型包括如下步骤:
步骤(1):首先对输入的训练集样本的拉曼光谱矩阵进行主成分分解,其后利用光谱主成分得分矩阵和精矿品位向量进行回归,其中对光谱矩阵X进行分解的公式如下:
Figure BDA0002437211100000051
公式中,QT为光谱矩阵X分解得到的得分矩阵,而P则为X分解得到的载荷矩阵,qk和pk分别为矩阵QT和P的第k列列向量,l为上述矩阵的列数,而EX则为应用偏最小二乘法拟合矩阵X所产生的偏差,矩阵符号T为求矩阵的转置;
步骤(2):基于光谱矩阵分解后所得到的得分矩阵QT及样本的精矿品位向量Y进行回归,从而得到模型的关联矩阵B计算如下:
B=(QTQ)-1QTY
步骤(3):利用定量分析模型对测试集样本的精矿品位进行预测,具体步骤为根据测试集样本的光谱矩阵X′以及模型的载荷矩阵P求出其得分矩阵Q′,从而最终得到其精矿品位预测值如下:
y=(X′)TB′
其中,B′=[(Q′)T(Q′)]-1(Q′)TY。
与现有技术相比,本发明的技术效果如下:
1)本发明通过所设计的光纤拉曼光谱系统能够在钨浮选中在线获取不同精矿品位泡沫层样本的拉曼光谱并建立定量校正模型,从而在线实时检测钨精矿品位,有效解决了钨浮选中精矿品位难以在线测定的问题,其测量精度能够达到工业现场要求。
2)本发明具有系统硬件简单、检测精度高、无需样本制备的特点。本发明所设计的光纤拉曼光谱系统通过准直透镜及长距离激光和收集光纤能够在开放环境下在线实时收集钨精选槽泡沫层矿物所产生的拉曼光谱。
3)系统特殊的硬件结构设计能够在一定程度上削弱现场强烈震动、浮选槽酸雾腐蚀和电磁干扰对仪器的影响,保证了光谱获取的准确和可靠性。
4)本发明针对多次测量并平均的方法准确得到各个精矿品位下的代表光谱,减小了矿物分布不均匀对定量预测模型精度的影响,所采用的平滑滤波、荧光背景扣除算法及基于高斯谱峰拟合的特征谱段筛选能够消除随机噪声、荧光散射信号和无信息波长的干扰,提高了光谱信号的代表性。
附图说明
图1是钨浮选精矿品位在线检测光纤拉曼系统的硬件结构示意图;
图2是精矿品位预测模型建立的流程图;
图3是在线检测钨浮选中精选槽泡沫层矿物品位的全流程框图;
图1中,1为功率为500mW,激光波长为785nm的激光发生器,2为微型制冷式拉曼光谱仪,3为激光入射光纤,4为收集光纤,5为复合型光谱光纤探头,6为聚焦透镜,7为钨精选槽泡沫层,8为溢流槽表层,9为钨精选槽,10为钨精选溢流槽,11为拉曼光谱系统与工业计算机连接数据光纤,12为工业计算机,13为光谱探头固定支架外保护套、14为光谱探头固定支架、15为光谱探头固定底座。
具体实施方法
下面结合附图和具体实施过程对本发明进行详细说明,所举实例只用于解释本发明,并非限定本发明范围。
钨浮选精矿品位在线检测光纤拉曼系统的硬件结构如图1所示,主要由功率500mW的785nm激光发生器1,微型制冷式拉曼光谱仪2,激光入射光纤3,收集光纤4,复合型拉曼光谱光纤探头5,光学准直透镜6,拉曼光谱系统与工业计算机连接数据光纤11,工业计算机12,光谱探头固定支架外保护套13、光谱探头固定支架14和光谱探头固定底座15组成。激光发生器1用于激发浮选槽泡沫层样本产生拉曼光谱。微型制冷式拉曼光谱仪2用于收集钨精选泡沫层样本拉曼散射。复合型拉曼光谱光纤探头5能够复用于传导激光和收集光谱。光学准直透镜用于实现激光聚焦距离的提升,从而提升拉曼光谱检测距离。工业计算机12用于在线构建钨精矿品位光谱定量模型。光谱探头固定支架外保护套13、光谱探头固定支架14和光谱探头固定底座15用于组合起来对复合型拉曼光谱光纤探头5进行固定及保护。钨精选槽9是钨粗选产生的矿浆与药物混合反应发生的场所,使得矿物粘合泡沫之上,浮到矿浆面并形成矿化泡沫,其表面是泡沫产生最多的区域,具有良好的隔离性。钨精选溢流槽10是盛放由刮板在浮选槽表面刮出的矿化泡沫的器皿,能够进一步提高矿物品位,为最终钨锭的冶炼提供原材料。光纤拉曼系统的复合型光谱探头连接着光学准直透镜6,由固定支架14使其垂直于钨精选泡沫层表面。500mw的785nm激光发生器持续产生稳定的平行激光由入射光纤3传导至复合型光谱探头5中并通过聚焦透镜进一步聚焦并射入钨精选泡沫层7表面粘附钨矿物颗粒的泡沫层中,最终激发含钨矿物产生拉曼散射并由复合光谱探头5中的收集光纤4传导入微型制冷式拉曼光谱仪2,并最终将光谱仪收集到的拉曼光谱发送至工业计算机12。
在浮选领域,钨精选槽泡沫层中含钨矿物的质量百分比也称为钨浮选精矿品位,本实施例通过以上系统实现钨浮选精矿品位在线检测包含以下步骤:
将光谱探头垂直于钨精选槽上方,探头与浮选槽泡沫层距离为15cm左右。同时,将光纤拉曼光谱系统中的制冷微型拉曼光谱仪放置于钨浮选现场的工控室,同时通过高性能光纤与钨精选槽上方的复合拉曼光谱探头相连接。实时采集钨浮选过程精选槽泡沫层样本和其拉曼光谱若干份,在采集拉曼光谱信号时,本实施例的单个光谱的积分时间为10s,实际实施中,光谱积分时间通过对制冷型拉曼光谱仪的光谱采集时间进行设置而改变。对连续采集到的三个拉曼光谱求取平均值,并将其作为钨精选的表征光谱。拉曼光纤光谱系统所采用的激光器波长为785nm,激光输出功率为100mW至500mW,光谱波长范围从150cm-1~2000cm-1,光谱分辨率大于10cm-1
对于采集的样本进行取样、烘干。基于XRF荧光光谱仪,采用标准XRF定量分析来离线确定各个样本的精矿品位精确值,从而作为钨浮选精矿品位的实时标定值,将其作为初始建模样本库。
将各个样本的拉曼光谱经过光谱仪模数转换并输出到工业计算机中存储,构成钨精选泡沫层样本的拉曼光谱库。
在测量过程中,由于拉曼散射信号较弱、仪器本身存在设计缺陷和采集过程中外界干扰的影响,从而采集到的拉曼光谱会叠加噪声信号,而噪声是一种无用信息,通常对有用信息的提取造成影响。为了消除噪声的影响,采用基于S-G滤波法对拉曼光谱进行平滑处理,S-G平滑滤波是基于多项式最小二乘拟合方法对平滑窗口内数据进行拟合,具体方法如下。
设x为采集到的钨精选泡沫层样本的拉曼光谱信号,z为平滑后的信号,其波数范围均为n,z相对于x的失真程度可以用二者之间的误差平方和表示,如公式(1)所示。
Figure BDA0002437211100000091
z的不平滑度可以用其差分平方和表示,如公式(2)所示。
Figure BDA0002437211100000092
引入惩罚因子,构造函数描述保真度与不平滑度之间的加权和如公式(3)所示。
Q=F+λR=||x-z||2+λ||Dz||2 (3)
公式(3)中,Dz即为矩阵z微分的结果,惩罚因子λ直接调节拟合数据的平滑程度,对函数Q应用最小二乘法,即可得到平滑后的拉曼光谱。
测量得到的拉曼光谱中会存在强荧光散射信号的干扰,造成整个光谱基线发生漂移,严重影响定量建模的准确性。为了校正光谱基线和扣除荧光背景的影响,基于airPLS方法对基线进行拟合,即通过调节惩罚因子和最小二乘迭代,改变拟合基线和原始光谱信号之间误差的平方和权重,从而逐步接近实际的基线形状,进而扣除荧光散射信号的影响。
获取含钨矿物的分子结构特征及其拉曼光谱表征结果,其拉曼光谱如图2所示。可以发现含钨矿物的拉曼光谱特征谱峰大多分布于1200-2000cm-1范围内,因此可以初步筛选出1000-2500cm-1作为含钨矿物的关键光谱信息区域。同时,基于含钨矿物的分子结构特征确定其特征谱峰位置,可知其特征谱峰主要位于1268cm-1,1337cm-1,1434cm-1,1499cm-1和1800cm-1。上述确定其特征谱峰位置是基于含钨矿物的分子结构特征来实现的,同时采用高斯谱峰拟合方法获取其特征谱峰的高斯数学模型,进而得到其特征谱峰宽度。具体来说,含钨矿物光谱特征峰是通过分析其矿物分子结构及化学键,从而确定其分子振动模式,最终得到其光谱特征峰信息。拉曼光谱高斯谱峰拟合方法通过已确定的含钨矿物光谱特征谱峰位置来实现。针对采集到的单个拉曼光谱而言,要拟合的拉曼光谱曲线由多个光谱特征峰组成,可以表示为:
Figure BDA0002437211100000101
其中,i为高斯谱峰的索引值,M为含钨矿物的拉曼特征谱峰个数,Rj为拉曼光谱上第j个采样点处的拉曼位移,
Figure BDA0002437211100000102
为拟合得到的拉曼光谱,Si,ωi和Rci分别为第i个高斯拟合谱峰的峰面积、半宽高和峰位置。
高斯谱峰拟合过程就是对高斯谱峰方程的优化求解过程,通常基于最小化均方误差方法使拟合曲线尽可能逼近原始光谱曲线,从而求解最优的Si和ωi,其寻优过程可以描述如下:
Figure BDA0002437211100000103
s.t.ωi>0,Si>0
基于高斯谱峰拟合方法确定其特征谱峰面积及半高宽,从而精确筛选其特征谱段,并减少无信息拉曼光谱波长点的干扰。
基于多个钨精矿样本预处理后的拉曼光谱及其精矿品位值构建钨精矿品位建模样本库,采用多元统计分析方法对拉曼光谱数据进行分析处理建立并建立钨精矿品位定量分析模型。其精矿品位值基于步骤二获取,所采用的的多元统计分析方法为偏最小二乘回归法。目的是获取钨精矿品位预测模型,从而用于钨浮选精矿品位的在线预测。基于图3所示的方法建立钨精矿品位定量分析模型,具体方法如下:
不同矿物品位拉曼代表光谱中存在无用的噪声信息,影响定量建模的精度,从而偏最小二乘(PLSR)建模方法消除消除噪声数据并得到精确的精矿品位定量预测模型。偏最小二乘(PLSR)建模方法对光谱阵X进行矩阵分解,同时利用精矿品位向量回归得到精矿品位定量模型。
首先,对光谱矩阵X进行分解,如下式所示。
Figure BDA0002437211100000111
分解后的矩阵T为得分矩阵,而P则为X分解得到的载荷矩阵,tk和pk分别为矩阵T和P的第k列列向量,l为上述矩阵的列数,而EX则为应用偏最小二乘法拟合矩阵X所产生的偏差。
基于指示函数法(IND)确定光谱阵X的主成分数,即得分向量阵T所含得分向量个数,如下式所示。
Figure BDA0002437211100000112
逐一计算f=1开始不同f取值下的HIND值,选取HIND极小值时的f即为主成分数,并获得其对应的f个得分向量矩阵T=[t1,t2,…,tf]与矿物品位对应的浓度向量y进行多元线性回归(MLR),即得到其偏最小二乘回归(PLSR)模型如下式所示。
Y=Tb+E
回归系数b的最小二乘解B的计算方法如下式所示。
B=(TTT)-1TTY
所建立的钨精选槽泡沫层矿物品位预测模型可用于钨浮选过程中在线获取精矿品位值,具体步骤如下:
步骤(1):在线采集钨精选槽泡沫层的实时拉曼光谱并完成对光谱的预处理,得到光谱阵X。
步骤(2):对光谱阵X分解得到载荷矩阵P,基于指示函数法确定其主成分数f,并计算其前f个得分向量矩阵T=[t1,t2,…,tf]。
步骤(3):基于已经建立的钨精选槽泡沫层矿物品位预测模型回归系数B及拟合残差E,计算得到实时矿物品位y=TB+E。

Claims (10)

1.一种用于钨浮选精矿品位在线检测的光纤拉曼系统,其特征在于,包括激光发生装置、光谱装置和工业计算机,所述的激光发生装置产生激光并竖直的照射至钨浮选槽内的浮选泡沫层上以激发泡沫层矿物产生拉曼光谱,所述的光谱装置采集浮选泡沫层被照射后产生的拉曼光谱,所述的工业计算机通信连接至光谱装置并接收采集的拉曼光谱信息。
2.根据权利要求1所述的系统,其特征在于,所述的激光发生装置和光谱装置通过光学装置来发射激光和收集光谱,所述的光学装置包括复合型拉曼光纤探头和光学准直透镜,所述的复合型拉曼光纤探头设置于浮选槽上方并竖直朝向泡沫层,所述的光学准直透镜固定于复合型拉曼光纤探头下,所述的激光发生装置通过复合型拉曼光纤探头朝向浮选泡沫层发射激光,所述的光谱装置通过光学准直透镜和复合型拉曼光纤探头采集拉曼光谱。
3.一种用于钨浮选精矿品位在线检测的方法,其特征在于,包括以下步骤:
步骤S1:设置采集钨浮选泡沫层拉曼光谱的系统;
步骤S2:在预设时间区间内,设定固定的光谱积分时间,从此时间区间内获得泡沫层样本的多个拉曼光谱采样结果,从而构成拉曼光谱集合;
步骤S3:持续采集一段时间内泡沫层样本的拉曼光谱,基于步骤S2构建多个拉曼光谱集合,同时人工采集与拉曼光谱集合相对应的钨精选泡沫层样本,基于ICP-MS方法获取其精矿品位值;
步骤S4:对各个集合中所包含的拉曼光谱进行预处理,滤除拉曼光谱的随机噪声及扣除其荧光散射背景信号;
步骤S5:对各个集合中预处理后的拉曼光谱求取其平均光谱作为代表性光谱,消除泡沫层矿物分布不均带来的检测误差;
步骤S6:基于含钨矿物的分子结构特征及其拉曼光谱表征方法获取泡沫层中含钨矿物的特征峰峰位信息;
步骤S7:结合步骤S6中所获取的特征峰峰位,基于高斯谱峰拟合方法获取含钨矿物特征谱峰的峰高、峰面积和谱峰区域信息,从而获取含钨矿物的特征谱段;
步骤S8:利用所述泡沫层样本的代表性拉曼光谱作为训练样本集合,提取其每个样本所对应光谱的特征谱段信息,基于偏最小二乘回归方法建立泡沫层拉曼光谱与钨浮选精矿品位间的回归模型;
步骤S9:对钨浮选中在线采集的泡沫层拉曼光谱基于所述步骤S3至所述步骤S7进行预处理,得到滤除随机噪声、扣除荧光散射背景和冗余光谱波长影响的含钨矿物拉曼光谱的特征谱段信息,其后将该特征谱段信息作为所述步骤S8中所建立回归模型的输入,从而获取钨浮选中的实时精矿品位。
4.根据权利要求3所述的方法,其特征在于,所述步骤S3中人工采样的泡沫层样本的精矿品位值基于ICP-MS方法离线测定。
5.根据权利要求3所述的方法,其特征在于,所述步骤S4中对泡沫层矿物拉曼光谱的预处理包括平滑滤波、荧光背景扣除和光谱归一化。
6.根据权利要求3所述的方法,其特征在于,所述步骤S6中含钨矿物的分子结构特征及其拉曼光谱表征方法是通过匹配含钨矿物分子键与拉曼光谱特征峰获取。
7.根据权利要求3所述的方法,其特征在于,所述步骤S7中高斯谱峰拟合方法是基于高斯谱峰的数学模型实现,数学模型公式为:
Figure FDA0002437211090000031
其中,R是拉曼位移,I是一定拉曼位移下的拉曼散射强度,I0为基线强度,当原始拉曼光谱经过基线校正后,I0的值为0,S为单个高斯谱峰的峰面积,W为高斯谱峰的半高宽,Rc为最大光谱散射强度所对应的谱峰处的拉曼位移,π为圆周率,因此,单个高斯峰由S、W和Rc三个参数决定。
8.根据权利要求3所述的方法,其特征在于,所述步骤S7中获取含钨矿物的特征谱段,是提取光谱的特征谱段基于高斯分峰拟合结果中特征谱峰的峰面积和半高宽,来确定其特征谱段范围,即多个特征谱峰所包含的拉曼光谱波长点的范围。
9.根据权利要求3所述的方法,其特征在于,所述步骤S8中基于偏最小二乘回归方法建立泡沫层拉曼光谱与钨浮选精矿品位间的回归模型,是基于偏最小二乘回归方法对钨浮选泡沫层样本的代表性拉曼光谱筛选特征谱段后的数据集合进行分析处理从而建立定量模型,包含以下三个步骤:
第一步,随机选取总样本中的70%作为训练集,而剩余30%样本作为测试集,将其划分成训练样本集和测试样本集两个光谱数据集;
第二步,基于训练样本集,采用偏最小二乘回归法建立钨精矿品位和拉曼光谱的对应回归模型;
第三步,通过测试样本对已建立的定量模型进行验证:将测试样本的拉曼光谱数据代入定量模型中进行预测和拟合分析,通过预测值与实际值之间的误差来判断模型性能,若两者之间的误差超过实际浮选过程中精矿品位检测所要求的误差,则通过增大建模样本数、调整偏最小二乘模型中的主成分个数的方式调整模型参数,直至误差符合要求,从而优化回归模型。
10.根据权利要求9所述的方法,其特征在于,所述步骤S8中第二步,采用偏最小二乘回归法建立钨精矿品位和拉曼光谱的对应回归模型包括如下步骤:
步骤(1):首先对输入的训练集样本的拉曼光谱矩阵进行主成分分解,其后利用光谱主成分得分矩阵和精矿品位向量进行回归,其中对光谱矩阵X进行分解的公式如下:
Figure FDA0002437211090000041
公式中,QT为光谱矩阵X分解得到的得分矩阵,而P则为X分解得到的载荷矩阵,qk和pk分别为矩阵QT和P的第k列列向量,l为上述矩阵的列数,而EX则为应用偏最小二乘法拟合矩阵X所产生的偏差,矩阵符号T为求矩阵的转置;
步骤(2):基于光谱矩阵分解后所得到的得分矩阵QT及样本的精矿品位向量Y进行回归,从而得到模型的关联矩阵B计算如下:
B=(QTQ)-1QTY
步骤(3):利用定量分析模型对测试集样本的精矿品位进行预测,具体步骤为根据测试集样本的光谱矩阵X以及模型的载荷矩阵P求出其得分矩阵Q′,从而最终得到其精矿品位预测值如下:
y=(X′)TB′
其中,B′=[(Q′)T(Q′)]-1(Q′)TY。
CN202010255653.XA 2020-04-02 2020-04-02 用于钨浮选精矿品位在线检测的光纤拉曼系统及方法 Pending CN111398243A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010255653.XA CN111398243A (zh) 2020-04-02 2020-04-02 用于钨浮选精矿品位在线检测的光纤拉曼系统及方法
PCT/CN2020/097212 WO2021196416A1 (zh) 2020-04-02 2020-06-19 用于钨浮选精矿品位在线检测的光纤拉曼系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010255653.XA CN111398243A (zh) 2020-04-02 2020-04-02 用于钨浮选精矿品位在线检测的光纤拉曼系统及方法

Publications (1)

Publication Number Publication Date
CN111398243A true CN111398243A (zh) 2020-07-10

Family

ID=71433579

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010255653.XA Pending CN111398243A (zh) 2020-04-02 2020-04-02 用于钨浮选精矿品位在线检测的光纤拉曼系统及方法

Country Status (2)

Country Link
CN (1) CN111398243A (zh)
WO (1) WO2021196416A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112161959A (zh) * 2020-09-25 2021-01-01 山东省科学院海洋仪器仪表研究所 一种基于荧光光谱全局重构的叶绿素浓度测量方法
CN112730373A (zh) * 2020-12-03 2021-04-30 北京信息科技大学 一种用于深度学习训练的拉曼光谱数据集分析方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114383985B (zh) * 2021-12-24 2023-09-15 鞍钢集团矿业有限公司 基于多参数特征融合的浮选精矿品位检测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1385691A (zh) * 2001-05-16 2002-12-18 昆明理工大学 一种宝石矿物的拉曼识别法
DE102005031140A1 (de) * 2005-07-04 2007-01-18 Erwin Kayser-Threde Gmbh Verfahren und Vorrichtung zur stofflichen Identifizierung von Oberflächen
JP2011043472A (ja) * 2009-08-24 2011-03-03 Agt Gem Laboratory:Kk 宝石物質の鑑別方法及び宝石物質の鑑別書
CN104655607A (zh) * 2015-01-29 2015-05-27 中南大学 一种浮选泡沫层硫化锑矿物品位检测方法
US9618651B2 (en) * 2012-09-26 2017-04-11 Panalytical Inc. Multi-sensor analysis of complex geologic materials
CN108844941A (zh) * 2018-05-30 2018-11-20 武汉工程大学 一种基于拉曼光谱和pca-hca的不同品位磷矿的鉴别和分类方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102495042B (zh) * 2011-12-09 2013-11-20 湖南大学 一种粉末混合物的拉曼光谱准确定量分析方法
KR20160041622A (ko) * 2014-10-08 2016-04-18 삼성전자주식회사 라만 분광 신호 다중 분석 방법
CN105300954B (zh) * 2015-09-02 2018-04-13 中南大学 一种锑矿浮选非均质泡沫层矿物的拉曼光谱表征方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1385691A (zh) * 2001-05-16 2002-12-18 昆明理工大学 一种宝石矿物的拉曼识别法
DE102005031140A1 (de) * 2005-07-04 2007-01-18 Erwin Kayser-Threde Gmbh Verfahren und Vorrichtung zur stofflichen Identifizierung von Oberflächen
JP2011043472A (ja) * 2009-08-24 2011-03-03 Agt Gem Laboratory:Kk 宝石物質の鑑別方法及び宝石物質の鑑別書
US9618651B2 (en) * 2012-09-26 2017-04-11 Panalytical Inc. Multi-sensor analysis of complex geologic materials
CN104655607A (zh) * 2015-01-29 2015-05-27 中南大学 一种浮选泡沫层硫化锑矿物品位检测方法
CN108844941A (zh) * 2018-05-30 2018-11-20 武汉工程大学 一种基于拉曼光谱和pca-hca的不同品位磷矿的鉴别和分类方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
傅强等: "《现代药物分离与分析技术》", 31 December 2017 *
李原鹰等: "基于拉曼光谱高斯分峰拟合的硫化锑含量检测方法", 《光谱学与光谱分析》 *
杨序纲等: "《拉曼光谱的分析与应用》", 30 November 2008 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112161959A (zh) * 2020-09-25 2021-01-01 山东省科学院海洋仪器仪表研究所 一种基于荧光光谱全局重构的叶绿素浓度测量方法
CN112730373A (zh) * 2020-12-03 2021-04-30 北京信息科技大学 一种用于深度学习训练的拉曼光谱数据集分析方法

Also Published As

Publication number Publication date
WO2021196416A1 (zh) 2021-10-07

Similar Documents

Publication Publication Date Title
CN111398243A (zh) 用于钨浮选精矿品位在线检测的光纤拉曼系统及方法
CN105092540B (zh) 一种食用油中重金属铅含量的快速高精度检测方法
CN104062257B (zh) 一种基于近红外光谱测定溶液中总黄酮含量的方法
CN109085136B (zh) 近红外漫反射光谱测定水泥生料氧化物成分含量的方法
CN102095744B (zh) 一种三能煤炭灰分在线监测装置
CN104655607B (zh) 一种浮选泡沫层硫化锑矿物品位检测方法
CN104897607A (zh) 便携式近红外光谱食品建模与快速检测一体化方法和系统
CN107703097B (zh) 利用近红外光谱仪构建快速预测原油性质的模型的方法
CN102410993B (zh) 基于激光诱导等离子体发射光谱标准化的元素测量方法
CN103884706B (zh) 原浆酒在线检测分级系统
CN108760647A (zh) 一种基于可见/近红外光谱技术的小麦带菌量线检测方法
CN105784672A (zh) 一种基于双树复小波算法的毒品检测仪标准化方法
AU6644896A (en) Biological fluid analysis using distance outlier detection
CN107860743A (zh) 利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用
CN105092436A (zh) 一种沉积物粒度光谱分析方法及装置
CN104977271A (zh) 一种红花醇沉过程有效成分近红外在线检测方法
CN105277531B (zh) 一种基于分档的煤质特性测量方法
CN110567894A (zh) 一种含水溶液pH值快速检测方法和检测装置
CN106442396A (zh) 基于近红外技术的甘蔗渣蔗糖含量快检方法
CN109001182A (zh) 封闭容器中酒精含量的拉曼光谱无损测定方法
CN103134764B (zh) 由透射红外光谱预测原油实沸点蒸馏曲线的方法
Cai et al. Online detection of concentrate grade in the antimony flotation process based on in situ Raman spectroscopy combined with a CNN-GRU hybrid model
CN116718613A (zh) 基于太赫兹谱的环氧玻璃纤维复合材料内部缺陷检测方法
CN104596982A (zh) 近红外漫反射光谱技术测定造纸法再造烟叶果胶的方法
CN105866065B (zh) 一种乌洛托品-醋酸溶液中乌洛托品含量分析方法

Legal Events

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

Application publication date: 20200710