CN101799394B - 一种水力旋流器溢流粒度分布的软测量方法 - Google Patents
一种水力旋流器溢流粒度分布的软测量方法 Download PDFInfo
- Publication number
- CN101799394B CN101799394B CN2010101020466A CN201010102046A CN101799394B CN 101799394 B CN101799394 B CN 101799394B CN 2010101020466 A CN2010101020466 A CN 2010101020466A CN 201010102046 A CN201010102046 A CN 201010102046A CN 101799394 B CN101799394 B CN 101799394B
- Authority
- CN
- China
- Prior art keywords
- hydrocyclone
- overflow
- soft
- overflow granularity
- value
- 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.)
- Expired - Fee Related
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 65
- 239000002245 particle Substances 0.000 title claims abstract description 31
- 238000000691 measurement method Methods 0.000 title claims abstract description 12
- 238000012549 training Methods 0.000 claims abstract description 62
- 238000005259 measurement Methods 0.000 claims abstract description 40
- 238000000034 method Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000000926 separation method Methods 0.000 claims description 11
- 230000007246 mechanism Effects 0.000 claims description 7
- 239000002893 slag Substances 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 3
- 230000001915 proofreading effect Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 14
- 235000019580 granularity Nutrition 0.000 description 83
- 230000006870 function Effects 0.000 description 23
- 238000001514 detection method Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 239000000126 substance Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 239000011362 coarse particle Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012773 Laboratory assay Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000002620 method output Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种水力旋流器溢流粒度分布的软测量方法,其包括以下步骤:1)设置一水力旋流器溢流粒度分布的软测量装置;2)确定先验知识表达式,将先验知识表达式作为先验知识输送给软测量模型训练模块;3)离线读取训练样本,输送到软测量模型训练模块内;4)软测量模型训练模块将先验知识和训练样本集,融合、训练得到水力旋流器溢流粒度分布软测量模型,并将该模型输送给溢流粒度分布测量模块;5)溢流粒度分布测量模块对水力旋流器溢流粒度分布软测量模型进行存储;6)调用溢流粒度分布测量模块中的水力旋流器溢流粒度分布软测量模型,结合实时采集的数据,计算出水力旋流器的溢流粒度分布,并通过监控计算机显示和存储。本发明构思巧妙,精确实用,可广泛用于实际测量过程中。
Description
技术领域
本发明涉及一种软测量方法,特别是关于一种水力旋流器溢流粒度分布的软测量方法。
背景技术
水力旋流器是一种分离混合物的设备,在化工、石油、选矿等领域应用十分广泛。水力旋流器的工作原理是离心沉降,利用离心力进行分级。由于在水力旋流器中,粗颗粒与细颗粒受到的离心力、向心力、流体曳力等大小不同,细颗粒由旋流器的上口溢流管排出,成为溢流产物,粗颗粒从下口底流管排出。在工业生产过程中,水力旋流器溢流产物在某特定粒度范围内的颗粒占颗粒总量的百分比作为反映产物溢流粒度分布的主要指标,称为水力旋流器的溢流粒度分布,记为β+μ(表示粒度大于μ的颗粒占颗粒总量的百分比)或β-μ(表示粒度小于μ的颗粒占颗粒总量的百分比)。实际工业中,通常用溢流产物中74μm以下颗粒占颗粒总量的百分比作为衡量产物粒度分布的重要指标,记为β-74。水力旋流器溢流粒度分布直接决定了后级浮选工艺的生产效率和产品质量,因此实时、准确地测量水力旋流器溢流产物的粒度分布具有重要意义。
水力旋流器溢流粒度分布的检测方法有离线化验和在线检测两种方法。离线化验为目前选矿生产广泛采用的方法,缺点为测量时间间隔过大,无法在线指导生产。
在线检测方法有在线粒度仪检测与软测量两种。在线粒度仪检测虽然可以在线测量,但是价格昂贵,而且测量间隔周期通常为数分钟,不适应闭环控制的需要。研究者们根据水力旋流器溢流粒度分布与水力旋流器工作状态的关系,提出了各种具体计算水力旋流器溢流粒度分布的经验公式。水力旋流器的工作状态包括入口处的给料浓度Cv、给料压力P、给料流量Q等,这些变量容易实现在线测量,因此,可以通过经验公式推算出溢流粒度分布。但是,这些经验公式的参数较难根据实际情况进行调整,灵活性较差,计算精度不能满足实际要求。
此外,还可以利用神经网络方法对水力旋流器的溢流粒度分布进行软测量,即选择合适的辅助变量,结合试验样本得到软测量算法,从而推算出水力旋流器的溢流粒度分布。但现有的神经网络软测量方法在实际运用过程中,都需要对水力旋流器出口处的溢流产物浓度进行实时测量,而在很多实际生产过程中,水力旋流器出口处的溢流浓度很难在线直接测量,这也限制了现有神经网络软测量方法的实际应用。
发明专利“测量水力旋流器固液分离溢流粒度分布的支持向量机方法”(专利号:200510086685.7,已授权)提出了一种溢流粒度分布的测量方法。该方法选择水力旋流器入口处给料的浓度、压力、流量和出口处的溢流浓度作为辅助变量,采用支持向量机方法推算出水力旋流器的溢流粒度分布。但是,该方法在运用过程中并没有解决需要测量溢流浓度的问题。
发明内容
针对上述问题,本发明的目的是提供一种无需测量溢流浓度,可实时在线测量,可靠性高,精确实用的水力旋流器溢流粒度分布的软测量方法。
为实现上述目的,本发明采取以下技术方案:一种水力旋流器溢流粒度分布的软测量方法,其包括以下步骤:1)设置一水力旋流器溢流粒度分布的软测量装置,所述软测量装置包括一渣浆泵,所述渣浆泵通过一管道连接一水力旋流器,所述水力旋流器的上口处设置有一溢流管,下口处设置有一底流管,在所述管道的所述水力旋流器的入口处设置有一浓度检测设备、一压强采集设备和一流量采集设备,所述浓度检测设备、压强采集设备和流量采集设备的输出端连接一数据采集设备,所述数据采集设备的输出端连接一带有显示屏的监控计算机;所述监控计算机内预设有一先验知识计算模块、一软测量模型训练模块和一溢流粒度分布测量模块;2)根据所述水力旋流器的具体型号以及实际工况,选定一表征所述水力旋流器工作状态与溢流粒度分布关系的经验公式,将经验公式表示成一种核函数的线性组合形式,并作为先验知识表达式输入到所述监控计算机中的所述先验知识计算模块内,由所述先验知识计算模块结合指定的参数,确定先验知识表达式,将先验知识表达式作为待融合的先验知识输送给所述软测量模型训练模块;3)离线读取训练样本,并对读取的训练样本集进行剔除野值、归一化处理后,输送到所述软测量模型训练模块;4)所述软测量模型训练模块将步骤2)中得到的待融合的先验知识和步骤3)中获得的训练样本集,融合、训练得到水力旋流器溢流粒度分布软测量模型,并将该模型输送给所述溢流粒度分布测量模块;5)所述溢流粒度分布测量模块对步骤4)中得到的水力旋流器溢流粒度分布软测量模型进行存储;6)调用所述溢流粒度分布测量模块中的水力旋流器溢流粒度分布软测量模型,结合所述数据采集设备实时输入的浓度值、压强值和流量值,计算出所述水力旋流器的溢流粒度分布,并通过所述监控计算机进行显示和存储。
所述步骤2)中,选定经验公式为:y=β-74=f(d50c)=p(x),将该经验公式表示成一种核函数的线性组合形式为:
式中,y为溢流粒度分布;β-74表示溢流产物中74μm以下颗粒占颗粒总量的比例;d50c=g(x)为校正分离粒度,校正分离粒度d50c是指一特定的粒度值,当一组颗粒的粒度均为d50c时,其进入所述水力旋流器后,从所述溢流管排出和从所述底流管排出的比例相同,均为50%,g(x)为解析形式;x=(Cv P Q)T为旋流器入口处给料的浓度值、压力值和流量值;d为一待定参数,其表示线性组合中核函数的总数;i=1、2、…d,其表示线性组合中第i个成分;γi为一待定参数,K(ti,x)为指定的核函数,ti为核函数中一待定基向量;基向量ti和待定参数d、γi均由所述先验知识计算模块进行迭代运算确定。
所述基向量ti和待定参数d、γi在所述先验知识计算模块中的具体求解步骤如下:
②令迭代步数k=0,记q(x)此时的表达函数q0(x)=0,则误差函数为:r0(x)=p(x)-q0(x);
③求tk+1,使误差函数 若|rk(tk+1)|<ε,则返回d=k,q(x)=qk(x),基向量ti的取值集合为Tk, 终止迭代;
④令Tk+1=Tk∪{tk+1},则构造出一(k+1)×(k+1)维矩阵Kk+1,Kk+1的第(i,j)个元素为K(ti,tj),进而构造出一(k+1)维向量p(T(k+1))=[p(t1),…,p(tk+1)]T;
⑤计算向量矩阵γ(k+1)=(Kk+1+μIk+1)-1p(T(k+1)),其中Ik+1为(k+1)维单位矩阵,令 其中γi (k+1)为向量矩阵γ(k+1)的第i个元素;
⑥令rk+1(x)=p(x)-qk+1(x),k=k+1,返回步骤③。
所述步骤4)中溢流粒度分布软测量模型的表达式为:
式中,y为溢流粒度分布软测量结果,αi *为一待定参数,参数d和基向量ti由所述步骤2)确定。
所述待定参数αi *通过求解以下优化公式得到:
式中,KXT为l×d维矩阵,其中第(i,j)个元素为K(ti,xj);KT的第(i,j)个元素为K(ti,tj);训练样本数量为l个,y=[y1,…,yl]T为训练样本输出向量;λ为一指定值;γ=[γ1,…,γi,…,γd]T由所述步骤2)确定,αi *为向量α*中的第i个元素。
所述步骤3)中读取的训练样本为:
xi=(Cv(i)P(i)Q(i))T,yi=β-74(i),i=1,2,…,l
式中,Cv(i)为第i份样品的入口浓度值;P(i)为第i份样品的入口压强值;Q(i)为第i份样品的入口流量值;β-74(i)为采样时刻第i份样品的溢流粒度分布;l为训练样本数量。
本发明由于采取以上技术方案,其具有以下优点:1、本发明设置的溢流粒度分布软测量装置,可实时在线检测水力旋流器入口处给料的浓度值、压强值和流量值,检测的数据通过输入监控计算机,由监控计算机内预制的模块实时运算出水力旋流器出口处的溢流粒度分布,因此,测量间隔时间短、效率高,能更方便地为生产提供依据。2、本发明在监控计算机内预制有先验知识计算模块、软测量模型训练模块和溢流粒度分布测量模块,测量过程中,只需要在线检测水力旋流器入口处给料的浓度值、压强值和流量值,就可以运算出水力旋流器的溢流粒度分布,因此,不需要测量水力旋流器出口处的溢流浓度值,解决了溢流浓度值很难在线测量的问题。3、本发明由于将水力旋器溢流粒度分布计算的经验公式做为先验知识,结合训练样本,融合、训练得到水力旋器溢流粒度分布的测量模型,因此,并不完全依赖于训练样本,可靠性要高于传统的完全依赖样本的学习方法,更适于实际溢流粒度分布的测量应用。本发明构思巧妙,精确实用,可广泛用于实际测量过程中。
附图说明
图1是本发明的软测量装置示意图
图2是本发明的软测量模块示意图
图3是本发明的测量流程示意图
图4是本发明用于实施例的测量结果示意图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明包括以下步骤:
1)如图1所示,设置一水力旋流器溢流粒度分布的软测量装置,该装置包括一渣浆泵1,渣浆泵1通过一管道2连接一水力旋流器3,水力旋流器3的上口处设置有一溢流管31,下口处设置有一底流管32。在管道2的水力旋流器3入口处设置有一浓度检测设备4、一压强检测设备5和一流量检测设备6,浓度检测设备4、压强检测设备5和流量检测设备6的输出端连接一数据采集设备7,数据采集设备7的输出端连接一带有显示屏的监控计算机8。如图2所示,监控计算机8内预设有一先验知识计算模块81、一软测量模型训练模块82和一溢流粒度分布测量模块83。
工作过程中,矿浆通过渣浆泵1泵入管道2中,在矿浆到达水力旋流器3的入口处时,通过浓度检测设备4、压强检测设备5和流量检测设备6分别对水力旋流器3入口处给料的浓度值、压强值和流量值进行测量,测量得到的数据输入到数据采集设备7,由数据采集设备7最后输入到监控计算机8中,进行溢流粒度分布的运算。在水流旋流器3的作用下,矿浆中的细颗粒经水力旋流器3上口的溢流管31排出,成为溢流产物,粗颗粒从水力旋流器3下口的底流管32排出。
2)根据水力旋流器3的具体型号以及实际工况,选定一表征水力旋流器3工作状态与溢流粒度分布关系的经验公式,将经验公式表示成一种核函数的线性组合形式,并作为先验知识表达式输入到监控计算机8中的先验知识计算模块81内,由先验知识计算模块81结合指定的参数,确定先验知识表达式,将先验知识表达式作为待融合的先验知识输送给软测量模型训练模块82。该步骤的具体过程如下:
对于不同型号的水力旋流器3,其入口处给料的浓度、压强、流量与溢流粒度分布关系的经验公式具有不同的表达形式。选定校正分离粒度d50c的经验公式d50c=g(x),然后通过溢流粒度分布公式β-74=f(d50c),计算溢流粒度分布。
其中,校正分离粒度d50c是指一特定的粒度值,当一组颗粒的粒度均为d50c时,其进入水力旋流器3后,从溢流管31排出和从底流管32排出的比例相同,均为50%。对于给定的水力旋流器,d50c取决于旋流器入口处给料的浓度值、压强值和流量值,即可以用g(x)表示,其中x=(Cv P Q)T为旋流器入口处给料的浓度值、压力值和流量值。β-74表示溢流产物中74μm以下颗粒占颗粒总量的比例。A.L.穆拉尔和G.V.杰根森提出了一个表格近似表示f(d50c),如表1所示,实际使用时可以通过插值计算得到溢流粒度分布β-74。
表1溢流粒度分布与校正分离粒度关系表
校正分离粒度d50c(μm) | 溢流粒度分布β-74 |
39.96 | 0.99 |
54.02 | 0.95 |
67.34 | 0.90 |
92.50 | 0.80 |
123.58 | 0.70 |
153.92 | 0.60 |
205.72 | 0.50 |
选定计算d50c的经验公式为g(x),通过插值计算即可确定计算溢流粒度分布的经验公式为:y=β-74=f(d50c)=p(x)。将该经验公式表示成一种核函数的线性组合形式为:
式中,y为溢流粒度分布,x为水力旋流器3入口处给料的浓度值、压强值与流量值;d为一待定参数,其表示线性组合中核函数的总数;i=1、2、…d,其表示线性组合中第i个成分;γi为一待定参数,K(ti,x)为指定的核函数,ti为核函数中一待定基向量;基向量ti和待定参数d、γi均由先验知识计算模块81进行迭代运算确定。
基向量ti和待定参数d、γi的具体求解步骤如下:
②令迭代步数k=0,记q(x)此时的表达函数q0(x)=0,则误差函数为:r0(x)=p(x)-q0(x);
③求tk+1,使误差函数 若|rk(tk+1)|<ε,则返回d=k,q(x)=qk(x),基向量ti的取值集合为Tk, 终止迭代;
④令Tk+1=Tk∪{tk+1},则构造出一(k+1)×(k+1)维矩阵Kk+1,Kk+1的第(i,j)个元素为K(ti,tj),进而构造出一(k+1)维向量p(T(k+1))=[p(t1),…,p(tk+1)]T;
⑤计算向量矩阵γ(k+1)=(Kk+1+μIk+1)-1p(T(k+1)),其中Ik+1为(k+1)维单位矩阵,令 其中γi (k+1)为向量矩阵γ(k+1)的第i个元素;
⑥令rk+1(x)=p(x)-qk+1(x),k=k+1,返回步骤③。
基向量ti和待定参数d、γi确定后,即可确定先验知识表达式(1),将选验知识表达式(1)作为待融合的先验知识输送给软测量模型训练模块82。
3)离线读取训练样本,并对读取的训练样本集进行剔除野值、归一化处理后,输送到软测量模型训练模块82中。该步骤的具体过程如下:
在旋流器3的溢流口31取矿浆样品,同时记录取样时刻的旋流器3入口处给料的浓度值、压力值和流量值,通过实验室化验得到溢流口矿浆样品的溢流粒度分布值β-74,并重复此过程,从而得到l个训练样本,其中第i份样品的入口浓度值记为Cv(i),第i份样品的入口压强值记为P(i),第i份样品的入口流量值记为Q(i),第i份样品的溢流粒度分布记为β-74(i),则Cv(i)、P(i)、Q(i)、β-74(i)构成一组训练样本,记为{xi,yi}:
xi=(Cv(i)P(i)Q(i))T,yi=β-74(i),i=1,2,…,l
将上述得到的训练样本集输入到软测量模型训练模块82中。
4)给定水力旋流器3溢流粒度分布软测量模型的一个核函数的线性组合形式,软测量模型训练模块82将步骤2)中得到的待融合的先验知识表达式(1)和步骤3)中获得的训练样本集{xi,yi},融合、训练得到水力旋流器溢流粒度分布软测量模型中的未知参数,最终确定水力旋流器溢流粒度分布的软测量模型,将该模型输送给溢流粒度分布测量模块83。该步骤的具体过程为:
给定水力旋流器3溢流粒度分布软测量模型的一个核函数的线性组合形式为:
式中,y=β-74为溢流粒度分布软测量结果,αi *一为待定参数,参数d和基向量ti由所述步骤2)确定。
将步骤2)中得到的待融合的先验知识表达式(1)中的参数d、γi和基向量ti与步骤3)中获得的训练样本集{xi,yi},融合、训练得到水力旋流器3溢流粒度分布软测量模型中的未知参数αi *,αi *的求解通过以下优化公式得到:
式中,KXT为l×d维矩阵,其中第(i,j)个元素为K(ti,xj);KT的第(i,j)个元素为K(ti,tj);训练样本数量为l个,y=[y1,…,yl]T为训练样本输出向量;λ为一指定值;γ=[γ1,…,γi,…,γd]T由所述步骤2)确定,αi *为向量α*中的第i个元素。
求得αi *后,水力旋流器溢流粒度分布软测量模型可以确定,将该模型输送给溢流粒度分布测量模块83。
5)溢流粒度分布测量模块83对步骤4)中得到的水力旋流器溢流粒度分布软测量模型进行存储。
6)调用溢流粒度分布测量模块83中的水力旋流器溢流粒度分布软测量模型,结合数据采集设备7实时输入的水力旋流器3入口处的给料浓度值、压强值和流量值,计算出水力旋流器3的溢流粒度分布,并通过监控计算机8进行显示和存储。在实际测量过程中,如果发现溢流粒度分布测量模块83中的溢流粒度分布软测量模型精度不能满足要求时,返回步骤2),再进行融合、训练得到新的溢流粒度分布软测量模型,并用该新模型替换原有的旧模型,继续进行溢流粒度分布的软测量。
如图3所示,在步骤6)中,测量溢流粒度分布的具体步骤如下:
(A)测量开始,决定使用方式:选择进行训练或者测量,若选择测量则转至(H),否则,转至(B)。
(B)选择合适的经验公式作为先验知识。
(C)指定参数:指定核函数K(t,x)、阈值ε、参数μ和λ的值。
(D)读取训练样本{xi,yi}。
(E)对训练样本进行预处理:剔除野值、归一化。
(F)进行训练。
(G)保存训练结果,训练过程结束,转至结尾(M)。
(H)读取训练结果表达式(2)。
(I)读取入口浓度、压强以及流量数据,得到表达式(2)中的输入量x。
(J)代入表达式(2)进行计算,得到溢流粒度分布软测量结果y。
(K)由监控计算机保存结果。
(L)由用户选择是否继续计算,若是,则转至(I),否则转至结尾(M)。
(M)结束。
实施例:该实施例所采用的数据为某选矿厂的实测数据,样本一共有520组,随机选取其中的50组作为训练样本,然后用所有数据检验训练得到结果。每组样本包括三维输入变量xi=(Cv(i)P(i)Q(i))T和输出变量yi=β-74(i)。利用本发明测量步骤如下:
1)选择进行训练;
2)经验公式包括两部分:校正分离粒度d50c的计算,通过校正分离粒度d50c计算溢流粒度分布β-74。实施例中采用如下公式计算d50c:
式中,δ为矿石密度、ρm为矿浆密度,δ可以通过实验室化验确定、ρm计算公式如下:
ρm=1+Cv(δ-1)
根据表1,并通过插值,可以得到经验公式如下。
y=β-74=f(d50c)=p(x)
3)选定核函数为高斯型核函数:
式中,‖·‖表示向量的范数,指定ε=0.05,μ=0.001,λ=0.1;
4)求取先验知识的表达式(1);
5)随机读取50组训练样本{xi,yi};
6)对训练样本进行预处理:判断训练样本是否都在正常范围内,是否存在异常数据需要剔除;
7)求取训练结果,保存参数γi、d、αi *和基向量ti;
8)利用上面的训练结果进行测量。对于某个待测样本,输入xi值,根据表达式(2)计算相应的yi值(读取数据、读取训练结果和显示结果的过程不再详述);
9)按照第8)步所述的方法计算所有样本的输出;
10)结束。
如图4所示,将采用本发明所测得的溢流粒度分布值和实验室化验所得的溢流粒度分布值进行绘图表示,横坐标为实验室化验所得溢流粒度分布真实值,纵坐标为采用本发明测得的溢流粒度分布值,两条直线为±0.02的误差线。利用本发明进行测量的误差主要指标如下:
误差最大值:9.19%
误差最小值:-14.30%
误差均值:0.024786%
误差标准差:1.62%
而直接采用经验公式计算,误差主要指标如下:
误差最大值:13.34%
误差最小值:-14.97%
误差均值:-0.085927%
误差标准差:2.072%
可见,采用本发明进行溢流粒度分布值测量,其测量精度有了很大提高。
上述实施例仅用于说明本发明,其中各部件的结构、设置位置及连接方式都是可以有所变化的,凡是在本发明技术方案的基础上进行的改进和等同变换,均不应排除在本发明的保护范围之外。
Claims (3)
1.一种水力旋流器溢流粒度分布的软测量方法,其包括以下步骤:
1)设置一水力旋流器溢流粒度分布的软测量装置,所述软测量装置包括一渣浆泵,所述渣浆泵通过一管道连接一水力旋流器,所述水力旋流器的上口处设置有一溢流管,下口处设置有一底流管,在所述水力旋流器的入口处设置有一浓度检测设备、一压强采集设备和一流量采集设备,所述浓度检测设备、压强采集设备和流量采集设备的输出端连接一数据采集设备,所述数据采集设备的输出端连接一带有显示屏的监控计算机;所述监控计算机内预设有一先验知识计算模块、一软测量模型训练模块和一溢流粒度分布测量模块;
2)根据所述水力旋流器的具体型号以及实际工况,选定一表征所述水力旋流器工作状态与溢流粒度分布关系的经验公式,将经验公式表示成一种核函数的线性组合形式,并将该经验公式输入到所述监控计算机中的所述先验知识计算模块内,由所述先验知识计算模块结合指定的参数之后形成先验知识表达式,将先验知识表达式作为待融合的先验知识输送给所述软测量模型训练模块;其中,选定经验公式为:y=β-74=f(d50c)=p(x),将该经验公式表示成一种核函数的线性组合形式为:
式中,y为溢流粒度分布;β-74表示溢流产物中74μm以下颗粒占颗粒总量的比例;d50c=g(x)为校正分离粒度,校正分离粒度d50c是指一特定的粒度值,当一组颗粒的粒度均为d50c时,其进入所述水力旋流器后,从所述溢流管排出和从所述底流管排出的比例相同,均为50%,g(x)为解析形式;x=(Cv P Q)T为旋流器入口处给料的浓度值、压力值和流量值;d为一待定参数,其表示线性组合中核函数的总数;i=1、2、…d,其表示线性组合中第i个成分;γi为一待定参数,K(ti,x)为指定的核函数,ti为核函数中一待定基向量;基向量ti和待定参数d、γi均由所述先验知识计算模块进行迭代运算确定;
3)离线读取训练样本,并对读取的训练样本集进行剔除野值、归一化处理后,输送到所述软测量模型训练模块;
4)所述软测量模型训练模块将步骤2)中得到的待融合的先验知识和步骤3)中获得的训练样本集,融合、训练得到水力旋流器溢流粒度分布软测量模型,并 将该模型输送给所述溢流粒度分布测量模块;其中,溢流粒度分布软测量模型的表达式为:
式中,KXT为l×d维矩阵,其中第(i,j)个元素为K(ti,xj);KT的第(i,j)个元素为K(ti,tj);训练样本数量为l个,y=[y1,…,yi]T为训练样本输出向量;λ为一指定值;γ=[γ1,…,γi,…,γd]T由所述步骤2)确定, 为向量α*中的第i个元素;
5)所述溢流粒度分布测量模块对步骤4)中得到的水力旋流器溢流粒度分布软测量模型进行存储;
6)调用所述溢流粒度分布测量模块中的水力旋流器溢流粒度分布软测量模型,结合所述数据采集设备实时输入的浓度值、压强值和流量值,计算出所述水力旋流器的溢流粒度分布,并通过所述监控计算机进行显示和存储。
2.如权利要求1所述的一种水力旋流器溢流粒度分布的软测量方法,其特征在于:所述基向量ti和待定参数d、γi在所述先验知识计算模块中的具体求解步骤如下:
②令迭代步数k=0,记q(x)此时的表达函数q0(x)=0,则误差函数为:r0(x)=p(x)-q0(x);
④令Tk+1=Tk∪{tk+1},则构造出一(k+1)×(k+1)维矩阵Kk+1,Kk+1的第(i,j)个元素为K(ti,tj),进而构造出一(k+1)维向量p(T(k+1))=[p(t1),…,p(tk+1)]T;
⑥令rk+1(x)=p(x)-qk+1(x),k=k+1,返回步骤③。
3.如权利要求1或2所述的一种水力旋流器溢流粒度分布的软测量方法,其特征在于:所述步骤3)中读取的训练样本为:
xi=(Cv(i) P(i) Q(i))T,yi=β-74(i),i=1,2,…,l
式中,Cv(i)为第i份样品的入口浓度值;P(i)为第i份样品的入口压强值;Q(i)为第i份样品的入口流量值;β-74(i)为采样时刻第i份样品的溢流粒度分布;l为训练样本数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101020466A CN101799394B (zh) | 2010-01-27 | 2010-01-27 | 一种水力旋流器溢流粒度分布的软测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101020466A CN101799394B (zh) | 2010-01-27 | 2010-01-27 | 一种水力旋流器溢流粒度分布的软测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101799394A CN101799394A (zh) | 2010-08-11 |
CN101799394B true CN101799394B (zh) | 2011-10-26 |
Family
ID=42595164
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101020466A Expired - Fee Related CN101799394B (zh) | 2010-01-27 | 2010-01-27 | 一种水力旋流器溢流粒度分布的软测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101799394B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102652925A (zh) * | 2012-04-26 | 2012-09-05 | 中冶南方工程技术有限公司 | 高炉喷煤中速磨制粉***煤粉粒度的测量*** |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103331205B (zh) * | 2013-06-04 | 2015-08-19 | 广西华锡集团股份有限公司 | 一种模糊闭路自整定动态水力旋流分级机 |
CN103350037B (zh) * | 2013-07-02 | 2015-06-17 | 太原科技大学 | 一种旋流器溢流粒度控制方法 |
CN104462816A (zh) * | 2014-12-09 | 2015-03-25 | 鞍钢集团矿业公司 | 一种确定水力旋流器分离粒度的方法 |
CN105117553A (zh) * | 2015-09-06 | 2015-12-02 | 长沙有色冶金设计研究院有限公司 | 一种水力旋流器溢流粒度的软测量方法 |
CN109829177A (zh) * | 2018-10-18 | 2019-05-31 | 太原理工大学 | 一种多工况环境下球磨机溢流粒度分布软测量方法 |
CN115270661B (zh) * | 2022-08-05 | 2024-07-12 | 金川集团股份有限公司 | 水力旋流器溢流细度软测量方法、装置、电子设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5248442A (en) * | 1989-10-20 | 1993-09-28 | Mintek | Method and apparatus for measuring shade of hydrocyclone underflow |
CN1598534A (zh) * | 2004-07-28 | 2005-03-23 | 东北大学 | 基于案例推理的磨矿***溢流粒度指标软测量方法 |
CN1769861A (zh) * | 2005-10-21 | 2006-05-10 | 清华大学 | 测量水力旋流器固液分离溢流粒度分布的支持向量机方法 |
-
2010
- 2010-01-27 CN CN2010101020466A patent/CN101799394B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5248442A (en) * | 1989-10-20 | 1993-09-28 | Mintek | Method and apparatus for measuring shade of hydrocyclone underflow |
CN1598534A (zh) * | 2004-07-28 | 2005-03-23 | 东北大学 | 基于案例推理的磨矿***溢流粒度指标软测量方法 |
CN1769861A (zh) * | 2005-10-21 | 2006-05-10 | 清华大学 | 测量水力旋流器固液分离溢流粒度分布的支持向量机方法 |
Non-Patent Citations (2)
Title |
---|
兰岚.水力旋流器控制方法的研究与应用.《矿冶》.2003,第12卷(第4期),19-22. * |
刘啸时等.选矿过程中水力旋流器溢流粒度分布软测量研究.《矿冶》.2005,第14卷(第1期),74-77. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102652925A (zh) * | 2012-04-26 | 2012-09-05 | 中冶南方工程技术有限公司 | 高炉喷煤中速磨制粉***煤粉粒度的测量*** |
Also Published As
Publication number | Publication date |
---|---|
CN101799394A (zh) | 2010-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101799394B (zh) | 一种水力旋流器溢流粒度分布的软测量方法 | |
CN109187534A (zh) | 水质检测方法及其水样图像识别装置 | |
CN101382556A (zh) | 一种基于数据驱动的煤泥浮选精煤灰份软测量方法 | |
CN101556234A (zh) | 气水两相流多参数测量方法和装置 | |
CN101441099A (zh) | 旋流浓集的稀疏气固两相流电容层析成像测量方法及装置 | |
CN1307415C (zh) | 基于案例推理的磨矿***溢流粒度指标软测量方法 | |
CN104459064A (zh) | 一种化学需氧量在线监控电路 | |
CN105301059B (zh) | 测定气液旋流液相停留时间分布的装置和方法 | |
CN104515562A (zh) | 一种多相流微压差测量装置及流量计量方法 | |
RU2344288C2 (ru) | Способ определения продуктивности группы скважин | |
CN103399134B (zh) | 一种基于输出观测器的污水cod软测量方法 | |
CN209027943U (zh) | 单井原油含水率在线称重测量*** | |
CN100570327C (zh) | 测量水力旋流器固液分离溢流粒度分布的支持向量机方法 | |
CN111141332A (zh) | 蒸馏酒摘酒过程的导流装置以及在线测量***和方法 | |
CN106404134A (zh) | 一种测量固液两相混合物中固体物的质量的方法及装置 | |
CN101650359A (zh) | 一种基于知识融合机器学习算法的煤热值软测量方法 | |
CN106198085A (zh) | 集成滤芯测试*** | |
CN110632267A (zh) | 一种对超出在线仪表适用量程的高浓度污水的连续测定***及方法 | |
CN109752507A (zh) | 一种岩心驱替实验用油气水计量装置 | |
CN109490143A (zh) | 酸洗液性能在线检测***及方法 | |
CN107607596B (zh) | 一种基于流动分析的硫酸盐法制浆过程纸浆卡伯值的在线检测方法与装置 | |
CN204988446U (zh) | 水表自动检测装置 | |
CN109932278A (zh) | 一种固定床反应器气相停留时间分布测定***及方法 | |
CN114527247A (zh) | 一种基于物联网的锅炉水质综合监测方法、***及设备 | |
CN206038337U (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111026 Termination date: 20220127 |