本发明属于电网安全稳定控制
技术领域：
，更为具体地讲，涉及一种基于多频段电力系统稳定调节器低频段结构(pss4b-l)抑制超低频振荡的方法。
背景技术：
近年来，实际电网中发生多次频率低于0.1hz的超低频振荡事件，2017年、锦苏直流进行孤岛实验期间，孤岛内部调速器失稳，导致系统出现频率约为0.07hz的频率异常波动情况。云南电网在2016年进行异步联网实验时，网内产生了频率为0.05hz，波动在49.9～50.1hz之间的振荡，振荡持续25分钟，严重威胁了电力系统安全运行。因此就超低频振荡的机理进行分析，提出有效的超低频振荡的控制策略对于电力系统的安全稳定运行具有十分重要的意义。针对超低频振荡的抑制问题，现有研究多是基于优化调速器参数来抑制超低频振荡，但该措施在抑制超低频振荡的同时会引起调速器调速性能的下降。而相关研究表明理论上可通过附加pss抑制超低频振荡，但传统pss在超低频段提供的阻尼较小，无法有效抑制超低频振荡。多频段结构的pss4b通过独立调节各频段的增益、相位、限幅以及滤波器参数，为各个频段的振荡提供合适的阻尼，解决了传统pss低频段抑制振荡能力不足的问题。鉴于此，本专利提出一种基于pss4b-l抑制超低频振荡的方法。技术实现要素：本发明的目的在于克服现有技术的不足，提供一种基于pss4b-l抑制超低频振荡的方法，通过在各种工况情形对pss4b-l各环节参数进行整定以确保其具备鲁棒性，再结合一类多变量阻尼比灵敏度进行pss4b-l的最优选址进行判定，以确定pss4b-l的最优安装位置，最终达到抑制超低频振荡的目的。为实现上述发明目的，本发明一种基于pss4b-l抑制超低频振荡的方法，其特征在于，包括以下步骤：(1)、依据水轮机组的运行工况，确定产生超低频振荡的场景集w，再对场景集内所有场景逐一进行辨识，形成振荡频率集f；其中，wp表示产生超低频振荡的第p个场景，fp表示第p个场景所对应的超低频振荡频率；(2)、基于多变量阻尼比相对灵敏度的pss4b-l选址判定(2.1)、计算第i种超低频振荡场景wi下，第k台发电机组配置pss4b-l的多变量阻尼比相对灵敏度其中，sk_i为场景为wi时第k台发电机组的多变量阻尼比相对灵敏度；ξi为系统的超低频振荡模式阻尼比，vk为第k台发电机组配置的pss4b-l输出，ξi0为ξi的给定值，vk0为变量vk对应的稳态值，r为采样点，θ1,θ2,…,θm为pss4b-l中各环节中的变量；对上式进行变形可得：(2.2)、采用中值法计算(2.3)、对进行离散化处理，再计算各采样点绝对值的平均值：其中，r为pss4b-l的输出曲线上的总采样点数；(2.4)、计算多变量阻尼比相对灵敏度的具体取值ak_i：(2.5)、重复步骤(2.1)-(2.4)依次计算场景集w内其余超低频振荡场景对pss4b-l的多变量阻尼比相对灵敏度的取值；(2.6)、对所有超低频振荡场景下的多变量阻尼比相对灵敏度的取值进行累加求均值处理，得到多种场景下的综合阻尼比灵敏度为：其中，ak-i为在场景i下的多变量阻尼比相对灵敏度的取值；(2.7)、将第k台发电机在多种场景下的综合阻尼比灵敏度与预设的阻尼比灵敏度阈值ψ比较，如果则在第k台发电机组处安装pss4b-l，否则不安装，并进入步骤(2.8)；(2.8)、重复步骤(2.1)-(2.7)，继续判断剩余台发电机组处是否安装pss4b-l，并在满足条件的发电机处安装上pss4b-l；(3)、设置pss4b-l带通环节参数其中，fd为带通环节的中心频率；r为带通的补偿系数，与带通宽度b相关；为了使各种工况情形下的振荡频率均在通频带内，fd和b满足如下关系：其中，fmin与fmax为振荡频率集f中频率的最小值与最大值，δ为频率裕量；(4)、设置pss4b-l超前滞后相位环节以及幅值限幅环节参数在多种工况情形，构建超低频振荡场景集的目标函数：其中，ui为模型中待优化参数，其包含超前滞后相位环节以及幅值限幅环节中所有待优化参数，u为其可行域；α(ui,wi)为在场景wi情形时，系统的待优化参数设置为ui时，超低频振荡模式所对应的实部，ξ'min(ui,wi)为除了超低频振荡模式外的其他振荡模式中阻尼最小值，ξ0为给定阻尼下限值；(5)、根据上述步骤(2)中各发电机组所确定的安装pss4b-l的位置，并结合步骤(3)、(4)设置的各参数，对系统进行模式识别，达到抑制系统的超低频振荡。本发明的发明目的是这样实现的：本发明基于pss4b-l抑制超低频振荡的方法，基于机组的运行工况确定超低频场景集，再计及上述场景集，通过设置pss4b-l带通环节、超前滞后环节、限幅环节相关变量的参数，结合一类多变量阻尼比灵敏度进行pss4b-l的最优选址进行判定，以确定pss4b-l的最优安装位置，最终达到抑制超低频振荡的目的。同时，本发明基于pss4b-l抑制超低频振荡的方法还应具有以下有益效果：(1)、本发明提出的pss4b-l参数设置方法充分考虑系统的多种工况情形，整定所得的pss4b-l参数具备较好的鲁棒性，在多种工况情形下均具有较好的抑制超低频振荡能力；(2)、依据多变量阻尼比灵敏度可识别超低频振荡模式的强相关机组，在灵敏度较大处安装pss4b-l具有更佳的控制效果。附图说明图1是本发明基于pss4b-l抑制超低频振荡的方法流程图；图2是pss4b-l的结构图；图3是西南某区域水电模型图；图4是配置pss4b-l前后系统部分机组的频率偏差的波动曲线；图5场景2情形下配置pss4b-l前后川江边#1机组频率偏差曲线；图6场景4情形下配置pss4b-l前后川江边#1机组频率偏差曲线。具体实施方式下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。实施例图1是本发明基于pss4b-l抑制超低频振荡的方法流程图。在本实施例中，将上述措施应用西南某区域水电进行实例分析，系统等值简化拓扑如图3所示，包含17台水电机组，为更加精细的反映系统机电暂态特性，发电机组均采用五阶模型。水电机组总装机容量约为1200mw，区域内负荷可忽略不计，水电机群产生的功率通过升压接入220kv电网，然后汇集至川九龙220kv变电站，最后经升压至500kv与主网进行连接实现外送。下面我们结合图1和图3，对本发明一种基于pss4b-l抑制超低频振荡的方法进行详细说明，具体包括以下步骤：s1、结合图3所示的水电模型拓扑图，依据水轮机组的运行工况，以不同的水锤效应时间常数值进行表征，通过调整部分水轮机组水锤效应时间常数值，形成如表1所示场景集；场景集川江边#1川穿偏桥#1川五一桥#1川沙坪#1川踏卡#1场景12.02.02.01.01.0场景21.02.03.02.04.0场景34.04.04.04.04.0场景42.52.53.03.04.0场景55.05.04.02.01.0场景65.05.05.05.05.0表1再对场景集内场景逐一进行辨识，形成表2所示的振荡频率集f。表2s2、基于多变量阻尼比灵敏度的pss4b-l选址判定确定产生超低频振荡的场景集w，再对场景集内所有场景逐一进行辨识，形成振荡频率集f；其中，wp表示产生超低频振荡的第p个场景，fp表示第p个场景所对应的超低频振荡频率；(2)、基于多变量阻尼比相对灵敏度的pss4b-l选址判定(2.1)、针对pss4b-l的安装位置的优选，采用一种多变量阻尼比灵敏度计算方法，即阻尼(超低频振荡模式)对pss4b-l输出的灵敏度，具体表现形式如下所示：其中，第i种超低频振荡场景wi下，第k台发电机组配置pss4b-l的多变量阻尼比灵敏度。为消除参数量纲对于灵敏度分析的影响，本发明引入相对灵敏度概念，即：计算第i种超低频振荡场景wi下，第k台发电机组配置pss4b-l的多变量阻尼比相对灵敏度其中，sk_i为场景为wi时第k台发电机组的多变量阻尼比相对灵敏度；ξi为系统的超低频振荡模式阻尼比，vk为第k台发电机组配置的pss4b-l输出，ξi0为ξi的给定值，vk0为变量vk对应的稳态值，r为采样点，θ1,θ2,…,θm为pss4b-l中各环节中的变量；对上式进行变形可得：分析上式可知，vk0/ξi0为常量，阻尼比相对灵敏度的大小与成反比，因此阻尼比相对灵敏度计算的关键在于的计算。s2.2、采用中值法计算针对的计算，为避免引起方向误差，采用中值法进行计算，即可表示为：上式计算中采取的两次摄动，对于pss4b-l各环节的全体变量同时进行扰动，再采用pss4b-l的输出作为全体变量的共同作用结果进行阻尼比灵敏度的计算，其反应了整体pss4b-l对于超低频振荡模式的影响，相较于单变量灵敏度更加合理。s2.3、上式计算所得实为一条曲线，为比较的大小，可对对进行离散化处理，再计算各采样点绝对值的平均值：其中，r为pss4b-l的输出曲线上的总采样点数；s2.4、计算多变量阻尼比相对灵敏度的具体取值ak：s2.5、重复步骤s2.1-s2.4依次计算场景集w内其余超低频振荡场景对pss4b-l的多变量阻尼比相对灵敏度的取值；s2.6、对所有超低频振荡场景下的多变量阻尼比相对灵敏度的取值进行累加求均值处理，得到多种场景下的综合阻尼比灵敏度为：其中，ak-i为在场景i下的多变量阻尼比相对灵敏度的取值；综合阻尼比灵敏度值越大，表明单位控制代价的阻尼效果愈佳，换言之在达到相同阻尼效果情形下，所付出的控制代价越小。即在综合阻尼比灵敏度较大处安装pss4b-l控制效果更佳。s2.7、将第k台发电机在多种场景下的综合阻尼比灵敏度与预设的阻尼比灵敏度阈值ψ比较，如果则在第k台发电机组处安装pss4b-l，否则不安装，并进入步骤s2.8；s2.8、重复步骤s2.1-s2.7，继续判断剩余台发电机组处是否安装pss4b-l，并在满足条件的发电机处安装上pss4b-l；在本实施例中，通过计算表2所示的场景集下系统各机组的综合阻尼比灵敏度，其具体计算结果如下面表3所示；机组综合灵敏度机组综合灵敏度川偏桥#32.959川沙坪#20.951川铁河厂#11.853川踏卡#10.916川江边#11.673川斜卡#30.867川偏桥#11.475川斜卡#10.657川偏桥#21.468川踏卡#20.628川江边#31.459川五一桥#10.565川江边#21.384川斜卡#20.589川铁河厂#21.366川五一桥#20.516川沙坪#11.249阀值1.500表3本实施例中，设置灵敏度阀值为1.5，综合上面表2分析可知，考虑在川偏桥#3、川铁河厂#1、川江边#1处安装pss4b-l。s3、设置pss4b-l带通环节参数ss4b-l带通环节的功能为实现中心频率差分滤波，使得中心频率点出输出幅值最大，相位为零。因此相关参数应当设置为：其中，fd为带通环节的中心频率；r为带通的补偿系数，与带通宽度b相关；为了使各种工况情形下的振荡频率均在通频带内，fd和b满足如下关系：其中，fmin与fmax为振荡频率集f中频率的最小值与最大值，δ为频率裕量,取为0.1；本实施例中，通过不断调节r值使得带通滤波器带宽为b，从而确定r的具体数值，带通环节相关变量的具体取值如表4所示。变量数值变量数值fd0.0379tl13.6399r1.1tl24.0039kl1231tl74.0039kl2231tl84.4043表4s4、设置pss4b-l超前滞后相位环节以及幅值限幅环节参数超前滞后相位环节以及幅值限幅环节的作用是通过设置合理的参数，使得其产生适当的附加转矩，抵消调速系统产生的负阻尼转矩，从而提升整个系统的阻尼，抑制超低频振荡。而且在提升超低频振荡模式时，需要避免恶化其他模式的阻尼。鉴于此，本发明充分考虑到系统存在的多种工况情形，构建超低频振荡场景集的目标函数，具体为：其中，ui为模型中待优化参数，其包含超前滞后相位环节以及幅值限幅环节中所有待优化参数，u为其可行域；α(ui,wi)为在场景wi情形时，系统的待优化参数设置为ui时，超低频振荡模式所对应的实部，ξ'min(ui,wi)为除了超低频振荡模式外的其他振荡模式中阻尼最小值，ξ0为给定阻尼下限值，取为0.06；在本实施例中，超前滞后相位环节以及幅值限幅环节变量的具体取值如表5所示；表5s5、根据上述步骤s2中各发电机组所确定的安装pss4b-l的位置，并结合步骤s3、s4设置的各参数，对系统进行模式识别，就可以达到抑制系统的超低频振荡。图4中给出了按照本发明所述控制策略前后系统部分机组的频率偏差的波动情况，其中包括了实施优化策略前机组的频率偏差曲线(虚线)，实施优化策略后机组的频率偏差曲线(实线)，图示结果证明了本发明所提措施可有效抑制超低频振荡现象。图5、图6给出了在场景2以及场景4情形下按照本发明所提出pss4b-l参数整定方法前后川江边#1机组的频率偏差情形。其中包括了实施策略前机组的频率偏差曲线(虚线)，实施策略后机组的频率偏差曲线(实线)。图示结果表明了本专利所提出的方法具备多种工况情形下的鲁棒性。尽管上面对本发明说明性的具体实施方式进行了描述，以便于本
技术领域：
的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本
技术领域：
的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。当前第1页12