一种基于Matlab软件计算石灰石浆液密度最优值算法模型的方法

技术领域

本发明属于湿法脱硫技术领域，具体属于一种基于Matlab软件计算石灰石浆液密度最优值算法模型的方法。

背景技术

随着全球经济的高速发展和工业化的不断推进，大气中二氧化硫排放量与日俱增，造成降水pH值下降，局部地方甚至形成酸雨，对人体健康和大气环境带来很大影响。为降低空气中二氧化硫排放量，提高空气质量，防止环境进一步恶化，燃煤发电厂必须安装烟气脱硫装置。。因此大力发展环保设备，清洁我们的生态环境也是我们的一个重要任务。目前，燃煤电厂应用最广泛的是石灰石/石膏湿法脱硫。

当前，湿法脱硫石灰石浆液配置方法主要有三种，分别是：外购石灰石粉，直接加水制浆；外购石灰石，利用湿式球磨机制备浆液及外购石灰石，利用干式球磨机制粉后加水制浆。

外购石灰石粉价格高，干磨系统投资大、电耗高，相比之下，湿磨系统投资小、成本低，占据火电厂脱硫石灰石浆液配置主流地位。

系统运行时，为控制装置运行所需的石灰石量，采用变频式称重皮带给料机控制进入湿式球磨机的石灰石量，根据进入磨机的石灰石量通过磨机进口补充水管道上流量调节阀控制进入磨机的过滤水量，磨制好的浆液由磨机出口自流到磨机浆液再循环箱，同时通过过滤水地坑泵将石膏滤液水打入石灰石浆液再循环箱稀释石灰石浆液，进入石灰石浆液再循环箱的滤液水流量大小通过电动调节阀控制，再经过磨机石灰石浆液再循环泵（1运1备）送入石灰石旋流器，旋流器溢流浆液达到满足系统出力及细度要求后，进入石灰石浆液罐，而底流浆液返回磨机进料口，参与湿磨机的重新磨制。非正常工况下，旋流站底流和溢流均返回磨机再循环浆液箱。储存在石灰石浆液罐的符合要求浆液密度的合格浆液，通过石灰石浆液泵送入吸收塔的石灰石供浆调节系统，去完成吸收塔内浆液PH的调整。

湿式球磨机制浆，石灰石与水的比例一般固定不变，制备的石灰石浆液通过过滤水地坑泵将石膏滤液水打到磨机再循环箱进行稀释后，再由石灰石浆液再循环泵打进石灰石浆液箱，过滤水至磨机再循环箱管道安装流量调节门，调节门仅根据再循环箱液位进行流量大小调节，若要改变石灰石浆液密度，就必须通过人为干预，在原烟气入口SO2浓度、机组负荷、原烟气流量等急剧升高，石灰石浆液供给流量达到设计值，而净烟气SO2浓度持续升高无法控制的情况下，再通过加大石灰石浆液密度的方式来控制净烟气SO2浓度时间根本来不及，造成环保指标超标。

发明内容

本发明的目的是提供一种基于Matlab软件计算石灰石浆液密度最优值算法模型的方法，要解决现有技术中石灰石浆液密度低而导致的净烟气SO

为实现上述目的，本发明采用如下技术方案：

一种基于Matlab软件计算石灰石浆液密度最优值算法模型的方法，其特征在于：包括以下步骤：

步骤一：采集历史数据：从DCS系统中导出原烟气SO

步骤二：模型训练：通过Matlab软件的BP神经网络模型作为实用模型进行训练，得出模型实例，利用模型实例对实际运行工况进行预测对比通过改变模型参数，反复训练、验证，直到模型实例与实际工况误差在允许范围内为止，得到最优值算法模型；

步骤三：嵌入DCS系统：从厂家获得DCS系统权限，将最优值算法模型嵌入DCS系统，根据某同一时刻的原烟气SO

步骤四：调整石灰石浆液密度：将实测石灰石浆液密度引入DCS系统，DCS系统将根据最优值算法模型的计算值与实测值的偏差，利用PID（微积分调节模块）调节磨机再循环箱过滤水调节阀的开度，控制再循环箱过滤水流量大小，调整石灰石浆液密度，最终达到密度实测值与公式计算值误差在允许范围。

进一步优选地，所述步骤二具体包括以下步骤：

步骤1，导入历史数据：在DCS系统历史数据中选取所需建模的数据，并导出为 csv格式，将导出的数据上传matlab 软件中；

步骤2，软件操作：依据需要进行分隔符、导入范围、及 str 类型数的选取，之后选择矩阵方式；

步骤3，数据预处理：神经网络以每列为一个样本，每行为单一变量，因此需对导入数据进行转置，在命令行窗口输入如下命令对数据进行输入和输出的转置：

命令：data=Assaydatabulian2';

取 2 到 14 行作为输入x，自动输出目标变量y，在命令行窗口输入如下命令对输入数据x进行选择：

命令：x_raw=data(2:14,:);y_raw=data(1,:);

步骤4，对输入数据x、目标变量y进行归一化处理：

在命令行窗口输入如下命令，对输入数据x、目标变量y进行归一化处理：

命令：[x,inputps]=mapminmax(x_raw);

[y,outputps]=mapminmax(y_raw);

归一化计算方法为：2*(x-xmin)/(xmax-xmin)+(-1)， xman表达x的最大数据，xmin表达x的最小数据；

步骤5，y值训练及创建神经网路模型：点击 Import 按钮，出现输入导出窗口，选择 x 作为输入数据，导入目标变量 y，软件自动进行y值训练，同时创建神经网路模型；

步骤6，对目标变量y进行反归一化处理：经神经网路模型计算后，需对输出数据目标变量y进行反归一化，以匹配系统实际输出在命令行窗口输入如下命令，对目标变量y进行反归一化：

命令:y_pre=mapminmax('reverse',y_tem,outputps);

y_pre 即为石灰石浆液最优值算法模型。

进一步地，所述步骤5还包括以下步骤：神经网络训练及结果查看，根据训练结果中实线实际计算结果与图中虚线预测结果是否接近判断训练结果是否接近理想效果，若两者夹角大于10°则进行重新训练直至达到理想效果。

此外，所述步骤四中DCS系统上设置绿、黄、红三种密度偏差指示灯，其中，

绿色表示偏差在±10kg/m

黄色表示偏差在±15kg/m

红色表示偏差在±20kg/m

更加优选地，绿色指示灯亮时：石灰石浆液箱正常工作；

黄色指示灯亮时：PID自动调节磨机再循环箱过滤水调节阀的开度，对石灰石浆液进行自动配比，直至绿色指示灯亮；

红色指示灯亮时：PID调节磨机再循环箱过滤水调节阀的开度对石灰石浆液进行自动配比在短时间内达不到“绿色”偏差范围，人工干预，手动配比直至绿色指示灯亮。

与现有技术相比本发明具有以下特点和有益效果：

本发明在石灰石浆液制备过程中可根据原烟气SO

附图说明

图1为本发明涉及的神经网络训练结果图示。

具体实施方式

为使本发明实现的技术手段、创新特征、达成目的与功效易于明白了解，下面对本发明进一步说明。

在此记载的实施例为本发明的特定的具体实施方式，用于说明本发明的构思，均是解释性和示例性的，不应解释为对本发明实施方式及本发明范围的限制。除在此记载的实施例外，本领域技术人员还能够基于本申请权利要求书和说明书所公开的内容采用显而易见的其它技术方案，这些技术方案包括采用对在此记载的实施例的做出任何显而易见的替换和修改的技术方案。

一种基于Matlab软件计算石灰石浆液密度最优值算法模型的方法，其特征在于：包括以下步骤：

步骤一：采集历史数据：从DCS系统中导出原烟气SO

步骤二：模型训练：通过Matlab软件的BP神经网络模型作为实用模型进行训练，得出模型实例，利用模型实例对实际运行工况进行预测对比通过改变模型参数，反复训练、验证，直到模型实例与实际工况误差在允许范围内为止，得到最优值算法模型；

具体包括以下步骤：

步骤1，导入历史数据：在DCS系统历史数据中选取所需建模的数据，并导出为 csv格式，将导出的数据上传matlab 软件中；

步骤2，软件操作：依据需要进行分隔符、导入范围、及 str 类型数的选取，之后选择矩阵方式；

步骤3，数据预处理：神经网络以每列为一个样本，每行为单一变量，因此需对导入数据进行转置，在命令行窗口输入如下命令对数据进行输入和输出的转置：

命令：data=Assaydatabulian2';

取 2 到 14 行作为输入x，自动输出目标变量y，在命令行窗口输入如下命令对输入数据x进行选择：

命令：x_raw=data(2:14,:);y_raw=data(1,:);

步骤4，对输入数据x、目标变量y进行归一化处理：

在命令行窗口输入如下命令，对输入数据x、目标变量y进行归一化处理：

命令：[x,inputps]=mapminmax(x_raw);

[y,outputps]=mapminmax(y_raw);

归一化计算方法为：2*(x-xmin)/(xmax-xmin)+(-1)， xman表达x的最大数据，xmin表达x的最小数据；

步骤5，y值训练及创建神经网路模型：点击 Import 按钮，出现输入导出窗口，选择 x 作为输入数据，导入目标变量 y，软件自动进行y值训练，同时创建神经网路模型，神经网络训练及结果查看，根据训练结果中实线实际计算结果与图中虚线预测结果是否接近判断训练结果是否接近理想效果，若两者夹角大于10°，利用matlab软件进行神经网络非线性拟合，经过反复训练、验证，确定模型实例，直至达到理想效果。

步骤6，对目标变量y进行反归一化处理：经神经网路模型计算后，需对输出数据目标变量y进行反归一化，以匹配系统实际输出在命令行窗口输入如下命令，对目标变量y进行反归一化：

命令:y_pre=mapminmax('reverse',y_tem,outputps);

y_pre 即为石灰石浆液最优值算法模型。

步骤三：嵌入DCS系统：从厂家获得DCS系统权限，将最优值算法模型嵌入DCS系统，根据某同一时刻的原烟气SO

步骤四：调整石灰石浆液密度：将实测石灰石浆液密度引入DCS系统，DCS系统将根据最优值算法模型的计算值与实测值的偏差，利用PID调节磨机再循环箱过滤水调节阀的开度，控制再循环箱过滤水流量大小，调整石灰石浆液密度，最终达到密度实测值与公式计算值误差在允许范围。

其中，DCS系统上设置绿、黄、红三种密度偏差指示灯，其中，

绿色表示偏差在±10kg/m

黄色表示偏差在±15kg/m

红色表示偏差在±20kg/m

绿色指示灯亮时：石灰石浆液箱正常工作；

黄色指示灯亮时：PID自动调节磨机再循环箱过滤水调节阀的开度，对石灰石浆液进行自动配比，直至绿色指示灯亮；

红色指示灯亮时：PID调节磨机再循环箱过滤水调节阀的开度对石灰石浆液进行自动配比在短时间内达不到“绿色”偏差范围，人工干预，手动配比直至绿色指示灯亮。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。