CN112865108A - 一种基于连续潮流仿真的电网自动电压控制仿真方法 - Google Patents

Abstract

本发明属于电网调度技术领域，尤其涉及一种基于连续潮流仿真的电网自动电压控制仿真方法。本发明方法利用电网实时运行的模型和数据，从整个系统的角度，通过仿真潮流计算形成电网数据和AVC策略的闭环控制，通过仿真数据分析对电网最佳无功电压调整方案给出建议和指导，以电压安全和优质为约束，以系统运行经济性为目标，连续闭环的进行电压的实时优化控制。通过连续循环执行本方法，实现多种典型日断面条件下电网无功电压自动控制仿真，为电网多种负荷状态下自动电压控制参数的设置提共参考和依据，优化无功电压自动控制策略，提高无功电压自动控制水平。

Description

一种基于连续潮流仿真的电网自动电压控制仿真方法

技术领域

本发明属于电网调度技术领域，尤其涉及一种基于连续潮流仿真的电网自动电压控制仿真方法。

背景技术

自动电压控制(以下简称AVC，Automatic Voltage Control)系统是实现输电网安全(提高电压稳定裕度)、经济(降低网络损耗)、优质(提高电压合格率)运行的重要手段。AVC系统架构在电网能量管理系统(以下简称EMS)之上，能够利用输电网实时运行数据，从输电网全局优化的角度科学决策出最佳的无功电压调整方案，自动下发给电厂、变电站以及下级电网调度机构执行。

在我国的电网系统中电压是电能质量的重要指标，合理的电压无功控制不仅可以提高电网的电压合格水平，还能够较好的减少电网的整体损耗。在实际应用中,随着电网的快速建设以及用电量的快速变化，电网调度往往不能及时的根据电网负荷变化调整无功控制参数，存在一些电压无功控制不合理的问题。

电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的数据，计算母线的电压、各元件的功率及网损，并对电网各处的运行状态进行评估。再根据计算得到的数据对电网系统的运行进行监测和优化。对于正在运行的电力系统，通过潮流计算可以评估当前系统中母线的电压、支路的功率等参数是否超限；如果出现异常，就应采取措施，调整运行方式。

发明内容

本发明目的是提出一种基于连续潮流仿真的电网自动电压控制仿真方法，以克服已有技术的缺点，利用电网实时运行的模型和数据，从整个系统的角度，通过仿真潮流计算形成电网数据和AVC策略的闭环控制。

本发明提出的基于连续潮流仿真的电网自动电压控制仿真方法，包括以下步骤：

(1)设定电网自动电压控制周期为T_c；

(2)从电网的电压自动控制的历史断面中，读取设定参考日的电网模型M和288点基态电网潮流断面F_i，形成以下两种仿真文件：F_m，i电网288点运行方式文件以及F_sch电网方式调整数据文件，其中：

288点电网运行方式文件F_m，i的表达式为：

F_m，i＝[i，i＝1，..288]{Ld_m，i，Un_m，i，Bs_m，i，Xf_m，i，Dev_m，i}；

其中，下标i为采样时刻，Ld_m，i、Un_m，i、Bs_m，i、Xf_m，i和Dev_m，i分别为与参考日的第i采样时刻相对应的电网模型中负荷Ld_m，i、发电机组Un_m，i、电网母线Bs_m，i、电网变压器Xf_m，i和设备Dev_m，i的信息数组，其中设备包括开关和变压器分接头；

电网负荷信息数组包含：负荷有功Ld_m，i，p和负荷无功Ld_m，i，q两部分，即：

Ld_m，i＝[Ld_m，i，p，Ld_m，i，q]；

发电机组信息数组包含：发电机组有功Un_m，i，p和发电机组无功Un_m，i，q两部分，即：

Un_m，i＝[Un_m，i，p，Un_m，i，q]；

电网母线信息数组包含：母线名称Bs_m，i，name和母线电压Bs_m，i，val两部分，即：

Bs_m，i＝[Bs_m，i，name，Bs_m，i，val]；

电网变压器信息数组包含：变压器有功Xf_m，i，p和变压器无功Xf_m，i，q两部分，即：

Xf_m，i＝[Xf_m，i，p，Xf_m，i，q]；

设备信息数组包含：开关状态Cb_m，i和变压器分接头Oltc_m，i两部分，即：

Dev_m，i＝[Cb_m，i，Oltc_m，i]；

其中，开关状态Cb_m，i包含：开关名称Cb_m，i，name和开关状态Cb_m，i，val，即：

Cb_m，i＝[Cb_m，i，name，Cb_m，i，val]；

变压器分接头Oltc_m，i包含：分接头名称Oltc_m，i，name和分接头挡位Oltc_m，i，val，即：

Oltc_m，i＝[Oltc_{m，i，name，}Oltc_m，i，val]；

F_sch电网方式调整数据文件的表达式如下：

F_sch＝{BS_sch}；

其中，BS_sch为电网母线中需调节母线的调整信息，设需调节母线的序号为n，n＝1......N，N表示文件F_sch中纪录的需调节母线的总数量，即：

BS_sch＝[n，n＝1，..N][bus_n，U_n，min，U_n，max，T_n，start，T_n，end]；

其中，bus_n为第n条记录的需调节母线的名称，U_n，min为第n条记录的需调节母线的电压下限值，U_n，max为第n条记录的需调节母线的电压上限值，T_n，start，T_n，end分别为第n条记录的需调节母线的调整生效的开始时刻和结束时刻；

(3)建立一个初始电网模型F_N：设定初始电网模型为参考日第0:00时刻的基态电网潮流断面F₀，包括电网发电机的机端电压、发电机有功，负荷有功、负荷无功、开关状态和变压器分接头挡位，具体过程如下：

(3-1)设定电网中可调节发电机组的序号为g，g＝1......G，G表示电网模型M中可调节发电机组的总数量，G_0，v和G_0，p分别为可调节发电机组的初始电压和初始有功：

G_0，v＝[g，g＝1，..G][V_0，g]；

G_0，p＝[g，g＝1，..G][P_0，g]；

其中，V_0，g和P_0，g分别为第g个可调节发电机组的初始电压和初始有功；

(3-2)设定电网中可调节负荷的序号为l，l＝1......L，L表示电网模型M中可调节负荷的总数量，L_0，p和L_0，q为可调节负荷的初始有功和初始无功：

L_0，p＝[l，l＝1，..L][P_0，l]；

L_0，q＝[l，l＝1，..L][Q_0，l]；

其中，P_0，l和Q_0，l分别为第l个可调节负荷的初始有功和初始无功；

(3-3)设定电网中可控开关的序号为c，c＝1......C，C表示电网模型M中可控开关的总数量，Cb₀为可控开关的初始状态数组：

Cb₀＝[c，c＝1，..C][Open_0，c]；

其中，Open_0，c为第c个可控开关的初始分合状态：

(3-4)设定电网中变压器可控分接头的序号为t，t＝1......T，T表示电网模型M中变压器可控分接头的总数量，Oltc₀为变压器可控分接头的初始挡位数组：

Oltc₀＝[t，t＝1，..T][Tap_0，t]；

其中，Tap_0，t为初始第t个变压器可控分接头的挡位：

(3-5)设定电网电压自动控制中可控制母线的序号为s，s＝1......S，S表示电网模型M中可控制母线的总数量，Bs₀为可控制母线的初始状态数组：

U_s，i，min和U_s，i，max分为与第s条可控制母线在第i个采样时刻相对应的初始上限值和初始下限值；

(3-6)根据以上描述，得到初始电网模型F_N如下：

F_N＝[G_0，v，G_0，p，L_0，p，L_0，q，Cb₀，Oltc₀，Bs₀]；

(4)在仿真启动的第0:00时刻，根据步骤(2)中的电网方式调整数据文件F_sch，对步骤(3-5)可控制母线的初始状态数组Bs₀进行更新，得到新的可控制母线电压状态数组Bs_0，new，包括以下步骤：

(4-1)对步骤(2)的电网方式调整数据文件F_sch进行检测，若电网方式调整数据文件F_sch中记录有需调节母线的信息，则进行步骤(4-2)，若电网方式调整数据文件F_sch中没有记录需调节母线的信息，则进行步骤(5)；

(4-2)根据步骤(2)中的电网方式调整数据文件F_sch，对步骤(3-5)的可控制母线的初始状态数组Bs₀进行更新，更新过程如下：

(4-2-1)从步骤(2)的电网方式调整文件F_sch中读取第n条需调节母线的名称bus_n和电压调整限值U_n，min，U_n，max，并将该需调节母线调整生效的开始时刻和结束时刻记为T_{n，i，start}，T_n，i，end；

(4-2-2)根据步骤(4-2-1)的需调节母线信息，对步骤(3-5)的可控制母线的初始状态数组Bs₀进行更新，方法如下：

(4-2-2-1)根据步骤(4-2-1)中得到的第n条需调节母线的名称bus_n，从步骤(3-5)的可控制母线的初始状态数组Bs₀中，检索到第n条需调节母线在数组Bs₀中的位置，将该位置记为s₀；

(4-2-2-2)用步骤(4-2-1)中读取的需调节母线的电压调整限值U_n，min，U_n，max、调整生效的开始时刻和结束时刻T_{n，i，start}，T_n，i，end和步骤(4-2-2-1)中得到的需调节母线序号s0，对步骤(3-5)中可控制母线初始状态数组Bs₀的U_n，min，U_n，max进行替换，得到时段T_{n，i，start}，T_n，i，end内可控制母线新的电压限值如下：

进而得到一个新的可控制母线电压状态数组Bs_0，new如下：

(5)在参考日第0:00时刻，对基态电网进行潮流计算，输出得到潮流信息数组F_0，out：

F_0，out＝[Bs_0，out，Xf_0.out，Cp_0，out，Un_0，out，Cb_0，out，Ln_0，out，Ld_0，out]；

其中，Bs_0，out为第0:00时刻电网母线信息，Xf_0.out为第0:00时刻电网变压器信息，Cp_0，out为第0:00时刻电网容抗器信息，Un_0，out为第0:00时刻电网发电机组信息，Cb_0，out为第0:00时刻电网开关信息，Ln_0，out为第0:00时刻电网线路信息，Ld_0，out为第0:00时刻电网负荷信息；

(6)根据步骤(5)的潮流计算结果，在第0:00时刻进行自动电压控制计算，输出得到电厂母线、容抗器、变压器分接头的相应控制参数A_0，strg：

A_0，strg＝[DVC_0，stg，SVC_0，strg]；

其中，DVC_0，stg为第0:00时刻的变电站控制参数，SVC_0，strg为第0:00时刻的电厂母线控制参数，其中：

变电站控制参数DVC_0，stg为：

DVC_0，stg＝[[x，x＝0，..X]Cp_0，x，strg，[y，y＝0，..Y]Oltc_0，y，strg]；

其中，x为容抗器控制参数记录序号，x＝0......X，X表示容抗器控制参数记录总数量，y为变压器分接头控制参数记录序号，y＝0......Y，Y为变压器分接头控制参数记录总数量，Cp_0，x，strg为第x条容抗器控制参数，Oltc_0，y，strg为第y条变压器分接头控制参数；

电厂母线控制参数SVC_0，strg为：

SVC_0，strg＝[[z，z＝0，..Z]Bs_0，z，strg]；

其中，z为电厂母线控制参数记录序号，z＝0......Z，Z表示电厂母线控制参数记录的总数量，Bs_0，z，strg为第z条电厂母线控制参数；

(7)根据步骤(6)的电厂母线、容抗器、变压器分接头的相应控制参数A_0，strg，对步骤(3)的初始电网模型F_N进行更新，并利用更新后的初始电网模型FN进行潮流计算，包括以下步骤：

(7-1)从步骤(6)的变电站控制参数DVC_0，stg中读取第x₀条容抗器控制参数

并对容抗器控制参数

进行判断，若容抗器控制参数

为投入运行，则设置与容抗器相对应的开关状态

为闭合，若容抗器控制参数

为退出运行，则设置与容抗器相对应的开关状态

为分开；根据上述开关状态，在步骤(3)的初始电网模型FN查找该开关在可控开关初始状态数组Cb₀中的位置，并将该位置记为c₀，对初始电网模型FN中的开关状态进行更新：

(7-2)从步骤(6)的变电站控制参数DVC_0，stg中读取第y₀条变压器分接头控制参数

得到变压器分接头对应的挡位值

在步骤(3)的初始电网模型F_N中查找第y₀条变压器分接头在初始挡位数组Oltc₀中的位置，并将该位置记为t₀，并对初始电网模型F_N中的变压器分接头挡位进行更新：

(7-3)从步骤(6)的电厂母线控制参数SVC_0，strg中读取第z₀条电厂母线控制参数

得到电厂母线设定电压目标值为

其中，

为电厂母线当前电压值，

为电厂母线设定电压目标值，计算得到电厂母线的电压调整量

对上述

进行判断，若

则判定电厂母线的电压调节为增电压，若

则判定电厂母线的电压调节为减电压；根据上述

判断结果，统计初始电网模型FN中该电厂可调节发电机组的个数，包括以下步骤：

(a)若

则统计初始电网模型F_N中第z₀条控制参数所在电厂可调节增电压的发电机组G_0，up的总数量U₁：

U_0，adjup＝[u₁，u₁＝1，..U₁]G_0，up；

设电厂母线的当前电压

与第u₁个电厂发电机组的机端电压

的比值为

根据上述比值

得到第0:00时刻第u₁个电厂发电机组的机端电压调整目标值为

在步骤(3)的初始电网模型F_N中查找第u₁个电厂发电机组在可调节发电机组初始状态数组G_0，v中的位置，并将该位置记为g₁，根据上述机端电压调整目标值

和位置g₁，对初始电网模型F_N中可调节发电机组的初始电压进行更新如下：

(b)若

则统计初始电网模型F_N中第z₀条控制参数所在电厂可调节减电压的发电机组G_0，down的总数量U₂：

U_0，adjdn＝[u₂，u₂＝1，..U₂]G_0，down；

设电厂母线的当前电压

与第u₂个电厂发电机组的机端电压

的比值为

根据上述比值

得到第0:00时刻第u₂个电厂发电机组的机端电压调整目标值为

在步骤(3)的初始电网模型F_N中查找第u₂个电厂发电机组在可调节发电机组初始状态数组G_0，v中的位置，并将该位置记为g₂，根据上述机端电压调整目标值

和位置g₂，对初始电网模型F_N中可调节发电机组的初始电压进行更新如下：

(7-4)利用上述更新后的初始电网模型FN进行潮流计算，输出得到潮流信息数组F′_0，out：

F′_0，out＝[Bs′_0，out，Xf′_0.out，Cp′_0，out，Un′_0，out，Cb′_0，out，Ln′_0，out，Ld′_0，out]；

(7-5)从步骤(7-4)的数组F′_0，out中，得到与自动电压控制设备相关的变压器分接头Oltc_strg，out数组、与容抗器对应的开关Cb_strg，cp数组和发电机组机端电压G_strg，out数组；

(8)在参考日第二采样时刻，读取相应的电网运行方式文件F_m，1，并在当前电网潮流结构中更新，包括以下步骤：

(8-1)使当前电网潮流计算模型中使与容抗器相连的开关以外的其他开关的运行为第二采样时刻F_m，1中状态，具体更新如下：

根据电网运行方式文件F_m，1，使与容抗器相连的开关以外的其他开关运行状态为：

Cb₀＝Cb_m，1；

使与容抗器相连的开关以及变压器分接头保持步骤(7-5)中的运行状态，即：

Cb₀＝Cb_strg，cp；

Oltc₀＝Oltc_strg，out；

(8-2)使电厂发电机组的机端电压保持步骤(7-5)中的运行状态，即：

G_0，v＝G_strg，out；

(8-3)根据电网运行方式文件F_m，1，使电厂发电机组有功、可控负荷有功、可控负荷无功为：

G_0，p＝Un_m，1，p；

Ld_0，p＝Ld_m，1，p；

Ld_0，q＝Ld_m，1，q；

(9)利用第二采样时刻调整后的电网结构，重复步骤(5)，进行潮流计算，输出得到潮流信息数组F_1，out：

F_1，out＝[Bs_1，out，Xf_1.out，Cp_1，out，Un_1，out，Cb_1，out，Ln_1，out，Ld_1，out]；

其中：Bs_1，out为第二采样时刻电网母线信息，Xf_1.out为第二采样时刻电网变压器信息，Cp_1，out为第二采样时刻电网容抗器信息，Un_1，out为第二采样时刻电网机组信息，Cb_1，out为第二采样时刻电网开关信息，Ln_1，out为第二采样时刻电网线路信息，Ld_1，out为第二采样时刻电网负荷信息；

(10)根据步骤(9)得到的潮流信息数组F_1，out，重复步骤(6)，计算得到第二采样时刻的自动电压控制参数数组A_1strg：

A_1，strg＝[DVC_1，stg，SVC_1，strg]；

DVC_1，stg＝[[x，x＝0，..X]Cp_1，x，strg，[y，y＝0，..Y]Oltc_1，y，strg]；

SVC_1，strg＝[[z，z＝0，..Z]Bs_1，z，strg]；

(11)用步骤(10)的自动电压控制参数数组A_1strg，对步骤(3)的初始电网模型F_N进行更新，重复步骤(7)，利用更新后的初始电网模型FN进行潮流计算，输出得到潮流信息数组F′_1，out；

(12)遍历所有采样时刻，重复步骤(8)～步骤(11)，输出仿真结果如下：

a、全天电网运行数据：

F_out＝[f，f＝1，..F][Bs_f，out，Xf_f.out，Cp_f，out，Un_f，out]；

其中，f为电网运行数据序号，f＝0......F，F表示电网运行数据的总数量，Bs_f，out为第f个断面电网运行数据中母线信息，Xf_f，out为第f个断面电网运行数据中主变信息，Cp_f，out为第f个断面电网运行数据中容抗器信息，Un_f，out为第f个断面电网运行数据中机组信息。

b、全天自动电压控制参数：

A_strg＝[k，k＝1，..K][Bs_k，stg，Cp_k，strg，Oltc_k，strg]；

其中，k为自动电压控制参数序号，k＝0......K，K表示自动电压控制参数的总数量。Bs_k，stg为第k条自动电压控制参数对应的电厂母线调整信息，Cp_k，strg为第k条自动电压控制参数对应的容抗器投切信息，Oltc_k，strg为第k条自动电压控制参数对应的主变分头调整信息；

实现基于连续潮流仿真的电网自动电压控制仿真。

本发明提出的基于连续潮流仿真的电网自动电压控制仿真方法，其优点是：

本发明基于连续潮流仿真的电网自动电压控制仿真方法，利用电网实时运行的模型和数据，从整个系统的角度，通过仿真潮流计算形成电网数据和AVC策略的闭环控制，通过仿真数据分析对电网最佳无功电压调整方案给出建议和指导，以电压安全和优质为约束，以系统运行经济性为目标，连续闭环的进行电压的实时优化控制。本发明方法可以通过连续循环执行，实现多种典型日断面条件下电网无功电压自动控制仿真，为电网多种负荷状态下自动电压控制参数的设置提共参考和依据，优化无功电压自动控制策略，提高无功电压自动控制水平。因此本发明方法克服了传统电网无功电压管理机制中存在的不足，有效地解决电网实际运行中的电压控制问题。本方法通过仿真潮流计算能够将电网的典型日断面和基于实际AVC程序给出的策略进行闭环模拟控制，通过对仿真输出结果数据的统计和分析，对实际电网无功电压控制中调压参数的设置给出积极建议和指导，有效地克服传统电网无功电压管理机制中存在的不足，确保电网安全运行，给客户带来更优质的用户体验。

附图说明

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

图2是本发明方法中步骤(7)的流程框图。

具体实施方式

本发明提出的基于连续潮流仿真的电网自动电压控制仿真方法，其流程框图如图1所示，包括以下步骤：

(1)设定电网自动电压控制周期为T_c；

(2)从电网的电压自动控制的历史断面中，读取设定参考日的电网模型M和288点基态电网潮流断面F_i，形成以下两种仿真文件：F_m，i电网288点运行方式文件以及F_sch电网方式调整数据文件，其中：

288点电网运行方式文件F_m，i的表达式为：

F_m，i＝[i，i＝1，..288]{Ld_m，i，Un_m，i，Bs_m，i，Xf_m，i，Dev_m，i}；

其中，下标i为采样时刻，Ld_m，i、Un_m，i、Bs_m，i、Xf_m，i和Dev_m，i分别为与参考日的第i采样时刻(本发明的一个实施例中每5分钟为一个采样时刻)相对应的电网模型中负荷Ld_m，i、发电机组Un_m，i、电网母线Bs_m，i、电网变压器Xf_m，i和设备Dev_m，i的信息数组，其中设备包括开关和变压器分接头；

电网负荷信息数组包含：负荷有功Ld_m，i，p和负荷无功Ld_m，i，q两部分，即：

Ld_m，i＝[Ld_m，i，p，Ld_m，i，q]；

发电机组信息数组包含：发电机组有功Un_m，i，p和发电机组无功Un_m，i，q两部分，即：

Un_m，i＝[Un_m，i，p，Un_m，i，q]；

电网母线信息数组包含：母线名称Bs_m，i，name和母线电压Bs_m，i，val两部分，即：

Bs_m，i＝[BS_m，i，name，Bs_m，i，val]；

电网变压器信息数组包含：变压器有功Xf_m，i，p和变压器无功Xf_m，i，q两部分，即：

Xf_m，i＝[Xf_m，i，p，Xf_m，i，q]；

设备信息数组包含：开关状态Cb_m，i和变压器分接头Oltc_m，i两部分，即：

Dev_m，i＝[Cb_m，i，Oltc_m，i]；

其中，开关状态Cb_m，i包含：开关名称Cb_m，i，name和开关状态Cb_m，i，val，即：

Cb_m，i＝[Cb_m，i，name，Cb_m，i，val]；

变压器分接头Oltc_m，i包含：分接头名称Oltc_m，i，name和分接头挡位Oltc_m，i，val，即：

Oltc_m，i＝[Oltc_m，i，name，Oltc_m，i，val]；

F_sch电网方式调整数据文件的表达式如下：

F_sch＝{BS_sch}；

其中，BS_sch为电网母线中需调节母线的调整信息，设需调节母线的序号为n，n＝1......N，N表示文件F_sch中纪录的需调节母线的总数量，即：

BS_sch＝[n，n＝1，..N][bus_n，U_n，min，U_n，max，T_n，start，T_n，end]；

其中，bus_n为第n条记录的需调节母线的名称，U_n，min为第n条记录的需调节母线的电压下限值，U_n，max为第n条记录的需调节母线的电压上限值，T_n，start，T_n，end分别为第n条记录的需调节母线的调整生效的开始时刻和结束时刻；

(3)建立一个初始电网模型F_N：设定初始电网模型为参考日第0:00时刻的基态电网潮流断面F₀，包括电网发电机的机端电压、发电机有功，负荷有功、负荷无功、开关状态和变压器分接头挡位，具体过程如下：

(3-1)设定电网中可调节发电机组的序号为g，g＝1......G，G表示电网模型M中可调节发电机组的总数量，G_0，v和G_0，p分别为可调节发电机组的初始电压和初始有功：

G_0，v＝[g，g＝1，..G][V_0，g]；

G_0，p＝[g，g＝1，..G][P_0，g]；

其中，V_0，g和P_0，g分别为第g个可调节发电机组的初始电压和初始有功；

(3-2)设定电网中可调节负荷的序号为l，l＝1......L，L表示电网模型M中可调节负荷的总数量，L_0，p和L_0，q为可调节负荷的初始有功和初始无功：

L_0，p＝[l，l＝1，..L][P_0，l]；

L_0，q＝[l，l＝1，..L][Q_0，l]；

其中，P_0，l和Q_0，l分别为第l个可调节负荷的初始有功和初始无功；

(3-3)设定电网中可控开关的序号为c，c＝1......C，C表示电网模型M中可控开关的总数量，Cb₀为可控开关的初始状态数组：

Cb₀＝[c，c＝1，..C][Open_0，c]；

其中，Open_0，c为第c个可控开关的初始分合状态：

(3-4)设定电网中变压器可控分接头的序号为t，t＝1......T，T表示电网模型M中变压器可控分接头的总数量，Oltc₀为变压器可控分接头的初始挡位数组：

Oltc₀＝[t，t＝1，..T][Tap_0，t]；

其中，Tap_0，t为初始第t个变压器可控分接头的挡位：

(3-5)设定电网电压自动控制中可控制母线的序号为s，s＝1......S，S表示电网模型M中可控制母线的总数量，Bs₀为可控制母线的初始状态数组：

U_s，i，min和U_s，i，max分为与第s条可控制母线在第i个采样时刻相对应的初始上限值和初始下限值；

(3-6)根据以上描述，得到初始电网模型FN如下：

FN＝[G_0，v，G_0，p，L_0，p，L_0，q，Cb₀，Oltc₀，Bs₀]；

(4)在仿真启动的第0:00时刻，根据步骤(2)中的电网方式调整数据文件F_sch，对步骤(3-5)可控制母线的初始状态数组Bs₀进行更新，得到新的可控制母线电压状态数组Bs_0，new，包括以下步骤：

(4-1)对步骤(2)的电网方式调整数据文件F_sch进行检测，若电网方式调整数据文件F_sch中记录有需调节母线的信息，则进行步骤(4-2)，若电网方式调整数据文件F_sch中没有记录需调节母线的信息，则进行步骤(5)；

(4-2)根据步骤(2)中的电网方式调整数据文件F_sch，对步骤(3-5)的可控制母线的初始状态数组Bs₀进行更新，更新过程如下：

(4-2-1)从步骤(2)的电网方式调整文件F_sch中读取第n条需调节母线的名称bus_n和电压调整限值U_n，min，U_n，max，并将该需调节母线调整生效的开始时刻和结束时刻记为T_{n，i，start}，T_n，i，end；

(4-2-2)根据步骤(4-2-1)的需调节母线信息，对步骤(3-5)的可控制母线的初始状态数组Bs₀进行更新，方法如下：

(4-2-2-1)根据步骤(4-2-1)中得到的第n条需调节母线的名称bus_n，从步骤(3-5)的可控制母线的初始状态数组Bs₀中，检索到第n条需调节母线在数组Bs₀中的位置，将该位置记为s₀；

(4-2-2-2)用步骤(4-2-1)中读取的需调节母线的电压调整限值U_n，min，U_n，max、调整生效的开始时刻和结束时刻T_{n，i，start}，T_n，i，end和步骤(4-2-2-1)中得到的需调节母线序号s₀，对步骤(3-5)中可控制母线初始状态数组Bs₀的U_n，min，U_n，max进行替换，得到时段T_{n，i，start}，T_n，i，end内可控制母线新的电压限值如下：

进而得到一个新的可控制母线电压状态数组Bs_0，new如下：

(5)在参考日第0:00时刻，对基态电网进行潮流计算，输出得到潮流信息数组F_0，out：

F_0，out＝[Bs_0，out，Xf_0.out，Cp_0，out，Un_0，out，Cb_0，out，Ln_0，out，Ld_0，out]；

其中，Bs_0，out为第0:00时刻电网母线信息，Xf_0.out为第0:00时刻电网变压器信息，Cp_0，out为第0:00时刻电网容抗器信息，Un_0，out为第0:00时刻电网发电机组信息，Cb_0，out为第0:00时刻电网开关信息，Ln_0，out为第0:00时刻电网线路信息，Ld_0，out为第0:00时刻电网负荷信息；

(6)根据步骤(5)的潮流计算结果，在第0:00时刻进行自动电压控制计算，输出得到电厂母线、容抗器、变压器分接头的相应控制参数A_0，strg：

A_0，strg＝[DVC_0，stg，SVC_0，strg]；

其中，DVC_0，stg为第0:00时刻的变电站控制参数，SVC_0，strg为第0:00时刻的电厂母线控制参数，其中：

变电站控制参数DVC_0，stg为：

DVC_0，stg＝[[x，x＝0，..X]Cp_0，x，strg，[y，y＝0，..Y]Oltc_0，y，strg]；

其中，x为容抗器控制参数记录序号，x＝0......X，X表示容抗器控制参数记录总数量，y为变压器分接头控制参数记录序号，y＝0......Y，Y为变压器分接头控制参数记录总数量，Cp_0，x，strg为第x条容抗器控制参数，Oltc_0，y，strg为第y条变压器分接头控制参数；

电厂母线控制参数SVC_0，strg为：

SVC_0，strg＝[[z，z＝0，..Z]Bs_0，z，strg]；

其中，z为电厂母线控制参数记录序号，z＝0......Z，Z表示电厂母线控制参数记录的总数量，Bs_0，z，strg为第z条电厂母线控制参数；

(7)根据步骤(6)的电厂母线、容抗器、变压器分接头的相应控制参数A_0，strg，对步骤(3)的初始电网模型F_N进行更新，并利用更新后的初始电网模型F_N进行潮流计算，其流程如图2所示，包括以下步骤：

(7-1)从步骤(6)的变电站控制参数DVC_0，stg中读取第x₀条容抗器控制参数

并对容抗器控制参数

进行判断，若容抗器控制参数

为投入运行，则设置与容抗器相对应的开关状态

为闭合，若容抗器控制参数

为退出运行，则设置与容抗器相对应的开关状态

为分开；根据上述开关状态，在步骤(3)的初始电网模型F_N查找该开关在可控开关初始状态数组Cb₀中的位置，并将该位置记为c₀，对初始电网模型F_N中的开关状态进行更新：

(7-2)从步骤(6)的变电站控制参数DVC_0，stg中读取第y₀条变压器分接头控制参数

得到变压器分接头对应的挡位值

在步骤(3)的初始电网模型F_N中查找第y₀条变压器分接头在初始挡位数组Oltc₀中的位置，并将该位置记为t₀，并对初始电网模型F_N中的变压器分接头挡位进行更新：

(7-3)从步骤(6)的电厂母线控制参数SVC_0，strg中读取第z₀条电厂母线控制参数

得到电厂母线设定电压目标值为

其中，

为电厂母线当前电压值，

为电厂母线设定电压目标值，计算得到电厂母线的电压调整量

对上述

进行判断，若

则判定电厂母线的电压调节为增电压，若

则判定电厂母线的电压调节为减电压；根据上述

判断结果，统计初始电网模型FN中该电厂可调节发电机组的个数，包括以下步骤：

(a)若

则统计初始电网模型FN中第z₀条控制参数所在电厂可调节增电压的发电机组G_0，up的总数量U1：

U_0，adjup＝[u₁，u₁＝1，..U₁]G_0，up；

设电厂母线的当前电压

与第u₁个电厂发电机组的机端电压

的比值为

根据上述比值

得到第0:00时刻第u₁个电厂发电机组的机端电压调整目标值为

在步骤(3)的初始电网模型F_N中查找第u₁个电厂发电机组在可调节发电机组初始状态数组G_0，v中的位置，并将该位置记为g₁，根据上述机端电压调整目标值

和位置g₁，对初始电网模型F_N中可调节发电机组的初始电压进行更新如下：

(b)若

则统计初始电网模型FN中第z₀条控制参数所在电厂可调节减电压的发电机组G_0，down的总数量U₂：

U_0，adjdn＝[u₂，u₂＝1，..U₂]G_0，down；

设电厂母线的当前电压

与第u₂个电厂发电机组的机端电压

的比值为

根据上述比值

得到第0:00时刻第u₂个电厂发电机组的机端电压调整目标值为

在步骤(3)的初始电网模型F_N中查找第u₂个电厂发电机组在可调节发电机组初始状态数组G_0，v中的位置，并将该位置记为g₂，根据上述机端电压调整目标值

和位置g₂，对初始电网模型F_N中可调节发电机组的初始电压进行更新如下：

(7-4)利用上述更新后的初始电网模型FN进行潮流计算，输出得到潮流信息数组F′_0，out：

F′_0，out＝[Bs′_0，out，Xf′_0.out，Cp′_0，out，Un′_0，out，Cb′_0，out，Ln′_0，out，Ld′_0，out]；

(7-5)从步骤(7-4)的数组F′_0，out中，得到与自动电压控制设备相关的变压器分接头Oltc_strg，out数组、与容抗器对应的开关Cb_strg，cp数组和发电机组机端电压G_strg，out数组；

(8)在参考日第二采样时刻，读取相应的电网运行方式文件F_m，1，并在当前电网潮流结构中更新，包括以下步骤：

(8-1)使当前电网潮流计算模型中使与容抗器相连的开关以外的其他开关的运行为第二采样时刻F_m，1中状态，具体更新如下：

根据电网运行方式文件F_m，1，使与容抗器相连的开关以外的其他开关运行状态为：

Cb₀＝Cb_m，1；

使与容抗器相连的开关以及变压器分接头保持步骤(7-5)中的运行状态，即：

Cb₀＝Cb_strg，cp；

Oltc₀＝Oltc_strg，out；

(8-2)使电厂发电机组的机端电压保持步骤(7-5)中的运行状态，即：

G_0，v＝G_strg，out；

(8-3)根据电网运行方式文件F_m，1，使电厂发电机组有功、可控负荷有功、可控负荷无功为：

G_0，p＝Un_m，1，p；

Ld_0，p＝Ld_m，1，p；

Ld_0，q＝Ld_m，1，q；

(9)利用第二采样时刻调整后的电网结构，重复步骤(5)，进行潮流计算，输出得到潮流信息数组F_1，out：

F_1，out＝[Bs_1，out，Xf_1.out，Cp_1，out，Un_1，out，Cb_1，out，Ln_1，out，Ld_1，out]；

其中：Bs_1，out为第二采样时刻电网母线信息，Xf_1.out为第二采样时刻电网变压器信息，Cp_1，out为第二采样时刻电网容抗器信息，Un_1，out为第二采样时刻电网机组信息，Cb_1，out为第二采样时刻电网开关信息，Ln_1，out为第二采样时刻电网线路信息，Ld_1，out为第二采样时刻电网负荷信息；

(10)根据步骤(9)得到的潮流信息数组F_1，out，重复步骤(6)，计算得到第二采样时刻的自动电压控制参数数组A_1strg：

A_1，strg＝[DVC_1，stg，SVC_1，strg]；

DVC_1，stg＝[[x，x＝0，..X]Cp_1，x，strg，[y，y＝0，..Y]O1tc_1，y，strg]；

SVC_1，strg＝[[z，z＝0，..Z]Bs_1，z，strg]；

(11)用步骤(10)的自动电压控制参数数组A_1strg，对步骤(3)的初始电网模型FN进行更新，重复步骤(7)，利用更新后的初始电网模型FN进行潮流计算，输出得到潮流信息数组F′_1，out；

(12)遍历所有采样时刻，重复步骤(8)～步骤(11)，输出仿真结果如下：

a、全天电网运行数据：

F_out＝[f，f＝1，..F][Bs_f，out，Xf_f.out，Cp_f，out，Un_f，out]；

其中，f为电网运行数据序号，f＝0......F，F表示电网运行数据的总数量，Bs_f，out为第f个断面电网运行数据中母线信息，Xf_f，out为第f个断面电网运行数据中主变信息，Cp_f，out为第f个断面电网运行数据中容抗器信息，Un_f，out为第f个断面电网运行数据中机组信息。

b、全天自动电压控制参数：

A_strg＝[k，k＝1，..K][Bs_k，stg，Cp_k，strg，Oltc_k，strg]；

其中，k为自动电压控制参数序号，k＝0......K，K表示自动电压控制参数的总数量。Bs_k，stg为第k条自动电压控制参数对应的电厂母线调整信息，Cp_k，strg为第k条自动电压控制参数对应的容抗器投切信息，Oltc_k，strg为第k条自动电压控制参数对应的主变分头调整信息；

实现基于连续潮流仿真的电网自动电压控制仿真。

Claims (1)

1.一种基于连续潮流仿真的电网自动电压控制仿真方法，其特征在于该方法包括以下步骤：

(1)设定电网自动电压控制周期为T_c；

(2)从电网的电压自动控制的历史断面中，读取设定参考日的电网模型M和288点基态电网潮流断面F_i，形成以下两种仿真文件：F_m，i电网288点运行方式文件以及F_sch电网方式调整数据文件，其中：

288点电网运行方式文件F_m，i的表达式为：

F_m，i＝[i，i＝1，..288]{Ld_m，i，Un_m，i，Bs_m，i，Xf_m，i，Dev_m，i}；

其中，下标i为采样时刻，Ld_m，i、Un_m，i、Bs_m，i、Xf_m，i和Dev_m，i分别为与参考日的第i采样时刻相对应的电网模型中负荷Ld_m，i、发电机组Un_m，i、电网母线Bs_m，i、电网变压器Xf_m，i和设备Dev_m，i的信息数组，其中设备包括开关和变压器分接头；

电网负荷信息数组包含：负荷有功Ld_m，i，p和负荷无功Ld_m，i，q两部分，即：

Ld_m，i＝[Ld_m，i，p，Ld_m，i，q]；

发电机组信息数组包含：发电机组有功Un_m，i，p和发电机组无功Un_m，i，q两部分，即：

Un_m，i＝[Un_m，i，p，Un_m，i，q]；

电网母线信息数组包含：母线名称Bs_m，i，name和母线电压Bs_m，i，val两部分，即：

Bs_m，i＝[Bs_m，i，name，Bs_m，i，val]；

电网变压器信息数组包含：变压器有功Xf_m，i，p和变压器无功Xf_m，i，q两部分，即：

Xf_m，i＝[Xf_m，i，p，Xf_m，i，q]；

设备信息数组包含：开关状态Cb_m，i和变压器分接头Oltc_m，i两部分，即：

Dev_m，i＝[Cb_m，i，Oltc_m，i]；

其中，开关状态Cb_m，i包含：开关名称Cb_m，i，name和开关状态Cb_m，i，val，即：

Cb_m，i＝[Cb_m，i，name，Cb_m，i，val]；

变压器分接头Oltc_m，i包含：分接头名称Oltc_m，i，name和分接头挡位Oltc_m，i，val，即：

Oltc_m，i＝[Oltc_m，i，name，Oltc_m，i，val]；

F_sch电网方式调整数据文件的表达式如下：

F_sch＝{BS_sch}；

其中，BS_sch为电网母线中需调节母线的调整信息，设需调节母线的序号为n，n＝1......N，N表示文件F_sch中纪录的需调节母线的总数量，即：

BS_sch＝[n，n＝1，..N][bus_n，U_n，min，U_n，max，T_n，start，T_n，end]；

其中，bus_n为第n条记录的需调节母线的名称，U_n，min为第n条记录的需调节母线的电压下限值，U_n，max为第n条记录的需调节母线的电压上限值，T_n，start，T_n，end分别为第n条记录的需调节母线的调整生效的开始时刻和结束时刻；

(3)建立一个初始电网模型F_N：设定初始电网模型为参考日第0：00时刻的基态电网潮流断面F₀，包括电网发电机的机端电压、发电机有功，负荷有功、负荷无功、开关状态和变压器分接头挡位，具体过程如下：

(3-1)设定电网中可调节发电机组的序号为g，g＝1......G，G表示电网模型M中可调节发电机组的总数量，G_0，v和G_0，p分别为可调节发电机组的初始电压和初始有功：

G_0，v＝[g，g＝1，..G][V_0，g]；

G_0，p＝[g，g＝1，..G][P_0，g]；

其中，V_0，g和P_0，g分别为第g个可调节发电机组的初始电压和初始有功；

(3-2)设定电网中可调节负荷的序号为l，l＝1...…L，L表示电网模型M中可调节负荷的总数量，L_0，p和L_0，q为可调节负荷的初始有功和初始无功：

L_0，p＝[l，l＝1，..L][P_0，l]；

L_0，q＝[l，l＝1，..L][Q_0，l]；

其中，P_0，l和Q_0，l分别为第l个可调节负荷的初始有功和初始无功；

(3-3)设定电网中可控开关的序号为c，c＝1......C，C表示电网模型M中可控开关的总数量，Cb₀为可控开关的初始状态数组：

Cb₀＝[c，c＝1，..C][Open_0，c]；

其中，Open_0，c为第c个可控开关的初始分合状态：

(3-4)设定电网中变压器可控分接头的序号为t，t＝1......T，T表示电网模型M中变压器可控分接头的总数量，Oltc₀为变压器可控分接头的初始挡位数组：

Oltc₀＝[t，t＝1，..T][Tap_0，t]；

其中，Tap_0，t为初始第t个变压器可控分接头的挡位：

(3-5)设定电网电压自动控制中可控制母线的序号为s，s＝1......S，S表示电网模型M中可控制母线的总数量，Bs₀为可控制母线的初始状态数组：

U_s，i，min和U_s，i，max分为与第s条可控制母线在第i个采样时刻相对应的初始上限值和初始下限值；

(3-6)根据以上描述，得到初始电网模型FN如下：

F_N＝[G_0，v，G_0，p，L_0，p，L_0，q，Cb₀，Oltc₀，Bs₀]；

(4)在仿真启动的第0：00时刻，根据步骤(2)中的电网方式调整数据文件F_sch，对步骤(3-5)可控制母线的初始状态数组Bs₀进行更新，得到新的可控制母线电压状态数组Bs_0，new，包括以下步骤：

(4-1)对步骤(2)的电网方式调整数据文件F_sch进行检测，若电网方式调整数据文件F_sch中记录有需调节母线的信息，则进行步骤(4-2)，若电网方式调整数据文件F_sch中没有记录需调节母线的信息，则进行步骤(5)；

(4-2)根据步骤(2)中的电网方式调整数据文件F_sch，对步骤(3-5)的可控制母线的初始状态数组Bs₀进行更新，更新过程如下：

(4-2-1)从步骤(2)的电网方式调整文件F_sch中读取第n条需调节母线的名称bus_n和电压调整限值U_n，oun，U_n，max，并将该需调节母线调整生效的开始时刻和结束时刻记为T_{n，i，start}，T_n，i，end；

(4-2-2)根据步骤(4-2-1)的需调节母线信息，对步骤(3-5)的可控制母线的初始状态数组Bs₀进行更新，方法如下：

(4-2-2-1)根据步骤(4-2-1)中得到的第n条需调节母线的名称bus_n，从步骤(3-5)的可控制母线的初始状态数组Bs₀中，检索到第n条需调节母线在数组Bs₀中的位置，将该位置记为s₀；

(4-2-2-2)用步骤(4-2-1)中读取的需调节母线的电压调整限值U_n，min，U_n，max、调整生效的开始时刻和结束时刻T_{n，i，start}，T_n，i，end和步骤(4-2-2-1)中得到的需调节母线序号s₀，对步骤(3-5)中可控制母线初始状态数组Bs₀的U_n，min，U_n，max进行替换，得到时段T_{n，i，start}，T_n，i，end内可控制母线新的电压限值如下：

进而得到一个新的可控制母线电压状态数组Bs_0，new如下：

(5)在参考日第0：00时刻，对基态电网进行潮流计算，输出得到潮流信息数组F_0，out：

F_0，out＝[Bs_0，out，Xf_0.out，Cp_0，out，Un_0，out，Cb_0，out，Ln_0，out，Ld_0，out]；

其中，Bs_0，out为第0：00时刻电网母线信息，Xf_0.out为第0：00时刻电网变压器信息，Cp_0，out为第0：00时刻电网容抗器信息，Un_0，out为第0：00时刻电网发电机组信息，Cb_0，out为第0：00时刻电网开关信息，Ln_0，out为第0：00时刻电网线路信息，Ld_0，out为第0：00时刻电网负荷信息；

(6)根据步骤(5)的潮流计算结果，在第0：00时刻进行自动电压控制计算，输出得到电厂母线、容抗器、变压器分接头的相应控制参数A_0，strg：

A_0，strg＝[DVC_0，stg，SVC_0，strg]；

其中，DVC_0，stg为第0：00时刻的变电站控制参数，SVC_0，strg为第0：00时刻的电厂母线控制参数，其中：

变电站控制参数DVC_0，stg为：

DVC_0，stg＝[[x，x＝0，..X]Cp_0，x，strg，[y，y＝0，..Y]Oltc_0，y，strg]；

其中，x为容抗器控制参数记录序号，x＝0......X，X表示容抗器控制参数记录总数量，y为变压器分接头控制参数记录序号，y＝0......Y，Y为变压器分接头控制参数记录总数量，Cp_0，x，strg为第x条容抗器控制参数，Oltc_0，y，strg为第y条变压器分接头控制参数；

电厂母线控制参数SVC_0，strg为：

SVC_0，strg＝[[z，z＝0，..Z]Bs_0，z，strg]；

其中，z为电厂母线控制参数记录序号，z＝0......Z，Z表示电厂母线控制参数记录的总数量，Bs_0，z，strg为第z条电厂母线控制参数；

(7)根据步骤(6)的电厂母线、容抗器、变压器分接头的相应控制参数A_0，strg，对步骤(3)的初始电网模型F_N进行更新，并利用更新后的初始电网模型F_N进行潮流计算，包括以下步骤：

(7-1)从步骤(6)的变电站控制参数DVC_0，stg中读取第x₀条容抗器控制参数

并对容抗器控制参数

进行判断，若容抗器控制参数

为投入运行，则设置与容抗器相对应的开关状态

为闭合，若容抗器控制参数

为退出运行，则设置与容抗器相对应的开关状态

为分开；根据上述开关状态，在步骤(3)的初始电网模型F_N查找该开关在可控开关初始状态数组Cb₀中的位置，并将该位置记为c₀，对初始电网模型F_N中的开关状态进行更新：

(7-2)从步骤(6)的变电站控制参数DVC_0，stg中读取第y₀条变压器分接头控制参数

得到变压器分接头对应的挡位值

在步骤(3)的初始电网模型F_N中查找第y₀条变压器分接头在初始挡位数组Oltc₀中的位置，并将该位置记为t₀，并对初始电网模型F_N中的变压器分接头挡位进行更新：

(7-3)从步骤(6)的电厂母线控制参数SVC_0，strg中读取第z₀条电厂母线控制参数

得到电厂母线设定电压目标值为

其中，

为电厂母线当前电压值，

为电厂母线设定电压目标值，计算得到电厂母线的电压调整量

对上述

进行判断，若

则判定电厂母线的电压调节为增电压，若

则判定电厂母线的电压调节为减电压；根据上述

判断结果，统计初始电网模型FN中该电厂可调节发电机组的个数，包括以下步骤：

(a)若

则统计初始电网模型F_N中第z₀条控制参数所在电厂可调节增电压的发电机组G_0，up的总数量U₁：

U_0，adjup＝[u₁，u₁＝1，..U₁]G_0，up；

设电厂母线的当前电压

与第u₁个电厂发电机组的机端电压

的比值为

根据上述比值

得到第0：00时刻第u₁个电厂发电机组的机端电压调整目标值为

在步骤(3)的初始电网模型F_N中查找第u₁个电厂发电机组在可调节发电机组初始状态数组G_0，v中的位置，并将该位置记为g₁，根据上述机端电压调整目标值

和位置g₁，对初始电网模型F_N中可调节发电机组的初始电压进行更新如下：

(b)若

则统计初始电网模型F_N中第z₀条控制参数所在电厂可调节减电压的发电机组G_0，down的总数量U₂：

U_0，adidn＝[u₂，u₂＝1，..U₂]G_0，down；

设电厂母线的当前电压

与第u₂个电厂发电机组的机端电压

的比值为

根据上述比值

得到第0：00时刻第u₂个电厂发电机组的机端电压调整目标值为

在步骤(3)的初始电网模型F_N中查找第u₂个电厂发电机组在可调节发电机组初始状态数组G_0，v中的位置，并将该位置记为g₂，根据上述机端电压调整目标值

和位置g₂，对初始电网模型F_N中可调节发电机组的初始电压进行更新如下：

(7-4)利用上述更新后的初始电网模型FN进行潮流计算，输出得到潮流信息数组F′_0，out：

F′_0，out＝[Bs′_0，out，Xf′_0.out，Cp′_0，out，Un′_0，out，Cb′_0，out，Ln′_0，out，Ld′_0，out]；

(7-5)从步骤(7-4)的数组F′_0，out中，得到与自动电压控制设备相关的变压器分接头Oltc_strg，out数组、与容抗器对应的开关Cb_strg，cp数组和发电机组机端电压G_strg，out数组；

(8)在参考日第二采样时刻，读取相应的电网运行方式文件F_m，1，并在当前电网潮流结构中更新，包括以下步骤：

(8-1)使当前电网潮流计算模型中使与容抗器相连的开关以外的其他开关的运行为第二采样时刻F_m，1中状态，具体更新如下：

根据电网运行方式文件F_m，1，使与容抗器相连的开关以外的其他开关运行状态为：

Cb₀＝Cb_m，1；

使与容抗器相连的开关以及变压器分接头保持步骤(7-5)中的运行状态，即：

Cb₀＝Cb_strg，cp；

Oltc₀＝Oltc_strg，out；

(8-2)使电厂发电机组的机端电压保持步骤(7-5)中的运行状态，即：

G_0，v＝G_strg，out；

(8-3)根据电网运行方式文件F_m，1，使电厂发电机组有功、可控负荷有功、可控负荷无功为：

G_0，p＝Un_m，1，p；

Ld_0，p＝Ld_m，1，p；

Ld_0，q＝Ld_m，1，q；

(9)利用第二采样时刻调整后的电网结构，重复步骤(5)，进行潮流计算，输出得到潮流信息数组F_1，out：

F_1，out＝[Bs_1，out，Xf_1.out，Cp_1，out，Un_1，out，Cb_1，out，Ln_1，out，Ld_1，out]；

其中：Bs_1，out为第二采样时刻电网母线信息，Xf_1.out为第二采样时刻电网变压器信息，Cp_1，out为第二采样时刻电网容抗器信息，Un_1，out为第二采样时刻电网机组信息，Cb_1，out为第二采样时刻电网开关信息，Ln_1，out为第二采样时刻电网线路信息，Ld_1，out为第二采样时刻电网负荷信息；

(10)根据步骤(9)得到的潮流信息数组F_1，out，重复步骤(6)，计算得到第二采样时刻的自动电压控制参数数组A_1strg：

A_1，strg＝[DVC_1，stg，SVC_1，strg]；

DVC_1，stg＝[[x，x＝0，..X]Cp_1，x，strg，[y，y＝0，..Y]Oltc_1，y，strg]；

SVC_1，strg＝[[z，z＝0，..Z]Bs_1，z，strg]；

(11)用步骤(10)的自动电压控制参数数组A_1strg，对步骤(3)的初始电网模型FN进行更新，重复步骤(7)，利用更新后的初始电网模型FN进行潮流计算，输出得到潮流信息数组F′_1，out；

(12)遍历所有采样时刻，重复步骤(8)～步骤(11)，输出仿真结果如下：

a、全天电网运行数据：

F_out＝[f，f＝1，..F][Bs_f，out，Xf_f.out，Cp_f，out，Un_f，out]；

其中，f为电网运行数据序号，f＝0......F，F表示电网运行数据的总数量，Bs_f，out为第f个断面电网运行数据中母线信息，Xf_f，out为第f个断面电网运行数据中主变信息，Cp_f，out为第f个断面电网运行数据中容抗器信息，Un_f，out为第f个断面电网运行数据中机组信息。

b、全天自动电压控制参数：

A_strg＝[k，k＝1，..K][Bs_k，stg，Cpk_，strg，Oltc_k，strg]；

其中，k为自动电压控制参数序号，k＝0......K，K表示自动电压控制参数的总数量。Bs_k，stg为第k条自动电压控制参数对应的电厂母线调整信息，Cp_k，strg为第k条自动电压控制参数对应的容抗器投切信息，Oltc_k，strg为第k条自动电压控制参数对应的主变分头调整信息；

实现基于连续潮流仿真的电网自动电压控制仿真。

Images