CN110350508B - Method for constructing unified fault recovery model of active power distribution network - Google Patents
Method for constructing unified fault recovery model of active power distribution network Download PDFInfo
- Publication number
- CN110350508B CN110350508B CN201910406777.0A CN201910406777A CN110350508B CN 110350508 B CN110350508 B CN 110350508B CN 201910406777 A CN201910406777 A CN 201910406777A CN 110350508 B CN110350508 B CN 110350508B
- Authority
- CN
- China
- Prior art keywords
- node
- load
- representing
- power
- constraint
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 54
- 238000011084 recovery Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000013486 operation strategy Methods 0.000 claims abstract description 18
- 238000013178 mathematical model Methods 0.000 claims abstract description 17
- 238000005192 partition Methods 0.000 claims description 33
- 238000005457 optimization Methods 0.000 claims description 22
- 230000005540 biological transmission Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000002040 relaxant effect Effects 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims 1
- 230000009466 transformation Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 3
- 230000004927 fusion Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 241001672018 Cercomela melanura Species 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/388—Islanding, i.e. disconnection of local power supply from the network
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
- Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
- Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a method for recovering a unified model of a fault of an active power distribution network by simultaneously considering reconstruction and islanding, which comprises the following steps: 1) acquiring power grid fault information and power distribution network information, and summarizing data required by model building; 2) determining an objective function of a fault recovery operation strategy; 3) constructing a unified mathematical model which comprises DistFlow power flow constraint and connectivity constraint using virtual power flow and comprises reconstruction and island division; 4) and solving the established unified mathematical model containing reconstruction and island division to obtain a final active power distribution network fault recovery operation strategy. The power distribution network fault recovery operation strategy simultaneously considers two strategies of network fault reconstruction and island division, considers controllability of load and distributed power supply output, and compared with the traditional power distribution network fault recovery, the fault recovery strategy is enabled to cover more possibilities in an optimized mode, the problem is converted into mixed integer linear programming, strategy generation time is saved, and the optimal solution can be obtained more accurately.
Description
Technical Field
The invention relates to a method for constructing a unified fault recovery model of an active power distribution network, and belongs to the field of operation and planning of the power distribution network.
Background
Nowadays, due to the access of distributed power sources and the popularization of flexible loads, a power distribution network is developed to a more controllable, more variable and more complex network, and students are also prompted to improve the fault recovery operation strategy of the traditional power distribution network. Meanwhile, after a subsequent line is automatically split when a network fault occurs, the peripheral load is driven by using the black start function of the distributed power supply, so that partial load is restored to power supply, namely, island division optimization. In the aspect of a model, active power distribution network fault reconstruction and island division are brought into a unified mathematical model for solving, and compared with the existing research, the integrated model is divided into two layers of operation, so that the optimal solution is easier to obtain.
In terms of the algorithm, for island division, a learner uses a knapsack Tree (TKP) tree algorithm and adds island fusion judgment after calculation is finished, island fusion is carried out after optimization by the method, the optimization result is modified, and the optimization result is determined to be influenced. However, the scholars use the mixed integer second-order cone programming, but set the island number as the number of the distributed power supplies, namely island fusion is not considered, and further optimization of the strategy is also limited. The invention solves the problem by adopting the radial constraint and the connectivity constraint suitable for the island division, establishes a unified model solution by using mixed integer second-order cone programming, is easier to obtain an optimal solution, and can save a power distribution network fault recovery operation strategy.
Disclosure of Invention
The invention aims to solve the problems and provides a method for constructing a fault recovery unified model of an active power distribution network.
In order to achieve the purpose, the method adopted by the invention is as follows: a method for constructing a unified fault recovery model of an active power distribution network comprises the following steps:
step 1) acquiring power grid fault information, a network topology structure, distribution network branches, a distributed power supply and switch state information, predicting load and fault time required by a fault period, and summarizing data required by model building;
step 2) determining a target function of a fault recovery operation strategy, considering the economy of a power grid, considering multi-period load shedding cost, considering the switch operation cost, and considering the load shedding cost and the weight grade of the load;
step 3) constructing a unified mathematical model containing reconstruction and island division, wherein the unified mathematical model contains DistFlow power flow constraint and connectivity constraint using virtual power flow;
and 4) solving the established unified mathematical model containing reconstruction and island division, and obtaining a final active power distribution network fault recovery operation strategy.
Preferably, in step 1), the input fault time period fault information is discretized, the fault time period is divided into n time intervals, and the corresponding predicted load is discretized in the same time scale, so as to prepare for establishing a multi-time-period dynamic unified mathematical model including reconstruction and island division.
Preferably, in step 2), the objective function of the fail-back operation strategy is determined by considering the load shedding cost and the weight level of the load while considering the load shedding cost, with the main objective of minimizing the load shedding amount and considering the switch operation cost. To ensure that the node load shedding amount is minimum, namely the recovery load amount is maximum, the predicted load is set as a constant, the load shedding amount of each node is set as an optimized variable, and meanwhile, the on-off state of a switch is set as a variable from 0 to 1, wherein 0 represents that the switch is off, and 1 represents that the switch is on. The target function considers the whole fault time interval, and the target function expression of the fault recovery operation strategy is as follows:
wherein, CopeRepresents the total cost, cloadRepresenting the load shedding cost factor, NTRepresenting the total number of time intervals of the fault period, NbusRepresenting the total number of nodes in the network, at representing the length of each time interval of the fault period,represents t time period, i node load shedding amount, lambdaiThe load weight of the ith node is obtained; c. CswitchRepresents a switch operation cost coefficient, and is an optimized variable of 0-1, alphapRepresenting the fault recovery of the No. p branchThe switch state of the complex time interval is 0-1 constant,representing the switching state before fault in branch No. p, NlineRepresenting the total number of branches in the network.
Preferably, in step 3), the constraint condition of the established unified mathematical model including reconstruction and islanding is as follows:
the method comprises the following steps of firstly, DistFlow power flow constraint, writing out a state equation of a system by taking active power, reactive power and a node voltage square as state variables of the system, and expressing the state equation as a constraint condition form. Considering the change of the patent in the aspect of network topology structure change during reconstruction and island division, a large M method is used for correction to obtain a corresponding node active power balance equation and a node reactive power balance equation:
wherein, the flow of the initial node to the final node in the branch information is specified to be a positive direction,represents the sum of the active power injected by the j node in the t period,representing the sum of the outgoing active power,representing the sum of the active line losses injected into the j node line at time t,representing lines with i node as the starting node and j node as the end node during t periodThe square term of the line current is,the j node predicts the active load representing the t period,represents the sum of all DG active power output of j nodes in the period t,representing the j node active load shedding amount in the t period,representing the sum of the reactive power injected at node j during time t,representing the sum of the reactive power flowing out,representing the sum of the reactive line losses injected into the j node line at time t,the reactive load is predicted at node j on behalf of time t,represents the sum of all DG reactive power outputs of j nodes in the t period,representing the reactive load throwing amount of the j node in the t period;
and a line pressure drop equilibrium equation, the following equation constraints and model conversion processing being applicable to any time period t, for each branch of each time period:
wherein M represents a value of largeThe constant number is a constant number,represents the square term of the voltage of the node i in the time period t (to eliminate the square term and achieve the linearization effect),representing the opening and closing condition of a branch circuit with an i node as a starting node and a j node as a tail node in a period t, wherein 0 or 1 represents that the circuit is opened, and 1 represents that the circuit is closed;
and secondly, limiting the branch power, wherein the rated capacity of the line transmission power needs to be considered in the topology changing process because the line transmission power of the actual power distribution network is limited. Secondly, the branch circuit state of opening and shutting decides whether line power can circulate, and when the branch circuit disconnection, branch circuit transmission power should be clamped as zero:
wherein,the active power and the reactive power circulating on the circuit in the time period t of the branch with i as the starting node and j as the last node are respectively represented, and the equation (6) restricts the circulating condition of the circuit by using a large M method, wherein M is a large number.
Node voltage quota restraint, need guarantee that electric wire netting node voltage is within the safety interval in the active power distribution network fault recovery period, voltage deviation does not exceed k% promptly, and general k gets 5 ~ 10:
and fourthly, DG capacity constraint, wherein the output of the distributed power supply connected into the power distribution network needs to be within a rated interval:
wherein,respectively representing the square terms of the active power and the reactive power output of the W-number distributed power supply in the t time period, S2 max DGSquare term representing rated capacity of load node, NDGRepresenting the total number of distributed power supplies;
current amplitude constraint, wherein the current of the line in each time period in the fault is required to be within a rated range, so that the current amplitude constraint is newly increased:
wherein,in the time period t, the square term of the current flowing through the line with i as a starting node and j as a final node is given;
sixthly, relaxing the secondary constraint into conical constraint:
radial and connectivity constraints are generated, in the process of dynamic reconfiguration and island division of a distribution network, a network structure needs to be kept radial, a virtual power flow model is built herein to standardize active power distribution network fault reconfiguration, the island division exists in an active power distribution network fault recovery optimization model, so that the distribution network is possibly divided into a plurality of areas, the number of root nodes is uncertain, the root nodes can be a transformer substation or a distributed power supply, the number of islands is also a part of optimization, the traditional radial and connectivity constraints cannot well adapt to the model built herein, and therefore the radiating and connectivity constraints of the island division suitable for the mixed integer second-order cone model are innovatively set herein:
wherein N isDGIs the total DG number, NsubIs the total number of substations, CqSelecting a variable for the root node, wherein the root node can be a DG (distributed generation) or a transformer substation for the 0-1 variable, and when the value of the variable C corresponding to the DG or the transformer substation is 1, the DG or the transformer substation is shown as the root node of one partition area, thenThe number of main nets additionally arranged outside the island is the number of the divided areas; a is one NlineLine, NDG+NsubThe 0-1 of the column optimizes the variable matrix.
Let NDG+Nsub=NrootThen in the active power distribution network fault recovery strategy, at most N occursrootA divided area including divided islands and main area left after division, Apq1 denotes a node for which the p-th node is selected as the q-th partition (q-th column), and if the distributed power supply is not selected as the root node of the partition, the corresponding C variable is set to zeroG corresponds to a node within the partition, i.e. a DG that is not selected does not partition out the partition. The expression (13) indicates that each node of the grid participating in the islanding optimization belongs to one and only one partition area. Equation (14) ensures that the root node corresponding to each partition is selected.
Therefore, in order to ensure radial constraint, namely, in the island division of the power grid, the number of the nodes of the power distribution network minus the number of the divided areas minus the number of the isolated nodes needs to be equal to the number of the running branches:
wherein,the number of main nets added outside the island; in the step a, columns with all zeros are removed, each column represents a partition area, so that connectivity in the area needs to be ensured during optimization, nodes in each partition area need to be connected with each other after the topology structure is changed, at this time, a model of "virtual power flow" is used to determine connectivity in an island, the load of a selected node in each partition area is set to be 1, power is injected from a root node, and if the connectivity in the partition area can be proved (a virtual power balance equation is satisfied):
wherein,for the virtual contribution of the root node (distributed power source),the line virtual power of i as the starting node and j as the end node,being k nodesAnd (4) virtual load. Compared with the common power flow, the virtual power flow greatly reduces the calculation amount and is more efficient.
The load shedding constraint is characterized in that loads in the active power distribution network can be divided into uncontrollable loads and controllable loads (flexible loads), the load shedding amount of the controllable loads is controllable, the uncontrollable loads have only two conditions, one is to shed all loads of the node, the other is not to execute the load shedding operation, and the load shedding amount has the following constraint requirements:
wherein, the subscript c represents the controllable load, the formula (16) is the load shedding constraint of the controllable load, the subscript uc represents the uncontrollable load, the formula (17) represents the load shedding constraint of the uncontrollable load, and z is a variable of 0-1.
As the optimization of the invention, the Mixed Integer Second Order Cone Problem (MISOCP) is solved, and an optimization result is obtained, so as to obtain a strategy for recovering unified operation according to dynamic information, network topology operation, DG output and load shedding in a fault time period, namely a fault recovery unified operation strategy of the active power distribution network.
Has the advantages that:
the method simultaneously considers two strategies of network fault reconstruction and island division, considers the controllability of load and distributed power supply output, enables the fault recovery strategy to be optimized and cover more possibilities compared with the traditional power distribution network fault recovery, converts the problem into mixed integer linear programming, saves the strategy generation time, and can more accurately obtain the optimal solution.
Detailed Description
The technical solution of the present invention is further described in detail with reference to the following specific embodiments.
The invention discloses a method for constructing a fault recovery unified model of an active power distribution network, which comprises the following steps of: the method comprises the following steps:
step 1) acquiring power grid fault information, a network topology structure, distribution network branches, a distributed power supply and switch state information, predicting load and fault time required in a fault period, and summarizing data required by model building. Discretizing the input fault information in the fault time interval, dividing the fault time interval into n time intervals, discretizing the corresponding predicted load in the same time scale and preparing for establishing a multi-time-interval dynamic reconstruction and island division model.
And 2) determining a target function of a fault recovery operation strategy, considering the economy of a power grid, considering multi-period load shedding cost, considering the switch operation cost, and considering the load shedding cost and the weight grade of the load.
And determining an objective function of the fault recovery operation strategy by taking the minimum load shedding amount as a main target and considering the switch operation cost, and considering the load shedding cost and the weight grade of the load. To ensure that the node load shedding amount is minimum, namely the recovery load amount is maximum, the predicted load is set as a constant, the load shedding amount of each node is set as an optimized variable, and meanwhile, the on-off state of a switch is set as a variable from 0 to 1, wherein 0 represents that the switch is off, and 1 represents that the switch is on. The target function considers the whole fault time interval, and the target function expression of the fault recovery operation strategy is as follows:
wherein, CopeRepresents the total cost, cloadRepresenting the load shedding cost factor, NTRepresenting the total number of time intervals of the fault period, NbusRepresenting the total number of nodes in the network, at representing the length of each time interval of the fault period,represents t time period, i node load shedding amount, lambdaiThe load weight of the ith node is obtained; c. CswitchRepresentative switch operating asThe coefficient, an optimization variable, alpha, of 0 to 1uThe switch state representing the fault recovery period of the u-th branch is a constant of 01,representing the switching state before fault in branch u, NlineRepresenting the total number of branches in the network.
Step 3) constructing a unified mathematical model containing reconstruction and island division, wherein the unified mathematical model contains DistFlow power flow constraint and connectivity constraint using virtual power flow; the established constraint conditions of the unified mathematical model containing reconstruction and island division are as follows:
the method comprises the following steps of firstly, DistFlow power flow constraint, writing out a state equation of a system by taking active power, reactive power and a node voltage square as state variables of the system, and expressing the state equation as a constraint condition form. Considering the change of the patent in the aspect of network topology structure change during reconstruction and island division, a large M method is used for correction to obtain a corresponding node active power balance equation and a node reactive power balance equation:
wherein, the flow of the initial node to the final node in the branch information is specified to be a positive direction,represents the sum of the active power injected by the j node in the t period,represents the sum of the active power flowing out from the j node in the t period,representing the sum of the active line losses injected into the j node line at time t,the square term of the line current of the line with the node i as the starting node and the node j as the end node in the period t is shown,the j node predicts the active load representing the t period,represents the sum of all DG active power output of j nodes in the period t,representing the j node active load shedding amount in the t period,representing the sum of the reactive power injected at node j during time t,representing the sum of the reactive power flowing out of node j during time t,representing the sum of the reactive line losses injected into the j node line at time t,the reactive load is predicted at node j on behalf of time t,represents the sum of all DG reactive power outputs of j nodes in the t period,representing the reactive load throwing amount of the j node in the t period;
and a line pressure drop equilibrium equation, the following equation constraints and model conversion processing being applicable to any time period t, for each branch of each time period:
wherein, M represents a constant with a large value,represents the square term of the voltage of the node i in the time period t (to eliminate the square term and achieve the linearization effect),and representing the opening and closing condition of the branch circuit with the node i as a starting node and the node j as a tail node in the period t, wherein 0 or 1, 0 represents that the circuit is opened, and 1 represents that the circuit is closed.
And secondly, limiting the branch power, wherein the rated capacity of the line transmission power needs to be considered in the topology changing process because the line transmission power of the actual power distribution network is limited. Secondly, the branch circuit state of opening and shutting decides whether line power can circulate, and when the branch circuit disconnection, branch circuit transmission power should be clamped as zero:
wherein,the active power and the reactive power circulating on the circuit in the time period t of the branch with i as the starting node and j as the last node are respectively represented, and the equation (6) restricts the circulating condition of the circuit by using a large M method, wherein M is a large number.
Node voltage quota restraint, need guarantee that electric wire netting node voltage is within the safety interval in the active power distribution network fault recovery period, voltage deviation does not exceed k% promptly, and general k gets 5 ~ 10:
and fourthly, DG capacity constraint, wherein the output of the distributed power supply connected into the power distribution network needs to be within a rated interval:
wherein,respectively representing the square terms of the active power and the reactive power output of the W-number distributed power supply in the t time period, S2 max DGSquare term representing rated capacity of load node, NDGRepresenting the total number of distributed power supplies;
current amplitude constraint, wherein the current of the line in each time period in the fault is required to be within a rated range, so that the current amplitude constraint is newly increased:
wherein,in the time period t, the square term of the current flowing through the line with i as the starting node and j as the end node,representing the branch current rating.
Sixthly, relaxing the secondary constraint into conical constraint:
radial and connectivity constraints are generated, in the process of dynamic reconfiguration and island division of a distribution network, a network structure needs to be kept radial, a virtual power flow model is built herein to standardize active power distribution network fault reconfiguration, the island division exists in an active power distribution network fault recovery optimization model, so that the distribution network is possibly divided into a plurality of areas, the number of root nodes is uncertain, the root nodes can be a transformer substation or a distributed power supply, the number of islands is also a part of optimization, the traditional radial and connectivity constraints cannot well adapt to the model built herein, and therefore the radiating and connectivity constraints of the island division suitable for the mixed integer second-order cone model are innovatively set herein:
wherein N isDGIs the total DG number, NsubIs the total number of substations, CqSelecting a variable for the root node, wherein the root node can be a DG (distributed generation) or a transformer substation for the 0-1 variable, and when the value of the variable C corresponding to the DG or the transformer substation is 1, the DG or the transformer substation is shown as the root node of one partition area, thenThe number of main nets additionally arranged outside the island is the number of the divided areas; a is one NlineLine, NDG+NsubThe 0-1 of the column optimizes the variable matrix.
Let NDG+Nsub=NrootAnd in the active power distribution network fault recovery strategy, the mostMultiple occurrence of NrootA divided area including divided islands and main area left after division, ApqBy 1, it means that the p-th node is selected as a node of the q-th partition area (q-th column), and if the distributed power supply is not selected as a root node of the partition area, the corresponding C variable is set to zero. The expression (13) indicates that each node of the grid participating in the islanding optimization belongs to one and only one partition area. Equation (14) ensures that the root node corresponding to each partition is selected.
Therefore, in order to ensure radial constraint, namely in the island division of the power grid, the number of the nodes of the power distribution network minus the number of the division areas needs to be equal to the number of the operating branches:
wherein,the number of main nets added outside the island; in the step a, columns with all zeros are removed, each column represents a partition area, so that connectivity in the area needs to be ensured during optimization, nodes in each partition area need to be connected with each other after the topology structure is changed, at this time, a model of "virtual power flow" is used to determine connectivity in an island, the load of a selected node in each partition area is set to be 1, power is injected from a root node, and if the connectivity in the partition area can be proved (a virtual power balance equation is satisfied):
wherein,as root node (distributed power supply)) The virtual output of (a) is determined,for a period t, i is taken as a starting node, j is taken as a line virtual power of an end node,is the virtual load of the k node. Compared with the common power flow, the virtual power flow greatly reduces the calculation amount and is more efficient.
The load shedding constraint is characterized in that loads in the active power distribution network can be divided into uncontrollable loads and controllable loads (flexible loads), the load shedding amount of the controllable loads is controllable, the uncontrollable loads have only two conditions, one is to shed all loads of the node, the other is not to execute the load shedding operation, and the load shedding amount has the following constraint requirements:
wherein, the subscript c represents the controllable load, the formula (16) is the load shedding constraint of the controllable load, the subscript uc represents the uncontrollable load, the formula (17) represents the load shedding constraint of the uncontrollable load, and z is a variable of 0-1.
And 4) solving the established unified mathematical model containing reconstruction and island division, and obtaining a final active power distribution network fault recovery operation strategy. And solving the Mixed Integer Second Order Cone Problem (MISOCP), obtaining an optimization result, and obtaining a dynamic information, network topology operation, DG output and load shedding strategy according to the fault time period, namely a fault recovery unified operation strategy of the active power distribution network.
The above description is only a preferred embodiment of the present invention, and the scope of the present invention is not limited to the above embodiment, but equivalent modifications or changes made by those skilled in the art according to the present disclosure should be included in the scope of the present invention as set forth in the appended claims.
Claims (4)
1. A method for constructing a unified model for fault recovery of an active power distribution network simultaneously considers two topological transformation strategies of reconstruction and islanding, and converts a non-convex nonlinear optimization problem into a solvable mixed integer second-order cone problem by setting various constraint conditions, and is characterized in that: the method comprises the following steps:
step 1) acquiring power grid fault information, a network topology structure, distribution network branches, a distributed power supply and switch state information, predicting load and fault time required by a fault period, and summarizing data required by model building;
step 2) determining a target function of the fault recovery operation strategy;
step 3) constructing a unified mathematical model containing reconstruction and island division, wherein the unified mathematical model contains DistFlow power flow constraint and connectivity constraint using virtual power flow; the constraint conditions of the established unified mathematical model containing the reconstruction and the island division are as follows:
firstly, DistFlow power flow constraint, taking active power, reactive power and node voltage square as state variables of a system, writing out a state equation of the system in a column, expressing the state equation as a constraint condition form, and correcting by using a large M method to obtain a corresponding node active power balance equation and a node reactive power balance equation:
wherein, the flow of the initial node to the final node in the branch information is specified to be a positive direction,representing j nodes in t periodThe sum of the active power is injected,represents the sum of the active power flowing out from the j node in the t period,representing the sum of the active line losses injected into the j node line at time t,the square term of the line current of the line with the node i as the starting node and the node j as the end node in the period t is shown,the j node predicts the active load representing the t period,represents the sum of all DG active power output of j nodes in the period t,representing the j node active load shedding amount in the t period,representing the sum of the reactive power injected at node j during time t,representing the sum of the reactive power flowing out of node j during time t,representing the sum of the reactive line losses injected into the j node line at time t,the reactive load is predicted at node j on behalf of time t,represents the sum of all DG reactive power outputs of j nodes in the t period,represents the reactive load rejection of the j node in the t period, NTRepresenting the total number of time intervals of the fault period;
and a line pressure drop equilibrium equation, the following equation constraints and model conversion processing being applicable to any time period t, for each branch of each time period:
wherein M represents a number having a large value,represents the squared term of the voltage at node i during the t period,representing the opening and closing condition of a branch circuit with an i node as a starting node and a j node as a tail node in a t period, wherein 0 or 1 represents that the circuit is opened and 1 represents that the circuit is closed;
secondly, the branch power limitation constraint needs to consider the rated capacity of the transmission power of the line in the topology changing process, and then the branch switching state determines whether the line power can circulate, and when the branch is switched off, the transmission power of the branch is clamped to zero:
wherein,respectively representing active power and reactive power circulating on a line in a time period t by taking i as a starting node and taking j as a final node, and constraining the circulation condition of the line by using a large M method in the formula (6), wherein M is a large number;
and limiting the node voltage, and ensuring that the node voltage of the power grid is within a safety interval within the fault recovery period of the active power distribution network, namely the voltage deviation is not more than k%, and k is 5-10:
and fourthly, DG capacity constraint, wherein the output of the distributed power supply connected into the power distribution network needs to be within a rated interval:
wherein,respectively representing the square terms of the active power and the reactive power output of the W-number distributed power supply in the t time period, S2 maxDGSquare term representing rated capacity of load node, NDGRepresenting the total number of distributed power supplies;
current amplitude constraint, wherein the current of the line in each time period in the fault is required to be within a rated range, so that the current amplitude constraint is newly increased:
wherein,in the time period t, the square term of the current flowing through the line with i as the starting node and j as the end node,representing the rated current of the branch circuit;
sixthly, relaxing the secondary constraint into conical constraint:
and seventhly, radial and connectivity constraint, namely setting island division radial and connectivity constraint applicable to a mixed integer second-order cone model:
wherein N isDGIs the total DG number, NsubIs the total number of substations, CqSelecting a variable for the root node, wherein the variable is a 0-1 variable, the root node can be a DG or a substation, and when the DG or the corresponding variable C of the substationqThe value is 1, and the DG or the transformer substation is a root node of a partition area; a is one NlineLine, NDG+Nsub0-1 optimization variable matrix of columns, NbusRepresenting the total number of nodes in the network;
let NDG+Nsub=NrootThen in the active power distribution network fault recovery strategy, at most N occursrootA divided area including divided islands and main area left after division, Apq1 indicates that the node p is selected to become a partition area q, if the distributed power supply is not selected to become a root node of the partition area, the corresponding C variable is set to zero, and according to the limitation, no load node is selected to become a node in the partition area corresponding to the unselected DG, that is, the unselected DG does not partition the partition area; the formula (13) represents that each node of the power grid participating in the islanding optimization belongs to one and only one partition area; equation (14) ensures that the root node corresponding to each partition region is selected;
therefore, in order to ensure radial constraint, namely in the island division of the power grid, the number of the nodes of the power distribution network minus the number of the division areas needs to be equal to the number of the operating branches:
wherein,the number of main nets added outside the island; in the step A, columns which are all zero are removed, each column represents a partition area, so that the connectivity in the area needs to be ensured during optimization, nodes in each partition area need to be mutually communicated after the topological structure is changed, at the moment, the connectivity in an island is determined by using a virtual power flow model, the load of a selected node in each partition area is set to be 1, power is injected from a root node, and if the power can be circulated and a virtual power balance equation is met, the connectivity in the partition areas can be proved:
wherein,is the virtual output of the root node, the distributed power supply,for a period t, i is taken as a starting node, j is taken as a line virtual power of an end node,is the virtual load of the k node;
and eighthly, load shedding constraint, namely the controllable load shedding amount is a controllable amount, and for an uncontrollable load, only two conditions are provided, one is to shed all loads of the node, and the other is not to execute the load shedding operation, so that the following constraint requirements are provided for the load shedding amount:
wherein, the subscript c represents controllable load, the formula (16) is load shedding constraint of the controllable load, the subscript uc represents uncontrollable load, the formula (17) represents load shedding constraint of the uncontrollable load, and z is a variable of 0-1;
and 4) solving the established unified mathematical model containing reconstruction and island division, and obtaining a final active power distribution network fault recovery operation strategy.
2. The method for constructing the unified model for fault recovery of the active power distribution network according to claim 1, wherein the unified model comprises the following steps: in the step 1), the input fault time interval fault information is discretized, the fault time interval is divided into n time intervals, the corresponding predicted load is discretized in the same time scale, and preparation is made for establishing a multi-time-interval dynamic unified mathematical model containing reconstruction and island division.
3. The method for constructing the unified model for fault recovery of the active power distribution network according to claim 1, wherein the unified model comprises the following steps: in the step 2), the predicted load is set as a constant, the load shedding amount of each node is set as an optimized variable, the on-off state of the switch is set as a variable of 0-1, 0 represents that the switch is off, 1 represents that the switch is on, and the objective function expression is as follows:
wherein, CopeRepresents the total cost, cloadRepresenting the load shedding cost factor, NTRepresenting the total number of time intervals of the fault period, NbusRepresenting the total number of nodes in the network, at representing the length of each time interval of the fault period,represents t time period, i node load shedding amount, lambdaiThe load weight of the ith node is obtained; c. CswitchRepresents a switch operation cost coefficient, and is an optimized variable of 0-1, alphauThe switch state representing the fault recovery time period of the u-th branch is a constant value of 0-1,representing the switching state before fault in branch u, NlineRepresenting the total number of branches in the network.
4. The method for constructing the unified model for failure recovery of the active power distribution network according to claim 1, wherein a mixed integer second order cone problem is solved, an optimization result is obtained, and a dynamic information, a network topology operation, a DG output and load shedding strategy according to a failure time period, namely a unified operation strategy for failure recovery of the active power distribution network, is obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910406777.0A CN110350508B (en) | 2019-05-16 | 2019-05-16 | Method for constructing unified fault recovery model of active power distribution network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910406777.0A CN110350508B (en) | 2019-05-16 | 2019-05-16 | Method for constructing unified fault recovery model of active power distribution network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110350508A CN110350508A (en) | 2019-10-18 |
CN110350508B true CN110350508B (en) | 2021-02-19 |
Family
ID=68174715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910406777.0A Expired - Fee Related CN110350508B (en) | 2019-05-16 | 2019-05-16 | Method for constructing unified fault recovery model of active power distribution network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110350508B (en) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110729770B (en) * | 2019-10-24 | 2022-07-29 | 北京交通大学 | Active power distribution network load fault recovery strategy optimization algorithm |
CN110912137B (en) * | 2019-10-25 | 2023-05-23 | 国网天津市电力公司电力科学研究院 | Flexible power distribution network operation domain model construction method considering alternating current power flow |
CN111027184B (en) * | 2019-11-21 | 2023-02-28 | 中国船舶重工集团公司第七一九研究所 | Ship power grid fault reconstruction convex optimization model considering reliability constraint |
CN111864791B (en) * | 2019-12-31 | 2022-04-01 | 武汉大学 | Decision table-based planned island online division method for power distribution network |
CN111682524A (en) * | 2020-05-28 | 2020-09-18 | 天津大学 | Multi-stage rapid reconstruction method for power distribution network |
CN111817345A (en) * | 2020-07-19 | 2020-10-23 | 国网山西省电力公司电力科学研究院 | Reconstruction method for power distribution network with distributed power supply after serious fault |
CN112186744B (en) * | 2020-09-16 | 2024-03-08 | 国网天津市电力公司 | Power supply recovery method suitable for power distribution network with distributed power supply and application |
CN112598159B (en) * | 2020-12-02 | 2022-04-19 | 广西大学 | Island fusion recovery method based on multi-source network reconstruction |
CN112510708B (en) * | 2020-12-15 | 2023-05-16 | 南方电网科学研究院有限责任公司 | Method and device for reconstructing radioactive power distribution network with variable number of connected areas |
CN112906169B (en) * | 2020-12-18 | 2022-03-22 | 天津大学 | Fault positioning method based on virtual stream |
CN112909942B (en) * | 2021-04-22 | 2023-09-26 | 国网宁夏电力有限公司电力科学研究院 | Active power distribution network fault recovery method and system based on hybrid solving strategy |
CN113572154B (en) * | 2021-07-08 | 2023-08-25 | 国网河北省电力有限公司电力科学研究院 | Power distribution network fault recovery method and device and electronic equipment |
CN113972693A (en) * | 2021-11-15 | 2022-01-25 | 国网江苏省电力有限公司常州供电分公司 | Power distribution network island dynamic division method |
CN114050573B (en) * | 2021-11-29 | 2024-07-02 | 杭州电力设备制造有限公司 | Fault recovery control method, device, equipment and medium for active power distribution network |
CN114336613A (en) * | 2022-01-04 | 2022-04-12 | 国网浙江省电力有限公司绍兴市上虞区供电公司 | Power distribution network reconstruction optimization method considering network segmentation uniformity |
CN116388185B (en) * | 2023-06-05 | 2023-09-12 | 昆明理工大学 | Active power distribution network fault processing and rapid self-healing method and system |
CN116780529B (en) * | 2023-06-30 | 2024-06-11 | 国网北京市电力公司 | Power distribution network fault recovery method, device, equipment and medium |
CN117154726A (en) * | 2023-11-01 | 2023-12-01 | 广东电网有限责任公司江门供电局 | Power distribution network reconstruction power restoration method and system based on island division |
CN117394353B (en) * | 2023-12-08 | 2024-05-14 | 国网天津市电力公司电力科学研究院 | Power distribution network load transferring and recovering method and device |
CN117810996B (en) * | 2024-03-01 | 2024-05-17 | 南京师范大学 | Active power distribution network fault recovery method considering island division and network reconstruction |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109120009A (en) * | 2018-09-11 | 2019-01-01 | 国网天津市电力公司电力科学研究院 | The active distribution network fault recovery method that meter and distributed generation resource power output change at random |
-
2019
- 2019-05-16 CN CN201910406777.0A patent/CN110350508B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110350508A (en) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110350508B (en) | Method for constructing unified fault recovery model of active power distribution network | |
Rouzbehi et al. | DC voltage control and power sharing in multiterminal DC grids based on optimal DC power flow and voltage-droop strategy | |
CN110266038A (en) | A kind of more virtual plant distributed coordination regulation methods | |
Sekhar et al. | An enhanced cuckoo search algorithm based contingency constrained economic load dispatch for security enhancement | |
CN104037765B (en) | The method of active power distribution network service restoration scheme is chosen based on improved adaptive GA-IAGA | |
CN112271727B (en) | Fault recovery method for flexible power distribution network containing flexible soft switch | |
CN103903055B (en) | Network reconstruction method based on all spanning trees of non-directed graph | |
El Helou et al. | Fully decentralized reinforcement learning-based control of photovoltaics in distribution grids for joint provision of real and reactive power | |
Jiang et al. | A hierarchical control structure for distributed energy storage system in DC micro-grid | |
CN107529678B (en) | Multi-objective optimization control method capable of expanding safety and stability constraint for alternating current and direct current power grid | |
CN107394775A (en) | A kind of active power distribution network sequential islet operation method of meter and load priority | |
CN107332277B (en) | Active power distribution network island operation method considering source load storage operation characteristics | |
CN112467748A (en) | Double-time-scale distributed voltage control method and system for three-phase unbalanced active power distribution network | |
CN107451743B (en) | Active power distribution network island operation method considering source load operation characteristics | |
CN113872213B (en) | Autonomous optimization control method and device for power distribution network voltage | |
Li et al. | Distributed control strategy for global economic operation and bus restorations in a hybrid AC/DC microgrid with interconnected subgrids | |
Taheri et al. | A new modified TLBO algorithm for placement of AVRs in distribution system | |
CN113097994A (en) | Power grid operation mode adjusting method and device based on multiple reinforcement learning agents | |
CN111860969A (en) | Power transmission network extension planning method based on reinforcement learning algorithm | |
Qin et al. | Tertiary control based on non-parametric model prediction for dc microgrid cluster | |
CN114784796A (en) | Multi-stage recovery method for flexible interconnected power distribution system based on multi-terminal SOP | |
CN115765035A (en) | Flexible power distribution network disturbance recovery method suitable for full-time dynamic reconstruction | |
CN116260201A (en) | Fault recovery method, device, equipment and storage medium for AC/DC hybrid power distribution network | |
Ling et al. | A distributed model for power system restoration based on ant colony optimization algorithm | |
Liu et al. | Optimal power flow of DC-grid based on improved PSO algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210219 |