发明名称 一种基于图像分析的废杂铜熔炼炉的控制方法及控制系统
摘要 本发明公开了一种基于图像分析的废杂铜熔炼炉的控制方法及控制系统,其通过采集废杂铜工频熔炼过程铜产品的彩色图像,提取其色调,快速估计铜产品的铜含量,能够及时判断熔炼产品是否达到设定值要求,解决铜含量离线检测时间滞后的问题,既能够显著减少离线检测时用于保温的电能消耗,又能够缩短每一炉废杂铜熔炼过程的时间;同时本发明通过实时采集工频熔炼炉内的烟尘图像,利用烟尘图像亮度指数和背景细节模糊度指数,实时估计引风机转速的设定值,进而自动控制引风机的转速。这样既可以减少熔炼炉热能流失、提高加热电能的利用效率,又可以降低引风机恒定高速运转造成的风机电耗,进一步节省电能。
申请公布号 CN103307896B 申请公布日期 2014.11.05
申请号 CN201310232883.4 申请日期 2013.06.09
申请人 浙江大学 发明人 宋执环;张宏伟
分类号 F27D19/00(2006.01)I 主分类号 F27D19/00(2006.01)I
代理机构 杭州天勤知识产权代理有限公司 33224 代理人 胡红娟
主权项 一种基于图像分析的废杂铜熔炼炉的控制方法,包括如下步骤:对于熔炼炉和引风机的停机控制:A1.在熔炼过程中,以采样间隔T1从熔炼炉中提取n次铜样本对应作为n个参考样本,并采集每个参考样本的样本图像,n为大于1的自然数;A2.对于每个参考样本,通过离线化验确定参考样本对应的铜品位,同时根据参考样本的样本图像,计算出参考样本对应的色调;A3.根据所有参考样本对应的铜品位和色调,拟合出熔炼炉内铜样本关于铜品位和色调的函数模型M1,其表达式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&rho;</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>i</mi></msub><mi>K</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mo>+</mo><mi>b</mi></mrow>]]></math><img file="FDA0000545942700000011.GIF" wi="581" he="179" /></maths>其中:ρ为待测样本的铜品位,K(u<sub>i</sub>,u)为u<sub>i</sub>与u的非线性核函数,u<sub>i</sub>为第i个参考样本对应的色调,u为待测样本对应的色调,b为偏置值,a<sub>i</sub>为u<sub>i</sub>对应的权重系数,i为自然数且1≤i≤n;A4.对于当前时刻,从熔炼炉中提取铜样本作为待测样本,采集待测样本的样本图像,并计算出待测样本对应的色调,进而将待测样本的色调代入函数模型M1中求得待测样本的铜品位;A5.若当前时刻待测样本的铜品位达到预定的铜品位指标,则使熔炼炉及引风机停机;若否,则保持熔炼炉及引风机工作状态;根据以下方法计算参考样本或待测样本对应的色调:首先,从参考样本或待测样本的样本图像中选取ROI;然后,根据ROI的RGB像素的灰度平均值,通过以下算式计算出参考样本或待测样本对应的色调:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&theta;</mi><mo>=</mo><mi>arccos</mi><mo>{</mo><mfrac><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><mrow><mo>(</mo><mi>Red</mi><mo>-</mo><mi>Green</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mi>Red</mi><mo>-</mo><mi>Blue</mi><mo>)</mo></mrow><mo>]</mo></mrow><msup><mrow><mo>[</mo><msup><mrow><mo>(</mo><mi>Red</mi><mo>-</mo><mi>Green</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mrow><mo>(</mo><mi>Red</mi><mo>-</mo><mi>Blue</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>Green</mi><mo>-</mo><mi>Blue</mi><mo>)</mo></mrow><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mfrac><mo>}</mo></mrow>]]></math><img file="FDA0000545942700000012.GIF" wi="1582" he="207" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>&theta;</mi></mtd><mtd><mi>if</mi></mtd><mtd><mi>Blue</mi><mo>&le;</mo><mi>Green</mi></mtd></mtr><mtr><mtd><mn>360</mn><mo>-</mo><mi>&theta;</mi></mtd><mtd><mi>if</mi></mtd><mtd><mi>Blue</mi><mo>></mo><mi>Green</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000545942700000013.GIF" wi="792" he="187" /></maths>其中:H为参考样本或待测样本对应的色调,Red、Green和Blue分别为参考样本或待测样本的样本图像ROI的红色通道像素的灰度平均值、绿色通道像素的灰度平均值和蓝色通道像素的灰度平均值;对于引风机转速的控制:B1.在熔炼过程中,以采样间隔T2采集m帧熔炼炉内的烟尘图像对应作为m帧参考图像,m为大于1的自然数;B2.对于每帧参考图像,通过人工经验估计确定参考图像对应的风机转速,同时计算出参考图像的烟尘图像亮度和背景细节模糊度组成参考图像的特征向量;B3.根据所有参考图像对应的风机转速和特征向量,拟合出熔炼炉内烟尘关于风机转速和特征向量的函数模型M2,其表达式如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&omega;</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&alpha;</mi><mi>j</mi></msub><mi>K</mi><mrow><mo>(</mo><msub><mi>U</mi><mi>j</mi></msub><mo>,</mo><mi>U</mi><mo>)</mo></mrow><mo>+</mo><mi>&beta;</mi></mrow>]]></math><img file="FDA0000545942700000021.GIF" wi="657" he="187" /></maths>其中:ω为待测图像对应的风机转速,K(U<sub>j</sub>,U)为U<sub>j</sub>与U的非线性核函数,U<sub>j</sub>为第j帧参考图像的特征向量,U为待测图像的特征向量,β为偏置值,α<sub>j</sub>为U<sub>j</sub>对应的权重系数,j为自然数且1≤j≤m;B4.对于当前时刻,采集熔炼炉内的烟尘图像作为待测图像,并计算出待测图像的烟尘图像亮度和背景细节模糊度组成待测图像的特征向量,进而将待测图像的特征向量代入函数模型M2中求得待测图像对应的风机转速;B5.根据待测图像对应的风机转速通过变频调速器控制引风机转速;根据以下方法计算参考图像或待测图像的烟尘图像亮度和背景细节模糊度:首先,从参考图像或待测图像中选取ROI;然后,根据以下算式计算出参考图像或待测图像的烟尘图像亮度和背景细节模糊度:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>L</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>255</mn><mo>&times;</mo><mi>X</mi><mo>&times;</mo><mi>Y</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>X</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>Y</mi></munderover><mi>G</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&times;</mo><mn>100</mn><mo>%</mo></mrow>]]></math><img file="FDA0000545942700000022.GIF" wi="944" he="179" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>BBI</mi><mo>=</mo><mfrac><msub><mi>n</mi><mi>edge</mi></msub><mi>N</mi></mfrac><mo>&times;</mo><mn>100</mn><mo>%</mo></mrow>]]></math><img file="FDA0000545942700000023.GIF" wi="508" he="153" /></maths>其中:L和BBI分别为烟尘图像亮度和背景细节模糊度,X和Y分别为参考图像或待测图像的ROI的横向尺寸和纵向尺寸,G(x,y)为参考图像或待测图像的ROI中第x行第y列像素的灰度值,N为参考图像或待测图像的ROI的像素总数,n<sub>edge</sub>为参考图像或待测图像的ROI的边缘像素总数,x为自然数且1≤x≤X,y为自然数且1≤y≤Y。
地址 310027 浙江省杭州市西湖区浙大路38号