[发明专利]一种基于轨迹灵敏度的风电配置比例快速计算方法

权利要求书

1.一种基于轨迹灵敏度的风电配置比例快速计算方法，其特征在于，所述方法包括如下步骤：

1）置迭代计算次数i=0，设置风电配置比例初值β₀，且0≤β≤1；

2）根据风电配置比例β_i生成运行方式，并采用潮流计算方法对含风电的系统进行潮流计算；

3）设置考核故障，进行暂态仿真，基于轨迹计算系统状态量x的轨迹灵敏度；

4）选定待分析状态变量x_j，并确定与x_j相关的n个轨迹灵敏度S₁、S₂…S_n，根据这n个轨迹灵敏度各时刻的瞬时值计算风电配置比例β_i下瞬时轨迹灵敏度范数TSN(t)；

5）判断所述迭代计算次数i是否等于0：若等于0则进行步骤6)；否则判断步骤4)中求取的TSN_i幅值是否发生突变，是则结束计算方法，否则进行步骤6)；

6）根据所述轨迹灵敏度调整风电配置比例β_i+1；

7）判断调整后的风电配置比例β_i+1是否大于1，是则进行步骤8），否则更新迭代计算次数i＝i+1并转向步骤2)；

8）计算结束，输出计算结果。

2.如权利要求1所述的快速计算方法，其特征在于，步骤1）所述风电配置比例初值β₀为0.5。

3.如权利要求1所述的快速计算方法，其特征在于，步骤3）所述计算系统状态量的轨迹灵敏度包括：

电力系统的动态过程用如下一组微分代数方程表示：

X·=F(X,y)---(1)]]>

0=g-(X,y)s(X,y)<0g+(X,y)s(X,y)>0---(2)]]>

其中X=xβ,]]>F=f0]]>

式中：x为系统动态状态量，y为系统代数量，β为风电配置比例，X为扩展后的状态量，F为扩展后的微分方程，g-、g+为系统运行模式变化前后的代数方程，s为系统运行模式的变化条件；当s(X,y)＝0时，系统运行模式发生变化；

将式(1)和式(2)对风电配置比例β求导得：

X·β=FXXβ+Fyyβ]]>

0=gX-Xβ+gy-yβgX+Xβ+gy+yβ---(3)]]>

式中：X_β、y_β即为系统轨迹灵敏度，F_X、F_y、为相应导函数；

求解式(1)-(3)，得到系统状态量的轨迹灵敏度。