发明名称 基于决策规则提取的选矿厂磨矿分级系统节能优化方法
摘要 本发明公开了一种基于决策规则提取的选矿厂磨矿分级系统节能优化方法,采用基于智能控制的模糊控制技术,对选矿厂磨矿分级系统中给矿量、返砂水量和排矿水量三个控制回路的设定值进行优化估计。该优化方法中决策规则集是在对生产过程记录数据使用面向数据的决策列联表构建方法基础上建立,完备决策规则集的增量更新是在对现场增量数据使用基于决策信息的增量更新和完备化方法的基础上实现。本发明所提方法避免了模糊控制传统方法的主观性和局限性,能够为选矿厂磨矿分级系统关键控制回路设定值估计提供可靠的优化结果,为实现磨矿分级系统的节能优化控制奠定基础。
申请公布号 CN105867117A 申请公布日期 2016.08.17
申请号 CN201610445622.4 申请日期 2016.06.20
申请人 西安电子科技大学 发明人 诸文智
分类号 G05B13/02(2006.01)I 主分类号 G05B13/02(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 徐文权
主权项 基于决策规则提取的选矿厂磨矿分级系统节能优化方法,其特征在于,包括以下步骤:1)根据选矿厂磨矿分级系统在人工操作状态下,利用模拟量采集模块采集记录数据信息,以给矿量M<sub>o</sub>、返砂水量W<sub>rs</sub>、排矿水量W<sub>ca</sub>、磨机电流I<sub>m</sub>及分级机电流I<sub>c</sub>组成模式属性集,分别以给矿量设定值<img file="FDA0001023343100000011.GIF" wi="116" he="58" />返砂水量设定值<img file="FDA0001023343100000012.GIF" wi="79" he="62" />及排矿水量设定值<img file="FDA0001023343100000013.GIF" wi="79" he="63" />组成决策属性集,按照系统预设时间t<sub>1</sub>进行数据采样记录,构成现场历史数据库;对已有历史数据库,应用统计预处理技术,获得采样时间t<sub>2</sub>=4t<sub>1</sub>的预选数据库,该数据库在t<sub>2</sub>时刻维度<img file="FDA0001023343100000014.GIF" wi="51" he="61" />上的取值为:<maths num="0001"><math><![CDATA[<mrow><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>t</mi><mn>1</mn></msub><msub><mi>t</mi><mn>2</mn></msub></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>t</mi><mn>2</mn></msub><mo>/</mo><msub><mi>t</mi><mn>1</mn></msub></mrow></munderover><msub><mi>F</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>&CenterDot;</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001023343100000015.GIF" wi="439" he="143" /></maths>其中,F<sub>i</sub>、<img file="FDA0001023343100000016.GIF" wi="781" he="86" />分别表示现场历史数据库和预选数据库D中的对应维度,F<sub>i</sub>(k·t<sub>1</sub>)表示现场历史数据库在k·t<sub>1</sub>时刻维度F<sub>i</sub>上的取值;而后,对数据实例数量为M预选数据库D进行归一化处理,数据实例x<sub>j</sub>在维度<img file="FDA0001023343100000017.GIF" wi="48" he="62" />上的归一化取值<img file="FDA0001023343100000018.GIF" wi="61" he="80" />通过下式进行计算:<maths num="0002"><math><![CDATA[<mrow><msubsup><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>j</mi><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>x</mi><mi>j</mi><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub></msubsup><mo>-</mo><msubsup><mi>V</mi><mi>min</mi><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub></msubsup></mrow><mrow><msubsup><mi>V</mi><mi>max</mi><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub></msubsup><mo>-</mo><msubsup><mi>V</mi><mi>min</mi><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub></msubsup></mrow></mfrac><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>M</mi></mrow>]]></math><img file="FDA0001023343100000019.GIF" wi="614" he="153" /></maths>式中,<img file="FDA00010233431000000110.GIF" wi="61" he="78" />表示数据实例x<sub>j</sub>在维度<img file="FDA00010233431000000111.GIF" wi="48" he="62" />上的原始取值,<img file="FDA00010233431000000112.GIF" wi="84" he="71" />和<img file="FDA00010233431000000113.GIF" wi="75" he="71" />分别表示M个数据实例在维度<img file="FDA00010233431000000114.GIF" wi="44" he="68" />上所取原始值的最大值和最小值;通过上述处理后,预选数据库D在模式属性集{M<sub>o</sub>,W<sub>rs</sub>,W<sub>ca</sub>,I<sub>m</sub>,I<sub>c</sub>}和决策属性集<img file="FDA00010233431000000115.GIF" wi="315" he="87" />上,每一维度的值域为[0,1];2)读取预选数据库D中已有M个数据实例,并组成初始数据集DS<sub>0</sub>;工业数据库继续存储数据,并形成m个数据规模固定为L的增量数据集,其中第k个增量数据集DS<sub>k</sub>表示为:DS<sub>k</sub>={x<sub>(k‑1)M+1</sub>,x<sub>(k‑1)M+2</sub>,…,x<sub>(k‑1)M+L</sub>}其中,k∈{1,2,…,m},x<sub>j</sub>和j∈{(k‑1)M+1,(k‑1)M+2,…,(k‑1)M+L}表示增量数据集DS<sub>k</sub>中的一个数据实例,M为增量数据集DS<sub>k</sub>起始位置标示,L为增量数据集DS<sub>k</sub>的实例数量,m为大于1的自然数,k为小于等于m的自然数;3)设定预选数据库D中给矿量M<sub>o</sub>、返砂水量W<sub>rs</sub>、排矿水量W<sub>ca</sub>、磨机电流I<sub>m</sub>、分级机电流I<sub>c</sub>、给矿量设定值<img file="FDA0001023343100000021.GIF" wi="119" he="64" />返砂水量设定值<img file="FDA0001023343100000022.GIF" wi="80" he="63" />及排矿水量设定值<img file="FDA0001023343100000023.GIF" wi="81" he="63" />的语言变量分别为MO、WRS、WCA、IM、IC、MO<sub>sv</sub>、WRS<sub>sv</sub>和WCA<sub>sv</sub>,其论域为[0,1];语言变量MO、WRS、WCA、IM、IC、MO<sub>sv</sub>、WRS<sub>sv</sub>和WCA<sub>sv</sub>的语言值均为{NB,NS,Z,PS,PB},NB、NS、Z、PS和PB分别代表负大、负小、零、正小和正大;以模式属性集{M<sub>o</sub>,W<sub>rs</sub>,W<sub>ca</sub>,I<sub>m</sub>,I<sub>c</sub>}分别与决策属性<img file="FDA0001023343100000024.GIF" wi="232" he="63" />和<img file="FDA0001023343100000025.GIF" wi="79" he="63" />组成给矿量设定数据集、返砂水量设定数据集和排矿水量设定数据集;4)根据步骤3)对选矿厂磨矿分级系统预选数据库所构建的语言变量,根据定义的语言值进行预选数据库D的维度扩展;如:根据语言变量MO的语言值,将预选数据库的维M<sub>o</sub>扩展成维<img file="FDA0001023343100000026.GIF" wi="131" he="70" />维<img file="FDA0001023343100000027.GIF" wi="130" he="70" />维<img file="FDA0001023343100000028.GIF" wi="106" he="68" />维<img file="FDA0001023343100000029.GIF" wi="96" he="70" />和维<img file="FDA00010233431000000210.GIF" wi="118" he="70" />并称维<img file="FDA00010233431000000211.GIF" wi="131" he="71" />维<img file="FDA00010233431000000212.GIF" wi="131" he="70" />维<img file="FDA00010233431000000213.GIF" wi="111" he="63" />维<img file="FDA00010233431000000214.GIF" wi="92" he="70" />和维<img file="FDA00010233431000000215.GIF" wi="94" he="63" />为维M<sub>o</sub>的扩展维;接着,根据维M<sub>o</sub>上语言值隶属度函数,以及预选数据库中每一实例在维M<sub>o</sub>上取值,确定每个实例在维<img file="FDA00010233431000000216.GIF" wi="125" he="62" />维<img file="FDA00010233431000000217.GIF" wi="128" he="71" />维<img file="FDA00010233431000000218.GIF" wi="107" he="71" />维<img file="FDA00010233431000000219.GIF" wi="97" he="71" />和维<img file="FDA00010233431000000220.GIF" wi="99" he="71" />上的隶属度值;语言变量MO、WRS、WCA、IM和IC采用论域为[0,1]的高斯型隶属度函数族,MO<sub>sv</sub>、WRS<sub>sv</sub>和WCA<sub>sv</sub>中的语言值均采用论域为[0,1]的等腰三角型隶属度函数族;采用相同方法,对预选数据库D中其它维进行扩展,并确定每个实例在给矿量设定数据集、返砂水量设定数据集和排矿水量设定数据集中对应扩展维上的隶属度值;5)使用面向数据的决策列联表构建方法,基于初始数据集DS<sub>0</sub>进行决策列联表的构建,根据论域空间在各模式上的决策信息分布进行划分,并将最终形成的决策列联表用于搜索相应控制回路设定值优化方案;分别构建出给矿量设定列联表、返砂水量设定列联表和排矿水量设定列联表;6)针对步骤5)所形成的给矿量设定列联表、返砂水量设定列联表和排矿水量设定列联表,使用基于决策信息的增量更新和完备化方法,基于增量数据集DS<sub>k</sub>构建并更新各决策规则集,并依据规则集先验知识推理计算出给矿量设定规则集、返砂水量设定规则集和排矿水量设定规则集的完备形式;7)选矿厂磨矿分级系统通过测量模块测量当前工况下的给矿量M<sub>o</sub>、返砂水量W<sub>rs</sub>、排矿水量W<sub>ca</sub>、磨机电流I<sub>m</sub>及分级机电流I<sub>c</sub>的运行值且分别由m′<sub>o</sub>、w′<sub>rs</sub>、w′<sub>ca</sub>、i′<sub>m</sub>和i′<sub>c</sub>来表征;根据步骤4)中论域选择为[0,1]的等腰三角形隶属度函数族,计算M′<sub>o</sub>、W′<sub>rs</sub>、W′<sub>ca</sub>、I′<sub>m</sub>和I′<sub>c</sub>分别在不同扩展维上的模糊隶属度值,其中M′<sub>o</sub>、W′<sub>rs</sub>、W′<sub>ca</sub>、I′<sub>m</sub>和I′<sub>c</sub>分别为运行值m′<sub>o</sub>、w′<sub>rs</sub>、w′<sub>ca</sub>、i′<sub>m</sub>和i′<sub>c</sub>对应的模糊语言值;选择各语言变量下获得最大隶属度值的语言值组合形成运行前件集PAS,该前件集为{M′<sub>o</sub>,W′<sub>rs</sub>,W′<sub>ca</sub>,I′<sub>m</sub>,I′<sub>c</sub>};8)搜索由步骤6)获得的规则集,如:给矿量设定规则集,匹配得到给矿量设定运行规则<img file="FDA0001023343100000031.GIF" wi="158" he="71" />该规则<img file="FDA0001023343100000032.GIF" wi="110" he="71" />的前件模式等于运行前件集PAS={M′<sub>o</sub>,W′<sub>rs</sub>,W′<sub>ca</sub>,I′<sub>m</sub>,I′<sub>c</sub>},运行规则<img file="FDA0001023343100000033.GIF" wi="114" he="71" />中规则后件为给矿量设定值语言变量MO<sub>sv</sub>语言值{NB,NS,Z,PS,PB}上对应的实例条件概率P<sub>lj</sub>,j=1,2,...,5;根据步骤4)中论域选择为[0,1]的等腰三角型隶属度函数族,并通过下式计算运行规则<img file="FDA0001023343100000034.GIF" wi="109" he="78" />对应给矿量设定值的去模糊化值<img file="FDA0001023343100000035.GIF" wi="138" he="62" /><maths num="0003"><math><![CDATA[<mrow><msubsup><mi>MO</mi><mrow><mi>s</mi><mi>v</mi></mrow><mi>l</mi></msubsup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>5</mn></munderover><msup><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>j</mi></msup><msub><mi>P</mi><mrow><mi>l</mi><mi>j</mi></mrow></msub></mrow>]]></math><img file="FDA0001023343100000036.GIF" wi="318" he="134" /></maths>上式中,<img file="FDA0001023343100000037.GIF" wi="315" he="63" />为语言变量MO<sub>sv</sub>对应语言值{NB,NS,Z,PS,PB}在论域[0,1]上定义的等腰三角型隶属度函数的中心;根据步骤1)中归一化公式,通过下式将论域为[0,1]的<img file="FDA0001023343100000041.GIF" wi="110" he="70" />变换回其原始值域内,并表示为给矿量运行设定值MO′<sub>sv</sub>:<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>MO</mi><mrow><mi>s</mi><mi>v</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><msubsup><mi>MO</mi><mrow><mi>s</mi><mi>v</mi></mrow><mi>l</mi></msubsup><mrow><mo>(</mo><msubsup><mi>V</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><msubsup><mi>M</mi><mi>o</mi><mrow><mi>s</mi><mi>v</mi></mrow></msubsup></msubsup><mo>-</mo><msubsup><mi>V</mi><mi>min</mi><msubsup><mi>M</mi><mi>o</mi><mrow><mi>s</mi><mi>v</mi></mrow></msubsup></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>V</mi><mi>min</mi><msubsup><mi>M</mi><mi>o</mi><mrow><mi>s</mi><mi>v</mi></mrow></msubsup></msubsup></mrow>]]></math><img file="FDA0001023343100000042.GIF" wi="684" he="95" /></maths>上式中,<img file="FDA0001023343100000043.GIF" wi="93" he="71" />和<img file="FDA0001023343100000044.GIF" wi="90" he="71" />分别表示预选数据D中数据实例在维度<img file="FDA0001023343100000045.GIF" wi="83" he="55" />上所取原始值的最大值和最小值;9)对返砂水量设定规则集和排矿水量设定规则集重复步骤8)的处理过程,获得当前工况下的返砂水量运行设定值WRS′<sub>sv</sub>和排矿水量运行设定值WCA′<sub>sv</sub>;如此得到的给矿量运行设定值MO′<sub>sv</sub>、返砂水量运行设定值WRS′<sub>sv</sub>和排矿水量运行设定值WCA′<sub>sv</sub>作为当前工况下选矿厂磨矿分级系统的给矿量、返砂水量和排矿水量的最佳设定值,提供给选矿厂磨矿分级系统的控制器进行设置,使选矿厂磨矿分级系统一直在最佳方式下运行。
地址 710065 陕西省西安市雁塔区太白南路2号