CN103701120A - 一种含风电场的大电网可靠性的评估方法 - Google Patents

一种含风电场的大电网可靠性的评估方法 Download PDF

Info

Publication number
CN103701120A
CN103701120A CN201310716801.3A CN201310716801A CN103701120A CN 103701120 A CN103701120 A CN 103701120A CN 201310716801 A CN201310716801 A CN 201310716801A CN 103701120 A CN103701120 A CN 103701120A
Authority
CN
China
Prior art keywords
wind
energy turbine
wind energy
turbine set
load
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310716801.3A
Other languages
English (en)
Other versions
CN103701120B (zh
Inventor
毛安家
高忠旭
黄昀思
陈得治
宋云亭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power University
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
North China Electric Power University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Electric Power Research Institute Co Ltd CEPRI, North China Electric Power University filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201310716801.3A priority Critical patent/CN103701120B/zh
Publication of CN103701120A publication Critical patent/CN103701120A/zh
Application granted granted Critical
Publication of CN103701120B publication Critical patent/CN103701120B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Control Of Eletrric Generators (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了属于电力系统可靠性评估技术领域的一种含风电场的大电网可靠性的评估方法。该方法首先建立风电场时序功率输出模型;然后将风电场时序功率输出模型和大电网系统结合,建立含风电场的大电网可靠性模型;建立含风电场的大电网可靠性评估的指标体系;最后运用序贯蒙特卡洛模拟法对含风电场的大电网的可靠性进行评估。本发明从电网侧以及风电场侧分别建立了相关可靠性模型及指标体系,用以衡量风电场的接入对大电网可靠性的影响,可以更充分反映系统的整体可靠性状况以及风电场对大电网充裕性的贡献。

Description

一种含风电场的大电网可靠性的评估方法
技术领域
本发明属于电力系统可靠性评估技术领域,特别涉及一种含风电场的大电网可靠性的评估方法。
背景技术
随着世界能源需求的日益增长和可再生资源的开发和利用,风能作为可再生清洁新能源越来越受到人们的重。风力发电由于技术发展成熟并且发电成本相对较低而成为常规能源最主要的一种替代形式。由于风能具有间歇性和随机性,在时间和空间上的分布具有很大的不均匀性,风电场的功率输出随风能资源、时间和空间的变化而波动,大规模风力发电并网后会对电力系统的安全稳定运行带来很大的影响,尤其是对大电网可靠性可靠性产生影响,传统的大电网可靠性评估模型和指标体系无法满足要求,需要提出一种新的含风电场的大电网可靠性的评估方法。
发明内容
针对上述现有技术的存在的问题,本发明提出一种含风电场的大电网可靠性的评估方法,其特征在于,该方法具体步骤如下:
步骤一:建立风电场时序功率输出模型;
步骤二:将步骤一的模型和大电网系统结合,建立含风电场的大电网可靠性模型,包括如下子步骤:
子步骤21)建立元件的序贯概率仿真模型;
子步骤22)建立常规发电机组的序贯概率仿真模型;
子步骤23)利用时序负荷曲线,建立大电网时序负荷模型;
子步骤24)建立含风电场的大电网时序负荷模型;
步骤三:建立含风电场的大电网可靠性评估的指标体系;
步骤四:读入风电场接入大电网后的整个系统的基本参数数据,设置仿真年数m;
步骤五:读入历史风速数据,运用ARMA模型进行模拟,仿真m年的风速数据,基于步骤一中风电场时序功率输出模型,求取风电场功率输出曲线;
步骤六:根据含风电场大电网可靠性模型子步骤21),对大电网中的元件状态持续时间抽样,得到相应的持续时间与故障时间,形成元件的状态序列;
步骤七:根据含风电场大电网可靠性模型子步骤22),形成发电机组的状态序列;
步骤八:将子步骤24)中得到的含风电场的大电网时序负荷看成是负的负荷,再与步骤六和步骤七得到的状态序列相结合,形成含风电场的系统状态序列;
步骤九:按时间顺序依选取步骤八中系统状态序列中的状态,然后依次对选取的每个状态进行判断,如果后评估的系统状态与前面评估过的系统状态一样,则直接读入前面的评估结果;如果不一样,则进入步骤十,对该状态进行分析;
步骤十:判断该状态下系统是否解列,如果不能够解列,则对整个系统进行分析;如果能够解列,则将整个系统分解成各个子系统,进入步骤十一,分别对每个子系统进行分析;
步骤十一:判断各个子系统发电总容量是否满足负荷,如果不满足,则削减负荷;
步骤十二:计算上面选取的系统状态下系统的直流潮流,以此来判断输电线路功率是否越限,若越限,则采用相应的负荷切除策略削减负荷;若不越限,则进入步骤十三;
步骤十三:计算各母线的切负荷量、系统总切负荷量,判断是否满足收敛判据,若满足则进入步骤十四;若不满足,则返回步骤九;
步骤十四:输出系统的可靠性指标。
所述步骤一风电场时序功率输出模型建立步骤如下:
11)录入风电机组的基本参数数据;所述基本参数数据包括风电机组的故障率、修复率,风电场风机台数,仿真年数;
12)统计大量历史风速数据,建立ARMA风速预测模型,并依据仿真年数模拟给出时序风速数据;
13)运用风电机组停运模型来模拟抽样得出风电机组的运行持续时间和故障恢复时间,建立风电机组的正常运行和故障停运的时序过程;
14)基于步骤12)中得到的时序风速数据,依据Jensen尾流损失模型,计算通过尾流效应后通过各风电机组转子处的风速;
15)依据风电机组的功率特性曲线,求取步骤14)中得到的某一风速数值的功率输出,并将风电场内各机组输出功率相叠加,建立风电场的时序功率输出模型。
所述步骤13)中风电机组停运模型为:假设风电机组正常持续时间和故障修复时间服从指数分布,风力发电机组的随机停运用正常运行持续时间和故障修复时间来描述,其中
正常持续运行时间t1为:
t 1 = - 1 λ 1 n ( R ( t ) ) = - 1 λ 1 n ( x ∈ ( 0 , 1 ] )
故障修复时间t2为:
t 2 = - 1 μ 1 nx ( x ∈ ( 0,1 ] )
平均持续工作时间为:
MTTF = 1 λ
平均修复时间为:
MTTR = 1 μ
元件的强迫停运率为:
for = λ λ + μ = MTTR MTTR + MTTF
其中,故障率λ和修复率μ为常数,R(t)为可靠度。
所述步骤15)中风电机组的功率特性曲线由分段函数表示:
P t = 0 0 &le; v < v ci ( A + Bv + Cv 2 ) P r v ci &le; v < v r P r v r &le; v < v co 0 v co &le; v
其中,vci、vr、vco分别表示风电机组的切入风速、额定风速和切出风速,Pr为风电机组的额定功率,A、B、C为风电机组的功率特性曲线参数,利用如下公式计算:
A = 1 ( V ci - V r ) 2 [ V ci ( V ci + V r ) - 4 ( V ci &times; V r ) [ V ci + V r 2 V r ] 3 ] B = 1 ( V ci - V r ) 2 [ 4 ( V ci + V r ) [ V ci + V r 2 V r ] 3 - ( 3 V ci + V r ) C = 1 ( V ci - V r ) 2 [ 2 - 4 [ V ci + V r 2 V r ] 3 ]
所述步骤21)中元件的序贯概率仿真模型为:按照时序,在一个时间跨度上进行模拟,通常假设元件的运行与修复过程服从指数分布:
D i = 1 &lambda; i 1 n R i
式中,Ri为[0,1]均匀分布的抽样值,若λi是某元件的故障率,则Di表示元件运行状态的持续时间;反之,若λi是某元件的修复率,则Di表示元件的故障的持续时间。
所述步骤三中的评估指标体系包括:
(1)切负荷频率
EFLC = 8760 T N
其中,N是有切负荷的状态总数,如果系统状态序列中连续几个系统状态均有切负荷,将其视为一个切负荷状态;
(2)切负荷持续时间
EDLC = &Sigma; i &Element; s t i T &times; 8760
其中,S是有切负荷的系统状态的集合;ti是系统状态i的持续时间,T是总模拟时间;
(3)电量不足期望值
EENS = 8760 T &Sigma; i &Element; s C i t i
其中,Ci是系统状态i的切负荷量;
(4)容量可信度
容量可信度Cc的定义为维持系统可靠性水平不变的情况下,风电场能够替换常规机组的容量与风电场容量的比值,用公式表示为:
C c = C equ C w &times; 100 %
其中,Cequ表示维持系统可靠性水平不变,风电场能够替换常规机组的容量,Cw为风电场的额定容量;
(5)风电场可中断能量效益
风电场可中断能量效益指标WGIEB反映风电场对大电网电量不足期望值的贡献,公式如下:
WGIEB = EENS bw - EENS aw C w
其中,EENSbw、EENSaw分别表示为大电网未加入风电场和加入风电场后的电量不足期望值,Cw为风电场的额定容量。
所述步骤十二中直流潮流表达式如下:
&Sigma; k &Element; &Omega; i A ik P k = P Gi - P Di i &Element; N P k = B k ( &theta; i - &theta; j ) k &Element; M
其中,Ωi为同节点i相连的支路集,Aik为节点i与支路k的关联系数,Pk为支路k潮流,PGi、PDi为节点i发电出力及负荷,N为系统节点集,Bk为支路电纳,θi、θj分别为母线i、j的相角;同时具有两类约束条件分别为发电机功率约束和输电线路潮流约束:
P G min &le; P Gi &le; P G max 0 &le; P k &le; P k max
其中,PGmax和PGmin分别为发电机有功功率输出的上下限,Pkmax为支路传输功率上限。
所述步骤十二中相应的负荷切除策略为采用最优负荷削减模型,其目标函数为最小总负荷削减量:
目标函数:
min &Sigma; i &Element; ND C i
约束条件:
T ( S ) = A ( S ) ( PG - PD + C ) &Sigma; i &Element; NG PG i + &Sigma; i &Element; ND C i = &Sigma; i &Element; ND PD i PG i max &le; PG i &le; PG i max ( i &Element; NG ) 0 &le; C i &le; PD i ( i &Element; NG ) | T k ( S ) | &le; T k max ( k &Element; L )
其中,T(S)表示状态S的有功潮流矢量,A(S)表示状态S的有功潮流和注入功率之间的关系矩阵,PG和PD分别表示发电输出和负荷功率矢量,C是负荷削减矢变量,PGi、PDi、Ci和Tk(S)分别是PG、PD、C和T(S)的元素;PGi min、PGi max和Tk max分别是PGi和Tk(S)的限值,NG、ND和L分别是系统发电母线、负荷母线以及支路的集合。
发明的有益效果:本发明从电网侧以及风电场侧分别建立了相关可靠性模型及指标体系,用以衡量风电场的接入对大电网可靠性的影响,可以更充分反映系统的整体可靠性状况以及风电场对大电网充裕性的贡献。
附图说明
图1为本发明提出的一种含风电场的大电网可靠性的评估方法的流程图;
图2为风电场时序功率输出模型建立方法的流程图;
图3为Jensen尾流损失模型;
图4风电机组的功率特性曲线;
图5为含风电场的大电网可靠性模型建立方法的流程图;
图6为常规发电机组不同运行状态之间的状态转移图;
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
本发明主要提出了一种含风电场的大电网可靠性的评估方法。这种评估方法的具体步骤如图1所示,该方法具体步骤如下:
步骤一:建立风电场时序功率输出模型;
如图2所示为风电场时序功率输出模型建立方法的流程图,具体步骤如下:
11)录入风电机组的基本参数数据;所述基本参数数据包括风电机组的故障率、修复率,风电场风机台数,仿真年数;
12)统计大量历史风速数据,建立ARMA风速预测模型,并依据仿真年数模拟给出时序风速数据;
13)运用风电机组停运模型来模拟抽样得出风电机组的运行持续时间和故障恢复时间,建立风电机组的正常运行和故障停运的时序过程;
风电机组停运模型:一般假设风电机组正常持续时间和故障修复时间服从指数分布,即故障率λ和修复率μ为常数,风力发电机组的随机停运可用正常运行持续时间和故障修复时间来描述。风电机组的故障率λ假定为常数后,则其与可靠度为指数函数关系:
R(t)=e-λt  (t≥0)
利用反函数变换法,得到正常持续运行时间t1:
x = R ( t ) = e - &lambda;t t 1 = - 1 &lambda; 1 n ( R ( t ) ) = - 1 &lambda; 1 nx ( x &Element; ( 0,1 ] )
同理,故障修复时间t2:
t 2 = - 1 &mu; 1 nx ( x &Element; ( 0,1 ] )
值得注意的是,平均持续工作时间MTTF
Figure BDA0000444051750000082
平均修复时间
Figure BDA0000444051750000083
元件的强迫停运率
Figure BDA0000444051750000084
用于计算每一风电机组的输出功率。
14)基于步骤12)中得到的时序风速数据,依据Jensen尾流损失模型,计算通过尾流效应后通过各风电机组转子处的风速;
如图3所示为Jensen尾流损失模型,具体为:
假定风电场是均匀场,X是两台机组之间的距离,R为叶轮半径,Rw为尾流半径,v0表示自然风速,vt表示通过前一个风机转子的风速,vx表示受到尾流效应后通过下一个风机转子的风速。根据动量理论:
&rho;&pi;R w 2 v x = &rho;&pi;R 2 v T + &rho;&pi; ( R w 2 - R 2 ) dR w dt = k w ( &sigma; G + &sigma; 0 ) dR w dx = dR w dt dt dx = k w ( &sigma; G + &sigma; 0 ) / v
其中,ρ表示空气密度,σG、σ0分别为风电机组产生的湍流和自言湍流的均方差,kw表示经验常数,v为平均风速,令K=kwG0)/v,定义为尾流下降系数;v0,vt与风力机组的推力系数CT之间的关系为:
v T = v 0 1 - G T
可推得:
v x = v 0 [ 1 - ( 1 - 1 - C T ) ( R R + KX ) 2 ]
得知自然风速v0,可以算出尾流效应后通过风电机组的转子处风速vx
15)依据风电机组的功率特性曲线,求取步骤14)中得到的某一风速数值的功率输出,并将风电场内各机组输出功率相叠加,建立风电场的时序功率输出模型。
如图4所示为风电机组的功率特性曲线图,该曲线可用如下分段函数表示:
P t = 0 0 &le; v < v ci ( A + Bv + Cv 2 ) P r v ci &le; v < v r P r v r &le; v < v co 0 v co &le; v
其中,vci、vr、vco分别表示风电机组的切入风速、额定风速和切出风速,Pr为风电机组的额定功率,A、B、C为风电机组的功率特性曲线参数,不同类型的风机其参数亦有所区别,可利用如下公式近似计算:
A = 1 ( V ci - V r ) 2 [ V ci ( V ci + V r ) - 4 ( V ci &times; V r ) [ V ci + V r 2 V r ] 3 ] B = 1 ( V ci - V r ) 2 [ 4 ( V ci + V r ) [ V ci + V r 2 V r ] 3 - ( 3 V ci + V r ) C = 1 ( V ci - V r ) 2 [ 2 - 4 [ V ci + V r 2 V r ] 3 ]
将风电场的模型与大电网系统相结合,按照在满足系统约束条件前提下充分利用风电的原则,基于序贯蒙特卡罗模拟法,研究风电接入电网的可靠性。
步骤二:将步骤一的模型和大电网系统结合,建立含风电场的大电网可靠性模型,如图5所示为模型的建立流程图,该模型的建立包括如下子步骤:
子步骤21)建立元件的序贯概率仿真模型;
元件指的是大电网中的元件,包括输电线、变压器等。元件的序贯概率仿真模型为:按照时序,在一个时间跨度上进行模拟,通常假设元件的运行与修复过程服从指数分布:
D i = 1 &lambda; i 1 n R i
式中,一个时间跨度通常指数十年;Ri为[0,1]均匀分布的抽样值,若λi是某元件的故障率,则Di表示元件运行状态的持续时间;反之,若λi是某元件的修复率,则Di表示元件的故障的持续时间。
子步骤22)建立常规发电机组的序贯概率仿真模型;
如图6所示为常规发电机组不同运行状态之间的状态转移图;假设不计机组完全停运状态与降额运行状态之间的状态转移,若常规发电机组初始运行状态处于状态1满额运行,则由状态1可进入状态2降额运行状态或状态3完全停运状态,其转移率分别为λ1和λ2,则状态1的持续时间可分为两种情况进行计算:
TTT 12 = - 1 &lambda; 1 1 n ( U 12 ) TTT 13 = - 1 &lambda; 2 1 n ( U 13 )
其中,U12、U13表示(0,1]之间的随机数。若TTT12小于TTT13,则状态1的持续时间为TTT12,当状态1持续时间结束后进入状态2;反之,若TTT12大于TTT13,则状态1的持续时间为TTT13,当状态1持续时间结束后进入状态3。在进入状态2或状态3后,其返回状态只有状态1,若进入状态2,则将λ1换成μ1即可;若进入状态3,则将λ2换成μ2即可。
子步骤23)利用时序负荷曲线,建立大电网时序负荷模型;
使用序贯蒙特卡罗模拟法时,为方便起见负荷模型利用时序负荷曲线,时间间隔为一小时,则某天的负荷曲线是由二十四个点构成的负荷曲线,假定一年为8760个小时,则年负荷曲线为由8760个负荷点构成的时序负荷曲线,基于此建立时序负荷模型。
子步骤24)建立含风电场的大电网时序负荷模型;
假设风速预测时间间隔为一小时,根据该小时的风速大小及方向,考虑尾流损失、机组故障率等计算出整个风电场在该时刻的功率输出,则风电场的年功率输出曲线由这8760个离散的功率输出值组成。
步骤三:建立含风电场的大电网可靠性评估的指标体系;
旧有的电力系统充裕性评估的指标体系不能够充分体现风电场对于大电网可靠性的贡献,因此,需从风电场角度出发建立相应的指标,衡量风电场接入后对大电网可靠性的影响,本文建立的主要评估指标体系包括:
(1)切负荷频率EFLC(occ/y)
EFLC = 8760 T N
其中,N是有切负荷的状态总数;如果系统状态序列中连续几个系统状态均有切负荷,将其视为一个切负荷状态。
(2)切负荷持续时间EDLC(h/y)
EDLC = &Sigma; i &Element; s t i T &times; 8760
其中,S是有切负荷的系统状态的集合;ti是系统状态i的持续时间,T是总模拟时间。
(3)电量不足期望值EENS(MW·H/y)
EENS = 8760 T &Sigma; i &Element; s C i t i
其中,Ci是系统状态i的切负荷量。
以上三个指标越小越好,越小表示系统越可靠。
(4)容量可信度(Cc
容量可信度的定义为维持系统可靠性水平不变的情况下,风电场可替换常规机组的容量与风电场容量的比值,用公式表示为:
R=f(C,L)=f(C+Cw,L+ΔL)=f(C+Cequ,L+ΔL)
其中,Cequ表示维持系统可靠性水平不变,风电场可替换常规机组的容量,Cw为风电场的额定容量,C、L分别表示原系统的发电容量和峰值负荷,ΔL表示接入风机和常规机组后,系统可额外承载的峰值负荷量。则根据定义,容量可信度的计算公式为:
C c = C equ C w &times; 100 %
容量可信度通常达到35%以上较好。
(5)风电场可中断能量效益(WGIEB)
风电场接入大电网后,可用可中断能量效益指标反映风电场对大电网电量不足期望值的贡献,公式如下:
WGIEB = EENS bw - EENS aw C w
其中,EENSbw、EENSaw分别表示为大电网未加入风电场和加入风电场后的电量不足期望值,Cw为风电场的额定容量。风电场可中断能量效益通常达到30以上较好。
步骤四:读入风电场接入大电网后的整个系统的基本参数数据,设置仿真年数m;
步骤五:读入历史风速数据,运用ARMA模型进行模拟,仿真m年的风速数据,基于步骤一中风电场时序功率输出模型,求取风电场功率输出曲线;
步骤六:根据含风电场大电网可靠性模型子步骤21),对大电网中的元件状态持续时间抽样,得到相应的持续时间与故障时间,形成元件的状态序列;
步骤七:根据含风电场大电网可靠性模型子步骤22),形成发电机组的状态序列;
步骤八:将子步骤24)中得到的含风电场的大电网时序负荷看成是负的负荷,再与步骤六和步骤七得到的状态序列相结合,形成含风电场的系统状态序列;
步骤九:按时间顺序依选取步骤八中系统状态序列中的状态,然后依次对选取的每个状态进行判断,如果后评估的系统状态与前面评估过的系统状态一样,则直接读入前面的评估结果;如果不一样,则进入步骤十,对该状态进行分析;
步骤十:判断该状态下系统是否解列,如果不能够解列,则对整个系统进行分析;如果能够解列,则将整个系统分解成各个子系统,进入步骤十一,分别对每个子系统进行分析;
步骤十一:判断各个子系统发电总容量是否满足负荷,如果不满足,则削减负荷;
步骤十二:计算上面选取的系统状态下系统的直流潮流,以此来判断输电线路功率是否越限,若越限,则采用相应的负荷切除策略削减负荷;若不越限,则进入步骤十三;
直流潮流表达式如下:
&Sigma; k &Element; &Omega; i A ik P k = P Gi - P Di i &Element; N P k = B k ( &theta; i - &theta; j ) k &Element; M
其中,Ωi为同节点i相连的支路集,Aik为节点i与支路k的关联系数,Pk为支路k潮流,PGi、PDi为节点i发电出力及负荷,N为系统节点集,Bk为支路电纳,θi、θj分别为母线i、j的相角;同时具有两类约束条件分别为发电机功率约束和输电线路潮流约束:
P G min &le; P Gi &le; P G max 0 &le; P k &le; P k max
其中,PGmax和PGmin分别为发电机有功功率输出的上下限,Pkmax为支路传输功率上限。
相应的负荷切除策略为采用最优负荷削减模型,其目标函数为最小总负荷削减量:
目标函数:
min &Sigma; i &Element; ND C i
约束条件:
T ( S ) = A ( S ) ( PG - PD + C ) &Sigma; i &Element; NG PG i + &Sigma; i &Element; ND C i = &Sigma; i &Element; ND PD i PG i max &le; PG i &le; PG i max ( i &Element; NG ) 0 &le; C i &le; PD i ( i &Element; NG ) | T k ( S ) | &le; T k max ( k &Element; L )
其中,T(S)表示状态S的有功潮流矢量,A(S)表示状态S的有功潮流和注入功率之间的关系矩阵,PG和PD分别表示发电输出和负荷功率矢量,C是负荷削减矢变量,PGi、PDi、Ci和Tk(S)分别是PG、PD、C和T(S)的元素;PGi min、PGi max和Tk max分别是PGi和Tk(S)的限值,NG、ND和L分别是系统发电母线、负荷母线以及支路的集合。
步骤十三:计算各母线的切负荷量、系统总切负荷量,判断是否满足收敛判据,若满足则进入步骤十四;若不满足,则返回步骤九;
步骤十四:输出系统的可靠性指标。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (8)

1.一种含风电场的大电网可靠性的评估方法,其特征在于,该方法具体步骤如下:
步骤一:建立风电场时序功率输出模型;
步骤二:将步骤一的模型和大电网系统结合,建立含风电场的大电网可靠性模型,包括如下子步骤:
子步骤21)建立元件的序贯概率仿真模型;
子步骤22)建立常规发电机组的序贯概率仿真模型;
子步骤23)利用时序负荷曲线,建立大电网时序负荷模型;
子步骤24)建立含风电场的大电网时序负荷模型;
步骤三:建立含风电场的大电网可靠性评估的指标体系;
步骤四:读入风电场接入大电网后的整个系统的基本参数数据,设置仿真年数m;
步骤五:读入历史风速数据,运用ARMA模型进行模拟,仿真m年的风速数据,基于步骤一中风电场时序功率输出模型,求取风电场功率输出曲线;
步骤六:根据含风电场大电网可靠性模型子步骤21),对大电网中的元件状态持续时间抽样,得到相应的持续时间与故障时间,形成元件的状态序列;
步骤七:根据含风电场大电网可靠性模型子步骤22),形成发电机组的状态序列;
步骤八:将子步骤24)中得到的含风电场的大电网时序负荷看成是负的负荷,再与步骤六和步骤七得到的状态序列相结合,形成含风电场的系统状态序列;
步骤九:按时间顺序依选取步骤八中系统状态序列中的状态,然后依次对选取的每个状态进行判断,如果后评估的系统状态与前面评估过的系统状态一样,则直接读入前面的评估结果;如果不一样,则进入步骤十,对该状态进行分析;
步骤十:判断该状态下系统是否解列,如果不能够解列,则对整个系统进行分析;如果能够解列,则将整个系统分解成各个子系统,进入步骤十一,分别对每个子系统进行分析;
步骤十一:判断各个子系统发电总容量是否满足负荷,如果不满足,则削减负荷;
步骤十二:计算上面选取的系统状态下系统的直流潮流,以此来判断输电线路功率是否越限,若越限,则采用相应的负荷切除策略削减负荷;若不越限,则进入步骤十三;
步骤十三:计算各母线的切负荷量、系统总切负荷量,判断是否满足收敛判据,若满足则进入步骤十四;若不满足,则返回步骤九;
步骤十四:输出系统的可靠性指标。
2.根据权利要求1所述的评估方法,其特征在于,所述步骤一风电场时序功率输出模型建立步骤如下:
11)录入风电机组的基本参数数据;所述基本参数数据包括风电机组的故障率、修复率,风电场风机台数,仿真年数;
12)统计大量历史风速数据,建立ARMA风速预测模型,并依据仿真年数模拟给出时序风速数据;
13)运用风电机组停运模型来模拟抽样得出风电机组的运行持续时间和故障恢复时间,建立风电机组的正常运行和故障停运的时序过程;
14)基于步骤12)中得到的时序风速数据,依据Jensen尾流损失模型,计算通过尾流效应后通过各风电机组转子处的风速;
15)依据风电机组的功率特性曲线,求取步骤14)中得到的某一风速数值的功率输出,并将风电场内各机组输出功率相叠加,建立风电场的时序功率输出模型。
3.根据权利要求2所述的评估方法,其特征在于,所述步骤13)中风电机组停运模型为:假设风电机组正常持续时间和故障修复时间服从指数分布,风力发电机组的随机停运用正常运行持续时间和故障修复时间来描述,其中
正常持续运行时间t1为:
t 1 = - 1 &lambda; 1 n ( R ( t ) ) = - 1 &lambda; 1 n ( x &Element; ( 0 , 1 ] )
故障修复时间t2为:
t 2 = - 1 &mu; 1 nx ( x &Element; ( 0,1 ] )
平均持续工作时间为:
MTTF = 1 &lambda;
平均修复时间为:
MTTR = 1 &mu;
元件的强迫停运率为:
for = &lambda; &lambda; + &mu; = MTTR MTTR + MTTF
其中,故障率λ和修复率μ为常数,R(t)为可靠度。
4.根据权利要求2所述的评估方法,其特征在于,所述步骤15)中风电机组的功率特性曲线由分段函数表示:
P t = 0 0 &le; v < v ci ( A + Bv + Cv 2 ) P r v ci &le; v < v r P r v r &le; v < v co 0 v co &le; v
其中,vci、vr、vco分别表示风电机组的切入风速、额定风速和切出风速,Pr为风电机组的额定功率,A、B、C为风电机组的功率特性曲线参数,利用如下公式计算:
A = 1 ( V ci - V r ) 2 [ V ci ( V ci + V r ) - 4 ( V ci &times; V r ) [ V ci + V r 2 V r ] 3 ] B = 1 ( V ci - V r ) 2 [ 4 ( V ci + V r ) [ V ci + V r 2 V r ] 3 - ( 3 V ci + V r ) C = 1 ( V ci - V r ) 2 [ 2 - 4 [ V ci + V r 2 V r ] 3 ]
5.根据权利要求1所述的评估方法,其特征在于,所述步骤21)中元件的序贯概率仿真模型为:按照时序,在一个时间跨度上进行模拟,通常假设元件的运行与修复过程服从指数分布:
D i = 1 &lambda; i 1 n R i
式中,Ri为[0,1]均匀分布的抽样值,若λi是某元件的故障率,则Di表示元件运行状态的持续时间;反之,若λi是某元件的修复率,则Di表示元件的故障的持续时间。
6.根据权利要求1所述的评估方法,其特征在于,所述步骤三中的评估指标体系包括:
(1)切负荷频率
EFLC = 8760 T N
其中,N是有切负荷的状态总数,如果系统状态序列中连续几个系统状态均有切负荷,将其视为一个切负荷状态;
(2)切负荷持续时间
EDLC = &Sigma; i &Element; s t i T &times; 8760
其中,S是有切负荷的系统状态的集合;ti是系统状态i的持续时间,T是总模拟时间;
(3)电量不足期望值
EENS = 8760 T &Sigma; i &Element; s C i t i
其中,Ci是系统状态i的切负荷量;
(4)容量可信度
容量可信度Cc的定义为维持系统可靠性水平不变的情况下,风电场能够替换常规机组的容量与风电场容量的比值,用公式表示为:
C c = C equ C w &times; 100 %
其中,Cequ表示维持系统可靠性水平不变,风电场能够替换常规机组的容量,Cw为风电场的额定容量;
(5)风电场可中断能量效益
风电场可中断能量效益指标WGIEB反映风电场对大电网电量不足期望值的贡献,公式如下:
WGIEB = EENS bw - EENS aw C w
其中,EENSbw、EENSaw分别表示为大电网未加入风电场和加入风电场后的电量不足期望值,Cw为风电场的额定容量。
7.根据权利要求1所述的评估方法,其特征在于,所述步骤十二中直流潮流表达式如下:
&Sigma; k &Element; &Omega; i A ik P k = P Gi - P Di i &Element; N P k = B k ( &theta; i - &theta; j ) k &Element; M
其中,Ωi为同节点i相连的支路集,Aik为节点i与支路k的关联系数,Pk为支路k潮流,PGi、PDi为节点i发电出力及负荷,N为系统节点集,Bk为支路电纳,θi、θj分别为母线i、j的相角;同时具有两类约束条件分别为发电机功率约束和输电线路潮流约束:
P G min &le; P Gi &le; P G max 0 &le; P k &le; P k max
其中,PGmax和PGmin分别为发电机有功功率输出的上下限,Pkmax为支路传输功率上限。
8.根据权利要求1所述的评估方法,其特征在于,所述步骤十二中相应的负荷切除策略为采用最优负荷削减模型,其目标函数为最小总负荷削减量:
目标函数:
min &Sigma; i &Element; ND C i
约束条件:
T ( S ) = A ( S ) ( PG - PD + C ) &Sigma; i &Element; NG PG i + &Sigma; i &Element; ND C i = &Sigma; i &Element; ND PD i PG i max &le; PG i &le; PG i max ( i &Element; NG ) 0 &le; C i &le; PD i ( i &Element; NG ) | T k ( S ) | &le; T k max ( k &Element; L )
其中,T(S)表示状态S的有功潮流矢量,A(S)表示状态S的有功潮流和注入功率之间的关系矩阵,PG和PD分别表示发电输出和负荷功率矢量,C是负荷削减矢变量,PGi、PDi、Ci和Tk(S)分别是PG、PD、C和T(S)的元素;PGi min、PGi max和Tk max分别是PGi和Tk(S)的限值,NG、ND和L分别是系统发电母线、负荷母线以及支路的集合。
CN201310716801.3A 2013-12-23 2013-12-23 一种含风电场的大电网可靠性的评估方法 Expired - Fee Related CN103701120B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310716801.3A CN103701120B (zh) 2013-12-23 2013-12-23 一种含风电场的大电网可靠性的评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310716801.3A CN103701120B (zh) 2013-12-23 2013-12-23 一种含风电场的大电网可靠性的评估方法

Publications (2)

Publication Number Publication Date
CN103701120A true CN103701120A (zh) 2014-04-02
CN103701120B CN103701120B (zh) 2015-08-12

Family

ID=50362565

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310716801.3A Expired - Fee Related CN103701120B (zh) 2013-12-23 2013-12-23 一种含风电场的大电网可靠性的评估方法

Country Status (1)

Country Link
CN (1) CN103701120B (zh)

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103996079A (zh) * 2014-05-29 2014-08-20 东南大学 一种基于条件概率的风电功率加权预测方法
CN104167734A (zh) * 2014-08-29 2014-11-26 重庆大学 基于时序模拟的并网型微网可靠性评估方法
CN104504606A (zh) * 2014-08-14 2015-04-08 贵州电力试验研究院 一种多类型电源分散并网的可信容量计算方法
CN105162111A (zh) * 2015-08-25 2015-12-16 国家电网公司 一种计及电网可靠性的电气设备检修运行方式的确定方法
CN105226650A (zh) * 2015-10-19 2016-01-06 重庆大学 基于微燃机-储能联合运行策略的微电网可靠性计算方法
CN105429129A (zh) * 2015-10-14 2016-03-23 中国电力科学研究院 一种考虑网络约束的间歇性能源发电容量置信度评估方法
CN105427005A (zh) * 2016-01-18 2016-03-23 重庆大学 一种风电场运行风险评估方法
CN105449667A (zh) * 2015-11-05 2016-03-30 国家电网公司 一种发输电系统可靠性预测方法
CN107035618A (zh) * 2017-06-12 2017-08-11 电子科技大学 一种基于尾流效应的风电场优化控制方法
CN107591833A (zh) * 2016-07-08 2018-01-16 华北电力大学(保定) 一种计及不同运行策略的微网可靠性评估方法
CN107633138A (zh) * 2017-09-21 2018-01-26 河海大学 基于实际出力特性的发电机频率模型建模方法
CN107742898A (zh) * 2017-10-17 2018-02-27 国网江苏省电力公司盐城供电公司 一种用于光伏多级接入的虚拟电厂源网协同能量管理系统
CN108053149A (zh) * 2018-01-05 2018-05-18 东南大学 一种考虑气象因素的光伏电站并网风险评估方法
CN108335004A (zh) * 2017-09-07 2018-07-27 广东石油化工学院 一种基于受阻电能相等的风力发电系统可靠性评价方法
CN108932568A (zh) * 2017-05-22 2018-12-04 北京金风科创风电设备有限公司 风力发电机组可靠性指标的分配方法和装置
CN109217294A (zh) * 2018-09-14 2019-01-15 大唐新疆清洁能源有限公司 风电场大型电网安全性评估方法
CN110472371A (zh) * 2019-09-06 2019-11-19 西安交通大学 一种基于恢复力的电力系统元件重要度的评估方法
CN111697572A (zh) * 2020-06-15 2020-09-22 西安交通大学 基于多阶段随机规划理论的电源及电力流结构优化方法
CN112380695A (zh) * 2020-11-13 2021-02-19 国网上海市电力公司 一种考虑实时风速的海上风电场集电系统可靠性获取方法
CN112434468A (zh) * 2020-12-02 2021-03-02 深圳供电局有限公司 一种配电网可靠性评估方法
CN113489066A (zh) * 2021-07-08 2021-10-08 华翔翔能科技股份有限公司 计及供需不确定性的含储能电网区间供电可靠性评估方法
WO2022016622A1 (zh) * 2020-07-22 2022-01-27 南京东博智慧能源研究院有限公司 真双极柔性直流输电系统故障下的自适应优化调控方法
CN115510677A (zh) * 2022-10-17 2022-12-23 华能山东发电有限公司 一种风电场发电能力评估方法及系统
CN116703249A (zh) * 2023-08-07 2023-09-05 南京师范大学 一种基于ckl风电容量预测的可信度分析方法
CN117277582A (zh) * 2023-10-24 2023-12-22 国网安徽省电力有限公司 一种基于大数据的电力企业运营监测分析系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101413981A (zh) * 2008-11-27 2009-04-22 上海交通大学 电力系统运行备用可靠性测试系统
CN101799681A (zh) * 2010-02-10 2010-08-11 刘文祥 智能电网
CN103078343A (zh) * 2012-11-05 2013-05-01 中国电力科学研究院 一种风电并网对电网暂态及中长期可靠性影响评估方法
CN103248051A (zh) * 2013-05-25 2013-08-14 南京南瑞集团公司 一种评估风电场功率波动影响电网运行安全风险的方法
US20130268131A1 (en) * 2012-04-09 2013-10-10 Clemson University Method and System for Dynamic Stochastic Optimal Electric Power Flow Control

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101413981A (zh) * 2008-11-27 2009-04-22 上海交通大学 电力系统运行备用可靠性测试系统
CN101799681A (zh) * 2010-02-10 2010-08-11 刘文祥 智能电网
US20130268131A1 (en) * 2012-04-09 2013-10-10 Clemson University Method and System for Dynamic Stochastic Optimal Electric Power Flow Control
CN103078343A (zh) * 2012-11-05 2013-05-01 中国电力科学研究院 一种风电并网对电网暂态及中长期可靠性影响评估方法
CN103248051A (zh) * 2013-05-25 2013-08-14 南京南瑞集团公司 一种评估风电场功率波动影响电网运行安全风险的方法

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103996079A (zh) * 2014-05-29 2014-08-20 东南大学 一种基于条件概率的风电功率加权预测方法
CN103996079B (zh) * 2014-05-29 2017-02-15 东南大学 一种基于条件概率的风电功率加权预测方法
CN104504606A (zh) * 2014-08-14 2015-04-08 贵州电力试验研究院 一种多类型电源分散并网的可信容量计算方法
CN104167734B (zh) * 2014-08-29 2016-04-27 重庆大学 基于时序模拟的并网型微网可靠性评估方法
CN104167734A (zh) * 2014-08-29 2014-11-26 重庆大学 基于时序模拟的并网型微网可靠性评估方法
CN105162111A (zh) * 2015-08-25 2015-12-16 国家电网公司 一种计及电网可靠性的电气设备检修运行方式的确定方法
CN105429129B (zh) * 2015-10-14 2020-07-07 中国电力科学研究院 一种考虑网络约束的间歇性能源发电容量置信度评估方法
CN105429129A (zh) * 2015-10-14 2016-03-23 中国电力科学研究院 一种考虑网络约束的间歇性能源发电容量置信度评估方法
CN105226650A (zh) * 2015-10-19 2016-01-06 重庆大学 基于微燃机-储能联合运行策略的微电网可靠性计算方法
CN105449667A (zh) * 2015-11-05 2016-03-30 国家电网公司 一种发输电系统可靠性预测方法
CN105449667B (zh) * 2015-11-05 2018-03-20 国家电网公司 一种发输电系统可靠性预测方法
CN105427005A (zh) * 2016-01-18 2016-03-23 重庆大学 一种风电场运行风险评估方法
CN107591833A (zh) * 2016-07-08 2018-01-16 华北电力大学(保定) 一种计及不同运行策略的微网可靠性评估方法
CN108932568A (zh) * 2017-05-22 2018-12-04 北京金风科创风电设备有限公司 风力发电机组可靠性指标的分配方法和装置
CN107035618A (zh) * 2017-06-12 2017-08-11 电子科技大学 一种基于尾流效应的风电场优化控制方法
CN107035618B (zh) * 2017-06-12 2019-03-01 电子科技大学 一种基于尾流效应的风电场优化控制方法
CN108335004A (zh) * 2017-09-07 2018-07-27 广东石油化工学院 一种基于受阻电能相等的风力发电系统可靠性评价方法
CN107633138A (zh) * 2017-09-21 2018-01-26 河海大学 基于实际出力特性的发电机频率模型建模方法
CN107633138B (zh) * 2017-09-21 2021-01-29 河海大学 基于实际出力特性的发电机频率模型建模方法
CN107742898A (zh) * 2017-10-17 2018-02-27 国网江苏省电力公司盐城供电公司 一种用于光伏多级接入的虚拟电厂源网协同能量管理系统
CN108053149A (zh) * 2018-01-05 2018-05-18 东南大学 一种考虑气象因素的光伏电站并网风险评估方法
CN109217294A (zh) * 2018-09-14 2019-01-15 大唐新疆清洁能源有限公司 风电场大型电网安全性评估方法
CN110472371A (zh) * 2019-09-06 2019-11-19 西安交通大学 一种基于恢复力的电力系统元件重要度的评估方法
CN111697572A (zh) * 2020-06-15 2020-09-22 西安交通大学 基于多阶段随机规划理论的电源及电力流结构优化方法
WO2022016622A1 (zh) * 2020-07-22 2022-01-27 南京东博智慧能源研究院有限公司 真双极柔性直流输电系统故障下的自适应优化调控方法
CN112380695A (zh) * 2020-11-13 2021-02-19 国网上海市电力公司 一种考虑实时风速的海上风电场集电系统可靠性获取方法
CN112380695B (zh) * 2020-11-13 2024-09-06 国网上海市电力公司 一种考虑实时风速的海上风电场集电系统可靠性获取方法
CN112434468A (zh) * 2020-12-02 2021-03-02 深圳供电局有限公司 一种配电网可靠性评估方法
CN113489066A (zh) * 2021-07-08 2021-10-08 华翔翔能科技股份有限公司 计及供需不确定性的含储能电网区间供电可靠性评估方法
CN115510677A (zh) * 2022-10-17 2022-12-23 华能山东发电有限公司 一种风电场发电能力评估方法及系统
CN115510677B (zh) * 2022-10-17 2024-01-09 华能山东发电有限公司 一种风电场发电能力评估方法及系统
CN116703249A (zh) * 2023-08-07 2023-09-05 南京师范大学 一种基于ckl风电容量预测的可信度分析方法
CN116703249B (zh) * 2023-08-07 2024-01-19 南京师范大学 一种基于ckl风电容量预测的可信度分析方法
CN117277582A (zh) * 2023-10-24 2023-12-22 国网安徽省电力有限公司 一种基于大数据的电力企业运营监测分析系统

Also Published As

Publication number Publication date
CN103701120B (zh) 2015-08-12

Similar Documents

Publication Publication Date Title
CN103701120B (zh) 一种含风电场的大电网可靠性的评估方法
Wen et al. A review on reliability assessment for wind power
Billinton et al. Multistate wind energy conversion system models for adequacy assessment of generating systems incorporating wind energy
CN103138267B (zh) 一种基于相对动能的电力系统暂态稳定实时紧急控制方法
CN103746370B (zh) 一种风电场可靠性建模方法
CN102831321B (zh) 一种基于蒙特卡洛方法的风电场风险评估方法
CN104682381B (zh) 大型风电场柔性直流输电系统可靠性计算方法
CN103904695B (zh) 一种基于mcs-pso的临近海岛多微网动态调度方法
CN105429129A (zh) 一种考虑网络约束的间歇性能源发电容量置信度评估方法
CN104077664B (zh) 一种风电储能发电系统的置信容量评估方法
CN104779611A (zh) 基于集中式和分布式双层优化策略的微电网经济调度方法
Mabel et al. Adequacy evaluation of wind power generation systems
CN102915471A (zh) 一种风电扰动在线安全预警的方法
CN105741025A (zh) 基于风电波动在线风险评估的预防控制方法
CN105226650A (zh) 基于微燃机-储能联合运行策略的微电网可靠性计算方法
CN109217294A (zh) 风电场大型电网安全性评估方法
CN101923685A (zh) 一种基于断线故障概率预测的电力切负荷决策系统及方法
Nivedh et al. Repowering of wind farms—a case study
CN105514990A (zh) 综合经济性与安全性的输电线路利用率提高平台及方法
Roy et al. A reliability assessment model of a wind farm for generation adequacy studies of wind integrated power system
Kyle et al. The transitional states of a floating wind turbine during high levels of surge
CN105703395A (zh) 一种风电消纳能力分析方法
CN102570447B (zh) 一种电网发展阶段划分方法
CN103886514A (zh) 一种考虑经济性与安全性的电网发展评估方法
Jiang et al. Research on nodal wind power values and optimal accommodation based on locational marginal price

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150812

Termination date: 20151223

EXPY Termination of patent right or utility model