CN103971299A - Power distribution network islanding method based on harmonic algorithm - Google Patents

Power distribution network islanding method based on harmonic algorithm Download PDF

Info

Publication number
CN103971299A
CN103971299A CN201410224887.2A CN201410224887A CN103971299A CN 103971299 A CN103971299 A CN 103971299A CN 201410224887 A CN201410224887 A CN 201410224887A CN 103971299 A CN103971299 A CN 103971299A
Authority
CN
China
Prior art keywords
isolated island
distribution network
load
power distribution
node
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
Application number
CN201410224887.2A
Other languages
Chinese (zh)
Other versions
CN103971299B (en
Inventor
汪沨
曾业运
陈春
曹一家
黄纯
董旭柱
黄小耘
许爱东
于力
彭飞进
雷金勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan University
CSG Electric Power Research Institute
Foshan Power Supply Bureau of Guangdong Power Grid Corp
Research Institute of Southern Power Grid Co Ltd
Original Assignee
Hunan University
Foshan Power Supply Bureau of Guangdong Power Grid Corp
Research Institute of Southern Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan University, Foshan Power Supply Bureau of Guangdong Power Grid Corp, Research Institute of Southern Power Grid Co Ltd filed Critical Hunan University
Priority to CN201410224887.2A priority Critical patent/CN103971299B/en
Publication of CN103971299A publication Critical patent/CN103971299A/en
Application granted granted Critical
Publication of CN103971299B publication Critical patent/CN103971299B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention relates to a power distribution network islanding method based on a harmonic algorithm. Aiming for largest recovered load switching value and smallest motion switching value, and considering constraint conditions such as active power balance and non-overload of line current in an island, the power distribution network islanding method based on the harmonic algorithm provided by the invention carries out 0-1 encoding for each node in the power distribution network, obtains a specific islanding scheme in combination with a breadth-first search algorithm, calculates an objective function value of each island, adjusts, evaluates and updates the islanding scheme by utilizing the harmonic algorithm, then carries out trend inspection on the current optimal solution, and finally gives a better islanding scheme. The method for realizing full search of a solution space to solve the power distribution network islanding problem based on the harmonic algorithm provided by the invention can utilize individual local information and group global information guidance algorithm for further search, is simple in principle and few in adjustable parameters, and can obtain the optimal islanding scheme in a relatively short period of time when the islanding problem is solved.

Description

A kind of power distribution network isolated island division methods based on harmony algorithm
Technical field
The present invention relates to a kind of power distribution network isolated island analytical approach, relate in particular to a kind of power distribution network isolated island division methods based on harmony algorithm.
Background technology
Power distribution network islet operation is the one abnormal method of operation of power distribution network after access distributed power source (DG).The division of power distribution network isolated island is exactly that higher level's electrical network of working as power distribution network outside breaks down, voltage to frequency is out-of-limit and vibrates in the situations such as step-out, according to position, the generating capacity etc. of electric network composition, DG, power distribution network off-the-line is become to one or several isolated subnet, to ensure the particularly power supply of critical load of loading, until fault is got rid of, system restores electricity.Under this method of operation, the certain customers of distribution system are only powered by DG, and the system that comprises this part user and DG is called as the isolated island in distribution system.
According to the relevant regulations of IEEE std.999-2000, after distribution system is stopped transport because of fault, for maintainer's personal safety, simultaneously in order to make power distribution network eliminate as early as possible fault, all DG in power distribution network must be out of service in short-term, forbids the generation of unplanned islet operation.
But because DG has independently-powered ability, after fault isolation, can utilize the power supply of DG recovered part important load and relevant load, thereby significantly improve the power supply reliability of system.To this, new standard IEEE std.1547-2003 encourages conscious islet operation.The in the situation that of system dead electricity, if can independently-powered isolated island can ensure the uninterrupted of important load power supply as power supply forms taking all kinds of DG, allow load as much as possible obtain service restoration simultaneously, can significantly promote the power supply reliability of distribution system.
How fully to effectively utilize limited distribution type electric energy, determine the best electric power isolated island taking DG as power center, and in time carry out islet operation and the synchronous operation with large electrical network, thereby bring into play to greatest extent the advantage of distributed power generation, ensure that the continued power of important load and many as far as possible loads is the key issues that realize the efficient operation of power distribution network isolated island.
Chinese patent " a kind of for consider interconnection switch containing distributed power generation distribution isolated island division methods ", patent publication No.: CN102983570A, disclose a kind of consider interconnection switch containing distributed power generation distribution isolated island division methods, based on distribution system graph model, obtain dividing new model containing DG distribution system isolated island; Build the fault graph model G of dead electricity region distribution system afterwards.For the distributed power generation power distribution network isolated island division methods that contains that builds the consideration interconnection switch that meets engineering demand provides theoretical foundation, propose the solution strategies containing the power distribution network isolated island partition problem of ring network structure taking this algorithm as basis, expanded the through engineering approaches range of application containing distributed power generation isolated island partitioning technology.
But it is not good that said method is processed for discrete optimization problem, and computing velocity is slow, is easily absorbed in local optimum, and in the time solving isolated island partition problem, the time is longer.
In view of above-mentioned defect, creator of the present invention has obtained this creation finally through long research and practice.
Summary of the invention
The object of the present invention is to provide a kind of power distribution network isolated island division methods based on harmony algorithm (HS algorithm, Harmony search), in order to overcome above-mentioned technological deficiency.
For achieving the above object, the invention provides a kind of power distribution network isolated island division methods based on harmony algorithm, it is characterized in that, its process is:
Step 1, distribution network system arranges the initialization of harmony algorithm parameter;
Step 2, distribution network system is to each dimension of each individuality in power distribution network initialization of encoding, produce N initial individual;
Step 3, distribution network system calculates the target function value of each individuality in initial harmony data base HM;
Step 4, distribution network system produces new explanation according to HS algorithmic rule, and calculates the target function value of new explanation;
Step 5, distribution network system upgrades harmony data base HM; Whether the target function value that judges the new explanation that above-mentioned steps 4 calculates is greater than in described harmony data base HM the target function value of poor solution, if be greater than, replaces the poorest solution in described harmony data base HM with new explanation;
Step 6, distribution network system judges whether to meet the end condition of harmony algorithm: actual iterations is less than the iterations arranging in above-mentioned steps 1; If meet, carry out following step 7; If do not meet, return to step 4, and carry out successively;
Step 7, distribution network system is obtaining after current optimum solution, and the each isolated island in isolated island scheme is carried out to trend verification;
Step 8, distribution network system output isolated island scheme, obtains splitting point.
Preferably, in above-mentioned steps 2, for each individuality, the cataloged procedure of its every one dimension is: for the node that is connected with DG, coding is initialized as 1; For the node that is connected with one-level or two stage loads, coding is initialized as 1; For the node that is connected with three stage loads, coding random initializtion is 0 or 1.
Preferably, in above-mentioned steps 3, first with breadth-first search, individuality is searched for, determine load and distributed power source DG that in splitting scheme, each isolated island comprises, judge whether each isolated island meets the requirement of constraint condition A and condition B, wherein, judge according to the following formula whether each isolated island meets constraint condition A and B
Wherein, Sump ifor the load restoration amount of single isolated island, m is the load number that this isolated island comprises, w jfor load weight, P jfor the active power of load; The basic ideas that above-mentioned constraint condition is judged as: in each isolated island scheme, may comprise n isolated island separated from one another, to single isolated island, judge whether this isolated island meets constraint condition A and B; If both meet simultaneously, calculate its load restoration amount, load restoration amount is the product of load active power and load level weights; If have at least one not meet, its load restoration amount set to 0.
Preferably, in above-mentioned steps 3, be calculated as follows the target function value of each individuality, the target function value that isolated island is divided is exactly the switch number that the load restoration amount sum of each isolated island in isolated island splitting scheme deducts action again,
max ( λ 1 Σ i = 1 n Sump i - λ 2 N breaker )
In formula, λ 1and λ 2for coefficient; Sump ifor the load restoration amount of single isolated island; N breakerthe number of switches that need to move for dividing isolated island.
Preferably, in above-mentioned steps 3, constraint condition A is: active power balance in isolated island;
The active power that in isolated island, generator sends must be greater than the active power of isolated island internal loading demand: is total to as shown in formula as following,
Σ G i ∈ V P G i > Σ L j ∈ V P L j
In formula: for distributed power source generated energy; for the load in isolated island.
Preferably, in above-mentioned steps 4, the coding value of each in definite this solution of a new explanation generation need; Wherein, in new explanation, the coding of a certain position x specifically determines that mode is:
Step 41, if x corresponding node is distributed power source DG, x indirect assignment is 1;
Step 42, if x corresponding node is important load, x assignment is 1; Then produce a random number rand between [0,1] 2if, rand 2< PAR, x accepts disturbance and becomes 0; Otherwise, remaining 1 constant, PAR is parameter preset, it is worth between 0-1;
Step 43, if x corresponding node is general load, produces a random number rand between [0,1] 1if, rand 1< HMCR, the value of x equals the coding of the x position of any one individuality in harmony data base HM, and HMCR is parameter preset, and it is worth between 0-1; Otherwise x random coded is 0 or 1; , then produce a random number rand between [0,1] then 2if, rand 2< PAR, x accepts disturbance change coding, and from being 1 to become 0 or be that 0 to become 1, PAR be parameter preset originally originally, it is worth between 0-1; Otherwise the value of x remains unchanged.
Preferably, in above-mentioned steps 7, the detailed process of each isolated island being carried out to trend verification is:
Step 71, selecting to have firm power in isolated island exports and has the distributed power source of certain frequency modulation and voltage modulation ability as the balance node of isolated island, other distributed power source and load, as PQ node, carry out trend calculating by Newton-Raphson method, and the trend that obtains whole isolated island distributes;
Step 72, according to calculation of tidal current, judges whether each branch road meets constraint condition C and D; If do not met, excise the interior current load level of isolated island minimum, the load of meritorious demand maximum, re-starts trend by above-mentioned steps 71 and calculates, until isolated island meets institute's Prescribed Properties.
Preferably, in above-mentioned steps 7, carry out trend checking according to constraint condition C and D, above-mentioned constraint condition C is: line current nonoverload; According to following formula,
I ij < I ij max
In formula: I ijfor flowing through the actual current of circuit; for flowing through the maximum allowed current of circuit.Preferably, above-mentioned constraint condition D is: busbar voltage is not out-of-limit; According to following formula,
U min<U i<U max
In formula: U ifor the virtual voltage of bus; U minand U maxfor the upper limit, the lower limit of busbar voltage permission.
Preferably, before this step 1 isolated island is divided, should determine the power distribution network region separating from large electrical network after fault isolation, and the scope of dividing using this as isolated island, power distribution network practical topology figure is reduced to a little and sideband has the non-directed graph of weights; Bus nodes is composed to the weights with the power of its connected load or distributed power source DG, then load and distributed power source DG node are deleted, power distribution network is further simplified to the node-branch road non-directed graph of node with weights.
Beneficial effect of the present invention is compared with prior art: the method that the present invention is based on harmony algorithm and realize the thorough search of power distribution network isolated island partition problem solution space, can utilize individual local message and colony's global information to instruct algorithm further to search for, and its principle is simple, adjustable parameter is few, easily realizes.Harmony algorithm can obtain in the short period of time optimum isolated island splitting scheme in the time solving isolated island partition problem.
In objective function of the present invention, one of target that switch number is also divided as isolated island.Islet operation is a kind of abnormal method of operation, and the time is shorter, and therefore, isolated island is divided and should be taken into full account after major network fault restoration, the isolated island after division and have a power failure load can rapid combination, obtain stable power supply; Step switch number is fewer, is more conducive to the service restoration after fault restoration.
Because being isolated island, important load divides the object of answering priority restores, according to the preferential power supply principle of important load, its coding is initialized as to 1 to be conducive to just find a preferably initial solution at the algorithm initial stage, reduce algorithm in early days because of the blind search number of times of individual random initializtion, thereby greatly save computing time.In algorithm search process, DG nodes encoding is not accepted disturbance, is always 1, and important load nodes encoding becomes 0 to accept disturbance compared with small probability, can avoid the generation of a lot of inferior solutions.
Brief description of the drawings
Fig. 1 is the power distribution network isolated island division methods that the present invention is based on harmony algorithm;
Fig. 2 is isolated island division result figure in the embodiment of the present invention.
Embodiment
Below in conjunction with accompanying drawing, technical characterictic and the advantage with other above-mentioned to the present invention are described in more detail.
The present invention is based on the power distribution network isolated island division methods of harmony algorithm, taking the on-load switch amount recovered at most and step switch minimum as target, consider the constraint conditions such as the interior active power balance of isolated island, line current nonoverload, each node in distribution network is carried out to 0-1 coding, obtain concrete isolated island splitting scheme in conjunction with breadth-first search, and calculate the target function value of each isolated island, utilize harmony algorithm that isolated island scheme is adjusted, evaluates, upgraded, then current optimum solution is carried out to trend inspection, finally provide preferably isolated island splitting scheme.
Refer to shown in Fig. 1, it is the process flow diagram that the present invention is based on the power distribution network isolated island division methods of harmony algorithm, and this detailed process is:
Step 1, distribution network system arranges the initialization of harmony algorithm parameter;
Harmony algorithm in the present invention comprises harmony data base size (HMS), the dimension M separating, harmony data base consideration probability (HMCR), fine setting probability (PAR) and iterations (Tmax).
In harmony algorithm, the parameter of optimization problem comprises objective function f (x), variable xi and gathers number N and each variable-value scope of x, variable; The parameter of HS algorithm comprises harmony data base size (HMS), the dimension M separating, harmony data base consideration probability (HMCR), fine setting probability (PAR), iterations (Tmax) and end condition.
Before this step 1 isolated island is divided, should determine the power distribution network region separating from large electrical network after fault isolation, and the scope of dividing using this as isolated island, power distribution network practical topology figure is reduced to a little and sideband has the non-directed graph of weights.Power distribution network is closed loop design, open loop operation, conventionally adopts radial structure.Therefore, power distribution network can be regarded as taking system side power supply as root node, taking the bus of connected load as interior point, and to load as leaf node, the tree taking the branch road of connection bus and bus, connection bus and load as branch; DG, as the power supply connecting system disperseing, can regard as and increase the leaf node with power supply natures.If, bus nodes is composed to the weights with the power of its connected load or DG, then load and DG node are deleted, power distribution network further can be simplified to the node-branch road non-directed graph of node with weights.
Step 2, distribution network system is to each dimension of each individuality in power distribution network initialization of encoding, produce N initial individual;
Individual dimension is exactly the node number in network, and the each node in distribution network is 0-1 coding (1 represents that this node is in isolated island, and 0 represents that this node is not in isolated island).
For each individuality, the coding initialization procedure of its every one dimension is: for the node that is connected with DG, coding is initialized as 1; For the node that is connected with one-level or two stage loads, coding is initialized as 1; For the node that is connected with three stage loads, coding random initializtion is 0 or 1.
Because being isolated island, important load divides the object of answering priority restores, according to the preferential power supply principle of important load, its coding is initialized as to 1 to be conducive to just find a preferably initial solution at the algorithm initial stage, reduce algorithm in early days because of the blind search number of times of individual random initializtion, thereby greatly save computing time.In algorithm search process, DG nodes encoding is not accepted disturbance, is always 1, and important load nodes encoding becomes 0 to accept disturbance compared with small probability, can avoid the generation of a lot of inferior solutions.
Step 3, distribution network system calculates the target function value of each individuality in initial harmony data base HM;
With breadth-first search, individuality is searched for, determined load and DG that in splitting scheme, each isolated island comprises, judge whether each isolated island meets the requirement of 4 constraint condition conditional A and condition B;
Judge by following formula (1) whether each isolated island meets constraint condition A and B,
Wherein, m is the load number that this isolated island comprises, and wj is load weight (load level is higher, and respective weights is larger), P jfor the active power of load.
The basic ideas that the present invention judges above-mentioned constraint condition as: in each isolated island scheme, may comprise n isolated island separated from one another, to single isolated island, judge whether this isolated island meets constraint condition A and B; If both meet simultaneously, calculate its load restoration amount, load restoration amount is the product of load active power and load level weights; If have at least one not meet, its load restoration amount set to 0.
Calculate the load restoration amount of each isolated island, so back-pushed-type (2) calculates this individual target function value.The target function value that isolated island is divided is exactly the switch number that the load restoration amount sum of each isolated island in isolated island splitting scheme deducts action again.
max ( &lambda; 1 &Sigma; i = 1 n Sump i - &lambda; 2 N breaker ) - - - ( 2 )
In formula, λ 1and λ 2for coefficient; Sump ifor the load restoration amount of single isolated island; N breakerthe number of switches that need to move for dividing isolated island.
In objective function of the present invention, one of target that switch number is also divided as isolated island.Islet operation is a kind of abnormal method of operation, and the time is shorter, and therefore, isolated island is divided and should be taken into full account after major network fault restoration, the isolated island after division and have a power failure load can rapid combination, obtain stable power supply; Step switch number is fewer, is more conducive to the service restoration after fault restoration.
Above-mentioned constraint condition A and B are respectively in the present invention:
Active power balance in A, isolated island
In system, reactive power is all generally to compensate in this locality, therefore, when isolated island is divided, only considers active power balance.The active power that in isolated island, generator sends must be greater than the active power of isolated island internal loading demand: as shown in following formula (3),
&Sigma; G i &Element; V P G i > &Sigma; L j &Element; V P L j - - - ( 3 )
In formula: for distributed power source generated energy; for the load in isolated island, G i∈ V represents that this constraint condition only considers generator in isolated island, L j∈ V represents that this constraint condition only considers isolated island internal loading.
In B, isolated island, at least contain the distributed power source of a firm power output.
Step 4, distribution network system produces new explanation according to HS algorithmic rule; Calculate the target function value of new explanation by above-mentioned formula (1) and formula (2);
Producing in the process of new explanation, arrange DG nodes encoding and do not accept disturbance, be always 1, important load nodes encoding becomes 0 to accept disturbance compared with small probability, suppresses the generation of inferior solution.
The generation of new explanation: for this isolated island division methods, produce a new explanation, only need to determine the coding (value 0 or 1) of each (M position altogether) in this solution.Wherein, in new explanation, the coding of a certain position x specifically determines that mode is:
Step 41, if x corresponding node is DG, x indirect assignment is 1;
Step 42, if x corresponding node is important load, x assignment is 1; Then produce a random number rand between [0,1] 2if, rand 2< PAR, x accepts disturbance and becomes 0; Otherwise, remain 1 constant; PAR is parameter preset, and it is worth between 0-1; Step 43, if x corresponding node is general load, produces a random number rand between [0,1] 1if, rand 1< HMCR, to equal the coding HMCR of the x position of any one individuality in harmony data base HM be parameter preset to the value of x, it is worth between 0-1; Otherwise x random coded is 0 or 1; , then produce a random number rand between [0,1] then 2if, rand 2< PAR, x accepts disturbance change coding, from being 1 to become 0 or be 0 to become 1 originally originally; Otherwise the value of x remains unchanged.
Step 5, distribution network system upgrades harmony data base HM; Whether the target function value that judges the new explanation that above-mentioned steps 4 calculates is greater than in described harmony data base HM the target function value of poor solution, if be greater than, replaces the poorest solution in described harmony data base HM with new explanation.
Step 6, distribution network system judges whether to meet the end condition of harmony algorithm: actual iterations is less than the iterations arranging in above-mentioned steps 1; If meet, carry out following step 7; If do not meet, return to step 4, and carry out successively.
Step 7, distribution network system is obtaining after current optimum solution, and the each isolated island in isolated island scheme is carried out to trend verification.Verify according to constraint condition C and D, if do not meet constraint condition C and D, excise successively isolated island internal loading grade minimum, the load of the meritorious demand maximum of load, until meet institute's Prescribed Properties.
Above-mentioned constraint condition C is: line current nonoverload; According to following formula (4),
I ij < I ij max - - - ( 4 )
In formula: I ijfor flowing through the actual current of circuit; for flowing through the maximum allowed current of circuit.Above-mentioned constraint condition D is: busbar voltage is not out-of-limit, according to following formula (5),
U min<U i<Umax (5)
In formula: U ifor the virtual voltage of bus; U minand U maxfor the upper limit, the lower limit of busbar voltage permission.
In this step, the detailed process of each isolated island being carried out to trend verification is:
Step 71, selecting to have firm power in isolated island exports and has the distributed power source of certain frequency modulation and voltage modulation ability as the balance node of isolated island, other distributed power source and load are as PQ node (node that active power and reactive power are fixing), carry out trend calculating by Newton-Raphson method, the trend that obtains whole isolated island distributes.
In said process, have multiplely if qualified, select the DG of fm capacity maximum.
The trend of described isolated island distributes, comprise flow through each branch road electric current and power, the voltage of each bus.
Step 72, according to calculation of tidal current, judges whether each branch road meets electric current nonoverload, whether each bus meets the not out-of-limit requirement of voltage.If do not met, excise the interior current load level of isolated island minimum, the load of meritorious demand maximum, re-starts trend by above-mentioned steps 71 and calculates, until isolated island meets institute's Prescribed Properties.
Step 8, distribution network system output isolated island scheme, obtains splitting point.Splitting point is exactly to carry out the switch that need to disconnect when isolated island is divided.
The present invention is based on harmony algorithm and realize the method for the thorough search of power distribution network isolated island partition problem solution space, can utilize individual local message and colony's global information to instruct algorithm further to search for, and its principle is simple, adjustable parameter is few, easily realizes.Harmony algorithm can obtain in the short period of time optimum isolated island splitting scheme in the time solving isolated island partition problem.
Followingly describe with simulative example:
Algorithm parameter HMS=5 is set, HMCR=0.9, PAR=0.1, Tmax=500000; In IEEE69 node distribution system, access DG as shown in Figure 2.At bus, 15,25,32,41,47,56,61,65 places add DG, are followed successively by DG1, DG2, DG3, DG4, DG5, DG6, DG7, DG8, and parameter is in table 1.Load in distribution system is classified by significance level: 1 stage load node is 6~9,13~15,20,22,33~35,40~42,48,52,54,55,58,62~64,68; 2 stage load nodes are 10~12,16~19,31,43~46,49,51,56,59; All the other are 3 stage loads.After external electrical network breaks down, disconnect and being connected of major network, divide the isolated island scheme that obtains with said method isolated island.
Table 1DG parameter
The foregoing is only preferred embodiment of the present invention, is only illustrative for invention, and nonrestrictive.Those skilled in the art is understood, and in the spirit and scope that limit, can carry out many changes to it in invention claim, amendment, and even equivalence, but all will fall within the scope of protection of the present invention.

Claims (10)

1. the power distribution network isolated island division methods based on harmony algorithm, is characterized in that, its process is:
Step 1, distribution network system arranges the initialization of harmony algorithm parameter;
Step 2, distribution network system is to each dimension of each individuality in power distribution network initialization of encoding, produce N initial individual;
Step 3, distribution network system calculates the target function value of each individuality in initial harmony data base HM;
Step 4, distribution network system produces new explanation according to HS algorithmic rule, and calculates the target function value of new explanation;
Step 5, distribution network system upgrades harmony data base HM; Whether the target function value that judges the new explanation that above-mentioned steps 4 calculates is greater than in described harmony data base HM the target function value of poor solution, if be greater than, replaces the poorest solution in described harmony data base HM with new explanation;
Step 6, distribution network system judges whether to meet the end condition of harmony algorithm: actual iterations is less than the iterations arranging in above-mentioned steps 1; If meet, carry out following step 7; If do not meet, return to step 4, and carry out successively;
Step 7, distribution network system is obtaining after current optimum solution, and the each isolated island in isolated island scheme is carried out to trend verification;
Step 8, distribution network system output isolated island scheme, obtains splitting point.
2. the power distribution network isolated island division methods based on harmony algorithm according to claim 1, is characterized in that, in above-mentioned steps 2, for each individuality, the cataloged procedure of its every one dimension is: for the node that is connected with DG, coding is initialized as 1; For the node that is connected with one-level or two stage loads, coding is initialized as 1; For the node that is connected with three stage loads, coding random initializtion is 0 or 1.
3. the power distribution network isolated island division methods based on harmony algorithm according to claim 1 and 2, it is characterized in that, in above-mentioned steps 3, first with breadth-first search, individuality is searched for, determine load and distributed power source DG that in splitting scheme, each isolated island comprises, judge whether each isolated island meets the requirement of constraint condition A and condition B, wherein, judge according to the following formula whether each isolated island meets constraint condition A and B
Wherein, Sump ifor the load restoration amount of single isolated island, m is the load number that this isolated island comprises, w jfor load weight, P jfor the active power of load; The basic ideas that above-mentioned constraint condition is judged as: in each isolated island scheme, may comprise n isolated island separated from one another, to single isolated island, judge whether this isolated island meets constraint condition A and B; If both meet simultaneously, calculate its load restoration amount, load restoration amount is the product of load active power and load level weights; If have at least one not meet, its load restoration amount set to 0.
4. the power distribution network isolated island division methods based on harmony algorithm according to claim 3, it is characterized in that, in above-mentioned steps 3, be calculated as follows the target function value of each individuality, the target function value that isolated island is divided is exactly the switch number that the load restoration amount sum of each isolated island in isolated island splitting scheme deducts action again
max ( &lambda; 1 &Sigma; i = 1 n Sump i - &lambda; 2 N breaker )
In formula, λ 1and λ 2for coefficient; Sump ifor the load restoration amount of single isolated island; N breakerthe number of switches that need to move for dividing isolated island.
5. the power distribution network isolated island division methods based on harmony algorithm according to claim 3, is characterized in that, in above-mentioned steps 3, constraint condition A is: active power balance in isolated island;
The active power that in isolated island, generator sends must be greater than the active power of isolated island internal loading demand: is total to as shown in formula as following,
&Sigma; G i &Element; V P G i > &Sigma; L j &Element; V P L j
In formula: for distributed power source generated energy; for the load in isolated island;
Above-mentioned constraint condition B is: the distributed power source DG that at least contains a firm power output in isolated island.
6. the power distribution network isolated island division methods based on harmony algorithm according to claim 1 and 2, is characterized in that, in above-mentioned steps 4, and the coding value of each in definite this solution of a new explanation generation need; Wherein, in new explanation, the coding of a certain position x specifically determines that mode is:
Step 41, if x corresponding node is distributed power source DG, x indirect assignment is 1;
Step 42, if x corresponding node is important load, x assignment is 1; Then produce a random number rand between [0,1] 2if, rand 2< PAR, x accepts disturbance and becomes 0; Otherwise, remaining 1 constant, PAR is parameter preset, it is worth between 0-1;
Step 43, if x corresponding node is general load, produces a random number rand between [0,1] 1if, rand 1< HMCR, the value of x equals the coding of the x position of any one individuality in harmony data base HM, and HMCR is parameter preset, and it is worth between 0-1; Otherwise x random coded is 0 or 1; , then produce a random number rand between [0,1] then 2if, rand 2< PAR, x accepts disturbance change coding, and from being 1 to become 0 or be that 0 to become 1, PAR be parameter preset originally originally, it is worth between 0-1; Otherwise the value of x remains unchanged.
7. the power distribution network isolated island division methods based on harmony algorithm according to claim 1, is characterized in that, in above-mentioned steps 7, the detailed process of each isolated island being carried out to trend verification is:
Step 71, selecting to have firm power in isolated island exports and has the distributed power source of certain frequency modulation and voltage modulation ability as the balance node of isolated island, other distributed power source and load, as PQ node, carry out trend calculating by Newton-Raphson method, and the trend that obtains whole isolated island distributes;
Step 72, according to calculation of tidal current, judges whether each branch road meets constraint condition C and D; If do not met, excise the interior current load level of isolated island minimum, the load of meritorious demand maximum, re-starts trend by above-mentioned steps 71 and calculates, until isolated island meets institute's Prescribed Properties.
8. the power distribution network isolated island division methods based on harmony algorithm according to claim 7, is characterized in that, in above-mentioned steps 7, carries out trend checking according to constraint condition C and D, and above-mentioned constraint condition C is: line current nonoverload; According to following formula,
I ij < I ij max
In formula: I ijfor flowing through the actual current of circuit; for flowing through the maximum allowed current of circuit.
9. the power distribution network isolated island division methods based on harmony algorithm according to claim 8, is characterized in that, above-mentioned constraint condition D is: busbar voltage is not out-of-limit; According to following formula,
U min<U i<U max
In formula: U ifor the virtual voltage of bus; U minand U maxfor the upper limit, the lower limit of busbar voltage permission.
10. the power distribution network isolated island division methods based on harmony algorithm according to claim 1 and 2, it is characterized in that, before this step 1 isolated island is divided, should determine the power distribution network region separating from large electrical network after fault isolation, and the scope of dividing using this as isolated island, power distribution network practical topology figure is reduced to a little and sideband has the non-directed graph of weights; Bus nodes is composed to the weights with the power of its connected load or distributed power source DG, then load and distributed power source DG node are deleted, power distribution network is further simplified to the node-branch road non-directed graph of node with weights.
CN201410224887.2A 2014-05-26 2014-05-26 A kind of power distribution network isolated island division methods based on harmony algorithm Active CN103971299B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410224887.2A CN103971299B (en) 2014-05-26 2014-05-26 A kind of power distribution network isolated island division methods based on harmony algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410224887.2A CN103971299B (en) 2014-05-26 2014-05-26 A kind of power distribution network isolated island division methods based on harmony algorithm

Publications (2)

Publication Number Publication Date
CN103971299A true CN103971299A (en) 2014-08-06
CN103971299B CN103971299B (en) 2016-09-28

Family

ID=51240752

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410224887.2A Active CN103971299B (en) 2014-05-26 2014-05-26 A kind of power distribution network isolated island division methods based on harmony algorithm

Country Status (1)

Country Link
CN (1) CN103971299B (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104485671A (en) * 2014-11-13 2015-04-01 国家电网公司 Voltage stability margin-based multi-microgrid system islanding method
CN104836333A (en) * 2015-05-13 2015-08-12 国家电网公司 Backup automatic switching apparatus monitoring method and apparatus
CN105119272A (en) * 2015-08-13 2015-12-02 浙江工业大学 Partition method of power distribution network island with electric vehicle battery swapping station based on hybrid algorithm
CN105372529A (en) * 2015-11-24 2016-03-02 国网甘肃省电力公司 Power grid islanding detection and anti-islanding protection method for large-scale wind power generation system
CN105406505A (en) * 2015-11-24 2016-03-16 国网甘肃省电力公司 Islanding state detection method for large-scale wind power
CN105486983A (en) * 2016-01-03 2016-04-13 国网江西省电力科学研究院 Fault-tolerance and distributed power supply contained power distribution network fault locating method
CN105978017A (en) * 2015-11-24 2016-09-28 国网甘肃省电力公司 Method for detection and protection of island with anti-theft and anti-tripping function aiming at large-scale wind power
CN107248754A (en) * 2017-07-03 2017-10-13 国家电网公司 A kind of addressing constant volume method of distributed power source
CN107565576A (en) * 2017-09-21 2018-01-09 国网福建省电力有限公司 A kind of active distribution network reactive Voltage Optimum method that more active management means are mutually coordinated
CN107909224A (en) * 2017-12-11 2018-04-13 国网山东省电力公司菏泽供电公司 The quick reconstruction method of power distribution network for removing infeasible solution based on harmony algorithm
CN110994612A (en) * 2019-12-25 2020-04-10 长沙理工大学 Power distribution network fault rapid recovery method based on network topology partition layering
CN111864791A (en) * 2019-12-31 2020-10-30 武汉大学 Decision table-based planned island online division method for power distribution network

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102868161A (en) * 2012-10-23 2013-01-09 四川大学 Optimization method of network variable structure with distributed type power supply distribution system
CN102983570A (en) * 2012-12-12 2013-03-20 天津市电力公司 Islanding method for power distribution network comprising distributed generation and considering interconnection switch
CN103001216A (en) * 2012-03-28 2013-03-27 河海大学 Method for quickly restoring power supply of power distribution network containing distributed power source

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103001216A (en) * 2012-03-28 2013-03-27 河海大学 Method for quickly restoring power supply of power distribution network containing distributed power source
CN102868161A (en) * 2012-10-23 2013-01-09 四川大学 Optimization method of network variable structure with distributed type power supply distribution system
CN102983570A (en) * 2012-12-12 2013-03-20 天津市电力公司 Islanding method for power distribution network comprising distributed generation and considering interconnection switch

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
M. MAHDAVI 等: "An improved harmony search algorithm for solving optimization problems", 《APPLIED MATHEMATICS AND COMPUTATION》 *
刘蓓 等: "一种配电网故障阶段式恢复方法", 《电力系统保护与控制》 *
李锰 等: "基于信息集成的配电网关键负荷保障", 《电工技术学报》 *
陈春 等: "基于基本环矩阵与改进和声搜索算法的配电网重构", 《电力系统自动化》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104485671B (en) * 2014-11-13 2017-01-11 国家电网公司 Voltage stability margin-based multi-microgrid system islanding method
CN104485671A (en) * 2014-11-13 2015-04-01 国家电网公司 Voltage stability margin-based multi-microgrid system islanding method
CN104836333A (en) * 2015-05-13 2015-08-12 国家电网公司 Backup automatic switching apparatus monitoring method and apparatus
CN104836333B (en) * 2015-05-13 2018-07-17 国家电网公司 Backup auto-activating device monitoring method and device
CN105119272A (en) * 2015-08-13 2015-12-02 浙江工业大学 Partition method of power distribution network island with electric vehicle battery swapping station based on hybrid algorithm
CN105978017B (en) * 2015-11-24 2018-10-26 国网甘肃省电力公司 There is the method for detecting and protecting isolated island of anti-jump function for large-scale wind power
CN105372529A (en) * 2015-11-24 2016-03-02 国网甘肃省电力公司 Power grid islanding detection and anti-islanding protection method for large-scale wind power generation system
CN105978017A (en) * 2015-11-24 2016-09-28 国网甘肃省电力公司 Method for detection and protection of island with anti-theft and anti-tripping function aiming at large-scale wind power
CN105372529B (en) * 2015-11-24 2019-03-19 国网甘肃省电力公司 For the power grid isolated island detection of large-scale wind power and anti-isolated island guard method
CN105406505A (en) * 2015-11-24 2016-03-16 国网甘肃省电力公司 Islanding state detection method for large-scale wind power
CN105486983A (en) * 2016-01-03 2016-04-13 国网江西省电力科学研究院 Fault-tolerance and distributed power supply contained power distribution network fault locating method
CN107248754A (en) * 2017-07-03 2017-10-13 国家电网公司 A kind of addressing constant volume method of distributed power source
CN107565576A (en) * 2017-09-21 2018-01-09 国网福建省电力有限公司 A kind of active distribution network reactive Voltage Optimum method that more active management means are mutually coordinated
CN107565576B (en) * 2017-09-21 2023-09-15 国网福建省电力有限公司 Reactive voltage optimization method for active power distribution network coordinated by multiple active management means
CN107909224A (en) * 2017-12-11 2018-04-13 国网山东省电力公司菏泽供电公司 The quick reconstruction method of power distribution network for removing infeasible solution based on harmony algorithm
CN110994612A (en) * 2019-12-25 2020-04-10 长沙理工大学 Power distribution network fault rapid recovery method based on network topology partition layering
CN110994612B (en) * 2019-12-25 2020-12-22 长沙理工大学 Power distribution network fault rapid recovery method based on network topology partition layering
CN111864791A (en) * 2019-12-31 2020-10-30 武汉大学 Decision table-based planned island online division method for power distribution network
CN111864791B (en) * 2019-12-31 2022-04-01 武汉大学 Decision table-based planned island online division method for power distribution network

Also Published As

Publication number Publication date
CN103971299B (en) 2016-09-28

Similar Documents

Publication Publication Date Title
CN103971299A (en) Power distribution network islanding method based on harmonic algorithm
Kumar et al. A review on overall control of DC microgrids
Muhtadi et al. Distributed energy resources based microgrid: Review of architecture, control, and reliability
Taher et al. Enhancing power grid synchronization and stability through time-delayed feedback control
Reddy et al. A review of Integration, Control, Communication and Metering (ICCM) of renewable energy based smart grid
Li et al. Controlled partitioning of a power network considering real and reactive power balance
Nguyen et al. Agent based restoration with distributed energy storage support in smart grids
Diaz et al. Fuzzy-logic-based gain-scheduling control for state-of-charge balance of distributed energy storage systems for DC microgrids
US20110106322A1 (en) Network connection manner of microgrid energy storage backup power source, and method for dispatching the same
Xiao et al. Online supervisory voltage control for grid interface of utility-level PV plants
CN106655253B (en) Single more microgrid regional dynamics division methods of three-phase
Wang et al. Unintentional islanding transition control strategy for three-/single-phase multimicrogrids based on artificial emotional reinforcement learning
Nemati et al. Optimization of microgrids short term operation based on an enhanced genetic algorithm
Sun et al. Implementation and CHIL testing of a microgrid control system
Shadabi et al. A decentralized non-linear dynamic droop control of a hybrid energy storage system bluefor primary frequency control in integrated AC-MTDC systems
Castro et al. Hierarchical optimisation strategy for energy scheduling and volt/var control in autonomous clusters of microgrids
Deng et al. Resilience-oriented dynamic distribution network with considering recovery ability of distributed resources
Yang et al. A novel partitioning method for the power grid restoration considering the support of multiple LCC-HVDC systems
Haiyun et al. A hierarchical control of microgrid based on droop controlled voltage source converter
CN108376997A (en) A kind of probabilistic active power distribution network isolated island division methods of consideration distributed generation resource
Peyghami-Akhuleh et al. Smart power management of DC microgrids in future milligrids
CN109412192A (en) From the back-to-back soft straight device operation method of energy storage multiterminal
Olvera et al. Using multi-terminal DC networks to improve the hosting capacity of distribution networks
Cai et al. Multistage bilevel planning model of energy storage system in urban power grid considering network reconfiguration
Saitoh et al. A proposal of emergency microgrid operation of distribution systems after large scale disasters

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant