汽轮机转子瞬态应力的检测方法与流程

本发明涉及汽轮机领域，尤其涉及一种汽轮机转子瞬态应力的检测方法。

背景技术：

汽轮机组中，转子的工作状况对汽轮机组运行安全以及寿命管理具有很重要的意义。而在汽轮机组的启动、停机等过程中，汽轮机组的主轴、叶片和叶轮等部件会因高速旋转、强扭矩和高温蒸汽而产生巨大的瞬态应力，威胁转子的正常工作。因此，对转子瞬态应力的分析和预测成为汽轮机组寿命管理的关键。但是，转子恶劣的工作环境导致其瞬态应力难以直接测量获得。

目前，对于汽轮机转子寿命的研究中，通常，以预测转子瞬态应力为研究核心，以有限元分析为基础，开发各种有限元分析软件，在绝大多数应用场景下具有较好的预测精度和建模复杂度。然而对大型汽轮机组多参数共同作用的复杂情况的预测效果不佳，建模过程繁琐耗时，且得到的结果准确度低，计算效率也不高。

技术实现要素：

本发明所要解决的技术问题是针对现有技术的不足，提供一种汽轮机转子瞬态应力的检测方法及一种存储介质。

本发明解决上述技术问题的技术方案如下：

一种汽轮机转子瞬态应力的检测方法，包括：

获取汽轮机组内影响汽轮机转子瞬态应力的特征参数；

建立瞬态应力预测模型，根据所述特征参数对所述瞬态应力预测模型进行训练，其中，所述瞬态应力预测模型包括两层结构，第一层结构由n个互不相同的机器学习模型组成，所述第一层结构中的n个所述机器学习模型的输出与第二层结构中的机器学习模型的输入连接，n≥2；

根据训练后的所述瞬态应力预测模型对汽轮机转子的瞬态应力进行检测。

本发明的有益效果是：本发明提供的检测方法，通过根据影响汽轮机转子瞬态应力的特征参数对具有两层结构的瞬态应力预测模型进行训练，提高了瞬态应力预测模型对复杂情况的预测能力，保证了充分学习数据间的信息，提高了多参数复杂情况下转子瞬态应力预测的精度，且不需要进行复杂建模。

本发明解决上述技术问题的另一种技术方案如下：

一种存储介质，所述存储介质中存储有指令，当计算机读取所述指令时，使所述计算机执行如上述技术方案所述的方法。

本发明附加的方面的优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明实践了解到。

附图说明

图1为本发明汽轮机转子瞬态应力的检测方法的实施例提供的流程示意图；

图2为本发明汽轮机转子瞬态应力的检测方法的其他实施例提供的瞬态应力预测模型结构示意图；

图3为本发明汽轮机转子瞬态应力的检测方法的其他实施例提供的数据集划分结构示意图。

具体实施方式

以下结合附图对本发明的原理和特征进行描述，所举实施例只用于解释本发明，并非用于限定本发明的范围。

如图1所示，为本发明汽轮机转子瞬态应力的检测方法的实施例提供的流程示意图，该检测方法包括：

s1，获取汽轮机组内影响汽轮机转子瞬态应力的特征参数。

需要说明的是，由于汽轮机组所涉及的特征参数较多，不同特征与转子应力的相关性存在差异，如果将相关性较低的特征加入模型中，相当于将噪声加入模型，不仅会降低模型精度，还会增加模型规模，加大计算量，造成过拟合。因此，特征参数应选取与转子瞬态应力的相关性较大或重要程度较高的汽轮机组的运行参数，例如，可以为主蒸汽压力、主蒸汽温度、机组转速、机组负荷、汽缸内表面金属温度、主蒸汽温度变化率等，这些参数均可在机组运行中直接测量获得。

s2，建立瞬态应力预测模型，根据特征参数对瞬态应力预测模型进行训练。

瞬态应力预测模型包括两层结构，第一层结构由n个互不相同的机器学习模型组成，第一层结构中的n个机器学习模型的输出与第二层结构中的机器学习模型的输入连接，n≥2。

需要说明的是，机器学习模型可以为lasso回归模型、弹性网络模型、随机森林模型、线性回归模型、循环神经网络rnn、长短时记忆网络lstm等。

其中，lasso回归模型、弹性网络模型、随机森林模型为较优的选择，这三个模型的预测精度较高。

例如，当n＝2时，第一层结构由2个互不相同的机器学习模型组成，可以为lasso回归模型和弹性网络模型，第二层结构可以由lasso回归模型组成。

又例如，当n＝3时，如图2所示，第一层结构1由3个互不相同的机器学习模型组成，可以为lasso回归模型、弹性网络模型和随机森林模型，第二层结构2可以由lasso回归模型组成。

需要说明的是，由于瞬态应力预测模型为两层结构，因此，在对其进行训练时，可以首先分别对第一层结构中的机器学习模型进行训练，分别得到训练结果后，将得到的这些训练结果与全部的训练集合并，作为新的训练集，再对第二层的机器学习模型进行训练，这样可以显著提高瞬态应力预测模型对复杂情况的预测精度。

应理解，在进行训练前，需要获取训练集。训练集可以根据特征参数获取。

例如，假设主蒸汽温度是与转子的瞬态压力相关度很高的特征参数，那么可以在预设的时间段内，隔预设间隔采集主蒸汽温度，就得到了一个数据集。可以根据以下方法对数据集进行划分：

根据交叉验证法将数据集分为n份，将其中1份作为测试集，其余作为训练集；

重复划分n次，得到n个划分后的数据集，每个数据集中都有n-1份训练集。

例如，如图3所示，假设n＝3，即第一层结构中有3个机器学习模型，那么可以使用3折交叉验证来划分数据集，可以在一定程度上降低过拟合的风险和减小方差。

首先，使用不重复抽样将数据集随机划分为3份，挑选其中的1份数据作为测试集3，其余2份作为训练集4，一共重复迭代3次，那么如图3所示，就能够得到3个不同的数据集，每个数据集中都有1份测试集3和2份训练集4，且2份训练集4的组合都是不同的，这样保证了每份子集都有机会作为测试集3和训练集4，提高预测的精度。

s3，根据训练后的瞬态应力预测模型对汽轮机转子的瞬态应力进行检测。

应理解，在对瞬态应力预测模型进行训练后，直接将检测到的汽轮机组的相关的运行数据输入到瞬态应力预测模型，就能快速的得到瞬态应力的分析结果，相比于传统方法，速度提高很多。

本实施例提供的检测方法，通过根据影响汽轮机转子瞬态应力的特征参数对具有两层结构的瞬态应力预测模型进行训练，提高了瞬态应力预测模型对复杂情况的预测能力，保证了充分学习数据间的信息，提高了多参数复杂情况下转子瞬态应力预测的精度，且不需要进行复杂建模。

可选地，在一些实施例中，还包括：

根据特征参数对瞬态应力预测模型进行预测。

可选地，在一些实施例中，根据特征参数对瞬态应力预测模型进行预测，具体包括：

根据交叉验证法将特征参数的数据集分为n份，将其中1份作为测试集，其余作为训练集；

重复划分n次，得到n个划分后的测试集；

将n个测试集一对一地输入到第一层结构中的n个机器学习模型中进行预测，得到n个初步预测结果；

对n个初步预测结果取均值，将均值和n个测试集输入到第二层结构中的机器学习模型中进行预测，得到预测结果。

可选地，在一些实施例中，获取汽轮机组内影响汽轮机转子瞬态应力的特征参数，具体包括：

获取汽轮机组运行时的全部运行参数；

根据预设特征选择模型对全部运行参数进行训练，根据训练结果确定全部运行参数的重要性，将重要性大于预设阈值的运行参数作为影响汽轮机转子瞬态应力的特征参数。

预设阈值可以根据实际需求设置。

优选地，还可以在计算完全部运行参数的重要性后，按重要性的大小进行排序，将排名在前x名的运行参数作为特征参数，x的值可以根据实际需求设置，这样也能够挑选出对汽轮机转子瞬态应力影响较大的参数，提高预测的精度。

应理解，用于特征重要性分析的模型可以采用随机森林模型、支持向量机、arima、人工神经网络、贝叶斯网络模型等，较优的模型是随机森林模型。

需要说明的是，目前在国内外相关研究领域以及国际相关工程应用中，对于建模特征参数的选取，仍然以物理模型分析为主，以转子设计的专家知识为辅，建模过程十分复杂，且涉及的参数过多。工况复杂时，有限元分析无法适应复杂的边界条件变化，从而无法得出问题的解。而本实施例利用特征重要性分析方法，通过分析数据对模型精度的影响程度，直接挖掘数据间的关联程度，挑选出适合建模的特征参数，有效避免了上述问题。

可选地，在一些实施例中，根据预设特征选择模型对全部运行参数进行训练，根据训练结果确定全部运行参数的重要性，具体包括：

将任一运行参数输入到预设特征选择模型中进行训练，得到第一训练结果；

根据第一训练结果与汽轮机转子瞬态应力进行误差计算，得到第一误差；

向运行参数中加入随机噪声；

将加入随机噪声后的运行参数重新输入到预设特征选择模型中进行训练，得到第二训练结果；

根据第二训练结果与汽轮机转子瞬态应力进行误差计算，得到第二误差；

根据第一误差、第二误差和预设重要性计算公式，计算得到运行参数的重要性；

重复计算其他运行参数的重要性，直到全部运行参数的重要性计算完成。

可选地，在一些实施例中，根据以下公式计算运行参数的重要性：

importace＝(|e1-e2|)/n

其中，importace为重要性，e1为第一误差，e2为第二误差，n为预设特征选择模型中决策树的数量。

可选地，在一些实施例中，预设特征选择模型为随机森林模型。

可选地，在一些实施例中，第一层结构由3个互不相同的机器学习模型组成，分别是lasso回归模型、弹性网络模型和随机森林模型。

可选地，在一些实施例中，第二层结构中的机器学习模型为lasso回归模型。

可以理解，在一些实施例中，可以包含如上述各实施例中的部分或全部步骤。

在本发明的其他实施例中，还提供一种存储介质，该存储介质中存储有指令，当计算机读取该指令时，使计算机执行如上述实施例中任一项所述的方法。

读者应理解，在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。

作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本发明实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。

以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。