CN115545501A - Multi-mode multi-target layered distributed comprehensive energy system economic dispatching method - Google Patents
Multi-mode multi-target layered distributed comprehensive energy system economic dispatching method Download PDFInfo
- Publication number
- CN115545501A CN115545501A CN202211257848.3A CN202211257848A CN115545501A CN 115545501 A CN115545501 A CN 115545501A CN 202211257848 A CN202211257848 A CN 202211257848A CN 115545501 A CN115545501 A CN 115545501A
- Authority
- CN
- China
- Prior art keywords
- layer
- region
- representing
- active
- formula
- 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 38
- 239000003795 chemical substances by application Substances 0.000 claims description 59
- 238000010248 power generation Methods 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000010586 diagram Methods 0.000 claims description 11
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 10
- 229910052799 carbon Inorganic materials 0.000 claims description 10
- 238000012216 screening Methods 0.000 claims description 10
- 239000000126 substance Substances 0.000 claims description 8
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 6
- 229910052739 hydrogen Inorganic materials 0.000 claims description 6
- 239000001257 hydrogen Substances 0.000 claims description 6
- 230000009194 climbing Effects 0.000 claims description 4
- 229910052731 fluorine Inorganic materials 0.000 claims description 4
- 125000001153 fluoro group Chemical group F* 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000012986 modification Methods 0.000 claims description 2
- 230000004048 modification Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides an economic dispatching method of a multi-mode multi-target layered distributed comprehensive energy system. Firstly, an analytic multi-mode multi-target method in the method is used for obtaining a multi-mode economic dispatching multi-target scheme with diversity. Secondly, the hierarchical distributed consistency method in the method utilizes hierarchical operation to quickly obtain an accurate economic dispatching decision scheme. The multi-mode multi-target layered distributed economic dispatching method for the comprehensive energy system can solve the problem of multi-mode economic dispatching of the large-scale comprehensive energy system, diversifies the decision scheme of economic dispatching, optimizes the stability of economic dispatching and improves the calculation speed and accuracy.
Description
Technical Field
The invention belongs to the field of power systems and comprehensive energy, and relates to an analytic multi-mode multi-target layered distributed method which is suitable for multi-mode multi-target layered distributed economic dispatching of a comprehensive energy system.
Background
The continuous expansion of the power market and the gradual expansion of the scale of the comprehensive energy system can cause the calculation speed of the economic dispatch of the comprehensive energy system to be slow and the information privacy to be poor. The economic dispatching of the comprehensive energy system is mostly aimed at a single target, and a multi-target method for the economic dispatching of the comprehensive energy is few. And moreover, multi-mode problems are not considered in economic dispatch in the research of the multi-target economic dispatch problem of the comprehensive energy system, so that the multi-target economic dispatch scheme lacks diversity and comprehensiveness.
Therefore, the method can keep the characteristics of multiple modes, can accelerate the speed of system processing problems and ensure the privacy and the robustness of the system in the multi-mode multi-target layered distributed economic dispatching problem of the comprehensive energy system, can also provide multi-mode selection for the dispatching of the comprehensive energy system, and ensures the stability and the safety of the economic dispatching.
Disclosure of Invention
The invention provides an economic dispatching method of a multi-mode multi-target layered distributed comprehensive energy system, which combines an analytic multi-mode multi-target method and a layered distributed consistency method for economic dispatching of the comprehensive energy system, increases the diversity and stability of economic dispatching, improves the calculation speed and precision, and provides rich decision schemes which can be replaced for decision makers; the steps in the using process are as follows:
step (1): the construction of an economic dispatching model of the comprehensive energy system comprises fire energy, water energy, wind energy, light energy and geothermal energy, and conforms to equality constraint and inequality constraint of economic dispatching by taking cost expense and carbon emission as multiple targets;
the objective function for the total power generation cost is:
wherein, the first and the second end of the pipe are connected with each other,total cost of power generation;representing the cost of the ith generating set at time t;representing the cost of the jth thermal power generating unit at the moment t;representing the cost of the kth hydroelectric generating set at time t;representing the cost of the first wind generating set at the time t;representing the cost of the r-th photovoltaic generator set at the time t;representing the cost of the s th geothermal energy generating set at the time t;representing the cost of the uth clean energy generator set at the time t; n is the number of the total generator sets; n is a radical of hydrogen TG The number of thermal power generator units; n is a radical of hydrogen HG The number of hydroelectric generator sets; n is a radical of hydrogen WG The number of wind driven generator sets; n is a radical of PG The number of photovoltaic generator sets; n is a radical of OG The number of geothermal energy generator sets; n is a radical of CleanG The number of generator sets is the clean energy; t is the statistical time of the objective function;for the ith generating set at the moment tAn amount of electricity;generating capacity of the jth thermal power generating set at the moment t;the generated energy of the kth hydroelectric generating set at the moment t;the generated energy of the first wind generating set at the moment t;generating capacity of the r-th photovoltaic generator set at the moment t;generating capacity of the s th geothermal energy generating set at the time t;generating capacity of the u clean energy generator set at the moment t; a is a j 、b j And c j The cost coefficient is the jth thermal power generating set; d k The cost coefficient of the kth hydroelectric generating set; e.g. of the type l The cost coefficient of the first wind generating set; g r The cost coefficient of the r photovoltaic generator set is set; m is s The cost coefficient of the s th geothermal energy generating set; kappa type u The cost coefficient of the u clean energy generator set is; cost of unit using clean energyReplacement of hydroelectric costsCost of wind powerCost of light energyAnd cost of geothermal energy
The objective function for total carbon emissions is:
wherein the content of the first and second substances,total carbon emissions;representing the carbon emission of the jth thermal power generating unit at the time t; alpha (alpha) ("alpha") j 、β j And gamma j The carbon emission coefficient of the jth thermal power generating unit is obtained; the equality constraint for power balancing is:
wherein, the first and the second end of the pipe are connected with each other,the total load demand value at time t;
the upper and lower limits of power are constrained as follows:
wherein the content of the first and second substances,the lower limit of the power generation amount of the jth thermal power generator set is set;the upper limit of the power generation amount of the jth thermal power generating set is set;the lower limit of the generating capacity of the kth hydroelectric generating set;the upper limit of the generating capacity of the kth hydroelectric generating set;the lower limit of the generating capacity of the first wind generating set;the upper limit of the generating capacity of the first wind generating set;the lower limit of the generating capacity of the r-th photovoltaic generator set;the upper limit of the generating capacity of the r-th photovoltaic generator set is set;the lower limit of the generating capacity of the s th geothermal energy generating set;the upper limit of the power generation amount of the s th geothermal energy power generation unit;
the ramp rate constraint of the thermal power generating unit is as follows:
wherein the content of the first and second substances,the power generation amount of the jth thermal power generator set at the moment t-1;the downward climbing speed of the jth thermal power generating unit is obtained;the upward climbing speed of the jth thermal power generating unit is obtained; t is 60 60 minutes;
the time constraints of the photovoltaic generator set are:
Wherein, R is the ratio of the values of the weighting factors; eta is a weight factor, and the numeric area of eta is more than or equal to 0 and less than or equal to 1; n is η Representing the number of segments for segmenting the value range of the weight factor;
and (3): converting the multi-target problem into a single-target problem by a linear weighting method;
wherein f is Total Is the total target value after linear weighting;a weighted target value is obtained for the jth thermal power unit line at the moment t;
and (4): setting an initial iteration time t =1;
and (5): inputting a predicted load value, and enabling the initial iteration step number k =0;
and (6): the comprehensive energy system is divided into four layers: the first layer is divided into an area A, an area B and an area C; the second layer is divided into a region Aa, a region Ab, a region Ac, a region Ba, a region Bb, a region Bc, a region Ca and a region Cb; the third layer is divided into a region Aa1, a region Aa2, a region Ab1, a region Ab2, a region Ac1, a region Ac2, a region Ba1, a region Ba2, a region Ba3, a region Bb1, a region Bb2, a region Bc1, a region Bc2, a region Ca1, a region Ca2, a region Cb1, a region Cb2, and a region Cb3; the fourth layer is a generator set in each area of the third layer; and (7): selecting leader and follower agents of each layer, wherein an area A of a first layer is a leader, and an area B and an area C are followers; the area Aa, the area Ba and the area Ca of the second layer are taken as leaders, and the other two-layer areas are taken as followers; the third layer area Aa1, the area Ab1, the area Ac1, the area Ba1, the area Bb1, the area Bc1, the area Ca1, and the area Cb1 are leaders, and the other third layer areas are followers; and (8): updating the consistency variable of the first layer of agents by combining the topological structure diagram of the agents and the formula (15) and the formula (16);
directed graph G = (V, E) of distributed network topology and adjacency matrix a = (a) of weight relationships between respective agents ij ) Forming a network topology structure of the multi-agent; laplace matrix L = [ L ] of directed graph G ij ]Comprises the following steps:
wherein l ii And l ij Is an element in the Laplace matrix L; a is ij Is an element in adjacency matrix a;
row random matrix H of agent ij Comprises the following steps:
the consistency variable λ is expressed as:
λ=[λ T ,λ Clean ] (12)
wherein λ represents a consistency variable of the integrated energy system; lambda [ alpha ] T Representing a consistency variable of the thermal generator set; lambda [ alpha ] Clean Representing a consistency variable of the clean energy generator set;
the consistency variable updating formulas of the follower and the leader are respectively as follows:
wherein N is IB Representing the number of agents; k represents the kth iteration;a consistency variable representing the nth agent at the kth iteration at the time t;a consistency variable representing the mth agent at the time t in the (k + 1) th iteration; τ represents a power balance adjustment factor of the distributed consistency method; delta P t (k) Representing the power deviation of the kth iteration at time t; h mn (k) A row random matrix representing an agent;
and (9): calculating the active output of each agent of the first layer according to a formula (17);
the active power formula of each unit is as follows:
wherein, the first and the second end of the pipe are connected with each other,the active output of the ith agent at the moment t in k +1 iteration is represented;the active power output of the jth thermal power generating unit at the moment t in k +1 iterations is represented;the active power output of the u-th clean energy generator set at the moment t in k +1 iterations is represented;the active upper limit of the mu region of the omega layer is expressed as the sum of the active upper limits of the lower sub-regions or the units divided by the region;representing a consistency variable of the ith agent at the moment t in k +1 iterations;representing a consistency variable of the jth thermal power generating unit at the moment t in k +1 iterations;representing a consistency variable of the u clean energy unit at the moment t in k +1 iterations;representing the number of thermal power generating units in the mu region of the omega layer;representing the number of clean energy generator sets of the mu region of the omega layer; u represents the total number of regions of the layer; w represents the total number of layers;
step (10): correcting the first layer active power output obtained in the step (9) according to a formula (18);
the active power correction formula is as follows:
wherein, the first and the second end of the pipe are connected with each other,representing the lower limit of the active output of the ith power generation intelligent agent;the upper limit of the active output of the ith power generation intelligent agent is represented;
step (11): solving an active power deviation value of the first layer according to a formula (19);
active power deviation Δ P t (k + 1) is:
step (12): judging whether the multi-mode active power deviation value of the first layer obtained in the step (11) meets the precision requirement or not
Wherein δ represents a power deviation allowable maximum value;representing a multi-modal power deviation allowance value;
if it is notEnabling the iteration step number k = k +1, returning to the step (8), and continuing to perform the iterative computation of the first layer;
if it is notStoring the obtained active output of the group of the intelligent agents, then performing another iteration step number k = k +1, returning to the step (8), and continuing to perform the iteration calculation of the first layer;
if it is notThe multi-mode active deviation value of the first layer meets the requirement, iteration is finished, and active outputs of a plurality of groups of intelligent bodies are obtained;
step (13): performing multi-mode judgment and decision congestion distance screening on the active power output of the multiple groups of intelligent generator sets obtained in the step (12) to obtain multi-mode active power output meeting the requirements of decision makers; the multi-modal judgment is as follows:
|f(x 1 )-f(x 2 )|<ε (21)
wherein x is 1 And x 2 Two sets of solutions for multimodal; epsilon represents the critical threshold for multiple modes;
the decision space crowding distance screening is as follows:
wherein the content of the first and second substances,a crowding distance representing a jth decision variable of the ith set of solutions; x is a radical of a fluorine atom i+1,j A jth decision variable representing the i +1 th set of solutions; x is a radical of a fluorine atom i,j A jth decision variable representing the ith set of solutions;a jth decision variable representing a maximum target value;a jth decision variable representing a minimum target value; x is the number of i-1,j A jth decision variable representing the i-1 th set of solutions;indicating the congestion distance of the i-th group solution; d represents the number of decision variables; ξ represents the minimum decision crowding distance allowed by multimodal;
step (14): inputting the multiple groups of multi-modal generated energy of the first layer obtained in the step (13) as load values into a second layer, and enabling the initial iteration step number k =0;
step (15): updating the consistency variable of the agent at the second layer by combining the topological structure diagram of the agent and the formula (15) and the formula (16);
step (16): working out the active output of each intelligent unit of the second layer according to a formula (17);
step (17): correcting the active power output obtained in the step (16) according to a formula (18);
step (18): solving the deviation value of the active power of the second layer according to a formula (19);
step (19): judging whether the multi-mode active power deviation value obtained in the step (18) meets the precision requirement or not;
if it is notMaking the iteration step number k = k +1, returning to the step (15), and continuing to perform the iterative computation of the second layer;
if it is notStoring the obtained active output of the group of the intelligent agents, then performing another iteration step number k = k +1, returning to the step (15), and continuing to perform iterative computation of the second layer;
if it is notThe multi-mode active deviation value of the second layer meets the requirement, iteration is finished, and active output of a plurality of groups of intelligent generator sets is obtained;
step (20): performing multi-mode judgment and decision congestion distance screening on the active output of the multiple groups of agents on the second layer obtained in the step (19) to obtain multi-mode active output meeting the requirements of a decision maker;
step (21): inputting the multi-group multi-modal power generation amount of the second layer obtained in the step (20) as a load value to obtain three layers, and enabling the initial iteration step number k =0;
step (22): updating the consistency variable of the agent at the third layer by combining the topological structure diagram of the agent and the formula (15) and the formula (16);
step (23): working out the active output of each intelligent machine set of the third layer according to a formula (17);
step (24): correcting the active power output obtained in the step (23) according to a formula (18);
step (25): solving the deviation value of the active power of the third layer according to a formula (19);
step (26): judging whether the multi-mode active power deviation value obtained in the step (25) meets the precision requirement or not;
if it is usedMaking the iteration step number k = k +1, returning to the step (22), and continuing to perform the iterative computation of the third layer;
if it is notStoring the obtained active output of the group of the intelligent agents, then performing another iteration step number k = k +1, returning to the step (22), and continuing to perform iterative computation of a third layer;
if it is notThen the multi-modal active deviation of the third layerThe value meets the requirement, iteration is finished, and the active output of a plurality of groups of intelligent generator sets is obtained;
step (27): performing multi-mode judgment and decision congestion distance screening on the active output of the multiple groups of agents on the third layer obtained in the step (26) to obtain multi-mode active output meeting the requirements of a decision maker;
step (28): inputting the multi-group multi-modal power generation amount of the third layer obtained in the step (27) as a load value to obtain four layers, and enabling the initial iteration step number k =0;
step (29): updating the consistency variable of each generator set on the fourth layer by combining the topological structure diagram of the agent and the formula (15) and the formula (16);
a step (30): working out the active power output of each generator set on the fourth layer according to a formula (17);
step (31): correcting the active power output obtained in the step (30) according to a formula (18); step (32): solving the deviation value of the active power of the fourth layer according to a formula (19);
step (33): judging whether the multi-mode active power deviation value obtained in the step (32) meets the precision requirement or not;
if it is usedMaking the iteration step number k = k +1, returning to the step (29), and continuing to perform the iterative computation of the fourth layer;
if it is usedStoring the obtained active power output of the group of generator sets, then performing another iteration step number k = k +1, returning to the step (29), and continuing to perform iterative calculation of the fourth layer;
if it is notThe multi-mode active deviation value of the fourth layer meets the requirement, iteration is finished, and active output of multiple groups of generator sets is obtained;
step (34): performing multi-mode judgment and decision congestion distance screening on the active power output of the multiple groups of generator sets on the fourth layer obtained in the step (33) to obtain a multi-mode active power output decision scheme meeting the requirements of a decision maker;
step (35): judging whether T < T is met, if so, enabling T = T +1, and turning to the step (5); if not, executing the next step;
step (36): judging whether eta <1 is satisfied, if yes, storing a multi-modal solution obtained when eta is taken as value, enabling eta = eta + R, and turning to the step (4); if not, executing the next step;
step (37): and reversely solving all the multi-mode solutions obtained by taking different eta values to obtain multi-mode multi-target values, wherein the multi-mode multi-target values obtained under all the eta values form a pareto frontier in a target space.
Compared with the prior art, the invention has the following advantages and effects:
(1) An analytic framework of the multi-mode multi-target layered distribution method is established for the economic dispatching problem of the comprehensive energy system.
(2) Because the constructed analytic multi-mode multi-target layered distribution method adopts layered scheduling and only exchanges the consistency variables of adjacent areas, the detailed operation information of each unit does not need to be exchanged, so that the calculation speed is accelerated and the information privacy is better guaranteed; and the problem of multi-modal characteristics in the dispatching of the comprehensive energy system is solved, and the diversity and comprehensiveness of an economic dispatching solution are increased.
(3) The obtained economic dispatching scheme of the comprehensive energy system has a plurality of mutually-replaced economic dispatching schemes, and when some sudden accidents or serious wind, light and water abandonment in partial areas are faced, decision makers have more alternative schemes, so that the stability and the safety of economic dispatching are improved.
Drawings
FIG. 1 is a schematic diagram of the hierarchical distribution of the process of the present invention.
FIG. 2 is a flow chart of an analytical multi-modal multi-target distributed consistency method of the present invention.
FIG. 3 is an overall flow chart of the method of the present invention.
Detailed Description
The invention provides an economic dispatching method of a multi-mode multi-target layered distributed comprehensive energy system, which is explained in detail by combining the attached drawings as follows:
FIG. 1 is a schematic diagram of the hierarchical distribution of the process of the present invention. The comprehensive energy system is divided into four layers: the first layer is divided into an area A, an area B and an area C; the second layer is divided into a region Aa, a region Ab, a region Ac, a region Ba, a region Bb, a region Bc, a region Ca and a region Cb; the third layer is divided into a region Aa1, a region Aa2, a region Ab1, a region Ab2, a region Ac1, a region Ac2, a region Ba1, a region Ba2, a region Ba3, a region Bb1, a region Bb2, a region Bc1, a region Bc2, a region Ca1, a region Ca2, a region Cb1, a region Cb2, and a region Cb3; the fourth layer is a generator set in each area of the third layer. Each layer elects a leader and a follower.
FIG. 2 is a flow chart of an analytical multi-modal multi-target distributed consistency method of the present invention. The method comprises the following specific steps:
step 1: converting the multi-target problem into a single-target problem according to a formula (9);
step 2: inputting a predicted load value, and enabling the initial iteration step number k =0;
and step 3: determining a leader and a follower of the agent, and updating a consistency variable of the agent by combining a topological structure diagram of the agent and a formula (15) and a formula (16);
and 4, step 4: working out the active output of each intelligent unit according to a formula (17);
and 5: correcting the active power output obtained in the step 4 according to a formula (18);
and 6: solving an active power deviation value according to a formula (19);
and 7: judging whether the multi-mode active power deviation value obtained in the step 6 meets the precision requirement, if so, judging whether the multi-mode active power deviation value meets the precision requirementEnabling the iteration step number k = k +1, and returning to the step 3; if it is usedStoring the obtained active output of the group of the intelligent agents, then iterating the step number k = k +1, and returning to the step 3; if it is notThe multi-mode active deviation value meets the requirement, iteration is finished, and active output of a plurality of groups of intelligent generator sets is obtained;
and step 8: and (4) performing multi-mode judgment and decision congestion distance screening on the active output of the plurality of groups of intelligent generator sets obtained in the step (7) to obtain a multi-mode active output scheduling scheme meeting the requirements of a decision maker.
Fig. 3 is an overall flow chart of the method of the present invention. The method comprises the following specific steps:
step 1: initializing the weighting factor, let eta =1/b η ;
Step 2: setting an initial iteration time t =1; and step 3: inputting a predicted load value, and enabling the initial iteration step number k =0;
and 4, step 4: dividing the system into four layers, wherein the first layer is divided into an area A, an area B and an area C; the second layer is divided into a region Aa, a region Ab, a region Ac, a region Ba, a region Bb, a region BC, a region Ca, and a region Cb; the third layer is divided into a region Aa1, a region Aa2, a region Ab1, a region Ab2, a region Ac1, a region Ac2, a region Ba1, a region Ba2, a region Ba3, a region Bb1, a region Bb2, a region Bc1, a region Bc2, a region Ca1, a region Ca2, a region Cb1, a region Cb2, and a region Cb3; the fourth layer is a generator set in each area of the third layer;
and 5: selecting leader and follower agents of each layer, wherein an area A of a first layer is a leader, and an area B and an area C are followers; the area Aa, the area Ba and the area Ca of the second layer are taken as leaders, and the other two-layer areas are taken as followers; the third layer area Aa1, the area Ab1, the area Ac1, the area Ba1, the area Bb1, the area Bc1, the area Ca1, and the area Cb1 are leaders, and the other third layer areas are followers;
step 6: carrying out optimization solution on the first layer of the system by using an analytic multi-mode multi-target distributed consistency method in FIG. 2;
and 7: inputting the multi-group multi-modal power generation quantities of the 3 areas of the first layer obtained in the step 6 into the second layer as load values respectively;
and 8: performing optimization solution on a second layer of the system by using the analytic multi-mode multi-target distributed consistency method in FIG. 2;
and step 9: inputting the multi-group multi-modal power generation quantities of the 8 areas of the second layer obtained in the step 8 into the third layer as load values respectively;
step 10: carrying out optimization solution on the third layer of the system by using the analytic multi-mode multi-target distributed consistency method in the figure 2;
step 11: inputting the multi-group multi-modal power generation quantities of the 18 areas in the third layer obtained in the step 10 into the fourth layer as load values, regarding each power generation unit in the fourth layer as an agent, and selecting a leader agent and a follower agent;
step 12: performing optimization solution on the fourth layer of the system by using the analytic multi-mode multi-target distributed consistency method in the figure 2;
step 13: judging whether T < T is met, if so, enabling T = T +1, and turning to the step 2; if not, executing the next step;
step 14: judging whether eta <1 is satisfied, if yes, storing a multi-modal solution obtained when eta is taken as value, enabling eta = eta + R, and turning to the step 2; if not, executing the next step; step 15: obtaining multi-modal multi-target values by reversely solving all multi-modal solutions obtained by taking different eta values; step 16: and forming pareto frontier in the target space by the multi-modal multi-target values obtained under all the eta values.
The above description is only a preferred embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by using the contents of the present specification and the accompanying drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.
Claims (1)
1. The economic dispatching method of the multi-mode multi-target layered distributed comprehensive energy system is characterized in that an analytic multi-mode multi-target method and a layered distributed consistency method are combined and used for economic dispatching of the comprehensive energy system, diversity and stability of economic dispatching are improved, computing speed and precision are improved, and a rich decision scheme capable of being replaced is provided for a decision maker; the steps in the using process are as follows:
step (1): the construction of an economic dispatching model of the comprehensive energy system comprises fire energy, water energy, wind energy, light energy and geothermal energy, and conforms to equality constraint and inequality constraint of economic dispatching by taking cost expense and carbon emission as multiple targets;
the objective function for the total power generation cost is:
wherein the content of the first and second substances,total cost of power generation;representing the cost of the ith generating set at time t;representing the cost of the jth thermal power generating unit at the moment t;representing the cost of the kth hydroelectric generating set at time t;indicates the first wind power generation at the time tThe cost of the unit;representing the cost of the r-th photovoltaic generator set at the time t;representing the cost of the s th geothermal energy generating set at the time t;representing the cost of the uth clean energy generator set at the time t; n is the number of the total generator sets; n is a radical of TG The number of thermal power generator units; n is a radical of HG The number of hydroelectric generator sets; n is a radical of WG The number of wind turbine generator sets; n is a radical of hydrogen PG The number of photovoltaic generator sets; n is a radical of hydrogen OG The number of geothermal energy generator sets; n is a radical of hydrogen CleanG The number of generator sets for clean energy; t is the statistical time of the objective function;generating capacity of an ith generating set at the moment t;generating capacity of the jth thermal power generating set at the moment t;the generating capacity of the kth hydroelectric generating set at the moment t;the generated energy of the first wind generating set at the moment t;generating capacity of the r-th photovoltaic generator set at the time t;generating capacity of the s th geothermal energy generating set at the time t;generating capacity of the u clean energy generator set at the moment t; a is j 、b j And c j The cost coefficient is the jth thermal power generating set; d is a radical of k The cost coefficient of the kth hydroelectric generating set; e.g. of a cylinder l The cost coefficient of the first wind generating set; g r The cost coefficient of the r-th photovoltaic generator set is obtained; m is s The cost coefficient of the s th geothermal energy generating set; k is a radical of u The cost coefficient of the u clean energy generator set is; cost of unit using clean energyReplacement of hydroelectric costsCost of wind powerCost of light energyAnd cost of geothermal energy
The objective function for total carbon emissions is:
wherein the content of the first and second substances,total carbon emissions;representing the carbon emission of the jth thermal power generating unit at the t moment; alpha (alpha) ("alpha") j 、β j And gamma j The carbon emission coefficient of the jth thermal power generating unit is obtained;
the equality constraint for power balancing is:
wherein, the first and the second end of the pipe are connected with each other,the total load demand value at the moment t;
the upper and lower limits of power are constrained as follows:
wherein, the first and the second end of the pipe are connected with each other,the lower limit of the power generation amount of the jth thermal power generator set is set;the upper limit of the generating capacity of the jth thermal power generating unit is set;the lower limit of the generating capacity of the kth hydroelectric generating set;the upper limit of the generating capacity of the kth hydroelectric generating set;for the first wind-driven generator setA lower limit of the electric quantity;the upper limit of the generating capacity of the first wind generating set;the lower limit of the generating capacity of the r-th photovoltaic generator set;the upper limit of the generating capacity of the r-th photovoltaic generator set is set;the lower limit of the generating capacity of the s th geothermal energy generating set;the upper limit of the power generation amount of the s th geothermal energy generating set;
the ramp rate constraint of the thermal power generating unit is as follows:
wherein, the first and the second end of the pipe are connected with each other,the power generation amount of the jth thermal power generator set at the moment t-1;the downward climbing speed of the jth thermal power generating unit is obtained;the upward climbing speed of the jth thermal power generating unit is obtained; t is a unit of 60 60 minutes;
the time constraints of the photovoltaic generator set are:
Wherein, R is the ratio of the values of the weighting factors; eta is a weight factor, and the numeric area of eta is more than or equal to 0 and less than or equal to 1; n is η The number of the segments for segmenting the value range of the weight factor is represented;
and (3): converting the multi-target problem into a single-target problem by a linear weighting method;
wherein, f Total Is the total target value after linear weighting;a weighted target value for the jth thermal power unit line at the time t;
and (4): setting an initial iteration time t =1;
and (5): inputting a predicted load value, and enabling the initial iteration step number k =0;
and (6): the comprehensive energy system is divided into four layers: the first layer is divided into an area A, an area B and an area C; the second layer is divided into a region Aa, a region Ab, a region Ac, a region Ba, a region Bb, a region Bc, a region Ca and a region Cb; the third layer is divided into a region Aa1, a region Aa2, a region Ab1, a region Ab2, a region Ac1, a region Ac2, a region Ba1, a region Ba2, a region Ba3, a region Bb1, a region Bb2, a region Bc1, a region Bc2, a region Ca1, a region Ca2, a region Cb1, a region Cb2, and a region Cb3; the fourth layer is a generator set in each area of the third layer; and (7): selecting leader and follower agents of each layer, wherein an area A of a first layer is a leader, and an area B and an area C are followers; the area Aa, the area Ba and the area Ca of the second layer are taken as leaders, and the other two-layer areas are taken as followers; a third layer region Aa1, a region Ab1, a region Ac1, a region Ba1, a region Bb1, a region Bc1, a region Ca1, and a region Cb1 are leaders, and the other third layer regions are followers; and (8): updating the consistency variable of the first layer of agents by combining the topological structure diagram of the agents and the formula (15) and the formula (16);
directed graph G = (V, E) of distributed network topologies and adjacency matrix a = (a) of weight relationships between respective agents ij ) Forming a network topology structure of the multi-agent; laplace matrix L = [ L ] of directed graph G ij ]Comprises the following steps:
wherein l ii And l ij Is an element in the Laplace matrix L; a is a ij Is an element in the adjacency matrix a;
row random matrix H of agent ij Comprises the following steps:
the consistency variable λ is expressed as:
λ=[λ T ,λ Clean ] (12)
wherein λ represents a consistency variable of the integrated energy system; lambda [ alpha ] T Representing a consistency variable of the thermal generator set; lambda Clean Representing a consistency variable of the clean energy generator set;
the consistency variable updating formulas of the follower and the leader are respectively as follows:
wherein N is IB Representing the number of agents; k represents the kth iteration;a consistency variable representing the nth agent at the kth iteration at the time t;representing the consistency variable of the mth agent at the time t in the (k + 1) th iteration; τ represents a power balance adjustment factor for the distributed consistency approach; delta P t (k) Representing the power deviation of the kth iteration at time t; h mn (k) A row random matrix representing an agent;
and (9): calculating the active output of each agent of the first layer according to a formula (17);
the active power formula of each unit is as follows:
wherein, the first and the second end of the pipe are connected with each other,denotes the ith agent at time t is at k +1 timesAn iterative active power output;the active power output of the jth thermal power generating unit at the moment t in k +1 iterations is represented;the active power output of the u clean energy generator set at the moment t in k +1 iterations is represented;the active upper limit of the mu region of the omega layer is expressed as the sum of the active upper limits of the lower sub-regions or the units divided by the region;representing a consistency variable of the ith agent at the moment t in k +1 iterations;representing a consistency variable of the jth thermal power generating unit at the moment t in k +1 iterations;representing a consistency variable of the u clean energy unit at the moment t in k +1 iterations;representing the number of thermal power generating units in the mu region of the omega layer;representing the number of clean energy generator sets of the mu region of the omega layer; u represents the total number of regions of the layer; w represents the total number of layers;
step (10): correcting the first layer active power output obtained in the step (9) according to a formula (18);
the active power modification formula is as follows:
wherein the content of the first and second substances,representing the lower limit of the active output of the ith power generation intelligent agent;representing the upper limit of the active output of the ith power generation agent;
step (11): solving an active power deviation value of the first layer according to a formula (19);
active power deviation Δ P t (k + 1) is:
step (12): judging whether the multi-mode active power deviation value of the first layer obtained in the step (11) meets the precision requirement or not
Wherein δ represents a power deviation allowable maximum value;representing a multi-modal power deviation allowance value;
if it is notEnabling the iteration step number k = k +1, returning to the step (8), and continuing to perform the iterative computation of the first layer;
if it is notStoring the obtained active output of the group of the intelligent agents, then performing another iteration step number k = k +1, returning to the step (8), and continuing to perform the iteration calculation of the first layer;
if it is usedThe multi-mode active deviation value of the first layer meets the requirement, iteration is finished, and active outputs of multiple groups of intelligent agents are obtained;
step (13): performing multi-mode judgment and decision congestion distance screening on the active power output of the multiple groups of intelligent generator sets obtained in the step (12) to obtain multi-mode active power output meeting the requirements of decision makers;
the multi-modal judgment is as follows:
|f(x 1 )-f(x 2 )|<ε (21) wherein x 1 And x 2 Two sets of solutions for multi-modal; epsilon represents the critical threshold for multiple modes;
the selection of the space crowding distance is as follows:
wherein the content of the first and second substances,a crowding distance representing a jth decision variable of the ith set of solutions; x is the number of i+1,j A jth decision variable representing the i +1 th set of solutions; x is a radical of a fluorine atom i,j A jth decision variable representing the ith set of solutions;a j-th decision variable representing a maximum target value;a j-th decision variable representing a minimum target value; x is a radical of a fluorine atom i-1,j A jth decision variable representing the i-1 th set of solutions;indicating the congestion distance of the i-th group solution; d represents the number of decision variables; ξ represents the minimum decision crowding distance allowed by multimodal;
step (14): inputting the multiple groups of multi-modal generated energy of the first layer obtained in the step (13) as load values into a second layer, and enabling the initial iteration step number k =0;
step (15): updating the consistency variable of the agent at the second layer by combining the topological structure diagram of the agent and the formula (15) and the formula (16);
step (16): working out the active output of each intelligent machine set of the second layer according to a formula (17);
step (17): correcting the active power output obtained in the step (16) according to a formula (18);
step (18): solving the deviation value of the active power of the second layer according to a formula (19);
step (19): judging whether the multi-mode active power deviation value obtained in the step (18) meets the precision requirement or not;
if it is usedEnabling the iteration step number k = k +1, returning to the step (15), and continuing to perform the iterative computation of the second layer;
if it is usedStoring the obtained active output of the group of the intelligent agents, then performing another iteration step number k = k +1, returning to the step (15), and continuing to perform iterative computation of the second layer;
if it is usedThe multi-mode active deviation value of the second layer meets the requirement, iteration is finished, and active output of multiple groups of intelligent generator sets is obtained;
step (20): performing multi-mode judgment and decision congestion distance screening on the active output of the multiple groups of agents on the second layer obtained in the step (19) to obtain multi-mode active output meeting the requirements of a decision maker;
step (21): inputting the multi-group multi-modal generating capacity of the second layer obtained in the step (20) as a load value to obtain three layers, and enabling the initial iteration step number k =0;
step (22): updating the consistency variable of the third layer of agents by combining the topological structure diagram of the agents and the formula (15) and the formula (16);
step (23): working out the active output of each intelligent unit of the third layer according to a formula (17);
step (24): correcting the active power output obtained in the step (23) according to a formula (18);
step (25): solving the deviation value of the active power of the third layer according to a formula (19);
step (26): judging whether the multi-mode active power deviation value obtained in the step (25) meets the precision requirement or not;
if it is notMaking the iteration step number k = k +1, returning to the step (22), and continuing to perform the iterative computation of the third layer;
if it is notStoring the obtained active output of the group of the intelligent agents, then performing another iteration step number k = k +1, returning to the step (22), and continuing to perform the iteration calculation of the third layer;
if it is notThe multi-mode active deviation value of the third layer meets the requirement, iteration is finished, and active output of a plurality of groups of intelligent generator sets is obtained;
step (27): performing multi-mode judgment and decision congestion distance screening on the active output of the multiple groups of intelligent agents on the third layer obtained in the step (26) to obtain multi-mode active output meeting the requirements of a decision maker;
step (28): inputting the multi-group multi-modal power generation amount of the third layer obtained in the step (27) as a load value to obtain four layers, and enabling the initial iteration step number k =0;
step (29): updating the consistency variable of each generator set on the fourth layer by combining the topological structure diagram of the agent and the formula (15) and the formula (16); a step (30): working out the active power output of each generator set on the fourth layer according to a formula (17);
step (31): correcting the active power output obtained in the step (30) according to a formula (18);
step (32): solving the deviation value of the active power of the fourth layer according to a formula (19);
step (33): judging whether the multi-mode active power deviation value obtained in the step (32) meets the precision requirement or not;
if it is notMaking the iteration step number k = k +1, returning to the step (29), and continuing to perform the iterative computation of the fourth layer;
if it is usedStoring the obtained active power output of the group of generator sets, then performing another iteration step number k = k +1, returning to the step (29), and continuing to perform iterative calculation of the fourth layer;
if it is notThe multi-mode active deviation value of the fourth layer meets the requirement, iteration is finished, and multiple groups of power generation are obtainedThe active power output of the unit;
step (34): performing multi-mode judgment and decision congestion distance screening on the active output of the multiple groups of generator sets on the fourth layer obtained in the step (33) to obtain a multi-mode active output decision scheme meeting the requirements of a decision maker;
step (35): judging whether T < T is met, if so, enabling T = T +1, and turning to the step (5); if not, executing the next step;
step (36): judging whether eta <1 is met, if yes, storing a multi-modal solution obtained when the value is eta, enabling eta = eta + R, and turning to the step (4); if not, executing the next step;
a step (37): and reversely solving all multi-modal solutions obtained by taking different eta values to obtain multi-modal multi-target values, wherein the multi-modal multi-target values obtained under all the eta values form a pareto frontier in a target space.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211257848.3A CN115545501A (en) | 2022-10-13 | 2022-10-13 | Multi-mode multi-target layered distributed comprehensive energy system economic dispatching method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211257848.3A CN115545501A (en) | 2022-10-13 | 2022-10-13 | Multi-mode multi-target layered distributed comprehensive energy system economic dispatching method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115545501A true CN115545501A (en) | 2022-12-30 |
Family
ID=84733649
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211257848.3A Pending CN115545501A (en) | 2022-10-13 | 2022-10-13 | Multi-mode multi-target layered distributed comprehensive energy system economic dispatching method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115545501A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115800276A (en) * | 2023-02-09 | 2023-03-14 | 四川大学 | Power system emergency scheduling method considering unit climbing |
-
2022
- 2022-10-13 CN CN202211257848.3A patent/CN115545501A/en active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115800276A (en) * | 2023-02-09 | 2023-03-14 | 四川大学 | Power system emergency scheduling method considering unit climbing |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106682810B (en) | Long-term operation method of cross-basin cascade hydropower station group under dynamic production of giant hydropower station | |
CN109492861B (en) | Method for decomposing medium-term electricity quantity trading plan of cascade hydropower station group | |
CN109740808B (en) | Wind-solar-water complementary power generation plan calculation method and system | |
CN109936170A (en) | Consider the honourable extreme misery complementation coordination optimization dispatching method of power supply flexibility nargin | |
CN107947244B (en) | New energy power station grid-connected power advanced control method | |
CN109347151A (en) | A kind of new energy participates in the sending end electric network source structural optimization method of peak regulation | |
CN111967666B (en) | Comprehensive cooperative scheduling system and scheduling method for river and lake water system | |
CN103106544B (en) | A kind of photovoltaic generation prognoses system based on T-S Fuzzy neutral net | |
CN107563538B (en) | Key water level control-oriented multi-target reservoir group scheduling optimization method under large power grid platform | |
Li et al. | Hierarchical multi-reservoir optimization modeling for real-world complexity with application to the Three Gorges system | |
CN111612292A (en) | Cascade hydropower station dispatching control system and method based on key water level control | |
CN114707767B (en) | New energy power system low-valley period adjustable peak power prediction method | |
CN115545501A (en) | Multi-mode multi-target layered distributed comprehensive energy system economic dispatching method | |
CN108092321A (en) | It is a kind of to consider probabilistic active power distribution network active reactive control method for coordinating | |
CN108649556A (en) | A kind of electric network emergency scene random optimization dispatching method | |
CN116526469A (en) | Long-term random dynamic scheduling method for water-wind-solar complementary system | |
CN115051388A (en) | Distribution robustness-based 'source-network-load-storage' two-stage scheduling optimization method | |
CN113255982A (en) | Medium-long term optimized scheduling method for wind-light-water complementary system | |
CN117526446A (en) | Wind-solar capacity double-layer optimization configuration method for cascade water-wind-solar multi-energy complementary power generation system | |
CN115545443A (en) | Multi-mode single-target layered distributed comprehensive energy system economic dispatching method | |
CN112686472B (en) | Power prediction method for distributed photovoltaic equivalent power station | |
CN113947291A (en) | Multi-mode distributed multi-target layered intelligent comprehensive energy system scheduling method | |
CN111612271B (en) | Clean energy delivery planning and operation optimization method considering clean energy utilization rate | |
CN111079066B (en) | Analysis method for two-objective competition relationship of power generation ecology of reservoir group | |
CN115239117A (en) | Classification and feature selection method based on multi-objective optimization |
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 |