一种考虑台风影响的海上风电场多状态可靠性获取方法与流程

[0001]
本发明涉及电力系统发电输电领域，尤其是涉及一种考虑台风影响的海上风电场多状态可靠性获取方法。


背景技术：

[0002]
近年来能源枯竭、环境污染问题日益严峻，而新能源发电具有无污染、可持续等特点，其应用前景极为广阔。风能是清洁的可再生能源，大力开发可再生能源是合理调整电源结构的需要，也是国内能源发展战略的重要组成部分。而海上风电场是风力发电的研究热点和未来风能利用的重要方向。
[0003]
相较陆上风电而言，海上风电具有三大优势：一是海上风速高，风速较为稳定，海上风力资源丰富，比陆地风力发电产能大；二是海上风电场靠近沿海用电负荷中心，可减少远距离输电成本；三是对环境负面影响较少，允许安装单机容量更大的风机，非常适合大规模开发。目前世界上最大的海上风力发电厂为英国的“伦敦阵列项目”，总装机容量为630mw。我国海上风力资源丰富，开发海上风力资源具有重要的现实意义。在建的23个海上风电项目分布在江苏、福建、浙江、广东、河北、辽宁和上海7个省(市)海域。其中江苏海域7个在建项目共计2107.34mw；福建海域6个在建项目共计1620mw；广东海域5个在建项目共计1500mw；辽宁海域2个在建项目共计599.85mw；河北海域1个在建项目300mw；浙江海域1个在建项目252mw；上海海域1个在建项目100mw。
[0004]
但是海上风电场的运行条件严酷，台风等极端恶劣的天气会使电气设备故障率增大，海上维修可进入性变差，从而严重影响海上风电场的可靠性水平。可靠性对于海上风电场的运行具有重要的意义。随着海上风电机组单机容量和海上风电场规模的不断增大，在考虑台风影响的条件下，对海上风电场寿命期内的可靠性的评估已成为亟须解决的关键问题之一。


技术实现要素：

[0005]
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种考虑台风影响的海上风电场多状态可靠性获取方法。
[0006]
本发明的目的可以通过以下技术方案来实现：
[0007]
一种考虑台风影响的海上风电场多状态可靠性获取方法，用于得到海上风电机组寿命周期内每个设定的单位时间内的可靠性参数数据，包括以下步骤：
[0008]
s1：获取海上风电机组的寿命周期，并获取寿命周期内风电机组每一单位时间的初始故障率和单台风机的初始风速数据；
[0009]
s2：获取寿命周期内每一单位时间内的台风关键参数，并根据台风对单台机组故障率的影响，对风电机组寿命周期内每一单位时间的初始故障率进行修正；
[0010]
s3：利用单台风机的初始风速数据和聚类算法对单台风机进行出力状态划分；
[0011]
s4：利用马尔可夫链模型，根据单台风机的出力状态和修正后的初始故障率，建立
寿命周期各单位时间的单台机组多状态模型；
[0012]
s5：组合单台机组多状态模型，建立海上风电场的多状态可靠性模型，获取海上风电机组寿命周期内每个单位时间内的可靠性参数数据，用于对海上风电场的可靠性评估。
[0013]
进一步地，所述的台风关键参数包括年台风出现频率、台风移动方向、台风移动速度、台风最小距离和台风中心气压差。
[0014]
进一步地，所述的台风对单台机组故障率的影响通过据batts风场模型以及风速与故障率关系公式确定。
[0015]
更进一步地，所述的batts风场模型的表达式为：
[0016][0017][0018][0019][0020]
其中，v
gx
为台风场最大梯度风速，k为经验常数，f为科氏系数，其值为f＝2ωsinψ，ω为地球自转角速度，ψ为纬度，为海平面10米高度处的平均最大风速，v
t
为台风移动速度，r
max
为台风最大风速半径，r为观测点到台风中心的距离，v
10,r
为距海平面10m处，与台风中心相距为r处的风速；
[0021]
所述的风速与故障率关系公式为：
[0022][0023][0024][0025]
其中，β
z
为风振系数，μ
z
为风压高度变化系数，μ
s
为体形系数，ρ为空气密度，λ
min
为切入风速时对应的故障率，v
ci
为切入风速，λ
max
为切出风速时对应的故障率，v
co
为切出风速。
[0026]
进一步地，所述的步骤s2中，对风电机组寿命周期内每一单位时间的初始故障率进行修正的表达式为：
[0027][0028]
其中，λ
i
为第i个单位时间的初始故障率，t
i
为第i个单位时间内受台风影响的时间，n
i
为第i个单位时间内的台风数量，t
iks
和t
ike
分别为第i个单位时间内第k个台风的起始时间和终止时间。
[0029]
进一步地，所述的步骤s3中，首先利用风速与功率关系，根据单台风机的初始风速
数据获取单台风机的输出功率，然后通过k-means聚类算法对单台风机输出功率场景聚类到多个离散的功率状态。
[0030]
更进一步地，所述的风速与功率关系的表达式为：
[0031][0032][0033]
其中，p为单台风机的输出功率，p
r
为单台风机的额定功率，v为风速，v
ci
为切入风速，v
r
为额定风速，v
co
为切出风速，a、b和c为模型参数。
[0034]
进一步地，所述的步骤s4中，首先利用马尔可夫链模型表征单台风电机组所有的工作运行状态，然后计算各状态之间的转移概率，所述的各状态之间的转移概率的计算式为：
[0035][0036]
其中，r
ij
为从状态i转移到状态j的转移概率，n
ij
从状态i转移到状态j的次数，t
i
为状态i在整个计算周期内出现的次数。
[0037]
进一步地，所述的步骤s5中，采用状态空间图组合单台机组多状态模型，建立海上风电场的多状态可靠性模型。
[0038]
进一步地，所述的每一单位时间的初始故障率为上一单位时间内修正的故障率，第一个单位时间的初始故障率通过文献数据或实际运行数据获得，所述的单位时间设定为一年。
[0039]
与现有技术相比，本发明具有以下优点：
[0040]
1)本发明的方法考虑台风对海上风电机组输出功率的影响以及风电机组故障率的累积影响，考虑风电机组随机故障下的风电场出力多状态模型，考虑因素多且符合实际风电机组故障情况，通过本发明方法获取的可靠性参数准确，可信度高；
[0041]
2)本发明方法所建立的考虑台风影响的海上风电场寿命期内的多状态可靠性模型，更加符合实际工程中海上风电场的出力模型，为海上风电场接入受端电网的可靠性评估提供了有效的计算参考标准。
附图说明
[0042]
图1为本发明中考虑台风影响海上风电多状态可靠性模型流程图；
[0043]
图2为本发明中不考虑风机随机故障单台风电机组状态空间图；
[0044]
图3为本发明中考虑风机随机故障后单台风电机组状态空间图(未标出非相邻状态之间的转移关系)；
[0045]
图4为本发明中考虑风机随机故障后n台风电机组状态空间图(未标出非相邻状态之间的转移关系)。
具体实施方式
[0046]
下面结合附图和具体实施例对本发明进行详细说明。显然，所描述的实施例是本发明的一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都应属于本发明保护的范围。
[0047]
实施例
[0048]
如图1所示，一种考虑台风影响的海上风电场多状态可靠性获取方法，用于得到海上风电机组寿命周期内每个设定的单位时间内的可靠性参数数据，包括以下步骤：
[0049]
s1：获取海上风电机组的寿命周期，并获取寿命周期内风电机组每一单位时间的初始故障率和单台风机的初始风速数据；
[0050]
s2：获取寿命周期内每一单位时间内的台风关键参数，并根据台风对单台机组故障率的影响，对风电机组寿命周期内每一单位时间的初始故障率进行修正；
[0051]
s3：利用单台风机的初始风速数据和聚类算法对单台风机进行出力状态划分；
[0052]
s4：利用马尔可夫链模型，根据单台风机的出力状态和修正后的初始故障率，建立寿命周期各单位时间的单台机组多状态模型；
[0053]
s5：组合单台机组多状态模型，建立海上风电场的多状态可靠性模型，获取海上风电机组寿命周期内每个单位时间内的可靠性参数数据，用于对海上风电场的可靠性评估。
[0054]
本发明所要解决的主要技术问题是考虑台风对海上风电机组输出功率的影响以及风电机组故障率的累积影响，考虑风电机组随机故障下的风电场出力多状态模型。本发明提出的一种考虑考虑台风影响的海上风电场寿命期内的多状态可靠性模型，主要包括以下几步骤：
[0055]
步骤1：通过海上风电单台机组一年的实测风速数据和风速与功率关系确定单台风机一年的输出功率，通过k-means聚类算法对单台风机出力进行状态划分。
[0056]
风速与功率关系公式为风速与功率关系的表达式为：
[0057][0058]
其中，p为单台风机的输出功率，p
r
为单台风机的额定功率，v为风速，v
ci
为切入风速，v
r
为额定风速，v
co
为切出风速，a、b和c为模型参数，计算公式如下：
[0059][0060]
根据一年风速数据和风速与功率间关系可得到一年单台风机的输出功率，采用k-means聚类算法将一年的单台风机输出功率场景聚类到多个离散的功率状态，其准则函数为：
[0061][0062]
式中，k表示聚类个数；x
i
表示属于第i的样本数据；c
i
表示第i个聚类中心。
[0063]
步骤2：确定海上风电机组第一年初的初始故障率和修复率等可靠性参数，根据batts风场模型以及风速与故障率间数学关系确定台风对单台机组风速和故障率的影响，并根据寿命期内每年台风关键参数的抽样修正风电机组寿命期内各年的初始故障率。
[0064]
确定海上风电机组第一年初的初始故障率和修复率等可靠性参数由文献数据或实际运行数据获得。batts风场模型，采用公式：
[0065][0066][0067][0068][0069]
其中，v
gx
为台风场最大梯度风速，k为经验常数，f为科氏系数，其值为f＝2ωsinψ，ω为地球自转角速度，ψ为纬度，为海平面10米高度处的平均最大风速，v
t
为台风移动速度，r
max
为台风最大风速半径，r为观测点到台风中心的距离，v
10,r
为距海平面10m处，与台风中心相距为r处的风速；
[0070]
台风关键参数包括年台风出现频率、台风移动方向、台风移动速度、台风最小距离和台风中心气压差。其中，年台风出现频率服从泊松分布，台风移动方向服从双正态分布，台风移动速度服从对数正态分布，台风最小距离服从均匀分布，台风中心气压差服从对数正态分布。
[0071]
风速与故障率关系公式为：
[0072]
[0073][0074][0075]
其中，β
z
为风振系数，μ
z
为风压高度变化系数，μ
s
为体形系数，ρ为空气密度，λ
min
为切入风速时对应的故障率，v
ci
为切入风速，λ
max
为切出风速时对应的故障率，v
co
为切出风速。
[0076]
对风电机组寿命周期内每一单位时间的初始故障率进行修正的表达式为：
[0077][0078]
其中，λ
i
为第i个单位时间的初始故障率，t
i
为第i个单位时间内受台风影响的时间，n
i
为第i个单位时间内的台风数量，t
iks
和t
ike
分别为第i个单位时间内第k个台风的起始时间和终止时间。
[0079]
步骤3：利用马尔可夫链模型表征单台风电机组所有的工作运行状态，确定寿命周期内各年的单台机组多状态模型。
[0080]
通过聚类后可以得到单台风电机组出力的markov模型，如图2所示。其中，s
i
表示聚类后的状态，m为聚类状态总数，r
ij
为从状态i到状态j的转移率为：
[0081][0082]
其中，r
ij
为从状态i转移到状态j的转移概率，n
ij
从状态i转移到状态j的次数，t
i
为状态i在整个计算周期内出现的次数。
[0083]
考虑风机随机故障后的单台风电机组状态空间图如图3，其中，λ为单台风机故障率，μ为单台风机修复率。根据状态空间图可以写出考虑风机随机故障的单台风机的转移状态概率矩阵t，求解线性方程组，便可得到图3中各状态对应的稳态概率：
[0084][0085]
其中，p为稳态概率矩阵且p＝{p1,p2,...,p
n
}，n为状态转移图中的状态总数。
[0086]
其中，虚线框中为单台风机出力为0的状态，可以进行状态合并。
[0087]
步骤4：组合单台机组多状态模型，建立海上风电场的多状态可靠性模型。
[0088]
考虑风机随机故障后的n台风电机组状态空间图如图4，其中，n为风机数量，λ为单台风机故障率，μ为单台风机修复率。根据状态空间图可以写出考虑随机故障的n台风机的转移状态概率矩阵t，求解线性方程组，便可得到图4中各状态对应的稳态概率。由于n台风机输出功率状态数较多，为了提高计算效率，同样根据上述k-值聚类算法对n台风机组成的风电场的输出功率进行聚类，确定最优的风电场状态划分数目，建立海上风电场的多状态
可靠性模型。
[0089]
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的工作人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。