CN108985561B - Active power distribution network island division method based on opportunity constraint - Google Patents
Active power distribution network island division method based on opportunity constraint Download PDFInfo
- Publication number
- CN108985561B CN108985561B CN201810588310.8A CN201810588310A CN108985561B CN 108985561 B CN108985561 B CN 108985561B CN 201810588310 A CN201810588310 A CN 201810588310A CN 108985561 B CN108985561 B CN 108985561B
- Authority
- CN
- China
- Prior art keywords
- power
- constraint
- node
- active power
- scene
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000010248 power generation Methods 0.000 claims abstract description 18
- 238000011084 recovery Methods 0.000 claims description 31
- 238000013486 operation strategy Methods 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
An active power distribution network island division method based on opportunity constraint comprises the following steps: according to the selected power distribution system, input line parameters, load levels, network topology connection relations, system operation voltage levels, branch active power limits, confidence parameters, controllable and uncontrollable distributed power access positions and capacities, a load prediction curve, system fault moments, reference voltages and reference power initial values; discretizing the power generation power of the uncontrollable distributed power supply to obtain a discrete probability distribution parameter of the power generation power of the uncontrollable distributed power supply; establishing an active power distribution network island division model based on opportunity constraint; converting system power flow constraint, node voltage probability constraint and branch active power probability constraint to obtain a mixed integer nonlinear model; solving the mixed integer nonlinear model by adopting an interior point method; and outputting a solving result. The invention can work out the island division strategy meeting different requirements, and ensures the safe operation of the system while ensuring the reliable power supply of part of important loads.
Description
Technical Field
The invention relates to an active power distribution network island division method. In particular to an active power distribution network island division method based on opportunity constraint.
Background
With the high-proportion wide access of distributed power supplies, the operation and scheduling modes of the power distribution system are changed profoundly and durably. Island operation is a special operation mode of an active power distribution system, and can supply power for important loads in the system by means of a distributed power supply under the condition of extreme faults, so that the power supply reliability and the power supply elasticity of the power distribution system are effectively improved.
However, after a large number of uncontrollable distributed power sources such as fans and photovoltaic power sources are connected to an active power distribution network, the operating characteristics of the uncontrollable distributed power sources are greatly influenced by the environment and have obvious randomness and volatility, and many problems such as node voltage out-of-limit, branch current overload and the like can be brought to the safe operation of the power distribution network. In the island operation process of the power distribution network, when the output of the uncontrollable distributed power supply is different from a predicted value, the controllable distributed power supply is limited in capacity and is often difficult to ensure power balance in an island region through self power regulation, so that a preset island strategy is difficult to normally operate, the problems of voltage out-of-limit and current overload in the power distribution network are more obvious, and the operation safety and the power supply reliability of the power distribution system are seriously influenced. Therefore, the influence of the output uncertainty on the island operation strategy is fully considered when the island division is carried out on the power distribution network containing the uncontrollable distributed power supply.
At present, the existing island division method considering the uncertainty of the distributed power output can be roughly divided into a robust optimization algorithm and a random programming method. The robust optimization algorithm replaces the probability exact distribution of the power generation power of the uncontrollable distributed power source with the uncertain set, takes the optimized operation strategy in the worst scene as the operation strategy of the system, and obtains a conservative result when the method is used for ensuring that the load in the system is supplied with power continuously as much as possible. The stochastic programming method is characterized in that uncertain information of the uncontrollable distributed power supply is mostly described in a scene probability distribution mode, when an island is divided by the stochastic programming method, an opportunity constraint model is mostly adopted, the model allows that an island operation strategy does not meet constraint conditions to a certain extent, but the probability of the constraint conditions being met is not lower than a certain confidence level. However, the existing active power distribution network island division method based on opportunity constraint mostly adopts a heuristic/intelligent algorithm + check two-stage algorithm, and the method has more iteration times, complicated calculation process, difficulty in finding a global optimal solution and certain limitation.
Therefore, a mathematical programming method capable of being directly solved is urgently needed, and is used for solving an active power distribution network island division model based on opportunity constraint, so that an island division strategy capable of meeting different requirements is worked out, the confidence level and the load recovery level of safe operation of a system are fully coordinated, and the system is enabled to operate safely while reliable power supply of part of important loads is ensured.
Disclosure of Invention
The invention aims to solve the technical problem of providing an active power distribution network island division method based on opportunity constraint, which can meet different requirements of an active power distribution network.
The technical scheme adopted by the invention is as follows: an active power distribution network island division method based on opportunity constraint comprises the following steps:
1) according to the selected power distribution system, input line parameters, load levels, network topology connection relations, system operation voltage levels, branch active power limits, confidence parameters, controllable and uncontrollable distributed power access positions and capacities, a load prediction curve, system fault moments, reference voltages and reference power initial values;
2) discretizing the power generation power of the uncontrollable distributed power supply to obtain a discrete probability distribution parameter of the power generation power of the uncontrollable distributed power supply;
3) establishing an active power distribution network island division model based on opportunity constraint, comprising the following steps: setting the maximum active load quantity recovered by a power distribution system in a period of time as an objective function, and respectively considering radial constraint, system power flow constraint, network reconstruction constraint, node voltage probability constraint, branch active power probability constraint, uncontrollable distributed power supply operation constraint and controllable distributed power supply operation constraint;
4) converting system power flow constraint, node voltage probability constraint and branch active power probability constraint to obtain a mixed integer nonlinear model;
5) solving the mixed integer nonlinear model obtained in the step 4) by adopting an interior point method;
6) and outputting the solving result of the step 5), including the power distribution system active load recovery quantity, the sectionalizing switch and interconnection switch states and the recovery load coefficient of each node which meet any confidence level.
The method for discretizing the power generation power of the uncontrollable distributed power supply in the step 2) comprises the following steps:
wherein, f (x) is a probability density function of the generated power of the uncontrollable distributed power supply; x is the generated power of the uncontrollable distributed power supply; m is a scene of discrete power generation power of the uncontrollable distributed power supply; p (m) is the generated power of the uncontrollable distributed power supply under the scene of m; f [ P (m)]Representing the probability that the power generation power of the uncontrollable distributed power source is P (m); q is a discretization step length; omegasIs a collection of scenes m.
The system power flow constraint in the step 3) is as follows:
in the formula, omegabIs a set of branches; pt,ji,m、Qt,ji,mRespectively the active power and the reactive power flowing through the m scene branch ji in the t time period; pt,i,m、Qt,i,mThe sum of active power and reactive power injected at the time t and the m scene node iSum of rates;respectively the active power and the reactive power consumed by the load on the node i in the t period; lambda [ alpha ]iFor the recovery of the load on node i, λi∈{0,1},λ i1 denotes the load recovery on node i, λ i0 indicates that the load on node i is not recovered;respectively the active power and the reactive power injected by the controllable distributed power supply at the time period t and the scene node m;respectively the active power and the reactive power injected by the uncontrollable distributed power supply on the m scene node i in the t time period; u. oft,i,mThe square of the voltage amplitude on the m scene node i in the time period t; i.e. it,ij,mThe current amplitude of the m scene branch ij is squared at the time t; rijIs the resistance, X, of branch ijijIs the reactance of branch ij.
The node voltage probability constraint in step 3) is as follows:
in the formula, Pr{. denotes the probability that an event holds; epsilon is a confidence parameter;respectively the upper and lower limits of the voltage amplitude; u. oft,i,mThe square of the voltage amplitude on the m scene node i in the time period t.
Step 3) the branch active power probability constraint is as follows:
in the formula, Pr{. represents a certain factProbability of a condition being established; epsilon is a confidence parameter;the upper limit and the lower limit of the active power of the branch circuit are respectively set; pt,ij,mThe active power flowing through the m scene branch ij in the time period t.
The transformation in the step 4) comprises the following steps:
introduction of the variable z from 0 to 1t,mAnd jointly converting system power flow constraint, node voltage probability constraint and branch active power probability constraint into a mixed integer nonlinear model:
in the formula, zt,mFor the introduced binary variable, zt,mThe expression 0 includes the scenario m, z when solving the optimal island operation strategy t,m1 represents that a scene m is not taken into account when solving the optimal island operation strategy; pit,mIs the probability of occurrence of the scene m at the time t; m represents a very large constant; omegabIs a set of branches; pt,ji,m、Qt,ji,mRespectively the active power and the reactive power flowing through the m scene branch ji in the t time period; pt,i,m、Qt,i,mRespectively being the sum of active power and reactive power injected at the t time period and the m scene node i;respectively the active power and the reactive power consumed by the load on the node i in the t period; lambda [ alpha ]iFor the recovery of the load on node i, λi∈{0,1},λ i1 denotes the load recovery on node i, λ i0 indicates that the load on node i is not recovered;respectively the active power and the reactive power injected by the controllable distributed power supply at the time period t and the scene node m;respectively the active power and the reactive power injected by the uncontrollable distributed power supply on the m scene node i in the t time period; u. oft,i,mThe square of the voltage amplitude on the m scene node i in the time period t; i.e. it,ij,mThe current amplitude of the m scene branch ij is squared at the time t; rijIs the resistance, X, of branch ijijReactance for branch ij; epsilon is a confidence parameter;respectively the upper and lower limits of the voltage amplitude;respectively the upper and lower limits of the active power of the branch.
The invention discloses an active power distribution network island division method based on opportunity constraint, which is based on solving the island division problem of an active power distribution network containing an uncontrollable distributed power supply, fully considering network topology constraint, system power flow constraint, node voltage probability constraint, branch active power probability constraint, uncontrollable distributed power supply operation constraint and controllable distributed power supply operation constraint, establishing an active power distribution network island division model based on opportunity constraint, solving by adopting an interior point method, and obtaining an island division strategy meeting node voltage and branch active power probability confidence level, wherein the mathematical essence of the model is a mixed integer nonlinear programming problem. The invention can work out the island division strategy meeting different requirements, fully coordinate the confidence level and the load recovery level of the safe operation of the system, and ensure the reliable power supply of part of important loads and ensure the safe operation of the system.
Drawings
FIG. 1 is a flow chart of an active power distribution network islanding method based on opportunity constraint according to the invention;
FIG. 2 is a diagram of an improved IEEE33 node algorithm;
FIG. 3 is a load prediction curve;
FIG. 4 is a probability density function distribution plot of the generated power of the photovoltaic unit at 9;
FIG. 5 is a discrete probability distribution graph of the power generated by the photovoltaic generator at 9;
FIG. 6 is a schematic diagram of an islanding strategy of scenario 1 when a system fails at 9;
FIG. 7 is a schematic diagram of an islanding strategy of scenario 2 when a system fails at 9;
FIG. 8 is a schematic diagram of an islanding strategy of scenario 3 when a system fails at 9;
FIG. 9 is a schematic diagram of the relationship between the load recovery level and the confidence level when the system fails at 9.
Detailed Description
The invention provides an active power distribution network islanding method based on opportunity constraint, which is described in detail below with reference to embodiments and drawings.
As shown in fig. 1, the method for dividing an island of an active power distribution network based on opportunity constraint of the present invention includes the following steps:
1) according to the selected power distribution system, input line parameters, load levels, network topology connection relations, system operation voltage levels, branch active power limits, confidence parameters, controllable and uncontrollable distributed power access positions and capacities, a load prediction curve, system fault moments, reference voltages and reference power initial values;
2) discretizing the power generation power of the uncontrollable distributed power supply to obtain a discrete probability distribution parameter of the power generation power of the uncontrollable distributed power supply; the method for discretizing the power generation power of the uncontrollable distributed power supply comprises the following steps:
wherein, f (x) is a probability density function of the generated power of the uncontrollable distributed power supply; x is the generated power of the uncontrollable distributed power supply; m is a scene of discrete power generation power of the uncontrollable distributed power supply; p (m) is the generated power of the uncontrollable distributed power supply under the scene of m; f [ P (m)]Representing the probability that the power generation power of the uncontrollable distributed power source is P (m); q is a discretization step length; omegasIs a collection of scenes m.
3) Establishing an active power distribution network island division model based on opportunity constraint, comprising the following steps: setting the maximum active load quantity recovered by a power distribution system in a period of time as an objective function, and respectively considering radial constraint, system power flow constraint, network reconstruction constraint, node voltage probability constraint, branch active power probability constraint, uncontrollable distributed power supply operation constraint and controllable distributed power supply operation constraint; wherein,
(1) the maximum recovery active load capacity of the power distribution system in a period of time is expressed as an objective function
In the formula, omegaτA set of power distribution system island operation times; omeganIs the collection of all nodes of the power distribution system; lambda [ alpha ]iFor the recovery of the load on node i, λi∈{0,1},λ i1 denotes the load recovery on node i, λ i0 indicates that the load on node i is not recovered;is the active load on node i during time t.
(2) The radial constraint is expressed as
αij=βij+βji,ij∈Ωb (3)
αij∈{0,1} (6)
0≤βij≤1,0≤βji≤1 (7)
In the formula, omegabRepresenting a collection of all branches of the power distribution system; omegarRepresenting a node set for supporting island voltage and frequency when the power distribution system operates in an island; alpha is alphaijRepresenting the open state of the switch on branch ij, alpha ij1 denotes switch closed, α ij0 indicates that the switch is open; beta is aijRepresents the relationship of node i and node j, β ij1 means that node j is the parent node of node i, otherwise βij=0。
(3) The system flow constraint is as follows:
in the formula, Pt,ji,m、Qt,ji,mRespectively the active power and the reactive power flowing through the m scene branch ji in the t time period; pt,i,m、Qt,i,mRespectively t period and m fieldThe sum of active power and reactive power injected at a scene node i;respectively the active power and the reactive power consumed by the load on the node i in the t period;respectively the active power and the reactive power injected by the controllable distributed power supply at the time period t and the scene node m;respectively the active power and the reactive power injected by the uncontrollable distributed power supply on the m scene node i in the t time period; u. oft,i,mThe voltage amplitude of the m scene node i is squared in the t time period; i.e. it,ij,mThe current amplitude of the m scene branch ij is squared at the time t; rijIs the resistance, X, of branch ijijIs the reactance of branch ij.
(4) The network reconfiguration constraint is expressed as
-Mαij≤Pt,ij,m≤Mαij (13)
-Mαij≤Qt,ij,m≤Mαij (14)
0≤it,ji,m≤Mαij (15)
In the formula, M represents a very large constant.
(5) The node voltage probability constraint is expressed as
In the formula, Pr{. denotes the probability that an event holds; epsilon is a confidence parameter;respectively the upper and lower limits of the voltage amplitude.
(6) The branch active power probability constraint is expressed as
(7) The controllable distributed power supply operation constraint is expressed as
In the formula,representing the capacity of the controllable distributed power supply on the node i;is the minimum power factor for the operation of the distributed power supply on node i.
(8) The uncontrollable distributed power supply operation constraint is expressed as
4) Converting system power flow constraint, node voltage probability constraint and branch active power probability constraint to obtain a mixed integer nonlinear model; the transformation comprises the following steps:
introduction of the variable z from 0 to 1t,mAnd jointly converting system power flow constraint, node voltage probability constraint and branch active power probability constraint into a mixed integer nonlinear model:
in the formula, zt,mFor the introduced binary variable, zt,mThe expression 0 includes the scenario m, z when solving the optimal island operation strategy t,m1 represents that a scene m is not taken into account when solving the optimal island operation strategy; pit,mIs the probability of occurrence of the scene m at the time t; m represents a very large constant.
5) Solving the mixed integer nonlinear model obtained in the step 4) by adopting an interior point method;
6) and outputting the solving result of the step 5), including the power distribution system active load recovery quantity, the sectionalizing switch and interconnection switch states and the recovery load coefficient of each node which meet any confidence level.
For the embodiment of the invention, firstly, the impedance value of a line element in an IEEE33 node system, the active power and the reactive power of a load element, the network topology connection relation, the confidence parameter, a load prediction curve and the distributed power supply parameter are input, the arithmetic structure is shown in figure 2, and the detailed parameters are shown in tables 1, 2 and 3; assuming that a permanent three-phase fault occurs between the branches 1-2 at 9 am, the fault isolation time is 1 hour; setting the reference voltage of the system to be 12.66kV and the reference power to be 1 MVA; to verify the validity of the method, the following 3 scenarios were used for analysis.
Scene 1: the island division strategy with the confidence level of 100 percent;
scene 2: the island division strategy corresponding to the confidence level of 90 percent;
scene 3: the island division strategy corresponds to the confidence level of 80%;
it is assumed that the uncontrollable distributed power sources accessed by the power distribution system are all photovoltaic units, and the probability density function of the generated power of the photovoltaic units is shown in fig. 4 at 9 am. The generated power of the photovoltaic unit is dispersed by the discretization method provided by the invention, and the obtained generated power discrete probability distribution is shown in figure 5. Solving is carried out by adopting the method, active power distribution network island division strategies under 3 scenes are shown in figures 6, 7 and 8, wherein solid nodes represent load recovery of the nodes, and hollow nodes represent load non-recovery of the nodes; the correlation between the load recovery amount and the confidence level is shown in fig. 9, and the load recovery conditions of different confidence levels are detailed in table 4; in order to verify the effectiveness of the invention, monte carlo tests are performed on the islanding strategy of 3 scenes, and the results are shown in table 5.
The computer hardware environment for executing the optimization calculation is Intel ICoreII5-3470CPU, the dominant frequency is 3.20GHz, and the memory is 4 GB; the software environment is the Windows 7 operating system.
As can be seen from fig. 6, 7 and 8, according to the invention, a corresponding islanding strategy of the active power distribution network can be formulated according to the node voltage and the probability confidence level that the branch active power constraint condition is established. When the confidence levels are different, the islanding strategies of the power distribution network are different, the load recovery amount is also different, the specific correlation between the load recovery level and the confidence level is shown in fig. 9, and the specific load recovery condition is shown in table 3. It can be seen that the active load recovery level of the power distribution system gradually decreases with the increase of the confidence level, because the higher the confidence level is, the more scenes the island operation strategy is suitable for are, and besides the scene that the generated power of the photovoltaic unit is higher, the island operation strategy must also meet the scene that the generated power of the photovoltaic unit is lower, so that the load recovery amount is reduced. To verify the effectiveness of the present invention, Monte Carlo tests were performed on 3 scenarios and the results are shown in Table 5. As can be seen from table 5, the active power distribution network islanding method based on opportunity constraint provided by the present invention can ensure that the obtained islanding operation strategy safely operates with a probability higher than a certain confidence level. The method can effectively coordinate the relationship between the confidence level of the safe operation of the system and the load recovery level, and formulate the corresponding island division strategy according to different requirements, thereby having important significance for the safe operation and reliable power supply of the power distribution system.
TABLE 1 IEEE33 node sample load access location and Power
TABLE 2 IEEE33 node exemplary line parameters
TABLE 3 distributed Power supply configuration
TABLE 4 load recovery scenarios with different confidence levels
TABLE 5 Monte Carlo test results
Claims (5)
1. An active power distribution network island division method based on opportunity constraint is characterized by comprising the following steps:
1) according to the selected power distribution system, input line parameters, load levels, network topology connection relations, system operation voltage levels, branch active power limits, confidence parameters, controllable and uncontrollable distributed power access positions and capacities, a load prediction curve, system fault moments, reference voltages and reference power initial values;
2) discretizing the power generation power of the uncontrollable distributed power supply to obtain a discrete probability distribution parameter of the power generation power of the uncontrollable distributed power supply;
3) establishing an active power distribution network island division model based on opportunity constraint, comprising the following steps: setting the maximum active load quantity recovered by a power distribution system in a period of time as an objective function, and respectively considering radial constraint, system power flow constraint, network reconstruction constraint, node voltage probability constraint, branch active power probability constraint, uncontrollable distributed power supply operation constraint and controllable distributed power supply operation constraint;
4) converting system power flow constraint, node voltage probability constraint and branch active power probability constraint to obtain a mixed integer nonlinear model; the transformation comprises the following steps:
introduction of the variable z from 0 to 1t,mAnd jointly converting system power flow constraint, node voltage probability constraint and branch active power probability constraint into a mixed integer nonlinear model:
in the formula, zt,mFor the introduced binary variable, zt,mThe expression 0 includes the scenario m, z when solving the optimal island operation strategyt,m1 represents that a scene m is not taken into account when solving the optimal island operation strategy; pit,mIs the probability of occurrence of the scene m at the time t; m represents a very large constant; omegabIs a set of branches; pt,ji,m、Qt,ji,mRespectively the active power and the reactive power flowing through the m scene branch ji in the t time period; pt,i,m、Qt,i,mRespectively being the sum of active power and reactive power injected at the t time period and the m scene node i;respectively the active power and the reactive power consumed by the load on the node i in the t period;λifor the recovery of the load on node i, λi∈{0,1},λi1 denotes the load recovery on node i, λi0 indicates that the load on node i is not recovered;respectively the active power and the reactive power injected by the controllable distributed power supply at the time period t and the scene node m;respectively the active power and the reactive power injected by the uncontrollable distributed power supply on the m scene node i in the t time period; u. oft,i,mThe square of the voltage amplitude on the m scene node i in the time period t; i.e. it,ij,mThe current amplitude of the m scene branch ij is squared at the time t; rijIs the resistance, X, of branch ijijReactance for branch ij; epsilon is a confidence parameter;respectively the upper and lower limits of the voltage amplitude;the upper limit and the lower limit of the active power of the branch circuit are respectively set;
5) solving the mixed integer nonlinear model obtained in the step 4) by adopting an interior point method;
6) and outputting the solving result of the step 5), including the power distribution system active load recovery quantity, the sectionalizing switch and interconnection switch states and the recovery load coefficient of each node which meet any confidence level.
2. The opportunity constraint-based active power distribution network islanding method according to claim 1, wherein the method for discretizing the power generation power of the uncontrollable distributed power source in the step 2) is as follows:
wherein, f (x) is a probability density function of the generated power of the uncontrollable distributed power supply; x is the generated power of the uncontrollable distributed power supply; m is a scene of discrete power generation power of the uncontrollable distributed power supply; p (m) is the generated power of the uncontrollable distributed power supply under the scene of m; f [ P (m)]Representing the probability that the power generation power of the uncontrollable distributed power source is P (m); q is a discretization step length; omegasIs a collection of scenes m.
3. The active power distribution network islanding method based on opportunity constraint according to claim 1, wherein the system power flow constraint in step 3) is as follows:
in the formula, omegabIs a set of branches; pt,ji,m、Qt,ji,mRespectively the active power and the reactive power flowing through the m scene branch ji in the t time period; pt,i,m、Qt,i,mRespectively being the sum of active power and reactive power injected at the t time period and the m scene node i;respectively the active power and the reactive power consumed by the load on the node i in the t period; lambda [ alpha ]iFor the recovery of the load on node i, λi∈{0,1},λi1 denotes the load recovery on node i, λi0 indicates that the load on node i is not recovered;respectively the active power and the reactive power injected by the controllable distributed power supply at the time period t and the scene node m;respectively the active power and the reactive power injected by the uncontrollable distributed power supply on the m scene node i in the t time period; u. oft,i,mThe square of the voltage amplitude on the m scene node i in the time period t; i.e. it,ij,mThe current amplitude of the m scene branch ij is squared at the time t; rijIs the resistance, X, of branch ijijIs the reactance of branch ij.
4. The opportunity constraint-based active power distribution network islanding method according to claim 1, wherein the node voltage probability constraint in step 3) is as follows:
5. The method for islanding an active power distribution network based on opportunity constraint according to claim 1, wherein the branch active power probability constraint in step 3) is as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810588310.8A CN108985561B (en) | 2018-06-08 | 2018-06-08 | Active power distribution network island division method based on opportunity constraint |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810588310.8A CN108985561B (en) | 2018-06-08 | 2018-06-08 | Active power distribution network island division method based on opportunity constraint |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108985561A CN108985561A (en) | 2018-12-11 |
CN108985561B true CN108985561B (en) | 2021-09-07 |
Family
ID=64540069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810588310.8A Expired - Fee Related CN108985561B (en) | 2018-06-08 | 2018-06-08 | Active power distribution network island division method based on opportunity constraint |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108985561B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109687469B (en) * | 2019-01-11 | 2023-04-28 | 国网天津市电力公司电力科学研究院 | Intelligent soft switching voltage control method for active power distribution network based on opportunity constraint planning |
CN110263391B (en) * | 2019-05-31 | 2023-05-26 | 广东电网有限责任公司 | Opportunity constraint-based intelligent soft switch planning method for active power distribution network |
CN112598159B (en) * | 2020-12-02 | 2022-04-19 | 广西大学 | Island fusion recovery method based on multi-source network reconstruction |
CN113011083B (en) * | 2021-02-25 | 2023-09-05 | 中国科学院电工研究所 | Island operation duration simulation evaluation method for comprehensive energy system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104868465A (en) * | 2014-02-26 | 2015-08-26 | 云南电力调度控制中心 | Power system grid structure reconfiguration and optimization method based on fuzzy chance constraint |
CN105958486A (en) * | 2016-06-15 | 2016-09-21 | 华北电力大学 | Power distribution network multi-period dynamic fault recovery method considering DG (Distributed Generation) output curve |
CN106921164A (en) * | 2017-04-05 | 2017-07-04 | 广东电网有限责任公司东莞供电局 | Mixed integer second-order cone planning method and system for distribution network voltage reactive power collaborative optimization |
-
2018
- 2018-06-08 CN CN201810588310.8A patent/CN108985561B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104868465A (en) * | 2014-02-26 | 2015-08-26 | 云南电力调度控制中心 | Power system grid structure reconfiguration and optimization method based on fuzzy chance constraint |
CN105958486A (en) * | 2016-06-15 | 2016-09-21 | 华北电力大学 | Power distribution network multi-period dynamic fault recovery method considering DG (Distributed Generation) output curve |
CN106921164A (en) * | 2017-04-05 | 2017-07-04 | 广东电网有限责任公司东莞供电局 | Mixed integer second-order cone planning method and system for distribution network voltage reactive power collaborative optimization |
Non-Patent Citations (3)
Title |
---|
Resilience-Oriented Critical Load Restoration Using Microgrids in Distribution Systems;Haixiang Gao等;《IEEE TRANSACTIONS ON SMART GRID》;20161130;第7卷(第6期);第2837-2848页 * |
基于机会约束规划含光伏发电的配电网故障恢复;陈昕玥 等;《电网技术》;20140131;第38卷(第1期);第99-106页 * |
考虑分布式发电波动性的有源配电网故障恢复策略;张磐 等;《电力***及其自动化学报》;20180131;第30卷(第1期);第115-120页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108985561A (en) | 2018-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108985561B (en) | Active power distribution network island division method based on opportunity constraint | |
Zhang et al. | Hierarchically-coordinated voltage/VAR control of distribution networks using PV inverters | |
Orfanos et al. | Transmission expansion planning of systems with increasing wind power integration | |
Shen et al. | Low-carbon electricity network transition considering retirement of aging coal generators | |
Vazquez et al. | A fully decentralized adaptive droop optimization strategy for power loss minimization in microgrids with PV-BESS | |
Dent et al. | Efficient secure AC OPF for network generation capacity assessment | |
Nassar et al. | A novel load flow algorithm for islanded AC/DC hybrid microgrids | |
Wang et al. | Optimal energy storage allocation for mitigating the unbalance in active distribution network via uncertainty quantification | |
Wang et al. | Risk assessment of distribution networks integrating large-scale distributed photovoltaics | |
Zhao et al. | An islanding partition method of active distribution networks based on chance-constrained programming | |
Galvani et al. | Optimal allocation of static synchronous series compensator (SSSC) in wind-integrated power system considering predictability | |
Abid et al. | Chaotic slime mould optimization algorithm for optimal load-shedding in distribution system | |
Parastvand et al. | Cyber-security constrained placement of FACTS devices in power networks from a novel topological perspective | |
Guo et al. | Short-term reliability assessment for islanded microgrid based on time-varying probability ordered tree screening algorithm | |
Davoudi et al. | Increasing penetration of distributed generation with meshed operation of distribution systems | |
Wei et al. | Reliability assessment for AC/DC hybrid distribution network with high penetration of renewable energy | |
Ivic et al. | An extended control strategy for weakly meshed distribution networks with soft open points and distributed generation | |
Ehsan et al. | Robust distribution system planning considering the uncertainties of renewable distributed generation and electricity demand | |
Yuan et al. | Dynamic equivalent model of VSC based on singular perturbation | |
Preece et al. | Probabilistic risk assessment of rotor angle instability using fuzzy inference systems | |
Taye et al. | A dynamic droop control for a DC microgrid to enhance voltage profile and proportional current sharing | |
Yang et al. | A novel partitioning method for the power grid restoration considering the support of multiple LCC-HVDC systems | |
Wang et al. | Unified load flow calculation for flexible interconnected AC-MTDC distribution system considering control strategy switching with high penetration of DGs | |
Nabi et al. | Fast distributed control design for DC linked microgrids | |
Davari et al. | A novel IGDT-based method to find the most susceptible points of cyberattack impacting operating costs of VSC-based microgrids |
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: 20210907 |