具体实施方式

下面结合附图对本申请作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本申请的保护范围。

本发明采用如下的技术方案。

如图1，一种基于支持向量回归机的水泥搅拌机能耗建模方法包括：

步骤1，采集水泥搅拌机的运行过程状态数据和稳定性数据。

具体地，步骤1中，运行过程状态数据包括：工作电流、阻抗、轴功率、桨叶排液量、压头、桨叶直径、搅拌转速、进料量、出料量、骨料平均粒径；稳定性数据包括：工作电压、工作频率。

水泥搅拌机的原始数据包括：水泥搅拌机运行过程中产生的主要电气量，包括但不限于电压、电流、频率、阻抗；核心生产工艺参数，包括但不限于轴功率、桨叶排液量、压头、桨叶直径、搅拌转速；实际操作的参数，包括但不限于进料量、出料量、骨料平均粒径。

其中，对于主要电气量参数而言，搅拌机由于接入电网，其电压、频率需要保持在标准值，维系产品和电网的稳定性，而阻抗和电流可以体现搅拌机的运行效率和过程，搅拌中骨料的材质、粒径和数量会影响搅拌机的阻抗，从而在电流的变化上有所体现，在启动中电流较大随搅拌过程进行电流不断减小趋于稳定，同时，当原料量增大负荷增大，工作电流也会上升；

对于核心生产工艺参数而言，轴功率、桨叶排液量、压头、桨叶直径、搅拌转速，这些参数互相存在耦合关系，都无法单独用来描述搅拌机的工作状态，桨叶排液量与桨叶直径、转速有关，搅拌消耗的轴功率和流体比重、桨叶转速等有关，一定功率下桨叶排液量和压头又可以通过改变桨叶直径和转速来调节等，因此在考虑这些因素时要综合考虑；

对于实际操作的参数而言，进料量、出料量、骨料平均粒径，这些因素都是可以方便人为调节，在运送搅拌原料时，可以调节搅拌的进料量来调节搅拌的负荷，改变骨料的粒径改变搅拌机运行过程的阻抗，在工作电流中有所体现，而因为粒径大多不规则，在获取数据时选择得到粒径的平均值作为数据，出料量可以体现搅拌机的运行效率。

步骤2，对运行过程状态数据进行预处理，得到参数矩阵A。

具体地，步骤2包括：

步骤2.1，从运行过程状态数据中提取各数据的特征值。

进一步，步骤2.1中，对于搅拌转速、进料量、出料量、桨叶直径，均取实际值作为特征值；对于轴功率、压头、桨叶直径、骨料平均粒径，均取平均值作为特征值；

根据工作电流，提取搅拌机的启动峰值电流和达到电流波动不超过5％的稳定阶段所用的电流稳定时间作为特征值；

提取桨叶排液量以及桨叶排液量的平均变化率作为特征值。

提取数据特征值时，首先工作电压和频率变化不大用于监控水泥搅拌机的稳定性，提取一段工作时间内的平均值作为特征；对工作电流的变化提取搅拌机启动时的峰值电流和达到电流波动不超过5％的稳定阶段所用时间来反映搅拌的负荷和搅拌的程度；桨叶直径、进料量、出料量用实际值，另外提取桨叶排液量及其平均变化率，对其他数据提取平均值作数据。

步骤2.2，对特征值数据，采用Z-score标准化方法，以如下关系式进行归一化处理：

式中，x表示某一项数据的特征值，表示某一项数据的平均值，σ表示某一项数据的标准差，x^*为归一化处理后得到的数据，其满足正态分布。

进一步，步骤2.2中，对于以实际值作为特征值的各数据，即搅拌转速、进料量、出料量不进行归一化处理；而其他数据间量纲不同，同一数据取值范围和变化范围相差过大，影响数据分析的结果，且因为数据间存在耦合情况，因此需要进行归一化处理。

步骤2.3，将归一化处理后的各数据特征值进行向量化处理，得到参数矩阵A。

进一步，步骤2.3中，向量化，即将归一化处理后的各数据特征值，按进料量、出料量、启动峰值电流、电流稳定时间、轴功率、桨叶排液量、压头、桨叶直径、搅拌转速、桨叶排液量变化率和骨料平均粒径的顺序依次放入参数矩阵A中。

步骤3，根据参数矩阵A，使用奇异值分解法得到新参数矩阵A^*；

具体地，步骤3包括：

步骤3.1，根据参数矩阵A构建第一方阵AA^T和第二方阵A^TA；

步骤3.2，求解第一方阵AA^T的第一特征值u_i和第一特征向量、第二方阵A^TA的第二特征值v_i和第二特征向量；

步骤3.3，利用第一特征向量构建维数为m×m的左奇异矩阵U，利用第二特征向量构建维数为n×n的右奇异矩阵V；

步骤3.4，以如下关系式计算奇异值矩阵的对角元δ_i，所得对角元δ_i按从大到小排列组成奇异值矩阵：

步骤3.5，选取k值，使得奇异值矩阵中前k列奇异值之和与全部奇异值之和之比大于等于95％；在右奇异矩阵V中取前k行组成新的右矩阵V^*；

值得注意的是，本发明优选实施例中，选取的k值使得奇异值矩阵中前k列奇异值之和与全部奇异值之和之比大于等于95％，这一比例数值的选取是一种非限制性的较优选择。

步骤3.6，以如下关系式计算得到新参数矩阵A^*：

A^*＝A(V^*)^T

新参数矩阵A^*是参数矩阵A的降维后结果。

步骤4，将稳定性数据加入新参数矩阵A^*，得到数据集B并作为自变量，将对应的水泥搅拌机消耗的总能耗值作为因变量，使用高斯径向基函数作为核函数，构建水泥搅拌机能耗模型。

具体地，步骤4包括：

步骤4.1，选择径向基函数作为核函数，构建支持向量回归机的模型，其中，所述径向基函数为高斯核函数，满足如下关系式：

式中，x和x′表示两个不同的向量，τ表示高斯滤波器带宽，决定着平滑程度，影响了核函数的局部作用范围。

步骤4.2，在新参数矩阵A^*中，加入工作电压和工作频率的特征值数据，得到数据集B作为自变量，即水泥搅拌机能耗模型的输入数据集；以水泥搅拌机的总能耗作为水泥搅拌机能耗模型的输出数据集；

步骤4.3，将输入数据集划分为训练集和测试集。

本发明优选实施例中，输入数据集按照7：3的比例划分为训练集和测试集，所采用的划分比例是一种非限制性的较优选择。

步骤4.4，使用训练集对水泥搅拌机能耗模型进行训练；

步骤4.5，使用测试集对训练好的模型进行测试，使得模型输出的总能耗的预测值和真实值之间的误差不大于5％。

进一步，步骤4.2中，提取一段工作时间内工作电压和工作频率的平均值作为特征值。

进一步，步骤4.4中，基于K折交叉验证法，将训练集分成K折；选其中一折数据作为验证集，其他K-1折数据作为训练集，进行训练和验证；然后，换另一折数据作为验证集，其他K-1折数据作为训练集；直到每一折数据都作过验证集后，得到训练好的模型。

实施例1。

按照图2所示的流程实施本发明提出的基于支持向量回归机的水泥搅拌机能耗建模方法，为了验证该建模方法得到的模型的有效性，收集部分水泥搅拌机的数据在Python中进行了建模，绘制水泥搅拌机能耗的预测值和真实值的拟合情况曲线图。图3给出了进行降维前水泥搅拌机能耗的预测值和真实值的拟合情况。图4给出了加入降维后水泥搅拌机能耗的预测值和真实值的拟合情况，由此可见对水泥搅拌机能耗情况预测精确。同时用均方根误差作为评价模型精确程度的标准，在加入降维处理前，模型的均方根误差为0.4737，在加入降维处理后，模型的均方根误差下降为0.4507，且模型搭建的速度更快，由此可见加入降维处理能有效提高模型的运行速度和预测精度。

本发明的有益效果在于，与现有技术相比：

1、能够有效分析典型工业负荷水泥搅拌机的能耗问题，方便分析各种参数对水泥搅拌机能耗产生的影响，方便电网实现工业负荷的优化调度问题。

2、选用奇异值分解避免了使用特征值分解时在大量数据特征情况下需要计算原始数据协方差矩阵的麻烦，降低了模型的复杂度，提高模型运算效率，方便在线式监控能耗情况。

3、使用支持向量回归机的方法，具有较好的鲁棒性，对未知数据拥有很好的泛化能力，尤其是在数据量较少的小样本情况下，相较于其他传统机器学习办法，具有更优的性能。

4、建立一套反映其能耗表征特性的模型，可以挖掘其可调潜力，提高电网对能量资源的综合利用和对该工业负荷的运行效率，引导用户方参与电网实时供需平衡。

本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述，但是本领域技术人员应该理解，以上实施示例仅为本发明的优选实施方案，详尽的说明只是为了帮助读者更好地理解本发明精神，而并非对本发明保护范围的限制，相反，任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。