风力发电机组的塔架载荷建模方法和装置与流程

1.本技术涉及风力发电机组领域，尤其涉及一种风力发电机组的塔架载荷建模方法和装置。


背景技术：

2.风力发电机组的设计正在向着深远海、尺寸和容量大型化的方向发展，相较于陆上风力发电机组，海上风力发电机组塔筒支撑系统具有更强的柔性、更加严酷的风浪流耦合服役环境，其服役性能与安全性对经济高效的海上风能资源开发至关重要，因而对其精细化分析与整体可靠性设计提出了更高要求。风力发电机组包括桨叶、机舱、塔体、基础等部件，且存在转矩控制、桨叶控制与气动耦合等电
‑
机结构非线性耦合效应，建立合理的海上风力发电机组统分析模型是进行整体可靠性分析与设计的基础。
3.与陆上风力发电机组不同的是，除气动载荷外海上风力发电机组还同时受到波浪和洋流载荷的影响(如图1所示)。在海上风电场的建设中，塔架
‑
基础结构的成本占总造价的比例较高，根据海上风电场不同的场区的水文特征、地质条件、使用要求，选用不同的基础结构形式，是保证海上风力发电机组基础稳定性、可靠性和经济性的关键。多样的结构形式、复杂的海洋环境条件等因素使得基础结构研发和设计成为海上风电领域的重点和难点，了解各类环境、载荷参数对不同基础结构的影响规律，便可根据环境、载荷条件对风力发电机组机位、基础形式作出合理的选择。
4.目前，针对风电机组“桨叶
‑
塔架
‑
基础”建立的一体化模型主要采用模态分析法、多体动力学法和有限元分析法。对于模态分析法，叶轮面经历的气动载荷仅简单的通过额定风速和修正系数来确定，随后依据转子和塔体的变形相容性条件实现运动方程的耦合，使用该方法建立的一体化载荷模型方法自由度较少、计算方便，但存在精度较低，代数推导繁杂，且难以考虑具有较大形变的几何非线性等问题。多体动力学法可同时考虑刚体位移和弹性变形，能够满足一定的模拟精度，然而多体动力学法对细部力学分析较低，而这对诸如疲劳载荷评估和可靠性设计至关重要。对于有限元法，风电机组基于有限元的建模分为获得精确局部响应的细部建模和获得整体响应的整体建模，获得各个部件的刚度、质量阻尼矩阵，并通过叶片
‑
塔筒的连接矩阵在轮毂处实现了叶片与塔筒的耦合，然而该一体化建模方法难以同时分析刚体运动和弹性变形问题及叶片的旋转效应。此外，上述一体化建模方法忽略了气动载荷、波浪载荷、洋流特征以及风力发电机组自身的控制器设计的相互影响等问题，这些因素都影响了一体化建模的精度，导致与风力发电机组支撑结构实际承受载荷情况有差异。


技术实现要素：

5.本技术提供一种风力发电机组的塔架载荷建模方法和装置。
6.具体地，本技术是通过如下技术方案实现的：
7.本技术实施例的第一方面，提供一种风力发电机组的塔架载荷建模方法，包括：
8.针对风力发电机组的多个运行工况，通过蒙特卡洛法仿真获得多个载荷数据集，多个所述运行工况与多个所述载荷数据集一一对应，各载荷数据集均包括塔架分别在对应运行工况的多个子工况下的载荷值，多个所述运行工况的至少一个工况参数的特征值的参考值大小不相等，每个运行工况的多个子工况对应包括该运行工况在不同时间段内的工况参数的特征值，每个运行工况的多个子工况的至少一个工况参数的特征值的变化范围不同，所述工况参数包括风浪流参数及所述风力发电机组的控制参数；
9.分析各载荷数据集中的载荷值的概率分布，所述概率分布用于表征各载荷值的概率大小；
10.根据各运行工况对应的载荷数据集中的载荷值的概率分布，确定能够代表该载荷数据集的有效载荷值，并确定与所述有效载荷值相对应的工况参数的特征值；
11.根据所述载荷数据集的有效载荷值及与所述有效载荷值相对应的工况参数的特征值，采用敏感性分析方法确定每个工况参数对所述塔架载荷的影响程度；
12.根据每个工况参数对所述塔架载荷的影响程度，确定多个所述工况参数中的有效影响工况参数；
13.根据多个所述载荷数据集的有效载荷值及对应的有效影响工况参数的特征值，通过多项式展开法对所述塔架载荷进行建模。
14.可选地，所述运行工况包括正常工况和极限工况。
15.可选地，所述针对风力发电机组的多个运行工况，通过蒙特卡洛法仿真获得多个载荷数据集，包括：
16.根据所述风力发电机组的各运行工况的工况参数的参考值，生成该运行工况的多个子工况序列，同一运行工况对应的多个子工况在时间序列上的各工况参数的值是随机变化的；
17.对各运行工况的多个所述子工况序列通过蒙特卡洛法进行仿真处理，获得塔架在风力发电机组的各子工况下的载荷力矩序列；
18.根据所述塔架在风力发电机组的各子工况下的载荷力矩序列，使用雨流计数法确定所述载荷力矩序列的等效损坏载荷值，获得各运行工况对应的载荷数据集。
19.可选地，所述针对风力发电机组的多个运行工况，通过蒙特卡洛法仿真获得多个载荷数据集，包括：
20.根据统计法对各载荷数据集中的载荷值的对称程度和集中程度进行分析，获得各载荷数据集中的各载荷值的数据集分布特征值；
21.根据各载荷数据集中的各载荷值的数据集分布特征值的对称程度和集中程度，确定该载荷数据集中的载荷值的分布特征；
22.判断各载荷数据集中的载荷值的分布特征是否满足预设分布特征，所述预设分布特征具有对称性，且集中程度大于或等于预设程度值；
23.当所述载荷数据集中的载荷值的分布特征满足所述预设分布特征时，停止对应运行工况下的蒙特卡洛法仿真；
24.当所述载荷数据集中的载荷值的分布特征不满足所述预设分布特征时，继续对应运行工况下的蒙特卡洛法仿真。
25.可选地，所述统计法为jarque
‑
bera统计法。
26.可选地，所述分析各载荷数据集中的载荷值的概率分布，包括：
27.根据非参数核密度估计方法，分别估计各载荷数据集中的载荷值的概率分布。
28.可选地，所述有效载荷值为所述载荷数据集中概率值最大的载荷值。
29.可选地，所述根据所述载荷数据集的有效载荷值及与所述有效载荷值相对应的工况参数的特征值，采用敏感性分析方法确定每个工况参数对所述塔架载荷的影响程度，包括：
30.根据所述载荷数据集的有效载荷值及所述工况参数的特征值之间的对应关系，采用敏感性分析方法确定每个工况参数对塔架载荷的影响程度。
31.可选地，所述根据多个所述载荷数据集的有效载荷值及对应的工况参数大小以及敏感性分析方法，确定每个工况参数对所述塔架载荷的影响程度，包括：
32.将多个所述载荷数据集均分成两组，获得两个载荷数据集组合；
33.根据两个所述载荷数据集组合，分别获得两个第一参数矩阵，其中，各第一参数矩阵的同一行为同一载荷数据集的有效载荷值对应的工况参数的特征值，同一列为不同载荷数据集的有效载荷值的同一类型的工况参数的特征值；
34.将两个所述第一参数矩阵的各相同列依次对换，分别获得该相同列对应的两个第二参数矩阵；
35.确定两个所述第一参数矩阵分别与对应的有效载荷值之间的函数关系；
36.根据所述函数关系，确定各相同列对应的两个所述第二参数矩阵分别对应的有效载荷值；
37.根据各相同列对应的其中一个所述第二参数矩阵对应的有效载荷值以及两个所述第一参数矩阵分别对应的有效载荷值，确定每个工况参数对所述塔架载荷的影响程度。
38.可选地，所述根据每个工况参数对所述塔架载荷的影响程度，确定多个所述工况参数中的有效影响工况参数，包括：
39.当所述影响程度大于预设程度阈值时，确定对应的工况参数为有效影响工况参数；或者，
40.确定所述影响程度靠前的预设数量的工况参数作为多个所述工况参数中的有效影响工况参数。
41.可选地，所述根据多个所述载荷数据集的有效载荷值及对应的有效影响工况参数的特征值，通过多项式展开法对所述塔架载荷进行建模，包括：
42.建立所述塔架载荷的多项式模型，所述多项式模型以所述有效影响工况参数的特征值为自变量，以所述有效载荷值为因变量；
43.将所述多个所述载荷数据集的有效载荷值及对应的有效影响工况参数的特征值分别代入所述多项式模型，获得所述多项式模型中每一项的系数。
44.可选地，所述风浪流参数包括风速、湍流强度、海浪高度、海浪周期、洋流速度中的至少一种，所述控制参数包括有功功率设定值、桨距角参考值、和转速设定值中的至少一种。
45.本技术实施例的第二方面，提供一种风力发电机组的塔架载荷建模装置，包括一个或多个处理器，用于实现第一方面中任一项所述的风力发电机组的塔架载荷建模方法。
46.本技术实施例的第三方面，提供一种计算机可读存储介质，其上存储有程序，该程
序被处理器执行时，实现第一方面中任一项所述的风力发电机组的塔架载荷建模方法。
47.根据本技术实施例提供的技术方案，采用蒙特卡洛法对多种运行工况进行仿真实验获得对应的塔架载荷数据集，克服了风力和水利随机性造成塔架载荷不确定性的影响，并分析各载荷数据集中的载荷值的概率分布以确定代表各载荷数据集的有效载荷值，再根据多个载荷数据集的有效载荷值及对应的工况参数的特征值，采用敏感性分析方法确定每个工况参数对塔架载荷的影响程度，并剔除对塔架载荷关联度较低的因素以简化建模的复杂度，最后根据多项式展开法实现对塔架载荷的精确建模，对进行整机可靠性分析与设计、优化单机控制策略和风电场运行维护策略，具有重要的工程应用价值。
48.应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本技术。
附图说明
49.此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本技术的实施例，并与说明书一起用于解释本技术的原理。
50.图1是目前海上风力发电机组塔架经历风
‑
浪
‑
流耦合的载荷工况示意图；
51.图2是本技术一示例性实施例示出的一种风力发电机组的塔架载荷建模方法的流程示意图；
52.图3是本技术一示例性实施例示出的一种同一运行工况下的不同子工况的风速值示意图；
53.图4是本技术一示例性实施例示出的一种塔架在风力发电机组的某一子工况下的载荷力矩序列图；
54.图5是本技术一示例性实施例示出的一种根据载荷数据集的有效载荷值及对应的有效影响工况参数的特征值，通过多项式展开法对塔架载荷进行建模的实现过程示意图；
55.图6是本技术一示例性实施例示出的一种风力发电机组的塔架载荷建模装置的结构示意图。
具体实施方式
56.这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本技术相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本技术的一些方面相一致的装置和方法的例子。
57.在本技术使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本技术。在本技术和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
58.应当理解，尽管在本技术可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本技术范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在
……
时”或“当
……
时”或“响应于确定”。
59.随着海上风力发电机组scada数据库的日益完善，有效利用风
‑
浪
‑
流和机组运行状态的检测数据能够提高估计风力发电机组塔架一体化载荷的精确度，分析风浪流数据和控制数据对塔架载荷的影响，建立合理的海上风力发电机组塔架载荷一体化精细模型不仅是进行整机可靠性分析与设计的基础，也是机组降载控制策略设计的重要参考，具有重要的工程应用价值。
60.海上风电机组大部件的设计寿命是25年，尽管在设计中的载荷参数能够满足25年的疲劳分析和静力分析，而且整个风力发电机组运行平稳，但却在使用几年后发生失效。问题不是出在静力和疲劳分析上，而是由于设计流程发生脱节导致的。脱节主要分为两部分，一方面是由于载荷计算与结构设计流程上发生了脱节，另一方面则发生在控制流程设计和结构设计发生脱节。特别是大型海上风力发电机组缺乏一种“风
‑
浪
‑
流
‑
控制”全耦合一体化高精度的解决方案。
61.对于此，本技术采用蒙特卡洛法对多种运行工况进行仿真实验获得对应的塔架载荷数据集，克服了风力和水利随机性造成塔架载荷不确定性的影响，并分析各载荷数据集中的载荷值的概率分布以确定代表各载荷数据集的有效载荷值，再根据多个载荷数据集的有效载荷值及对应的工况参数的特征值，采用敏感性分析方法确定每个工况参数对塔架载荷的影响程度，并剔除对塔架载荷关联度较低的因素以简化建模的复杂度，最后根据多项式展开法实现对塔架载荷的精确建模，对进行整机可靠性分析与设计、优化单机控制策略和风电场运行维护策略，具有重要的工程应用价值。
62.下面结合附图，对本技术的风力发电机组的塔架载荷建模方法和装置进行详细说明。在不冲突的情况下，下述的实施例及实施方式中的特征可以相互组合。
63.图2是本技术一示例性实施例示出的一种风力发电机组的塔架载荷建模方法的流程示意图；本技术实施例的风力发电机组的塔架载荷建模方法的执行主体可以包括风力发电机组，也可以为其他具有数据处理能力的设备，如计算机。
64.参见图2，本技术实施例提供的一种风力发电机组的塔架载荷建模方法可包括步骤s21～s25。
65.其中，在s21中、针对风力发电机组的多个运行工况，通过蒙特卡洛法仿真获得多个载荷数据集，多个运行工况与多个载荷数据集一一对应，各载荷数据集均包括塔架分别在对应运行工况的多个子工况下的载荷值，多个运行工况的至少一个工况参数的特征值的参考值大小不相等，每个运行工况的多个子工况对应包括该运行工况在不同时间段内的工况参数的特征值，每个运行工况的多个子工况的特征值的至少一个工况参数的变化范围不同，工况参数包括风浪流参数及风力发电机组的控制参数。
66.本技术实施例中，运行工况可包括正常工况和极限工况，同时考虑正常工况和极限工况，能够有效克服风力和水利变化的随机性造成塔架载荷的不确定性。在其他实施例中，运行工况可仅包括正常工况。一个运行工况对应一个载荷数据集，每个载荷数据集均包括塔架分别在风力发电机组的多个子工况下的载荷值，即每个运行工况对应多个子工况。示例性的，多个运行工况包括工况a、工况b和工况c，工况a对应的子工况包括子工况a1、子工况a2和子工况a3、
…
、子工况an，工况b对应的子工况包括子工况b1、子工况b2、子工况b3、
…
、子工况bn，工况c对应的子工况包括子工况c1、子工况c2、子工况c3、
…
、子工况an。
67.示例性地，针对工况a，通过蒙特卡洛法仿真获得载荷数据集a11，其中，载荷数据集a11包括塔架在分别子工况a1、a2、a3、
…
、an下的载荷值，一个运行工况对应一个载荷数据集，一个子工况对应一个载荷值。
68.本技术实施例的同一运行工况及子工况包括相同类型的工况参数，不同运行工况的工况参数的类型完全相同。同一运行工况的不同子工况的同一种工况参数的特征值的变化范围的参考值大小相等，示例性地，对于同一运行工况，不同子工况的工况参数的特征值是在该子工况所属的运行工况的对应工况参数的特征值的参考值的上下波动获得的。
69.示例性地，同一运行工况包括相同类型的工况参数，并且，同一种类型的工况参数的特征值的参考值相等，如，风速的参考值为15m/s，湍流强度的参考值为10％，海浪高度的参考值为1m，海浪周期的变化范围的参考值为10s。
70.同一运行工况的多个子工况分别是该运行工况中对应时间段内的工况参数的特征值大小，在每一个时间段内的子工况表现出来的效果是不一样的，即每个运行工况的多个子工况的至少一个工况参数的特征值的变化范围不同，如图3所示，子工况1、子工况2、子工况3、子工况4、
…
、子工况n对应的风速值分别为同一运行工况的不同时间段内的风速值。另外，子工况1、子工况2、子工况3及子工况4、
…
、子工况n对应的其他工况参数的特征值也为对应时间段内的。
71.需要说明的是，若运行工况为正常工况，则该运行工况对应的多个子工况也分别为正常工况；若运行工况为极限工况，则该运行工况对应的多个子工况也分别为极限工况。风浪流参数可包括平均风速、湍流强度、海浪高度、海浪周期、洋流速度中的至少一种，但不限于此；控制参数可包括有功功率设定值、桨距角参考值和转速设定值中的至少一种，但不限于此。
72.示例性的，在一些实施例中，风浪流参数包括风速、湍流强度、海浪高度、海浪周期和洋流速度，控制参数包括有功功率设定值和转速设定值，同时考虑到多种风浪流参数及多种控制参数，能够提高建模的精确度和准确性。
73.其中，运行工况及子工况包括相同类型的工况参数是指运行工况及子工况均包括相同类型的风浪流参数及相同类型的风力发电机组的控制参数，如，运行工况及子工况均包括风速、湍流强度、海浪高度、海浪周期、洋流速度、有功功率设定值、桨距角度参考值和转速设定值这些工况参数。
74.同一运行工况对应的多个子工况的具有相同的工况参数，但同一运行工况对应的子工况在时间序列上的各工况参数的值具有随机性，具体表现为在时间序列上工况的幅值、频率波动变化的特征具有任意性。
75.一种针对风力发电机组的多个运行工况，通过蒙特卡洛法仿真获得多个载荷数据集的实现过程可包括但不限于如下步骤：
76.(1)、根据风力发电机组的各运行工况的工况参数的参考值，生成相同该运行工况的多个子工况序列，同一运行工况对应的多个子工况在时间序列上的各工况参数的值是随机变化的。
77.即同一运行工况对应的多个子工况序列虽然具有相同的工况参数，但在时间序列上具有随机性。
78.例如，可以根据下述表1和表2分别给出的风力发电机组的正常工况和极限工况的
参考值以及表3给出的风力发电机组的控制参数设计，随机生成k个具有相同风浪流参数及风力发电机组的控制参数的子工况序列，可选地，各子工况序列的时长为10分钟，当然，各子工况序列的时长也可设置为其他时长大小。各子工况序列的时长大小可以相等，也可以部分不相等，具体可根据需要设置各子工况序列的时长大小。
79.表1
80.风浪流正常参数参数设定范围自由度风速[m/s][3:2:25]12湍流强度[％][1:3:15]5海浪高度[m][0:2:6]3海浪周期[s][10:4:24]6洋流速度[m/s][0:0.5:2]4发电机组转速[rpm][5.0:0.5:7.5]6
[0081]
表2
[0082][0083]
表3
[0084]
控制参数数值额定功率[mw]15叶片半径[m]120切入、额定、切出风速[m/s]3、10.6、25额定转速[rpm]7.56塔基、塔顶直径[m]10、6.5塔基、塔顶厚度[mm]100、30
[0085]
需要说明书的是，对于表1中风速的参数设定范围[3:2:25]，3表示最小风速，25表示最大风速，2表示多个运行工况的取值间隔，其他参数的参数设定范围的含义类似，不再赘述。比如，风速采用了3m/s、5m/s、7m/s、9m/s、11m/s、13m/s、15m/s、17m/s、19m/s、21m/s、23m/s、25m/s这12个风速值，自由度则为12。
[0086]
示例性的，运行工况1
‑
12的风速分别为3m/s、5m/s、7m/s、9m/s、11m/s、13m/s、15m/
s、17m/s、19m/s、21m/s、23m/s、25m/s，每个运行工况的多个子工况的风速则在对应风速值周围上下变化，不同子工况的风速的变化范围不同，例如，对于运行工况1，包括子工况11和12。
[0087]
(2)、对各运行工况的多个子工况序列通过蒙特卡洛法进行仿真处理，获得塔架在风力发电机组的各运行工况下的载荷力矩序列。
[0088]
对上述实施例中的各运行工况的k个子工况序列通过蒙特卡洛法进行仿真处理，获得塔架在风力发电机组的各运行工况下的k个载荷力矩序列。
[0089]
需要说明的是，若子工况序列为正常工况，则对应的载荷力矩序列为疲劳载荷力矩序列；若子工况序列为极限工况，则对应的载荷力矩序列为极限载荷力矩序列。
[0090]
(3)、根据塔架在风力发电机组的各运行工况下的载荷力矩序列，使用雨流计数法确定载荷力矩序列的等效损坏载荷值，获得各个运行工况对应的载荷数据集。
[0091]
雨流计数法的公式如下：
[0092][0093]
公式(1)中，del为等效损坏载荷值，n
i
和分别为通过雨流计数法计算得到的循环数和循环幅值，t为载荷时间历程评估周期，f为给定的正弦载荷频率，m为材料s
‑
n曲线中的斜率，m是由部件材料属性决定。
[0094]
在风力发电机组经历子工况后，塔架会受到振动力的作用，这样来回的振动会对塔架造成损伤，如图4所示，为塔架在风力发电机组的某一子工况下的载荷力矩序列图。通过雨流计数法来描述这样的损伤，这个值可以代表塔架的损坏程度，即等效损坏载荷。
[0095]
一种通过蒙特卡洛法仿真获得塔架多个载荷数据集的实现过程可包括但不限于如下步骤：
[0096]
(1)、根据统计法对各载荷数据集中的载荷值的对称程度和集中程度进行分析，获得各载荷数据集中的各载荷值的数据集分布特征值；
[0097]
可使用偏度表征对称程度，使用峰度表征集中程度。
[0098]
在一些实施例中，使用jarque
‑
bera(总体分布特性检验)统计法对各载荷数据集中的载荷值的偏度和峰度进行分析，其中，jarque
‑
bera统计法的计算公式为：
[0099][0100]
公式(2)中，jb为数据集分布特征值，c
sk
和c
sh
分别为载荷数据集的偏度系数和峰度系数，可选地，偏度系数c
sk
为0，峰度系数大于等于3，del
i
为载荷数据集的第i个载荷值，
d
i
e为载荷数据集的平均值，为载荷数据集的方差，k为各载荷数据集中的载荷值的数量。
[0101]
在一些实施例中，也可以使用其他统计法对各载荷数据集中的载荷值的偏度和峰度进行分析。
[0102]
(2)、根据各载荷数据集中的各载荷值的数据集分布特征值的对称程度和集中程度，确定该载荷数据集中的载荷值的分布特征；
[0103]
(3)、判断各载荷数据集中的载荷值的分布特征是否满足预设分布特征，预设分布特征具有对称性，且集中程度大于或等于预设程度值(即预设分布特征具有尖峰)；当载荷数据集中的载荷值的分布特征满足预设分布特征时，进入步骤(4)；当载荷数据集中的载荷值的分布特征不满足预设分布特征时，进入步骤(5)；
[0104]
(4)、停止对应运行工况下的蒙特卡洛法仿真；
[0105]
(5)、继续对应运行工况下的蒙特卡洛法仿真。
[0106]
例如，当载荷数据集中的载荷值的分布特征满足预设分布特征时，停止该载荷数据集对应的运行工况下的蒙特卡洛法仿真；当载荷数据集中的载荷值的分布特征不满足预设分布特征时，继续该载荷数据集对应的运行工况下的蒙特卡洛法仿真。
[0107]
直到所有载荷数据集中的载荷值的分布特征均满足预设分布特征的要求。
[0108]
在s22中、分析各载荷数据集中的载荷值的概率分布，概率分布用于表征各载荷值的概率大小。
[0109]
可根据非参数核密度估计方法，分别估计各载荷数据集中的载荷值的概率分布，非参数核密度估计方法不依赖于先验知识，准确度高。
[0110]
非参数核密度估计方法的核密度函数表达式为：
[0111][0112]
公式(3)中，del
i
(i＝1，2，...，k)为各载荷数据集中的载荷值，f
h
(x)是内核函数(高斯核函数是处理非线性问题最常用的函数之一，其核方程为h是参数带宽，x为各载荷数据集中的最小载荷值至最大载荷值之间的任意一个。
[0113]
对于各运行工况，可以获得相应的载荷数据集del
n
＝[del
n1
，del
n2
，...，del
nk
]，其中，del
nk
表示基于第n个运行工况的蒙特卡洛法仿真实验获得的第k个载荷值。对于塔架载荷的载荷数据集del
n
，通过(3)中所示的非参数核密度估计方法获得载荷数据集的概率分布。
[0114]
应当理解地，也可基于其他方法分析各载荷数据集中的载荷值的概率分布。
[0115]
在s23中、根据各载荷数据集中的载荷值的概率分布，确定能够代表该载荷数据集的有效载荷值，并确定与所述有效载荷值相对应的工况参数的特征值。
[0116]
其中，载荷数据集的有效值与工况参数的特征值的具有对应关系。
[0117]
在一些实施例中，各载荷数据集的有效载荷值的数量为一个，例如，各载荷数据集
的有效载荷值为各载荷数据集中概率值最大的载荷值，基于s22确定出的概率分布的分布特征可以得到各载荷数据集中最大概率的载荷值。
[0118]
随后，各载荷数据集中最大概率的载荷值作为与对应运行工况下的有效载荷值，即：
[0119][0120]
公式(4)中，v
n，
ti
n
，h
sea_n
，f
sea_n
，v
sea_n
，p
set_n
，β
ref_n
，ω
set_n
分别表示运行工况n的平均风速、湍流强度、海浪高度、海浪周期、洋流速度、有功功率设定值、桨距角参考值、转速设定值。
[0121]
在另外一些实施例中，各载荷数据集的有效载荷值的数量为多个，如，各载荷数据集的有效载荷值为各载荷数据集中概率大小大于预设概率阈值的载荷值，预设概率阈值的大小可以根据需要设置，如80％、85％或90％等等。
[0122]
在s24中、根据载荷数据集的有效载荷值及与有效载荷值相对应的与工况参数的特征值，采用敏感性分析方法确定每个工况参数对塔架载荷的影响程度。
[0123]
具体可根据载荷数据集的有效载荷值与工况参数的特征值之间的对应关系，采用敏感性分析方法确定每个工况参数对塔架载荷的影响程度。
[0124]
可根据载荷数据集的有效载荷值与工况参数的特征值之间的对应关系，并通过敏感性分析方法，确定每个工况参数对塔架载荷的影响程度，即使用敏感性分析方法确定每个工况参数对塔架载荷的影响程度，从而剔除对塔架载荷关联度较低的工况参数，简化塔架载荷建模的复杂度。
[0125]
下面，介绍敏感性分析方法确定每个工况参数对塔架载荷的影响程度的实现过程：
[0126]
(1)、将多个载荷数据集均分成两组，获得两个载荷数据集组合；
[0127]
(2)、根据两个载荷数据集组合，分别获得两个第一参数矩阵，其中，各第一参数矩阵的同一行为同一载荷数据集的有效载荷值对应的工况参数的特征值，同一列为不同载荷数据集的有效载荷值的同一类型的工况参数的特征值；
[0128]
(3)、将两个第一参数矩阵的各相同列依次对换，分别获得该相同列对应的两个第二参数矩阵；
[0129]
(4)、确定两个第一参数矩阵分别与对应的有效载荷值之间的函数关系；
[0130]
(5)、根据函数关系，确定各相同列对应的两个第二参数矩阵分别对应的有效载荷值；
[0131]
(6)、根据各相同列对应的其中一个第二参数矩阵对应的有效载荷值以及两个第一参数矩阵分别对应的有效载荷值，确定每个工况参数对塔架载荷的影响程度。
[0132]
示例性的，通过步骤s21～s23，获得一个(2k，m)的随机数矩阵，其中，2k为样本容量，即2k个载荷数据集，m为各载荷数据集的工况参数的个数，即风浪流参数和风力发电机组的控制参数的个数，定义两个第一参数矩阵矩阵(a和b)，每个第一参数矩阵包含一半的样本。
[0133][0134]
公式(5)中，x1、x2、
…
、x
m
分别为1
‑
m的工况参数，a和b分别包含1
‑
k的载荷数据集的工况参数、(k 1)
‑
(k k)的载荷数据集的工况参数。
[0135]
定义第二参数矩阵c
i
和d
i
，为了充分利用第一参数矩阵a和b，其中，c
i
由将a的第i列分配到b的第i列，d
i
通过将b的第i列分配到a的第i列。
[0136][0137]
将a、b、c
i
和d
i
分别作为工况参数的输入，得到不同组的(k，1)的塔架载荷向量(即上述函数关系)，即：
[0138][0139]
公式(7)中，y为有效载荷值。
[0140]
然后，计算各工况参数对塔架载荷的一阶敏感性系数s
i
：
[0141][0142]
公式(8)中，s
i
为第i个工况参数对塔架载荷的一阶敏感性系数，x
‑
i
代表除了x
i
所有输出变量的向量。
[0143]
每个工况参数对塔架载荷的一阶敏感性系数即为每个工况参数对塔架载荷的影响程度。
[0144]
可基于方差的敏感性分析方法确定每个工况参数对塔架载荷的影响程度，也可基于其他算法的敏感性分析方法确定每个工况参数对塔架载荷的影响程度。
[0145]
在s25中、根据每个工况参数对塔架载荷的影响程度，确定多个工况参数中的有效影响工况参数。
[0146]
在完成敏感性分析后，为了降低塔架载荷建模的复杂程度，去除对塔架载荷影响程度较低的因素。在一些实施例中，当影响程度大于预设程度阈值时，确定对应的工况参数为有效影响工况参数，从而去除对塔架载荷影响程度较低的因素，降低塔架载荷建模的复杂度。其中，预设程度阈值的大小可以根据需要设置，如预设程度阈值可以为20％、25％或者30％等等。
[0147]
在另外一些实施例中，确定影响程度靠前的预设数量的工况参数作为多个工况参数中的有效影响工况参数，从而去除对塔架载荷影响程度较低的因素，降低塔架载荷建模的复杂度。其中，预设数量可以根据需要设置，如，工况参数的数量为7个，预设数量可以为6个或5个等等。
[0148]
在s26中、根据多个载荷数据集的有效载荷值及对应的有效影响工况参数的特征值，通过多项式展开法对塔架载荷进行建模。
[0149]
参见图5，一种根据多个载荷数据集的有效载荷值及对应的有效影响工况参数的特征值，通过多项式展开法对塔架载荷进行建模的实现过程可包括如下步骤：
[0150]
s51、建立塔架载荷的多项式模型，多项式模型以有效影响工况参数的特征值为自变量，以有效载荷值为因变量；
[0151]
s52、将多个载荷数据集的有效载荷值及对应的有效影响工况参数的特征值分别代入多项式模型，获得多项式模型中每一项的系数。
[0152]
可使用任意多项式混沌展开的方法建立塔架载荷输出的多项式模型，取决于在事件ω、代数λ、概率度量γ下的随机独立的矢量输入x(包含风浪流参数和控制参数，假设矢量输入有n个工况参数)，作为随机多元正交多项式ψ
i
(x)和确定性系数c
i
的线性组合：
[0153][0154]
公式(9)中，系数c
i
通过使用正交特性获得，展开式中m个项的数量取决于输入的工况参数的个数n和展开式阶数p，即m＝(n p)！/(n！p！)。
[0155]
ψ
i
为一元正交或正交多项式ψ
i
(x)的乘积：
[0156][0157]
公式(11)中，指数(α)是一个m
×
n的矩阵，包含展开式项i中参数编号j的响应阶数。以二元二次方程为例，当ψ
i
(x)＝x2 xy y2时：x和y表示为输入的工况参数。
[0158]
参见图6，本技术实施例还提供一种风力发电机组的塔架载荷建模装置，包括一个或多个处理器，用于实现上述实施例中的风力发电机组的塔架载荷建模方法。
[0159]
风力发电机组的塔架载荷建模装置的实施例可以应用在风力发电机组上。装置实施例可以通过软件实现，也可以通过硬件或者软硬件结合的方式实现。以软件实现为例，作为一个逻辑意义上的装置，是通过其所在风力发电机组的处理器将非易失性存储器中对应的计算机程序指令读取到内存中运行形成的。从硬件层面而言，如图6所示，为本技术风力发电机组的塔架载荷建模装置所在风力发电机组的一种硬件结构图，除了图6所示的处理器、内部总线、内存、网络接口、以及非易失性存储器之外，实施例中装置所在的风力发电机组通常根据该风力发电机的实际功能，还可以包括其他硬件，对此不再赘述。
[0160]
本技术实施例还提供一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时，实现上述实施例中的风力发电机组的塔架载荷建模方法。
[0161]
所述计算机可读存储介质可以是前述任一实施例所述的风力发电机组的内部存储单元，例如硬盘或内存。所述计算机可读存储介质也可以是风力发电机的外部存储设备，例如所述设备上配备的插接式硬盘、智能存储卡(smart media card，smc)、sd卡、闪存卡(flash card)等。进一步的，所述计算机可读存储介质还可以既包括风力发电机组的内部存储单元也包括外部存储设备。所述计算机可读存储介质用于存储所述计算机程序以及所述风力发电机组所需的其他程序和数据，还可以用于暂时地存储已经输出或者将要输出的
数据。
[0162]
以上所述仅为本技术的较佳实施例而已，并不用以限制本技术，凡在本技术的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本技术保护的范围之内。