Power distribution network operation envelope determination method and device, storage medium and computer device

Abstract

The application discloses a power distribution network operation envelope determination method and device, a storage medium and computer equipment, and comprises the following steps: constructing a power distribution network linearization safety constraint model based on power grid attribute information; determining a high-dimensional feasible region of a target power distribution network based on the power distribution network linearization safety constraint model; projecting the high-dimensional feasible region to an active-reactive power plane of each distributed power source access node in the power distribution network to obtain a two-dimensional operation envelope, determining a Cartesian product space of multi-node power combinations based on the two-dimensional operation envelope of the node; determining an extreme power vertex based on each power vertex in the two-dimensional operation envelope to obtain a power test point set, performing operation constraint verification on the Cartesian product space based on the power test point set, and determining a safe power point based on the verification result; determining a translation center and a scaling coefficient, translating and scaling the two-dimensional operation envelope of the node based on the translation center and the scaling coefficient, and determining a final operation envelope of the power distribution network based on the translation and scaling results.

Description

Technical Field

[0001] This invention relates to the field of power distribution network safety technology, and in particular to a method, apparatus, storage medium, and computer equipment for determining the operating envelope of a power distribution network. Background Technology

[0002] With the large-scale integration of distributed photovoltaic and other new energy sources into the distribution network, the randomness and volatility of their output pose a severe challenge to the safe operation of the distribution network. As a novel control method, the operating envelope maps the global security constraints of the distribution network to the active-reactive (PQ) power security boundary of each distributed energy access node, allowing each node to autonomously adjust its output within the boundary, achieving local autonomous operation, and effectively reducing the communication overhead and computational pressure of centralized dispatching.

[0003] Currently, the operating envelope of the distribution network is usually determined based on human experience. However, this manual method is time-consuming and labor-intensive, and due to the varying skill levels of the staff, errors in the operating envelope determination can occur. Summary of the Invention

[0004] This invention provides a method, apparatus, storage medium, and computer equipment for determining the operating envelope of a power distribution network, which mainly improves the efficiency and accuracy of determining the operating envelope of a power distribution network.

[0005] According to a first aspect of the present invention, a method for determining the operating envelope of a distribution network is provided, comprising: Obtain the power grid attribute information of the target distribution network, and construct a linearized security constraint model of the distribution network based on the power grid attribute information, wherein the power grid attribute information includes topology, line parameters and distributed power generation equipment capacity; The high-dimensional feasible region of the target distribution network is determined based on the linearized security constraint model of the distribution network. The high-dimensional feasible region is projected onto the active-reactive power plane of each distributed power access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power access node, and the Cartesian product space of the power combination of multiple nodes is determined based on the two-dimensional operating envelope of each distributed power access node. Based on each power vertex in the two-dimensional operating envelope, extreme power vertices are determined, and each extreme power vertex constitutes a power test point set. Based on the power test point set, global safe operation constraints are verified on the Cartesian product space. Based on the global safe operation constraints verification results, safe power points of distributed power access nodes that satisfy global safe operation constraints are determined in the Cartesian product space. The translation center and scaling factor of each distributed power source access node are dynamically determined. Based on the scaling factor and the translation center, the safe power point of the corresponding distributed power source access node is translated and scaled. Based on the translation and scaling results, the final operating envelope of the target distribution network is determined.

[0006] Optionally, a linearized security constraint model for the distribution network is constructed based on the power grid attribute information, including: Based on the power grid attribute information, an augmented impedance matrix is ​​determined, and the node-branch admittance correlation matrix of the target distribution network is determined. Based on the augmented impedance matrix and the admittance correlation matrix, a linear mapping relationship between the square of the node voltage and power injection at the distributed generation access node is determined. ,in, The square vector of node voltages For the augmented impedance matrix (unit: ), Inject vectors into the active power of each distributed power source access node. Injecting reactive power vectors into each distributed power source access node, For admittance correlation matrix, It is a constant vector; Determine the node power constraints for each distributed power source access node: , ,in, and The first Active power injection and reactive power injection of each distributed power source access node and The first The lower and upper limits of active power injection for each distributed power source access node. and The first The lower and upper limits of reactive power injection for each distributed power source access node; Determine the node voltage constraints for each distributed power source access node: ,in, and These are the lower and upper safety limits for node voltages, respectively. For the first The squared value of the voltage of each distributed power supply access node; Determine the root node power constraints of the target distribution network: , , ,in, and These are the active and reactive power of the root node of the target distribution network, respectively. and These represent the total active power and total reactive power of the target distribution network's fixed loads, respectively. The rated apparent power of the root node transformer of the target distribution network is m, and the total number of nodes in the target distribution network where distributed power sources are connected is m. Determine the line power flow constraints of the target distribution network: ,in, and The first in the target distribution network are respectively Line 1 Active power and reactive power flow of the segment For the first Line 1 Segmentation angle, For the first Thermal stability capacity of the line; The linearized security constraint model of the distribution network is determined based on the linear mapping relationship between the square of the node voltage and the power injection, the node power constraint, the node voltage constraint, the root node power constraint, and the line power flow constraint.

[0007] Optionally, determining the high-dimensional feasible region of the target distribution network based on the linearized security constraint model of the distribution network includes: Substituting the node voltage constraint into the linear mapping relationship between the node voltage square and power injection yields a linear inequality constraint for the node active power and node reactive power of each distributed power access node. Based on the node power constraints, the root node power constraints, and the line power flow constraints, a constraint coefficient matrix is ​​determined. Based on the constraint coefficient matrix and the linear inequality constraints, the high-dimensional feasible region of the target distribution network is determined. ,in, For high-dimensional feasible domains, For the constraint coefficient matrix, The decision variable vector of active and reactive power of the target distribution network. To constrain the vector on the right.

[0008] Optionally, determining the Cartesian product space of the power combination of multiple nodes based on the two-dimensional running envelope of each of the distributed power access nodes includes: For each distributed power access node's two-dimensional operating envelope, a power vertex to be added is randomly determined, and the distance of the power vertex to be added relative to the surface point of the two-dimensional operating envelope is determined. The center distance between the surface point of the two-dimensional operating envelope and the center point of the two-dimensional operating envelope is determined. Based on the distance of ... For each new two-dimensional running envelope, a new power vertex to be added is randomly determined, and a new overshoot distance of the new power vertex to be added relative to the surface point of the new two-dimensional running envelope is determined. A new center distance between the surface point of the new two-dimensional running envelope and the center point of the new two-dimensional running envelope is determined. Based on the new overshoot distance and the new center distance, a new increment rate of the new power vertex to be added is determined. If the new increment rate is greater than a preset increment rate threshold, the new power vertex to be added is added as another new vertex to the new vertex set of the new two-dimensional running envelope to obtain another new vertex set. Based on the other new vertex set, another new two-dimensional running envelope is determined. In this way, new vertices are continuously added to the two-dimensional running envelope until the vertex addition termination condition is met, thus obtaining the extended two-dimensional running envelope corresponding to each two-dimensional running envelope. Based on the extended two-dimensional operating envelope, a power combination set for the envelope of each distributed power access node is determined, and the Cartesian product space is determined based on the power combination set.

[0009] Optionally, determining the extreme power vertex based on each power vertex in the two-dimensional running envelope includes: The constraint coefficient matrix is ​​decomposed according to the distributed power access node to obtain the power contribution function of each distributed power access node for each constraint: ,in, , These are the constraint coefficient matrices. The Middle The row corresponds to the first The active and reactive power coefficients of each distributed power supply access node. For the first The distributed power access node connects to the first... Power contribution of the constraints; Substitute each power vertex in the two-dimensional running envelope into the power contribution function to obtain the power contribution of the corresponding power vertex, and determine the extreme power vertex among each power vertex based on the power contribution.

[0010] Optionally, a global safety operation constraint verification is performed on the Cartesian product space based on the power test point set, and a safe power point of the distributed power access node that satisfies the global safety operation constraint is determined in the Cartesian product space based on the global safety operation constraint verification result, including: For each power test point in the power test point set, determine whether the power test point meets the global safe operation constraints: ,in, The total number of power test points in the power test point set. For the first The combination of active and reactive power at each power test point; If the power test point satisfies the global safe operation constraint, then the power test point is determined to be a safe power point in the Cartesian product space.

[0011] Optionally, the translation center and scaling factor of each distributed power access node are dynamically determined, including: Determine the scaling factor constraints for each distributed power access node: ,in, For the first The weighting coefficient of each distributed power access node, For the first Scaling factor of each distributed power access node This represents the total number of nodes connected to the distributed power source. Determine the constraint on the sum of the two-dimensional running envelope areas of each distributed power access node after translation and scaling: ,in, For the first The two-dimensional operating envelope area of ​​a distributed power access node; The spatial constraints of the Cartesian product space determined by the two-dimensional running envelope of each distributed power access node after translation and scaling are as follows: ,in, The power point in the Cartesian product space after translation and scaling; The translation center and scaling factor of each distributed power access node are determined based on the scaling factor constraint, the sum of the two-dimensional running envelope areas constraint, and the spatial constraints of the Cartesian product space.

[0012] According to a second aspect of the present invention, a distribution network operation envelope determination device is provided, comprising: A construction unit is used to acquire the power grid attribute information of the target distribution network and construct a linearized security constraint model of the distribution network based on the power grid attribute information. The power grid attribute information includes topology, line parameters and distributed power generation equipment capacity. A determining unit is used to determine the high-dimensional feasible region of the target distribution network based on the linearized security constraint model of the distribution network. The projection unit is used to project the high-dimensional feasible domain onto the active-reactive power plane of each distributed power access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power access node, and to determine the Cartesian product space of the power combination of multiple nodes based on the two-dimensional operating envelope of each distributed power access node. The constraint verification unit is used to determine the extreme power vertex based on each power vertex in the two-dimensional running envelope, form a power test point set by each extreme power vertex, perform global safe operation constraint verification on the Cartesian product space based on the power test point set, and determine the safe power point of the distributed power access node that satisfies the global safe operation constraint in the Cartesian product space based on the global safe operation constraint verification result. The scaling and translation unit is used to dynamically determine the translation center and scaling factor of each distributed power access node, and to translate and scale the safe power point of the corresponding distributed power access node based on the scaling factor and the translation center, and to determine the final operating envelope of the target distribution network based on the translation and scaling results.

[0013] According to a third aspect of the present invention, a computer-readable storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the above-described method for determining the operating envelope of a power distribution network.

[0014] According to a fourth aspect of the present invention, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the above-described method for determining the operating envelope of a power distribution network.

[0015] According to the present invention, a method, apparatus, storage medium, and computer equipment for determining the operating envelope of a distribution network are provided. Compared with the current method of determining the operating envelope of a distribution network through manual experience, the present invention obtains the topology, line parameters, and distributed power generation capacity of the distribution network, constructs a linearized security constraint model, and determines a high-dimensional feasible region. By integrating the physical properties of the power grid, the nonlinear equations describing the safe operation of the power grid are transformed into a linear convex polyhedral feasible region. This transforms a complex physical problem that is difficult to handle directly into a linear problem that can be solved efficiently using mathematical programming, reducing the complexity of subsequent calculations. By projecting the high-dimensional feasible region onto the active-reactive power plane of each node, a two-dimensional operating envelope is obtained. Based on this, a Cartesian product space of multi-node power combinations is constructed, decomposing the complex high-dimensional security region, in which all nodes are interconnected, into each distributed power generation node. The independent two-dimensional operating range, constructing a Cartesian product space, exhaustively enumerates all possible node combination states. This enables comprehensive analysis of the security of single nodes and combined nodes during the determination of the distribution network operating envelope, thereby improving the accuracy of the distribution network operating envelope determination. By selecting extreme vertices for safety operation constraint verification, the time and resources consumed by verifying all vertices in the operating envelope can be avoided, thus improving verification efficiency and saving verification resources. By dynamically determining the translation center and scaling factor, each distributed power source can have the largest possible power adjustment freedom while ensuring that no limit is exceeded in all extreme cases, avoiding excessively small envelope settings and maximizing the absorption capacity of new energy sources and the flexibility of the power grid. The determination of the distribution network operating envelope in this invention does not require manual intervention, thereby improving the efficiency and accuracy of the distribution network operating envelope determination. Attached Figure Description

[0016] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this application, illustrate exemplary embodiments of the invention and, together with their description, serve to explain the invention and do not constitute an undue limitation thereof. In the drawings: Figure 1 A flowchart of a method for determining the operating envelope of a power distribution network according to an embodiment of the present invention is shown; Figure 2 A flowchart of another method for determining the operating envelope of a power distribution network provided by an embodiment of the present invention is shown; Figure 3 This diagram illustrates the structure of a power distribution network operation envelope determination device according to an embodiment of the present invention. Figure 4 This invention provides a schematic diagram of another distribution network operation envelope determination device according to an embodiment of the present invention. Figure 5 A schematic diagram of the physical structure of a computer device provided in an embodiment of the present invention is shown. Detailed Implementation

[0017] The present invention will be described in detail below with reference to the accompanying drawings and embodiments. It should be noted that, unless otherwise specified, the embodiments and features described in the present application can be combined with each other.

[0018] Currently, the operating envelope of the distribution network is determined based on human experience. However, this manual method is time-consuming and labor-intensive, and due to the varying skill levels of the staff, errors in the operating envelope determination can occur.

[0019] To address the aforementioned problems, embodiments of the present invention provide a method for determining the operating envelope of a distribution network, such as... Figure 1 As shown, the method includes: 101. Obtain the power grid attribute information of the target distribution network, and construct a linearized security constraint model of the distribution network based on the power grid attribute information. The power grid attribute information includes topology, line parameters and distributed power generation equipment capacity.

[0020] The topology includes the node-branch connection relationship of the distribution network, specifying which nodes are power supply nodes, load nodes, or tie nodes; the line parameters include, but are not limited to, the resistance and reactance parameters of the branches of the distribution network and the thermal stability current limit of the lines; the capacity of distributed power generation equipment includes, but is not limited to, the maximum active power output, minimum active power output, maximum reactive power output, minimum reactive power output, and rated apparent power of the distributed power generation access nodes of the distribution network.

[0021] In this embodiment of the invention, the linear mapping relationship between the square of the node voltage and the power injection is determined based on the power grid attribute information. At the same time, the node power constraints, node voltage constraints, root node (the connection point between the distribution network and the transmission network) power constraints and line power flow constraints of each access node in the distribution network are determined. The above constraint information and the linear mapping relationship between the square of the node voltage and the power injection are called the linearized security constraint model of the distribution network.

[0022] 102. Determine the high-dimensional feasible region of the target distribution network based on the linearized security constraint model of the distribution network.

[0023] In this embodiment of the invention, the linear mapping relationship between the square of node voltage and power injection, node power constraints, node voltage constraints, root node power constraints, and line power flow constraints are integrated to obtain a set of thread inequalities with the combination of active and reactive power of nodes as independent variables. Based on this set of thread inequalities, a large number of active and reactive power combinations of distributed power sources connected to nodes that satisfy the set of thread inequalities can be determined. The region formed by the points of these large number of active and reactive power combinations is used as the high-dimensional feasible region of the target distribution network. This embodiment of the invention obtains the topology, line parameters, and distributed power source capacity of the distribution network, constructs a linearized security constraint model, and determines the high-dimensional feasible region. By integrating the physical properties of the power grid, the nonlinear equations describing the safe operation of the power grid are transformed into a linear convex polyhedral feasible region. This transforms a complex physical problem that is difficult to handle directly into a linear problem that can be solved efficiently using mathematical programming, reducing the complexity of subsequent calculations.

[0024] 103. Project the high-dimensional feasible region onto the active-reactive power plane of each distributed power source access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power source access node, and determine the Cartesian product space of the power combination of multiple nodes based on the two-dimensional operating envelope of each distributed power source access node.

[0025] In this embodiment of the invention, for each distributed power source access node i in the target distribution network, a two-dimensional Cartesian coordinate system is established with active power as the horizontal axis and reactive power as the vertical axis, i.e., the active-reactive power plane. A polyhedral projection algorithm is used to eliminate all node variables except node i in the high-dimensional feasible region (high-dimensional convex polyhedron). After projection, a closed polygonal region is obtained on the two-dimensional plane of node i, which is the two-dimensional operational envelope of node i. Further, if the two-dimensional operational envelope of the i-th distributed power source node is a set... Assuming there are N distributed generation nodes in the distribution network, then there are N such sets. This invention constructs a high-dimensional power combination space as the Cartesian product space by performing Cartesian product operations on the two-dimensional operational envelopes of all nodes. This Cartesian product space contains all possible power combinations of all nodes. In this embodiment, the high-dimensional feasible region is projected onto the active-reactive power plane of each node to obtain the two-dimensional operational envelope. Based on this, a Cartesian product space for multi-node power combinations is constructed. This decomposes the complex, interconnected high-dimensional security domain into an independent two-dimensional operational range for each distributed power node. Constructing the Cartesian product space exhaustively enumerates all possible node combination states. This enables comprehensive analysis of the security of single nodes and combined nodes during the determination of the distribution network operational envelope, thereby improving the accuracy of the distribution network operational envelope determination.

[0026] 104. Determine the extreme power vertices based on each power vertex in the two-dimensional operating envelope, and form a power test point set by each extreme power vertex. Perform global safe operation constraint verification on the Cartesian product space based on the power test point set. Based on the global safe operation constraint verification result, determine the safe power point of the distributed power access node that satisfies the global safe operation constraint in the Cartesian product space.

[0027] In this embodiment of the invention, all vertices on the boundary of each two-dimensional operating envelope are identified and extracted. These vertices are called power vertices. Among these vertices, the point with the greatest impact on the distribution network security constraints (such as power constraints, voltage constraints, and line flow constraints) is selected as the extreme power vertex. The extracted extreme power vertices are used to construct a power test point set for global verification. The extreme power vertices of each node are combined. Although the Cartesian product space contains all possible power combinations, in order to reduce the amount of computation, this embodiment of the invention uses extreme power vertices as test points for testing. Each test point represents a multi-node power combination state. It is determined whether each test point in the power test point set satisfies the power constraints, voltage constraints, and line flow constraints of the distribution network. If it does, the test point is regarded as a safe power point. Safe power points are retained in the Cartesian product space, and non-safe power points and their associated non-safe regions are removed. The set of retained power points is the safe power point that satisfies the global safe operation constraints. These points constitute the safety kernel of multi-node collaborative operation, ensuring that the distribution network can remain safe and stable no matter how distributed power sources are combined and operated at these points. The embodiments of the present invention simplify a complex continuous space into the verification of a finite number of discrete points (extreme power vertices), avoiding the resources and time consumed by verifying a large number of vertices.

[0028] 105. Dynamically determine the translation center and scaling factor of each distributed power access node, translate and scale the safe power point of the corresponding distributed power access node based on the translation center and scaling factor, and determine the final operating envelope of the target distribution network based on the translation and scaling results.

[0029] In this embodiment of the invention, the predicted or historical average output value of each distributed power generation access node at the current moment is used as a benchmark, and the optimal center point coordinates are determined as the translation center in conjunction with the current voltage deviation of the power grid. Based on the current stability margin of the power grid, such as voltage stability margin and thermal stability margin, the scaling factor is dynamically calculated. For example, when the power grid is under heavy load or in a vulnerable state, the scaling factor is reduced to compress the adjustable range of the node to ensure absolute safety. When the power grid is under light load or in a robust state, the scaling factor is increased to allow the distributed power generation to have greater adjustment space. Therefore, the translation center and scaling factor can be determined in the above manner. Geometric transformations are then performed on the set of safe power points using the translation center and scaling factor. For example, the safe power points are translated according to the translation center, and the region formed by the translated safe power points is scaled using the scaling factor. A set operation is then performed on all the scaled and translated safe power points of each node. These transformed discrete points are connected using a convex hull algorithm or a polygon fitting algorithm to form a closed two-dimensional geometric region. This region is the final operating envelope of the target distribution network at that moment. The final operating envelope intuitively defines the power range boundary within which all distributed power sources can operate safely and stably under the current dynamic environment. This embodiment of the invention, by dynamically determining the translation center and scaling factor, allows each distributed power source to have the largest possible power adjustment freedom while ensuring that no limit is exceeded under any extreme conditions. This avoids excessively small envelope settings, maximizing the absorption capacity of new energy sources and the flexibility of the power grid. Furthermore, the determination of the distribution network operating envelope in this invention does not require manual intervention, thereby improving the efficiency and accuracy of the determination.

[0030] According to the present invention, a method for determining the operating envelope of a distribution network, compared with the current method of determining the operating envelope of a distribution network through manual experience, the present invention obtains the topology, line parameters, and distributed power generation capacity of the distribution network, constructs a linearized security constraint model, and determines a high-dimensional feasible region. By integrating the physical properties of the power grid, the nonlinear equations describing the safe operation of the power grid are transformed into a linear convex polyhedral feasible region. This transforms a complex physical problem that is difficult to handle directly into a linear problem that can be solved efficiently using mathematical programming, reducing the complexity of subsequent calculations. By projecting the high-dimensional feasible region onto the active-reactive power plane of each node, a two-dimensional operating envelope is obtained. Based on this, a Cartesian product space of multi-node power combinations is constructed, decomposing the complex, interconnected high-dimensional security domain into an independent two-dimensional operating region for each distributed power generation node. The Cartesian product space is constructed to exhaustively enumerate all possible node combinations, enabling comprehensive analysis of the security of single nodes and combined nodes during the determination of the distribution network operating envelope, thereby improving the accuracy of the distribution network operating envelope determination. By selecting extreme vertices for safety operation constraint verification, the time and resources consumed by verifying all vertices in the operating envelope can be avoided, thus improving verification efficiency and saving verification resources. By dynamically determining the translation center and scaling factor, each distributed power source can have the largest possible power adjustment freedom while ensuring that no limit is exceeded in all extreme cases, avoiding excessively small envelope settings and maximizing the absorption capacity of new energy sources and the flexibility of the power grid. The determination of the distribution network operating envelope in this invention does not require manual intervention, thereby improving the efficiency and accuracy of the determination of the distribution network operating envelope.

[0031] Furthermore, to better illustrate the process of determining the operating envelope of the distribution network described above, as a refinement and extension of the above embodiments, this invention provides another method for determining the operating envelope of the distribution network, such as... Figure 2 As shown, the method includes: 201. Obtain the power grid attribute information of the target distribution network, and construct a linearized security constraint model of the distribution network based on the power grid attribute information. The power grid attribute information includes topology, line parameters and distributed power generation equipment capacity.

[0032] In this embodiment of the invention, to improve the efficiency of determining the safe operation envelope of the distribution network, it is first necessary to construct a linearized safety constraint model of the distribution network. Based on this, step 201 specifically includes: determining the augmented impedance matrix based on the power grid attribute information, determining the node-branch admittance correlation matrix of the target distribution network, and determining the linear mapping relationship between the square of the node voltage and power injection of the distributed power access node based on the augmented impedance matrix and the admittance correlation matrix. ,in, The square vector of node voltages For the augmented impedance matrix (unit: ), Inject vectors into the active power of each distributed power source access node. Injecting reactive power vectors into each distributed power source access node, For admittance correlation matrix, Let the vector be a constant; determine the node power constraints for each distributed power source access node: , ,in, and The first Active power injection and reactive power injection of each distributed power source access node and The first The lower and upper limits of active power injection for each distributed power source access node. and The first The lower and upper limits of reactive power injection for each distributed power source access node are determined; the node voltage constraints for each distributed power source access node are also determined. ,in, and These are the lower and upper safety limits for node voltages, respectively. For the first The squared value of the voltage of each distributed power supply access node; Determine the root node power constraints of the target distribution network: , , ,in, and These are the active and reactive power of the root node of the target distribution network, respectively. and These represent the total active power and total reactive power of the target distribution network's fixed loads, respectively. Let m be the rated apparent power of the root node transformer of the target distribution network, and m be the total number of nodes in the target distribution network where distributed generation is connected; determine the line power flow constraints of the target distribution network: ,in, and The first in the target distribution network are respectively Line 1 Active power and reactive power flow of the segment For the first Line 1 Segmentation angle, For the first The thermal stability capacity of each line; the linearized security constraint model of the distribution network is determined based on the linear mapping relationship between the square of the node voltage and the power injection, the node power constraint, the node voltage constraint, the root node power constraint and the line power flow constraint.

[0033] Specifically, a linear power flow model of the distribution network is established using linearized branch power flow equations based on the square of voltage magnitude. The node-branch correlation matrix and admittance matrix are determined based on the grid attribute information of the distribution network. The linear mapping relationship between the square of node voltage and power injection is derived based on the matrix, as shown in the following expression:

[0034] in, The square vector of node voltages (unit: ), dimension ( (Total number of nodes in the distribution network) For the augmented impedance matrix (unit: ), combining network topology and line parameters, with dimensions of , Inject vectors into the active power of each node (unit: (Inflow into the power grid is positive, outflow is negative, dimension is...) , Inject vectors for reactive power at each node (unit: ), dimension ; Admittance Independence Matrix (unit: ), dimension , It is a constant vector, determined by the fixed load and network parameters, with dimension 1. .

[0035] The node power constraints are determined as follows:

[0036]

[0037] in, and The first Active and reactive power injection at each node , The first The lower and upper limits of active power injection at each node. , The first The lower limit and upper limit of reactive power injection at each node. , The number of controllable nodes.

[0038] The node voltage constraints are determined as follows:

[0039] in, , These are the lower and upper safety limits for node voltages, respectively. For the first The squared voltage values ​​of each node are substituted into the formula for the linear mapping relationship between the squared voltage of the conducting node and power injection, transforming it into a function of... , Linear inequality constraints.

[0040] The root node power constraints are determined as follows:

[0041]

[0042]

[0043] in, , These are the active and reactive power of the root node (the connection point between the distribution network and the transmission network). , These represent the total active power and total reactive power of the fixed loads in the distribution network. Rated apparent power of the root node transformer (unit: The power of the root node is constrained to not exceed the device capacity.

[0044] Determine the power flow constraints of the line: A piecewise linearization method is used to approximate the thermal stability constraints of the line. The first of the lines The segmental linear constraints are shown below:

[0045] in, , The first Line 1 The active and reactive currents of the segment, For the first Line 1 Segmentation angle, For the first Thermal stability capacity of each line (unit: ), , This represents the number of segments.

[0046] Finally, the linear mapping relationship between the square of the node voltage and the power injection with the above constraints is used as the linearized security constraint model for the distribution network.

[0047] 202. Determine the high-dimensional feasible region of the target distribution network based on the linearized security constraint model of the distribution network.

[0048] For embodiments of the present invention, in order to determine the feasible region of the distribution network, it is first necessary to determine the high-dimensional feasible region of the distribution network. Based on this, step 202 specifically includes: substituting the node voltage constraints into the linear mapping relationship between the node voltage square and power injection to obtain linear inequality constraints on the node active power and node reactive power of each distributed power access node; determining the constraint coefficient matrix based on the node power constraints, the root node power constraints, and the line power flow constraints; and determining the high-dimensional feasible region of the target distribution network based on the constraint coefficient matrix and the linear inequality constraints. ,in, For high-dimensional feasible domains, For the constraint coefficient matrix, The decision variable vector of active and reactive power of the target distribution network. To constrain the vector on the right.

[0049] Specifically, by integrating the various constraints and the linear mapping relationship between the square of the node voltage and the power injection determined in the preceding steps, the following linear inequality is obtained:

[0050] The above formula can be simplified as follows:

[0051] The region consisting of the points corresponding to the active and reactive power combinations that satisfy this linear inequality is taken as the high-dimensional feasible region of the target distribution network.

[0052] 203. Project the high-dimensional feasible region onto the active-reactive power plane of each distributed power source access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power source access node, and determine the Cartesian product space of the power combination of multiple nodes based on the two-dimensional operating envelope of each distributed power source access node.

[0053] For each embodiment of the present invention, the high-dimensional feasible region is projected onto the active-reactive power plane of the corresponding node to obtain the two-dimensional operating envelope of the node as shown below:

[0054] in, For the two-dimensional running envelope of node i, For the first The active-reactive power combination of each node is projected to preserve convexity, therefore... For a convex polygon, it can be determined through a vertex set. Precise description For the first The number of vertices in an envelope.

[0055] Next, it is necessary to determine the Cartesian product space of the power combination of multiple nodes based on the two-dimensional running envelope of each node. Based on this, step 203 specifically includes: for the two-dimensional running envelope of each distributed power access node, randomly determine the power vertex to be added, and determine the distance of the power vertex to be added relative to the surface point of the two-dimensional running envelope, determine the center distance between the surface point of the two-dimensional running envelope and the center point of the two-dimensional running envelope, determine the increment rate of the power vertex to be added based on the distance and the center distance, if the increment rate is greater than a preset increment rate threshold, then add the power vertex to be added as a new vertex to the vertex set of the two-dimensional running envelope to obtain a new vertex set, and determine a new two-dimensional running envelope based on the new vertex set; for each new two-dimensional running envelope, randomly determine a new power vertex to be added again, and determine the new power vertex to be added... The new overshoot distance of a point relative to the surface point of the new two-dimensional operating envelope is used to determine the new center distance between the surface point of the new two-dimensional operating envelope and the center point of the new two-dimensional operating envelope. Based on the new overshoot distance and the new center distance, the new increment rate of the new power vertex to be added is determined. If the new increment rate is greater than a preset increment rate threshold, the new power vertex to be added is added as another new vertex to the new vertex set of the new two-dimensional operating envelope to obtain another new vertex set. Based on the other new vertex set, another new two-dimensional operating envelope is determined. This process is repeated to continuously add new vertices to the two-dimensional operating envelope until the vertex addition termination condition is met, thus obtaining the extended two-dimensional operating envelope corresponding to each two-dimensional operating envelope. Based on the extended two-dimensional operating envelope, the power combination set of each distributed power access node envelope is determined, and the Cartesian product space is determined based on the power combination set.

[0056] Specifically, an adaptive vertex enumeration algorithm is used to solve the projected envelope. The algorithm achieves accurate vertex enumeration through initialization search, surface dilation, and increment rate verification. The increment rate is determined as follows:

[0057] in, The out-of-range distance of the new search vertex (the power vertex to be added) relative to the 2D running envelope surface. The distance between the two-dimensional running envelope surface and the center point of the envelope. For dimensionless increment rates, a preset increment rate threshold is set according to actual needs. ,when When a new search vertex is identified, it is determined to be a new vertex and added to the vertex set corresponding to the two-dimensional running envelope. This process is repeated to continuously identify new search vertices and check if the increment rate of each new search vertex meets the requirement. If so, it is added to the vertex set until the number of vertices in the vertex set meets the requirement, or until the number of determinations of new search vertices reaches the requirement. Finally, the final determined vertex set constitutes the extended two-dimensional running envelope. Finally, a Cartesian product space containing the power points of each extended two-dimensional running envelope is constructed. As shown below:

[0058] Here, the Cartesian product space is the set of power combinations of all single-node envelopes, and its number of vertices is the product of the number of vertices of each single-node envelope. There is an outer inclusion relationship. , For single-node operation envelope, Run the envelope for the i-th node, where m is the total number of nodes.

[0059] 204. Determine the extreme power vertices based on each power vertex in the two-dimensional operating envelope, and form a power test point set by each extreme power vertex. Perform global safe operation constraint verification on the Cartesian product space based on the power test point set. Based on the global safe operation constraint verification result, determine the safe power point of the distributed power access node that satisfies the global safe operation constraint in the Cartesian product space.

[0060] In this embodiment of the invention, to improve the detection efficiency of power points, it is first necessary to determine the extreme power vertices at each power vertex in the two-dimensional running envelope. Based on this, step 204 specifically includes: decomposing the constraint coefficient matrix according to the distributed power access nodes to obtain the power contribution function of each distributed power access node to each constraint. ,in, , These are the constraint coefficient matrices. The Middle The row corresponds to the first The active and reactive power coefficients of each distributed power supply access node. For the first The distributed power access node connects to the first... The power contribution of each constraint is determined by substituting each power vertex in the two-dimensional running envelope into the power contribution function to obtain the power contribution of the corresponding power vertex, and the extreme power vertex is determined among each power vertex based on the power contribution.

[0061] Specifically, The constraint coefficient matrix A in the matrix is ​​decomposed by nodes, and the first node is defined as follows: constraints ( ) for the first The contribution function of each node is shown below:

[0062] in, , These are the constraint coefficient matrices. No. The corresponding line in the row The coefficients of active and reactive power at each node, For the first The node is the first The power contribution of the constraint reflects the degree to which the nodal power affects the constraint. For the _th Article constraint, if the first The envelope vertices of each node Make the contribution function If the maximum value is reached, then that vertex is the [number]th vertex. The node is the first The extreme vertices (extreme power vertices) of the constraint are denoted as The extreme power vertex satisfies the following formula:

[0063] Wherein, the index of the extreme vertex is such that the first... The node power states most likely to exceed the constraints are identified. Further, an extreme vertex index matrix is ​​constructed. ,element Indicates the first The first constraint The extreme vertex indices corresponding to each node are used to calculate the inner product matrix for fast extreme vertex localization, as shown below:

[0064] in, For the first The constraint corresponds to the first The corresponding node is the 1st The inner product of the vertices is obtained by solving... This yields the indices of the extreme vertices corresponding to each constraint. Then, the extreme vertex combinations are deduplicated to obtain a unique set of test points. Each test point in this test point set Corresponding to a set of extreme vertex combinations, , Let be the i-th extreme vertex.

[0065] Furthermore, after determining the set of test points, it is necessary to perform a global safety operation constraint verification on the Cartesian product space based on the set of test points. Therefore, the method includes: for each power test point in the power test point set, determining whether the power test point satisfies the global safety operation constraints. ,in, The total number of power test points in the power test point set. For the first The combination of active and reactive power at each power test point; if the power test point satisfies the global safe operation constraint, then the power test point is determined to be a safe power point in the Cartesian product space.

[0066] Specifically, verify whether all test points satisfy the following global constraints:

[0067] in, The number of test points after deduplication is much smaller than the number of vertices in the Cartesian product, reducing the verification complexity from exponential to linear. If a test point satisfies the global constraints, then the test point is taken as a safe power point.

[0068] 205. Determine the scaling factor constraint for each distributed power access node, the sum of the two-dimensional running envelope areas of each distributed power access node after translation and scaling, and the spatial constraints of the Cartesian product space.

[0069] Specifically, to maximize the sum of weighted scaling factors while ensuring fairness across all nodes, the scaling factor constraints for anomaly determination are as follows:

[0070] in, For the first The weight coefficients of each node can be set according to the importance of the node, satisfying the following conditions: , For the first Scaling factor of each distributed power access node Let be the total number of distributed power supply access nodes. Maximizing the sum of the operating envelope areas of each node and maximizing the power regulation space, the resulting constraint on the sum of the two-dimensional operating envelope areas of each distributed power supply access node after translation and scaling is shown below:

[0071] in, For the first The two-dimensional operating envelope area of ​​a distributed power access node can be calculated using the following formula: ,

[0072] in, For vertex identifiers, For the first The number of vertices in an envelope.

[0073] The optimization constraints are set so that the Cartesian product after scaling and translation is contained within the original feasible region. The spatial constraints of the Cartesian product space determined by the two-dimensional running envelope of each distributed power access node after translation and scaling, as determined by this requirement, are as follows:

[0074] in, Let represent the active and reactive power of each distributed power access node after translation and scaling. Combining extreme vertex characteristics, it can be simplified as follows:

[0075] in, The test point is the power point in the Cartesian product space after scaling and translation. It is obtained by scaling and translating the extreme vertices as shown below:

[0076]

[0077] in, The active power translation center of extreme vertex i, The scaling factor for extreme vertex i, The reactive power translation center of extreme vertex i The vertex power combination before translation and scaling.

[0078] 206. Based on the scaling factor constraint, the sum of the two-dimensional running envelope area constraint, and the spatial constraint of the Cartesian product space, determine the translation center and scaling factor of each distributed power access node.

[0079] Specifically, based on the scaling factor constraint, the sum of the two-dimensional running envelope areas constraint, and the spatial constraint of the Cartesian product space, a scaling factor that simultaneously satisfies all of the above constraints is determined. Peace Relocation Center .

[0080] 207. Based on the translation center and scaling factor, the safe power point of the corresponding distributed power access node is translated and scaled, and the final operating envelope of the target distribution network is determined based on the translation and scaling results.

[0081] Specifically, after determining the translation center and scaling factor, the safe power point is translated and scaled according to the translation center and scaling factor, and the final operating envelope of the target distribution network is constructed from the translated and scaled safe power point. :

[0082] in, For the safe power point after translation and scaling, For the first Scaling factor for each node , To translate the center, The safe power point before translation and scaling. , The geometric center of node i before translation and scaling. The running envelope of node i before translation and scaling.

[0083] According to another method for determining the operating envelope of a distribution network provided by this invention, compared with the current method of determining the operating envelope of a distribution network through manual experience, this invention obtains the topology, line parameters, and distributed power generation capacity of the distribution network, constructs a linearized security constraint model, and determines a high-dimensional feasible region. By integrating the physical properties of the power grid, the nonlinear equations describing the safe operation of the power grid are transformed into a linear convex polyhedral feasible region. This transforms a complex physical problem that is difficult to handle directly into a linear problem that can be efficiently solved using mathematical programming, reducing the complexity of subsequent calculations. By projecting the high-dimensional feasible region onto the active-reactive power plane of each node, a two-dimensional operating envelope is obtained. Based on this, a Cartesian product space of multi-node power combinations is constructed, decomposing the complex, interconnected high-dimensional security domain into an independent two-dimensional operating envelope for each distributed power generation node. The Cartesian product space, constructed within the range of nodes, exhaustively enumerates all possible node combinations. This enables comprehensive analysis of the security of single and combined nodes during the determination of the distribution network operating envelope, thereby improving the accuracy of the distribution network operating envelope determination. By selecting extreme vertices for safety constraint verification, the time and resources consumed by verifying all vertices in the operating envelope are avoided, thus improving verification efficiency and saving verification resources. Dynamically determining the translation center and scaling factor allows each distributed power source to have the maximum possible power adjustment freedom while ensuring that no limit is exceeded in any extreme case, avoiding excessively small envelope settings and maximizing the absorption capacity of new energy sources and the flexibility of the power grid. Furthermore, the determination of the distribution network operating envelope in this invention does not require manual intervention, thereby improving the efficiency and accuracy of the determination.

[0084] Furthermore, as Figure 1In specific implementation, embodiments of the present invention provide a distribution network operation envelope determination device, such as... Figure 3 As shown, the device includes: a construction unit 31, a determination unit 32, a projection unit 33, a constraint verification unit 34, and a scaling and translation unit 35.

[0085] The construction unit 31 can be used to obtain the power grid attribute information of the target distribution network and construct a linearized security constraint model of the distribution network based on the power grid attribute information. The power grid attribute information includes topology, line parameters and distributed power generation equipment capacity.

[0086] The determining unit 32 can be used to determine the high-dimensional feasible region of the target distribution network based on the linearized security constraint model of the distribution network.

[0087] The projection unit 33 can be used to project the high-dimensional feasible domain onto the active-reactive power plane of each distributed power access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power access node, and determine the Cartesian product space of the power combination of multiple nodes based on the two-dimensional operating envelope of each distributed power access node.

[0088] The constraint verification unit 34 can be used to determine extreme power vertices based on each power vertex in the two-dimensional operating envelope, form a power test point set by each extreme power vertex, perform global safe operation constraint verification on the Cartesian product space based on the power test point set, and determine the safe power point of the distributed power access node that satisfies the global safe operation constraint in the Cartesian product space based on the global safe operation constraint verification result.

[0089] The scaling and translation unit 35 can be used to dynamically determine the translation center and scaling factor of each distributed power access node, translate and scale the safe power point of the corresponding distributed power access node based on the scaling factor and the translation center, and determine the final operating envelope of the target distribution network based on the translation and scaling results.

[0090] In specific application scenarios, in order to construct a linearized security constraint model for the distribution network, the construction unit 31 can be used to determine the augmented impedance matrix based on the power grid attribute information, determine the node-branch admittance correlation matrix of the target distribution network, and determine the linear mapping relationship between the square of the node voltage and power injection of the distributed generation access node based on the augmented impedance matrix and the admittance correlation matrix. ,in, The square vector of node voltages For the augmented impedance matrix (unit: ), Inject vectors into the active power of each distributed power source access node. Injecting reactive power vectors into each distributed power source access node, For admittance correlation matrix, Let the vector be a constant; determine the node power constraints for each distributed power source access node: , ,in, and The first Active power injection and reactive power injection of each distributed power source access node and The first The lower and upper limits of active power injection for each distributed power source access node. and The first The lower and upper limits of reactive power injection for each distributed power source access node are determined; the node voltage constraints for each distributed power source access node are also determined. ,in, and These are the lower and upper safety limits for node voltages, respectively. For the first The squared voltage values ​​of each distributed power source access node; determining the root node power constraints of the target distribution network: , , ,in, and These are the active and reactive power of the root node of the target distribution network, respectively. and These represent the total active power and total reactive power of the target distribution network's fixed loads, respectively. Let m be the rated apparent power of the root node transformer of the target distribution network, and m be the total number of nodes in the target distribution network where distributed generation is connected; determine the line power flow constraints of the target distribution network: ,in, and The first in the target distribution network are respectively Line 1 Active power and reactive power flow of the segment For the first Line 1 Segmentation angle, For the first The thermal stability capacity of each line; the linearized security constraint model of the distribution network is determined based on the linear mapping relationship between the square of the node voltage and the power injection, the node power constraint, the node voltage constraint, the root node power constraint and the line power flow constraint.

[0091] In specific application scenarios, in order to determine the high-dimensional feasible region of the target distribution network, such as Figure 4 As shown, the determining unit 32 includes a power injection module 321 and a feasible region determining module 322.

[0092] The power injection module 321 can be used to substitute the node voltage constraint into the linear mapping relationship between the node voltage square and power injection to obtain a linear inequality constraint on the node active power and node reactive power of each distributed power access node.

[0093] The feasible region determination module 322 can be used to determine a constraint coefficient matrix based on the node power constraints, the root node power constraints, and the line power flow constraints, and to determine the high-dimensional feasible region of the target distribution network based on the constraint coefficient matrix and the linear inequality constraints. ,in, For high-dimensional feasible domains, For the constraint coefficient matrix, The decision variable vector of active and reactive power of the target distribution network. To constrain the vector on the right.

[0094] In specific application scenarios, in order to determine the Cartesian product space, the projection unit 33 includes an expansion module 331 and a space determination module 332.

[0095] The expansion module 331 can be used to randomly determine the power vertex to be added for each distributed power access node's two-dimensional operating envelope, determine the distance of the power vertex to be added relative to the surface point of the two-dimensional operating envelope, determine the center distance between the surface point of the two-dimensional operating envelope and the center point of the two-dimensional operating envelope, determine the increment rate of the power vertex to be added based on the distance and the center distance, and if the increment rate is greater than a preset increment rate threshold, then the power vertex to be added is added as a new vertex to the vertex set of the two-dimensional operating envelope to obtain a new vertex set, and a new two-dimensional operating envelope is determined based on the new vertex set.

[0096] The expansion module 331 can also be used to randomly determine new power vertices to be added for each new two-dimensional running envelope, determine the new overshoot distance of the new power vertices to be added relative to the surface points of the new two-dimensional running envelope, determine the new center distance between the surface points of the new two-dimensional running envelope and the center point of the new two-dimensional running envelope, determine the new increment rate of the new power vertices to be added based on the new overshoot distance and the new center distance, and if the new increment rate is greater than a preset increment rate threshold, then the new power vertices to be added are added as new vertices to the new vertex set of the new two-dimensional running envelope to obtain another new vertex set, and another new two-dimensional running envelope is determined based on the other new vertex set. In this way, new vertices are continuously added to the two-dimensional running envelope until the vertex addition termination condition is met, so as to obtain the expanded two-dimensional running envelope corresponding to each two-dimensional running envelope.

[0097] The space determination module 332 can be used to determine the power combination set of the envelope of each distributed power access node based on the extended two-dimensional running envelope, and to determine the Cartesian product space based on the power combination set.

[0098] In specific application scenarios, in order to determine the extreme power vertex, the constraint verification unit 34 includes a decomposition module 341 and a solution module 342.

[0099] The decomposition module 341 can be used to decompose the constraint coefficient matrix according to the distributed power access node to obtain the power contribution function of each distributed power access node to each constraint: ,in, , These are the constraint coefficient matrices. The Middle The row corresponds to the first The active and reactive power coefficients of each distributed power supply access node. For the first The distributed power access node connects to the first... Power contribution under constraints.

[0100] The solution module 342 can be used to substitute each power vertex in the two-dimensional running envelope into the power contribution function to obtain the power contribution of the corresponding power vertex, and determine the extreme power vertex among each power vertex based on the power contribution.

[0101] In specific application scenarios, in order to perform global safe operation constraint verification on the Cartesian product space, the constraint verification unit 34 further includes: a judgment module 343 and a determination module 344.

[0102] The judgment module 343 can be used to determine whether each power test point in the power test point set meets the global safe operation constraints. ,in, The total number of power test points in the power test point set. For the first The combination of active and reactive power at each power test point.

[0103] The determination module 344 can be used to determine that the power test point is a safe power point in the Cartesian product space if the power test point satisfies the global safe operation constraint.

[0104] In specific application scenarios, in order to determine the translation center and scaling factor of each distributed power access node, the scaling translation unit 35 can be used to determine the scaling factor constraint of each distributed power access node: ,in, For the first The weighting coefficient of each distributed power access node, For the first Scaling factor of each distributed power access node Given the total number of distributed power source access nodes; determine the constraint on the sum of the two-dimensional running envelope areas of each distributed power source access node after translation and scaling: ,in, For the first The two-dimensional operating envelope area of ​​each distributed power access node; the spatial constraints of the Cartesian product space determined by the two-dimensional operating envelope of each distributed power access node after translation and scaling: ,in, For the power points in the Cartesian product space after translation and scaling; determine the translation center and scaling factor of each distributed power access node based on the scaling factor constraint, the sum of the two-dimensional running envelope area constraint, and the spatial constraints of the Cartesian product space.

[0105] It should be noted that other corresponding descriptions of the functional modules involved in the distribution network operation envelope determination device provided in this embodiment of the invention can be found in [reference needed]. Figure 1 The corresponding description of the method shown will not be repeated here.

[0106] Based on the above, Figure 1Correspondingly, this embodiment of the invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, performs the following steps: obtaining grid attribute information of a target distribution network; constructing a linearized security constraint model of the distribution network based on the grid attribute information, wherein the grid attribute information includes topology, line parameters, and distributed power generation equipment capacity; determining the high-dimensional feasible region of the target distribution network based on the linearized security constraint model of the distribution network; projecting the high-dimensional feasible region onto the active-reactive power plane of each distributed power generation access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power generation access node, and based on the two-dimensional operating envelope of each distributed power generation access node... The row envelope determines the Cartesian product space of the power combination of multiple nodes; based on each power vertex in the two-dimensional operating envelope, extreme power vertices are determined, and each extreme power vertex constitutes a power test point set. Based on the power test point set, the global safe operation constraint verification is performed on the Cartesian product space. Based on the global safe operation constraint verification result, the safe power point of the distributed power access node that satisfies the global safe operation constraint is determined in the Cartesian product space; the translation center and scaling factor of each distributed power access node are dynamically determined respectively, and the safe power point of the corresponding distributed power access node is translated and scaled based on the scaling factor and the translation center. Based on the translation and scaling results, the final operating envelope of the target distribution network is determined.

[0107] Based on the above, Figure 1 The method shown and as Figure 3 The embodiment of the device shown in the invention also provides a physical structure diagram of a computer device, such as... Figure 5As shown, the computer device includes: a processor 41, a memory 42, and a computer program stored in the memory 42 and executable on the processor. Both the memory 42 and the processor 41 are mounted on a bus 43. When the processor 41 executes the program, it performs the following steps: acquiring grid attribute information of the target distribution network; constructing a linearized security constraint model of the distribution network based on the grid attribute information, wherein the grid attribute information includes topology, line parameters, and distributed power generation equipment capacity; determining the high-dimensional feasible region of the target distribution network based on the linearized security constraint model; projecting the high-dimensional feasible region onto the active-reactive power plane of each distributed power generation access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power generation access node, and based on... The two-dimensional operating envelope of each distributed power access node determines the Cartesian product space of the power combination of multiple nodes; extreme power vertices are determined based on each power vertex in the two-dimensional operating envelope, and a power test point set is formed by each extreme power vertex. Global safety operation constraint verification is performed on the Cartesian product space based on the power test point set. Based on the global safety operation constraint verification result, the safe power point of the distributed power access node that satisfies the global safety operation constraint is determined in the Cartesian product space; the translation center and scaling factor of each distributed power access node are dynamically determined respectively, and the safe power point of the corresponding distributed power access node is translated and scaled based on the scaling factor and the translation center. The final operating envelope of the target distribution network is determined based on the translation and scaling results.

[0108] Through the technical solution of this invention, the topology, line parameters, and distributed power generation capacity of the distribution network are obtained to construct a linearized security constraint model and determine a high-dimensional feasible region. By integrating the physical properties of the power grid, the nonlinear equations describing the safe operation of the power grid are transformed into a linear convex polyhedral feasible region. This transforms a complex physical problem that is difficult to handle directly into a linear problem that can be efficiently solved using mathematical programming, reducing the complexity of subsequent calculations. By projecting the high-dimensional feasible region onto the active-reactive power plane of each node, a two-dimensional operating envelope is obtained. Based on this, a Cartesian product space of multi-node power combinations is constructed. The complex high-dimensional security region, in which all nodes are interconnected, is decomposed into an independent two-dimensional operating range for each distributed power generation node. The construction of the Cartesian product space exhaustively enumerates all... The possible node combination states enable comprehensive analysis of the security of single nodes and combined nodes during the determination of the distribution network operating envelope, thereby improving the accuracy of the distribution network operating envelope determination. By selecting extreme vertices for safety operation constraint verification, the time and resources consumed by verifying all vertices in the operating envelope can be avoided, thus improving verification efficiency and saving verification resources. By dynamically determining the translation center and scaling factor, each distributed power source can have the largest possible power adjustment freedom while ensuring that no limit is exceeded in all extreme cases, avoiding excessively small envelope settings and maximizing the absorption capacity of new energy sources and the flexibility of the power grid. The determination of the distribution network operating envelope in this invention does not require manual intervention, thereby improving the efficiency and accuracy of the determination of the distribution network operating envelope.

[0109] It is obvious to those skilled in the art that the modules or steps of the present invention described above can be implemented using general-purpose computing devices. They can be centralized on a single computing device or distributed across a network of multiple computing devices. Optionally, they can be implemented using computer-executable program code, thereby storing them in a storage device for execution by a computing device. In some cases, the steps shown or described can be performed in a different order than those presented herein, or they can be fabricated as separate integrated circuit modules, or multiple modules or steps can be fabricated as a single integrated circuit module. Thus, the present invention is not limited to any particular combination of hardware and software.

[0110] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for determining the operating envelope of a power distribution network, characterized in that, include: Obtain the power grid attribute information of the target distribution network, and construct a linearized security constraint model of the distribution network based on the power grid attribute information, wherein the power grid attribute information includes topology, line parameters and distributed power generation equipment capacity; The high-dimensional feasible region of the target distribution network is determined based on the linearized security constraint model of the distribution network. The high-dimensional feasible region is projected onto the active-reactive power plane of each distributed power access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power access node, and the Cartesian product space of the power combination of multiple nodes is determined based on the two-dimensional operating envelope of each distributed power access node. Based on each power vertex in the two-dimensional operating envelope, extreme power vertices are determined, and each extreme power vertex constitutes a power test point set. Based on the power test point set, global safe operation constraints are verified on the Cartesian product space. Based on the global safe operation constraints verification results, safe power points of distributed power access nodes that satisfy global safe operation constraints are determined in the Cartesian product space. The translation center and scaling factor of each distributed power source access node are dynamically determined. Based on the scaling factor and the translation center, the safe power point of the corresponding distributed power source access node is translated and scaled. Based on the translation and scaling results, the final operating envelope of the target distribution network is determined.

2. The method according to claim 1, characterized in that, Based on the aforementioned power grid attribute information, a linearized security constraint model for the distribution network is constructed, including: Based on the power grid attribute information, an augmented impedance matrix is ​​determined, and the node-branch admittance correlation matrix of the target distribution network is determined. Based on the augmented impedance matrix and the admittance correlation matrix, a linear mapping relationship between the square of the node voltage and power injection at the distributed generation access node is determined. ,in, The square vector of node voltages For the augmented impedance matrix (unit: ), Inject vectors into the active power of each distributed power source access node. Injecting reactive power vectors into each distributed power source access node, For admittance correlation matrix, It is a constant vector; Determine the node power constraints for each distributed power supply access node: , ,in, and The first Active power injection and reactive power injection of each distributed power source access node and The first The lower and upper limits of active power injection for each distributed power source access node. and The first The lower and upper limits of reactive power injection for each distributed power source access node; Determine the node voltage constraints for each distributed power source access node: ,in, and These are the lower and upper safety limits for node voltages, respectively. For the first The squared value of the voltage of each distributed power supply access node; Determine the root node power constraints of the target distribution network: , , ,in, and These are the active and reactive power of the root node of the target distribution network, respectively. and These represent the total active power and total reactive power of the target distribution network's fixed loads, respectively. The rated apparent power of the root node transformer of the target distribution network is m, and the total number of nodes in the target distribution network where distributed power sources are connected is m. Determine the line flow of the target distribution network Flow constraints: ,in, and The first in the target distribution network Line 1 Active power and reactive power flow of the segment For the first Line 1 Segmentation angle, For the first Thermal stability capacity of the line; The linearized security constraint model of the distribution network is determined based on the linear mapping relationship between the square of the node voltage and the power injection, the node power constraint, the node voltage constraint, the root node power constraint, and the line power flow constraint.

3. The method according to claim 2, characterized in that, The high-dimensional feasible region of the target distribution network is determined based on the linearized security constraint model of the distribution network, including: Substituting the node voltage constraint into the linear mapping relationship between the node voltage square and power injection yields a linear inequality constraint for the node active power and node reactive power of each distributed power access node. Based on the node power constraints, the root node power constraints, and the line power flow constraints, a constraint coefficient matrix is ​​determined. Based on the constraint coefficient matrix and the linear inequality constraints, the high-dimensional feasible region of the target distribution network is determined. ,in, For high-dimensional feasible domains, For the constraint coefficient matrix, The decision variable vector of active and reactive power of the target distribution network. To constrain the vector on the right.

4. The method according to claim 1, characterized in that, The Cartesian product space of the power combination of multiple nodes is determined based on the two-dimensional running envelope of each of the distributed power access nodes, including: For each distributed power access node's two-dimensional operating envelope, a power vertex to be added is randomly determined, and the distance of the power vertex to be added relative to the surface point of the two-dimensional operating envelope is determined. The center distance between the surface point of the two-dimensional operating envelope and the center point of the two-dimensional operating envelope is determined. Based on the distance of ... For each new two-dimensional running envelope, a new power vertex to be added is randomly determined, and a new overshoot distance of the new power vertex to be added relative to the surface point of the new two-dimensional running envelope is determined. A new center distance between the surface point of the new two-dimensional running envelope and the center point of the new two-dimensional running envelope is determined. Based on the new overshoot distance and the new center distance, a new increment rate of the new power vertex to be added is determined. If the new increment rate is greater than a preset increment rate threshold, the new power vertex to be added is added as another new vertex to the new vertex set of the new two-dimensional running envelope to obtain another new vertex set. Based on the other new vertex set, another new two-dimensional running envelope is determined. In this way, new vertices are continuously added to the two-dimensional running envelope until the vertex addition termination condition is met, thus obtaining the extended two-dimensional running envelope corresponding to each two-dimensional running envelope. Based on the extended two-dimensional operating envelope, a power combination set for the envelope of each distributed power access node is determined, and the Cartesian product space is determined based on the power combination set.

5. The method according to claim 3, characterized in that, Determining extreme power vertices based on each power vertices in the two-dimensional running envelope includes: The constraint coefficient matrix is ​​decomposed according to the distributed power access node to obtain the power contribution function of each distributed power access node for each constraint: ,in, , These are the constraint coefficient matrices. The Middle The row corresponds to the first The active and reactive power coefficients of each distributed power supply access node. For the first The distributed power access node connects to the first... Power contribution of the constraints; Substitute each power vertex in the two-dimensional running envelope into the power contribution function to obtain the power contribution of the corresponding power vertex, and determine the extreme power vertex among each power vertex based on the power contribution.

6. The method according to claim 3, characterized in that, Based on the power test point set, a global safety operation constraint verification is performed on the Cartesian product space. Based on the global safety operation constraint verification result, safe power points of distributed power access nodes that satisfy the global safety operation constraints are determined in the Cartesian product space, including: For each power test point in the power test point set, determine whether the power test point meets the global safe operation constraints: ,in, The total number of power test points in the power test point set. For the first The combination of active and reactive power at each power test point; If the power test point satisfies the global safe operation constraint, then the power test point is determined to be a safe power point in the Cartesian product space.

7. The method according to claim 3, characterized in that, The translation center and scaling factor of each distributed power access node are dynamically determined, including: Determine the scaling factor constraints for each distributed power access node: ,in, For the first The weighting coefficient of each distributed power access node, For the first Scaling factor of each distributed power access node, This represents the total number of nodes connected to the distributed power source. Determine the constraint on the sum of the two-dimensional running envelope areas of each distributed power access node after translation and scaling: ,in, For the first The two-dimensional operating envelope area of ​​a distributed power access node; The spatial constraints of the Cartesian product space determined by the two-dimensional running envelope of each distributed power access node after translation and scaling are as follows: ,in, The power point in the Cartesian product space after translation and scaling; The translation center and scaling factor of each distributed power access node are determined based on the scaling factor constraint, the sum of the two-dimensional running envelope areas constraint, and the spatial constraints of the Cartesian product space.

8. A device for determining the operating envelope of a power distribution network, characterized in that, include: A construction unit is used to acquire the power grid attribute information of the target distribution network and construct a linearized security constraint model of the distribution network based on the power grid attribute information. The power grid attribute information includes topology, line parameters and distributed power generation equipment capacity. A determining unit is used to determine the high-dimensional feasible region of the target distribution network based on the linearized security constraint model of the distribution network. The projection unit is used to project the high-dimensional feasible domain onto the active-reactive power plane of each distributed power access node in the target distribution network to obtain the two-dimensional operating envelope of the corresponding distributed power access node, and to determine the Cartesian product space of the power combination of multiple nodes based on the two-dimensional operating envelope of each distributed power access node. The constraint verification unit is used to determine the extreme power vertex based on each power vertex in the two-dimensional running envelope, form a power test point set by each extreme power vertex, perform global safe operation constraint verification on the Cartesian product space based on the power test point set, and determine the safe power point of the distributed power access node that satisfies the global safe operation constraint in the Cartesian product space based on the global safe operation constraint verification result. The scaling and translation unit is used to dynamically determine the translation center and scaling factor of each distributed power access node, and to translate and scale the safe power point of the corresponding distributed power access node based on the scaling factor and the translation center, and to determine the final operating envelope of the target distribution network based on the translation and scaling results.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

10. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.

Images