发明名称 一种高炉铁水硫含量预报方法
摘要 本发明提供一种高炉铁水硫含量预报方法,包括参数选取、数据预处理、预测算法设计、结果输出四个步骤。以硫含量短期均值、硫含量中期均值、硫含量长期均值、前次铁水对应的炉渣碱度、前次铁水含硅量、入炉焦炭S含量、入炉煤粉S含量作为铁水含硫预报的输入变量,利用高炉形成铁水的化学反应过程,结合RBF神经网络,预报下一次铁水的含硫量,在炉况波动小时依靠均线系统,在炉况波动大时自动加入炉渣碱度及前次铁水含硅量,以较少的自变量,最大限度地消除炉况波动对硫含量的影响,获得较好的预报精度,提高高炉铁水硫含量的预报命中率,实现对生铁中硫含量的有效控制,避免出现生铁中硫含量过高的情况。
申请公布号 CN103320559A 申请公布日期 2013.09.25
申请号 CN201310287621.8 申请日期 2013.07.10
申请人 鞍钢股份有限公司 发明人 孙鹏;车玉满;李连成;郭天永;姚硕;陈国一
分类号 C21B5/00(2006.01)I;C21B7/24(2006.01)I 主分类号 C21B5/00(2006.01)I
代理机构 鞍山华惠专利事务所 21213 代理人 赵长芳
主权项 1.一种高炉铁水硫含量预报方法,其特征在于,以硫含量短期均值、硫含量中期均值、硫含量长期均值、前次铁水对应的炉渣碱度、前次铁水含硅量、入炉焦炭S含量、入炉煤粉S含量作为铁水含硫预报的输入变量,利用高炉形成铁水的化学反应过程,结合RBF神经网络,预报下一次铁水的含硫量,在炉况波动小时依靠均线系统,在炉况波动大时自动加入炉渣碱度及前次铁水含硅量,以较少的自变量,获得较好的预报精度,提高高炉铁水硫含量的预报命中率;包括参数选取、数据预处理、预测算法设计、结果输出四个步骤;参数选取:参数选取采用能够反映高炉铁水中硫含量受高炉化学变化影响的变量,包括硫含量短期均值、硫含量中期均值、硫含量长期均值、前次铁水对应的炉渣碱度、前次铁水含硅量、入炉焦炭硫含量及入炉煤粉硫含量;结合高炉冶炼周期的不同选取和计算不同的时间周期内硫含量的均值,且计算时间周期可以调整;设铁水硫序列为:S={S<sub>1</sub>,S<sub>2</sub>,S<sub>3</sub>,,,,,,S<sub>i</sub>,S<sub>i十1</sub>}其中:S<sub>i</sub>为当前铁水硫含量,%;S<sub>i+1</sub>为要预报的下次铁水硫含量,%;设S<sub>AVG</sub>为铁水硫序列的平均值,<maths num="0001"><![CDATA[<math><mrow><msub><mi>S</mi><mi>AVG</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>i</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>i</mi></munderover><msub><mi>S</mi><mi>j</mi></msub></mrow></math>]]></maths>S<sub>j</sub>为铁水中硫含量;(1)硫含量短期均值:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>S</mi><mi>i</mi><mi>n</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><msub><mi>n</mi><mi>n</mi></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>-</mo><msub><mi>n</mi><mi>n</mi></msub><mo>-</mo><mn>1</mn></mrow><mi>i</mi></munderover><msub><mi>S</mi><mi>j</mi></msub></mrow></math>]]></maths>其中:n<sub>n</sub>为短期均值周期,次;参考值为平均冶炼周期1/4;(2)硫含量中期均值:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>S</mi><mi>i</mi><mi>m</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><msub><mi>n</mi><mi>m</mi></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>-</mo><msub><mi>n</mi><mi>m</mi></msub><mo>-</mo><mn>1</mn></mrow><mi>i</mi></munderover><msub><mi>S</mi><mi>j</mi></msub></mrow></math>]]></maths>其中:n<sub>m</sub>为中期均值周期,次;参考值为平均冶炼周期1/2;(3)硫含量长期均值:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>S</mi><mi>i</mi><mi>l</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><msub><mi>n</mi><mi>l</mi></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>-</mo><msub><mi>n</mi><mi>l</mi></msub><mo>-</mo><mn>1</mn></mrow><mi>i</mi></munderover><msub><mi>S</mi><mi>j</mi></msub></mrow></math>]]></maths>其中:n<sub>l</sub>为长期均值周期,次;参考值为平均冶炼周期;均值系统在炉况稳定时是预测下次铁水硫含量的影响因素;数据预处理:数据预处理利用反正切和反余切函数,体现前次铁水对应的炉渣碱度和前次铁水含硅量作为对炉况波动大时的影响参数,炉渣碱度与铁水硫含量的结果有明显的负相关;(1)炉渣碱度在炉况变动加大时是决定铁水含硫量的重要因素:设炉渣碱度序列为:T={T<sub>1</sub>,T<sub>2</sub>,T<sub>3</sub>,,,,,T<sub>i</sub>}其中:T<sub>i</sub>为前次铁水硫对应炉渣碱度;设T<sub>AVR</sub>为炉渣碱度的平均值<img file="FDA00003487547400023.GIF" wi="191" he="138" />(2)前次铁水含硅量:在炉况变动加大时前次铁水含硅量是决定铁水含硫量的重要因素,特别在硫含量变化剧烈时,硫和硅有明显的负相关,而在硫均值附近硫和硅没有明显的相关性;设铁水含硅量序列为:Si={Si<sub>1</sub>,Si<sub>2</sub>,Si<sub>3</sub>,,,,,Si<sub>i</sub>Si<sub>i+1</sub>}其中:Si<sub>i</sub>为前次铁水含硅量,%;设Si<sub>AVG</sub>为铁水含硅量的平均值:<maths num="0005"><![CDATA[<math><mrow><msub><mi>Si</mi><mi>AVG</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>i</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>i</mi></munderover><msub><mi>Si</mi><mi>j</mi></msub></mrow></math>]]></maths>(3)前次铁水对应的炉渣碱度及前次铁水含硅量影响因素分别为:<img file="FDA00003487547400032.GIF" wi="557" he="105" /><maths num="0006"><![CDATA[<math><mrow><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>Si</mi></msubsup><mo>=</mo><msub><mi>K</mi><mi>Si</mi></msub><mo>*</mo><mfrac><mn>2</mn><mi>&pi;</mi></mfrac><mo>*</mo><mi>Arc</mi><mi>cot</mi><mo>(</mo><mfrac><mrow><msub><mi>Si</mi><mi>i</mi></msub><mo>-</mo><mi>S</mi></mrow><msub><mi>Si</mi><mi>A</mi></msub></mfrac><mo>*</mo></mrow></math>]]></maths>其中:<img file="FDA00003487547400034.GIF" wi="84" he="76" />为前次铁水对应的炉渣碱度影响因素;<img file="FDA00003487547400035.GIF" wi="84" he="74" />为前次铁水含硅量影响因素;K<sub>T</sub>前次铁水对应的炉渣碱度影响因素权重;K<sub>Si</sub>前次铁水含硅量影响因素权重;(4)入炉焦炭S含量、入炉煤粉S含量对铁水硫含量的影响因素分别为:<maths num="0007"><![CDATA[<math><mrow><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>Sj</mi></msubsup><mo>=</mo><msub><mi>K</mi><mi>Sj</mi></msub><mo>*</mo><msub><mi>Sj</mi><mi>i</mi></msub></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>Sm</mi></msubsup><mo>=</mo><msub><mi>K</mi><mi>Sm</mi></msub><mo>*</mo><msub><mi>Sm</mi><mi>i</mi></msub></mrow></math>]]></maths>其中:<img file="FDA00003487547400038.GIF" wi="88" he="83" />为入炉焦炭S含量对铁水硫含量影响因素;<img file="FDA00003487547400039.GIF" wi="86" he="79" />为入炉煤粉S含量对铁水硫含量影响因素;K<sub>Sj</sub>为入炉焦炭S含量对铁水硫含量影响因素权重;K<sub>Sm</sub>为入炉焦炭S含量对铁水硫含量影响因素权重;(5)含硫量短期、中期及长期影响因素由以下三个多项式表示:<maths num="0009"><![CDATA[<math><mrow><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>A</mi></msubsup><mo>=</mo><msub><mi>K</mi><mi>n</mi></msub><mo>*</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>S</mi><mi>i</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>m</mi></msub><mo>*</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>S</mi><mi>i</mi><mi>m</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>l</mi></msub><mo>*</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>S</mi><mi>i</mi><mi>l</mi></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>其中:<img file="FDA000034875474000311.GIF" wi="78" he="80" />为含硫量均线系统影响因素;K<sub>n</sub>短期均线影响因素权重;K<sub>m</sub>中期均线影响因素权重;K<sub>l</sub>长期均线影响因素权重;预测算法设计:预测算法设计根据高炉冶炼的化学原理结合高炉冶炼的物理过程预报下次铁水硫含量;其中,硫含量短期均值、硫含量中期均值、硫含量长期均值、前次铁水对应的炉渣碱度是基于高炉化学反应来考虑的,而前次铁水含硅量、入炉焦炭硫含量、入炉煤粉硫含量则是根据高炉冶炼过程中温度和原燃料<img file="FDA000034875474000312.GIF" wi="184" he="55" />考虑的;硫预报函数为:<maths num="0010"><![CDATA[<math><mrow><msub><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mrow><mo>(</mo><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>A</mi></msubsup><mo>+</mo><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>T</mi></msubsup><mo>+</mo><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>Si</mi></msubsup><mo>+</mo><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>Sj</mi></msubsup><mo>+</mo><msubsup><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>Sm</mi></msubsup><mo>)</mo></mrow><mo>/</mo></mrow></math>]]></maths><maths num="0011"><![CDATA[<math><mrow><mo>[</mo><msub><mi>K</mi><mi>n</mi></msub><mo>+</mo><msub><mi>K</mi><mi>m</mi></msub><mo>+</mo><msub><mi>K</mi><mi>l</mi></msub><mo>+</mo><msub><mi>K</mi><mi>Sj</mi></msub><mo>+</mo><msub><mi>K</mi><mi>Sm</mi></msub><mo>+</mo></mrow></math>]]></maths><maths num="0012"><![CDATA[<math><mrow><msub><mi>K</mi><mi>T</mi></msub><mo>*</mo><mi>Arc</mi><mi>tan</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>T</mi><mi>i</mi></msub><mo>-</mo><msub><mi>T</mi><mi>AVR</mi></msub></mrow><msub><mi>T</mi><mi>AVR</mi></msub></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>Si</mi></msub><mo>*</mo><mi>Arc</mi><mi>cot</mi><mo>(</mo><mfrac><mrow><msub><mi>Si</mi><mi>i</mi></msub><mo>-</mo><mi>S</mi></mrow><msub><mi>Si</mi><mi>A</mi></msub></mfrac><mo>]</mo></mrow></math>]]></maths><maths num="0013"><![CDATA[<math><mrow><mo>=</mo><mo>[</mo><msub><mi>K</mi><mi>n</mi></msub><mo>*</mo><mrow><mo>(</mo><msubsup><mrow><msub><mi>S</mi><mi>i</mi></msub><mo>+</mo><mi>S</mi></mrow><mi>i</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>m</mi></msub><mo>*</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>S</mi><mi>i</mi><mi>m</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>l</mi></msub><mo>*</mo><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>S</mi><mi>i</mi><mi>l</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>Sj</mi></msub><mo>*</mo><msub><mi>Sj</mi><mi>i</mi></msub><mo>+</mo><msub><mi>K</mi><mi>Sm</mi></msub><mo>*</mo><msub><mi>Sm</mi><mi>i</mi></msub><mo>+</mo></mrow></math>]]></maths><maths num="0014"><![CDATA[<math><mrow><msub><mi>K</mi><mi>T</mi></msub><mo>*</mo><mfrac><mn>2</mn><mi>&pi;</mi></mfrac><mo>*</mo><mi>Arc</mi><mi>tan</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>T</mi><mi>i</mi></msub><mo>-</mo><msub><mi>T</mi><mi>AVR</mi></msub></mrow><msub><mi>T</mi><mi>AVR</mi></msub></mfrac><mo>)</mo></mrow><mo>*</mo><msub><mi>T</mi><mi>i</mi></msub><mo>+</mo><msub><mi>K</mi><mi>Si</mi></msub><mo>*</mo><mfrac><mn>2</mn><mi>&pi;</mi></mfrac><mo>*</mo><mi>Arc</mi><mi>cot</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>Si</mi><mi>i</mi></msub><mo>-</mo><msub><mi>Si</mi><mi>AVR</mi></msub></mrow><msub><mi>Si</mi><mi>AVR</mi></msub></mfrac><mo>)</mo></mrow><mo>*</mo><msub><mi>Si</mi><mi>i</mi></msub><mo>]</mo><mo>/</mo></mrow></math>]]></maths><maths num="0015"><![CDATA[<math><mrow><mo>[</mo><msub><mi>K</mi><mi>n</mi></msub><mo>+</mo><msub><mi>K</mi><mi>m</mi></msub><mo>+</mo><msub><mi>K</mi><mi>l</mi></msub><mo>+</mo><msub><mi>K</mi><mi>Sj</mi></msub><mo>+</mo><msub><mi>K</mi><mi>Sm</mi></msub><mo>+</mo></mrow></math>]]></maths><maths num="0016"><![CDATA[<math><mrow><msub><mi>K</mi><mi>T</mi></msub><mo>*</mo><mi>Arc</mi><mi>tan</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>T</mi><mi>i</mi></msub><mo>-</mo><msub><mi>T</mi><mi>AVR</mi></msub></mrow><msub><mi>T</mi><mi>AVR</mi></msub></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>Si</mi></msub><mo>*</mo><mi>Arc</mi><mi>cot</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>Si</mi><mi>i</mi></msub><mo>-</mo><msub><mi>Si</mi><mi>AVR</mi></msub></mrow><msub><mi>Si</mi><mi>AVR</mi></msub></mfrac><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>预测算法中的权重系数K<sub>n</sub>,K<sub>m</sub>,K<sub>l</sub>,K<sub>T</sub>,K<sub>Si</sub>,K<sub>Sj</sub>,K<sub>Sm</sub>,采用RBF神经网络训练获得;结果输出:结果输出以曲线和数字形式显示在操作屏上,指导炉长对相关参数进行调节,达到稳定炉温使炉况顺行的目的。
地址 114021 辽宁省鞍山市铁西区鞍钢厂区内