一种基于区域惯性中心的电力系统惯量评估方法及系统与流程

本发明涉及电力系统安全稳定分析领域，尤其涉及一种基于区域惯性中心的电力系统惯量评估方法及系统。

背景技术：

1、惯量评估是一种确定系统的等效惯量的方式，旨在了解系统整体在受到扰动后抵抗频率变化的能力，可以指导电力系统的运行、控制和规划。现有惯量评估主要分为离线评估和在线评估，离线评估主要基于大扰动频率事件，用于事故后系统惯量水平分析；

2、电力系统惯性中心是指当电力系统受到外部扰动时，能反映系统整体惯量特性的等效节点，体现了系统的稳定性和频率响应特性。惯性中心一般通过系统的发电机参数、负荷特性和系统运行状态来确定，可基于惯量中心频率对系统进行频率控制或惯量评估。现有技术中，是直接将系统中所有网络节点作为采样点，根据所有网络节点信息进行惯量评估，然而，电气距离相近的节点在同一时间测到的频率频率偏差和频率变化率差异较小，现有技术直接使用所有节点信息进行惯性评估会带来惯量评估所需频率采样点过多、数据计算量过大的情况，影响评估效率。

技术实现思路

1、有鉴于此，本发明提供一种基于区域惯性中心的电力系统惯量评估方法及系统，通过将基于惯性中心的评估方式应用于分区后的小型区域，可减少惯量评估所需频率采样点，进而优化惯量评估结果，有效提高评估效率，保证电网运维经济性和安全性。

2、本发明实施例解决其技术问题所采用的技术方案是：

3、一种基于区域惯性中心的电力系统惯量评估方法，包括：

4、步骤s1，基于聚类算法进行发电机分区；

5、步骤s2，计及网络节点频率与发电机频率的线性关系，以步骤s1的发电机分区结果为基础进行网络节点分区；

6、步骤s3，基于惯性中心频率评估区域惯量；

7、步骤s4，获取惯量的空间分布特征及进行系统整体惯量误差计算。

8、较优地，所述步骤s1包括：

9、步骤s11，在系统导纳矩阵的基础上考虑发电机暂态电抗形成增广矩阵ys：

10、

11、其中，系统中包含m+n个网络节点，m为发电机端节点数目，n为其他节点数目；g指示发电机内节点，g∈[1,m]；m指示发电机端外节点，m∈[1,m]；n指示系统中的其他节点，n∈[m+1,m+n]；ygg为发电机内节点自导纳；ygm为发电机内节点与发电机外节点互导纳；ymm为发电机外节点自导纳；ymn为发电机外节点和其他节点互导纳；ynn为其他节点自导纳；所述增广矩阵ys中置零位置代表节点间没有直接连接；

12、步骤s12，对所述增广矩阵ys进行kron简化，收缩为m×m的矩阵yg，对所述矩阵yg求逆得到阻抗增广矩阵zg，记阻抗增广矩阵zg中第i行第j列的元素为z(i,j)，进一步建立各发电机内节点间表示电气距离的相关度矩阵d，所述相关度矩阵d中第i行第j列元素d(i,j)，其中：

13、d(i,j)＝z(i,i)+z(j,j)-z(i,j)-z(j,i),i∈[1,m],j∈[1,m]

14、步骤s13，使用k-means聚类方法，预设分区数a，根据所述相关度矩阵d中的电气距离d(i,j)进行发电机分区，形成a个子区域。

15、较优地，所述步骤s2包括：

16、步骤s21，推导系统中网络节点频率与发电机频率的线性关系，其中，网络节点k与发电机j之间的关联系数fk/fj越大表示关联越紧密；

17、步骤s22，根据所述关联系数关联网络节点与发电机，得出各个所述网络节点所在的所述子分区；具体是，为所述网络节点选取出最大关联系数所对应的发电机，令所述最大关联系数所对应的发电机所在的子分区作为所述网络节点对应的子分区；

18、所述步骤s21包括：

19、步骤s211，首先建立系统增广导纳矩阵对应网络方程：

20、

21、其中，in为发电机向系统注入的电流向量，in＝{in-1,in-2,in-3,...,in-m}；ue为发电机内电势节点的电压向量ue＝{ue-1,ue-2,ue-3,...,ue-m}；unet为网络节点电压向量unet＝{unet-1,unet-2,unet-3,...,unet-m,unet-m+1,unet-m+2,unet-m+3,...,unet-m+n}；y1、y2、y3、y4表示增广导纳矩阵ys中的部分矩阵：

22、y1＝ygg

23、y2＝[-ygn 0]

24、y3＝[-ygm 0]t

25、

26、步骤s212，用发电机内电势节点电压表示网络节点电压，即：

27、unet＝rue

28、其中，r表示所述发电机内电势节点电压与所述网络节点电压之间的关联矩阵r：

29、

30、步骤s213，若所述网络节点k处发生故障，将所述关联矩阵r的表达式展开并提出与所述网络节点k相关的部分，即第k行的所有元素，k∈[1,m+n]，将故障点电压表示为各发电机内电势，经简化后将电压角度推广至频率关系，表示为：

31、

32、其中，j∈[1,m]；fk表示所述网络节点k的频率；fj为发电机j的频率；vj为所述发电机j的电压幅值；rkj为所述关联矩阵r中第k行第j列的元素，表示所述网络节点k和所述发电机j之间电压关系；

33、步骤s214，根据步骤s213计算出所述关联系数fk/fj。

34、较优地，所述步骤s3包括：

35、步骤s31，基于惯性中心概念，确定各子区域的惯性中心：

36、

37、其中，令单个子区域中的网络节点总数为m’+n’，所述单个子区域中发电机数目为m’、其他节点数目为n’；d(i,j)表示所述单个子区域中网络节点i与发电机j之间的电气距离；hj为所述单个子区域中所述发电机j的惯性时间常数；sj为所述单个子区域中所述发电机j的额定容量；对应最小f值的节点i作为所述单个子区域的惯性中心，通过本步骤，得到a个fa；

38、步骤s32，在不同位置施加扰动测量惯性中心频率，进而评估各个所述单个子区域惯量ha：

39、

40、其中，△p为所述单个子区域的有功功率缺额；d△fs/dt为所述单个子区域中的一个网络节点扰动发生时刻的频率变化率。

41、较优地，所述步骤s4包括：

42、步骤s41，分析系统惯量的空间分布特征，综合子区域的空间位置与子区域的区域惯量得出系统惯量的地理分布状况，进而确定出系统惯量薄弱环节；

43、步骤s41之后，所述步骤s4还包括：

44、步骤s42，根据各个所述单个子区域惯量ha进行系统惯量的误差计算；

45、所述步骤s42包括：

46、步骤s421，将各个所述单个子区域惯量ha线性相加得到系统惯量hc、并进行误差计算，得出总误差εc：

47、

48、

49、

50、其中，ht表示实际的系统惯量，j表示系统中的发电机，sj是发电机j的额定容量；hj为第j台发电机的惯性时间常数；

51、或者，

52、步骤s422，计算各个所述单个子区域的区域惯量hca、并进行误差计算，得出误差εca：

53、

54、

55、其中，hta表示实际的所述单个子区域的惯量，j表示所述单个子区域中的发电机，sj是发电机j的额定容量；

56、步骤s423，计算出总误差εc：

57、

58、本发明还提供一种基于区域惯性中心的电力系统惯量评估系统，包括：

59、电机分区模块，用于基于聚类算法进行发电机分区；

60、节点分区模块，用于计及网络节点频率与发电机频率的线性关系，以所述电机分区模块的发电机分区结果为基础进行网络节点分区；

61、区域惯量评估模块，用于基于惯性中心频率评估区域惯量；

62、分析模块，用于获取惯量的空间分布特征及进行系统整体惯量误差计算。

63、较优地，所述电机分区模块用于：

64、在系统导纳矩阵的基础上考虑发电机暂态电抗形成增广矩阵ys：

65、

66、其中，系统中包含m+n个网络节点，m为发电机端节点数目，n为其他节点数目；g指示发电机内节点，g∈[1,m]；m指示发电机端外节点，m∈[1,m]；n指示系统中的其他节点，n∈[m+1,m+n]；ygg为发电机内节点自导纳；ygm为发电机内节点与发电机外节点互导纳；ymm为发电机外节点自导纳；ymn为发电机外节点和其他节点互导纳；ynn为其他节点自导纳；所述增广矩阵ys中置零位置代表节点间没有直接连接；

67、对所述增广矩阵ys进行kron简化，收缩为m×m的矩阵yg，对所述矩阵yg求逆得到阻抗增广矩阵zg，记阻抗增广矩阵zg中第i行第j列的元素为z(i,j)，进一步建立各发电机内节点间表示电气距离的相关度矩阵d，所述相关度矩阵d中第i行第j列元素d(i,j)，其中：

68、d(i,j)＝z(i,i)+z(j,j)-z(i,j)-z(j,i),i∈[1,m],j∈[1,m]

69、使用k-means聚类方法，预设分区数a，根据所述相关度矩阵d中的电气距离d(i,j)进行发电机分区，形成a个子区域。

70、较优地，所述节点分区模块用于：

71、推导单元，用于推导系统中网络节点频率与发电机频率的线性关系，其中，网络节点k与发电机j之间的关联系数fk/fj越大表示关联越紧密；

72、节点分区单元，用于根据所述关联系数关联网络节点与发电机，得出各个所述网络节点所在的所述子分区；具体是，为所述网络节点选取出最大关联系数所对应的发电机，令所述最大关联系数所对应的发电机所在的子分区作为所述网络节点对应的子分区；

73、所述节点分区单元具体用于：

74、首先建立系统增广导纳矩阵对应网络方程：

75、

76、其中，in为发电机向系统注入的电流向量，in＝{in-1,in-2,in-3,...,in-m}；ue为发电机内电势节点的电压向量ue＝{ue-1,ue-2,ue-3,...,ue-m}；unet为网络节点电压向量unet＝{unet-1,unet-2,unet-3,...,unet-m,unet-m+1,unet-m+2,unet-m+3,...,unet-m+n}；y1、y2、y3、y4表示增广导纳矩阵ys中的部分矩阵：

77、y1＝ygg

78、y2＝[-ygn 0]

79、y3＝[-ygm 0]t

80、

81、用发电机内电势节点电压表示网络节点电压，即：

82、unet＝rue

83、其中，r表示所述发电机内电势节点电压与所述网络节点电压之间的关联矩阵r：

84、

85、若所述网络节点k处发生故障，将所述关联矩阵r的表达式展开并提出与所述网络节点k相关的部分，即第k行的所有元素，k∈[1,m+n]，将故障点电压表示为各发电机内电势，经简化后将电压角度推广至频率关系，表示为：

86、

87、其中，j∈[1,m]；fk表示所述网络节点k的频率；fj为发电机j的频率；vj为所述发电机j的电压幅值；rkj为所述关联矩阵r中第k行第j列的元素，表示所述网络节点k和所述发电机j之间电压关系；

88、计算出所述关联系数fk/fj。

89、较优地，所述区域惯量评估模块用于：

90、基于惯性中心概念，确定各子区域的惯性中心：

91、

92、其中，令单个子区域中的网络节点总数为m’+n’，所述单个子区域中发电机数目为m’、其他节点数目为n’；d(i,j)表示所述单个子区域中网络节点i与发电机j之间的电气距离；hj为所述单个子区域中所述发电机j的惯性时间常数；sj为所述单个子区域中所述发电机j的额定容量；对应最小f值的节点i作为所述单个子区域的惯性中心，通过本步骤，得到a个fa；

93、在不同位置施加扰动测量惯性中心频率，进而评估各个所述单个子区域惯量ha：

94、

95、其中，△p为所述单个子区域的有功功率缺额；d△fs/dt为所述单个子区域中的一个网络节点扰动发生时刻的频率变化率。

96、较优地，所述分析模块用于：

97、分析系统惯量的空间分布特征，综合子区域的空间位置与子区域的区域惯量得出系统惯量的地理分布状况，进而确定出系统惯量薄弱环节；

98、所述分析模块，还用于根据各个所述单个子区域惯量ha进行系统惯量的误差计算；

99、所述分析模块进行误差计算具体是：

100、将各个所述单个子区域惯量ha线性相加得到系统惯量hc、并进行误差计算，得出总误差εc：

101、

102、

103、

104、其中，ht表示实际的系统惯量，j表示系统中的发电机，sj是发电机j的额定容量；hj为第j台发电机的惯性时间常数；

105、或者，

106、计算各个所述单个子区域的区域惯量hca、并进行误差计算，得出误差εca：

107、

108、

109、其中，hta表示实际的所述单个子区域的惯量，j表示所述单个子区域中的发电机，sj是发电机j的额定容量；

110、计算出总误差εc：

111、

112、由上述技术方案可知，本发明实施例提供的基于区域惯性中心的电力系统惯量评估方法及系统，首先基于聚类算法进行发电机分区；计及网络节点频率与发电机频率的线性关系，以发电机分区结果为基础进行网络节点分区；基于惯性中心频率评估区域惯量；获取惯量的空间分布特征及进行系统整体惯量误差计算。本发明通过将基于惯性中心的评估方式应用于分区后的小型区域，可减少惯量评估所需频率采样点，进而优化惯量评估结果，有效提高电网运维经济性和安全性。