CN107069801A - A kind of power distribution network isolated island division methods based on minimum Custom interruption cost - Google Patents
A kind of power distribution network isolated island division methods based on minimum Custom interruption cost Download PDFInfo
- Publication number
- CN107069801A CN107069801A CN201710144090.5A CN201710144090A CN107069801A CN 107069801 A CN107069801 A CN 107069801A CN 201710144090 A CN201710144090 A CN 201710144090A CN 107069801 A CN107069801 A CN 107069801A
- Authority
- CN
- China
- Prior art keywords
- node
- isolated island
- load
- distribution network
- power source
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000005611 electricity Effects 0.000 claims abstract description 35
- 238000000638 solvent extraction Methods 0.000 claims abstract description 31
- 238000011144 upstream manufacturing Methods 0.000 claims description 7
- 238000009472 formulation Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 abstract description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004146 energy storage Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011234 economic evaluation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- GOLXNESZZPUPJE-UHFFFAOYSA-N spiromesifen Chemical compound CC1=CC(C)=CC(C)=C1C(C(O1)=O)=C(OC(=O)CC(C)(C)C)C11CCCC1 GOLXNESZZPUPJE-UHFFFAOYSA-N 0.000 description 1
- 230000009466 transformation Effects 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
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of power distribution network isolated island division methods based on minimum Custom interruption cost, this method comprises the following steps：Step 10）Build dead electricity regional distribution network simplified model after failure；Step 20）According to dead electricity region actual conditions, different classes of Subscriber Unit loss of outage function is determined；Step 30）Determine exert oneself situation and it is expected that the power off time of dead electricity moment distributed power source；Step 40）Set up to reduce isolated island partitioning model of the Custom interruption cost as target；Step 50）Feasible zone is searched for, isolated island partitioning model is solved, obtains isolated island division methods, isolated island division is carried out according to the isolated island division methods.On the premise of the various operation constraintss of power distribution network are met, this method reduction Custom interruption cost.
Description
Technical field
The invention belongs to power distribution network steady-state analysis field, it is related to a kind of power distribution network isolated island division methods, more particularly, to
A kind of power distribution network isolated island division methods based on minimum Custom interruption cost.
Background technology
Islet operation pattern is a kind of new method of operation that distribution system is introduced after distributed power generation.In this state
Under, power consumer is only by distributed power source (abbreviation DG in text) power supply.In this case, make rational planning for the connecing of distributed power source
Access point, access capacity have obvious action to reducing power distribution network loss of outage, improving system reliability.However, it is contemplated that by
Caused by failure there is hidden danger in unplanned isolated island to grid equipment safety, personal safety, while departing from distributed power source after major network
Power supply quality be difficult to ensure that, require not allowing the appearance of islet operation pattern in IEEE929-2000 standards.With distribution
The maturation of generation technology and popularization, IEEE1547-2003 standards no longer forbid the presence of islet operation pattern, encourage supply of electric power
Business realizes islet operation by technological means：All distributed power sources are out of service in short-term after failure occurs, after Fault Isolation
Using the power supply of distributed power source recovered part load.
For islet operation pattern, how to exert oneself power-balance, circuit between workload demand meeting distributed power source
On the basis of the requirement such as trend constraint, it is considered to recover user's confession while the controllability and significance level of load as much as possible
Electricity, be in recent years more popular research topic simultaneously, in terms of engineering actual angle, isolated island divides and is used as the one of fault recovery
The means of kind, are timely responded to quickly handling very necessary, thus the rapidity of isolated island splitting scheme is also that appraisal procedure is good and bad
Key factor.The purpose that isolated island is divided is reduction loss of outage, and the target of isolated island partitioning model is used as using Custom interruption cost
So that the effect that isolated island is divided is more directly perceived notable.
The content of the invention
Technical problem：The present invention provides a kind of power distribution network isolated island division methods based on minimum Custom interruption cost, full
On the premise of the sufficient various operation constraintss of power distribution network, Custom interruption cost is reduced.
Technical scheme：In order to solve the above technical problems, the embodiment of the present invention uses following technical scheme：
A kind of power distribution network isolated island division methods based on minimum Custom interruption cost, this method comprises the following steps：
Step 10) build dead electricity regional distribution network simplified model after failure；
Step 20) according to dead electricity region actual conditions, determine different classes of Subscriber Unit loss of outage function；
Step 30) determine exert oneself situation and it is expected that the power off time of dead electricity moment distributed power source；
Step 40) set up to reduce isolated island partitioning model of the Custom interruption cost as target；
Step 50) search feasible zone, isolated island partitioning model is solved, isolated island division methods are obtained, according to the isolated island division side
Method carries out isolated island division.
As preference, the step 10) include：By the transformer that may be broken down in actual power distribution network, feeder line section
And switch element is regarded as the structure node that fault rate is definite value, load is 0；It is 0, load that customer charge is regarded as into fault rate
For the load bus of definite value, power distribution network is converted into the topological structure being made up of the node comprising fault rate and information on load.
As preference, the step 20) in, according to dead electricity region all types of user history loss of outage statistics, intend
Close all types of user unit loss of outage function for the region.
As preference, the step 40) specifically include：Based on the step 10) build dead electricity regional distribution network letter
Change model, set up to reduce Custom interruption cost as target, using power, line current and node voltage as the isolated island of constraints
Partitioning model；
It is described shown as the function such as formula (1) of target to reduce Custom interruption cost：
Wherein, loss is total Custom interruption cost, and n is dead electricity Area Node number, tfFor estimated power off time, LiFor i-th
The load value of node, fk(t) it is corresponding unit loss of outage function, x when i-th of node institute on-load is kth type loadiFor
Whether i-th of node is included into isolated island, xi=1 represents that the node is not resumed, and power failure causes economic loss, xi=0 represents the section
Point is included into isolated island, and power failure does not cause economic loss to it；
The constraints is included：
Set up the power constraint as shown in formula (2)：
Wherein, PGiRepresent the active power output of i-th of distributed power source in the range of isolated island, LjFor node j active demand；
Set up the line current constraint as shown in formula (3)：
Ii.max<Ii.NFormula (3)
Wherein, Ii.maxThe maximum load current born for i-th line road or transformer；Ii.NFor i-th line road or transformer
Load current value；
Set up the node voltage constraint as shown in formula (4)：
0.95UN<Uk<1.05UNFormula (4)
Wherein, UkFor k-th of node voltage value, UNFor k-th of node load voltage value.
As preference, the step 50) in, search feasible zone includes：
Node where power supply is the center of circle in a distributed manner, along all directions of network topology, with dead electricity moment distributed power source
Exert oneself as radius, according to the constraint formulations as shown in formula (5), search is possible to be included into the node of isolated island, forms isolated island division
The feasible zone of model：
Wherein, path l between node j node g where node i to distributed power sourcei,gOn a bit, LjFor node j's
Active demand, G exerts oneself for distributed power source.
As preference, the step 50) in, solve isolated island partitioning model process as follows：
Step 501) find the dead electricity moment and exert oneself the distributed power source of minimum, determine the distribution according to the feasible zone that searches
The isolated island of formula power supply divides feasible zone scope；
Step 502) definition node weight, average weight index, accumulative load index；
Step 503) the maximum node of average weight is chosen, by the node and its to node shortest path where distributed power source
All nodes on footpath, i.e. upstream node all includes isolated island；If multiple node average weights are identical, accumulative load is chosen
Maximum node carries out above-mentioned processing；
Step 504) isolated island division target and distributed power source residual capacity are updated, update the whole for being included into isolated island node
The average weight of downstream node and accumulative information on load, by accumulative load being averaged more than the node of distributed power source residual capacity
Weight index and accumulative load index are reset；
Step 505) judge whether the accumulative load of each node in the range of isolated island division feasible zone has been cleared, such as
It is that then the distributed power source completes isolated island division work, into step 506), otherwise return to step 503)；
Step 506) raw information of all nodes for being included into isolated island is all reset, distributed power source residue is held
Amount is regarded as negative load；Then judge whether that the isolated island for completing all distributed power sources is divided, in this way, that is, complete step 506)
Isolated island partitioning model solution procedure, otherwise return to step 501) isolated island division is carried out to remaining distributed power source.
Beneficial effect：The present invention compared with prior art, with advantages below：The inventive method is directed to it is determined that position connects
Enter the power distribution network of distributed power source, the Custom interruption cost that reduction distribution network failure is caused is divided by isolated island, by power distribution network
Fail-safe analysis combines with economic evaluation.Meanwhile, the derivation algorithm time complexity of the inventive method is low, meets actual
To the rapidity demand of isolated island division work in engineering.The inventive method first divides Subscriber Unit loss of outage as isolated island
During distribution network load node weight index.The higher load of unit loss of outage is even more important, extensive in distributed power source
It should be considered first during multiple power supply.Compared with prior art, by introducing Subscriber Unit loss of outage so that distribution
Load bus in net possesses changeable weight.This make it that under different power failure situations isolated island divides recovery load bus and powered
Precedence, compare fixed weight, more tally with the actual situation.
Brief description of the drawings
Fig. 1 is the method flow schematic diagram of the embodiment of the present invention；
Fig. 2 is structural representation after the PG＆E69 Node power distribution system failures transformed；
Fig. 3 is isolated island partitioning model solution schematic flow sheet in the embodiment of the present invention；
Fig. 4 is the PG＆E69 Node power distribution system isolated island division result schematic diagrames of transformation.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, below in conjunction with drawings and Examples pair
The present invention is in depth described in detail.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in restriction invention.
As shown in figure 1, a kind of power distribution network isolated island division methods based on minimum Custom interruption cost of the present invention, following step
Suddenly：
Step 10) build dead electricity regional distribution network simplified model after failure.
Step 10) specifically include：By the transformer that may be broken down in actual power distribution network, feeder line section and switch element
It is regarded as the structure node that fault rate is definite value, load is 0；It is the load that 0, load is definite value that customer charge is regarded as into fault rate
Node, power distribution network is converted into the topological structure being made up of the node comprising fault rate and information on load.
Step 20) according to dead electricity region actual conditions, determine different classes of Subscriber Unit loss of outage function.
Step 20) in, according to dead electricity region all types of user history loss of outage statistics, fitting is each for the region
Class Subscriber Unit loss of outage function.The history loss of outage includes direct economic loss and indirect loss.
Step 30) determine exert oneself situation and it is expected that the power off time of dead electricity moment distributed power source.
Distributed power source is constrained due to the primary energy of conversion by natural conditions, present the fluctuation exerted oneself with it is random
Property；And in isolated island partition problem, it is believed that exerted oneself within a certain period of time when distributed power source is equipped with energy storage device constant, thus
It needs to be determined that the situation of exerting oneself of dead electricity moment distributed power source is as the power constraints in isolated island partitioning model solution procedure,
Meanwhile, Custom interruption cost and power failure duration are closely related, need to determine estimated power off time before isolated island partitioning model is solved.
Step 40) set up to reduce isolated island partitioning model of the Custom interruption cost as target.
Step 40) specifically include：Based on the step 10) the dead electricity regional distribution network simplified model that builds, set up to drop
Low Custom interruption cost is target, the isolated island partitioning model by constraints of power, line current and node voltage；
It is described shown as the function such as formula (1) of target to reduce Custom interruption cost：
Wherein, loss is total Custom interruption cost, and n is dead electricity Area Node number, tfFor estimated power off time, LiFor i-th
The load value of node, fk(t) it is corresponding unit loss of outage function, x when i-th of node institute on-load is kth type loadiFor
Whether i-th of node is included into isolated island, xi=1 represents that the node is not resumed, and power failure causes economic loss, xi=0 represents the section
Point is included into isolated island, and power failure does not cause economic loss to it；
The reactive requirement of load is met by local compensation, therefore power constraint considers that DG active power output is with bearing in isolated island
Active balance between lotus demand.Set up the power constraint as shown in formula (2)：
Wherein, PGiRepresent the active power output of i-th of distributed power source in the range of isolated island, LjFor node j active demand；
Set up the line current constraint as shown in formula (3)：
Ii.max<Ii.NFormula (3)
Wherein, Ii.maxThe maximum load current born for i-th line road or transformer；Ii.NFor i-th line road or transformer
Load current value；
Set up the node voltage constraint as shown in formula (4)：
0.95UN<Uk<1.05UNFormula (4)
Wherein, UkFor k-th of node voltage value, UNFor k-th of node load voltage value.
Step 50) search feasible zone, isolated island partitioning model is solved, isolated island division methods are obtained, according to the isolated island division side
Method carries out isolated island division.
The step 50) in, search feasible zone includes：
Node where power supply is the center of circle in a distributed manner, along all directions of network topology, with dead electricity moment distributed power source
Exert oneself as radius, according to the constraint formulations as shown in formula (5), search is possible to be included into the node of isolated island, forms isolated island division
The feasible zone of model：
Wherein, path l between node j node g where node i to distributed power sourcei,gOn a bit, LjFor node j's
Active demand, G exerts oneself for distributed power source.
As shown in figure 3, it is as follows to solve isolated island partitioning model process：
Step 501) find the dead electricity moment and exert oneself the distributed power source of minimum, determine the distribution according to the feasible zone that searches
The isolated island of formula power supply divides feasible zone scope；
Step 502) definition node weight, average weight index, accumulative load index；
Correspondence power off time tfKth class Subscriber Unit loss of outage functional value be i-th node weights ωi, description is such as
Under：
ωi=fk(tf)
Average weight index ΩiIt is defined as follows：
Accumulative load index XiIt is defined as follows：
Wherein, JiRepresent the set of all nodes on the shortest path of node where from i-th of node to DG (not including the
I node), LjRepresent the active demand of j-th of node；ωjRepresent the weight of j-th of node.LiRepresent the active need of node i
Ask.
Step 503) the maximum node of average weight is chosen, by the node and its to node shortest path where distributed power source
All nodes on footpath, i.e. upstream node all includes isolated island；If multiple node average weights are identical, i.e., average weight is more
Individual node gets equal maximum, then chooses the maximum node of accumulative load and carry out above-mentioned processing；
Step 504) isolated island division target and distributed power source residual capacity are updated, update the whole for being included into isolated island node
The average weight of downstream node and accumulative information on load, by accumulative load being averaged more than the node of distributed power source residual capacity
Weight index and accumulative load index are reset；
Step 505) judge whether the accumulative load of each node in the range of isolated island division feasible zone has been cleared, such as
It is that then the distributed power source completes isolated island division work, into step 506), otherwise return to step 503)；
Step 506) raw information of all nodes for being included into isolated island is all reset, distributed power source residue is held
Amount is regarded as negative load；Then judge whether that the isolated island for completing all distributed power sources is divided, in this way, that is, complete step 506)
Isolated island partitioning model solution procedure, otherwise return to step 501) isolated island division is carried out to remaining distributed power source.
The method of above-described embodiment passes through step 40), unit power failure is introduced in the object function of isolated island partitioning model and is damaged
Function is lost, as the determination foundation of load changeable weight, the purpose of reduction user's loss is realized.Meanwhile, above-mentioned implementation
The method of example passes through step 20), it is not necessary to structure is simplified to power distribution network repeatedly and carries out tree-shaped search and iteration, isolated island is improved and draws
The computational efficiency divided.
A specific embodiment is enumerated below.
As shown in Fig. 2 in the present embodiment, choosing the PG＆E69 nodes changed as distribution example.Node 5,19,
32nd, 36,52,65 access distributed power source.Failure occurs between position shown in Fig. 2, i.e. node 2 and 3, downstream after Fault Isolation
System whole dead electricity.
Simplified model of distribution networks is as follows：
It is fixed that the transformer that may be broken down in actual power distribution network, feeder line section and switch element are regarded as into fault rate
Value, load are 0 structure node, and it is the load bus that 0, load is definite value that customer charge is regarded as into fault rate, and power distribution network is turned
The topological structure being made up of the node comprising fault rate and information on load is turned to, as shown in Figure 2.
According to dead electricity region all types of user history loss of outage statistics, all types of user unit of the fitting for the region
Loss of outage function, the solution for follow-up isolated island partitioning model.The history loss of outage include direct economic loss and
Connect loss.Distributed power source is constrained due to the primary energy of conversion by natural conditions, present the fluctuation exerted oneself with it is random
Property；And in isolated island partition problem, it is believed that exerted oneself within a certain period of time when distributed power source is equipped with energy storage device constant, thus
It needs to be determined that the situation of exerting oneself of dead electricity moment distributed power source is as the power constraints in isolated island partitioning model solution procedure,
Meanwhile, Custom interruption cost and power failure duration are closely related, need to determine estimated power off time before isolated island partitioning model is solved.
In the present embodiment, have a power failure 6 DG of the moment exert oneself be followed successively by 250kW, 400kW, 40kW, 50kW, 1300kW,
100kW, the expection power off time for the moment generation that has a power failure is 1 hour.Describe for convenience, the load in the region is divided into 3 classes, recognized
It is respectively 1 $/kWh, 10 $/kWh, 100 $/kWh for the correspondence unit loss of outage of 1 hour, then corresponds to what is this time had a power failure
Load bus information is as shown in table 1.
Table 1PG＆E69 Node power distribution system load bus information (load power units：kW)
Powered to recover the big load of loss of outage as far as possible, while recovering load as much as possible, based on dead electricity area
Domain simplified model of distribution networks, sets up to reduce Custom interruption cost as target, using power, line current and node voltage as constraint
The isolated island partitioning model of condition, it is described to be described as follows using reducing Custom interruption cost as the function of target：
Power constraint is described as follows：
Line current constraint is as follows：
Ii.max<Ii.N
Node voltage constraint is as follows：
0.95UN<Uk<1.05UN
For a certain distributed power source, it is as follows that the isolated island centered on it divides feasible zone search procedure：
Node where power supply is the center of circle in a distributed manner, along all directions of network topology, with dead electricity moment distributed power source
Exert oneself as radius, search for the node for being possible to be included into isolated island according to following constraint formulations, formation isolated island partitioning model can
Row domain：
In the present embodiment, by taking the DG3 accessed at node 32 as an example, the institute for meeting above formula statement is found in all directions
There is node, then the corresponding feasible zones of DG3 are { 29,30,31,32,33,34,35 }
Isolated island partitioning model solution procedure is as follows：
(1) find the dead electricity moment to exert oneself the distributed power source of minimum, the distributed electrical is determined according to the feasible zone that searches
The isolated island in source divides feasible zone scope.In the present embodiment, isolated island division is carried out to DG3 first；
(2) definition node weight, average weight index, accumulative load index：
In the present embodiment, 9*8 isolated island Matrix dividing I record cast solution procedurees are set up, matrix is as follows：
Wherein, l is that isolated island divides target, is initially ∑ Li`fk(tf)；S is DG residual capacity, is initially fault moment
DG's exerts oneself；2nd~8 row the 1st, 2 rows record the average weight index and accumulative load of 7 alternative load buses and refer to respectively
Mark；The identifier that whether is included into of 3rd~9 row the 1st row 7 alternative load buses of record, be initially 1 (1 expression is not included into,
It can still cause damage；0 represents to be included into, and does not result in loss；0) identifier of node is where giving tacit consent to DG；Remaining 7*7 matrixes
Topological relation between area identification load bus, RmjIdentify the annexation between node m and node j, Rmj=0 represents node m
It is indirect between node j to be connected or do not connect；The angle of node, R where from DGmj=1 represents that node m saves for node j father
Point, Rmj=-1 represents the child node that node m is node j.The then initial isolated island Matrix dividing I centered on DG3 of the present embodiment
It is as follows：
(3) the maximum node of average weight is chosen, by the node and its on node shortest path where distributed power source
All nodes, i.e. upstream node all includes isolated island, if multiple node average weights are identical, chooses accumulative load maximum
Node carry out above-mentioned processing.For the present embodiment, node 35 first is selected, is recalled successively according to node annexation and found
Its upstream node 33,34；
(4) update isolated island and divide target and distributed power source residual capacity, will be included into isolated island whole nodes it is flat
Equal weight index and accumulative load index are reset, and are updated the average weight for the whole downstream nodes for being included into isolated island node and are added up
Information on load, it is as follows：
Wherein, l is that isolated island divides target, is initially ∑ Li`fk(tf)；L ' divides target, Ω for the isolated island after updatingkFor section
The corresponding average weight indexs of point k, XkFor the corresponding accumulative load indexs of node k；S is DG residual capacity, is initially failure
Moment DG's exerts oneself；S ' is the DG residual capacities after updating；ΩupFor the node k corresponding average weight index of each node in upstream,
Ω’upFor the Ω after renewalup, XupFor the node k corresponding accumulative load index of each node in upstream, X 'upFor the X after renewalup,
ΩdownFor the node k corresponding average weight index of each node in downstream, Ω 'downFor the Ω after renewaldown, XdownFor node k's
The corresponding accumulative load index of each node in downstream；X’downX after respectively updatingdown。
It is important to note that for selected node k, its father node i is unique, and node i may be also
In the presence of other child nodes, it is also required to update average weight and accumulative information on load for these branch nodes.
Meanwhile, the average weight index and accumulative load that accumulative load is more than to the node of distributed power source residual capacity refer to
SD zero, is not considered further that.In the present embodiment, for DG3 first time isolated island divide after, isolated island Matrix dividing I be updated to as
Shown in lower：
(5) judge whether that the accumulative load of each node has been cleared, in this way, then the distributed power source completes isolated island stroke
The division of labor is made, into step (6), otherwise return to step (3)；
(6) raw information of all nodes for being included into isolated island is all reset, prevents to repeat when isolated island crosses
Calculate, influence isolated island division result；Distributed power source residual capacity is regarded as negative load, after judge whether to complete and be distributed
The isolated island of formula power supply is divided, and in this way, that is, completes the solution procedure of whole isolated island partitioning model, otherwise return to step (1) is to residue
Distributed power source carry out isolated island division.In the present embodiment, node 32, then be considered as by residue 0.5kW surplus after DG3 is divided
Load be -0.5kW, unit loss of outage be 10 $/kWh node, hereafter successively using DG4, DG6, DG1, DG2 and DG5 as
Center carries out isolated island division work, and here is omitted.
After initial isolated island division terminates, initial isolated island is verified according to isolated island partitioning model constraints, obtained most
Whole isolated island division result, as shown in Figure 4.It is $ 13888.9 that this isolated island, which divides obtained loss of outage, compares and does not carry out orphan
There being $ 70502 loss of outage when island is divided is greatly reduced.Meanwhile, isolated island divides target and is intuitively demonstrated by economic index,
The quality of isolated island splitting scheme is set to be compared at a glance.
Above-described embodiment is only the preferred embodiment of the present invention, it should be pointed out that：For the ordinary skill of the art
For personnel, under the premise without departing from the principles of the invention, some improvement and equivalent substitution can also be made, these are to the present invention
Claim be improved with the technical scheme after equivalent substitution, each fall within protection scope of the present invention.
Claims (6)
1. a kind of power distribution network isolated island division methods based on minimum Custom interruption cost, it is characterised in that this method includes as follows
Step：
Step 10) build dead electricity regional distribution network simplified model after failure；
Step 20) according to dead electricity region actual conditions, determine different classes of Subscriber Unit loss of outage function；
Step 30) determine exert oneself situation and it is expected that the power off time of dead electricity moment distributed power source；
Step 40) set up to reduce isolated island partitioning model of the Custom interruption cost as target；
Step 50) search feasible zone, isolated island partitioning model is solved, isolated island division methods is obtained, enters according to the isolated island division methods
Row isolated island is divided.
2. the power distribution network isolated island division methods according to claim 1 based on minimum Custom interruption cost, it is characterised in that
The step 10) include：The transformer that may be broken down in actual power distribution network, feeder line section and switch element are regarded as failure
The structure node that rate is definite value, load is 0；It is the load bus that 0, load is definite value that customer charge is regarded as into fault rate, will
Power distribution network is converted into the topological structure being made up of the node comprising fault rate and information on load.
3. the power distribution network isolated island division methods according to claim 1 or 2 based on minimum Custom interruption cost, its feature exists
In the step 20) in, according to dead electricity region all types of user history loss of outage statistics, fitting is each for the region
Class Subscriber Unit loss of outage function.
4. the power distribution network isolated island division methods according to claim 1 or 2 based on minimum Custom interruption cost, its feature exists
In the step 40) specifically include：Based on the step 10) the dead electricity regional distribution network simplified model that builds, set up to reduce
Custom interruption cost is target, the isolated island partitioning model by constraints of power, line current and node voltage；
It is described shown as the function such as formula (1) of target to reduce Custom interruption cost：
Wherein, loss is total Custom interruption cost, and n is dead electricity Area Node number, tfFor estimated power off time, LiFor i-th of node
Load value, fk(t) it is corresponding unit loss of outage function, x when i-th of node institute on-load is kth type loadiFor i-th
Whether node is included into isolated island, xi=1 represents that the node is not resumed, and power failure causes economic loss, xi=0 represents the node quilt
Isolated island is included, power failure does not cause economic loss to it；
The constraints is included：
Set up the power constraint as shown in formula (2)：
Wherein, PGiRepresent the active power output of i-th of distributed power source in the range of isolated island, LjFor node j active demand；
Set up the line current constraint as shown in formula (3)：
Ii.max<Ii.NFormula (3)
Wherein, Ii.maxThe maximum load current born for i-th line road or transformer；Ii.NFor i-th line road or the volume of transformer
Determine current value；
Set up the node voltage constraint as shown in formula (4)：
0.95UN<Uk<1.05UNFormula (4)
Wherein, UkFor k-th of node voltage value, UNFor k-th of node load voltage value.
5. the power distribution network isolated island division methods according to claim 1 or 2 based on minimum Custom interruption cost, its feature exists
In the step 50) in, search feasible zone includes：
Node where power supply is the center of circle in a distributed manner, along all directions of network topology, is exerted oneself with dead electricity moment distributed power source
For radius, according to the constraint formulations as shown in formula (5), search is possible to be included into the node of isolated island, forms isolated island partitioning model
Feasible zone：
Wherein, path l between node j node g where node i to distributed power sourcei,gOn a bit, LjFor the active of node j
Demand, G exerts oneself for distributed power source.
6. the power distribution network isolated island division methods according to claim 1 or 2 based on minimum Custom interruption cost, its feature exists
In the step 50) in, solve isolated island partitioning model process as follows：
Step 501) find the dead electricity moment and exert oneself the distributed power source of minimum, determine the distributed electrical according to the feasible zone that searches
The isolated island in source divides feasible zone scope；
Step 502) definition node weight, average weight index, accumulative load index；
Step 503) the maximum node of average weight is chosen, by the node and its on node shortest path where distributed power source
All nodes, i.e. upstream node all includes isolated island；If multiple node average weights are identical, accumulative load is chosen maximum
Node carry out above-mentioned processing；
Step 504) isolated island division target and distributed power source residual capacity are updated, update the whole downstreams for being included into isolated island node
The average weight of node and accumulative information on load, accumulative load are more than the average weight of the node of distributed power source residual capacity
Index and accumulative load index are reset；
Step 505) judge whether the accumulative load of each node in the range of isolated island division feasible zone has been cleared, in this way,
Then the distributed power source completes isolated island division work, into step 506), otherwise return to step 503)；
Step 506) raw information of all nodes for being included into isolated island is all reset, distributed power source residual capacity is regarded
Make negative load；Then the isolated island for judging whether to complete all distributed power sources is divided, in this way, that is, completes step 506) isolated island
Partitioning model solution procedure, otherwise return to step 501) isolated island division is carried out to remaining distributed power source.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710144090.5A CN107069801B (en) | 2017-03-13 | 2017-03-13 | A kind of power distribution network isolated island division methods based on minimum Custom interruption cost |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710144090.5A CN107069801B (en) | 2017-03-13 | 2017-03-13 | A kind of power distribution network isolated island division methods based on minimum Custom interruption cost |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107069801A true CN107069801A (en) | 2017-08-18 |
| CN107069801B CN107069801B (en) | 2019-09-24 |
Family
ID=59621688
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710144090.5A Expired - Fee Related CN107069801B (en) | 2017-03-13 | 2017-03-13 | A kind of power distribution network isolated island division methods based on minimum Custom interruption cost |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107069801B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107611963A (en) * | 2017-09-01 | 2018-01-19 | 东南大学 | A kind of PMU Optimal Configuration Methods for considering loss of outage coefficient |
| CN110518572A (en) * | 2019-07-12 | 2019-11-29 | 国网浙江省电力有限公司杭州供电公司 | A kind of power distribution network isolated island division methods based on minimum load loss |
| CN111626603A (en) * | 2020-05-26 | 2020-09-04 | 中国石油化工股份有限公司 | Oil field ordered power utilization system with genetic algorithm and ordered power utilization distribution method thereof |
| CN112465235A (en) * | 2020-12-01 | 2021-03-09 | 国网浙江杭州市富阳区供电有限公司 | Power failure interval prediction method for reducing electric quantity loss |
| CN112989531A (en) * | 2021-02-22 | 2021-06-18 | 国网安徽省电力有限公司合肥供电公司 | Method and system for automatically generating platform-to-platform project planning scheme |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101986496B (en) * | 2010-09-29 | 2012-08-08 | 天津大学 | Method for forming islands of distribution system containing distributed power supply |
| CN103078391A (en) * | 2013-01-10 | 2013-05-01 | 长兴县供电局 | Power distribution network power supply power restoration method based on photovoltaic power generation system |
| CN105098762A (en) * | 2014-05-09 | 2015-11-25 | 国家电网公司 | Island dividing method for power distribution network having distributed power sources |
| CN106340899A (en) * | 2016-09-21 | 2017-01-18 | 东南大学 | Dynamic islanding method for reducing user outage loss target |
-
2017
- 2017-03-13 CN CN201710144090.5A patent/CN107069801B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101986496B (en) * | 2010-09-29 | 2012-08-08 | 天津大学 | Method for forming islands of distribution system containing distributed power supply |
| CN103078391A (en) * | 2013-01-10 | 2013-05-01 | 长兴县供电局 | Power distribution network power supply power restoration method based on photovoltaic power generation system |
| CN105098762A (en) * | 2014-05-09 | 2015-11-25 | 国家电网公司 | Island dividing method for power distribution network having distributed power sources |
| CN106340899A (en) * | 2016-09-21 | 2017-01-18 | 东南大学 | Dynamic islanding method for reducing user outage loss target |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107611963A (en) * | 2017-09-01 | 2018-01-19 | 东南大学 | A kind of PMU Optimal Configuration Methods for considering loss of outage coefficient |
| CN107611963B (en) * | 2017-09-01 | 2019-08-20 | 东南大学 | A kind of PMU Optimal Configuration Method considering loss of outage coefficient |
| CN110518572A (en) * | 2019-07-12 | 2019-11-29 | 国网浙江省电力有限公司杭州供电公司 | A kind of power distribution network isolated island division methods based on minimum load loss |
| CN111626603A (en) * | 2020-05-26 | 2020-09-04 | 中国石油化工股份有限公司 | Oil field ordered power utilization system with genetic algorithm and ordered power utilization distribution method thereof |
| CN111626603B (en) * | 2020-05-26 | 2023-08-04 | 中国石油化工股份有限公司 | Oil field ordered electricity utilization system with genetic algorithm and ordered electricity utilization distribution method thereof |
| CN112465235A (en) * | 2020-12-01 | 2021-03-09 | 国网浙江杭州市富阳区供电有限公司 | Power failure interval prediction method for reducing electric quantity loss |
| CN112989531A (en) * | 2021-02-22 | 2021-06-18 | 国网安徽省电力有限公司合肥供电公司 | Method and system for automatically generating platform-to-platform project planning scheme |
| CN112989531B (en) * | 2021-02-22 | 2022-09-23 | 国网安徽省电力有限公司合肥供电公司 | Method and system for automatically generating platform-to-platform project planning scheme |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107069801B (en) | 2019-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107069801A (en) | A kind of power distribution network isolated island division methods based on minimum Custom interruption cost | |
| CN102983629B (en) | Auxiliary decision-making method for on-line power system restoration | |
| Zhou et al. | Distributed secondary control for islanded microgrids with mobile emergency resources | |
| CN104820865B (en) | Intelligent distribution network fault recovery intelligent optimization method based on graph theory | |
| CN104537580B (en) | The construction method of distribution network structure based on K means clusters | |
| CN107017622A (en) | The multiple faults multiple target of distribution containing DG combined optimization repairing recovery policy is asked for | |
| CN105762790A (en) | Evaluation index algorithm for AC and DC distribution network power supply mode selection | |
| CN102801790A (en) | Microgrid energy management system based on IEC (International Electrotechnical Commission) 61850 and multi-agent method | |
| CN110350511B (en) | Power distribution network fault recovery method based on transfer power affinity | |
| CN109598377B (en) | AC/DC hybrid power distribution network robust planning method based on fault constraint | |
| CN107622360A (en) | A kind of critical circuits recognition methods for considering subjective and objective factor | |
| Zand et al. | Sensitivity analysis index to determine the optimal location of multi-objective UPFC for improvement of power quality parameters | |
| CN111079972A (en) | Method, device and medium for planning reliability of active power distribution network | |
| CN106655253A (en) | Single-/three-phase multi-micro-grid region dynamic partitioning method | |
| CN105958472B (en) | A kind of power distribution network operation optimizing system | |
| CN103106561B (en) | Power system load recovers aid decision-making system and decision making method thereof | |
| CN107017635A (en) | A kind of adaptive islanding control method and system | |
| CN108092267A (en) | A kind of power distribution network access planning system and method based on intelligent body | |
| CN107294086B (en) | The service restoration method realized based on network equivalent and parallelization | |
| CN110472867A (en) | Power distribution network Network adaptation assessment processing method and device | |
| CN115549075A (en) | Power supply recovery method and system for power distribution network containing micro-grid | |
| Huang et al. | Multi-objective Day-ahead Polarity Switching Optimization for Voltage Unbalance in Bipolar DC Distribution Networks | |
| CN105576658A (en) | Method for analyzing economy of power distribution network reconstruction project on basis of wiring mode | |
| CN111092453A (en) | Power grid key line identification method for multipoint access photovoltaic power supply | |
| CN109066675A (en) | A kind of topological arrangement method for power distribution network wide area self-healing system |
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 | ||
| CB02 | Change of applicant information |
Address after: No. 15, Qinghe Xiaoying East Road, Haidian District, Beijing Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant after: STATE GRID CORPORATION OF CHINA Applicant after: SOUTHEAST University Address before: No. 15, Qinghe Xiaoying East Road, Haidian District, Beijing Applicant before: China Electric Power Research Institute Applicant before: State Grid Corporation of China Applicant before: Southeast University |
|
| CB02 | Change of applicant information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180905 Address after: No. 15, Qinghe Xiaoying East Road, Haidian District, Beijing Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant after: STATE GRID CORPORATION OF CHINA Applicant after: SOUTHEAST University Applicant after: STATE GRID ZHEJIANG ELECTRIC POWER Co.,Ltd. Address before: No. 15, Qinghe Xiaoying East Road, Haidian District, Beijing Applicant before: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant before: State Grid Corporation of China Applicant before: Southeast University |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190924 Termination date: 20200313 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |