CN105610156B - 一种多线路并发合环方法 - Google Patents

Abstract

本发明公开了一种多线路并发合环方法，包括：获取主配网模型，并对主配网模型进行拼接；对拼接后的主配网模型进行网络拓扑分析，通过确定母线为节点，母线间的连接设备为支路，建立节点‑支路的电网模型；获取各节点的节点电压，各节点的节点电压包括各节点的节点电压相位及节点电压幅值；获取并发合环导纳矩阵，并发合环导纳矩阵包括常规节点导纳矩阵及并发合环增量矩阵；基于并发合环导纳矩阵、各节点的节点电压计算各合环开关的合环电流；判断合环电流是否越保护限额，如果否，则将主配网模型中的各合环开关同时闭合。由此，大大提高了负荷转移效率，降低了负荷转移时间，并且，进一步避免了诱发大面积停电的发生。

Description

一种多线路并发合环方法

技术领域

本发明涉及电网调度技术领域，更具体地说，涉及一种多线路并发合环方法。

背景技术

目前，随着城市电网负荷密度加大，电网环网化率与电缆化率越来越高，手拉手供电，多回路供电等已经成为配网设计的常态。当电网中的线路发生故障时，通过配电自动化系统进行线路合环，可以使得负荷快速转供到对侧线路，从而有效避免停电。

当变电站面临全停风险(如暴雨导致变电站内水位上升、一条母线检修而另一条母线发生故障等)时，必须在尽可能短的时间内把负荷快速转移出去，现有技术中应对这种情况的技术方案为逐条线路依次合环，但是发明人发现，这种方案存在负荷转移效率低，转移时间过长，容易诱发大面积停电的问题。

综上所述，现有技术中存在负荷转移效率低、转移时间过长及容易诱发大面积停电的问题。

发明内容

本发明的目的是提供一种多线路并发合环方法，以解决现有技术中存在的负荷转移效率低、转移时间过长及容易诱发大面积停电的问题。

为了实现上述目的，本发明提供如下技术方案：

一种多线路并发合环方法，包括：

获取主配网模型，并对所述主配网模型进行拼接；

对拼接后的所述主配网模型进行网络拓扑分析，并通过确定母线为节点，母线间的连接设备为支路，建立节点-支路的电网模型；

获取各节点的节点电压，其中，所述各节点的节点电压包括各节点的节点电压相位及节点电压幅值；

获取并发合环导纳矩阵，所述并发合环导纳矩阵包括常规节点导纳矩阵及并发合环增量矩阵；

基于所述并发合环导纳矩阵、所述各节点的节点电压计算各合环开关的合环电流；

判断所述合环电流是否越保护限额，如果否，则将所述主配网模型中的各合环开关同时闭合。

优选的，获取并发合环导纳矩阵，包括：

获取所述并发合环导纳矩阵，所述并发合环导纳矩阵包括常规节点导纳矩阵及并发合环增量矩阵，其中，所述并发合环增量矩阵表示为下式：

其中，ΔY表示所述并发合环增量矩阵，ΔY_ii表示节点i的自导纳变化量，ΔY_jj表示节点j的自导纳变化量，ΔY_ij表示节点i和节点j之间的互导纳变化量，ΔY_ji表示节点j与节点i之间的互导纳变化量，节点j为与节点i拟合环的节点。

优选的，分别按照下列公式计算ΔY_ii、ΔY_jj、ΔY_ij及ΔY_ji：

ΔY_jj＝y_ij

ΔY_ij＝-y_ij

ΔY_ji＝-y_ij

其中，y_ij为节点i和节点j之间支路的导纳，n为与节点i拟合环的节点j的总数量。

优选的，获取ΔY_ii、ΔY_jj、ΔY_ij及ΔY_ji，包括：

如果所述节点j与所述节点i位于同一个变电站内，则y_ij取为0。

优选的，基于所述并发合环导纳矩阵、所述各节点的节点电压计算各合环开关的合环电流，并确定该合环电流为合环潮流，包括：

按照下列公式计算及

其中，Y'为所述并发合环导纳矩阵，为任意合环开关的合环电流，及分别对应分布于该合环开关两侧的节点i的节点电压及节点j的节点电压；

按照下列公式计算任意合环开关的合环电流：

其中，为节点i及节点j的电压差。

优选的，获取各节点的节点电压相位，包括：

将各节点的节点电压幅值的标幺值取为1，利用Givens正交变换法，基于各支路对应的支路有功潮流或者各节点对应的节点有功注入计算各节点的节点电压相位。

优选的，获取各节点的节点电压幅值，包括：

利用部分节点的节点电压幅值量测获取该部分节点的节点电压幅值，基于其余部分节点的零注入电流虚部量测获取该部分节点的节点电压幅值。

优选的，所述基于其余部分节点的零注入电流虚部量测获取该部分节点的节点电压幅值，包括：

对于所述其余部分节点中的任意节点i，按照下列公式计算该结点i的节点电压幅值：

其中，为节点i的注入电流相量，U_j为节点j的节点电压幅值，θ_ij为节点i和节点j之间的电压相位差，Y_ij为所述常规节点导纳矩阵中第i行第j列元素，g_ij为Y_ij的实部，b_ij为Y_ij的虚部，N为节点的总数量，img(·)为取虚部函数。

本发明提供的一种多线路并发合环方法，包括：获取主配网模型，并对所述主配网模型进行拼接；对拼接后的所述主配网模型进行网络拓扑分析，并通过确定母线为节点，母线间的连接设备为支路，建立节点-支路的电网模型；获取各节点的节点电压，其中，所述各节点的节点电压包括各节点的节点电压相位及节点电压幅值；获取并发合环导纳矩阵，所述并发合环导纳矩阵包括常规节点导纳矩阵及并发合环增量矩阵；基于所述并发合环导纳矩阵、所述各节点的节点电压计算各合环开关的合环电流；判断所述合环电流是否越保护限额，如果否，则将所述主配网模型中的各合环开关同时闭合。通过本申请公开的上述技术特征，首先基于并发合环导纳矩阵及各节点的节点电压计算出各合环开关的合环电流，进而在判断出合环电流并未越保护限额时直接将各合环开关同时闭合，由此，在变电站面临全停风险时，能够在尽可能短的时间内将负荷快速转移出去，与现有技术中逐条线路依次合环相比，本申请通过上述方案实现了多条线路的并发合环，大大提高了负荷转移效率，降低了负荷转移时间，并且，进一步避免了诱发大面积停电的发生。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

图1为本发明实施例提供的一种多线路并发合环方法的流程图；

图2为本发明实施例提供的合环前两条线路的合环等值电路；

图3为本发明实施例提供的合环时两条线路的合环等值电路；

图4为本发明实施例提供的合环后反向电压单独作用时的两条线路的合环等值电路；

图5为本发明实施例提供的合环后将反向电压替换成一等值电流时两条线路的合环等值电路。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

请参阅图1，其示出了本发明实施例提供的一种多线路并发合环方法的流程图，可以包括以下步骤：

S11：获取主配网模型，并对主配网模型进行拼接。

需要说明的是，配电网中负荷的大小和分布位置，将影响到配电网中电压的分布，使得合环开关两侧存在一定的电势差，当进行合环操作(即将合环开关进行闭合操作)时，在上述电势差的作用下，环网中将产生环流，导致配电网的潮流分布发生变化，因此，网络结构及负荷分布的不对称性所引起的合环开关两侧的电势差是合环潮流(对应于合环电流)产生的主要原因。

由于合环潮流与合环点所在供区内环路阻抗的分布有关，而环路阻抗受主变阻抗影响较大，因此，为了准确计算配电网中各线路的环路阻抗，需要对主配网模型进行拼接。其中，主网模型在输出时通常可以采用符合IEC61970的CIME格式，配电网模型通常可以只包括10kV出线模型，在拼接时需要依据变电站名称与调度编号对主配网模型进行拼接。现有技术中对于上述拼接已有较成熟的技术方案，在此不再赘述。

S12：对拼接后的主配网模型进行网络拓扑分析，并通过确定母线为节点，母线间的连接设备为支路，建立节点-支路的电网模型。

需要说明的是，网络拓扑分析是解决电力系统分析计算问题的基础，其主要任务是处理开关信息的变化，形成新的电网结构，其实质是把母线-开关之间描述的物理模型通过分析转化为节点-支路之间描述的数学模型。具体来说，在进行网络拓扑分析时，采用矩阵法，把电力系统物理模型中的连接关系抽象成图论中的图，即把物理模型中的节点或母线抽象成节点，把开关或支路抽象成图中的边，这样图中的节点与节点之间、节点与边之间的关系可以通过邻接矩阵和关联矩阵表示出来，邻接矩阵反映节点对之间是否通过边直接连接，关联矩阵是反映节点和边之间的连接关系。本发明实施例中，将进行拼接后的主配网模型中的母线抽象成节点，即负荷节点；而母线之间的连接设备抽象成支路，其中，母线之间的连接设备可以包括主变电器、线路间的联络开关、电缆支路、架空线路、开关站等。本申请中，支路可以分为联络支路与非联络支路，联络支路是指具备合环功能的支路，非联络支路是指不具备合环功能的支路。

S13：获取各节点的节点电压，其中，各节点的节点电压包括各节点的节点电压相位及节点电压幅值。

S14：获取并发合环导纳矩阵，并发合环导纳矩阵包括常规节点导纳矩阵及并发合环增量矩阵。

需要说明的是，常规节点导纳矩阵是现有技术中的运行方式所对应的导纳矩阵，在此不再赘述。而并发合环增量矩阵是由多线路并发合环引起的增量，以进一步实现对于多线路并发合环所产生的合环电流的计算。

S15：基于并发合环导纳矩阵、各节点的节点电压计算各合环开关的合环电流。

S16：判断合环电流是否越保护限额，如果是，则执行步骤S17，如果否，则执行步骤S18。

其中，判断合环电流是否越保护限额即为判断合环电流是否大于保护限额，保护限额可以由工作人员根据实际需要进行确定，如果合环电流大于保护限额，说明合环电流过大，如果此时合环的话可能会导致整条线路停电，因此，此时不进行并发合环，而如果合环电流不大于保护限额，说明合环电流符合要求，可以进行并发合环。

S17：结束操作。

S18：将主配网模型中的各合环开关同时闭合。

将主配网模型中的各合环开关同时闭合，即实现多条线路的并发合环，不同于现有技术中的逐条线路依次合环，本申请中多条线路的合环是同时进行的。

通过本申请公开的上述技术特征，首先基于并发合环导纳矩阵及各节点的节点电压计算出各合环开关的合环电流，进而在判断出合环电流并未越保护限额时直接将各合环开关同时闭合，由此，在变电站面临全停风险时，能够在尽可能短的时间内将负荷快速转移出去，与现有技术中逐条线路依次合环相比，本申请通过上述方案实现了多条线路的并发合环，大大提高了负荷转移效率，降低了负荷转移时间，并且，进一步避免了诱发大面积停电的发生。

上述实施例提供的一种多线路并发合环方法中，获取各节点的节点电压相位，可以包括：

将各节点的节点电压幅值的标幺值取为1，利用Givens正交变换法，基于各支路对应的支路有功潮流或者各节点对应的节点有功注入计算各节点的节点电压相位。

在实际供、配电网中，合环潮流的大小与合环开关两侧的电压相量差关系密切。为求得合环开关两侧的电压相量差(简称电压差)，除了需要知道合环开关两侧的电压幅值外，还需知道其相位。由于实际供、配电网中电压相位一般没有测量，不能直接获取，因此，本申请中利用电网中有功和无功相对解耦实现对于合环开关两侧的电压相量差的计算，且有功量测的精度及可靠性均相对较好，因此，采用有功量测估计合环开关两侧对应的节点的节点电压相位，可减小对基础数据及量测质量的依赖性，提高合环电流算法的可靠性及计算速度。

具体来说，本申请中在对各节点的节点电压相位进行估计时，只考虑有功功率量测(支路量测和注入量测)和节点电压相位之间的关系，而忽略其节点电压幅值的变化，因此，将节点电压幅值的标幺值均取为1。而估计节点电压相位的状态量为各节点的节点电压相位，量测量为支路有功潮流或节点有功注入。

其中，对于与节点i及节点j对应的支路有功潮流量测，有：

P_ij＝B_ij(θ_i-θ_j)

i为任意支路的始端节点对应的编号，j为该支路的末端节点对应的编号，P_ij为该支路由始端节点i流向末端节点j的有功功率，B_ij为常规节点导纳矩阵第i行第j列元素的虚部，θ为对应节点的节点电压相位。

对于节点i的有功注入量测P_i，有：

n为与节点i对应的支路的总数量，即可以为以节点i为始端节点的支路的总数量。

由上述两个公式可以看出，节点的相位估计问题的量测方程均为状态量的线性函数，故相位估计问题为加权线性最小二乘估计问题，其数值求解不需要迭代，具有很高的可靠性。

需要说明的是，加权线性最小二乘估计问题在数学上可描述为如下的极小化问题：

min J(x)＝[z-h(x)]^TW[z-h(x)]

其中，x为状态向量，n维；z为量测向量，m维；h为量测方程向量，m维；W为量测权重矩阵，m×m维，一般为正定对角矩阵。一般而言，量测值的测量精度决定了其在目标函数中权因子的大小，精度越高，所赋予的权因子越大。实际量测的精度主要取决于测量仪表本身；伪量测的误差一般较大；虚拟量测则绝对准确，不存在误差。

求解加权线性最小二乘估计问题的方法主要有正则方程法及正交变换法，其中，正交变换法相对于正则方程法来说，具有极好的数值稳定性，而正交变换方法主要有Givens正交变换法及Household正交变换法，对于大规模稀疏矩阵，相对于Household正交变换法，Givens正交变换具有更高的效率，对于实际的大规模电网，加权线性相位估计问题的加权量测对应矩阵为大规模稀疏矩阵，因此，采用Givens正交变换法比Household正交变换法具有更高的执行效率，故本发明实施例中利用Givens正交变换法求解加权线性相位估计问题，以获得节点的节点电压相位。

上述实施例提供的一种多线路并发合环方法中，获取各节点的节点电压幅值，可以包括：

利用部分节点的节点电压幅值量测获取该部分节点的节点电压幅值，基于其余部分节点的零注入电流虚部量测获取该部分节点的节点电压幅值。

采用叠加原理计算合环潮流，首先需要计算合环前合环开关两侧的电压相量差，由上述实施例求得合环开关两侧的节点电压相位，还需要求得合环开关两侧的节点电压幅值，才能求得电压相量差。对于主配网模型中的母线联络开关或母线分段开关，开关两侧一般均有电压幅值量测，且电压幅值量测的可靠性及精度一般均较高，因此可以直接求得开关两侧的电压相量差。而对于元件开关，其母线侧一般有电压幅值量测，但元件侧一般没有电压幅值量测，为求得此类开关两侧的电压相量差，需对其电压幅值进行估计。

具体来说，对于具有N个节点的电网模型，如果网络结构和网络元件参数已知，则可得网络方程：

其中，Y为N×N阶的常规节点导纳矩阵，为N×1维节点电压列矢量，为N×1维节点注入电流列矢量。

对于节点i，有：

其中，为节点i的注入电流相量，U_j为节点j的节点电压幅值，θ_ij为节点i和节点j之间的相位差，Y_ij为常规节点导纳矩阵的第i行第j列元素，g_ij为Y_ij的实部，b_ij为Y_ij的虚部。需要说明的是，节点i与节点j为对应于同一支路的两个节点。

若节点i为零注入节点，即上述的其余部分节点中的任意节点，其注入电流相量为0，则有：

因此，有：

由此可见，对于零注入节点i，可以添加零注入量测：

式中，img(·)为取虚部函数。

由于各节点的节点电压相位已经求得，因此系数b_ijcosθ_ij+g_ijsinθ_ij为已知量，进而可以求得各节点的节点电压幅值。

由此，可以利用节点电压幅值量测及零注入电流虚部量测构建电压幅值估计器，对各节点的节点电压幅值进行估计。显然，由这两类量测构建的电压幅值估计器为线性估计器。

若电网模型中所有零注入节点均不配置电压幅值量测，所有非零注入节点均配置有电压幅值量测，则上述电压幅值估计器可以满足可观测性要求。对于实际输、供电网，非零注入节点(电源点或负荷点)一般配置有电压幅值量测，零注入节点则部分配置有电压幅值量测，部分未配置有电压幅值量测，且配置电压幅值量测的比例一般较高，对于少数未配置电压幅值量测的节点，可添加虚拟的零注入电流伪量测，并以较小的权重参与电压幅值估计，以保证电压幅值估计器的可观测性。由此，不但满足可观测性要求，而且具有一定的冗余度。

上述实施例提供的一种多线路并发合环方法中，获取并发合环导纳矩阵，可以包括：

获取并发合环导纳矩阵，并发合环导纳矩阵包括常规节点导纳矩阵及并发合环增量矩阵，其中，并发合环增量矩阵表示为下式：

其中，ΔY表示并发合环增量矩阵，ΔY_ii表示节点i的自导纳变化量，ΔY_jj表示节点j的自导纳变化量，ΔY_ij表示节点i和节点j之间的互导纳变化量，ΔY_ji表示节点j与节点i之间的互导纳变化量，节点j为与节点i拟合环的节点。需要说明的是，此处的节点i和节点j与上述实施例中的节点i和节点j并不矛盾，这是因为，当节点i和节点j对应的同一支路为具有合环开关的支路时，节点j则为与节点i拟合环的节点。

需要说明的是，与两条线路之间合环电流的计算不同的是，本申请中的并发合环导纳矩阵Y'，与现有技术中的常规节点导纳矩阵不同在于其包含了与节点i拟合环的所有节点之间的导纳值，即包含了一个并发合环增量矩阵，并且同一个变电站内的两个拟合环的节点间的导纳值为0，即Y'＝Y+ΔY，其中Y为常规节点导纳矩阵，ΔY表示并发合环增量矩阵。

上述ΔY的表达式中，ΔY_ii为节点i的自导纳变化量：

其中，y_ij为节点i和与节点i拟合环的节点j之间支路的导纳，n为与节点i拟合环的节点j的总数量。

ΔY_ij为节点i和与节点i拟合环的节点j之间的互导纳变化量：

ΔY_ij＝-y_ij

ΔY_jj为与节点i拟合环的节点j的自导纳变化量：

ΔY_jj＝y_ij

ΔY_ji为与节点i拟合环的节点j与节点i之间的互导纳变化量：

ΔY_ji＝-y_ij

需要说明的是，如果节点i和节点j位于同一个变电站内，y_ij＝0，对应的，ΔY_ii、ΔY_jj、ΔY_ij及ΔY_ji均为0。

上述实施例提供的一种多线路并发合环方法中，基于并发合环导纳矩阵、各节点的节点电压计算各合环开关的合环电流，并确定该合环电流为合环潮流，可以包括：

按照下列公式计算及

其中，Y'为并发合环导纳矩阵，为任意合环开关的合环电流，及分别对应分布于该合环开关两侧的节点i的节点电压及节点j的节点电压；需要说明的是，此处的节点i和节点j与上述实施例中的节点i和节点j并不矛盾，这是因为，当节点i和节点j对应的同一支路为具有合环开关的支路时，节点j则为与节点i拟合环的节点，对应的节点i和节点j即为分布于该合环开关两侧的节点。

按照下列公式计算任意合环开关的合环电流：

其中，为节点i及节点j的电压差。

具体来说，多条线路之间的并发合环可以看作多对两条线路之间的并发合环，与两条线路之间的合环不同的是其合环导纳矩阵的变化，可以将每个合环开关的合环过程全部利用两条线路合环中的计算步骤分别计算合环电流。

其中，对两条线路之间的合环电流的计算，需要对合环点(即合环开关)以外的电网模型进行等值，根据叠加原理把电网模型中每条线路上的潮流值都由合环前的初始值和合环操作引起的增量两部分组成。设化简后的两条线路的合环等值电路如图2至图5所示。在合环之前，节点i与节点j之间相当于串联一个电压源如图2所示；合环时，节点i与节点j之间短路，相当于在节点i与j之间又串联一个反向电压源将节点i与节点j之间的电势差拉至零，如图3所示。根据叠加定理，合环后各线路上的电势值分为两个部分：一部分是合环前合环开关两侧电压差和有源网络中的电源共同起作用的结果，对应于合环前的状态；另一部分是单独作用引起的结果，如图4所示。由于原本节点i与节点j之间开路，节点i与节点j之间流过的电流为零，合环时流过合环点的合环电流全部由引起。根据替换定理将反向电压源替换成一等值电流源其大小正是要求的合环电流，如图5所示。

对于具有N个节点的电网模型，如果网络结构和网络元件参数已知，则可得网络方程：

其中，Y为N×N阶常规节点导纳矩阵，是N×1维节点电压列矢量，是N×1维节点注入电流列矢量。

将上述改写成增量型式，则有：

由图5可知，等值电流源在节点i处注入电流在节点j处注入电流故有：

即：

由于多条线路并发合环时，位于同一个变电站内两个节点之间进行合环时，合环电流比逐条线路依次合环时的合环电流要大得多，在多条线路并发合环时，需要排除同一个变电站内两个拟合环的节点之间进行合环时对其它线路合环之间的影响。由于与两条线路之间合环计算不同的是，此处的Y矩阵为并发合环导纳矩阵，即Y'，因此，上式可以改写为：

而Y'已由上述实施例中获得，因此，对该式进行求解，可得

由于合环前合环开关两侧的电压差为合环后合环开关两侧的电压差为0，故有：

由此可求得合环电流：

在多条线路进行并发合环时，由于对应的多个合环开关同时被闭合，其合环前的电压差互不影响，合环后的电压差也均为0，故可以对每个合环开关分别计算合环电流。

举例说明，假设某变电站AB有两段10kV母线，I段母线与II段母线分列运行，I段母线位于节点A，II段母线位于节点B，与I段母线相联的手拉手馈线四条，分别通过北埠216线与节点B合环，双流218线与节点C合环，凌家319线与节点D合环，里龙613线与节点E合环，黄少608线与节点F合环。

采用依次合环时，即按照黄沙608线—〉北埠216线—〉凌家319线—〉里龙613线—〉双流218线的顺序进行合环时，合环电流分别为表1；采用本发明实施例提供的一种多条线路并发合环的方法对其进行并发合环，其合环电流分别为表2；

可知，节点A与节点B因位于同一个变电站，采用并发合环时，其合环电流极剧增大，而其它环路合环电流减少，为减少节点B对节点A的合环影响，对双流218线剔除并发合环后其它线路进行并发合环计算结果如表3所示。

可知，如果排除双流218线进行并发合环时，其它线路合环潮流比包含双流218线进行并发合环的潮流值变小。所以，在进行并发合环时，应避免位于同一变电站内的不同节点之间的线路合环对其它线路合环的影响。

表1

表2

表3

合环时间	合环开关	环路	环路电流(Amp)
				0:33:36	rdf:ID＝"284113195"	凌家319线	316.257956
0:33:36	rdf:ID＝"284043074"	里龙613线	204.675050
				0:33:36	rdf:ID＝"284113141"	北埠216线	103.668838
0:33:36	rdf:ID＝"284042771"	黄沙608线	164.779634

对所公开的实施例的上述说明，使本领域技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (6)

1.一种多线路并发合环方法，其特征在于，包括：

获取主配网模型，并对所述主配网模型进行拼接；

对拼接后的所述主配网模型进行网络拓扑分析，并通过确定母线为节点，母线间的连接设备为支路，建立节点-支路的电网模型；

获取各节点的节点电压，其中，所述各节点的节点电压包括各节点的节点电压相位及节点电压幅值；

获取并发合环导纳矩阵，所述并发合环导纳矩阵包括常规节点导纳矩阵及并发合环增量矩阵；

基于所述并发合环导纳矩阵、所述各节点的节点电压计算各合环开关的合环电流；

判断所述合环电流是否越保护限额，如果否，则将所述主配网模型中的各合环开关同时闭合；

获取并发合环导纳矩阵，包括：

获取所述并发合环导纳矩阵，所述并发合环导纳矩阵包括常规节点导纳矩阵及并发合环增量矩阵，其中，所述并发合环增量矩阵表示为下式：

其中，ΔY表示所述并发合环增量矩阵，ΔY_ii表示节点i的自导纳变化量，ΔY_jj表示节点j的自导纳变化量，ΔY_ij表示节点i和节点j之间的互导纳变化量，ΔY_ji表示节点j与节点i之间的互导纳变化量，节点j为与节点i拟合环的节点；

分别按照下列公式计算ΔY_ii、ΔY_jj、ΔY_ij及ΔY_ji：

ΔY_jj＝y_ij

ΔY_ij＝-y_ij

ΔY_ji＝-y_ij

其中，y_ij为节点i和节点j之间支路的导纳，n为与节点i拟合环的节点j的总数量。

2.根据权利要求1所述的方法，其特征在于，获取ΔY_ii、ΔY_jj、ΔY_ij及ΔY_ji，包括：

如果所述节点i和所述节点j位于同一个变电站内，则y_ij取为0。

3.根据权利要求1所述的方法，其特征在于，基于所述并发合环导纳矩阵、所述各节点的节点电压计算各合环开关的合环电流，并确定该合环电流为合环潮流，包括：

按照下列公式计算及

其中，Y'为所述并发合环导纳矩阵，为任意合环开关的合环电流，及分别对应分布于该合环开关两侧的节点i的节点电压及节点j的节点电压；

按照下列公式计算任意合环开关的合环电流：

其中，为节点i及节点j的电压差。

4.根据权利要求1所述的方法，其特征在于，获取各节点的节点电压相位，包括：

将各节点的节点电压幅值的标幺值取为1，利用Givens正交变换法，基于各支路对应的支路有功潮流或者各节点对应的节点有功注入计算各节点的节点电压相位。

5.根据权利要求1所述的方法，其特征在于，获取各节点的节点电压幅值，包括：

利用部分节点的节点电压幅值量测获取该部分节点的节点电压幅值，基于其余部分节点的零注入电流虚部量测获取该部分节点的节点电压幅值。

6.根据权利要求5所述的方法，其特征在于，所述基于其余部分节点的零注入电流虚部量测获取该部分节点的节点电压幅值，包括：

对于所述其余部分节点中的任意节点i，按照下列公式计算该节 点i的节点电压幅值：

其中，为节点i的注入电流相量，U_j为节点j的节点电压幅值，θ_ij为节点i和节点j之间的电压相位差，Y_ij为所述常规节点导纳矩阵中第i行第j列元素，g_ij为Y_ij的实部，b_ij为Y_ij的虚部，N为节点的总数量，img(·)为取虚部函数。