CN102169077A - 湿式磨矿过程溢流粒度指标混合智能软测量方法 - Google Patents
湿式磨矿过程溢流粒度指标混合智能软测量方法 Download PDFInfo
- Publication number
- CN102169077A CN102169077A CN2010106097356A CN201010609735A CN102169077A CN 102169077 A CN102169077 A CN 102169077A CN 2010106097356 A CN2010106097356 A CN 2010106097356A CN 201010609735 A CN201010609735 A CN 201010609735A CN 102169077 A CN102169077 A CN 102169077A
- Authority
- CN
- China
- Prior art keywords
- case
- granularity
- soft
- hydrocyclone
- reasoning
- 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.)
- Granted
Links
Images
Landscapes
- Crushing And Grinding (AREA)
- Feedback Control In General (AREA)
Abstract
一种湿式磨矿过程溢流粒度指标混合智能软测量方法,属于自动化测量技术领域。该装置包括球磨机、水力旋流器、泵池、底流泵、阀门、流量计、浓度计、压力计、磨机新给矿的皮带、给矿水管路、水力旋流器给料管、溢流管、数据采集器和计算机。混合智能软测量方法包括以下步骤:(1)辅助变量的选择,(2)样本数据的取得,(3)基于案例推理的磨矿粒度软测量,(4)基于神经网络的磨矿粒度软测量,(5)基于案例推理的可信度因子求解,(6)基于专家规则推理的磨矿粒度最终求解。本发明的优点:在磨矿过程正常运行期间,能够根据过程的实时数据估计出磨矿粒度值,相对误差小、可信度高,是具有很高实用价值、低成本的粒度计量手段。
Description
技术领域
本发明属于自动化测量技术领域,特别涉及一种湿式磨矿过程溢流粒度指标混合智能软测量方法。
背景技术
在选矿行业中,由球磨机和水力旋流器组成的湿式磨矿过程广泛用于将矿石研磨至选矿工艺要求的粒度范围内,粒度过大或过小均对后续的选别作业都会产生不利影响,因此水力旋流器溢流粒度(也称磨矿粒度、溢流粒度)是衡量磨矿过程运行品质的重要指标。目前,常规的磨矿粒度的检测方法有两种:一是人工取样离线化验,在实验室人工测量;另一种是使用粒度的检测设备——粒度计进行在线测量,前一种方法的不足在于:1.人工操作时人为因素影响大,测量结果的客观性差;2.测量的时间间隔长,测量结果反馈的时间也长,因此得到的信息对操作人员缺乏指导意义;第二种方法不足之处在于虽然能够得到比较客观的测量结果,但粒度计价格昂贵,我国多数选矿厂难以配备,而且容易堵塞,现场维护的工作量很大。另外,粒度计分析的时间一般也得几分钟,因此也难以实现真正意义的粒度在线测量。目前实现粒度在线测量的最有效方法是采用软测量技术。现有的粒度软测量方法主要有基于神经网络(NeuralNetwork,NN)的方法(球磨机磨矿过程溢流粒度指标软测量方法,发明专利,专利号:ZL 03133951.4)和基于案例推理(Case-Based Reasoning,CBR)的方法(基于案例推理的磨矿过程溢流粒度指标软测量方法,发明专利,专利号:ZL 200410021565.4)。采用NN进行粒度软测量算法复杂,模型难于训练和实现,模型校正困难且校正效果不明显。而基于CBR的粒度软测量方法具有实现简单,易于维护,具有较好的自学习能力,特别适合工况缓慢变化的工业对象。但是CBR模型逼近非线性函数的能力有限,而NN却能以任意精度逼近非线性函数,可实现更复杂的输入-输出非线性映射关系。因此以上粒度软测量方法各有优缺点和不同的适用范围。由于磨矿过程具有时变,多工况性,因此以上基于单一技术和模型的粒度软测量方法可信度不高,因而难以满足实际生产和过程控制与优化的要求。
发明内容
为了解决现有湿式磨矿过程溢流粒度测量方法的不足,本发明提供一种基于案例推理和神经网络的湿式磨矿过程溢流粒度指标混合智能软测量方法,通过常规在线测量仪表提供的辅助变量测量参数,给出当前水力旋流器溢流粒度的估计值,为磨矿生产过程的优化操作和优化运行提供关键指标。
本发明方法由硬件平台及测量软件组成,其中硬件平台核心包括球磨机、水力旋流器、泵池、底流泵、阀门、流量计、浓度计、压力计、磨机新给矿的皮带、给矿水管路、水力旋流器给料管、溢流管、数据采集器和计算机。其硬件的联接是球磨机的输入端与磨机新给矿的皮带、给矿水管路以及水力旋流器的沉砂入口相接,球磨机输出端与泵池相接,泵池连接底流泵和泵池加水管,底流泵通过水力旋流器给料管与水力旋流器相接,溢流管连接水力旋转器的出口,数据采集器连接流量计、压力计和密度计,数据采集器通过通信总线连接计算机,阀门安装在给矿水管路和泵池加水管的入口处,最终产品从水力旋流器溢流口排出。详细结构如图1所示。以球磨机和水力旋流器组成的闭路磨矿过程,其测量仪表包括:
一个流量计,用于在线测量水力旋流器给矿矿浆流量LG,安装在水力旋流器给料管上;
一个压力计,用于在线测量水力旋流器给矿矿浆压力PG,安装在水力旋流器的给料管上;
两个密度计,用于在线测量水力旋流器给矿矿浆浓度DG和水力旋流器溢流矿浆浓度DO,分别安装在水力旋流器给料管上和溢流管上。
该磨矿过程同时配置了控制计算机(分布式计算机控制***(DCS)、或可编程逻辑控制器(PLC)、或工业控制计算机(IPC))或数据采集器、用以采集测量仪表的信号。
本发明的软测量软件既可以运行在计算机控制***的监控计算机上,也可以运行于独立的计算机上,该软件通过与控制计算机或数据采集器进行通讯,获得实时的过程数据,并给出磨矿粒度的估计结果。
本发明的湿式磨矿过程溢流粒度混合智能软测量模型结构如图2所示,软测量流程如图3所示,方法包括以下步骤:(1)辅助变量的选择,(2)样本数据的取得,(3)基于案例推理的磨矿粒度软测量,(4)基于神经网络的磨矿粒度软测量,(5)基于案例推理的可信度因子求解,(6)基于专家规则推理的磨矿粒度最终求解。
(1)辅助变量的选择
根据对球磨机-水力旋流器组成的磨矿过程机理进行分析,本发明选择磨矿粒度软测量的辅助变量包括:
水力旋流器给矿矿浆浓度DG;
水力旋流器给矿矿浆压力PG;
水力旋流器给矿矿浆流量LG;
水力旋流器溢流矿浆浓度DO。
(2)样本数据的取得
在设备承受能力之内,在覆盖正常操作范围并稍大于正常操作的范围内,给出一组独立变量(新给矿量OF、磨机入口加水流量WF、泵池加水流量WB)的不同的设定值的组合,形成下面的设定值集合
Ssetp={[OF,i,WF,i,WB,i]|i=1,L,m}
其中m为该集合内的元素数量,每个元素包含一个[OF,i,WF,i,WB,i]的三元组。将该设定值集合的每一个元素依次施加于磨矿过程之上,每加一次元素[OF,i,WF,i,WB,i]之后待磨矿过程进入稳态,在水力旋流器溢流处由人工采集样本,送实验室测量样本的粒度数值S3,i,同时记录采样的时刻Ti,作为本次采样的时间变量,并同时记录水力旋流器给矿浓度DG,i、水力旋流器给矿压力PG,i、水力旋流器给矿流量LG,i以及水力旋流器溢流矿浆浓度DO,i。待本次采样及记录完成后,对磨矿过程施加下一元素。由此,可以得到如下的数据集合,即成为软测量样本数据的训练集
{[DG,i,PG,i,LG,i,DO,i]|i=1,L,m}→{[S3,i]|i=1,L,m}
其中,→符号左边的变量为磨矿粒度软测量模型的输入变量(也就是辅助变量),→符号右边的变量为软测量模型的主导变量即磨矿粒度值。
(3)基于案例推理的磨矿粒度软测量
采用发明人已有成果“基于案例推理的磨矿过程溢流粒度指标软测量方法”进行磨矿粒度软测量,得到磨矿粒度的估计值S1。该方法已获得了国家发明专利,其专利号为:ZL200410021565.4。该方法包括:案例表示、案例库初始案例的取得、案例推理、案例存储与维护等实施步骤。
(4)基于神经网络的磨矿粒度软测量
采用发明人已有成果“球磨机磨矿过程溢流粒度指标软测量方法”进行磨矿粒度软测量,得到磨矿粒度的估计值S2。该方法已获得国家发明专利,其专利号为:ZL 03133951.4。该方法包括:辅助变量的选择、训练数据的取得、神经网络软测量模型的学习和使用等步骤。
(5)基于案例推理的可信度因子求解
采用案例推理技术,根据案例推理粒度软测量的粒度估计值输出S1和神经网络粒度软测量的粒度估计值输出S2,并结合软测量模型的辅助变量信息OC={DG,PG,LG,DO},求解该工况下的粒度估计S1,S2的可信度因子u1,u2。基于案例推理的可信度因子u1,u2的具体求解步骤如下:
(A)案例表示
粒度软测量模型的可信度因子求解***以一定的结构进行组织并以案例的形式存于案例数据库中。案例库中每条案例由案例工况描述和案例解组成,案例工况描述F即为软测量模型的辅助变量信息{DG,PG,LG,DO}以及软测量模型的粒度估计输出S1,S2。案例解J就是该工况下的粒度软测量估计值S1,S2的可信度因子J={u1,u2}。另外,为了便于案例检索与匹配及其它案例操作,在案例库表中再增加时间属性,其中时间为案例获得时间。所以对案例库中的案例进行如下案例表示:
为了运算方便,简记为
(B)案例库初始案例的取得
通过工业试验的方式按照如下数据配对获得样本数据建立初始案例库
{[Ti,DG,i,PG,i,FG,i,DO,i,SC,i,SN,i]|i=1,L,m}→{[u1,i,u2,i]|i=1,L,m}
其中可信度因子定义如下:
显然0≤uj,i≤1。
通常,磨矿分级作业要分期处理多种矿石,每一种矿石的硬度、可磨性是不同的,应该针对不同的矿石的可磨性差异分别建立其对应的案例数据库。
(C)案例推理
本发明可信度因子求解的案例推理采用可信度因子求解方法,其基本流程框图如图4所示。左列为案例检索与匹配及案例重用的流程,右列为案例评价与修正的流程。其详细步骤如下:
(a)初始化
进行所有变量的初始化。
(b)案例推理软测量模型和神经网络软测量模型是否以对粒度进行估计,如果否,等待;如果是,则转至(c)。
(c)是否进行可信度因子求解,如果是,则转至(d),进行案例检索与匹配及案例重用的过程;如果否,则转至(k),进行案例评价与修正的过程。
步骤(d)至(j)为案例检索与匹配及案例重用流程,采用案例检索策略中的最近相邻策略。
(d)选择矿石可磨性
不同的矿石具有不同的可磨性,其对应的案例数据库也是不同的,所以选择矿石可磨性也就是选择案例数据库。
(e)读取当前工况描述
也就是读取案例推理软测量模型和神经网络软测量模型各自的粒度估计输出S1,S2以及当前工况描述参数DG,PG,FG,DO。
(f)案例相似度计算
设磨矿过程当前运行工况为MGK,定义MGK的工况描述为F=(f1,f2,f3,f4,f5,f6),MGK的解为JGK。定义案例库中案例为C1,C2ΛCn,其中案例Ck(k=1,L,n)的工况描述为Fk=(f1,k,f2,k,f3,k,f4,k,f5,k,f6,k),Ck的解为Jk。
那么当前工况描述MGK的描述特征fi(i=1,L,6)与案例Ck(k=1,L,n)的案例描述特征fi,k的相似度函数为:
当前工况描述MGK和案例Ck(k=1,L,n)的相似度函数为:
(g)确定相似度阈值
(h)案例检索与匹配
从案例库中挑选出SIM(MGK,Ck)≥SIMyz的所有案例作为匹配案例并依次按“相似度”、“时间”(案例存储时间)值降序排序。
(i)案例重用
一般情况下案例库中不存在与当前工况描述完全匹配的案例,因而检索出的匹配工况的解并不能直接作为当前工况的解,这就需要对检索得到的相似案例进行重用,具体方法如下:
设匹配案例集为
其中m为匹配案例数,SIMi为匹配案例与当前工况描述的案例相似度数值。则当前工况描述的解为
(j)显示并保存可信度因子求解结果
案例重用后就可以把可信度因子求解结果显示在人机界面上,并把当前工况描述、时间、可信度因子等数据保存到相关历史数据库中,以供案例修正与评价及其它操作使用。
步骤(k)至(T)为案例评价与修正流程。
(k)选择矿石可磨性
(l)读取磨矿粒度化验值和化验取样时间:
记粒度化验值为S3,化验取样时间为TA。
(m)记录检索
在历史数据库中检索“时间”属性值T与化验取样时间TA最接近的数据纪录,记检索到的数据纪录为CT
(n)求取实际的粒度软测量可信度因子
分别根据下式求取工况MGK下案例推理粒度软测量和神经网络粒度软测量的实际可信度因子u1,A,u2,A。
(o)求估计可信度因子和实际可信度因子之差通过下式求取
(p)是否满足精度要求
若max{Δu1,Δu2}≤uΔ(uΔ为精度合格标准,这里uΔ=0.04),则说明案例推理可信度因子求解模型的精度合格,不需要进行案例修正;否则,若max{Δu1,Δu2}>uΔ,则说明精度不合格,需要进行案例修正来改善模型精度。
(q)输入与取样时间对应的工况描述
将FT=(f1 T,f2 T,f3 T,f4 T,f5 T,f6 T)再次输入到可信度因子求解***中,重新进行下述(r)~(t)的操作。
(r)案例相似度计算
与步骤(f)一样
(s)确定相似度阈值:
与步骤(g)一样。
(t)案例检索与匹配
与步骤(h)一样。
(u)案例修正
从匹配案例中挑选出具有最大相似度SIMmax的案例并确定其个数Num。
若SIMmax<0.95,则直接把下述案例Cnew加入到案例库中。
Cnew:{[f1 T,f2 T,f3 T,f4 T,f5 T,f6 T]→[u1,A,u2,A]}
若SIMmax≥0.95,则将Cnew取代原案例库中与FT=(f1 T,f2 T,f3 T,f4 T,f5 T,f6 T)具有最大相似度SIMmax的案例。
(v)案例保存
把修正后的案例库中案例保存起来,以供后续案例推理运算使用。
(w)结束
(6)基于专家规则推理的磨矿粒度最终求解
可信度因子u1,u2求解完毕后,由粒度最终求解***对u1,u2的数值进行分析,根据其相互之间的大小情况,给出最终的粒度软测量输出SE。专家推理***的知识表达方式选用广泛使用的产生式规则,即用如下的形式表示:前提为可信度因子u1,u2,结论就是所求解的最终的粒度软测量输出SE,具体求解SE的推理规则如下所示,共6条规则,即R1,R2,R3,R4,R5,R6:
本发明的优点:在磨矿过程正常运行期间,能够根据过程的实时数据估计出磨矿粒度值,相对误差小、可信度高,是具有很高实用价值、低成本的粒度计量手段。
附图说明
图1磨矿过程的流程、测量仪表及计算机配置图;
图2本发明的湿式磨矿过程溢流粒度混合智能软测量模型结构图;
图3本发明湿式磨矿过程溢流粒度指标混合智能软测量方法流程框图;
图4本发明的基于案例推理的可信度因子求解流程框图;
图中:1球磨机,2水力旋流器,3泵池,4底流泵,5阀门,6浓度计,7压力计,8流量计,9数据采集器,10计算机。
具体实施方式
以某赤铁矿选矿厂的磨矿分级工段为实施案例,该选矿生产磨矿过程的具体工艺如图1所示,其中,实线箭头表示物流(原矿,水和矿浆),虚线表示信号流,
水力旋流器给矿压力——PG
水力旋流器给矿流量——LG
水力旋流器给矿浓度——DG
水力旋流器溢流浓度——DO
水力旋流器溢流浓度——DO
软测量模型的辅助变量信息集——OC
案例推理软测量模型的粒度估计——S1
神经网络软测量模型的粒度估计——S2
磨矿粒度的实际化验值——S3
磨矿粒度的(最终)软测量值——SE
案例推理软测量模型粒度估计的可信度因子——u1
神经网络软测量模型粒度估计的可信度因子——u2
浓度计——DT
流量计——FT
压力计——PT
原矿石和一定比例的水送入磨机进行研磨,研磨后的矿浆排入泵池,同时在泵池入口补加一定量的水。泵池内的矿浆由底流泵打入水力旋流器进行分级,细粒级矿浆从水力旋流器溢流口排出进入选别工序,粗粒级矿浆返回磨机再磨。
该装置包括球磨机、水力旋流器、泵池、底流泵、阀门、流量计、浓度计、压力计、磨机新给矿的皮带、给矿水管路、水力旋流器给料管、溢流管、数据采集器和计算机;
球磨机型号为Φ3200mm×3500mm,有效容积25.3m3,筒体转速18.5r/min,最大装球量54吨。
水力旋流器为聚胺脂水力旋流器,筒体内径φ350mm,给矿口尺寸80×60mm,溢流管直径φ100mm,溢流管长度295mm,沉砂嘴直径φ45~55mm,锥角20度,入口压力0.05~0.31Mpa。
首先,按照本说明书的要求安装如下的测量仪表,包括:
核子密度计测量水力旋流器给矿矿浆浓度DG和水力旋流器溢流矿浆浓度OC;
电磁流量计测量水力旋流器给矿流量LG;
电容式压力变送器测量水力旋流器给矿压力PG。
该磨矿过程的正常工作范围为:
球磨机处理量——65±10吨/(台·小时)
磨机内矿浆浓度——70%~80%
水力旋流器溢流浓度——35%~40%
水力旋流器溢流粒度——72%~82%(-200目含量)
水力旋流器给矿压力——0.13MPa~0.17MPa
水力旋流器给矿浓度——40%~55%
介质填充率——38%~42%
球磨机的输入端与磨机新给矿的皮带、给矿水管路以及水力旋流器的沉砂入口相接,球磨机输出端与泵池相接,泵池连接底流泵和泵池加水管,底流泵通过水力旋流器给料管与水力旋流器相接,溢流管连接水力旋转器的出口,一个流量计和一个压力计均安装在水力旋流器给料管上;两个密度计分别安装在水力旋流器给料管上和溢流管上,数据采集器连接流量计、压力计和密度计,数据采集器通过通信总线连接计算机,阀门安装在给矿水管路和泵池加水管的入口处。
以美国Rockwell的Controllogix5K实现基础控制回路的自动控制,并读取上述的过程数据。在监控计算机中以Rockwell的RSView32软件实现监控人机界面。
本发明的湿式磨矿过程溢流粒度混合智能软测量模型结构如图2所示,软测量程序用RSView32提供的VBA应用软件编制。软测量软件在单独的优化计算机上运行,该计算机上装有RSLinx通讯程序负责与Controllogix5K***和上位监控机进行数据通讯,RSLinx与软测量程序之间通过DDE方式进行双向通讯。粒度软测量的案例库和相关数据库采用MicrosoftAccess数据库实现,相关复杂计算通过VBA调用Matlab程序实现。
按照本说明书所述的实现方法选择辅助变量和样本选择,然后采用发明人已有成果,即基于案例推理的磨矿过程溢流粒度指标软测量方法(发明专利,专利号:ZL 200410021565.4)和球磨机磨矿过程溢流粒度指标软测量方法(发明专利,专利号:ZL 03133951.4),分别对基于案例推理的磨矿粒度软测量模型和基于神经网络的磨矿粒度软测量模型进行软件实现。
然后采用本说明书所述方法对可信度因子求解模型进行软件实现。采用工业实验的方法建立可信度因子求解模型的案例库,共86组案例。根据专家经验,案例工况描述特征的权值确定为:
水力旋流器给矿浓度特征权值——ω1=0.15
水力旋流器给矿压力特征权值——ω2=0.05
水力旋流器给矿流量特征权值——ω3=0.05
水力旋流器溢流浓度特征权值——ω4=0.15
案例推理模型的粒度估计输出特征权值——ω5=0.3
神经网络模型的粒度估计输出特征权值——ω6=0.3进行案例修正时,为可信度因子的精度合格标准uΔ确定为0.04,即uΔ=0.04。
最后采用本说明书所述方法对基于规则推理的磨矿粒度最终求解模型进行软件实现。
本实施例在下位DCS***中设置有定时器,每隔10分钟触发一次上位RSView32中的相关标签,从而调用一次磨矿粒度软测量程序,进行一次该时间工况下的粒度软测量,并及时把磨矿粒度软测量结果输出到优化控制计算机的***管理画面,为操作者或相关技术工人进行参数优化和监督操作提供一个很好的参考、指导作用,同时把软测量结果保存到相关历史数据库中。
湿式磨矿过程溢流粒度指标混合智能软测量装置的测量方法,按如下步骤进行:如图3所示,
步骤一、辅助变量的选择
根据对球磨机-水力旋流器组成的磨矿过程机理进行分析,本发明所选择磨矿粒度软测量的辅助变量包括:
水力旋流器给矿矿浆浓度DG;
水力旋流器给矿矿浆压力PG;
水力旋流器给矿矿浆流量LG;
水力旋流器溢流矿浆浓度DO;
步骤二、样本数据的取得
设定独立变量:新给矿量OF、磨机入口加水流量WF、泵池加水流量WB的设定值组合,形成下面的设定值集合
Ssetp={[OF,i,WF,i,WB,i]|i=1,L,m}
其中m为该集合内的元素数量,每个元素包含一个[OF,i,WF,i,WB,i]的三元组,将该设定值集合的每一个元素依次施加于磨矿过程之上,每加一次元素[OF,i,WF,i,WB,i]之后待磨矿过程进入稳态,在水力旋流器溢流处由人工采集样本,送实验室测量样本的粒度数值S3,i,同时记录采样的时刻Ti,作为本次采样的时间变量,并同时记录水力旋流器给矿浓度DG,i、水力旋流器给矿压力PG,i、水力旋流器给矿流量LG,i以及水力旋流器溢流矿浆浓度DO,i,待本次采样及记录完成后,对磨矿过程施加下一元素,由此,得到如下的数据集合,即成为软测量样本数据的训练集
{[DG,i,PG,i,LG,i,DO,i]|i=1,L,m}→{[S3,i]|i=1,L,m}
其中,→符号左边的变量为磨矿粒度软测量模型的输入变量也就是辅助变量,→符号右边的变量为软测量模型的主导变量即磨矿粒度值;
步骤三、基于案例推理的磨矿粒度软测量
采用基于案例推理的磨矿过程溢流粒度指标软测量方法进行磨矿粒度软测量,得到磨矿粒度的估计值S1,该方法包括:案例表示、案例库初始案例的取得、案例推理、案例存储与维护;
步骤四、基于神经网络的磨矿粒度软测量
采用球磨机磨矿过程溢流粒度指标软测量方法进行磨矿粒度软测量,得到磨矿粒度的估计值S2,该方法包括:辅助变量的选择、训练数据的取得、神经网络软测量模型的学习和使用;
步骤五、基于案例推理的可信度因子求解
采用案例推理技术,根据案例推理粒度软测量的粒度估计值输出S1和神经网络粒度软测量的粒度估计值输出S2,并结合软测量模型的辅助变量信息OC={DG,PG,LG,DO},求解该工况下的粒度估计S1,S2的可信度因子u1,u2;
步骤六、基于专家规则推理的磨矿粒度最终求解
可信度因子求解完毕后,由粒度最终求解***对u1,u2的数值进行分析,根据u1,u2相互之间的大小情况,给出最终的粒度软测量输出SE,专家推理***的知识表达方式选用广泛使用的产生式规则,即用如下的形式表示:前提为可信度因子u1,u2,结论就是所求解的最终的粒度软测量输出SE,具体求解SE的推理规则如下所示,共6条规则,即R1,R2,R3,R4,R5,R6:
所述的步骤五中基于案例推理的可信度因子u1,u2的具体求解步骤如下:
(A)案例表示
粒度软测量模型的可信度因子求解***以一定的结构进行组织并以案例的形式存于案例数据库中,案例库中每条案例由案例工况描述和案例解组成,案例工况描述F即为软测量模型的辅助变量信息{DG,PG,LG,DO}以及软测量模型的粒度估计输出S1,S2;案例解J就是该工况下的粒度软测量估计值S1,S2的可信度因子J={u1,u2},在案例库表中再增加时间属性,其中时间为案例获得时间,所以对案例库中的案例进行如下案例表示:
简记为
(B)案例库初始案例的取得
通过工业试验的方式按照如下数据配对获得样本数据建立初始案例库
{[Ti,DG,i,PG,i,FG,i,DO,i,SC,i,SN,i]|i=1,L,m}→{[u1,i,u2,i]|i=1,L,m}
其中可信度因子定义如下:
针对矿石的可磨性差异分别建立案例数据库;
(C)案例推理
可信度因子求解的案例推理采用可信度因子求解方法实现。
步骤C所述的可信度因子求解方法,按如下步骤进行:
(a)初始化
进行所有变量的初始化;
(b)案例推理软测量模型和神经网络软测量模型是否以对粒度进行估计,如果否,等待;如果是,则转至(c);
(c)是否进行可信度因子求解,如果是,则转至(d),进行案例检索与匹配及案例重用的过程;如果否,则转至(k),进行案例评价与修正的过程;
步骤(d)至(j)为案例的检索与匹配及案例重用流程,采用案例检索策略中的最近相邻策略;
(d)选择矿石可磨性
按照不同的矿石的可磨性,选择矿石可磨性的案例数据库;
(e)读取当前工况描述
读取案例推理软测量模型和神经网络软测量模型各自的粒度估计输出S1,S2以及当前工况描述参数DG,PG,FG,DO;
(f)案例相似度计算
设磨矿过程当前运行工况为MGK,定义MGK的工况描述为F=(f1,f2,f3,f4,f5,f6),MGK的解为JGK,定义案例库中案例为C1,C2ΛCn,其中案例Ck(k=1,L,n)的工况描述为Fk=(f1,k,f2,k,f3,k,f4,k,f5,k,f6,k),Ck的解为Jk;
那么当前工况描述MGK的描述特征fi(i=1,L,6)与案例Ck(k=1,L,n)的案例描述特征fi,k的相似度函数为:
当前工况描述MGK和案例Ck(k=1,L,n)的相似度函数为:
(g)确定相似度阈值
(h)案例检索与匹配
从案例库中挑选出SIM(MGK,Ck)≥SIMyz的所有案例作为匹配案例并依次按“相似度”、“时间”(案例存储时间)值降序排序;
(i)案例重用
对检索得到的相似案例进行重用,具体方法如下:
设匹配案例集为
其中m为匹配案例数,SIMi为匹配案例与当前工况描述的案例相似度数值,则当前工况描述的解为
(j)显示并保存可信度因子求解结果
案例重用后把可信度因子求解结果显示在人机界面上,并把当前工况描述、时间、可信度因子等数据保存到相关历史数据库中,以供案例修正与评价及其它操作使用;
步骤(k)至(T)为案例评价与修正流程;
(k)选择矿石可磨性
(l)读取磨矿粒度化验值和化验取样时间:
记粒度化验值为S3,化验取样时间为TA;
(m)记录检索
在历史数据库中检索“时间”属性值T与化验取样时间TA最接近的数据纪录,记检索到的数据纪录为CT
(n)求取实际的粒度软测量可信度因子
分别根据下式求取工况MGK下案例推理粒度软测量和神经网络粒度软测量的实际可信度因子u1A,u2A;
(o)求估计可信度因子和实际可信度因子之差通过下式求取
(p)是否满足精度要求
若max{Δu1,Δu2}≤uΔ(uΔ为精度合格标准,这里uΔ=0.04),则说明案例推理可信度因子求解模型的精度合格,不需要进行案例修正;否则,若max{Δu1,Δu2}>uΔ,则说明精度不合格,需要进行案例修正来改善模型精度;
(q)输入与取样时间对应的工况描述
将FT=(f1 T,f2 T,f3 T,f4 T,f5 T,f6 T)再次输入到可信度因子求解模型中,重新进行下述(r)~(t)的操作;
(r)案例相似度计算
与步骤(f)一样;
(s)确定相似度阈值:
与步骤(g)一样;
(t)案例检索与匹配
与步骤(h)一样;
(u)案例修正
从匹配案例中挑选出具有最大相似度SIMmax的案例并确定其个数Num;
若SIMmax<0.95,则直接把下述案例Cnew加入到案例库中;
Cnew:{[f1 T,f2 T,f3 T,f4 T,f5 T,f6 T]→[u1,A,u2,A]}
若SIMmax≥0.95,则将Cnew取代原案例库中与FT=(f1 T,f2 T,f3 T,f4 T,f5 T,f6 T)具有最大相似度SIMmax的案例;
(v)案例保存
把修正后的案例库中案例保存起来,以供以后案例推理运算使用;
(w)结束。
本发明即湿式磨矿过程溢流粒度软测量方法在磨矿过程正常运行期间,能够根据过程的实时数据估计出磨矿粒度值,相对误差不超过2.5%,成为一个具有很高实用价值、低成本的粒度计量手段。
Claims (4)
1.一种湿式磨矿过程溢流粒度指标混合智能软测量装置,其特征在于:该装置包括球磨机、水力旋流器、泵池、底流泵、阀门、流量计、浓度计、压力计、磨机新给矿的皮带、给矿水管路、水力旋流器给料管、溢流管、数据采集器和计算机;球磨机的输入端与磨机新给矿的皮带、给矿水管路以及水力旋流器的沉砂入口相接,球磨机输出端与泵池相接,泵池连接底流泵和泵池加水管,底流泵通过水力旋流器给料管与水力旋流器相接,溢流管连接水力旋转器的出口,一个流量计和一个压力计均安装在水力旋流器给料管上;两个密度计分别安装在水力旋流器给料管上和溢流管上,数据采集器连接流量计、压力计和密度计,数据采集器通过通信总线连接计算机,阀门安装在给矿水管路和泵池加水管的入口处。
2.采用权利要求1所述的湿式磨矿过程溢流粒度指标混合智能软测量装置的测量方法,其特征在于:按如下步骤进行:
步骤一、辅助变量的选择
水力旋流器给矿矿浆浓度DG,
水力旋流器给矿矿浆压力PG,
水力旋流器给矿矿浆流量LG,
水力旋流器溢流矿浆浓度DO,
步骤二、样本数据的取得
设定独立变量:新给矿量OF、磨机入口加水流量WF、泵池加水流量WB的设定值组合,形成下面的设定值集合
Ssetp={[OF,i,WF,i,WB,i]|i=1,L,m}
其中m为该集合内的元素数量,每个元素包含一个[OF,i,WF,i,WB,i]的三元组,将该设定值集合的每一个元素依次施加于磨矿过程之上,每加一次元素[OF,i,WF,i,WB,i]之后待磨矿过程进入稳态,在水力旋流器溢流处由人工采集样本,送实验室测量样本的粒度数值S3,i,同时记录采样的时刻Ti,作为本次采样的时间变量,并同时记录水力旋流器给矿浓度DG,i、水力旋流器给矿压力PG,i、水力旋流器给矿流量LG,i以及水力旋流器溢流矿浆浓度DO,i,待本次采样及记录完成后,对磨矿过程施加下一元素,由此,得到如下的数据集合,即成为软测量样本数据的训练集
{[DG,i,PG,i,LG,i,DO,i]|i=1,L,m}→{[S3,i]|i=1,L,m}
其中,→符号左边的变量为磨矿粒度软测量模型的输入变量也就是辅助变量,→符号右边的变量为软测量模型的主导变量即磨矿粒度值;
步骤三、基于案例推理的磨矿粒度软测量
采用基于案例推理的磨矿过程溢流粒度指标软测量方法进行磨矿粒度软测量,得到磨矿粒度的估计值S1,该方法包括:案例表示、案例库初始案例的取得、案例推理、案例存储与维护;
步骤四、基于神经网络的磨矿粒度软测量
采用球磨机磨矿过程溢流粒度指标软测量方法进行磨矿粒度软测量,得到磨矿粒度的估计值S2,该方法包括:辅助变量的选择、训练数据的取得、神经网络软测量模型的学习和使用;
步骤五、基于案例推理的可信度因子求解
采用案例推理技术,根据案例推理粒度软测量的粒度估计值输出S1和神经网络粒度软测量的粒度估计值输出S2,并结合软测量模型的辅助变量信息OC={DG,PG,LG,DO},求解该工况下的粒度估计S1,S2的可信度因子u1,u2;
步骤六、基于专家规则推理的磨矿粒度最终求解
可信度因子求解完毕后,由粒度最终求解***对u1,u2的数值进行分析,根据u1,u2相互之间的大小情况,给出最终的粒度软测量输出SE,专家推理***的知识表达方式选用广泛使用的产生式规则,即用如下的形式表示:前提为可信度因子u1,u2,结论就是所求解的最终的粒度软测量输出SE,具体求解SE的推理规则如下所示,共6条规则,即R1,R2,R3,R4,R5,R6:
3.根据权利要求2所述的湿式磨矿过程溢流粒度指标混合智能软测量装置的测量方法,其特征在于:所述的步骤五中基于案例推理的可信度因子u1,u2的具体求解步骤如下:
(A)案例表示
粒度软测量模型的可信度因子求解***以一定的结构进行组织并以案例的形式存于案例数据库中,案例库中每条案例由案例工况描述和案例解组成,案例工况描述F即为软测量模型的辅助变量信息{DG,PG,LG,DO}以及软测量模型的粒度估计输出S1,S2;案例解J就是该工况下的粒度软测量估计值S1,S2的可信度因子J={u1,u2},在案例库表中再增加时间属性,其中时间为案例获得时间,所以对案例库中的案例进行如下案例表示:
简记为
(B)案例库初始案例的取得
通过工业试验的方式按照如下数据配对获得样本数据建立初始案例库
{[Ti,DG,i,PG,i,FG,i,DO,i,SC,i,SN,i]|i=1,L,m}→{[u1,i,u2,i]|i=1,L,m}
其中可信度因子定义如下:
针对矿石的可磨性差异分别建立案例数据库;
(C)案例推理
可信度因子求解的案例推理采用可信度因子求解方法实现。
4.根据权利要求3所述的湿式磨矿过程溢流粒度指标混合智能软测量装置的测量方法,其特征在于:步骤C所述的可信度因子求解方法,按如下步骤进行:
(a)初始化
进行所有变量的初始化;
(b)案例推理软测量模型和神经网络软测量模型是否以对粒度进行估计,如果否,等待;如果是,则转至(c);
(c)是否进行可信度因子求解,如果是,则转至(d),进行案例检索与匹配及案例重用的过程;如果否,则转至(k),进行案例评价与修正的过程;
步骤(d)至(j)为案例的检索与匹配及案例重用流程,采用案例检索策略中的最近相邻策略;
(d)选择矿石可磨性
按照不同的矿石的可磨性,选择矿石可磨性的案例数据库;
(e)读取当前工况描述
读取案例推理软测量模型和神经网络软测量模型各自的粒度估计输出S1,S2以及当前工况描述参数DG,PG,FG,DO;
(f)案例相似度计算
设磨矿过程当前运行工况为MGK,定义MGK的工况描述为F=(f1,f2,f3,f4,f5,f6),MGK的解为JGK,定义案例库中案例为C1,C2ΛCn,其中案例Ck(k=1,L,n)的工况描述为Fk=(f1,k,f2,k,f3,k,f4,k,f5,k,f6,k),Ck的解为Jk;
那么当前工况描述MGK的描述特征fi(i=1,L,6)与案例Ck(k=1,L,n)的案例描述特征fi,k的相似度函数为:
当前工况描述MGK和案例Ck(k=1,L,n)的相似度函数为:
(g)确定相似度阈值
设SIMmax为所有上述求得的相似度的最大值,即:阈值SIMyz可通过下式确定:
(h)案例检索与匹配
从案例库中挑选出SIM(MGK,Ck)≥SIMyz的所有案例作为匹配案例并依次按“相似度”、“时间”(案例存储时间)值降序排序;
(i)案例重用
对检索得到的相似案例进行重用,具体方法如下:
设匹配案例集为
其中m为匹配案例数,SIMi为匹配案例与当前工况描述的案例相似度数值,则当前工况描述的解为
(j)显示并保存可信度因子求解结果
案例重用后把可信度因子求解结果显示在人机界面上,并把当前工况描述、时间、可信度因子等数据保存到相关历史数据库中,以供案例修正与评价及其它操作使用;
步骤(k)至(T)为案例评价与修正流程;
(k)选择矿石可磨性
(l)读取磨矿粒度化验值和化验取样时间:
记粒度化验值为S3,化验取样时间为TA;
(m)记录检索
在历史数据库中检索“时间”属性值T与化验取样时间TA最接近的数据纪录,记检索到的数据纪录为CT
(n)求取实际的粒度软测量可信度因子
分别根据下式求取工况MGK下案例推理粒度软测量和神经网络粒度软测量的实际可信度因子u1A,u2A;
(o)求估计可信度因子和实际可信度因子之差
通过下式求取
(p)是否满足精度要
若max{Δu1,Δu2}≤uΔ,uΔ为精度合格标准,,则说明案例推理可信度因子求解模型的精度合格,不需要进行案例修正;否则,若max{Δu1,Δu2}>uΔ,则说明精度不合格,需要进行案例修正来改善模型精度;
(q)输入与取样时间对应的工况描述
将FT=(f1 T,f2 T,f3 T,f4 T,f5 T,f6 T)再次输入到可信度因子求解模型中,重新进行下述(r)~(t)的操作;
(r)案例相似度计算
与步骤(f)一样;
(s)确定相似度阈值:
与步骤(g)一样;
(t)案例检索与匹配
与步骤(h)一样;
(u)案例修正
从匹配案例中挑选出具有最大相似度SIMmax的案例并确定其个数Num;
若SIMmax<0.95,则直接把下述案例Cnew加入到案例库中;
Cnew:{[f1 T,f2 T,f3 T,f4 T,f5 T,f6 T]→[u1,A,u2,A]}
若SIMmax≥0.95,则将Cnew取代原案例库中与FT=(f1 T,f2 T,f3 T,f4 T,f5 T,f6 T)具有最大相似度SIMmax的案例;
(v)案例保存
把修正后的案例库中案例保存起来,以供以后案例推理运算使用;
(w)结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010609735 CN102169077B (zh) | 2010-12-28 | 2010-12-28 | 湿式磨矿过程溢流粒度指标混合智能软测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010609735 CN102169077B (zh) | 2010-12-28 | 2010-12-28 | 湿式磨矿过程溢流粒度指标混合智能软测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102169077A true CN102169077A (zh) | 2011-08-31 |
CN102169077B CN102169077B (zh) | 2013-04-17 |
Family
ID=44490318
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010609735 Active CN102169077B (zh) | 2010-12-28 | 2010-12-28 | 湿式磨矿过程溢流粒度指标混合智能软测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102169077B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103331205A (zh) * | 2013-06-04 | 2013-10-02 | 广西华锡集团股份有限公司 | 一种模糊闭路自整定动态水力旋流分级机 |
CN103969412A (zh) * | 2014-04-13 | 2014-08-06 | 北京工业大学 | 一种基于群决策案例推理的溶解氧浓度软测量方法 |
CN105117553A (zh) * | 2015-09-06 | 2015-12-02 | 长沙有色冶金设计研究院有限公司 | 一种水力旋流器溢流粒度的软测量方法 |
CN107096611A (zh) * | 2017-06-06 | 2017-08-29 | 淮南矿业(集团)有限责任公司 | 一种选煤厂介质制备装置及介质制备方法 |
CN109190226A (zh) * | 2018-08-24 | 2019-01-11 | 东北大学 | 一种磨矿***溢流粒度指标软测量方法 |
CN109635337A (zh) * | 2018-11-13 | 2019-04-16 | 中国矿业大学 | 一种基于块增量随机配置网络的工业过程软测量建模方法 |
CN111026057A (zh) * | 2019-12-16 | 2020-04-17 | 中国建材检验认证集团北京天誉有限公司 | 一种基于vba和深度学习的水泥烧成工艺模型孪生方法及*** |
CN111444580A (zh) * | 2020-03-18 | 2020-07-24 | 武汉理工大学 | 一种大型环件轧制用环坯推演设计方法及*** |
CN115270661A (zh) * | 2022-08-05 | 2022-11-01 | 金川集团股份有限公司 | 水力旋流器溢流细度软测量方法、装置、电子设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1144631A (ja) * | 1997-07-28 | 1999-02-16 | Nikkiso Co Ltd | 湿式系粒度分布測定装置の試料循環器 |
CN1525153A (zh) * | 2003-09-12 | 2004-09-01 | 东北大学 | 球磨机磨矿***溢流粒度指标软测量方法 |
CN1598534A (zh) * | 2004-07-28 | 2005-03-23 | 东北大学 | 基于案例推理的磨矿***溢流粒度指标软测量方法 |
CA2646327A1 (en) * | 2006-03-21 | 2007-09-27 | Exxonmobil Research And Engineering Company | Application of abnormal event detection (aed) technology to polymers process |
-
2010
- 2010-12-28 CN CN 201010609735 patent/CN102169077B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1144631A (ja) * | 1997-07-28 | 1999-02-16 | Nikkiso Co Ltd | 湿式系粒度分布測定装置の試料循環器 |
CN1525153A (zh) * | 2003-09-12 | 2004-09-01 | 东北大学 | 球磨机磨矿***溢流粒度指标软测量方法 |
CN1598534A (zh) * | 2004-07-28 | 2005-03-23 | 东北大学 | 基于案例推理的磨矿***溢流粒度指标软测量方法 |
CA2646327A1 (en) * | 2006-03-21 | 2007-09-27 | Exxonmobil Research And Engineering Company | Application of abnormal event detection (aed) technology to polymers process |
Non-Patent Citations (2)
Title |
---|
《***仿真学报》 20071231 周平等 基于案例推理的磨矿粒度软测量及其软件实现 全文 1 第19卷, 第23期 * |
周平等: "基于案例推理的磨矿粒度软测量及其软件实现", 《***仿真学报》, vol. 19, no. 23, 31 December 2007 (2007-12-31) * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103331205A (zh) * | 2013-06-04 | 2013-10-02 | 广西华锡集团股份有限公司 | 一种模糊闭路自整定动态水力旋流分级机 |
CN103331205B (zh) * | 2013-06-04 | 2015-08-19 | 广西华锡集团股份有限公司 | 一种模糊闭路自整定动态水力旋流分级机 |
CN103969412A (zh) * | 2014-04-13 | 2014-08-06 | 北京工业大学 | 一种基于群决策案例推理的溶解氧浓度软测量方法 |
CN103969412B (zh) * | 2014-04-13 | 2015-11-11 | 北京工业大学 | 一种基于群决策案例推理的溶解氧浓度软测量方法 |
CN105117553A (zh) * | 2015-09-06 | 2015-12-02 | 长沙有色冶金设计研究院有限公司 | 一种水力旋流器溢流粒度的软测量方法 |
CN107096611A (zh) * | 2017-06-06 | 2017-08-29 | 淮南矿业(集团)有限责任公司 | 一种选煤厂介质制备装置及介质制备方法 |
CN109190226A (zh) * | 2018-08-24 | 2019-01-11 | 东北大学 | 一种磨矿***溢流粒度指标软测量方法 |
CN109190226B (zh) * | 2018-08-24 | 2023-01-31 | 东北大学 | 一种磨矿***溢流粒度指标软测量方法 |
CN109635337A (zh) * | 2018-11-13 | 2019-04-16 | 中国矿业大学 | 一种基于块增量随机配置网络的工业过程软测量建模方法 |
CN109635337B (zh) * | 2018-11-13 | 2021-04-27 | 中国矿业大学 | 一种基于块增量随机配置网络的工业过程软测量建模方法 |
CN111026057A (zh) * | 2019-12-16 | 2020-04-17 | 中国建材检验认证集团北京天誉有限公司 | 一种基于vba和深度学习的水泥烧成工艺模型孪生方法及*** |
CN111026057B (zh) * | 2019-12-16 | 2020-12-15 | 中国建材检验认证集团北京天誉有限公司 | 一种基于vba和深度学习的水泥烧成工艺模型孪生方法及*** |
CN111444580A (zh) * | 2020-03-18 | 2020-07-24 | 武汉理工大学 | 一种大型环件轧制用环坯推演设计方法及*** |
CN115270661A (zh) * | 2022-08-05 | 2022-11-01 | 金川集团股份有限公司 | 水力旋流器溢流细度软测量方法、装置、电子设备及介质 |
CN115270661B (zh) * | 2022-08-05 | 2024-07-12 | 金川集团股份有限公司 | 水力旋流器溢流细度软测量方法、装置、电子设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102169077B (zh) | 2013-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102169077B (zh) | 湿式磨矿过程溢流粒度指标混合智能软测量方法 | |
CN1307415C (zh) | 基于案例推理的磨矿***溢流粒度指标软测量方法 | |
CN100401213C (zh) | 选矿过程综合生产指标智能优化控制方法 | |
CN101244403B (zh) | 一种磨矿分级过程优化控制方法 | |
CN103412489B (zh) | 一种磨矿粒度在线预报***及方法 | |
CN104778378B (zh) | 一种油气田产量递减影响因素分析方法 | |
CN101382556B (zh) | 一种基于数据驱动的煤泥浮选精煤灰份软测量方法 | |
CN105278520A (zh) | 基于t-kprm复杂工业过程运行状态评价方法及应用 | |
CN104062953B (zh) | 金湿法冶金过程运行状态在线评价方法 | |
CN103530818B (zh) | 一种基于brb***的供水管网建模方法 | |
CN109190226A (zh) | 一种磨矿***溢流粒度指标软测量方法 | |
CN110322014A (zh) | 一种基于bp神经网络的水泥成品比表面积预测方法 | |
CN109242194A (zh) | 一种基于混合模型的浓密机底流浓度预测方法 | |
CN100394163C (zh) | 球磨机磨矿***溢流粒度指标软测量方法 | |
CN104199389B (zh) | 金湿法冶金置换过程的优化方法 | |
CN101954309B (zh) | 一种磨矿分级过程的抗扰动控制装置及其方法 | |
CN104978484B (zh) | 基于模糊预测模型的选矿厂磨矿过程矿浆浓度检测方法 | |
CN102357395A (zh) | 一种磨矿分级过程的模糊监督控制装置及方法 | |
CN100570327C (zh) | 测量水力旋流器固液分离溢流粒度分布的支持向量机方法 | |
CN108536979A (zh) | 一种基于浓密机机理模型的底流浓度预测方法 | |
CN110222825A (zh) | 一种水泥成品比表面积预测方法及*** | |
CN104199285A (zh) | 一种湿法冶金金氰化浸出过程浸出率预测方法 | |
CN106485326A (zh) | 一种矿石破碎生产过程中的硬度检测方法 | |
CN105867117A (zh) | 基于决策规则提取的选矿厂磨矿分级***节能优化方法 | |
CN106886154B (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 |