CN103390098B - 一种脱硫计算方法及其应用的系统 - Google Patents

一种脱硫计算方法及其应用的系统 Download PDF

Info

Publication number
CN103390098B
CN103390098B CN201310268262.1A CN201310268262A CN103390098B CN 103390098 B CN103390098 B CN 103390098B CN 201310268262 A CN201310268262 A CN 201310268262A CN 103390098 B CN103390098 B CN 103390098B
Authority
CN
China
Prior art keywords
sigma
desulfurization
stove
data
dosage
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
Application number
CN201310268262.1A
Other languages
English (en)
Other versions
CN103390098A (zh
Inventor
郑毅
朱诺南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan Iron Steel Engineering Technology Group Automation Co ltd
Wuhan Iron and Steel Group Corp
Original Assignee
Wuhan Iron Steel Engineering Technology Group Automation Co ltd
Wuhan Iron and Steel Group Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan Iron Steel Engineering Technology Group Automation Co ltd, Wuhan Iron and Steel Group Corp filed Critical Wuhan Iron Steel Engineering Technology Group Automation Co ltd
Priority to CN201310268262.1A priority Critical patent/CN103390098B/zh
Publication of CN103390098A publication Critical patent/CN103390098A/zh
Application granted granted Critical
Publication of CN103390098B publication Critical patent/CN103390098B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Refinement Of Pig-Iron, Manufacture Of Cast Iron, And Steel Manufacture Other Than In Revolving Furnaces (AREA)

Abstract

本发明涉及金属冶炼技术领域，特别涉及一种脱硫计算方法及其应用的系统。脱硫计算方法包括采用多元线性回归的方法计算出脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3。系统包括控制器、一级机与二级机。一级机与控制器连接，一级机用于显示脱硫过程中的各项参数以及用于输入控制指令修正本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值。二级机与一级机连接，二级机用于计算本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值及存储计算所需的数据。控制器与一级机连接，控制器用于控制脱硫设备进行备料以及控制脱硫装置进行脱硫。本发明提供的脱硫计算方法及其应用的系统提高了脱硫的成功率以及脱硫的工作效率。

Description

一种脱硫计算方法及其应用的系统
技术领域
本发明涉及金属冶炼技术领域，特别涉及一种脱硫计算方法及其应用的系统。
背景技术
在炼钢领域，一般都有铁水预处理装置，它是关系到钢材质量好坏的重要环节。铁水预处理的主要工艺目标是：脱硫、脱磷和脱硅。而脱硫是铁水预处理过程的最重要的环节。降低铁水中的硫含量，对提高钢的质量有重要的影响。目前，由日本开发KR法为代表的机械搅拌脱硫和喷吹法脱硫，作为两种主要的铁水脱硫手段，在我国炼钢行业中得到广泛地应用。而这两种脱硫法的关键是要准确地根据铁水的重量、铁水的温度、铁水原始含硫量(前硫含量)、铁水成品含硫量(目标硫含量)、脱硫剂成分，得出需要添加的脱硫剂量、需要搅拌的时间、搅拌的速度。传统的KR法计算脱硫剂的重量是利用预先建立好的一个表，通过查表的方法来计算出脱硫剂的加入量。该表也叫做KR脱硫剂加入量计算表，人工查表的方法固然简单，却存在以下问题：一、人工查表的数据本身就是一个经验积累的数据，只能作为加入脱硫剂的参考，数据并不是准确的。二、人工查表的数据没有考虑到搅拌器的使用次数对脱硫的影响，得到结果误差较大。
发明内容
本发明所要解决的技术问题是提供一种能够准确的计算脱硫过程中所需的脱硫剂量、搅拌速度及搅拌时间参数的脱硫计算方法及其应用的系统。
为解决上述技术问题，本发明提供了一种脱硫计算方法，一种铁水脱硫计算方法，包括：根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉铁水原始硫值数据X3及本炉铁水目标硫值数据X4，采用多元线性回归的方法分别计算出脱硫剂量Y1。根据所述本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉脱硫剂量数据X6，采用多元线性回归的方法计算出搅拌速度Y2。根据所述本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉搅拌速度数据X7，采用多元线性回归的方法计算出搅拌时间Y3。
其中，所述根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉铁水原始硫值数据X3及本炉铁水目标硫值数据X4，采用多元线性回归的方法计算出脱硫剂量Y1包括：根据公式
Y1＝a0+a1*X1+a2*X2+a3*X3+a4*X4及公式
Σ Y 1 = n * a 0 + a 1 * Σ X 1 + a 2 * Σ X 2 + a 3 * Σ X 3 + a 4 * Σ X 4 Σ X 1 Y 1 = a 0 * Σ X 1 + a 1 * Σ X 1 2 + a 2 * Σ X 1 X 2 + a 3 * Σ X 1 X 3 + a 4 * Σ X 1 X 4 Σ X 2 Y 1 = a 0 * Σ X 2 + a 1 * Σ X 1 X 2 + a 2 * Σ X 2 2 + a 3 * Σ X 2 X 3 + a 4 * Σ X 2 X 4 Σ X 3 Y 1 = a 0 * Σ X 3 + a 1 * Σ X 1 X 3 + a 2 * Σ X 2 X 3 + a 3 * Σ X 3 2 + a 4 * Σ X 3 X 4 Σ X 4 Y 1 = a 0 * Σ X 4 + a 1 * Σ X 1 X 4 + a 2 * Σ X 2 X 4 + a 3 * Σ X 3 X 4 + a 4 * Σ X 4 2
计算出脱硫剂量Y1的值。所述n为每炉中成功脱硫的历史次数，所述a0、a1、a2、a3及a4为需要计算出的参数，所述∑X1、∑X2、∑X3、∑X4及所述∑Y1分别为对应量的n次历史数据之和，所述∑X1X2、∑X1X3、∑X2X4、∑X3X4、∑X1Y1、∑X2Y1、∑X3Y1及所述∑X4Y1分别为对应量的n次历史数据乘积之和。
进一步地，所述根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉脱硫剂量数据X6，采用多元线性回归的方法计算出搅拌速度Y2包括：根据公式
Y2＝b0+b1*X1+b2*X2+b3*X5+b4*X6及公式
Σ Y 2 = n * b 0 + b 1 * Σ X 1 + b 2 * Σ X 2 + b 3 * Σ X 5 + b 4 * Σ X 6 Σ X 1 Y 2 = b 0 * Σ X 1 + b 1 * ΣX 1 2 + b 2 * Σ X 1 X 2 + b 3 * Σ X 1 X 5 + b 4 * Σ X 1 X 6 Σ X 2 Y 2 = b 0 * Σ X 2 + b 1 * Σ X 1 X 2 + b 2 * Σ X 2 2 + b 3 * Σ X 2 X 5 + b 4 * Σ X 2 X 6 Σ X 5 Y 2 = b 0 * Σ X 5 + b 1 * Σ X 1 X 5 + b 2 * Σ X 2 X 5 + b 3 * Σ X 5 2 + b 4 * Σ X 5 X 6 Σ X 6 Y 2 = b 0 * Σ X 6 + b 1 * Σ X 1 X 6 + b 2 * Σ X 2 X 6 + b 3 * Σ X 5 X 6 + b 4 * Σ X 6 2
计算出搅拌速度Y2的值。本炉脱硫剂量数据X6与Y1的值相等；所述n为每炉中成功脱硫的历史次数，所述b0、b1、b2、b3及b4为需要计算出的参数，所述∑X1、∑X2、∑X5、∑X6及所述∑Y2分别为对应量的n次历史数据之和，所述∑X1X2、∑X1X5、∑X2X6、∑X5X6、∑X1Y2、∑X2Y2、∑X5Y5及所述∑X6Y2分别为对应量的n次历史数据乘积之和。
进一步地，所述根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉搅拌速度数据X7，采用多元线性回归的方法计算出搅拌时间Y3包括：根据公式
Y3＝c0+c1*X1+c2*X2+c3*X5+c4*X7及公式
Σ Y 3 = n * c 0 + c 1 * Σ X 1 + c 2 * Σ X 2 + c 3 * Σ X 5 + c 4 * Σ X 7 Σ X 1 Y 3 = c 0 * Σ X 1 + c 1 * Σ X 1 2 + c 2 * Σ X 1 X 2 + c 3 * Σ X 1 X 5 + c 4 * Σ X 1 X 7 Σ X 2 Y 3 = c 0 * Σ X 2 + c 1 * Σ X 1 X 2 + c 2 * Σ X 2 2 + c 3 * Σ X 2 X 5 + c 4 * Σ X 2 X 7 Σ X 5 Y 3 = c 0 * Σ X 5 + c 1 * Σ X 1 X 5 + c 2 * Σ X 2 X 5 + c 3 * Σ X 5 2 + c 4 * Σ X 5 X 7 Σ X 7 Y 3 = c 0 * Σ X 7 + c 1 * Σ X 1 X 7 + c 2 * Σ X 2 X 7 + c 3 * Σ X 5 X 7 + c 4 * Σ X 7 2
计算出搅拌时间Y3的值。所述n为每炉中成功脱硫的历史次数，所述c0、c1、c2、c3及c4为需要计算出的参数，所述∑X1、∑X2、∑X5、∑X7及所述∑Y3分别为对应量的n次历史数据之和，所述∑X1X2、∑X1X5、∑X2X7、∑X5X7、∑X1Y3、∑X2Y3、∑X5Y3及所述∑X7Y3分别为对应量的n次历史数据乘积之和。
本发明还提供了一种铁水脱硫计算方法所应用的系统，包括一级机、二级机及控制器。所述一级机与所述控制器连接，且所述一级机用于显示脱硫过程中的各项参数以及用于输入控制指令修正本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值。所述二级机与所述一级机连接，所述二级机用于存储KR脱硫剂加入量计算表、每炉中成功脱硫的历史数据及其它相关数据，并计算本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值。所述控制器与所述一级机连接，所述控制器用于根据脱硫剂量Y1的值控制脱硫设备进行备料，根据脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值控制脱硫装置进行脱硫。
进一步地，所述二级机包括：数据存储单元，用于存储KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据，并将KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据发送给数据处理单元。数据处理单元，用于确定本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值。
进一步地，所述数据存储单元包括：存储模块，用于存储KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据。发送模块，用于将存储模块中存储的所述KR脱硫剂加入量计算表及每炉中成功脱硫的所述历史数据发送给所述数据处理单元。
进一步地，所述数据处理单元包括：接收模块，用于接收所述发送模块传来的所述KR脱硫剂加入量计算表及每炉中成功脱硫的所述历史数据，并将所述KR脱硫剂加入量计算表及所述历史数据发送给处理模块。处理模块，根据所述KR脱硫剂加入量计算表及每炉中成功脱硫的所述历史数据查询出本次的脱硫剂量Y1或利用权利要求4所述的多元线性回归的方法分别计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值；并将本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值通过所述一级机发送给所述控制器，以及存储到所述存储模块中。
进一步地，所述处理模块包括：查询模块，如果本次脱硫工作之前没有任何脱硫工作，则查询模块根据KR脱硫剂加入量表可以查询出本次脱硫的脱硫剂量Y1的值；并将本次脱硫的脱硫剂量Y1的值通过所述一级机发送给所述控制器，以及存储到所述存储模块中。计算模块，如果本次脱硫工作之前有脱硫工作，利用权利要求4所述的多元线性回归的方法分别计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值；并将本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值通过所述一级机发送给所述控制器，以及存储到所述存储模块中。
本发明提供的铁水脱硫计算方法及其应用的系统，采用每炉中成功脱硫的历史数据进行多元线性回归计算，得到本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3，使得本次脱硫的所需的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的参数更加的准确，提高了脱硫的成功率；同时，二级机利用本发明提供的铁水脱硫计算方法，自动计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的参数，并将这些参数提供给控制器用于控制脱硫设备的运转，实现了铁水脱硫过程中的自动化控制，提高了铁水脱硫的工作效率。
,附图说明
图1为本发明实施例提供的铁水温度与脱硫率的关系图。
图2为本发明实施例提供的脱硫剂与脱硫率的关系图。
图3为本发明实施例提供的搅拌时间与铁水硫含量关系图。
图4为本发明实施例提供的初始硫含量与脱硫率关系图。
图5为本发明实施例提供的系统内部软件单元示意图。
具体实施方式
在脱硫领域，脱硫率与铁水温度、脱硫剂、搅拌时间及搅拌速度等因素密切相关。铁水温度与脱硫率是成正比的线性关系。参见图1，脱硫剂加入量与脱硫率的关系曲线在一定范围内是成正比的线性关系，但到了一定程度，就到达饱和状态。参见图2，搅拌时间随铁水硫含量的增加而延长，在保证脱硫目标的前提下，搅拌时间与铁水硫含量有关，通常脱硫反应进行10min左右可达到最佳脱硫效果。搅拌速度对脱硫率的影响，脱硫时加快搅拌头旋转速度，就加快了脱硫反应的传质，也加快了脱硫反应速度。参见图3，随着铁水初始硫含量的增加，脱硫率越高，这是因为铁水初始硫含量越高，硫的活度越大，脱硫反应越容易进行，脱硫率越高；在初始硫大于0.030％后，脱硫率基本稳定。
基于以上分析，本发明实施例提供了一种脱硫计算方法，包括：
步骤10：根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉铁水原始硫值数据X3及本炉铁水目标硫值数据X4，采用多元线性回归的方法计算出脱硫剂量Y1，具体为：根据公式
Y1＝a0+a1*X1+a2*X2+a3*X3+a4*X4及公式
Σ Y 1 = n * a 0 + a 1 * Σ X 1 + a 2 * Σ X 2 + a 3 * Σ X 3 + a 4 * Σ X 4 Σ X 1 Y 1 = a 0 * Σ X 1 + a 1 * Σ X 1 2 + a 2 * Σ X 1 X 2 + a 3 * Σ X 1 X 3 + a 4 * Σ X 1 X 4 Σ X 2 Y 1 = a 0 * Σ X 2 + a 1 * Σ X 1 X 2 + a 2 * Σ X 2 2 + a 3 * Σ X 2 X 3 + a 4 * Σ X 2 X 4 Σ X 3 Y 1 = a 0 * Σ X 3 + a 1 * Σ X 1 X 3 + a 2 * Σ X 2 X 3 + a 3 * Σ X 3 2 + a 4 * Σ X 3 X 4 Σ X 4 Y 1 = a 0 * Σ X 4 + a 1 * Σ X 1 X 4 + a 2 * Σ X 2 X 4 + a 3 * Σ X 3 X 4 + a 4 * Σ X 4 2
计算出脱硫剂量Y1的值。其中，n为每炉中成功脱硫的历史次数，a0、a1、a2、a3及a4为需要计算出的参数。需要说明的是，历史铁水温度数据之和∑X1、历史铁水重量数据之和∑X2、历史铁水原始硫值数据之和∑X3、历史铁水目标硫值数据之和∑X4及历史脱硫剂量之和∑Y1是指在本次脱硫之前的一次或多次脱硫过程中留下的数据，例如，∑X1为若干次历史成功脱硫过程中的铁水温度数据之和。
步骤20：根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉脱硫剂量数据X6，采用多元线性回归的方法计算出搅拌速度Y2，具体为：根据公式
Y2＝b0+b1*X1+b2*X2+b3*X5+b4*X6及公式
Σ Y 2 = n * b 0 + b 1 * Σ X 1 + b 2 * Σ X 2 + b 3 * Σ X 5 + b 4 * Σ X 6 Σ X 1 Y 2 = b 0 * Σ X 1 + b 1 * ΣX 1 2 + b 2 * Σ X 1 X 2 + b 3 * Σ X 1 X 5 + b 4 * Σ X 1 X 6 Σ X 2 Y 2 = b 0 * Σ X 2 + b 1 * Σ X 1 X 2 + b 2 * Σ X 2 2 + b 3 * Σ X 2 X 5 + b 4 * Σ X 2 X 6 Σ X 5 Y 2 = b 0 * Σ X 5 + b 1 * Σ X 1 X 5 + b 2 * Σ X 2 X 5 + b 3 * Σ X 5 2 + b 4 * Σ X 5 X 6 Σ X 6 Y 2 = b 0 * Σ X 6 + b 1 * Σ X 1 X 6 + b 2 * Σ X 2 X 6 + b 3 * Σ X 5 X 6 + b 4 * Σ X 6 2
计算出搅拌速度Y2的值。其中，X6为本炉脱硫剂量数据Y1，需要说明的是，历史搅拌器使用次数数据之和∑X5、历史脱硫剂量数据之和∑X6及历史搅拌速度之和∑Y2是本次脱硫之前的一次或多次脱硫过程中留下的数据。n为每炉中成功脱硫的历史次数，b0、b1、b2、b3及b4为需要计算出的参数，∑X1、∑X2、∑X5及∑X6分别为对应量的n次历史数据之和，∑X1X2、∑X1X5、∑X2X6及∑X5X6分别为对应量的n次历史数据乘积之和。
步骤30：根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉搅拌速度数据X7，采用多元线性回归的方法计算出搅拌时间Y3，具体为：根据公式
Y3＝c0+c1*X1+c2*X2+c3*X5+c4*X7及公式
Σ Y 3 = n * c 0 + c 1 * Σ X 1 + c 2 * Σ X 2 + c 3 * Σ X 5 + c 4 * Σ X 7 Σ X 1 Y 3 = c 0 * Σ X 1 + c 1 * Σ X 1 2 + c 2 * Σ X 1 X 2 + c 3 * Σ X 1 X 5 + c 4 * Σ X 1 X 7 Σ X 2 Y 3 = c 0 * Σ X 2 + c 1 * Σ X 1 X 2 + c 2 * Σ X 2 2 + c 3 * Σ X 2 X 5 + c 4 * Σ X 2 X 7 Σ X 5 Y 3 = c 0 * Σ X 5 + c 1 * Σ X 1 X 5 + c 2 * Σ X 2 X 5 + c 3 * Σ X 5 2 + c 4 * Σ X 5 X 7 Σ X 7 Y 3 = c 0 * Σ X 7 + c 1 * Σ X 1 X 7 + c 2 * Σ X 2 X 7 + c 3 * Σ X 5 X 7 + c 4 * Σ X 7 2
计算出搅拌时间Y3的值。X7为本炉搅拌速度数据Y2，n为每炉中成功脱硫的历史次数，c0、c1、c2、c3及c4为需要计算出的参数，∑X1、∑X2、∑X5及∑X7分别为对应量的n次历史数据之和，所述∑X1X2、∑X1X5、∑X2X7及∑X5X7分别为对应量的n次历史数据乘积之和。
步骤40：每进行一次成功脱硫，则会将铁水温度数据X1、铁水重量数据X2、铁水原始硫值数据X3、铁水目标硫值数据X4、脱硫剂量Y1搅拌速度、Y2及搅拌时间Y3的值记录到历史数据库中。
本发明实施例还提供了一种上述脱硫方法所应用的系统，参见图5，该系统包括控制器、一级机及二级机。控制器通过数据导线与一级机连接，一级机与二级机通过以太网连接。二级机设置有：数据存储单元，用于存储KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据。数据处理单元，用于确定本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值。数据存储单元包括：存储模块，用于存储KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据。发送模块，用于将存储模块中存储的KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据发送给数据处理单元。数据处理单元包括：接收模块，用于接收发送模块传来的KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据，并将KR脱硫剂加入量计算表及历史数据发送给处理模块。处理模块，根据KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据查询出本次的脱硫剂量Y1或利用本发明实施例提供的多元线性回归的方法分别计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值；并将本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值通过一级机发送给控制器，以及存储到存储模块中。处理模块包括：查询模块，如果本次脱硫工作之前没有任何脱硫工作，则查询模块根据KR脱硫剂加入量表可以查询出本次脱硫的脱硫剂量Y1的值；并将本次脱硫的脱硫剂量Y1的值通过一级机发送给控制器，以及存储到存储模块中。计算模块，如果本次脱硫工作之前有脱硫工作，利用权利要求4所述的多元线性回归的方法分别计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值；并将本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值通过一级机发送给所述控制器，以及存储到存储模块中。控制器用于根据脱硫剂量Y1的值控制脱硫设备进行备料，根据脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值控制脱硫装置进行脱硫。一级机用于显示脱硫过程中的各项参数以及对控制器发号施令以及用于输入控制指令修正本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值。本实施例中，一级机为HMI设备，二级机为服务器或计算机，控制器为PLC或单片机。
下面以一次脱硫工作为例，对本发明提供的脱硫计算方法所应用的系统的工作原理进行说明，二级机中的存储模块存储KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据，若本次脱硫工作之前没有任何脱硫工作，那么仅存储有KR脱硫剂加入量计算表。在进行本次脱硫工作之前，二级机中的计算模块通过本实施例介绍的多元线性回归的方法，计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值，并由一级机把得到的这些参数(本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值)传给控制器。特别的，如果本次脱硫工作之前没有任何脱硫工作，则查询模块根据KR脱硫剂加入量表可以查询出本次脱硫的脱硫剂量Y1的值，并人工判断输入搅拌速度Y2及搅拌时间Y3的值，并由一级机把得到的这些参数(本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值)传给控制器。控制器根据脱硫剂量Y1的值控制脱硫设备进行备料，备料完成以后，控制器根据脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值控制脱硫装置进行脱硫，当脱硫工作完成以后，控制器将本次脱硫工作使用的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3及搅拌器使用次数X5的值作为历史数据通过一级机存储到二级机中的存储模块。控制器控制脱硫装置进行脱硫的过程中，脱硫进度、脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3及搅拌器使用次数X5等参数通过一级机显示，供工作人员查看。同时，工作人员可以通过一级机设备输入控制指令修正本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值。
本发明实施例具有以下有益效果：
1、采用每炉中成功脱硫的历史数据进行多元线性回归计算，得到本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3，使得本次脱硫的所需的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的参数更加的准确，提高了脱硫的成功率。
2、二级机利用本发明提供的铁水脱硫计算方法，自动计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的参数，并将这些参数提供给控制器用于控制脱硫设备的运转，实现了铁水脱硫过程中的自动化控制，提高了铁水脱硫的工作效率。
3、硬件结构简单、取材方面，便于在工业上广泛实施。
最后所应说明的是，以上具体实施方式仅用以说明本发明的技术方案而非限制，尽管参照实例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种铁水脱硫计算方法，特征在于，包括：
根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉铁水原始硫值数据X3及本炉铁水目标硫值数据X4，采用多元线性回归的方法计算出脱硫剂量Y1；
根据所述本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉脱硫剂量数据X6，采用多元线性回归的方法计算出搅拌速度Y2；
根据所述本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉搅拌速度数据X7，采用多元线性回归的方法计算出搅拌时间Y3；
其中，
所述根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉铁水原始硫值数据X3及本炉铁水目标硫值数据X4，采用多元线性回归的方法计算出脱硫剂量Y1包括：根据公式
Y1＝a0+a1*X1+a2*X2+a3*X3+a4*X4及公式
ΣY 1 = n * a 0 + a 1 * ΣX 1 + a 2 * ΣX 2 + a 3 * ΣX 3 + a 4 * ΣX 4 ΣX 1 Y 1 = a 0 * ΣX 1 + a 1 * ΣX 1 2 + a 2 * ΣX 1 X 2 + a 3 * ΣX 1 X 3 + a 4 * ΣX 1 X 4 ΣX 2 Y 1 = a 0 * ΣX 2 + a 1 * ΣX 1 X 2 + a 2 * ΣX 2 2 + a 3 * ΣX 2 X 3 + a 4 * ΣX 2 X 4 ΣX 3 Y 1 = a 0 * ΣX 3 + a 1 * ΣX 1 X 3 + a 2 * ΣX 2 X 3 + a 3 * ΣX 3 2 + a 4 * ΣX 3 X 4 ΣX 4 Y 1 = a 0 * ΣX 4 + a 1 * ΣX 1 X 4 + a 2 * ΣX 2 X 4 + a 3 * ΣX 3 X 4 + a 4 * ΣX 4 2
计算出脱硫剂量Y1的值；
所述n为每炉中成功脱硫的历史次数，所述a0、a1、a2、a3及a4为需要计算出的参数，所述∑X1、∑X2、∑X3、∑X4及所述∑Y1分别为对应量的n次历史数据之和，所述∑X1X2、∑X1X3、∑X2X4、∑X3X4、∑X1Y1、∑X2Y1、∑X3Y1及所述∑X4Y1分别为对应量的n次历史数据乘积之和。
2.根据权利要求1所述的铁水脱硫计算方法，其特征在于，所述根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉脱硫剂量数据X6，采用多元线性回归的方法计算出搅拌速度Y2包括：根据公式
Y2＝b0+b1*X1+b2*X2+b3*X5+b4*X6及公式
ΣY 2 = n * b 0 + b 1 * ΣX 1 + b 2 * ΣX 2 + b 3 * ΣX 5 + b 4 * ΣX 6 ΣX 1 Y 2 = b 0 * ΣX 1 + b 1 * ΣX 1 2 + b 2 * ΣX 1 X 2 + b 3 * ΣX 1 X 5 + b 4 * ΣX 1 X 6 ΣX 2 Y 2 = b 0 * ΣX 2 + b 1 * ΣX 1 X 2 + b 2 * ΣX 2 2 + b 3 * ΣX 2 X 5 + b 4 * ΣX 2 X 6 ΣX 5 Y 2 = b 0 * ΣX 5 + b 1 * ΣX 1 X 5 + b 2 * ΣX 2 X 5 + b 3 * ΣX 5 2 + b 4 * ΣX 5 X 6 ΣX 6 Y 2 = b 0 * ΣX 6 + b 1 * ΣX 1 X 6 + b 2 * ΣX 2 X 6 + b 3 * ΣX 5 X 6 + b 4 * ΣX 6 2
计算出搅拌速度Y2的值；
本炉脱硫剂量数据X6与Y1的值相等；所述n为每炉中成功脱硫的历史次数，所述b0、b1、b2、b3及b4为需要计算出的参数，所述∑X1、∑X2、∑X5、∑X6及所述∑Y2分别为对应量的n次历史数据之和，所述∑X1X2、∑X1X5、∑X2X6、∑X5X6、∑X1Y2、∑X2Y2、∑X5Y5及所述∑X6Y2分别为对应量的n次历史数据乘积之和。
3.根据权利要求2所述的铁水脱硫计算方法，其特征在于，所述根据本炉铁水温度数据X1、本炉铁水重量数据X2、本炉搅拌器使用次数数据X5及本炉搅拌速度数据X7，采用多元线性回归的方法计算出搅拌时间Y3包括：根据公式
Y3＝c0+c1*X1+c2*X2+c3*X5+c4*X7及公式
ΣY 3 = n * c 0 + c 1 * ΣX 1 + c 2 * ΣX 2 + c 3 * ΣX 5 + c 4 * ΣX 7 ΣX 1 Y 3 = c 0 * ΣX 1 + c 1 * ΣX 1 2 + c 2 * ΣX 1 X 2 + c 3 * ΣX 1 X 5 + c 4 * ΣX 1 X 7 ΣX 2 Y 3 = c 0 * ΣX 2 + c 1 * ΣX 1 X 2 + c 2 * ΣX 2 2 + c 3 * ΣX 2 X 5 + c 4 * ΣX 2 X 7 ΣX 5 Y 3 = c 0 * ΣX 5 + c 1 * ΣX 1 X 5 + c 2 * ΣX 2 X 5 + c 3 * ΣX 5 2 + c 4 * ΣX 5 X 7 ΣX 7 Y 3 = c 0 * ΣX 7 + c 1 * ΣX 1 X 7 + c 2 * ΣX 2 X 7 + c 3 * ΣX 5 X 7 + c 4 * ΣX 7 2
计算出搅拌时间Y3的值；
本炉搅拌速度数据X7与Y2的值相等；所述n为每炉中成功脱硫的历史次数，所述c0、c1、c2、c3及c4为需要计算出的参数，所述∑X1、∑X2、∑X5、∑X7及所述∑Y3分别为对应量的n次历史数据之和，所述∑X1X2、∑X1X5、∑X2X7、∑X5X7、∑X1Y3、∑X2Y3、∑X5Y3及所述∑X7Y3分别为对应量的n次历史数据乘积之和。
4.一种权利要求3所述的铁水脱硫计算方法所应用的系统，其特征在于，包括一级机、二级机及控制器；
所述一级机与所述控制器连接，且所述一级机用于显示脱硫过程中的各项参数以及用于输入控制指令修正本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值；
所述二级机与所述一级机连接，所述二级机用于存储KR脱硫剂加入量计算表、每炉中成功脱硫的历史数据及其它相关数据，并计算本次脱硫的脱硫剂量Y1、搅拌速度Y2及搅拌时间Y3的值；
所述控制器与所述一级机连接，所述控制器用于根据脱硫剂量Y1的值控制脱硫设备进行备料，根据脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值控制脱硫装置进行脱硫。
5.根据权利要求4所述的系统，其特征在于，所述二级机包括：
数据存储单元，用于存储KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据，并将KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据发送给数据处理单元；
数据处理单元，用于确定本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值。
6.根据权利要求5所述的系统，其特征在于，所述数据存储单元包括：
存储模块，用于存储KR脱硫剂加入量计算表及每炉中成功脱硫的历史数据；
发送模块，用于将存储模块中存储的所述KR脱硫剂加入量计算表及每炉中成功脱硫的所述历史数据发送给所述数据处理单元。
7.根据权利要求6所述的系统，其特征在于，所述数据处理单元包括：
接收模块，用于接收所述发送模块传来的所述KR脱硫剂加入量计算表及每炉中成功脱硫的所述历史数据，并将所述KR脱硫剂加入量计算表及所述历史数据发送给处理模块；
处理模块，根据所述KR脱硫剂加入量计算表及每炉中成功脱硫的所述历史数据查询出本次的脱硫剂量Y1或利用权利要求4所述的多元线性回归的方法分别计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值；并将本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值通过所述一级机发送给所述控制器，以及存储到所述存储模块中。
8.根据权利要求7所述的系统，其特征在于，所述处理模块包括：
查询模块，如果本次脱硫工作之前没有任何脱硫工作，则查询模块根据KR脱硫剂加入量表可以查询出本次脱硫的脱硫剂量Y1的值；并将本次脱硫的脱硫剂量Y1的值通过所述一级机发送给所述控制器，以及存储到所述存储模块中；
计算模块，如果本次脱硫工作之前有脱硫工作，利用权利要求4所述的多元线性回归的方法分别计算出本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值；并将本次脱硫的脱硫剂量Y1、搅拌速度Y2、搅拌时间Y3的值通过所述一级机发送给所述控制器，以及存储到所述存储模块中。
CN201310268262.1A 2013-06-28 2013-06-28 一种脱硫计算方法及其应用的系统 Expired - Fee Related CN103390098B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310268262.1A CN103390098B (zh) 2013-06-28 2013-06-28 一种脱硫计算方法及其应用的系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310268262.1A CN103390098B (zh) 2013-06-28 2013-06-28 一种脱硫计算方法及其应用的系统

Publications (2)

Publication Number Publication Date
CN103390098A CN103390098A (zh) 2013-11-13
CN103390098B true CN103390098B (zh) 2016-07-06

Family

ID=49534367

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310268262.1A Expired - Fee Related CN103390098B (zh) 2013-06-28 2013-06-28 一种脱硫计算方法及其应用的系统

Country Status (1)

Country Link
CN (1) CN103390098B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104815545B (zh) * 2015-04-23 2017-03-15 西安西热电站信息技术有限公司 一种实时计算湿法脱硫系统脱硫剂及副产品流量的方法
CN106503413B (zh) * 2015-08-31 2019-06-18 上海梅山钢铁股份有限公司 一种精确计算铁水脱硫用镁粉量的方法
CN106544471B (zh) * 2016-10-18 2018-05-29 北京首钢自动化信息技术有限公司 一种kr法脱硫冶炼过程搅拌转速的控制方法
CN111382393B (zh) * 2020-04-29 2024-01-09 鞍钢股份有限公司 一种铁水预处理脱硫粉剂计算方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2309006A1 (en) * 2009-09-21 2011-04-13 Korea Institute of Industrial Technology Desulfurizing agent and method for manufacturing the same
CN103074461A (zh) * 2012-10-17 2013-05-01 南京梅山冶金发展有限公司 一种适用于炼钢铁水脱硫模型调度装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2309006A1 (en) * 2009-09-21 2011-04-13 Korea Institute of Industrial Technology Desulfurizing agent and method for manufacturing the same
CN103074461A (zh) * 2012-10-17 2013-05-01 南京梅山冶金发展有限公司 一种适用于炼钢铁水脱硫模型调度装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
KR脱硫自动控制系统的应用;陈帅等;《钢铁研究》;20130610;第41卷(第3期);正文第45-49页 *
基于神经网络的铁水KR脱硫预报模型;但斌斌等;《计算机应用软件》;20110115;第28卷(第1期);正文第2.1-2.3节 *
炼钢厂铁水预处理一键脱硫控制系统;王熊;《科技传播》;20120308(第3期);正文第91-92页 *
铁水预处理炉外喷粉脱硫经验模型;杨辉合等;《世界钢铁》;20100315(第2期);正文第14页 *

Also Published As

Publication number Publication date
CN103390098A (zh) 2013-11-13

Similar Documents

Publication Publication Date Title
CN103390098B (zh) 一种脱硫计算方法及其应用的系统
CN103737878B (zh) 一种注塑缺陷在线修正方法及注塑机
MX2019014559A (es) Metodo y dispositivo de correccion de datos de mapa.
CN102766728A (zh) 钢包精炼炉精炼过程钢水硫含量实时预测的方法及装置
Wang et al. Prediction of calcium yield during calcium treatment process performed in steelmaking using neural network
CN104252565B (zh) 一种自动添加拼板中镀金引线的方法
CN103205531A (zh) 一种钢包自动增氮装置及方法
PL1791440T3 (pl) Sposób i instalacja do produkcji wyrobów cukierniczych
CN109280746A (zh) 一种kr脱硫搅拌的动态控制方法
CN108886935B (zh) 管道施肥系统施肥过程中施肥方案的确定方法、装置及系统
CN203745920U (zh) 一种磁悬液浓度智能控制装置
CN111832786A (zh) 一种基于滑动误差的光伏发电超短期功率预测方法和装置
CN110144426A (zh) 用于高炉上料配料组罐的控制方法、系统及可读存储介质
CN110760642B (zh) 一种铁水脱硫过程中脱硫剂投入量的控制方法
CN100507018C (zh) 一种确定rh精炼过程中吹氧量和冷材加入量的方法
CN102690923B (zh) 一种转炉副枪过程碳含量预报方法
CN105070164B (zh) 面向虚拟对象单回路控制系统实验装置
CN205903881U (zh) 一种化肥助剂的连续生产装置
CN107885084A (zh) 基于rbf与adhdp的天然气吸收塔脱硫过程控制方法
CN206078980U (zh) 一种饲料混合系统
CN103744402B (zh) 汽油管道调合调度批次自动切换方法
Kehr Project QI-Digital-Quo vadis Use Case H2 Refueling Station?
CN204689632U (zh) 一种液态物料自动计量装罐装置
CN202968595U (zh) 转炉煤气回收双限控制系统
CN110322256A (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

Granted publication date: 20160706

Termination date: 20170628

CF01 Termination of patent right due to non-payment of annual fee