CN108964124B - Distributed power supply grid-connected optimization configuration method for counting and electricity price response in incremental power distribution network - Google Patents

Distributed power supply grid-connected optimization configuration method for counting and electricity price response in incremental power distribution network Download PDF

Info

Publication number
CN108964124B
CN108964124B CN201810793769.1A CN201810793769A CN108964124B CN 108964124 B CN108964124 B CN 108964124B CN 201810793769 A CN201810793769 A CN 201810793769A CN 108964124 B CN108964124 B CN 108964124B
Authority
CN
China
Prior art keywords
distribution network
power
power distribution
electricity
incremental
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.)
Active
Application number
CN201810793769.1A
Other languages
Chinese (zh)
Other versions
CN108964124A (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.)
State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch
State Grid Corp of China SGCC
Southeast University
State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch
State Grid Corp of China SGCC
Southeast University
State Grid Jiangsu Electric Power 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 State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch, State Grid Corp of China SGCC, Southeast University, State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Jiangsu Electric Power Co ltd Suqian Power Supply Branch
Priority to CN201810793769.1A priority Critical patent/CN108964124B/en
Publication of CN108964124A publication Critical patent/CN108964124A/en
Application granted granted Critical
Publication of CN108964124B publication Critical patent/CN108964124B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H02J3/383
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/14Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2310/00The network for supplying or distributing electric power characterised by its spatial reach or by the load
    • H02J2310/50The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads
    • H02J2310/56The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads characterised by the condition upon which the selective controlling is based
    • H02J2310/62The condition being non-electrical, e.g. temperature
    • H02J2310/64The condition being economic, e.g. tariff based load management
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
    • Y02B70/3225Demand response systems, e.g. load shedding, peak shaving
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S20/00Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20End-user application control systems
    • Y04S20/222Demand response systems, e.g. load shedding, peak shaving

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a distributed power supply grid-connected optimization configuration method for counting and responding to electricity prices in an incremental power distribution network, relates to a distributed power supply grid-connected optimization configuration method for counting and responding to electricity prices, is a multi-objective nonlinear programming method, can be applied to the location determination and the peak-valley electricity price determination of a distributed power supply of the power distribution network, and belongs to the field of power system planning. The method comprises the steps of obtaining the wind speed and the sunlight intensity of the incremental distribution network area
Figure 100004_DEST_PATH_IMAGE001
The method comprises the steps of (1) data distribution on a full time sequence of each time period, and constructing a wind power and photovoltaic output model; building each load node in incremental distribution network
Figure 509671DEST_PATH_IMAGE001
Distributing in full time sequence of each time interval, and constructing a load demand response model of the node; comprehensively considering the benefits of main parties of the incremental power distribution network market, and constructing objective functions of the main parties; the genetic algorithm is improved, so that the algorithm can solve the optimization problem related to the invention, the position and capacity distribution of the distributed power supply on each node of the power distribution network is obtained according to the result obtained by the genetic algorithm, and the value of the peak-valley time-of-use electricity price is obtained.

Description

Distributed power supply grid-connected optimization configuration method for counting and electricity price response in incremental power distribution network
Technical Field
The invention discloses a distributed power supply grid-connected optimization configuration method for counting and responding to electricity prices in an incremental power distribution network, relates to a distributed power supply grid-connected optimization configuration method for counting the electricity price responses, is a multi-objective nonlinear programming method, can be applied to site selection and volume determination of a distributed power supply of the power distribution network and establishment of peak-valley electricity prices, and belongs to the field of power system programming.
Background
As society develops, the demand for energy is increasing. In recent years, with the rapid development of scientific technology, various new energy sources are effectively developed and utilized, and the distributed power generation technology is gradually applied to some power distribution networks. The distributed power supply is used as a new energy source mainly developed by renewable energy sources, has the characteristics of randomness and uncontrollable performance, and can be scientifically connected to the distributed power supply, so that the situation that the current energy sources are increasingly tense can be relieved, and the electric energy quality of a power distribution network area can be effectively improved. Therefore, the first problem of accessing the distributed power supply to the power distribution network is where and when to access a large-capacity power supply, and the problems to be considered include the economic efficiency, reliability and the like of the power distribution network, and the problems should be reasonably considered and planned before the distributed power supply is connected to the power distribution network.
The electricity price response is a demand response behavior which mainly aims at stimulating users to participate in power grid peak shaving and guides the users to scientifically and reasonably use electricity. The electricity price response can realize that the power consumption is optimized, reduce the peak-valley difference, increase new forms of energy and consume, this not only needs the power company to formulate reasonable electricity price to guide the user to shift peak load, considers that distributed power source itself just has the chronology in the distribution network, the event still should combine reasonable planning overall arrangement of distributed power generation and configuration optimization to reduce the great electricity sales income of power supplier loss that the user that the demand elasticity is big loses the load and cause, just so can realize distributed power source and utilize and accomplish the overall planning of multiple target and compromise most rationally. For the problem of configuration optimization of a distributed power supply in a power distribution network, the following defects or shortcomings mainly exist in previous related researches:
1. distribution characteristics of the distributed power supply on a full time sequence are not considered, so that a certain error is generated in an optimization result of DG in power distribution network site selection and constant volume.
2. Some researches only consider the capacity configuration under the fixed-point condition or only consider the address selection problem with fixed capacity, and although the difficulty of the problem and the complexity of the algorithm are reduced to a certain extent, the application scene is greatly limited, and the obtained solution is not accurate enough.
3. Aiming at the problem of site selection optimization of DGs in the power distribution network, the selected intelligent optimization algorithm needs to be properly improved so as to accelerate convergence and not lose feasibility and reliability of solution.
4. Most studies only consider the analysis and research of a single optimization objective, and in consideration of actual needs, a plurality of optimization objectives are usually processed for decision-making and application.
5. Partial research carries out double-layer optimization on demand response and optimal configuration or scheduling of DGs, and the global optimal possibility of the solution in the complete set is reduced to a certain extent.
The operation mode of the power grid company in the incremental power distribution network is changed from 'purchasing electricity price difference' to 'admitting income' (general purchasing and general selling → logistics), the charging mode is changed from the state network company for replacing the power plant and charging the electricity selling company for uniformly charging the electricity to the electricity selling company by the user, and the electricity selling company respectively pays the electricity distribution fee to the power grid enterprise and pays the electricity to the electricity generating enterprise. Therefore, in the current distributed power optimization problem of the incremental power distribution network, a three-way game of a power grid enterprise, a user and a power generation enterprise needs to be considered, and multi-party target optimization needs to be constructed so as to provide multi-party decision for planning or prepare for deep multi-party cooperative planning.
Disclosure of Invention
In order to overcome various problems and difficulties caused by the fact that the distributed power supply is connected into a power distribution network to fix the volume and select the site, the invention provides a distributed power supply grid-connected optimization configuration method considering electricity price response in an incremental power distribution network.
The invention is realized by adopting the following technical scheme:
the invention discloses a distributed power supply grid-connected optimization configuration method for counting and responding to electricity prices in an incremental power distribution network, which comprises the following steps of:
step S1, acquiring the wind speed and the sunlight intensity of the incremental distribution network area
Figure DEST_PATH_IMAGE001
The method comprises the steps of (1) data distribution on a full time sequence of each time period, and constructing a wind power and photovoltaic output model;
step S2, constructing each load node in the incremental distribution network
Figure 706036DEST_PATH_IMAGE001
Distribution in each time period and full time sequence is realized, differences of different types of industries in response to power price change are considered, and acquisition nodes are constructed according to similar historical data of demand response of other power distribution networks
Figure 623176DEST_PATH_IMAGE001
The method comprises the steps of (1) establishing a load demand response model by using a time-interval full-time-sequence crossed elastic matrix;
step S3, the power supply main body aims at minimizing the load peak-valley difference of the power distribution network, minimizing the total investment cost of the power distribution network and minimizing the network loss of the power distribution network, and the power selling main body aims at maximizing the power selling income and constructs a multi-objective optimization function;
step S4, improving the genetic algorithm to solve, according to a plurality of objective functions in the power supply main body, establishing the weight based on the variance, calculating the comprehensive fitness of the individual, synthesizing the external objective function, namely synthesizing the single objective function of the power selling main body, obtaining the final total fitness, and adopting the replication operator, the self-adaptive crossover operator and the mutation operator of the elite strategy to accelerate the convergence;
s5, decoding according to the result of the genetic algorithm to obtain the distributed power distribution point position and the configuration capacity;
the above-mentioned
Figure 402913DEST_PATH_IMAGE001
Typically 24 or 96.
Compared with other research methods of the same type, the method of the invention has the advantages that:
1. considering full time sequence distribution of distributed power supply output;
2. solving a combined planning problem of site selection, volume fixing and electricity price response of distributed power generation by adopting an improved genetic algorithm;
3. multi-target optimization is carried out by considering multi-party benefits of the incremental power distribution network;
the DG optimal configuration method considering the electricity price response, which is constructed by the invention, can provide better reference for the plan and decision of incremental power distribution service development.
Drawings
Fig. 1 is a flowchart of a distributed power grid-connected optimization configuration method for counting and responding to electricity prices in an incremental power distribution network.
FIG. 2 is a flow chart of an improved genetic algorithm in accordance with the present invention.
Detailed Description
The invention will be further described with reference to fig. 1 and 2.
The distributed power supply grid-connected optimization configuration method for counting and responding to electricity prices in the incremental power distribution network comprises the following steps of:
in the following examples, reference will be made to the description above
Figure 860440DEST_PATH_IMAGE001
The time periods are set to 24 hours;
And S1, acquiring data distribution of wind speed and sunlight intensity of the incremental distribution network area in a full 24-hour time sequence, and constructing a wind power and photovoltaic output model.
S1.1, acquiring wind speed data of a 24-hour full time sequence in a typical day in the incremental power distribution network from a meteorological office, and establishing the following wind power generation output model for a single distributed wind power generator:
Figure 100002_DEST_PATH_IMAGE002
(9)
in the formula (9), the reaction mixture is,
Figure DEST_PATH_IMAGE003
is the wind speed of the distribution network area for a certain period of time,
Figure 100002_DEST_PATH_IMAGE004
cutting into the wind speed for the wind power generator,
Figure DEST_PATH_IMAGE005
is the rated wind speed, and is,
Figure 100002_DEST_PATH_IMAGE006
in order to cut out the wind speed,
Figure DEST_PATH_IMAGE007
is the rated power of a single distributed wind power generator.
S1.2, acquiring the 24-hour full-time illumination intensity distribution of the incremental distribution network in a typical day from a meteorological office, and constructing the following photovoltaic output model for a single distributed photovoltaic cell:
Figure 100002_DEST_PATH_IMAGE008
(10)
in the formula (10), the compound represented by the formula (10),
Figure DEST_PATH_IMAGE009
the active power output of a single photovoltaic cell is provided,
Figure 100002_DEST_PATH_IMAGE010
the size of the area of the battery is,
Figure DEST_PATH_IMAGE011
it is the intensity of the light that is irradiated,
Figure 100002_DEST_PATH_IMAGE012
the photoelectric conversion efficiency.
And S2, constructing the distribution of all load nodes in the incremental power distribution network in a full time sequence of 24 hours, and constructing a load demand response model of the nodes. In practice, the subscriber unit power factor is high and the power in calculating the bus bar nodes only takes into account the demand response of the active power. Acquiring the average electricity price of the power distribution network area
Figure DEST_PATH_IMAGE013
The following requirement elastic model is established according to the distribution data of each node in the full time sequence of 24 hours:
Figure 100002_DEST_PATH_IMAGE014
(11)
in the formula (11), the reaction mixture is,
Figure DEST_PATH_IMAGE015
the elastic coefficient of the electricity quantity in the m time period to the electricity price change in the n time period,
Figure 100002_DEST_PATH_IMAGE016
for a certain historical average electricity price,
Figure DEST_PATH_IMAGE017
in order to correspond to the electric quantity demand of the m-th time period under the electricity price,
Figure 100002_DEST_PATH_IMAGE018
Figure DEST_PATH_IMAGE019
electricity consumption in m periods respectivelyThe amount variation and the electricity price in the n period are compared with
Figure 111030DEST_PATH_IMAGE016
The change that occurred.
When analyzing the influence factors of different user power demands, the invention mainly considers the following two aspects:
1. analyzing the bearing capacity of the user to the electricity price change, namely the magnitude of the self-elastic coefficient;
2. according to different industrial characteristics and electricity utilization habits, the capacity of transferring the electricity consumption in the peak time period to the electricity consumption in the valley time period, namely the size of the cross elasticity coefficient, of a user is analyzed when the peak-valley electricity price is analyzed.
The demand elasticity of different types of industries for 24 hours can be approximately calculated according to historical data of other power systems or other similar power distribution areas, so the demand elasticity coefficient of the current power distribution area can be calculated by the following formula:
Figure 100002_DEST_PATH_IMAGE020
(12)
in the formula (12), the reaction mixture is,
Figure DEST_PATH_IMAGE021
at the basic electricity price for the k-th industry
Figure 995809DEST_PATH_IMAGE016
The weight of the load occupied by the lower period,
Figure 100002_DEST_PATH_IMAGE022
the corresponding spring constant.
In summary, the electric quantity variation of the full time sequence of the user at a certain node of the incremental distribution network can be obtained by the following formula:
Figure DEST_PATH_IMAGE023
(13)
the electricity consumption of the user t period after the electricity price response is carried out is as follows:
Figure 100002_DEST_PATH_IMAGE024
(14)
and step S3, constructing a multi-party objective function in consideration of different interest points of power generation and power supply parties.
From the grid company's perspective, assuming that the distributed power operations are all invested in full by the power supply company, the objective function may be considered:
Figure DEST_PATH_IMAGE025
(15)
in the formula (15), the first objective function is that the total loss of the network is minimum, the second objective function is that the peak-valley difference of the access public network node is minimum, the third objective function is that the total investment of the power grid company is minimum, wherein the first two items are the total investment of wind power and photovoltaic,
Figure 100002_DEST_PATH_IMAGE026
for the cost of the generation of the unit active power of the wind power,
Figure DEST_PATH_IMAGE027
in order to reduce the power generation cost per unit area of the photovoltaic cell,
Figure 100002_DEST_PATH_IMAGE028
and
Figure DEST_PATH_IMAGE029
the output of a single fan and the output of a photovoltaic cell are respectively, the third term is the distributed generation benefit,
Figure 100002_DEST_PATH_IMAGE030
in order to increase the price of the power transmission and distribution of the power distribution network,
Figure DEST_PATH_IMAGE031
and (4) distributing the electric quantity participating in market trading for the t period. From a power plant perspective, an objective function of the power selling income maximization of the whole network can be considered:
Figure 100002_DEST_PATH_IMAGE032
(16)
Figure DEST_PATH_IMAGE033
for the time period t, the price of selling electricity,
Figure 100002_DEST_PATH_IMAGE034
the electric quantity in the first time period under a certain node,
Figure DEST_PATH_IMAGE035
the total number of nodes in the incremental power distribution network.
The following inequality constraint constraints are considered:
Figure 100002_DEST_PATH_IMAGE036
(17)
in the formula (17), the first term is constrained to be that the voltage is not out of limit, the second term is constrained to be that the distributed generation permeability is limited, the ratio of the total capacity of the accessed DGs to the total capacity of the power distribution network does not exceed the permeability, and the third term is limited to be that the electricity price does not exceed the maximum electricity price
Figure DEST_PATH_IMAGE037
And not lower than the lowest electricity price
Figure 100002_DEST_PATH_IMAGE038
And step S4, improving the genetic algorithm to solve, and acquiring the distributed power distribution point position and the configuration capacity of the distributed power. Assuming that the rated power and the battery area of a single distributed wind power and photovoltaic are fixed, so that the number of the wind power and photovoltaic configurations on each node of the distribution network area is obtained, the following method can be adopted for coding:
Figure DEST_PATH_IMAGE039
(18)
in the formula (18), the reaction mixture,
Figure 100002_DEST_PATH_IMAGE040
Figure DEST_PATH_IMAGE041
are respectively the first
Figure 100002_DEST_PATH_IMAGE042
The number of fans and photovoltaic devices installed on each node,
Figure DEST_PATH_IMAGE043
Figure 100002_DEST_PATH_IMAGE044
Figure DEST_PATH_IMAGE045
for the peak period, the usual period, and the valley period of the region, the electricity prices can be generally distributed as follows with respect to the period division: peaks (8: 00-11: 00, 16: 00-21: 00), plateaus (7: 00-8: 00, 11: 00-16: 00, 21: 00-23: 00), valleys (23: 00-7: 00).
The genetic algorithm solution steps are as follows:
y1, setting the population size N and setting the maximum variation rate
Figure 100002_DEST_PATH_IMAGE046
And minimum rate of variation
Figure DEST_PATH_IMAGE047
Setting the maximum crossing rate
Figure 100002_DEST_PATH_IMAGE048
And minimum crossing rate
Figure DEST_PATH_IMAGE049
Maximum number of iterations;
y2, forming an initial population, wherein the gene X adopts integer coding to represent the number of distributed power supplies, and the gene P adopts real number coding;
y3, calculating the value of each objective function. Firstly, considering an objective function of a power grid company level, because the objective functions are a plurality of objective functions in the company and no interest dispute exists among the objective functions, under the condition of seeking overall optimum, the following method is adopted to calculate the weight and the comprehensive fitness of each objective function:
and S4.1, carrying out load flow on the power distribution network represented by each individual, and calculating the value of each objective function of all the individuals. And (3) adopting a normalization method to scale each objective function value to an interval as the fitness of a single objective function, and assuming that
Figure 100002_DEST_PATH_IMAGE050
And
Figure DEST_PATH_IMAGE051
for the fitness of the maximum and minimum population of the current algebra, the normalization formula is:
Figure 100002_DEST_PATH_IMAGE052
(19)
s4.2, calculating the variance of the fitness values of all individuals on each objective function, wherein the larger the weight of the objective function with the large variance is, the larger the search area of the Pareto optimal solution is, so that the possibility that the subsequent solution is too concentrated, the larger the loss of the optimal solution is avoided:
Figure DEST_PATH_IMAGE053
(20)
in the formula (20), the reaction mixture is,
Figure 100002_DEST_PATH_IMAGE054
is as follows
Figure DEST_PATH_IMAGE055
The weight of the respective objective function is,
Figure 100002_DEST_PATH_IMAGE056
is as follows
Figure 276225DEST_PATH_IMAGE055
The fitness of each objective function is calculated,
Figure DEST_PATH_IMAGE057
the number of objective functions is considered for the grid company. Therefore, the total fitness of the genetic algorithm individual multi-objective function for the power grid company is as follows:
Figure 100002_DEST_PATH_IMAGE058
(21)
in the formula (21), the compound represented by the formula,
Figure DEST_PATH_IMAGE059
is the fitness of the kth objective function. Similarly, the generating company fitness function is calculated in the same way.
And S4.3, considering the cooperation level of the comprehensive power generation company and the power supply company, and having agreed weight (for example, 1: 1) on each target, wherein the final fitness of each individual is as follows:
Figure 100002_DEST_PATH_IMAGE060
(22)
s4.4, copying operation. Adopting individual elite strategy to select the one with the maximum fitness
Figure DEST_PATH_IMAGE061
Directly enters the next generation operation as an individual elite. The rest of the
Figure 100002_DEST_PATH_IMAGE062
And selecting the wheel disc probability according to the fitness of the individual.
And S4.5, performing crossover operation. The current population is paired, and the fitness is maximum
Figure 872464DEST_PATH_IMAGE061
Individuals did not participate in crossover, which was performed as follows:
Figure DEST_PATH_IMAGE063
(23)
in the formula (23), the number of distributed power supplies adopts integer intersection, and the electricity price adopts real number intersection. The improved crossover rate was as follows:
Figure 100002_DEST_PATH_IMAGE064
(24)
in the formula (24), T is the current population genetic algebra, T is the maximum genetic algebra, the cross rate is reduced along with the increase of the genetic algebra, and the convergence of the algorithm is accelerated. For each pair of genes, the roulette method was used, and if it was determined that the crossing was performed, the crossing was performed according to equation (23), and if the crossing was not performed, the alleles were directly introduced into the next generation.
And S4.6, mutation operation, wherein the individual elite does not participate in the mutation operation. Each gene is subjected to mutation operation by adopting a wheel disc method according to the following improved self-adaptive mutation rate:
Figure DEST_PATH_IMAGE065
(25)
in the formula (25), the reaction mixture,
Figure 100002_DEST_PATH_IMAGE066
for the purpose of the current average fitness measure,
Figure DEST_PATH_IMAGE067
in order to reach the current maximum fitness, the formula realizes that the variation rate is reduced along with the genetic algebra and the population average fitness is improved, so that the algorithm convergence is accelerated.
S4.7, if the maximum iteration number T is reached, ending the algorithm; if the number of iterations T is not the maximum number of iterations T, the process returns to step S4.1.
And step S5, selecting the optimal individual or individuals, acquiring the position and number information of the distributed power supply and the peak-valley bisection time interval electricity price, and calculating each objective function value to facilitate decision making.
The method is applied to the optimization configuration problem of the distributed power supply accessed to the power distribution network, meanwhile, the electricity price demand response of a user side is comprehensively considered, a demand response model of the user for electricity price change is established, so that a combined planning model of the distributed power supply and the electricity price decision is established, the multi-objective nonlinear optimization problem aiming at the minimum network loss, the minimum power grid load loss, the maximum comprehensive benefit, the minimum peak-valley difference and the like is considered, the solution is solved by improving a genetic algorithm, the optimization configuration scheme of the distributed power supply is obtained, and the peak-valley time-sharing electricity price is obtained at the same time.
The above-mentioned embodiments describe the present invention in further detail, and some parameters and functions are instantiated, so that equivalent substitutions can be made in practical applications, and suitable parameters can be selected according to specific situations.

Claims (5)

1. A distributed power supply grid-connected optimization configuration method for considering electricity price response in an incremental power distribution network is characterized by comprising the following steps:
step S1, acquiring the wind speed and the sunlight intensity of the incremental distribution network area
Figure DEST_PATH_IMAGE002
Establishing a wind power and photovoltaic output model by data distribution on a time interval full time sequence;
step S2, constructing each load node in the incremental distribution network
Figure DEST_PATH_IMAGE002A
The method comprises the steps that the time intervals are distributed in a full time sequence, the difference of different types of industries in response to power price change is considered, and a cross elastic matrix for acquiring the full time sequence of each time interval of nodes is constructed according to similar historical data of demand response of other power distribution networks to establish a load demand response model;
step S3, the power supply main body aims at minimizing the load peak-valley difference of the power distribution network, minimizing the total investment cost of the power distribution network and minimizing the network loss of the power distribution network, and the power selling main body aims at maximizing the power selling income and constructs a multi-objective optimization function;
step S4, improving the genetic algorithm to solve, according to a plurality of objective functions in the power supply main body, establishing the weight based on the variance, calculating the comprehensive fitness of the individual, synthesizing the external objective function, namely synthesizing the single objective function of the power selling main body, obtaining the final total fitness, and adopting the replication operator, the self-adaptive crossover operator and the mutation operator of the elite strategy to accelerate the convergence;
s5, decoding according to the result of the genetic algorithm to obtain the distributed power distribution point position and the configuration capacity;
of the nodes constructed in step S2
Figure DEST_PATH_IMAGE002AA
The elements of the cross elastic matrix of the full time sequence of each period are obtained by elastically weighting and summing requirements of different industries, so that a calculation formula of the requirement elastic coefficient of the current power distribution network area is obtained:
Figure DEST_PATH_IMAGE004
(1)
in the formula (1), the reaction mixture is,
Figure DEST_PATH_IMAGE006
at the basic electricity price for the k-th industry
Figure DEST_PATH_IMAGE008
Lower part
Figure DEST_PATH_IMAGE010
The weight of the load taken up by the time period,
Figure DEST_PATH_IMAGE012
for the purpose of the corresponding elastic coefficient,
Figure DEST_PATH_IMAGE014
the demand elastic coefficient of the electricity quantity of a certain node at the m moment to the electricity price change at the n moment is obtained;
the electric quantity change of a user full time sequence under a certain node of the incremental power distribution network is obtained by the following formula:
Figure DEST_PATH_IMAGE016
the electricity consumption of the user t period after the electricity price response is carried out is as follows:
Figure DEST_PATH_IMAGE018
Figure DEST_PATH_IMAGE020
in order to be the historical average electricity prices,
Figure DEST_PATH_IMAGE022
in order to correspond to the electric quantity demand of the m-th time period under the electricity price,
Figure DEST_PATH_IMAGE024
Figure DEST_PATH_IMAGE026
the electricity consumption change in m periods and the electricity price in n periods are respectively compared with
Figure DEST_PATH_IMAGE008A
The change that occurred;
the multi-objective optimization function is constructed in step S3 as follows:
Figure DEST_PATH_IMAGE028
(2)
Figure DEST_PATH_IMAGE030
(3)
the formula (2) is a multi-objective function considered by the power supply main body, wherein the first objective function is the minimum total loss of the network, the second objective function is the minimum peak-to-valley difference of the nodes accessed to the public network,
Figure DEST_PATH_IMAGE032
the third objective function is the minimum total investment of the power distribution network for a full time sequence of electric quantity, wherein the first two items are the total investment sum of wind power and photovoltaic,
Figure DEST_PATH_IMAGE034
for the cost of the generation of the unit active power of the wind power,
Figure DEST_PATH_IMAGE036
in order to reduce the power generation cost per unit area of the photovoltaic cell,
Figure DEST_PATH_IMAGE038
and
Figure DEST_PATH_IMAGE040
the output of a single fan and the output of a photovoltaic cell are respectively, the third term is the distributed generation benefit,
Figure DEST_PATH_IMAGE042
for the time period t, the price of selling electricity,
Figure DEST_PATH_IMAGE044
in order to increase the price of the power transmission and distribution of the power distribution network,
Figure DEST_PATH_IMAGE046
the electric quantity participating in market trading is distributed generation at the t time period; the formula (3) is an objective function considered by the electricity selling main body and is an objective function for maximizing the income of electricity selling in the whole network;
step S4 is carried out on the grid-connected optimal configuration of the distributed power supply, meanwhile, the demand response of the user to the peak-valley electricity price is considered, a comprehensive planning model is established, and the genetic algorithm gene sequence is as follows:
Figure DEST_PATH_IMAGE048
(4)
Figure DEST_PATH_IMAGE050
Figure DEST_PATH_IMAGE052
are respectively the first
Figure DEST_PATH_IMAGE054
The number of fans and photovoltaic devices installed on each node,
Figure DEST_PATH_IMAGE056
Figure DEST_PATH_IMAGE058
Figure DEST_PATH_IMAGE060
the electricity prices of the peak time period, the ordinary time period and the valley time period of the region are obtained;
step S4, optimizing the multi-target function in the power supply main body by applying a genetic algorithm, and determining the weight of each function by adopting the fitness variance contribution based on the target function in order to enlarge the search area of the Pareto optimal solution:
Figure DEST_PATH_IMAGE062
(5)
Figure DEST_PATH_IMAGE064
(6)
in the formula (5), the reaction mixture is,
Figure DEST_PATH_IMAGE066
is as follows
Figure DEST_PATH_IMAGE068
The weight of the respective objective function is,
Figure DEST_PATH_IMAGE070
is as follows
Figure DEST_PATH_IMAGE068A
The fitness of each objective function is calculated,
Figure DEST_PATH_IMAGE072
the number of the functions of a certain target is shown in the formula (6)
Figure DEST_PATH_IMAGE074
Is as follows
Figure DEST_PATH_IMAGE068AA
The fitness of the individual objective functions is,
Figure DEST_PATH_IMAGE076
the comprehensive fitness of the multi-objective function in the power supply main body is obtained.
2. The distributed power grid-connected optimization configuration method for the counting and electricity price response in the incremental power distribution network according to claim 1, characterized in that:
Figure DEST_PATH_IMAGE002AAA
the value is 24 or 96.
3. The distributed power grid-connected optimization configuration method for the neutralization electricity price response in the incremental power distribution network according to claim 1, wherein the optimization problem is solved according to the genetic algorithm applied in the step S4, a replication operator based on an elite strategy is adopted, and the maximum fitness is selected in each generation
Figure DEST_PATH_IMAGE078
Directly enters the next generation operation as the individual elite without mutation and cross operation.
4. The distributed power grid-connected optimization configuration method for counting and responding to electricity prices in incremental power distribution network according to claim 1, wherein the step S4 is executed according toSolving the optimization problem by applying a genetic algorithm, adopting a self-adaptive crossover operator, wherein T is the current population genetic algebra, T is the maximum genetic algebra,
Figure DEST_PATH_IMAGE080
in order to maximize the rate of cross-over,
Figure DEST_PATH_IMAGE082
minimum crossover rate:
Figure DEST_PATH_IMAGE084
5. the distributed power grid-connected optimization configuration method for the counting and electricity price response in the incremental power distribution network according to claim 1, wherein the optimization problem is solved according to the genetic algorithm applied in the step 4, an adaptive mutation operator is adopted,
Figure DEST_PATH_IMAGE086
Figure DEST_PATH_IMAGE088
for the purpose of the current average fitness measure,
Figure DEST_PATH_IMAGE090
in order to arrive at the current maximum fitness,
Figure DEST_PATH_IMAGE080A
and
Figure DEST_PATH_IMAGE092
the maximum and minimum mutation rates are respectively, T is the current population genetic algebra, and T is the maximum genetic algebra.
CN201810793769.1A 2018-07-19 2018-07-19 Distributed power supply grid-connected optimization configuration method for counting and electricity price response in incremental power distribution network Active CN108964124B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810793769.1A CN108964124B (en) 2018-07-19 2018-07-19 Distributed power supply grid-connected optimization configuration method for counting and electricity price response in incremental power distribution network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810793769.1A CN108964124B (en) 2018-07-19 2018-07-19 Distributed power supply grid-connected optimization configuration method for counting and electricity price response in incremental power distribution network

Publications (2)

Publication Number Publication Date
CN108964124A CN108964124A (en) 2018-12-07
CN108964124B true CN108964124B (en) 2021-03-16

Family

ID=64496328

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810793769.1A Active CN108964124B (en) 2018-07-19 2018-07-19 Distributed power supply grid-connected optimization configuration method for counting and electricity price response in incremental power distribution network

Country Status (1)

Country Link
CN (1) CN108964124B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109886524A (en) * 2018-12-26 2019-06-14 东南大学 A kind of edge cooperated computing method that distributed generation resource responds tou power price power generation
CN109861295B (en) * 2019-04-08 2021-06-08 国网河北省电力有限公司经济技术研究院 Multi-target new energy planning method and device
CN111259315B (en) * 2020-01-17 2022-04-19 福州大学 Decentralized scheduling method of multi-subject coordinated pricing mode
CN111967639B (en) * 2020-07-06 2023-04-07 国网河北省电力有限公司保定供电分公司 Power distribution network optimization method and system considering wind-solar power output and price demand response
CN112085291A (en) * 2020-09-18 2020-12-15 杭州市电力设计院有限公司余杭分公司 Genetic algorithm, system and device for planning and optimizing power distribution network
CN112510684B (en) * 2020-11-13 2023-10-13 安徽电力交易中心有限公司 Interruptible load configuration method for adjacent incremental power distribution network
CN117709554B (en) * 2024-02-04 2024-04-16 天津龙海科技有限公司 Energy scheduling method and system combining heat storage with electromagnetic heater

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7629705B2 (en) * 2006-10-20 2009-12-08 General Electric Company Method and apparatus for operating electrical machines
CN101847891A (en) * 2009-10-30 2010-09-29 深圳市科陆电子科技股份有限公司 Method for monitoring distributed power supply network incorporation in distribution network terminal
CN101834431A (en) * 2010-05-21 2010-09-15 深圳市科陆电子科技股份有限公司 Method for preventing grid-connected island effect of distributed power source
CN105071443A (en) * 2015-07-17 2015-11-18 上海电机学院 Distributed power supply grid connection control method
CN106300323B (en) * 2016-10-11 2019-03-08 耿天侃 Distributed generation resource power grid
CN106532769B (en) * 2016-11-25 2019-04-30 云南电网有限责任公司电力科学研究院 Optimization method based on grid-connected probabilistic Demand-side tou power price
CN107612016B (en) * 2017-08-08 2020-01-14 西安理工大学 Planning method of distributed power supply in power distribution network based on maximum voltage correlation entropy

Also Published As

Publication number Publication date
CN108964124A (en) 2018-12-07

Similar Documents

Publication Publication Date Title
CN108964124B (en) Distributed power supply grid-connected optimization configuration method for counting and electricity price response in incremental power distribution network
Wang et al. Incentivizing energy trading for interconnected microgrids
Diekerhof et al. Hierarchical distributed robust optimization for demand response services
Nojavan et al. Selling price determination by electricity retailer in the smart grid under demand side management in the presence of the electrolyser and fuel cell as hydrogen storage system
Yang et al. A multi-objective stochastic optimization model for electricity retailers with energy storage system considering uncertainty and demand response
Luo et al. Distributed peer-to-peer energy trading based on game theory in a community microgrid considering ownership complexity of distributed energy resources
Senemar et al. Dynamic structural sizing of residential energy hubs
CN110350523A (en) Multi-energy complementation Optimization Scheduling based on demand response
Lu et al. Interactive game vector: A stochastic operation-based pricing mechanism for smart energy systems with coupled-microgrids
Wang et al. Bargaining-based energy trading market for interconnected microgrids
Jordehi Two-stage stochastic programming for risk-aware scheduling of energy hubs participating in day-ahead and real-time electricity markets
CN112036934A (en) Quotation method for participation of load aggregators in demand response considering thermoelectric coordinated operation
CN112800658A (en) Active power distribution network scheduling method considering source storage load interaction
Qiu et al. Mean-field multi-agent reinforcement learning for peer-to-peer multi-energy trading
Amicarelli et al. Optimization algorithm for microgrids day-ahead scheduling and aggregator proposal
Dahlin et al. Scheduling flexible nonpreemptive loads in smart-grid networks
CN110571795A (en) arrangement method of energy storage unit in high-wind-force penetration power system
Dorahaki et al. A sharing economy model for a sustainable community energy storage considering end-user comfort
Hayati et al. A two-stage stochastic optimization scheduling approach for integrating renewable energy sources and deferrable demand in the spinning reserve market
Chreim et al. Recent sizing, placement, and management techniques for individual and shared battery energy storage systems in residential areas: A review
Nammouchi et al. Robust opportunistic optimal energy management of a mixed microgrid under asymmetrical uncertainties
CN113988471A (en) Multi-objective optimization method for micro-grid operation
Kong et al. Independence enhancement of distributed generation systems by integrating shared energy storage system and energy community with internal market
CN112598175A (en) Watershed type virtual power plant scheduling method
Kang et al. Multi-objective sizing and real-time scheduling of battery energy storage in energy-sharing community based on reinforcement learning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant