一种风电机组的疲劳寿命评估方法、系统、设备及介质与流程

1.本发明涉及风力发电机组发电技术领域，具体涉及一种风电机组的疲劳寿命评估方法、系统、设备及介质。


背景技术：

2.随着风电行业的发展，机组功率等级和尺寸的加大对风电机组整体性能和机组安全性的要求也越来越高，其中风电机组的疲劳寿命就是机组重要的安全性指标。目前风电机组的疲劳寿命通过机器学习方法预测得到或通过仿真计算得到。
3.现有技术中申请号为cn2020102580882,名称为《一种风电叶片剩余疲劳寿命评估方法》的专利申请文件公开内容如图1所示，主要应用于在役机组叶片技改项目，通过评估技改前后的疲劳损伤验证技改后延长了叶片的剩余疲劳寿命。主要是针对已知叶片疲劳损伤过大后通过优化控制策略对叶片进行控制改造，是一种验证方法，对于原叶片已经造成的疲劳损伤并未明确方法；并且控制策略的改变可能增加机组其他部件的疲劳损伤。该发明使用的是仿真计算与有限单元法相结合的方法计算疲劳损伤，计算量过大才能完全覆盖机组实际状态。
4.现有技术中申请号为cn2020100786467,名称为《一种风电机组疲劳寿命的快速预测的方法》的专利申请文件公开内容如图2所示，主要是应用机器学习方法预测风电机组疲劳寿命，通过对雨流计算的疲劳损伤数据集训练、验证与测试，得到疲劳损伤的机器模型。该方法虽然无需耗费大量的计算资源及时间，但受限于机器计算模型与特征量的选取，预测结果随机性比较强。
5.近年来随着装机容量的不断增大，风电机组关键部件，如叶片、塔筒、传动链等所承受的载荷范围越来越大，机组倒塔、叶片折断等安全事故率也逐渐增加， 风力发电机组的寿命与安全问题日益凸显。风电机组疲劳载荷的数据主要来源于试验测试和仿真计算两种方法，但因为机组的设计寿命一般为20年，而测试和仿真计算采样时间短，通过此两种方法所获取的载荷数据直接用于疲劳分析会出现较大的分析误差和无法预计的困难。


技术实现要素：

6.为了解决现有技术中由于试验测试和仿真计算的采样时间短以及外部环境条件限制导致获取的载荷数据用于疲劳分析会出现较大的分析误差和困难的问题，本发明提出了一种风电机组的疲劳寿命评估方法，包括：基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值；基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱；基于所述计算周期的疲劳载荷谱结合s-n曲线进行评估得到在役风电机组的剩余疲劳寿命。
7.优选的，所述基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值，包括：基于所述在役风电机组的运行状态进行区间划分，得到各运行状态下的各风速区间；基于所述各运行状态下的各风速区间对所述风电场外部气象参数进行采样，得到各运行状态下各风速区间内的载荷样本；基于所述各运行状态下各风速区间内的载荷样本进行仿真计算，得到各运行状态下各风速区间内第一单位时长的载荷时序值；其中，所述在役风电机组的运行状态包括：正常发电状态、停机状态和启停机制动状态。
8.优选的，所述基于所述在役风电机组的运行状态进行区间划分，得到各风速区间，包括：当在役风电机组的运行状态为正常发电状态时，基于风电机组的切入风速以设定间隔等间隔递增至切出风速进行区间划分，得到正常发电状态下的各风速区间；当在役风电机组的运行状态为停机状态时，基于初始风速以设定间隔等间隔递增至目标风速进行区间划分，得到停机状态下的各风速区间；当在役风电机组的运行状态为启停机制动状态时，基于初始风速以设定间隔等间隔递增至目标风速进行区间划分，得到启停机制动状态下的各风速区间。
9.优选的，所述基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱，包括：基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值进行雨流计算，得到载荷幅值、载荷幅值极大值和载荷幅值极小值；基于在役风电机组计算周期内第二单位时长的平均风速频率分布结合各幅值区间，得到各风速区间内第一单位时长内各幅值区间的循环作用次数；基于所述各风速区间内第一单位时长内各幅值区间的循环作用次数结合所述载荷幅值得到第一单位时长疲劳载荷谱；基于在役风电机组计算周期内第二单位时长的平均风速频率分布计算得到各风速区间内计算周期内各幅值区间的循环作用次数；基于所述载荷幅值结合所述各风速区间内计算周期内各幅值区间的循环作用次数得到计算周期疲劳载荷谱；其中，所述各幅值区间是基于载荷幅值极大值和载荷幅值极小值进行区间划分得到的。
10.优选的，所述载荷幅值极大值按下述方法得到，包括：基于在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值进行雨流计算，得到各运行状态下各风速区间内的载荷幅值、载荷均值以及相应的循环作用次数和雨流矩阵；基于所述雨流矩阵中的载荷幅值做双参数威尔布分布，得到载荷幅值的概率密度函数；
基于所述载荷幅值的概率密度函数进行积分得到载荷幅值累积频率；基于所述载荷幅值累积频率得到载荷幅值的极大值。
11.优选的，所述载荷幅值极小值按下述方法得到，包括：基于所述载荷幅值和所述相应的循环作用次数得到s-n曲线，并基于所述载荷幅值划分为各个等级的载荷；基于所述s-n曲线结合线性损伤理论，计算得到各个等级的载荷产生的疲劳损伤和各个等级的载荷产生的疲劳损伤密度；基于所述各个等级的载荷产生的疲劳损伤和各个等级的载荷产生的疲劳损伤密度确定载荷幅值极小值。
12.优选的，所述各个等级的载荷产生的疲劳损伤的计算式如下所示：
13.式中，为载荷谱第k级载荷产生的疲劳损伤；为第k级载荷对应的循环次数；为第k级载荷；为与材料性质相关的拟合参数；为第k级载荷m次方。
14.优选的，所述各个等级的载荷产生的疲劳损伤密度的计算式如下所示：
15.式中，为载荷谱第k级载荷产生的疲劳损伤；为第k级载荷对应的循环次数；为第k级载荷产生的疲劳损伤密度；为第k级载荷；为第k级载荷m次方，为级数。
16.优选的，所述计算周期内各幅值区间的循环作用次数按下式计算：
17.式中，为载荷幅值对应的循环次数；为i区间风速的分布频率；为发电状态下i风速区间中载荷幅值对应的循环次数；为停机状态下i风速区间中载荷幅值对应的循环次数；为启停机制动状态下i风速区间中载荷幅值对应的循环次数；为正常发电状态的年频率；a为计算周期；b为第三单位时长；为风速区间i的风速；为区间序号。
18.优选的，所述运行状态的年频率是基于数据采集与监视控制系统的数据统计得到的。
19.优选的，所述风电场外部气象参数包括在役风电机组计算周期内的第二单位时长平均风速频率分布、第二单位时长湍流强度、平均温度、平均湿度、平均气压、平均风剪切、平均风转向和平均入流角。
20.基于同一发明构思，本发明还提出了一种风电机组的疲劳寿命评估系统，包括：仿真计算模块，用于基于获取的风电场外部气象参数和在役风电机组实际运行状
态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值；统计外推模块，用于基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱；寿命评估模块，用于基于所述计算周期的疲劳载荷谱结合s-n曲线进行评估得到在役风电机组的剩余疲劳寿命。
21.再一方面，本发明还提出了一种计算机设备，包括：一个或多个处理器；处理器，用于执行一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行时，实现上述的一种风电机组的疲劳寿命评估方法。
22.再一方面，本发明还提出了一种计算机可读存储介质，其上存有计算机程序，所述计算机程序被执行时，实现上述的一种风电机组的疲劳寿命评估方法。
23.与现有技术相比，本发明的有益效果为：一种风电机组的疲劳寿命评估方法、系统、设备及介质，包括：基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值；基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱；基于所述计算周期的疲劳载荷谱结合s-n曲线进行评估得到在役风电机组的剩余疲劳寿命；本发明考虑了风电机组实际运行外部环境条件和实际运行状态，准确的还原了在役风电机组实际疲劳载荷单位时长的时间历程；本法利用统计外推法准确描述了载荷分布的变化趋势并填补了数据库的空缺，减小疲劳分析过程产生的误差和难度；本发明无需高精度的传感器，只需要scada运行数据，使得更经济可靠。
附图说明
24.图1为现有技术的一种风电叶片剩余疲劳寿命评估方法流程图；图2为现有技术的一种风电机组疲劳寿命的快速预测的方法流程图；图3为本发明的一种风电机组的疲劳寿命评估方法流程图。
具体实施方式
25.为了解决由于测试或仿真时间短及外部环境条件限制所带来的问题，统计分析的方法被广泛应用于机组疲劳载荷分析研究中，我们称之为“统计外推法”。统计外推法是根据以往的经验以及对问题本身的认知程度，假设变量符合某一特定的函数，来估算总体的概率分布，由于分布函数的参数未知，因此要通过样本值来估计，样本值的选取具有随机性与普遍性。统计外推法不仅能够描述载荷分布的变化趋势，填补数据库的空缺，还能够基于样本对数据进行平滑的外推。本发明基于此法，对载荷雨流计数的幅值、均值和外部风速分布进行外推，获取在役风电机组的1h疲劳载荷谱和年疲劳载荷谱以评估在役机组的剩余疲劳寿命。为了更好地理解本发明，下面结合说明书附图和实施例对本发明的内容做进一步的说明。
26.实施例 1：一种风电机组的疲劳寿命评估方法，具体过程如图3所示，包括：步骤1，基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值；步骤2，基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱；步骤3，基于所述计算周期的疲劳载荷谱结合s-n曲线进行评估得到在役风电机组的剩余疲劳寿命。
27.本发明中计算周期、第一单位时长、第二单位时长和第三单位时长均为设定时长，且第一单位时长、第二单位时长和第三单位时长在计算周期范围内；本实施例以计算周期为一年，第一单位时长为1h，第二单位时长为10min，第三单位时长为24h为例进行介绍。
28.步骤1中的，基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值，具体包括：基于获取的风电场历史气象数据，统计出在役风电机组计算周期内第二单位时长的平均风速频率分布，并拟合成双参数威布尔分布；统计出在役机组计算周期内第二单位时长的湍流强度，并将湍流强度平均分成设定个数的等级，以备后续载荷仿真计算使用，本实施例中设定为6个等级；统计出在役机组计算周期内平均温度、湿度、气压、风剪切、风转向与入流角，以备后续载荷仿真计算使用。
29.将在役风电机组运行状态分为正常发电状态、停机状态和启停机制动状态，每种状态时间历程设置为第二单位时长。基于在役风电机组scada数据，统计出三种状态年频率，分别为正常发电状态年频率，停机状态年频率与启停机制动状态年频率，且。
30.基于在役风电机组的运行状态进行区间划分，得到各运行状态下的各风速区间：当运行状态为正常发电状态时，根据风电机组切入风速、切出风速将风速以设定间隔进行区间划分，本实施例中设定间隔为0.5m/s，即划分为：{[+0.5)、[+0.5+1.0)、
…
[-0.5,)}；当运行状态为停机状态或启停机制动状态时，基于初始风速以设定间隔等间隔递增至目标风速进行区间划分，本实施例中初始风速设定为0m/s，目标风速设定为50m/s，间隔设定为0.5m/s，即划分为：{[0 0.5)、[0.5 1.0)、
…
[49.5,50)}。
[0031]
基于各运行状态下的各风速区间确定载荷样本：载荷样本的采样时间为10min，采样频率为50hz。由于风速波动大，应用仿真软件计算时，每个风速区间取6个样本，每个样本的平均风速为区间风速的中值，湍流强度分别为6个湍流强度等级，其余气象参数取自风电场历史气象数据统计的年平均值，基于各运行状态下各风速区间内的载荷样本进行仿真计算，得到各运行状态下各风速区间内第一单位时长的载荷时序值。
[0032]
本发明基于统计外推法考虑了影响风电机组疲劳载荷的风电场真实外部风参数与机组实际运行状态，准确地还原了在役机组实际疲劳载荷第一单位时长的时间历程，对特定场址的在役机组的疲劳寿命尤其精确。
[0033]
步骤2中的，基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱，具体包括：对各运行状态下各风速区间的载荷时序值进行雨流计算，得到大小不同的疲劳载荷幅值x、载荷均值y以及相应的循环作用次数n，从而得到雨流矩阵n（x,y）。
[0034]
根据三维雨流矩阵的幅值统计载荷做双参数威布尔分布,其概率密度函数的计算式如下所示：
[0035]
式中，x为载荷幅值，f（x）为载荷幅值的概率密度函数，为形状参数，为尺度参数。
[0036]
均值统计载荷做正态分布的统计，其概率密度函数的计算式如下所示：
[0037]
式中，为载荷均值，为载荷均值的概率密度函数，为尺度参数，为位置参数。
[0038]
对载荷幅值与载荷均值进行相关性分析，检验验证均幅值相互独立，其联合概率密度函数的计算式如下所示：
[0039]
式中，x为载荷幅值，y为载荷均值，为载荷幅值和载荷均值的联合概率密度函数，为形状参数，为尺度参数，为尺度参数，为位置参数。
[0040]
载荷幅值的概率密度函数进行积分可得幅值累积频率，如下所示：
[0041]
式中，p为载荷幅值的累计频率，x为载荷幅值，为形状参数，为尺度参数。
[0042]
根据幅值累计频率可得到载荷幅值极大值如下所示：
[0043]
式中，为载荷幅值极大值，为尺度参数，p为载荷幅值的累计频率,取99.5%。
[0044]
由于风电机组载荷分布尾部较大的载荷对部件的疲劳损伤起主导作用，为使拟合的概率分布的尾部更接近实际载荷数据，根据载荷幅值和相应的循环作用次数得到s-n曲线，并将载荷幅值划分为各个等级的载荷。根据s-n曲线和线性损伤理论可得到各个等级的载荷产生的疲劳损伤和各个等级的载荷产生的疲劳损伤密度。
[0045]
各个等级的载荷产生的疲劳损伤的计算式如下所示：
[0046]
式中，为第k级载荷；为载荷谱第k级载荷产生的疲劳损伤；为第k级载荷对应的循环次数；n为循环次数；为第k级载荷m次方；为与材料性质相关的拟合参数。
[0047]
各个等级的载荷产生的疲劳损伤密度的计算式如下所示：
[0048]
式中，为第k级载荷；为第k级载荷产生的疲劳损伤密度；为载荷谱第k级载荷产生的疲劳损伤；为第k级载荷对应的循环次数；n为循环次数；为第k级载荷m次方；为级数。
[0049]
工程应用中定义为，可得到载荷幅值极小值。因此仅对大于载荷幅值极小值的载荷幅值的雨流计数进行拟合。
[0050]
根据载荷幅值极大值和载荷幅值极小值将载荷幅值划分100等分，载荷幅值范围为[
ꢀꢀ
]，取三种运行状态中的最小值，取三种运行状态中的最大值。根据幅值概率分布可得出各个风速区间内第一单位时长内各幅值区间的循环作用的次数，即可得出风电机组各运行状态下各风速区间内第一单位时长的疲劳载荷谱。
[0051]
根据风电场年平均风速概率分布，计算周期内各幅值区间的循环作用次数按下式计算：
[0052]
式中，为载荷幅值对应的循环次数；为i区间风速的分布频率；为风速区间i的风速；为发电状态下i风速区间中载荷幅值对应的循环次数；为停机状态下i风速区间中载荷幅值对应的循环次数；为启停机制动状态下i风速区间中载荷幅值对应的循环次数；为正常发电状态的年频率；a为计算周期；b为第三单位时长；为级数。
[0053]
由载荷幅值与载荷幅值对应的循环次数可得到风电机组计算周期疲劳载荷谱。本发明基于基于在役机组实际疲劳载荷第一单位时长的时间历程外推至计算周期疲劳载荷谱；而且本发明无需增加高精度传感器，经济又可靠，外推方法简便，计算量少，可在风电场中大规模推行。
[0054]
步骤3中的，基于所述计算周期的疲劳载荷谱结合s-n曲线进行评估得到在役风电机组的剩余疲劳寿命，包括：基于计算周期的疲劳载荷谱结合材料的s-n曲线进行对比，即可得到剩余疲劳寿命。
[0055]
具体包括：
材料的s-n曲线包括评估的疲劳载荷谱，证明在使用寿命期内；假设每年的气象、机组运行状态一致，年疲劳载荷谱中的循环数乘以机组运行年数即可得出机组当年的疲劳载荷谱。
[0056]
本发明利用s-n曲线和miner线性累计损伤理论从而估算出在役机组剩余疲劳寿命，进而保障风电机组安全运行，也为风电机组后评估提供可靠的新方法，对风电机组的延寿与技改具有指导意义；关键部件载荷包括叶片载荷、塔筒载荷与传动系统载荷。
[0057]
准确的疲劳寿命评估对于风电机组部件结构设计及运行维护安全性至关重要。通过本发明准确评估机组各部件疲劳寿命，可在部件发生灾难性故障或失效前进行维修或更换，确保风电机组长期安全稳定运行。
[0058]
实施例 2：一种风电机组的疲劳寿命评估系统，包括：仿真计算模块，用于基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值；统计外推模块，用于基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱；寿命评估模块，用于基于所述计算周期的疲劳载荷谱结合s-n曲线进行评估得到在役风电机组的剩余疲劳寿命。
[0059]
仿真计算模块，包括：风速区间子模块，用于基于所述在役风电机组的运行状态进行区间划分，得到各运行状态下的各风速区间；载荷样本子模块，用于基于所述各运行状态下的各风速区间对所述风电场外部气象参数进行采样，得到各运行状态下各风速区间内的载荷样本；载荷时序值子模块，用于基于所述各运行状态下各风速区间内的载荷样本进行仿真计算，得到各运行状态下各风速区间内第一单位时长的载荷时序值；其中，所述在役风电机组的运行状态包括：正常发电状态、停机状态和启停机制动状态。
[0060]
风速区间子模块，具体用于：当在役风电机组的运行状态为正常发电状态时，基于风电机组的切入风速以设定间隔等间隔递增至切出风速进行区间划分，得到正常发电状态下的各风速区间；当在役风电机组的运行状态为停机状态时，基于初始风速以设定间隔等间隔递增至目标风速进行区间划分，得到停机状态下的各风速区间；当在役风电机组的运行状态为启停机制动状态时，基于初始风速以设定间隔等间隔递增至目标风速进行区间划分，得到启停机制动状态下的各风速区间。
[0061]
统计外推模块，具体用于：基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值进行雨流计算，得到载荷幅值、载荷幅值极大值和载荷幅值极小值；基于在役风电机组计算周期内第二单位时长的平均风速频率分布结合各幅值区间，得到各风速区间内第一单位时长内各幅值区间的循环作用次数；
基于所述各风速区间内第一单位时长内各幅值区间的循环作用次数结合所述载荷幅值得到第一单位时长疲劳载荷谱；基于在役风电机组计算周期内第二单位时长的平均风速频率分布计算得到各风速区间内计算周期内各幅值区间的循环作用次数；基于所述载荷幅值结合所述各风速区间内计算周期内各幅值区间的循环作用次数得到计算周期疲劳载荷谱；其中，所述各幅值区间是基于载荷幅值极大值和载荷幅值极小值进行区间划分得到的。
[0062]
统计外推模块还包括载荷幅值极大值子模块，载荷幅值极大值子模块具体用于：基于在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值进行雨流计算，得到各运行状态下各风速区间内的载荷幅值、载荷均值以及相应的循环作用次数和雨流矩阵；基于所述雨流矩阵中的载荷幅值做双参数威尔布分布，得到载荷幅值的概率密度函数；基于所述载荷幅值的概率密度函数进行积分得到载荷幅值累积频率；基于所述载荷幅值累积频率得到载荷幅值的极大值。
[0063]
统计外推模块还包括载荷幅值极小值子模块，载荷幅值极小值子模块具体用于：基于所述载荷幅值和所述相应的循环作用次数得到s-n曲线，并基于所述载荷幅值划分为各个等级的载荷；基于所述s-n曲线结合线性损伤理论，计算得到各个等级的载荷产生的疲劳损伤和各个等级的载荷产生的疲劳损伤密度；基于所述各个等级的载荷产生的疲劳损伤和各个等级的载荷产生的疲劳损伤密度确定载荷幅值极小值。
[0064]
载荷幅值极小值子模块中的各个等级的载荷产生的疲劳损伤的计算式如下所示：
[0065]
式中，为载荷谱第k级载荷产生的疲劳损伤；为第k级载荷对应的循环次数；为第k级载荷；为第k级载荷m次方；为与材料性质相关的拟合参数。
[0066]
载荷幅值极小值子模块中的各个等级的载荷产生的疲劳损伤密度的计算式如下所示：
[0067]
式中，为载荷谱第k级载荷产生的疲劳损伤；为第k级载荷对应的循环次数；为第k级载荷产生的疲劳损伤密度；为第k级载荷；为第k级载荷m次方，为级数。
[0068]
统计外推模块中的计算周期内各幅值区间的循环作用次数按下式计算：
[0069]
式中，为载荷幅值对应的循环次数；为i区间风速的分布频率；为发电状态下i风速区间中载荷幅值对应的循环次数；为停机状态下i风速区间中载荷幅值对应的循环次数；为启停机制动状态下i风速区间中载荷幅值对应的循环次数；为正常发电状态的年频率；为运行状态的年频率；a为计算周期，设定为365；b为第三单位时长，设定为24；为风速区间i的风速；为区间序号。
[0070]
统计外推模块中的运行状态的年频率是基于数据采集与监视控制系统的数据统计得到的。
[0071]
风电场外部气象参数包括在役风电机组计算周期内的第二单位时长平均风速频率分布、第二单位时长湍流强度、平均温度、平均湿度、平均气压、平均风剪切、平均风转向和平均入流角。
[0072]
实施例 3：基于同一发明构思，本发明再一个实施例中，提供了一种计算机设备，该计算机设备包括处理器以及存储器，所述存储器用于存储计算机程序，所述计算机 程序包括程序指令，所述处理器用于执行所述计算机存储介质存储的程序指令。 处理器可能是中央处理单元(central processing unit，cpu)，还可以是其他通 用处理器、数字信号处理器(digital signal processor、dsp)、专用集成电路 (applicationspecific integrated circuit，asic)、现成可编程门阵列 (field-programmable gatearray，fpga)或者其他可编程逻辑器件、分立门或者 晶体管逻辑器件、分立硬件组件等，其是终端的计算核心以及控制核心，其适于 实现一条或一条以上指令，具体适于加载并执行计算机存储介质内一条或一条以 上指令从而实现相应方法流程或相应功能；本发明实施例所述的处理器可以用于 执行一种风电机组的疲劳寿命评估方法的步骤。
[0073]
实施例 4：基于同一发明构思，本发明再一个实施例中，本发明还提供了一种存储介质， 具体为计算机可读存储介质（memory），所述计算机可读存储介质是计算机设备 中的记忆设备，用于存放程序和数据。可以理解的是，此处的计算机可读存储介质既可以包括计算机设备中的内置存储介质，当然也可以包括计算机设备所支持 的扩展存储介质。计算机可读存储介质提供存储空间，该存储空间存储了终端的 操作系统。并且，在该存储空间中还存放了适于被处理器加载并执行的一条或一 条以上的指令，这些指令可以是一个或一个以上的计算机程序（包括程序代码）。 需要说明的是，此处的计算机可读存储介质可以是高速 ram 存储器，也可以是非不稳定的存储器(non-volatile memory)，例如至少一个磁盘存储器。可由处理 器加载并执行计算机可读存储介质中存放的一条或一条以上指令，以实现上述实施例中一种风电机组的疲劳寿命评估方法的相应步骤。
[0074]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、cd-rom、光学存储器等）上实施的计算机程序产
品的形式。
[0075]
本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。
[0076]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。
[0077]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的步骤。
[0078]
以上仅为本发明的实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均包含在发明待批的本发明的权利要求范围之内。

技术特征：
1.一种风电机组的疲劳寿命评估方法，其特征在于，包括：基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值；基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱；基于所述计算周期的疲劳载荷谱结合s-n曲线进行评估得到在役风电机组的剩余疲劳寿命。2.根据权利要求1所述方法，其特征在于，所述基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值，包括：基于所述在役风电机组的运行状态进行区间划分，得到各运行状态下的各风速区间；基于所述各运行状态下的各风速区间对所述风电场外部气象参数进行采样，得到各运行状态下各风速区间内的载荷样本；基于所述各运行状态下各风速区间内的载荷样本进行仿真计算，得到各运行状态下各风速区间内第一单位时长的载荷时序值；其中，所述在役风电机组的运行状态包括：正常发电状态、停机状态和启停机制动状态。3.根据权利要求2所述方法，其特征在于，所述基于所述在役风电机组的运行状态进行区间划分，得到各风速区间，包括：当在役风电机组的运行状态为正常发电状态时，基于风电机组的切入风速以设定间隔等间隔递增至切出风速进行区间划分，得到正常发电状态下的各风速区间；当在役风电机组的运行状态为停机状态时，基于初始风速以设定间隔等间隔递增至目标风速进行区间划分，得到停机状态下的各风速区间；当在役风电机组的运行状态为启停机制动状态时，基于初始风速以设定间隔等间隔递增至目标风速进行区间划分，得到启停机制动状态下的各风速区间。4.根据权利要求1所述方法，其特征在于，所述基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱，包括：基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值进行雨流计算，得到载荷幅值、载荷幅值极大值和载荷幅值极小值；基于在役风电机组计算周期内第二单位时长的平均风速频率分布结合各幅值区间，得到各风速区间内第一单位时长内各幅值区间的循环作用次数；基于所述各风速区间内第一单位时长内各幅值区间的循环作用次数结合所述载荷幅值得到第一单位时长疲劳载荷谱；基于在役风电机组计算周期内第二单位时长的平均风速频率分布计算得到各风速区间内计算周期内各幅值区间的循环作用次数；基于所述载荷幅值结合所述各风速区间内计算周期内各幅值区间的循环作用次数得到计算周期疲劳载荷谱；其中，所述各幅值区间是基于载荷幅值极大值和载荷幅值极小值进行区间划分得到
的。5.根据权利要求4所述方法，其特征在于，所述载荷幅值极大值按下述方法得到，包括：基于在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值进行雨流计算，得到各运行状态下各风速区间内的载荷幅值、载荷均值以及相应的循环作用次数和雨流矩阵；基于所述雨流矩阵中的载荷幅值做双参数威尔布分布，得到载荷幅值的概率密度函数；基于所述载荷幅值的概率密度函数进行积分得到载荷幅值累积频率；基于所述载荷幅值累积频率得到载荷幅值的极大值。6.根据权利要求5所述方法，其特征在于，所述载荷幅值极小值按下述方法得到，包括：基于所述载荷幅值和所述相应的循环作用次数得到s-n曲线，并基于所述载荷幅值划分为各个等级的载荷；基于所述s-n曲线结合线性损伤理论，计算得到各个等级的载荷产生的疲劳损伤和各个等级的载荷产生的疲劳损伤密度；基于所述各个等级的载荷产生的疲劳损伤和各个等级的载荷产生的疲劳损伤密度确定载荷幅值极小值。7.根据权利要求6所述方法，其特征在于，所述各个等级的载荷产生的疲劳损伤的计算式如下所示：式中，为载荷谱第k级载荷产生的疲劳损伤；为第k级载荷对应的循环次数；为第k级载荷；为与材料性质相关的拟合参数；为第k级载荷m次方。8.根据权利要求6所述方法，其特征在于，所述各个等级的载荷产生的疲劳损伤密度的计算式如下所示：式中，为载荷谱第k级载荷产生的疲劳损伤；为第k级载荷对应的循环次数；为第k级载荷产生的疲劳损伤密度；为第k级载荷；为第k级载荷m次方，为级数。9.根据权利要求4所述方法，其特征在于，所述计算周期内各幅值区间的循环作用次数按下式计算：式中，为载荷幅值对应的循环次数；为i区间风速的分布频率；为发电状态下i风速区间中载荷幅值对应的循环次数；为停机状态下i风速区间中载
荷幅值对应的循环次数；为启停机制动状态下i风速区间中载荷幅值对应的循环次数；为正常发电状态的年频率；a为计算周期；b为第三单位时长；为风速区间i的风速；为区间序号。10.根据权利要求9所述方法，其特征在于，所述运行状态的年频率是基于数据采集与监视控制系统的数据统计得到的。11.根据权利要求1所述方法，其特征在于，所述风电场外部气象参数包括在役风电机组计算周期内的第二单位时长平均风速频率分布、第二单位时长湍流强度、平均温度、平均湿度、平均气压、平均风剪切、平均风转向和平均入流角。12.一种风电机组的疲劳寿命评估系统，其特征在于，包括：仿真计算模块，用于基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值；统计外推模块，用于基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱；寿命评估模块，用于基于所述计算周期的疲劳载荷谱结合s-n曲线进行评估得到在役风电机组的剩余疲劳寿命。13.一种计算机设备，其特征在于，包括：一个或多个处理器；处理器，用于执行一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行时，实现如权利要求求1至 11中任一所述的一种风电机组的疲劳寿命评估方法。14.一种计算机可读存储介质，其特征在于，其上存有计算机程序，所述计算机程序被执行时，实现如权利要求 1 至 11 中任一所述的一种风电机组的疲劳寿命评估方法。

技术总结
一种风电机组的疲劳寿命评估方法、系统、设备及介质，包括：基于获取的风电场外部气象参数和在役风电机组实际运行状态进行仿真计算，得到在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值；基于所述在役风电机组各运行状态下各风速区间内第一单位时长的载荷时序值利用统计外推法计算得到计算周期的疲劳载荷谱；基于所述计算周期的疲劳载荷谱结合S-N曲线进行评估得到在役风电机组的剩余疲劳寿命；本发明考虑了风电机组实际运行外部环境条件和实际运行状态，准确的还原了在役风电机组实际疲劳载荷单位时长的时间历程；本法利用统计外推法准确描述了载荷分布的变化趋势并填补了数据库的空缺，减小疲劳分析过程产生的误差和难度。程产生的误差和难度。程产生的误差和难度。


技术研发人员：龚利策 秦世耀 王安庆 付德义 王瑞明 薛扬 孔令行 赵娜 贾海坤 李松迪
受保护的技术使用者：中国电力科学研究院有限公司
技术研发日：2023.07.27
技术公布日：2023/8/28