A parallel computing method for power distribution network commercial load cluster security domain

Abstract

The application discloses a kind of parallel computing methods of distribution network commercial load cluster security domain, comprising: in each low-voltage distribution network, only considering the operation domain of power flowing into each low-voltage distribution network from medium-voltage bus in low-voltage distribution network constraint parallel computing;In each possible topology of medium-voltage distribution network, the operation domain of power flowing into the low-voltage distribution network connected to each medium-voltage bus is parallel calculated;The security domain of power flowing into the low-voltage distribution network connected to each medium-voltage bus is calculated;The security domain of power flowing into the low-voltage distribution network connected to each medium-voltage bus is refined into the security domain of distribution network commercial load cluster.The application is based on the analysis of distribution network topology, and the calculation of distribution network commercial load cluster security domain is divided into multiple steps, and each step is divided into multiple parallel executable subtasks, which realizes the rapid definition of distribution network commercial load cluster security domain, and supports distribution network commercial load cluster to participate in demand side response.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of power distribution network safety domain analysis and calculation, and particularly relates to a parallel computing method for a commercial load cluster safety domain of a power distribution network. BACKGROUND

[0002] With the rapid growth of renewable energy penetration in power systems, the demand for real-time balancing of power systems is also increasing. In recent years, as an important resource for maintaining the real-time balance of power grid generation and consumption, the demand-side response capability of commercial loads has received increasing attention. Commercial loads contain some flexible loads, and adjusting their power within a certain range will not cause significant impact. Therefore, commercial loads can respond to load increase / decrease instructions from the power dispatch center to assist in maintaining the real-time balance of generation and consumption. Considering the thermal stability constraints of the distribution network, it is necessary to define the safety domain of the commercial load cluster in the state space composed of the power adjustable range of all commercial loads, and limit the power adjustment range of each commercial load within the safety domain to ensure that the distribution network after participating in the demand-side response still meets the N-1 safety operation criterion. The N-1 safety operation criterion means that after any one of the expected accident sets occurs and the power supply is transferred, the power of each feeder section of the distribution network does not exceed the limit.

[0003] According to the granularity from coarse to fine, the safety domain of the distribution network can be divided into the substation perspective, the main transformer perspective, the feeder perspective and the load perspective. The existing distribution network safety domain calculation method adopts the boundary point simulation method, which determines the boundary of the safety domain by simulating the power flow of each working point. Considering the amount of calculation, the existing method can only calculate the safety domain at the first three perspectives. Refining the safety domain to the load perspective will lead to an exponential increase in the amount of calculation, which in turn leads to a long time-consuming for calculating the safety domain, affecting the real-time performance. However, the safety domain at the first three perspectives is too coarse in granularity, which cannot effectively constrain the power adjustment range of each commercial load, and cannot be applied to the scenario of demand-side response of commercial load clusters. SUMMARY

[0004] The present application constructs a parallel computing method for the safety domain of a commercial load cluster of a distribution network to solve the problem that the calculation of the safety domain of the distribution network cannot balance granularity and speed.

[0005] The present application discloses a parallel computing method for the safety domain of a commercial load cluster of a distribution network, which comprises the following steps:

[0006] S1: parallel computing in each low-voltage distribution network, in each low-voltage distribution network, according to the adjustable range of the power of each commercial load and the maximum carrying capacity of each feeder section, the operating domain of the power of each low-voltage feeder section is calculated in order from downstream to upstream, and then the operating domain of the power flowing into the low-voltage distribution network from the medium-voltage bus is obtained only considering the constraints of the low-voltage distribution network;

[0007] S2: calculating in parallel in each possible topology of the medium voltage distribution network, in each possible topology of the medium voltage distribution network, according to the operating domain of the power flowing from each medium voltage bus into the connected low voltage distribution network and the maximum current carrying capacity of each feeder section, in the order from downstream to upstream, the operating domain of the power of each medium voltage feeder section in the topology;

[0008] S3: calculating in parallel in each possible topology of the medium voltage distribution network, in each possible topology of the medium voltage distribution network, according to the operating domain of the power of each medium voltage feeder section, in the order from upstream to downstream, determining the operating domain of the power flowing from each medium voltage bus into the connected low voltage distribution network in the topology;

[0009] S4: taking the overlapping part of the operating domain of the power flowing from each medium voltage bus into the connected low voltage distribution network in all possible topologies of the medium voltage distribution network as the security domain of the power flowing from each medium voltage bus into the connected low voltage distribution network;

[0010] S5: replacing each power variable of the power flowing from each medium voltage bus into the connected low voltage distribution network in the security domain obtained in step S4 with the sum of the power variables of all commercial loads in the connected low voltage distribution network of the medium voltage bus, and refining the security domain of the power flowing from each medium voltage bus into the connected low voltage distribution network into the security domain of the commercial load cluster of the distribution network.

[0011] Step S1 further comprises:

[0012] S1.1: pointing a pointer to the most downstream feeder section in the low voltage distribution network;

[0013] S1.2: if the feeder section pointed by the pointer is the most downstream feeder section, calculating the operating domain of the power thereof according to its maximum current carrying capacity and the adjustable range of the power of the downstream adjacent commercial load; if the feeder section pointed by the pointer is not the most downstream feeder section, calculating the operating domain of the power thereof according to its maximum current carrying capacity, the adjustable range of the power of the downstream adjacent commercial load and the operating domain of the power of the downstream adjacent feeder section;

[0014] S1.3: judging whether the feeder section pointed by the pointer is the most upstream feeder section in the low voltage distribution network, if yes, marking the task status as completed and turning to step S1.4, otherwise moving the pointer by one feeder section upstream and returning to step S1.2;

[0015] S1.4: taking the operating domain of the power of the most upstream feeder section of the low voltage distribution network as the operating domain of the power flowing from the medium voltage bus into the low voltage distribution network only considering the constraints of the low voltage distribution network.

[0016] Preferably, in step S1.2, if the feeder section pointed by the pointer is the most downstream feeder section, the operating domain of the power of the feeder section pointed by the pointer S Line(LV) is:

[0017]

[0018] wherein, Sdownstreamadjacentload denotes the adjustable range of the power of the downstream adjacent load of the feeder segment pointed by the pointer, and denote the upper limit and the lower limit of the adjustable range S D,Load , respectively; Smax denotes the maximum load flow of the feeder segment pointed by the pointer; min(a, b) denotes the minimum of a and b;

[0019] If the feeder segment pointed by the pointer is not the most downstream feeder segment, the operating range S Line(LV) of the power of the feeder segment pointed by the pointer is:

[0020]

[0021] wherein, Sdownstreamadjacentfeeder denotes the operating range of the power of the downstream adjacent feeder segment of the feeder segment pointed by the pointer, and denote the upper limit and the lower limit of the operating range S D,Line(LV) , respectively.

[0022] Preferably, the possible topologies of the medium voltage distribution network are the topologies of the distribution network formed by reconfiguration after any one of the set of expected faults occurs.

[0023] The step S2 further comprises:

[0024] S2.1: pointing the pointer to the most downstream medium voltage feeder segment in the current topology;

[0025] S2.2: if the feeder segment pointed by the pointer is the most downstream feeder segment, calculating the operating range of the power of the feeder segment pointed by the pointer according to its maximum load flow and the operating range of the power of the downstream adjacent medium voltage bus feeding the connected low voltage distribution network; if the feeder segment pointed by the pointer is not the most downstream feeder segment, calculating the operating range of the power of the feeder segment pointed by the pointer according to its maximum load flow, the operating range of the power of the downstream adjacent medium voltage bus feeding the connected low voltage distribution network and the operating range of the power of the downstream adjacent feeder segment;

[0026] S2.3: if the feeder segment pointed by the pointer is the most upstream feeder segment in the current topology, marking the task state as completed and ending the flow; otherwise, moving the pointer one feeder segment upstream in the current topology and going to step S2.2.

[0027] Preferably, in step S2.2, if the feeder segment pointed by the pointer is the most downstream feeder segment, the operating range S Line(MV) of the power of the feeder segment pointed by the pointer is:

[0028]

[0029] In the formula, This indicates the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar of the feeder segment pointed to by the pointer. and Representing the runtime domain S respectively D,Bus The upper and lower limits; This indicates the maximum current carrying capacity of the feeder segment pointed to by the pointer;

[0030] If the feeder segment pointed to by the pointer is not the downstream feeder segment, the operating domain S of the power of the feeder segment pointed to by the pointer is... Line(MV) for:

[0031]

[0032] in, This indicates the operating domain of the power of the downstream adjacent feeder segment to which the pointer points. and Representing the runtime domain S respectively D,Line(MV) The upper and lower limits.

[0033] Step S3 further includes:

[0034] S3.1: Point the pointer to the upstream medium-voltage feeder segment in the current topology;

[0035] S3.2: If the feeder segment pointed to by the pointer is the most downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar is determined according to its power operating domain; if the feeder segment pointed to by the pointer is not the most downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar and the power of its downstream adjacent feeder segment are determined according to its power operating domain.

[0036] S3.3: Determine whether the feeder segment pointed to by the pointer is the downstream feeder segment in the current topology. If so, mark the task status as completed and end the process; otherwise, move the pointer one feeder segment downstream in the current topology and return to step S3.2.

[0037] Preferably, in step S3.2, if the feeder segment pointed to by the pointer is the most downstream feeder segment, the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage bus is:

[0038]

[0039] In the formula, S U,Bus This represents a vector consisting of the power flowing into the connected low-voltage distribution network from each medium-voltage busbar upstream of the feeder segment pointed to by the pointer; S Line(MV) This indicates the power of the feeder segment pointed to by the pointer; This indicates the operating domain of the power of the feeder segment pointed to by the pointer, where i represents the topology number and f... i (S Line(MV) ,S U,Bus ) indicates that S Line(MV) and S U,Bus A function of variables and Its upper and lower limits; S D,Bus This indicates the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar of the feeder segment pointed to by the pointer;

[0040] If the feeder segment pointed to by the pointer is not the downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar and the power of its downstream adjacent feeder segment is:

[0041]

[0042] In the formula, S D,Line(MV) This indicates the power of the adjacent feeder segment downstream of the feeder segment pointed to by the pointer, f. i (S D,Line(MV) +S D,Bus ,S U,Bus ) indicates that f i (S Line(MV) ,S U,Bus S in ) Line(MV) Replace all with S D,Line(MV) +S D,Bus The function that is formed afterward.

[0043] Step S4 further includes:

[0044] Connect the operating domains of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in all topologies:

[0045]

[0046] Where i represents the topology number, n represents the number of all possible topologies for the medium-voltage distribution network, and S Bus This represents the power variable flowing into the connected low-voltage distribution network from each medium-voltage bus. g represents the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in topology i. i (S Bus ) represents S Bus A linear combination of the variables in the equation. and They represent g respectively i (S Bus The upper and lower limits of );

[0047] The simplex method is used to eliminate redundant constraints, thus obtaining the safety domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus:

[0048] h min ≤h(S Bus )≤h max

[0049] Wherein, h(S) Bus ) is S Bus A linear combination of the variables in h max and h min For h(S) Bus The upper and lower limits of ).

[0050] Furthermore, in step S5, the safety domain h of the power flowing into the connected low-voltage distribution network from each medium-voltage bus is defined. min ≤h(S Bus )≤h max In the middle, S Bus =[S1,S2,…,S m Replace with:

[0051]

[0052] Obtain the security domain of the commercial load cluster in the distribution network:

[0053] h min ≤h(S load )≤h max ;

[0054] Among them, S j Let S be the power variable flowing from medium-voltage bus j into the connected low-voltage distribution network, m be the number of medium-voltage buses, and S be the power variable. j,k Let l be the power variable of the kth commercial load in the low-voltage distribution network connected to the medium-voltage bus j. j The number of commercial loads in the low-voltage distribution network connected to the medium-voltage bus j.

[0055] The beneficial effects of this invention are as follows:

[0056] This invention addresses the problem that current point-by-point simulation-based distribution network security domain calculation methods struggle to balance security domain granularity and computational speed. It proposes a parallel distribution network security domain calculation method that can quickly perform fine-grained distribution network security domain calculations from a load perspective, enabling the definition of the security domain for commercial load clusters in the distribution network and supporting their participation in demand-side response. Attached Figure Description

[0057] Figure 1 This is a flowchart of the parallel computing method for the security domain of a commercial load cluster in a distribution network according to an embodiment of the present invention;

[0058] Figure 2This is a schematic diagram of a distribution network example used in an embodiment of the present invention to illustrate the parallel computing method for the security domain of a commercial load cluster in a distribution network.

[0059] Figure 3 This is a schematic diagram of all possible topologies of a medium-voltage distribution network according to an embodiment of the present invention. Detailed Implementation

[0060] The following embodiments are provided to enable those skilled in the art to more fully understand the present invention, but do not limit the invention in any way.

[0061] like Figure 1 As shown, the parallel computing method for the security domain of a commercial load cluster in a distribution network according to the present invention includes the following steps:

[0062] S1: Calculations are performed in parallel in each low-voltage distribution network. In each low-voltage distribution network, the operating domain of the power of each low-voltage feeder segment is calculated in order from downstream to upstream, based on the adjustable range of the power of each commercial load and the maximum current carrying capacity of each feeder segment. This yields the operating domain of the power flowing into the low-voltage distribution network from the medium-voltage bus when only considering the constraints of the low-voltage distribution network.

[0063] S2: Perform calculations in parallel across all possible topologies of the medium-voltage distribution network. In each possible topology of the medium-voltage distribution network, based on the operating domain of the power flowing from each medium-voltage bus into the connected low-voltage distribution network and the maximum current carrying capacity of each feeder segment, calculate the operating domain of the power of each medium-voltage feeder segment in the topology from downstream to upstream.

[0064] S3: Perform calculations in parallel across all possible topologies of the medium-voltage distribution network. In each possible topology of the medium-voltage distribution network, based on the operating domain of the power of each medium-voltage feeder segment, determine the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in the topology in the order from upstream to downstream.

[0065] S4: Take the overlapping part of the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in all possible topologies of the medium-voltage distribution network, and use it as the safety domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus.

[0066] S5: Replace the power variable flowing into the connected low-voltage distribution network from each medium-voltage bus in the security domain obtained in step S4 with the sum of the power variables of all commercial loads in the low-voltage distribution network connected to the medium-voltage bus, thereby refining the security domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus into the security domain of the commercial load cluster of the distribution network.

[0067] In step S1, within each low-voltage distribution network, the operating domain of the power of each low-voltage feeder segment is calculated in order from downstream to upstream, thereby obtaining the operating domain of the power flowing into the low-voltage distribution network from the medium-voltage bus when only considering the constraints of the low-voltage distribution network. Specifically, this includes:

[0068] S1.1 directs the pointer to the downstream feeder segment in this low-voltage distribution network;

[0069] S1.2 If the feeder segment pointed to by the pointer is the most downstream feeder segment, then its power operating range is calculated based on its maximum current carrying capacity and the adjustable range of the power of its downstream adjacent commercial load; if the feeder segment pointed to by the pointer is not the most downstream feeder segment, then its power operating range is calculated based on its maximum current carrying capacity, the adjustable range of the power of its downstream adjacent commercial load and the power operating range of its downstream adjacent feeder segment.

[0070] S1.3 If the feeder segment pointed to by the pointer is the upstream feeder segment in the low-voltage distribution network, then mark the "task status" as "complete"; otherwise, move the pointer upstream by one feeder segment.

[0071] S1.4 Repeat S1.2 and S1.3 until the "Task Status" is "Completed";

[0072] S1.5 The operating domain of the power of the upstream feeder segment of the low-voltage distribution network shall be taken as the operating domain of the power flowing into the low-voltage distribution network from the medium-voltage bus when only considering the constraints of the low-voltage distribution network.

[0073] Step S1.2 specifically includes:

[0074] If the feeder segment pointed to by the pointer is the most downstream feeder segment, let the adjustable range of the downstream adjacent commercial load power of the feeder segment pointed to by the pointer be... The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer Line(LV) The operating domain is

[0075]

[0076] If the feeder segment pointed to by the pointer is not the downstream feeder segment, let the adjustable range of the downstream adjacent commercial load power of the feeder segment pointed to by the pointer be... The operating domain of the power of the downstream adjacent feeder segment pointed to by the pointer is The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer Line(LV) The operating domain is

[0077]

[0078] in, and It is S D,Load The upper and lower limits of the adjustable range, and It is S D,Line(LV)The upper and lower limits of the operating domain are given by min(a,b), which represents the minimum value between a and b.

[0079] In steps S2 and S3, each possible topology of the medium-voltage distribution network refers to the distribution network topology formed by reconstruction after any one of the anticipated faults in the medium-voltage distribution network occurs.

[0080] In step S2, in each possible topology of the medium-voltage distribution network, the operating domain of the power of each medium-voltage feeder segment in that topology is calculated in order from downstream to upstream, specifically including:

[0081] S2.1 points the pointer to the downstream medium-voltage feeder segment in the current topology;

[0082] S2.2 If the feeder segment pointed to by the pointer is the most downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage bus is determined according to its power operating domain; if the feeder segment pointed to by the pointer is not the most downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage bus and the power of its downstream adjacent feeder segment are determined according to its power operating domain.

[0083] S2.3 If the feeder segment pointed to by the pointer is the upstream feeder segment in the current topology, then mark the "task status" as "complete"; otherwise, move the pointer one feeder segment upstream in the current topology.

[0084] S2.4 Repeat S2.2 and S2.3 until the "Task Status" is "Completed".

[0085] Step S2.2 specifically includes:

[0086] Let the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage bus of the feeder segment pointed to by the pointer be... The operating range of power for the downstream adjacent feeder segment (if any) pointed to by the pointer is: The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer Line(MV) The operating domain is

[0087]

[0088] If the feeder segment pointed to by the pointer is the most downstream feeder segment, let the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar of the feeder segment pointed to by the pointer be... The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer Line(MV) The operating domain is

[0089]

[0090] If the feeder segment pointed to by the pointer is not the downstream feeder segment, let the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage bus of the feeder segment pointed to by the pointer be... The operating domain of the power of the downstream adjacent feeder segment pointed to by the pointer is The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer Line(MV) The operating domain is

[0091]

[0092] in, and For S D,Bus The upper and lower limits of the operating domain, and For S D,Line(MV) The upper and lower limits of the operating domain.

[0093] In step S3, in each possible topology of the medium-voltage distribution network, the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in that topology is determined in order from upstream to downstream. Specifically, this includes:

[0094] S3.1 points the pointer to the upstream medium-voltage feeder segment in the current topology;

[0095] S3.2 If the feeder segment pointed to by the pointer is the most downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar is determined according to its power operating domain; if the feeder segment pointed to by the pointer is not the most downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar and the power of its downstream adjacent feeder segment are determined according to its power operating domain.

[0096] S3.3 If the feeder segment pointed to by the pointer is the downstream feeder segment in the current topology, then mark the "task status" as "complete"; otherwise, move the pointer one feeder segment downstream in the current topology.

[0097] S3.4 Repeat S3.2 and S3.3 until the "Task Status" is "Completed".

[0098] Step S3.2 specifically includes:

[0099] Let S be the vector of power flowing into the connected low-voltage distribution network from each medium-voltage bus upstream of the feeder segment pointed to by the pointer. U,Bus The power of the feeder segment pointed to by the pointer is S. Line(MV) The operating domain of the power of the feeder segment pointed to by the pointer is The operating range of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar and the power of its downstream adjacent feeder segment (if any) is as follows:

[0100]

[0101] Where i is the index of the current topology, f i (S Line(MV) ,S U,Bus ) as S Line(MV) and S U,Bus A function of variables and For its upper and lower limits, S D,Bus S represents the power flowing into the connected low-voltage distribution network from the adjacent medium-voltage bus downstream of the feeder segment indicated by the pointer. D,Line(MV) f represents the power of the feeder segment downstream of the feeder segment pointed to by the pointer. i (S D,Line(MV) +S D,Bus ,S U,Bus ) for f i (S Line(MV) ,S U,Bus S in ) Line(MV) Replace all with S D,Line(MV) +S D,Bus The function that is formed afterward.

[0102] In step S4, the overlapping portion of the operating domains of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in all possible topologies of the medium-voltage distribution network is taken as the safe domain of the power flowing into the connected low-voltage distribution network. Specifically, this includes:

[0103] Connect the operating domains of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in all topologies:

[0104]

[0105] Where i is the topology index, n is the number of all possible topologies for the medium-voltage distribution network, and S Bus The power variables flowing into the connected low-voltage distribution network from each medium-voltage busbar are represented. For the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in topology i calculated in step S3, g i (S Bus ) is S Bus A linear combination of the variables in the equation. and For g i (S Bus The upper and lower limits of ).

[0106] The simplex method is used to eliminate redundant constraints, thus obtaining the safety domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus:

[0107] h min ≤h(SBus )≤h max

[0108] Wherein, h(S) Bus ) is S Bus A linear combination of the variables in h max and h min For h(S) Bus The upper and lower limits of ).

[0109] In step S5, the power variable flowing into the connected low-voltage distribution network from each medium-voltage bus in the safety domain obtained in step S4 is replaced with the sum of the power variables of all commercial loads in the low-voltage distribution network connected to the medium-voltage bus. Specifically, this includes:

[0110] In h min ≤h(S Bus )≤h max In the middle,

[0111] S Bus =[S1,S2,…,S m ]

[0112] Replace with

[0113]

[0114] Obtain the security domain of the commercial load cluster of the distribution network

[0115] h min ≤h(S load )≤h max

[0116] Among them, S j Let S be the power variable flowing from medium-voltage bus j into the connected low-voltage distribution network, m be the number of medium-voltage buses, and S be the power variable. j,k Let l be the power variable of the kth commercial load in the low-voltage distribution network connected to the medium-voltage bus j. j The number of commercial loads in the low-voltage distribution network connected to the medium-voltage bus j. Specific Implementation

[0118] The present invention provides a parallel computing method for the security domain of a commercial load cluster in a distribution network, comprising:

[0119] The first step is to perform calculations in parallel in each low-voltage distribution network. In each low-voltage distribution network, based on the adjustable range of each commercial load power and the maximum current carrying capacity of each feeder segment, the operating domain of the power of each low-voltage feeder segment is calculated in order from downstream to upstream. This yields the operating domain of the power flowing into the low-voltage distribution network from the medium-voltage bus when only considering the constraints of the low-voltage distribution network.

[0120] Figure 2The distribution network shown comprises three low-voltage distribution networks, and calculations are performed simultaneously across these networks. Taking low-voltage distribution network 1 as an example, the pointer is directed to the downstream feeder segment "Commercial Load 1_2 - Commercial Load 1_3". Based on the adjustable power range and maximum current carrying capacity of the adjacent downstream commercial load 1_3, its power operating range is calculated. Let the adjustable power range of commercial load 1_3 be... The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer (1,2),(1,3) The operating domain is

[0121]

[0122] Here, min(a,b) represents the minimum value between a and b.

[0123] Point the pointer to feeder segment "Commercial Load 1_1 - Commercial Load 1_2". Based on the adjustable power range of its downstream adjacent Commercial Load 1_2, the operating power range of its downstream adjacent feeder segment "Commercial Load 1_2 - Commercial Load 1_3", and its maximum current carrying capacity, calculate its power operating range. Let the adjustable power range of Commercial Load 1_2 be... The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer (1,1),(1,2) The operating domain is

[0124]

[0125] Point the pointer to feeder segment "Medium Voltage Busbar 1 - Commercial Load 1_1". Based on the adjustable power range of its downstream adjacent commercial load 1_1, the operating power range of its downstream adjacent feeder segment "Commercial Load 1_1 - Commercial Load 1_2", and its maximum current carrying capacity, calculate its power operating range. Let the adjustable power range of commercial load 1_1 be... The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer (1,0),(1,1) The operating domain is

[0126]

[0127] At this point, the pointer has pointed to the upstream feeder segment, and the calculation task in low-voltage distribution network 1 is complete. The power S of this feeder segment is... (1,0),(1,1) The operating domain is the operating domain of the power flowing into low-voltage distribution network 1 from medium-voltage bus 1, considering only the constraints of the low-voltage distribution network. During the above calculation, similar calculations are also performed simultaneously in low-voltage distribution networks 2 and 3 to obtain the operating domains of the power flowing into low-voltage distribution network 2 from medium-voltage bus 2 and the power flowing into low-voltage distribution network 3 from medium-voltage bus 3, considering only the constraints of the low-voltage distribution network.

[0128] The second step is to perform calculations in parallel across all possible topologies of the medium-voltage distribution network. In each possible topology of the medium-voltage distribution network, based on the operating domain of the power flowing from each medium-voltage bus into the connected low-voltage distribution network and the maximum current carrying capacity of each feeder segment, the operating domain of the power of each medium-voltage feeder segment is calculated in order from downstream to upstream.

[0129] exist Figure 2 In the distribution network shown, there are four possible topologies for the medium-voltage distribution network, such as... Figure 3 As shown. Calculations are performed simultaneously in these four topologies. Taking topology 1 as an example, the pointer is directed to the downstream feeder segment "Medium-voltage bus 1 - Medium-voltage bus 2". Based on the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage bus 1 and the maximum current carrying capacity of the feeder segment pointed to by the pointer, the operating domain of the power in the feeder segment pointed to by the pointer is calculated. Let the operating domain of the power flowing into the connected low-voltage distribution network from medium-voltage bus 1 be... The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer 1,2(MV) The operating domain is

[0130]

[0131] Point the pointer to feeder segment "Medium-voltage busbar 2 - Medium-voltage busbar 3". Based on the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar 2, the operating domain of the power from the downstream adjacent feeder segment "Medium-voltage busbar 1 - Medium-voltage busbar 2", and the maximum current carrying capacity of the feeder segment pointed to by the pointer, calculate the operating domain of the power of the feeder segment pointed to by the pointer. Let the operating domain of the power flowing into the connected low-voltage distribution network from medium-voltage busbar 2 be... The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer 2,3(MV) The operating domain is

[0132]

[0133] Point the pointer to feeder segment "Medium-voltage busbar 3 - Substation". Based on the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar 3, the operating domain of the power from the downstream adjacent feeder segment "Medium-voltage busbar 2 - Medium-voltage busbar 3", and the maximum current carrying capacity of the feeder segment pointed to by the pointer, calculate the operating domain of the power of the feeder segment pointed to by the pointer. Let the operating domain of the power flowing into the connected low-voltage distribution network from medium-voltage busbar 3 be... The maximum current carrying capacity of the feeder segment pointed to by the pointer is The power S of the feeder segment pointed to by the pointer 3,4(MV) The operating domain is

[0134]

[0135] At this point, the pointer has pointed to the upstream feeder segment, and the calculation task in topology 1 is complete. During the above calculation, similar calculations are also being performed simultaneously in topologies 2, 3, and 4.

[0136] The third step is to perform calculations in parallel across all possible topologies of the medium-voltage distribution network. In each possible topology of the medium-voltage distribution network, based on the operating domain of the power of each medium-voltage feeder segment, the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage busbar is determined in order from upstream to downstream.

[0137] Taking Topology 1 as an example, the pointer is directed to the upstream feeder segment "Medium-voltage Busbar 3 - Substation". Based on the operating domain of the power of the feeder segment indicated by the pointer, the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar 3 and the power of the downstream adjacent feeder segment "Medium-voltage Busbar 2 - Medium-voltage Busbar 3" is determined. Let the operating domain of the feeder segment "Medium-voltage Busbar 3 - Substation" be... The operating range of the power flowing into the low-voltage distribution network from medium-voltage bus 3 and the power of the feeder section "medium-voltage bus 2 - medium-voltage bus 3" is as follows:

[0138]

[0139] Point the pointer to feeder segment "Medium-voltage busbar 2 - Medium-voltage busbar 3". Based on the operating domain of the power of the feeder segment indicated by the pointer, determine the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar 2 and the power of the downstream adjacent feeder segment "Medium-voltage busbar 1 - Medium-voltage busbar 2".

[0140]

[0141] Point the pointer to the feeder segment "Medium Voltage Busbar 1 - Medium Voltage Busbar 2". Based on the operating domain of the power of the feeder segment indicated by the pointer, determine the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium voltage busbar 1.

[0142]

[0143] At this point, the pointer has pointed to the downstream feeder segment, and the calculation task in topology 1 is complete. During the above calculations, similar calculations are also being performed simultaneously in topologies 2, 3, and 4.

[0144] The fourth step is to take the overlapping part of the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in all possible topologies of the medium-voltage distribution network, and use it as the safety domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus.

[0145] From Figure 3 The operating domains of the power flowing into the connected low-voltage distribution network from each of the four topologies shown are combined as follows:

[0146]

[0147] The simplex method is used to eliminate redundant constraints, thus obtaining the safety domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus:

[0148] h min ≤h(S1,S2,S3)≤h max

[0149] The fifth step is to replace the power variable flowing into the connected low-voltage distribution network from each medium-voltage bus in the obtained security domain with the sum of the power variables of all commercial loads in the connected low-voltage distribution network, thereby refining the security domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus into the security domain of the commercial load cluster of the distribution network.

[0150] In h min ≤h(S1,S2,S3)≤h max In the middle, replace S1 with S 1,1 +S 1,2 +S 1,3 S2 is replaced with S 2,1 +S 2,2 +S 2,3 S3 is replaced with S 3,1 +S 3,2 +S 3,3 The security domain of the commercial load cluster of the distribution network is obtained.

[0151] h min ≤h(S 1,1 +S 1,2 +S 1,3 ,S 2,1 +S 2,2 +S 2,3 ,S 3,1 +S 3,2 +S 3,3 )≤h max .

[0152] Based on the analysis of the distribution network topology, this invention decomposes the calculation of the security domain of the commercial load cluster of the distribution network into multiple steps, and then decomposes each step into multiple sub-tasks that can be executed in parallel, thereby realizing the rapid definition of the security domain of the commercial load cluster of the distribution network and supporting the participation of the commercial load cluster of the distribution network in demand-side response.

[0153] The above are merely preferred embodiments of the present invention. The scope of protection of the present invention is not limited to the above embodiments. All technical solutions falling within the scope of the present invention's concept are within the scope of protection of the present invention. It should be noted that for those skilled in the art, any improvements and modifications made without departing from the principles of the present invention should be considered within the scope of protection of the present invention.

Claims

1. A parallel computing method for the security domain of commercial load clusters in a distribution network, characterized in that, The method includes the following steps: S1: Calculations are performed in parallel in each low-voltage distribution network. In each low-voltage distribution network, the operating domain of the power of each low-voltage feeder segment is calculated in order from downstream to upstream, based on the adjustable range of the power of each commercial load and the maximum current carrying capacity of each feeder segment. This yields the operating domain of the power flowing into the low-voltage distribution network from the medium-voltage bus when only considering the constraints of the low-voltage distribution network. S2: Perform calculations in parallel across all possible topologies of the medium-voltage distribution network. In each possible topology of the medium-voltage distribution network, based on the operating domain of the power flowing from each medium-voltage bus into the connected low-voltage distribution network and the maximum current carrying capacity of each feeder segment, calculate the operating domain of the power of each medium-voltage feeder segment in the topology from downstream to upstream. S3: Perform calculations in parallel across all possible topologies of the medium-voltage distribution network. In each possible topology of the medium-voltage distribution network, based on the operating domain of the power of each medium-voltage feeder segment, determine the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in the topology in the order from upstream to downstream. S4: Take the overlapping part of the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in all possible topologies of the medium-voltage distribution network, and use it as the safety domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus. S5: Replace the power variable flowing into the connected low-voltage distribution network from each medium-voltage bus in the security domain obtained in step S4 with the sum of the power variables of all commercial loads in the low-voltage distribution network connected to the medium-voltage bus, and refine the security domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus into the security domain of the commercial load cluster of the distribution network.

2. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 1, characterized in that, Step S1 further includes: S1.1: Point the pointer to the downstream feeder segment in the low-voltage distribution network; S1.2: If the feeder segment pointed to by the pointer is the most downstream feeder segment, then its power operating range is calculated based on its maximum current carrying capacity and the adjustable range of the power of its downstream adjacent commercial load; if the feeder segment pointed to by the pointer is not the most downstream feeder segment, then its power operating range is calculated based on its maximum current carrying capacity, the adjustable range of the power of its downstream adjacent commercial load and the power operating range of its downstream adjacent feeder segment. S1.3: Determine whether the feeder segment pointed to by the pointer is the upstream feeder segment in the low-voltage distribution network. If so, mark the task status as completed and proceed to step S1.

4. Otherwise, move the pointer upstream by one feeder segment and return to step S1.

2. S1.4: The operating domain of the power of the upstream feeder segment of the low-voltage distribution network is taken as the operating domain of the power flowing into the low-voltage distribution network from the medium-voltage bus when only considering the constraints of the low-voltage distribution network.

3. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 2, characterized in that, In step S1.2, if the feeder segment pointed to by the pointer is the most downstream feeder segment, the power S of the feeder segment pointed to by the pointer is... Line(LV) The operating domain is: In the formula, This indicates the adjustable range of power for the downstream adjacent commercial load of the feeder segment pointed to by the pointer. and These represent the adjustable range S respectively. D,Load The upper and lower limits; This indicates the maximum current carrying capacity of the feeder segment pointed to by the pointer; min(a,b) represents the minimum value between a and b; If the feeder segment pointed to by the pointer is not the downstream feeder segment, the operating domain S of the power of the feeder segment pointed to by the pointer is... Line(LV) for: in, This indicates the operating domain of the power of the downstream adjacent feeder segment to which the pointer points. and Representing the runtime domain S respectively D,Line(LV) The upper and lower limits.

4. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 1, characterized in that, Each possible topology of the medium-voltage distribution network refers to the distribution network topology formed by reconfiguration after any one of the anticipated faults in the medium-voltage distribution network occurs.

5. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 1, characterized in that, Step S2 further includes: S2.1: Point the pointer to the downstream medium-voltage feeder segment in the current topology; S2.2: If the feeder segment pointed to by the pointer is the most downstream feeder segment, then its operating range is calculated based on its maximum current carrying capacity and the operating range of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage bus. If the feeder segment pointed to by the pointer is not the most downstream feeder segment, then its operating range is calculated based on its maximum current carrying capacity, the operating range of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage bus, and the operating range of the power from its downstream adjacent feeder segment. S2.3: If the feeder segment pointed to by the pointer is the upstream feeder segment in the current topology, then mark the task status as completed and end the process; otherwise, move the pointer one feeder segment upstream in the current topology and proceed to step S2.

2.

6. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 5, characterized in that, In step S2.2, if the feeder segment pointed to by the pointer is the most downstream feeder segment, the operating domain S of the power of the feeder segment pointed to by the pointer is... Line(MV) for: In the formula, This indicates the operating domain of the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar of the feeder segment pointed to by the pointer. and Representing the runtime domain S respectively D,Bus The upper and lower limits; This indicates the maximum current carrying capacity of the feeder segment pointed to by the pointer; If the feeder segment pointed to by the pointer is not the downstream feeder segment, the operating domain S of the power of the feeder segment pointed to by the pointer is... Line(MV) for: in, This indicates the operating domain of the power of the downstream adjacent feeder segment to which the pointer points. and Representing the runtime domain S respectively D,Line(MV) The upper and lower limits.

7. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 1, characterized in that, Step S3 further includes: S3.1: Point the pointer to the upstream medium-voltage feeder segment in the current topology; S3.2: If the feeder segment pointed to by the pointer is the most downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar is determined according to its power operating domain; if the feeder segment pointed to by the pointer is not the most downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar and the power of its downstream adjacent feeder segment are determined according to its power operating domain. S3.3: Determine whether the feeder segment pointed to by the pointer is the downstream feeder segment in the current topology. If so, mark the task status as completed and end the process; otherwise, move the pointer one feeder segment downstream in the current topology and return to step S3.

2.

8. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 7, characterized in that, In step S3.2, if the feeder segment pointed to by the pointer is the most downstream feeder segment, the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage bus is: f i min ≤f i (S D,Bus ,S U,Bus )≤f i max ; In the formula, S U,Bus This represents a vector consisting of the power flowing into the connected low-voltage distribution network from each medium-voltage busbar upstream of the feeder segment pointed to by the pointer; S Line(MV) Indicates the power of the feeder segment pointed to by the pointer; f i min ≤f i (S Line(MV) ,S U,Bus )≤f i max This indicates the operating domain of the power of the feeder segment pointed to by the pointer, where i represents the topology number and f... i (S Line(MV) ,S U,Bus ) indicates that S Line(MV) and S U,Bus f is a function of variables. i max and f i min Its upper and lower limits; S D,Bus This indicates the power flowing into the connected low-voltage distribution network from the downstream adjacent medium-voltage busbar of the feeder segment pointed to by the pointer; If the feeder segment pointed to by the pointer is not the downstream feeder segment, then the operating domain of the power flowing into the connected low-voltage distribution network from its downstream adjacent medium-voltage busbar and the power of its downstream adjacent feeder segment is: f i min ≤f i (S D,Line(MV) +S D,Bus ,S U,Bus )≤f i max ; In the formula, S D,Line(MV) This indicates the power of the adjacent feeder segment downstream of the feeder segment pointed to by the pointer, f. i (S D,Line(MV) +S D,Bus ,S U,Bus ) indicates that f i (S Line(MV) ,S U,Bus S in ) Line(MV) Replace all with S D,Line(MV) +S D,Bus The function that is formed afterward.

9. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 1, characterized in that, Step S4 further includes: Connect the operating domains of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in all topologies: Where i represents the topology number, n represents the number of all possible topologies for the medium-voltage distribution network, and S Bus This represents the power variable flowing into the connected low-voltage distribution network from each medium-voltage bus. g represents the operating domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus in topology i. i (S Bus ) represents S Bus A linear combination of the variables in the equation. and They represent g respectively i (S Bus The upper and lower limits of ); The simplex method is used to eliminate redundant constraints, thus obtaining the safety domain of the power flowing into the connected low-voltage distribution network from each medium-voltage bus: h min ≤h(S Bus )≤h max Wherein, h(S) Bus ) is S Bus A linear combination of the variables in h max and h min For h(S) Bus The upper and lower limits of ).

10. The parallel computing method for the security domain of commercial load clusters in a distribution network according to claim 1, characterized in that, In step S5, the safety domain h of the power flowing into the connected low-voltage distribution network from each medium-voltage bus is determined. min ≤h(S Bus )≤h max In the middle, S Bus =[S1,S2,…,S m Replace with: Obtain the security domain of the commercial load cluster in the distribution network: h min ≤h(S load )≤h max ; Among them, S j Let S be the power variable flowing from medium-voltage bus j into the connected low-voltage distribution network, m be the number of medium-voltage buses, and S be the power variable. j,k Let l be the power variable of the kth commercial load in the low-voltage distribution network connected to the medium-voltage bus j. j The number of commercial loads in the low-voltage distribution network connected to the medium-voltage bus j.

Images