CN112615380A - 一种基于母线电压历史数据分析的逆调压控制方法 - Google Patents

Abstract

本发明属于电力系统自动电压控制技术领域，特别涉及一种基于母线电压历史数据分析的逆调压控制方法的方法。本方法在电压控制中，对于投入自动电压控制系统闭环控制的电厂和变电站的母线电压值进行统计，以最近10～30日的电压历史数据作为统计分析的依据，结合电厂机组、变电站容抗器等无功设备的运行情况，自动计算母线电压上限，一方面满足了母线电压的日常运行区间，另一方面实现了逆调压；并且通过在线自动滚动计算，也可以自动适应电网运行方式的变化。本发明应用后能够实现系统的逆调压，在低谷时段压低电压值上限,在高峰时段提高电压值上限,并将计算的控制限值应用于自动电压控制系统，实现逆调压，以提高电力系统的运行稳定性。

Description

一种基于母线电压历史数据分析的逆调压控制方法

技术领域

本发明属于电力系统自动电压控制技术领域，特别涉及一种基于母线电压历史数据分析的逆调压控制方法的方法。

背景技术

自动电压控制(以下简称AVC，Automatic Voltage Control)系统是实现电网安全(提高电压稳定裕度)、经济(降低网络损耗)、优质(提高电压合格率)运行的重要手段。AVC系统架构在电网能量管理系统(EMS)之上，能够利用电网实时运行数据，从电网全局优化的角度科学决策出最佳的无功电压调整方案，自动下发给电厂、以及下级电网调度机构执行。AVC系统以电压安全和优质为约束，以系统运行经济性为目标，连续闭环地进行电压的实时优化控制，实现了无功电压协调控制方案的在线生成、实时下发、闭环自动控制等一整套分析、决策、控制，以及再分析，再决策、再控制的无功电压实时追踪控制问题，能够有效地克服传统的电网无功电压控制手段存在的不足，提高电网安全稳定经济运行的水平。孙宏斌、张伯明、郭庆来在《基于软分区的全局电压优化控制系统设计》(电力系统自动化,2003年，第27卷第8期，16-20页)中说明了大电网自动电压控制的体系结构。

目前省级电网调度中心的AVC系统采用协调二级电压控制方法(CSVC)实现对某个包含发电厂的区域电网的电压控制，二级控制为分区解耦的控制策略计算，其将电网自动划分为解耦的各个分区，每个分区内选择若干中枢母线，对每个分区分别计算分区内包括电厂在内的各种无功资源的控制策略，以追随该分区中枢母线的优化目标。郭庆来，孙宏斌，张伯明在《协调二级电压控制的研究》(电力系统自动化，2005年12月，V29N23，pp.19-24)中提出了一种协调二级电压控制(CSVC)模型，该模型在优先考虑中枢母线电压偏差最小的前提下，利用多余的控制自由度保证本区域发电机运行在无功裕度更大、出力更均衡的状态。该模型所涉及到的各个变量的具体物理含义可以从图1所示的二次规划计算电厂模型中直观看出，其中，Q_g1、Q_g2、Q_g3、Q_g4分别表示发电机1、2、3、4的当前无功出力，V_g1、V_g2、V_g3、V_g4分别表示发电机1、2、3、4机端母线当前电压，V_p1、V_p2分别表示中枢母线1和2的当前电压，V_H1、V_H2分别表示发电机高压侧母线1和2的当前电压，C_g为中枢母线无功电压灵敏度矩阵和C_vg为高压侧母线无功电压灵敏度矩阵，满足：

ΔV_p＝C_gΔQ_g (1.1)

ΔV_H＝C_vgΔQ_g (1.2)

其中，ΔQ_g为发电机无功出力的调节量，ΔV_p为机组对中枢母线的电压调整量，ΔV_H为机组对电厂高压侧母线的电压调整量。区域中的中枢母线一般为预先人工选择指定。该模型采用了二次规划模型来计算电厂的控制，二次规划模型的目标函数如下：

其中，ΔQ_g为二次规划的优化变量，表示发电机无功出力的调节量；

和

分别表示无功下限和无功上限；

表示区域中的中枢母线设定电压；W_p和W_q为权重系数，两者一般均为0.1到1之间，α为增益系数，一般为0.1到1之间；式中目标函数第一部分体现了调整发电机无功使得中枢母线达到三级控制给出的目标值。同时，为了实现增大发电机无功裕度并使之出力更加均衡的目的，在目标函数的第二部分引入了无功裕度向量Θ_g，其第i个分量为:

其中，

为发电机i的当前无功出力，

为发电机i的无功出力的调节量，

为发电机i的无功最大值，

为发电机i的无功最小值，i为发电机的序号，||Θ_g||²为无功裕度向量集，||Θ_g||²在二次规划目标函数中，可以保证一方面增加控制发电机的无功裕度，另一方面促使各台控制发电机向无功出力更加均衡的方向发展。

省调AVC的二级电压控制模块在满足安全约束条件的情况下来求解极小化问题，这些约束包括：

其中

和

分别表示中枢母线电压下限和中枢母线电压上限；

和

分别表示发电机无功下限和无功上限；

和

分别表示发电机高压侧母线的电压下限、电压上限和允许的单步最大调整量。

上述内容涉及到无功电压灵敏度矩阵C_g和C_vg的计算。孙宏斌，张伯明，相年德在《准稳态的灵敏度分析方法》(中国电机工程学报，1999年4月V19 N4，pp.9-13)中提出了准稳态灵敏度方法，与常规的静态的灵敏度分析方法不同，准稳态灵敏度方法考虑了电力系统准稳态的物理响应，计及系统控制前后新旧稳态间的总变化，有效提高了灵敏度分析的精度。该方法基于电力系统的PQ解耦模型，当发电机安装有自动电压调节器(AVR)时，可认为该发电机节点为PV节点；而当发电机装有自动无功功率调节(AQR)或自动功率因数调节(APFR)时，可认为该发电机节点与普通负荷节点相同均为PQ节点。此外，将负荷电压静特性考虑成节点电压的一次或二次曲线。这样所建立的潮流模型就自然地将这些准稳态的物理响应加以考虑，从而基于潮流模型计算出的灵敏度即为准稳态的灵敏度。省调AVC中，C_g和C_vg均采用准稳态的灵敏度。

电力系统中各级电压偏差值应不超过允许范围。行业标准DL/T 1773-2017《电力系统电压和无功电力技术导则》中规定:对于500(330)kV母线，正常运行方式时、最高运行电压不得超过系统额定电压的+10％，最低运行电压不应影响电力系统同步稳定、发电厂用电的正常使用及下一级电压的调节。对于220kV母线.正常运行方式时.电压允许偏差为系统额定电压的0～+10％。事故运行方式时为系统额定电压的-5％～+10％。对于110-35kv母线，正常运行方式时为相应系统额定电压的-3％～+7％事故后为额定电压的±10％。电力系统调压是指为使电力系统中各电压中枢点运行电压保持在规定允许范围之内所采取的技术措施，主要调压方式分为逆调压，恒调压和顺调压三种。

(1)逆调压。逆调压是指在最大负荷时段，提高系统电压中枢点电压至105％倍标准电压以补偿线路上增加的电压损失，最小负荷时段降低中枢点电压至标准电压以防止受端电压过高的电压调整方式。逆调压的作用是为了使电压偏差符合用电设备端电压的要求，逆调压的范围宜为额定电压的0～+5％。

(2)恒调压。常调压是指无论负荷如何变动，系统电压中枢点电压基本保持不变的电压调整方式，一般保持中枢点电压在102％—105％倍额定电压。通常适用于负荷变动小或者线路上电压损耗小的情况。

(3)顺调压。顺调压是指在最大负荷时适当降低中枢点电压，但不低于102.5％倍额定电压，最小负荷时适当加大中枢点电压的电压调整方式，但不高于107.5％倍额定电压。通常适用于出线线路不太长，负荷变化不大的情况。

在输电网自动电压控制过程中，通常在系统负荷的低谷时段对变电站母线电压进行逆调压控制，在保证电网电压合格同时，适度降低母线电压运行水平，一方面提高设备绝缘器件的运行寿命，另一方面也可以预防故障引起的过电压。在AVC中实现逆调压的方法是在低谷时段降低的母线电压上限值，在平峰或高峰时段提高母线的电压上限值，从而实现低谷时段电压低于其他时段的电压值。目前常规逆调压控制中的低谷、平峰、高峰时段母线电压的上限值，是人工根据经验提前设定的，是固定不变的。由于没有考虑电网中无功资源配置的实际情况，以及电网运行方式的变化，人工设定的各时段上限值存在不合理的情况，影响自动电压控制的效果。

发明内容

本发明的目的是提出一种基于母线电压历史数据分析的逆调压控制方法，对已有的自动电压控制方法进行改进，以电压历史数据作为统计分析的依据，根据电压数值、计划限值、以及对应的发电机\容抗器等控制设备的调节能力，在低谷时段压低电压值上限,在高峰时段提高电压值上限,并将计算的控制限值应用于自动电压控制系统，实现逆调压，以提高电力系统的稳定性。

本发明提出的基于母线电压历史数据分析的逆调压控制方法，包括以下步骤：

(1)从电力系统的自动电压控制系统中，在每日0时到来时，读取电力系统中的每条母线的名称，形成母线名称集合B：B_i＝|B₁ B₂ B₃ B₄ ... B_I|，i＝1,...,I，B为母线名称集合，下标i为母线序号，B_i为第i个母线；读取电力系统中发电机组的名称，形成发电机组名称的集合G：G_j＝|G₁ G₂ G₃ G₄ ... G_J|，j＝1,...,J，G为发电机组名称集合，下标j为发电机组序号，G_j为第j个发电机组；读取电力系统中容抗器的名称，形成容抗器的名称集合C，C_k＝|C₁ C₂ C₃ C₄ ... C_K|，k＝1,...,K，C为容抗器名称集合，下标k为容抗器序号，C_k为第k个容抗器；

(2)从电力系统中读取母线电压设定计划值以及母线、发电机组和容抗器的历史值，设定一天内的采样点为n＝1,...,N、历史天数m＝1,...,M，历史天数m为从当前0点前的M天中的一天，发电机组对母线的灵敏度和变电站低压母线对高压母线的灵敏度，容抗器容量，进行数据整合形成数据集合，包括以下步骤：

(2-1)根据母线集合B，读取母线电压设定计划值，形成如下所有母线电压设定计划值的集合：

V^max、V^min分别为母线电压设定计划值的上限集合和下限集合，

分别为第i个母线在一天内的第n个采样点的电压设定计划值的上限数据值和电压设定计划值的下限数据值。

(2-2)根据母线名称集合B，读取第i个母线的前M天每天采样点N个点的电压历史值，形成数据集合：

V_i ^m-n为第i个母线的第m天第n个电压历史采样点，遍历母线序号i，形成所有母线的电压历史数据集合V：V＝|V₁ V₂ V₃ V₄ ... V_I|；

(2-3)根据发电机组集合G，读取第j个发电机组的前M天每天采样点N个点的无功历史值，形成数据集合：

为第j个发电机组的第m天第n个无功历史采样点，遍历发电机组序号j，形成所有发电机组的无功历史数据集合：

Q＝|Q₁ Q₂ Q₃ Q₄ ... Q_J|；

(2-4)根据容抗器集合C，读取第k个容抗器的前M天每天采样点N个点的运行状态历史值，形成数据集合：

为第k个发电机组的第m天第n个运行状态历史采样点，遍历容抗器序号k，形成所有容抗器的运行状态历史数据集合：D＝|D₁ D₂ D₃ D₄ ... D_K|；

(2-5)根据发电机组集合G和母线集合B，读取发电机组对母线的灵敏度信息，形成集合如下：

其中，S_GV为发电机组对母线的灵敏度矩阵，S_ji为第j个发电机组对第i个母线的灵敏度；

(2-6)读取变电站低压母线对高压母线的灵敏度，形成变电站低压母线对高压母线的灵敏度集合S_CV：

其中，S_CV为变电站低压母线对高压母线的灵敏度矩阵，S_ii为第w个低压母线对第r个高压母线的灵敏度；w为变电站低压母线个数，r为变电站高压母线个数；

(2-7)根据容抗器集合C，读取第k个容抗器的容量

形成所有容抗器的容量集合：

(3)根据步骤(2)得到的数据集合，计算发电厂控制母线、变电站控制母线的逆调压电压上限，步骤如下：

(3-1)读取母线集合的母线B_i，对母线B_i进行判断，若B_i为发电厂母线，则执行步骤(3-2),若B_i为变电站母线，则执行步骤(3-3)；

(3-2)根据母线集合的母线B_i，从发电机组集合G中读取与母线B_i相连接的发电机组，形成集合G_G：G_f＝|G₁ G₂ G₃ G₄ ... G_F|、f＝1,...,F，其中，F为与母线B_i相连接的发电机组的个数，G_f为第f个发电机组，并根据相连接的发电机组集合G_G分别从发电机组的无功历史数据集合Q中获取该发电机组的无功历史值Q_f，同时从母线B_i电压设定计划值集合V^max、V^min中分别读取母线B_i的设定计划值上限

设定计划值下限

根据上述读取到的集合数据，计算母线B_i的逆调压电压上限V_i ^smax，包括以下步骤：

(3-2-1)设定母线B_i电压向下调节能力V_id＝0、电压向上调节能力V_iu＝0，计算母线B_i的电压向下调节能力V_id、电压向上调节能力V_iu，步骤如下：

(3-2-1-1)设定发电机组G_f可减少无功Q_fs＝0、发电机组G_f可减少无功的采样计数点Z₁＝0，记发电机组G_f的最小无功出力为Q_min，从发电机组的无功历史值Q_f中读取

其中，上标m-n为第m天的第n个采样点、

为第m天的第n个点的无功值，将发电机组G_f的第m天第n个点的历史值

与发电机组最小无功出力Q_min进行比较，若

则使发电机组G_f的可减少无功

使采样计数点Z₁＝Z₁+1；遍历历史天数m，将所有得到的发电机组G_f的可减少无功Q_fs进行累加，并将所有发电机组G_f可减少无功的采样计数点Z₁的累加；根据累加值，计算发电机组G_f可减少无功值的平均值Q_avg＝Q_fs/Z₁；读取发电机组G_f对母线B_i的灵敏度S_GV，根据灵敏度S_GV和发电机组G_f可减少无功平均值Q_avg计算得到母线B_i电压向下调节能力V_adjd＝Q_avg×S_GV×V_x,V_x为安全系数，根据电压向下调节能力V_adjd，计算得到母线B_i的电压向下调节能力累加：V_id＝V_id+V_adjd；

(3-2-1-2)设定发电机组G_f的可增加无功Q_fu＝0、发电机组G_f可增加无功的采样计数点Z₂＝0，将发电机组G_f的最大无功出力记为Q_max，从发电机组G_f的无功历史值Q_f中读取

将发电机组的第m天第n个点的历史值

与发电机组最大无功出力进行比较，若

则使

使Z₂＝Z₂+1；遍历历史天数m，将所有得到的发电机组的可增加无功Q_fu进行累加，并将发电机组可增加无功的采样计数点Z₂进行累加，根据累加值计算得到发电机组G_f可增加无功平均值Q_pju＝Q_fu/Z₂；读取发电机组G_f对母线B_i的灵敏度S_GV，根据灵敏度S_GV和发电机组G_f可增加无功平均值Q_avg，计算得到母线B_i电压向上调节能力V_adju＝Q_pju×S_GV×V_x,V_x为安全系数；据电压向上调节能力V_adju，计算得到母线B_i的电压向上调节能力累加：V_iu＝V_iu+V_adju；

(3-2-1-3)对与母线B_i相连接的发电机组集合G_G，遍历发电机组个数f的取值，重复步骤(3-2-1-1)～步骤(3-2-1-2)，得到所有与母线B_i相连接的发电机组对母线B_i总的电压向下调节能力V_id和电压向上调节能力V_iu；

(3-2-2)设定母线B_i的电压值V_s＝0，母线B_i电压计数点Z₃＝0，从步骤(2-2)的母线B_i的电压历史值V_i中读取V_i ^m-n，将母线B_i第m天第n个点的历史值V_i ^m-n与母线B_i的电压设定计划值

比较，若

并且

则使母线B_i电压V_s＝V_s+V_i ^m-n，并使计数点Z₃＝Z₃+1；遍历历史天数m，对得到的母线B_i的所有电压值V_s进行累加，对所有母线B_i的电压计数点Z₃进行累加，根据得到的累加数据计算得到母线B_i的电压平均值V_ipj＝V_s/Z₃；

(3-2-3)对采样点n进行判断，若采样点n为电力系统负荷低谷时刻，则根据(3-2-2)计算得到的母线B_i的电压平均值V_ipj和(3-2-1-3)计算得到的母线B_i总的向下调节能力V_id，计算得到母线B_i的逆调压电压上限

并读取母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

若采样点n的时刻为电力系统负荷平峰或高峰时刻，则根据步骤(3-2-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-2-1-3)计算得到的母线B_i总的电压向上调节能力V_iu，计算母线B_i的逆调压电压上限

并读取母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

(3-3)根据母线集合的母线B_i，从容抗器集合C中读取与母线B_i相连接的容抗器C_C＝|G₁ G₂ G₃ G₄ ... G_O|，其中o＝1,...,O，O为相连接的容抗器的个数，根据相连接的容抗器集合C_C，从容抗器G_o的无功历史数据集合D中获取容抗器的历史值D_o，从母线B_i电压设定计划值集合V^max、V^min中读取母线B_i的设定计划值上限

设定计划值下限

根据上述读取到的集合数据，计算母线B_i的逆调压电压上限V_i ^smax，包括以下步骤：

(3-3-1)设定母线B_i的电压向下调节能力V_id＝0、电压向上调节能力V_iu＝0，计算母线B_i的电压向下调节能力V_id、电压向上调节能力V_iu，步骤如下：

(3-3-1-1)设定母线B_i可增加无功容量Q_vh＝0、母线B_i可增加无功容量的采样计数点Z₄＝0、母线B_i可减无功容量Q_vl＝0、母线B_i可减无功容量的采样计数点Z₅＝0,从步骤(2-4)容抗器G_o的运行状态历史值D_o中读取

并从步骤(2-7)容抗器的容量集合C^rmx中读取

对G_o的容量

进行判断：若

则判定G_o为电容器，并进而根据电容器G_o的运行状态的历史值

进行判断，若C_o为已投入运行的电容器，则使

Q_vl＝Q_vl+1，若C_o为热备用的电容器，则使

使Z₄＝Z₄+1；若

为负，则判定G_o为电抗器，并进而根据电抗器G_o的运行状态的历史值

进行判断，若C_o为已投入运行的电抗器，则使

Z₄＝Z₄+1，若C_o为热备用的电抗器，则使

Q_vl＝Q_vl+1；遍历历史天数m,得到母线B_i所有可增加无功容量Q_vh的累加值、母线B_i所有可增加无功容量的采样计数点Z₄的累加值、母线B_i所有可减无功容量Q_vl的累加值以及母线B_i所有可减无功容量的采样计数点Z₅的累加值；

(3-3-1-2)根据步骤(3-3-1-1)的母线B_i所有可增加无功容量Q_vh的累加值、B_i所有可增加无功容量的采样计数点Z₄的累加值，计算母线B_i的可增加无功容量平均值Q_pjh＝Q_vh/Z₄，根据母线B_i所有可减无功容量Q_vl的累加、母线B_i所有可减无功容量的采样计数点Z₅的累加值，计算得到母线B_i的可减少无功容量平均值Q_pjl＝Q_vl/Z₅；从步骤(2-5)读取变电站低压母线对高压母线灵敏度S_VV，根据变电站低压母线对高压母线灵敏度S_VV和可减少无功容量平均值Q_pjl、可增加无功容量平均值Q_pjh，分别计算得到母线B_i的电压向下调节能力V_tjd＝Q_pjl×S_VV×V_x和电压向上调节能力V_tju＝Q_pjh×S_VV×V_x，V_x为安全系数；根据计算得到的母线B_i的电压向下调节能力V_tjd、电压向上调节能力V_tju，分别计算得到母线B_i的累加电压向下调节能力V_id＝V_id+V_tjd和电压向上调节能力V_iu＝V_iu+V_tju；

(3-3-1-3)对与母线B_i相连接的容抗器集合C_C，遍历容抗器个数o的取值，重复步骤(3-2-1-1)～步骤(3-2-1-2)，得到与母线B_i的所有连接的容抗器对母线B_i的总向下调节能力V_id、总电压向上调节能力V_iu；

(3-3-2)设定母线B_i电压V_s＝0，母线B_i电压采样点数Z₃＝0，从母线B_i的电压历史值V_i中读取V_i ^m-n，将母线B_i第m天第n个点的历史值V_i ^m-n与母线B_i的电压设定计划值

比较，若

且

则使母线B_i的电压V_s＝V_s+V_i ^m-n、使Z₃＝Z₃+1；遍历历史天数m，得到母线B_i的所有电压V_s的累加值、母线B_i电压采样点数Z₃的累加值，根据累加值计算得到母线B_i的电压平均值V_ipj＝V_s/Z₃；

(3-3-3)对采样点n进行判断，若采样点n为电力系统负荷低谷时刻，则根据步骤(3-3-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-3-1-3)计算得到的母线B_i总的向下调节能力V_id，计算得到母线B_i的逆调压电压上限

并进而根据母线B_i的电压设定计划值

对逆调压电压上限进行如下修正：若

则使

若

则使

若采样点n的时刻为电力系统负荷平峰或高峰时刻，则根据步骤(3-3-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-3-1-3)计算得到的母线B_i总的电压向上调节能力V_iu，计算得到母线B_i的逆调压电压上限

并进而根据母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

(3-4)遍历步骤(1)的母线集合B中的所有母线，遍历采样点n，重复步骤(3-1)～步骤(3-3)，计算并形成母线B_i的逆调压电压上限集合：

为第i个母线第n个采样点的逆调压电压上限值；

(4)设定自动电压控制周期为T_c，在每个自动电压控制周期T_c到来时，记到来时刻为t₀，计算自动电压控制发电厂无功控制指令，包括以下步骤：

(4-1)建立一个自动电压控制优化模型，该优化模型的目标函数的表达式如下：

其中，V_p表示电力系统区域中枢母线的当前电压，

表示电力系统区域中枢母线设定电压，C_g为中枢母线无功电压灵敏度矩阵，ΔQ_g表示发电机无功出力的调节量，W_p和W_q为发电机组的权重系数，取值范围为(0.1-10)，α为增益系数，取值范围为0.1～10，Θ_g为无功裕度向量，Θ_g的第i个分量为：

其中，Q_gi为发电机i的当前无功出力，ΔQ_gi为发电机i的无功出力的调节量，

为发电机i的无功最大值，

为发电机i的无功最小值，i为发电机的序号，||Θ_g||²为无功裕度向量集，以上数据量分别从电力系统中自动电压控制系统中获取；

(4-2)确定自动电压控制优化模型的约束条件，表达式如下：

其中，C₀、C₁、C₂、C₃分别为限值约束方程，C_vg为高压侧母线无功电压灵敏度矩阵，从电力系统中读取，V_H表示发电机高压侧母线的当前电压，

和

分别表示中枢母线电压下限和中枢母线电压上限；

和

分别表示发电机无功下限和无功上限；

和

分别表示发电机高压侧母线的电压下限、电压上限和允许的单步电压最大调整量，以上数据量分别从电力系统中自动电压控制系统中获取；

(4-3)根据(2)定义的一天内的采样点n＝1,...,N，计算采样点时间间隔t_t＝24×60÷N，通过当前到来时刻t₀及采样点时间间隔t_t，计算得到当前采样点n₀＝t₀÷t_t；从步骤(3-4)逆调压电压上限集合V^smax中，获取与当前到来t₀时刻相对应的逆调压电压上限

并以该

分别取代步骤(4-2)中限值约束方程C₁中的

和限值约束方程C₂中的

得到更新后的自动电压控制优化模型；

(4-4)采用二次规划计算方法，求解步骤(4-3)的优化模型，得到发电机无功出力的调节量ΔQ_g，并将该发电机无功出力的调节量ΔQ_g下发到发电厂的各自动电压控制电厂子站，实现基于母线电压历史数据分析的逆调压控制。

本发明提出的基于母线电压历史数据分析的逆调压控制方法，其优点是：

本发明的基于母线电压历史数据分析的逆调压控制方法，在电压控制中，根据变电站在最近时期内的母线电压运行的实际范围，结合电厂机组、变电站容抗器等无功设备的运行情况，自动计算母线电压上限，一方面满足了母线电压的日常运行区间，另一方面实现了逆调压；并且通过在线自动滚动计算，也可以自动适应电网运行方式的变化。本发明应用后能够实现系统的逆调压，提高了电力系统的运行稳定性。

附图说明

图1是本发明方法涉及的二次规划计算电厂模型示意图。

图2是本发明方法的流程框图。

具体实施方式

本发明提出的基于母线电压历史数据分析的逆调压控制方法，其流程框图如图2所示，包括以下步骤：

(1)从电力系统的自动电压控制系统中，在每日0时到来时，读取电力系统中的每条母线的名称，形成母线名称集合B：B_i＝|B₁ B₂ B₃ B₄ ... B_I|，i＝1,...,I，B为母线名称集合，下标i为母线序号，B_i为第i个母线；读取电力系统中发电机组的名称，形成发电机组名称的集合G：G_j＝|G₁ G₂ G₃ G₄ ... G_J|，j＝1,...,J，G为发电机组名称集合，下标j为发电机组序号，G_j为第j个发电机组；读取电力系统中容抗器的名称，形成容抗器的名称集合C，C_k＝|C₁ C₂ C₃ C₄ ... C_K|，k＝1,...,K，C为容抗器名称集合，下标k为容抗器序号，C_k为第k个容抗器；

(2)从电力系统中读取母线电压设定计划值以及母线、发电机组和容抗器的历史值，设定一天内的采样点为n＝1,...,N、历史天数m＝1,...,M，一般m取15分钟为一个采样点，一天共96点，历史天数m为从当前0点前的M天中的一天，一般M取30，发电机组对母线的灵敏度和变电站低压母线对高压母线的灵敏度，容抗器容量，进行数据整合形成数据集合，包括以下步骤：

(2-1)根据母线集合B，读取母线电压设定计划值，形成如下所有母线电压设定计划值的集合：

V^max、V^min分别为母线电压设定计划值的上限集合和下限集合，

分别为第i个母线在一天内的第n个采样点的电压设定计划值的上限数据值和电压设定计划值的下限数据值。

(2-2)根据母线名称集合B，读取第i个母线的前M天每天采样点N个点的电压历史值，形成数据集合：

V_i ^m-n为第i个母线的第m天第n个电压历史采样点，遍历母线序号i，形成所有母线的电压历史数据集合V：V＝|V₁ V₂ V₃ V₄ ... V_I|；

(2-3)根据发电机组集合G，读取第j个发电机组的前M天每天采样点N个点的无功历史值，形成数据集合：

为第j个发电机组的第m天第n个无功历史采样点，遍历发电机组序号j，形成所有发电机组的无功历史数据集合：

Q＝|Q₁ Q₂ Q₃ Q₄ ... Q_J|；

(2-4)根据容抗器集合C，读取第k个容抗器的前M天每天采样点N个点的运行状态历史值，形成数据集合：

为第k个发电机组的第m天第n个运行状态历史采样点，遍历容抗器序号k，形成所有容抗器的运行状态历史数据集合：D＝|D₁ D₂ D₃ D₄ ... D_K|；

(2-5)根据发电机组集合G和母线集合B，读取发电机组对母线的灵敏度信息，形成集合如下：

其中，S_GV为发电机组对母线的灵敏度矩阵，S_ji为第j个发电机组对第i个母线的灵敏度；

(2-6)读取变电站低压母线对高压母线的灵敏度，形成变电站低压母线对高压母线的灵敏度集合S_CV：

其中，S_CV为变电站低压母线对高压母线的灵敏度矩阵，S_ii为第w个低压母线对第r个高压母线的灵敏度；w为变电站低压母线个数，r为变电站高压母线个数；

(2-7)根据容抗器集合C，读取第k个容抗器的容量

形成所有容抗器的容量集合：

(3)根据步骤(2)得到的数据集合，计算发电厂控制母线、变电站控制母线的逆调压电压上限，步骤如下：

(3-1)读取母线集合的母线B_i，对母线B_i进行判断，若B_i为发电厂母线，则执行步骤(3-2),若B_i为变电站母线，则执行步骤(3-3)；

(3-2)根据母线集合的母线B_i，从发电机组集合G中读取与母线B_i相连接的发电机组，形成集合G_G：G_f＝|G₁ G₂ G₃ G₄ ... G_F|、f＝1,...,F，其中，F为与母线B_i相连接的发电机组的个数，G_f为第f个发电机组，并根据相连接的发电机组集合G_G分别从发电机组的无功历史数据集合Q中获取该发电机组的无功历史值Q_f，同时从母线B_i电压设定计划值集合V^max、V^min中分别读取母线B_i的设定计划值上限

设定计划值下限

根据上述读取到的集合数据，计算母线B_i的逆调压电压上限V_i ^smax，包括以下步骤：

(3-2-1)设定母线B_i电压向下调节能力V_id＝0、电压向上调节能力V_iu＝0，计算母线B_i的电压向下调节能力V_id、电压向上调节能力V_iu，步骤如下：

(3-2-1-1)设定发电机组G_f可减少无功Q_fs＝0、发电机组G_f可减少无功的采样计数点Z₁＝0，记发电机组G_f的最小无功出力为Q_min，本发明的一个实施例中，G_f的最小无功出力为Q_min＝20，从发电机组的无功历史值Q_f中读取

其中，上标m-n为第m天的第n个采样点、

为第m天的第n个点的无功值，将发电机组G_f的第m天第n个点的历史值

与发电机组最小无功出力Q_min进行比较，若

则使发电机组G_f的可减少无功

使采样计数点Z₁＝Z₁+1；遍历历史天数m，将所有得到的发电机组G_f的可减少无功Q_fs进行累加，并将所有发电机组G_f可减少无功的采样计数点Z₁的累加；根据累加值，计算发电机组G_f可减少无功值的平均值Q_avg＝Q_fs/Z₁；读取发电机组G_f对母线B_i的灵敏度S_GV，根据灵敏度S_GV和发电机组G_f可减少无功平均值Q_avg计算得到母线B_i电压向下调节能力V_adjd＝Q_avg×S_GV×V_x,V_x为安全系数，本发明的一个实施例中，V_x的取值为0.5，根据电压向下调节能力V_adjd，计算得到母线B_i的电压向下调节能力累加：V_id＝V_id+V_adjd；

(3-2-1-2)设定发电机组G_f的可增加无功Q_fu＝0、发电机组G_f可增加无功的采样计数点Z₂＝0，将发电机组G_f的最大无功出力记为Q_max，本发明的一个实施例中Q_max的取值为100，从发电机组G_f的无功历史值Q_f中读取

将发电机组的第m天第n个点的历史值

与发电机组最大无功出力进行比较，若

则使

使Z₂＝Z₂+1；遍历历史天数m，将所有得到的发电机组的可增加无功Q_fu进行累加，并将发电机组可增加无功的采样计数点Z₂进行累加，根据累加值计算得到发电机组G_f可增加无功平均值Q_pju＝Q_fu/Z₂；读取发电机组G_f对母线B_i的灵敏度S_GV，根据灵敏度S_GV和发电机组G_f可增加无功平均值Q_avg，计算得到母线B_i电压向上调节能力V_adju＝Q_pju×S_GV×V_x,V_x为安全系数，在本发明的一个实施例中V_x的取值为0.5；根据电压向上调节能力V_adju，计算得到母线B_i的电压向上调节能力累加：V_iu＝V_iu+V_adju；

(3-2-1-3)对与母线B_i相连接的发电机组集合G_G，遍历发电机组个数f的取值，重复步骤(3-2-1-1)～步骤(3-2-1-2)，得到所有与母线B_i相连接的发电机组对母线B_i总的电压向下调节能力V_id和电压向上调节能力V_iu；

(3-2-2)设定母线B_i的电压值V_s＝0，母线B_i电压计数点Z₃＝0，从步骤(2-2)的母线B_i的电压历史值V_i中读取V_i ^m-n，将母线B_i第m天第n个点的历史值V_i ^m-n与母线B_i的电压设定计划值

比较，若

并且

则使母线B_i电压V_s＝V_s+V_i ^m-n，并使计数点Z₃＝Z₃+1；遍历历史天数m，对得到的母线B_i的所有电压值V_s进行累加，对所有母线B_i的电压计数点Z₃进行累加，根据得到的累加数据计算得到母线B_i的电压平均值V_ipj＝V_s/Z₃；

(3-2-3)对采样点n进行判断，若采样点n为电力系统负荷低谷时刻，则根据(3-2-2)计算得到的母线B_i的电压平均值V_ipj和(3-2-1-3)计算得到的母线B_i总的向下调节能力V_id，计算得到母线B_i的逆调压电压上限

并读取母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

若采样点n的时刻为电力系统负荷平峰或高峰时刻，则根据步骤(3-2-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-2-1-3)计算得到的母线B_i总的电压向上调节能力V_iu，计算母线B_i的逆调压电压上限

并读取母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

(3-3)根据母线集合的母线B_i，从容抗器集合C中读取与母线B_i相连接的容抗器C_C＝|G₁ G₂ G₃ G₄ ... G_O|，其中o＝1,...,O，O为相连接的容抗器的个数，根据相连接的容抗器集合C_C，从容抗器G_o的无功历史数据集合D中获取容抗器的历史值D_o，从母线B_i电压设定计划值集合V^max、V^min中读取母线B_i的设定计划值上限

设定计划值下限

根据上述读取到的集合数据，计算母线B_i的逆调压电压上限V_i ^smax，包括以下步骤：

(3-3-1)设定母线B_i的电压向下调节能力V_id＝0、电压向上调节能力V_iu＝0，计算母线B_i的电压向下调节能力V_id、电压向上调节能力V_iu，步骤如下：

(3-3-1-1)设定母线B_i可增加无功容量Q_vh＝0、母线B_i可增加无功容量的采样计数点Z₄＝0、母线B_i可减无功容量Q_vl＝0、母线B_i可减无功容量的采样计数点Z₅＝0,从步骤(2-4)容抗器G_o的运行状态历史值D_o中读取

并从步骤(2-7)容抗器的容量集合C^rmx中读取

对G_o的容量

进行判断：若

则判定G_o为电容器，并进而根据电容器G_o的运行状态的历史值

进行判断，若C_o为已投入运行的电容器，则使

Q_vl＝Q_vl+1，若C_o为热备用的电容器，则使

使Z₄＝Z₄+1；若

为负，则判定G_o为电抗器，并进而根据电抗器G_o的运行状态的历史值

进行判断，若C_o为已投入运行的电抗器，则使

Z₄＝Z₄+1，若C_o为热备用的电抗器，则使

遍历历史天数m,得到母线B_i所有可增加无功容量Q_vh的累加值、母线B_i所有可增加无功容量的采样计数点Z₄的累加值、母线B_i所有可减无功容量Q_vl的累加值以及母线B_i所有可减无功容量的采样计数点Z₅的累加值；

(3-3-1-2)根据步骤(3-3-1-1)的母线B_i所有可增加无功容量Q_vh的累加值、B_i所有可增加无功容量的采样计数点Z₄的累加值，计算母线B_i的可增加无功容量平均值Q_pjh＝Q_vh/Z₄，根据母线B_i所有可减无功容量Q_vl的累加、母线B_i所有可减无功容量的采样计数点Z₅的累加值，计算得到母线B_i的可减少无功容量平均值Q_pjl＝Q_vl/Z₅；从步骤(2-5)读取变电站低压母线对高压母线灵敏度S_VV，根据变电站低压母线对高压母线灵敏度S_VV和可减少无功容量平均值Q_pjl、可增加无功容量平均值Q_pjh，分别计算得到母线B_i的电压向下调节能力V_tjd＝Q_pjl×S_VV×V_x和电压向上调节能力V_tju＝Q_pjh×S_VV×V_x，V_x为安全系数，本发明的一个实施例中为0.5；根据计算得到的母线B_i的电压向下调节能力V_tjd、电压向上调节能力V_tju，分别计算得到母线B_i的累加电压向下调节能力V_id＝V_id+V_tjd和电压向上调节能力V_iu＝V_iu+V_tju；

(3-3-1-3)对与母线B_i相连接的容抗器集合C_C，遍历容抗器个数o的取值，重复步骤(3-2-1-1)～步骤(3-2-1-2)，得到与母线B_i的所有连接的容抗器对母线B_i的总向下调节能力V_id、总电压向上调节能力V_iu；

(3-3-2)设定母线B_i电压V_s＝0，母线B_i电压采样点数Z₃＝0，从母线B_i的电压历史值V_i中读取V_i ^m-n，将母线B_i第m天第n个点的历史值V_i ^m-n与母线B_i的电压设定计划值

比较，若

且

则使母线B_i的电压V_s＝V_s+V_i ^m-n、使Z₃＝Z₃+1；遍历历史天数m，得到母线B_i的所有电压V_s的累加值、母线B_i电压采样点数Z₃的累加值，根据累加值计算得到母线B_i的电压平均值V_ipj＝V_s/Z₃；

(3-3-3)对采样点n进行判断，若采样点n为电力系统负荷低谷时刻，则根据步骤(3-3-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-3-1-3)计算得到的母线B_i总的向下调节能力V_id，计算得到母线B_i的逆调压电压上限

并进而根据母线B_i的电压设定计划值

对逆调压电压上限进行如下修正：若

则使

若

则使

若采样点n的时刻为电力系统负荷平峰或高峰时刻，则根据步骤(3-3-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-3-1-3)计算得到的母线B_i总的电压向上调节能力V_iu，计算得到母线B_i的逆调压电压上限

并进而根据母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

(3-4)遍历步骤(1)的母线集合B中的所有母线，遍历采样点n，重复步骤(3-1)～步骤(3-3)，计算并形成母线B_i的逆调压电压上限集合：

为第i个母线第n个采样点的逆调压电压上限值；

(4)设定自动电压控制周期为T_c，在每个自动电压控制周期T_c到来时，记到来时刻为t₀，计算自动电压控制发电厂无功控制指令，包括以下步骤：

(4-1)建立一个自动电压控制优化模型，该优化模型的目标函数的表达式如下：

其中，V_p表示电力系统区域中枢母线的当前电压，

表示电力系统区域中枢母线设定电压，C_g为中枢母线无功电压灵敏度矩阵，ΔQ_g表示发电机无功出力的调节量，W_p和W_q为发电机组的权重系数，取值范围为(0.1-10)，α为增益系数，取值范围为0.1～10，Θ_g为无功裕度向量，Θ_g的第i个分量为：

其中，

为发电机i的当前无功出力，

为发电机i的无功出力的调节量，

为发电机i的无功最大值，

为发电机i的无功最小值，i为发电机的序号，||Θ_g||²为无功裕度向量集，以上数据量分别从电力系统中自动电压控制系统中获取；

(4-2)确定自动电压控制优化模型的约束条件，表达式如下：

其中，C₀、C₁、C₂、C₃分别为限值约束方程，C_vg为高压侧母线无功电压灵敏度矩阵，从电力系统中读取，V_H表示发电机高压侧母线的当前电压，

和

分别表示中枢母线电压下限和中枢母线电压上限；

和

分别表示发电机无功下限和无功上限；

和

分别表示发电机高压侧母线的电压下限、电压上限和允许的单步电压最大调整量，以上数据量分别从电力系统中自动电压控制系统中获取；

(4-3)根据(2)定义的一天内的采样点n＝1,...,N，计算采样点时间间隔t_t＝24×60÷N，通过当前到来时刻t₀及采样点时间间隔t_t，计算得到当前采样点n₀＝t₀÷t_t；从步骤(3-4)逆调压电压上限集合V^smax中，获取与当前到来t₀时刻相对应的逆调压电压上限

并以该

分别取代步骤(4-2)中限值约束方程C₁中的

和限值约束方程C₂中的

得到更新后的自动电压控制优化模型；

(4-4)采用二次规划计算方法，求解步骤(4-3)的优化模型，得到发电机无功出力的调节量ΔQ_g，并将该发电机无功出力的调节量ΔQ_g下发到发电厂的各自动电压控制电厂子站，实现基于母线电压历史数据分析的逆调压控制。

Claims (1)

1.一种基于母线电压历史数据分析的逆调压控制方法，其特征在于，该方法包括以下步骤：

(1)从电力系统的自动电压控制系统中，在每日0时到来时，读取电力系统中的每条母线的名称，形成母线名称集合B：B_i＝|B₁ B₂ B₃ B₄...B_I|，i＝1,...,I，B为母线名称集合，下标i为母线序号，B_i为第i个母线；读取电力系统中发电机组的名称，形成发电机组名称的集合G：G_j＝|G₁ G₂ G₃ G₄...G_J|，j＝1,...,J，G为发电机组名称集合，下标j为发电机组序号，G_j为第j个发电机组；读取电力系统中容抗器的名称，形成容抗器的名称集合C，C_k＝|C₁ C₂C₃ C₄...C_K|，k＝1,...,K，C为容抗器名称集合，下标k为容抗器序号，C_k为第k个容抗器；

(2)从电力系统中读取母线电压设定计划值以及母线、发电机组和容抗器的历史值，设定一天内的采样点为n＝1,...,N、历史天数m＝1,...,M，历史天数m为从当前0点前的M天中的一天，发电机组对母线的灵敏度和变电站低压母线对高压母线的灵敏度，容抗器容量，进行数据整合形成数据集合，包括以下步骤：

(2-1)根据母线集合B，读取母线电压设定计划值，形成如下所有母线电压设定计划值的集合：

V^max、V^min分别为母线电压设定计划值的上限集合和下限集合，

分别为第i个母线在一天内的第n个采样点的电压设定计划值的上限数据值和电压设定计划值的下限数据值。

(2-2)根据母线名称集合B，读取第i个母线的前M天每天采样点N个点的电压历史值，形成数据集合：

V_i ^m-n为第i个母线的第m天第n个电压历史采样点，遍历母线序号i，形成所有母线的电压历史数据集合V：V＝|V₁ V₂ V₃ V₄...V_I|；

(2-3)根据发电机组集合G，读取第j个发电机组的前M天每天采样点N个点的无功历史值，形成数据集合：

为第j个发电机组的第m天第n个无功历史采样点，遍历发电机组序号j，形成所有发电机组的无功历史数据集合：

Q＝|Q₁ Q₂ Q₃ Q₄...Q_J|；

(2-4)根据容抗器集合C，读取第k个容抗器的前M天每天采样点N个点的运行状态历史值，形成数据集合：

为第k个发电机组的第m天第n个运行状态历史采样点，遍历容抗器序号k，形成所有容抗器的运行状态历史数据集合：D＝|D₁ D₂ D₃ D₄...D_K|；

(2-5)根据发电机组集合G和母线集合B，读取发电机组对母线的灵敏度信息，形成集合如下：

其中，S_GV为发电机组对母线的灵敏度矩阵，S_ji为第j个发电机组对第i个母线的灵敏度；

(2-6)读取变电站低压母线对高压母线的灵敏度，形成变电站低压母线对高压母线的灵敏度集合S_CV：

其中，S_CV为变电站低压母线对高压母线的灵敏度矩阵，S_ii为第w个低压母线对第r个高压母线的灵敏度；w为变电站低压母线个数，r为变电站高压母线个数；

(2-7)根据容抗器集合C，读取第k个容抗器的容量

形成所有容抗器的容量集合：

(3)根据步骤(2)得到的数据集合，计算发电厂控制母线、变电站控制母线的逆调压电压上限，步骤如下：

(3-1)读取母线集合的母线B_i，对母线B_i进行判断，若B_i为发电厂母线，则执行步骤(3-2),若B_i为变电站母线，则执行步骤(3-3)；

(3-2)根据母线集合的母线B_i，从发电机组集合G中读取与母线B_i相连接的发电机组，形成集合G_G：G_f＝|G₁ G₂ G₃ G₄...G_F|、f＝1,...,F，其中，F为与母线B_i相连接的发电机组的个数，G_f为第f个发电机组，并根据相连接的发电机组集合G_G分别从发电机组的无功历史数据集合Q中获取该发电机组的无功历史值Q_f，同时从母线B_i电压设定计划值集合V^max、V^min中分别读取母线B_i的设定计划值上限

设定计划值下限

根据上述读取到的集合数据，计算母线B_i的逆调压电压上限V_i ^smax，包括以下步骤：

(3-2-1)设定母线B_i电压向下调节能力V_id＝0、电压向上调节能力V_iu＝0，计算母线B_i的电压向下调节能力V_id、电压向上调节能力V_iu，步骤如下：

(3-2-1-1)设定发电机组G_f可减少无功Q_fs＝0、发电机组G_f可减少无功的采样计数点Z₁＝0，记发电机组G_f的最小无功出力为Q_min，从发电机组的无功历史值Q_f中读取

其中，上标m-n为第m天的第n个采样点、

为第m天的第n个点的无功值，将发电机组G_f的第m天第n个点的历史值

与发电机组最小无功出力Q_min进行比较，若

则使发电机组G_f的可减少无功

使采样计数点Z₁＝Z₁+1；遍历历史天数m，将所有得到的发电机组G_f的可减少无功Q_fs进行累加，并将所有发电机组G_f可减少无功的采样计数点Z₁的累加；根据累加值，计算发电机组G_f可减少无功值的平均值Q_avg＝Q_fs/Z₁；读取发电机组G_f对母线B_i的灵敏度S_GV，根据灵敏度S_GV和发电机组G_f可减少无功平均值Q_avg计算得到母线B_i电压向下调节能力V_adjd＝Q_avg×S_GV×V_x,V_x为安全系数，根据电压向下调节能力V_adjd，计算得到母线B_i的电压向下调节能力累加：V_id＝V_id+V_adjd；

(3-2-1-2)设定发电机组G_f的可增加无功Q_fu＝0、发电机组G_f可增加无功的采样计数点Z₂＝0，将发电机组G_f的最大无功出力记为Q_max，从发电机组G_f的无功历史值Q_f中读取

将发电机组的第m天第n个点的历史值

与发电机组最大无功出力进行比较，若

则使

使Z₂＝Z₂+1；遍历历史天数m，将所有得到的发电机组的可增加无功Q_fu进行累加，并将发电机组可增加无功的采样计数点Z₂进行累加，根据累加值计算得到发电机组G_f可增加无功平均值Q_pju＝Q_fu/Z₂；读取发电机组G_f对母线B_i的灵敏度S_GV，根据灵敏度S_GV和发电机组G_f可增加无功平均值Q_avg，计算得到母线B_i电压向上调节能力V_adju＝Q_pju×S_GV×V_x,V_x为安全系数；据电压向上调节能力V_adju，计算得到母线B_i的电压向上调节能力累加：V_iu＝V_iu+V_adju；

(3-2-1-3)对与母线B_i相连接的发电机组集合G_G，遍历发电机组个数f的取值，重复步骤(3-2-1-1)～步骤(3-2-1-2)，得到所有与母线B_i相连接的发电机组对母线B_i总的电压向下调节能力V_id和电压向上调节能力V_iu；

(3-2-2)设定母线B_i的电压值V_s＝0，母线B_i电压计数点Z₃＝0，从步骤(2-2)的母线B_i的电压历史值V_i中读取V_i ^m-n，将母线B_i第m天第n个点的历史值V_i ^m-n与母线B_i的电压设定计划值

比较，若

并且

则使母线B_i电压V_s＝V_s+V_i ^m-n，并使计数点Z₃＝Z₃+1；遍历历史天数m，对得到的母线B_i的所有电压值V_s进行累加，对所有母线B_i的电压计数点Z₃进行累加，根据得到的累加数据计算得到母线B_i的电压平均值V_ipj＝V_s/Z₃；

(3-2-3)对采样点n进行判断，若采样点n为电力系统负荷低谷时刻，则根据(3-2-2)计算得到的母线B_i的电压平均值V_ipj和(3-2-1-3)计算得到的母线B_i总的向下调节能力V_id，计算得到母线B_i的逆调压电压上限

并读取母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

若采样点n的时刻为电力系统负荷平峰或高峰时刻，则根据步骤(3-2-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-2-1-3)计算得到的母线B_i总的电压向上调节能力V_iu，计算母线B_i的逆调压电压上限

并读取母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

(3-3)根据母线集合的母线B_i，从容抗器集合C中读取与母线B_i相连接的容抗器C_C＝|G₁G₂ G₃ G₄...G_O|，其中o＝1,...,O，O为相连接的容抗器的个数，根据相连接的容抗器集合C_C，从容抗器G_o的无功历史数据集合D中获取容抗器的历史值D_o，从母线B_i电压设定计划值集合V^max、V^min中读取母线B_i的设定计划值上限

设定计划值下限

根据上述读取到的集合数据，计算母线B_i的逆调压电压上限V_i ^smax，包括以下步骤：

(3-3-1)设定母线B_i的电压向下调节能力V_id＝0、电压向上调节能力V_iu＝0，计算母线B_i的电压向下调节能力V_id、电压向上调节能力V_iu，步骤如下：

(3-3-1-1)设定母线B_i可增加无功容量Q_vh＝0、母线B_i可增加无功容量的采样计数点Z₄＝0、母线B_i可减无功容量Q_vl＝0、母线B_i可减无功容量的采样计数点Z₅＝0,从步骤(2-4)容抗器G_o的运行状态历史值D_o中读取

并从步骤(2-7)容抗器的容量集合C^rmx中读取

对G_o的容量

进行判断：若

则判定G_o为电容器，并进而根据电容器G_o的运行状态的历史值

进行判断，若C_o为已投入运行的电容器，则使

Q_vl＝Q_vl+1，若C_o为热备用的电容器，则使

使Z₄＝Z₄+1；若

为负，则判定G_o为电抗器，并进而根据电抗器G_o的运行状态的历史值

进行判断，若C_o为已投入运行的电抗器，则使

Z₄＝Z₄+1，若C_o为热备用的电抗器，则使

Q_vl＝Q_vl+1；遍历历史天数m,得到母线B_i所有可增加无功容量Q_vh的累加值、母线B_i所有可增加无功容量的采样计数点Z₄的累加值、母线B_i所有可减无功容量Q_vl的累加值以及母线B_i所有可减无功容量的采样计数点Z₅的累加值；

(3-3-1-2)根据步骤(3-3-1-1)的母线B_i所有可增加无功容量Q_vh的累加值、B_i所有可增加无功容量的采样计数点Z₄的累加值，计算母线B_i的可增加无功容量平均值Q_pjh＝Q_vh/Z₄，根据母线B_i所有可减无功容量Q_vl的累加、母线B_i所有可减无功容量的采样计数点Z₅的累加值，计算得到母线B_i的可减少无功容量平均值Q_pjl＝Q_vl/Z₅；从步骤(2-5)读取变电站低压母线对高压母线灵敏度S_VV，根据变电站低压母线对高压母线灵敏度S_VV和可减少无功容量平均值Q_pjl、可增加无功容量平均值Q_pjh，分别计算得到母线B_i的电压向下调节能力V_tjd＝Q_pjl×S_VV×V_x和电压向上调节能力V_tju＝Q_pjh×S_VV×V_x，V_x为安全系数；根据计算得到的母线B_i的电压向下调节能力V_tjd、电压向上调节能力V_tju，分别计算得到母线B_i的累加电压向下调节能力V_id＝V_id+V_tjd和电压向上调节能力V_iu＝V_iu+V_tju；

(3-3-1-3)对与母线B_i相连接的容抗器集合C_C，遍历容抗器个数o的取值，重复步骤(3-2-1-1)～步骤(3-2-1-2)，得到与母线B_i的所有连接的容抗器对母线B_i的总向下调节能力V_id、总电压向上调节能力V_iu；

(3-3-2)设定母线B_i电压V_s＝0，母线B_i电压采样点数Z₃＝0，从母线B_i的电压历史值V_i中读取V_i ^m-n，将母线B_i第m天第n个点的历史值V_i ^m-n与母线B_i的电压设定计划值

比较，若

且

则使母线B_i的电压V_s＝V_s+V_i ^m-n、使Z₃＝Z₃+1；遍历历史天数m，得到母线B_i的所有电压V_s的累加值、母线B_i电压采样点数Z₃的累加值，根据累加值计算得到母线B_i的电压平均值V_ipj＝V_s/Z₃；

(3-3-3)对采样点n进行判断，若采样点n为电力系统负荷低谷时刻，则根据步骤(3-3-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-3-1-3)计算得到的母线B_i总的向下调节能力V_id，计算得到母线B_i的逆调压电压上限

并进而根据母线B_i的电压设定计划值

对逆调压电压上限进行如下修正：若

则使

若

则使

若采样点n的时刻为电力系统负荷平峰或高峰时刻，则根据步骤(3-3-2)计算得到的母线B_i的电压平均值V_ipj和步骤(3-3-1-3)计算得到的母线B_i总的电压向上调节能力V_iu，计算得到母线B_i的逆调压电压上限

并进而根据母线B_i的电压设定计划值

对母线B_i的逆调压电压上限进行如下修正：若

则使

若

则使

(3-4)遍历步骤(1)的母线集合B中的所有母线，遍历采样点n，重复步骤(3-1)～步骤(3-3)，计算并形成母线B_i的逆调压电压上限集合：

为第i个母线第n个采样点的逆调压电压上限值；

(4)设定自动电压控制周期为T_c，在每个自动电压控制周期T_c到来时，记到来时刻为t₀，计算自动电压控制发电厂无功控制指令，包括以下步骤：

(4-1)建立一个自动电压控制优化模型，该优化模型的目标函数的表达式如下：

其中，V_p表示电力系统区域中枢母线的当前电压，

表示电力系统区域中枢母线设定电压，C_g为中枢母线无功电压灵敏度矩阵，ΔQ_g表示发电机无功出力的调节量，W_p和W_q为发电机组的权重系数，取值范围为(0.1-10)，α为增益系数，取值范围为0.1～10，Θ_g为无功裕度向量，Θ_g的第i个分量为：

其中，

为发电机i的当前无功出力，

为发电机i的无功出力的调节量，

为发电机i的无功最大值，

为发电机i的无功最小值，i为发电机的序号，||Θ_g||²为无功裕度向量集，以上数据量分别从电力系统中自动电压控制系统中获取；

(4-2)确定自动电压控制优化模型的约束条件，表达式如下：

其中，C₀、C₁、C₂、C₃分别为限值约束方程，C_vg为高压侧母线无功电压灵敏度矩阵，从电力系统中读取，V_H表示发电机高压侧母线的当前电压，

和

分别表示中枢母线电压下限和中枢母线电压上限；

和

分别表示发电机无功下限和无功上限；

和

分别表示发电机高压侧母线的电压下限、电压上限和允许的单步电压最大调整量，以上数据量分别从电力系统中自动电压控制系统中获取；

(4-3)根据(2)定义的一天内的采样点n＝1,...,N，计算采样点时间间隔t_t＝24×60÷N，通过当前到来时刻t₀及采样点时间间隔t_t，计算得到当前采样点n₀＝t₀÷t_t；从步骤(3-4)逆调压电压上限集合V^smax中，获取与当前到来t₀时刻相对应的逆调压电压上限

并以该

分别取代步骤(4-2)中限值约束方程C₁中的

和限值约束方程C₂中的

得到更新后的自动电压控制优化模型；

(4-4)采用二次规划计算方法，求解步骤(4-3)的优化模型，得到发电机无功出力的调节量ΔQ_g，并将该发电机无功出力的调节量ΔQ_g下发到发电厂的各自动电压控制电厂子站，实现基于母线电压历史数据分析的逆调压控制。

Images