CN112766532A - DG planning method based on improved mixed integer differential evolution algorithm - Google Patents
DG planning method based on improved mixed integer differential evolution algorithm Download PDFInfo
- Publication number
- CN112766532A CN112766532A CN202010999265.2A CN202010999265A CN112766532A CN 112766532 A CN112766532 A CN 112766532A CN 202010999265 A CN202010999265 A CN 202010999265A CN 112766532 A CN112766532 A CN 112766532A
- Authority
- CN
- China
- Prior art keywords
- power
- constraint
- distribution network
- individuals
- power distribution
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000035772 mutation Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003034 coal gas Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Primary Health Care (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (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 DG planning method based on an improved mixed integer differential evolution algorithm, which comprises the following steps: establishing a line loss minimum objective function of power distribution network reconstruction; establishing constraint conditions of power distribution network reconstruction, including power flow constraint, voltage constraint, DG power constraint and branch power constraint; and solving by adopting an improved mixed integer differential evolution algorithm to obtain the DG configuration of the power distribution network which meets the constraint condition and has the minimum line loss. The method and the system can reasonably plan the DG capacity accessed to the power distribution network, so that the power distribution network frame is more reasonably planned under the condition of accessing a large number of DGs.
Description
Technical Field
The invention relates to distributed energy (DG) capacity planning of an access distribution network, in particular to a DG planning method based on an improved mixed integer differential evolution algorithm.
Background
In recent years, with the gradual depletion of conventional energy sources, particularly petroleum, coal and natural gas, on the earth and the growing attention on the global problem of environmental protection and energy conservation, and the continuous development of power science and technology, the Distributed Generation (DG) technology mainly based on renewable energy sources such as wind power Generation and photovoltaic power Generation is becoming mature and becoming a hot point of research at home and abroad. However, the access of a large number of DG to the distribution network will change the structure and operation control manner of the conventional distribution network greatly. The research on reconstruction and fault recovery of the power distribution network containing multiple DGs is less, and especially, the technical method how to reasonably plan the island effect of the DGs is less in ensuring the power supply reliability of important loads and ensuring the economic operation, so that the reasonable planning of the DG capacity is very important and practical.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a DG planning method based on an improved mixed integer differential evolution algorithm, which can reasonably plan the DG capacity accessed to a power distribution network, so that the power distribution network frame planning is more reasonable under the condition of accessing a large amount of DGs.
The purpose of the invention is realized by the following technical scheme.
The DG planning method based on the improved mixed integer differential evolution algorithm comprises the following processes:
the method comprises the following steps: establishing a line loss minimum objective function of power distribution network reconstruction;
step two: establishing constraint conditions of power distribution network reconstruction, including power flow constraint, voltage constraint, DG power constraint and branch power constraint;
step three: and solving by adopting an improved mixed integer differential evolution algorithm to obtain the DG configuration of the power distribution network which meets the constraint condition and has the minimum line loss.
In the first step, the line loss minimum objective function is:
in the formula: r isijIs the resistance of branch i-j; pij、QijRespectively the active power and the reactive power flowing through each time period at the tail ends of the branches i-j; vijIs the node voltage at the end of branch i-j; n is the number of nodes in the power distribution network; k is a radical ofijThe state variable of the switch for branch i-j isDiscrete amount, 0 for open, 1 for closed; min F is the minimum function of line loss.
And step two, calculating the constraint condition of the power distribution network reconstruction according to the following formula:
and (3) power flow constraint:
in the formula: pi、QiRespectively input active power and reactive power of a node i; pDGi、QDGiRespectively injecting active power and reactive power into the node i for the DG; pDi、QDiRespectively the active power and the reactive power of the load at the node i; u shapei、UjVoltages of nodes i and j, respectively; gij、BijAnd deltaijRespectively the conductance, susceptance and phase angle difference of the branches i-j; n is the number of nodes in the power distribution network;
voltage constraint:
in the formula: n is the number of nodes in the power distribution network; u shapei,minAnd Ui,maxRespectively representing the lower limit and the upper limit of the allowable voltage of the node i; i isijIs the current of branch i-j; i isij,maxThe upper limit of the current for branch i-j;
DG power constraint:
in the formula: pDGi、QDGiRespectively injecting active power and reactive power into the node i for the DG; pDG,max、PDG,minThe maximum value and the minimum value of the active power of the DG are respectively; qDG,max、QDG,minMaximum and minimum reactive power of DG respectively;
branch power constraint:
The solving process by adopting the improved mixed integer differential evolution algorithm in the third step is as follows:
step 1: initializing parameter population size NP (total number of individuals), wherein each individual is a D-dimensional vector, a scaling factor F and a cross factor CR;
step 2: randomly generating an initial population P in a decision space;
step 3: calculating the fitness of each individual, namely an objective function value;
step 6: crossing individuals in the population;
step 7: and (3) selecting the individuals in the population: when new vector individualsIs more than the target vector individualWhen better, reserveTo the next generation population, otherwise target vector individualsStill remain in the population, once again as the next generation parent vector;
step 8: and if the maximum iteration number is reached or the error requirement is met, obtaining the optimal DG configuration result, and if not, returning to Step 3.
Compared with the prior art, the technical scheme of the invention has the following beneficial effects:
the method for planning the DGs based on the improved mixed integer differential evolution algorithm promotes the economic operation of the power distribution network under the condition that a large amount of DGs are accessed into the power distribution network, and ensures that the planning configuration of the DGs is more reasonable and more adaptive.
Drawings
Fig. 1 is a topological structure diagram of a 10kV power distribution network.
Detailed Description
For a better understanding of the present invention, reference is made to the following description taken in conjunction with the accompanying drawings and examples.
The DG planning method based on the improved mixed integer differential evolution algorithm comprises the following processes:
the method comprises the following steps: establishing a line loss minimum objective function of power distribution network reconstruction, as follows:
in the formula: r isijIs the resistance of branch i-j; pij、QijRespectively the active power and the reactive power flowing through each time period at the tail ends of the branches i-j; vijIs the node voltage at the end of branch i-j; n is the number of nodes in the power distribution network; k is a radical ofijThe state variable of the switch for branch i-j isDiscrete amount, 0 for open, 1 for closed; minF is a minimum function of line loss.
Step two: and establishing constraint conditions of power distribution network reconstruction, including power flow constraint, voltage constraint, DG power constraint and branch power constraint.
The power flow constraint is calculated according to the following formula:
in the formula: pi、QiRespectively input active power and reactive power of a node i; pDGi、QDGiRespectively injecting active power and reactive power into the node i for the DG; pDi、QDiRespectively the active power and the reactive power of the load at the node i; u shapei、UjVoltages of nodes i and j, respectively; gij、BijAnd deltaijRespectively the conductance, susceptance and phase angle difference of the branches i-j; and n is the number of nodes in the power distribution network.
The voltage constraint is calculated as follows:
in the formula: n is the number of nodes in the power distribution network; u shapei,minAnd Ui,maxRespectively representing the lower limit and the upper limit of the allowable voltage of the node i; i isijIs the current of branch i-j; i isij,maxThe upper limit for the current in branch i-j.
The DG power constraint is calculated as follows:
in the formula: pDGi、QDGiRespectively injecting active power and reactive power into the node i for the DG; pDG,max、PDG,minThe maximum value and the minimum value of the active power of the DG are respectively; qDG,max、QDG,minThe maximum and minimum values of the reactive power of DG, respectively.
The branch power constraint is calculated according to the following formula:
Step three: and solving by adopting an improved mixed integer differential evolution algorithm to obtain the DG configuration of the power distribution network which meets the constraint condition and has the minimum line loss.
Solving using the modified mixed integer differential evolution algorithm generally follows the following steps:
step 1: initializing parameter population size NP (total number of individuals), wherein each individual is a D-dimensional vector, a scaling factor F and a cross factor CR;
step 2: randomly generating an initial population P in a decision space;
step 3: calculating the fitness of each individual, namely an objective function value;
step 6: crossing individuals in the population;
step 7: and (3) selecting the individuals in the population: when new vector individualsIs more than the target vector individualWhen better, reserveTo the next generation population, otherwise target vector individualsStill remain in the population, once again as the next generation parent vector;
step 8: and if the maximum iteration number is reached or the error requirement is met, obtaining the optimal DG configuration result, and if not, returning to Step 3.
Example (b):
the calculation example takes an IEEE 33 node system as an example, network reconstruction calculation is carried out by using the DG planning method based on the improved mixed integer differential evolution algorithm and taking the minimum line loss as an objective function, and meanwhile, the lowest node voltage when the line loss is minimum is output. The method comprises the following steps:
the method comprises the following steps: and establishing a reconstructed line loss minimum objective function of the power distribution network, as shown in formula (1).
Step two: and establishing constraint conditions for power distribution network reconstruction, as shown in the formulas (2) to (6).
Step three: the specific solving process of the improved differential evolution algorithm comprises the following steps:
step 1: parameters population size NP (total number of individuals) are initialized, each individual being a D-dimensional vector, a scaling factor F, a crossover factor CR. (ii) a
Step 2: randomly generating an initial population P in a decision space;
step 3: calculating the fitness of each individual, namely an objective function value;
step 6: crossing individuals in the population;
step 7: and (3) selecting the individuals in the population: when new vector individualsIs more than the target vector individualWhen better, reserveTo the next generation population, otherwise target vector individualsStill remain in the population, once again as the next generation parent vector;
step 8: if the maximum iteration number is reached or the error requirement is met, obtaining the optimal DG configuration result, otherwise, returning to Step 3: .
The power distribution network is provided with 33 nodes and 37 branches, wherein 5 of the nodes are interconnection switches, the rated voltage is 12.66kV, the total load is 5048.26kVA + j2547.32kVar, and the topological structure of the power distribution network is shown in figure 1.The DG installation location and capacity are shown in table 1. Initial grid loss without distributed generation is 0.02023 p.u., lowest node voltage 0.9132p.u. parameter set: population size NP is 5, maximum number of generations of variation Gmax60, the scaling factor F is 0.5 and the crossover factor CR is 0.5.
TABLE 1 installation location, type and Capacity
While the present invention has been described in terms of its functions and operations with reference to the accompanying drawings, it is to be understood that the invention is not limited to the precise functions and operations described above, and that the above-described embodiments are illustrative rather than restrictive, and that various changes and modifications may be effected therein by one skilled in the art without departing from the scope or spirit of the invention as defined by the appended claims.
Claims (4)
1. A DG planning method based on an improved mixed integer differential evolution algorithm is characterized by comprising the following processes:
the method comprises the following steps: establishing a line loss minimum objective function of power distribution network reconstruction;
step two: establishing constraint conditions of power distribution network reconstruction, including power flow constraint, voltage constraint, DG power constraint and branch power constraint;
step three: and solving by adopting an improved mixed integer differential evolution algorithm to obtain the DG configuration of the power distribution network which meets the constraint condition and has the minimum line loss.
2. The DG planning method based on the improved mixed integer differential evolution algorithm of claim 1, wherein the line loss minimum objective function in the first step is:
in the formula: r isijIs the resistance of branch i-j; pij、QijRespectively the active power and the reactive power flowing through each time period at the tail ends of the branches i-j; vijIs the node voltage at the end of branch i-j; n is the number of nodes in the power distribution network; k is a radical ofijThe state variable of the branch i-j switch is 0-1 discrete quantity, wherein 0 represents opening, and 1 represents closing; minF is a minimum function of line loss.
3. The DG planning method based on the improved mixed integer differential evolution algorithm of claim 1, wherein the constraint conditions for the reconstruction of the distribution network in the second step are calculated according to the following formula:
and (3) power flow constraint:
in the formula: pi、QiRespectively input active power and reactive power of a node i; pDGi、QDGiRespectively injecting active power and reactive power into the node i for the DG; pDi、QDiRespectively the active power and the reactive power of the load at the node i; u shapei、UjVoltages of nodes i and j, respectively; gij、BijAnd deltaijRespectively the conductance, susceptance and phase angle difference of the branches i-j; n is the number of nodes in the power distribution network;
voltage constraint:
in the formula: n is the number of nodes in the power distribution network; u shapei,minAnd Ui,maxRespectively representing the lower limit and the upper limit of the allowable voltage of the node i; i isijIs the current of branch i-j; i isij,maxThe upper limit of the current for branch i-j;
DG power constraint:
in the formula: pDGi、QDGiRespectively injecting active power and reactive power into the node i for the DG; pDG,max、PDG,minThe maximum value and the minimum value of the active power of the DG are respectively; qDG,max、QDG,minMaximum and minimum reactive power of DG respectively;
branch power constraint:
4. The DG planning method based on the improved mixed integer differential evolution algorithm of claim 1, wherein the solving process adopting the improved mixed integer differential evolution algorithm in the third step is as follows:
step 1: initializing parameter population size NP (total number of individuals), wherein each individual is a D-dimensional vector, a scaling factor F and a cross factor CR;
step 2: randomly generating an initial population P in a decision space;
step 3: calculating the fitness of each individual, namely an objective function value;
step 6: crossing individuals in the population;
step 7: and (3) selecting the individuals in the population: when new vector individualsIs more than the target vector individualWhen better, reserveTo the next generation population, otherwise target vector individualsStill remain in the population, once again as the next generation parent vector;
step 8: and if the maximum iteration number is reached or the error requirement is met, obtaining the optimal DG configuration result, and if not, returning to Step 3.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010999265.2A CN112766532A (en) | 2020-09-22 | 2020-09-22 | DG planning method based on improved mixed integer differential evolution algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010999265.2A CN112766532A (en) | 2020-09-22 | 2020-09-22 | DG planning method based on improved mixed integer differential evolution algorithm |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112766532A true CN112766532A (en) | 2021-05-07 |
Family
ID=75693084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010999265.2A Pending CN112766532A (en) | 2020-09-22 | 2020-09-22 | DG planning method based on improved mixed integer differential evolution algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112766532A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116187723A (en) * | 2023-04-26 | 2023-05-30 | 佰聆数据股份有限公司 | Resource scheduling method and device applied to distribution line loss reduction scene |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103972927A (en) * | 2014-05-26 | 2014-08-06 | 武汉大学 | Integrated control method for transforming microgrid containing photovoltaic/stored energy generating system from connected grid to isolated grid |
CN105719196A (en) * | 2016-01-18 | 2016-06-29 | 天津大学 | Active power distribution network pressure reactive power control method based on intelligent soft normally open point |
-
2020
- 2020-09-22 CN CN202010999265.2A patent/CN112766532A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103972927A (en) * | 2014-05-26 | 2014-08-06 | 武汉大学 | Integrated control method for transforming microgrid containing photovoltaic/stored energy generating system from connected grid to isolated grid |
CN105719196A (en) * | 2016-01-18 | 2016-06-29 | 天津大学 | Active power distribution network pressure reactive power control method based on intelligent soft normally open point |
Non-Patent Citations (2)
Title |
---|
肖镇铭;: "基于粒子群算法的含分布式发电配电网规划", 云南电力技术 * |
邹必昌;孙洪斌;龚庆武;刘栋;陈道君;: "含分布式发电的改进混合整数差分算法的配电网重构", 电力系统保护与控制 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116187723A (en) * | 2023-04-26 | 2023-05-30 | 佰聆数据股份有限公司 | Resource scheduling method and device applied to distribution line loss reduction scene |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111342461B (en) | Power distribution network optimal scheduling method and system considering dynamic reconfiguration of network frame | |
CN110348048B (en) | Power distribution network optimization reconstruction method based on consideration of heat island effect load prediction | |
CN103310065B (en) | Meter and the intelligent network distribution reconstructing method of distributed power generation and energy-storage units | |
CN109599894B (en) | DG grid-connected optimization configuration method based on improved genetic algorithm | |
Tang et al. | Study on day-ahead optimal economic operation of active distribution networks based on Kriging model assisted particle swarm optimization with constraint handling techniques | |
CN109038545B (en) | Power distribution network reconstruction method based on differential evolution invasive weed algorithm | |
Azimi et al. | Multiobjective daily Volt/VAr control in distribution systems with distributed generation using binary ant colony optimization | |
CN108631278B (en) | The Optimal Configuration Method of breaker and fault current limiter in a kind of looped network formula direct-current micro-grid | |
Li et al. | Active distribution network active and reactive power coordinated dispatching method based on discrete monkey algorithm | |
CN112766532A (en) | DG planning method based on improved mixed integer differential evolution algorithm | |
You et al. | Optimal topology of LVAC in a rural village using Water Cycle Algorithm | |
Diaaeldin et al. | A novel reconfiguration methodology of radial distribution systems for power loss minimization using expanded invasive weed optimization | |
CN117335384A (en) | Power distribution network fault recovery reconstruction method based on ant colony algorithm | |
Cho et al. | A selecting method of optimal load on time varying distribution system for network reconfiguration considering DG | |
CN114844051A (en) | Reactive power supply optimal configuration method and terminal for active power distribution network | |
Siqing et al. | Distribution network fault restoration based on improved adaptive genetic algorithm | |
Wang et al. | Reconfiguration of distribution network with dispersed generators based on improved forward-backward sweep method | |
CN111740419A (en) | Active power distribution network fault recovery method based on differential evolution algorithm | |
Dong et al. | Distribution Network Reconfiguration Method with Distributed Generators Based on an Improved Shuffled Frog Leaping Algorithm | |
Doagou-Mojarrad et al. | Probabilistic interactive fuzzy satisfying generation and transmission expansion planning using fuzzy adaptive chaotic binary PSO algorithm | |
CN110994665A (en) | Distributed photovoltaic multi-point access low-voltage distribution network site selection method | |
Li et al. | Combining Differential Evolution Algorithm with biogeography-based optimization algorithm for reconfiguration of distribution network | |
Liu et al. | Multi-objective flexible planning of transmission network considering generation and load uncertainties | |
Bai | Reconfiguration performance of the urban power distribution system based on the genetic-ant colony fusion algorithm | |
Zhengji et al. | Research on self-healing function of smart distribution network based on adaptive genetic algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210507 |