基于机组出力调整的小干扰稳定辅助决策计算W.docVIP

基于机组出力调整的小干扰稳定辅助决策计算 鲍颜红1,2，徐伟2，徐泰山2，郑伟2，方勇杰2 （1．华北电力大学电气与电子工程学院，北京市 102206； 2．国网电力科学研究院，南京市 210003） 摘要：关键词： 中图分类号：TM7120 引言Arnoldi 算法 (IRAM) 可使用基于稀疏技术的代数方程求解而应用到非常大型的系统[6]。通过IRAM算法可以提供给定阻尼比和频率范围和给定点附近的两种特征值计算方式。无论采用哪种计算方式，计算结果当中都会包含多个模式。因此，需要设计特征值匹配算法确定机组出力调整后与感兴趣的模式相对应的模式。另外，虽然IRAM算法可以充分利用矩阵的稀疏性，但特征值计算所用时间及内存与状态矩阵的维数密切相关，随着互联电网的规模不断增大，IRAM算法的计算速度仍有可能无法满足在线应用的要求。 1.2 灵敏度计算技术 根据模式阻尼对机组有功的灵敏度可以确定机组的出力调整性能和调整方向。状态矩阵中每个元素都可以表示为节点电压和注入电流的函数，因此计算模式阻尼对机组有功的灵敏度需要先求出节点电压和注入电流对运行参数（机组有功）的偏导数，再求取状态矩阵元素对节点电压和注入电流的偏导数，从而得到状态矩阵元素对运行参数的偏导数，最后计算出特征值对参数的灵敏度[7]。另外一种求取灵敏度的方法是摄动法，该方法不需要求取上述的大量导数，仅针对每一个运行参数求取摄动后的特征值。无论采用哪一种方法都需要额外的计算时间。 对于单个模式，需要计算模式阻尼对多个机组有功的灵敏度，若存在多个弱阻尼模式则计算量更大，因此，通过灵敏度来计算小干扰稳定辅助决策并不适用于在线应用。 1.3 并行计算技术 随着集群技术的发展，用高性能计算技术解决大规模电力系统的计算问题，被认为是实现大规模、复杂电网在线分析计算的有效方法。可以将给定范围的特征值计算转化为不同移位点的特征值计算，然后将计算任务下发到不同CPU（对于多核CPU则是不同的核）进行计算，然后将计算结果进行合并。 在小干扰稳定辅助决策计算中包含多个模式时，可以按机组出力调整方案来并行计算控制措施[8]，也可以按模式来并行计算控制措施。前者在每次迭代过程中都需要进行给定范围的特征值计算；后者在计算单个模式控制措施的迭代过程中需要进行给定点的特征值计算。与给定范围的特征值计算相比，给定点的特征值计算相对较快，本文采用按模式并行的计算方案。此时，需要在计算各模式的控制措施前确定各机组的出力调整方向，在得到各模式的控制措施后还需要将控制措施进行合并，然后再对合并后的控制措施进行校核。 2 单模式的辅助决策计算 2.1 确定参与控制的机组 参与因子反映了状态变量与特征根之间的关联程度。高的参与因子是机组出力调整有效提高模式阻尼的必要条件。参与因子虽然无法精确度量模式阻尼对机组有功出力的灵敏度，但参与因子可由特征值计算结果直接给出，不需要额外的计算量。因此，可将参与因子作为机组出力调整的近似性能指标，设定参与因子门槛值，选取参与因子大于的机组参与控制。 2.2 确定机组出力调整方向 参与因子无法给出机组出力调整方向的信息。参与区域振荡的机组可以分为两群，对应于实际系统的送端区域和受端区域。一般情况下，将送端区域的机组减出力、受端区域的机组增出力可以提高模式的阻尼。因此，可将参与控制的机组分为两群，分别确定两群机组的出力调整方向。 对于某一个复数特征根，其模态向量为，在状态变量和上观察相应的过渡过程时，其振幅大小的比值等于，相对的振荡相位差为。因此，可以根据模态的相角来划分相互振荡的两群。以模态向量中幅值最大的元素作为参考，计算各机组转子角（或转速）对应的模态相角；设置角度门槛值，的绝对值在范围内的机组为S群，的绝对值在范围内的机组为A群。仿真实验表明，将设为可保证分群结果的正确性。 将给定的功率摄动量（）平均分摊到S群机组，按参与因子由小到大的顺序在A群中寻找出力可增加的机组平衡出力，计算摄动后的模式阻尼；采用同样的方法计算A群机组功率摄动后的模式阻尼；根据摄动前后的模式阻尼、和确定两群机组的出力调整方向。 2.2 特征值匹配算法 机组出力调整措施实施前后同一个振荡模式的机组分群情况，特别是参与因子较大机组的分群情况基本保持一致。因此，采用如下的特征值匹配算法来确定控制后的模式阻尼。 设定频率门槛值，待匹配模式的频率为，根据模式的频率筛选出的模式；设待匹配模式的模态向量中幅值最大的元素对应的机组为，将机组对应的模态相角设为0，重新计算其它机组对应的模态相角；根据模态相角将模式中参与因子大于的机组分成两群；确定筛选出的模式集中与待匹配模式的分群情况相同的模式，即控制后模式。 2.3 计算流程 考虑计算速度和按模式并行的计算方案，本文采用灵敏度法确定机组调整量，具体步骤如下： 1）将机组出力调