CN113256045B - Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty - Google Patents

Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty Download PDF

Info

Publication number
CN113256045B
CN113256045B CN202010773898.1A CN202010773898A CN113256045B CN 113256045 B CN113256045 B CN 113256045B CN 202010773898 A CN202010773898 A CN 202010773898A CN 113256045 B CN113256045 B CN 113256045B
Authority
CN
China
Prior art keywords
power
gas
energy
heat
maximum
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
CN202010773898.1A
Other languages
Chinese (zh)
Other versions
CN113256045A (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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN202010773898.1A priority Critical patent/CN113256045B/en
Publication of CN113256045A publication Critical patent/CN113256045A/en
Application granted granted Critical
Publication of CN113256045B publication Critical patent/CN113256045B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • G06Q10/06375Prediction of business process outcome or impact based on a proposed change
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • 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
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty. Because the CCG method has better calculation efficiency in processing the robust optimization model, the method adopts the CCG method to solve. And finally, carrying out example simulation by using a commercial solver Gurobi to obtain a two-stage robust optimization scheduling strategy considering the combined thermoelectric demand response, verifying that the method can better treat the uncertainty in the system, and simultaneously can reduce the operation cost of the park and promote the consumption of new energy.

Description

Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty
Technical Field
The invention belongs to the technical field of optimization operation of a comprehensive energy system, and particularly relates to a park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty.
Background
The variance of the actual output and the predicted output of the wind power is far larger than the traditional load variance, and the traditional power system operation mode is not enough to ensure the reliability of the system. Due to physical limitations of an electric power system (such as climbing limitation of capacity of a conventional generator and a power transmission line), wind power reduction frequently occurs, so that the wind power utilization rate is low, and the enthusiasm of wind power investment is inhibited in the long run. When a large-scale renewable energy source is connected into a power system, the output of the renewable energy source is difficult to be accurately predicted due to strong randomness and intermittence of the renewable energy source, and in planning and scheduling of the power system, the error of the output predicted value of the renewable energy source causes non-negligible influence, in other words, the random fluctuation of renewable energy source power generation reduces the supply flexibility of the energy source system. Therefore, the traditional deterministic optimal scheduling method cannot meet the requirements, and the uncertainty of renewable energy sources needs to be considered on the basis of the original deterministic optimal scheduling method.
Stochastic optimization and robust optimization are typical methods for solving the problem of optimization with uncertainty. Compared with random optimization, the robust optimization has the advantages of simple data acquisition, high solving speed, suitability for solving large-scale uncertainty problems and the like, and is widely applied to processing uncertainty problems. At present, many scholars research the application of robust optimization in an integrated energy system, but the existing literature does not fully consider the influence of the energy coupling relation of various devices in the system and the response of a load side demand side on the uncertainty of renewable energy in the process of utilizing the robust optimization to process the wind-solar uncertainty of the integrated energy system.
The comprehensive energy system model can adaptively adjust the output of the unit to adapt to the change of the generated energy of the renewable energy source through the energy conversion relation among the devices, and the running safety of the system is ensured. The combined thermoelectric demand response can fully utilize the coupling relation among various devices, further improve the capability of the system for dealing with the uncertainty of renewable energy sources, reduce the abandoned wind and abandoned light and increase the permeability of the renewable energy sources. Therefore, on the basis of the existing research, a two-stage adjustable robust optimization model of the park comprehensive energy system based on the combined heat and power demand response is further considered, and the two-stage adjustable robust optimization model has important significance for the optimization operation of the comprehensive energy system.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty, a two-stage adjustable robust optimization model considering combined thermoelectric demand response is established to process uncertainty problems caused by wind and light output in a park, a dual-layer max-min problem in the model is converted into a single-layer max problem through a dual theory, and a CCG (column and constraint generation) method is adopted to solve, so that the dispatching result is more practical while the solving speed is improved, and the consumption of new energy is further promoted.
In order to solve the technical problems, the invention adopts the technical scheme that: in order to consider the uncertainty of renewable energy sources and aim at the problem of day-ahead economic dispatching of the park integrated energy system, the park integrated energy system day-ahead economic dispatching method considering the uncertainty of wind and light is provided, and a two-stage adjustable robust model considering the combined heat and power demand response is established. The day-ahead economic dispatching of the park comprehensive energy system aims at a basic scene with wind power and photovoltaic output values as predicted values, and when uncertainty occurs in operation, the park comprehensive energy system can adaptively and safely redistribute a generator set, a heat supply unit, P2G equipment, energy storage equipment and energy exchange between a park and a superior network. The optimized scheduling decision completely accords with the idea of a two-stage adjustable robust model. In other words, the start-stop state of the gas turbine and the cogeneration unit is determined in the first stage under the basic scene, and the worst unsafe scene of the system is searched when uncertainty occurs in the second stage. The two-stage adjustable robust scheduling model provided by the invention assumes that the start-stop state of the unit is a first-stage variable, namely when the uncertain variable of the system changes in the fluctuation interval of the uncertain variable, the start-stop state of the unit is kept unchanged, because the physical characteristics of most generator sets limit that the start-stop state of the unit cannot be changed rapidly under the uncertain condition.
A park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty comprises the following steps:
(1) Determining the specific composition of the multi-energy park, including the introduced new energy form and the specific equipment composition;
(2) Establishing models of various energy conversion equipment in the park;
(3) Establishing a demand response model;
(4) On the premise of meeting system safety constraints, establishing a two-stage adjustable robust model considering combined heat and power demand response by taking the minimum running cost of a basic scene as a target function;
(5) Obtaining an abstract expression of a two-stage adjustable robust optimization model of the park comprehensive energy system day-ahead economic dispatching;
(6) Establishing a maximum and minimum subproblem of worst scene identification;
(7) Solving a two-stage adjustable robust model considering combined heat and power demand response by using a CCG method;
(8) Inputting energy access, new energy output data, equipment parameters and operating parameters of the park integrated energy system, and solving a two-stage robust optimization model of the park integrated energy system day-ahead economic dispatching considering wind-solar uncertainty by adopting a commercial solver Gurobi to obtain a dispatching strategy of the park integrated energy system day-ahead economic dispatching.
Further, the specific composition of the park integrated energy system in the step (1) is as follows:
(1) The new energy form of the park integrated energy system is as follows: wind power and photovoltaic power generation;
(2) The energy conversion equipment introduced into the comprehensive energy system of the park comprises: electricity-to-gas equipment, an electric boiler, a gas turbine, a cogeneration unit and gas/heat storage equipment.
The energy conversion equipment models in the park are as follows;
(2.1) model of electric gas-converting apparatus
Figure SMS_1
Figure SMS_2
In the formula: t is scheduling time; m is an index of the electric-to-gas equipment;
Figure SMS_3
respectively gas making power, consumed power and electricity-to-gas efficiency, L, of an electricity-to-gas (P2G) facility HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure SMS_4
The minimum and maximum pneumatic power of the mth station P2G.
(2.2) electric boiler model
Figure SMS_5
Figure SMS_6
In the formula: t is scheduling time; n is an electric boiler index;
Figure SMS_7
and &>
Figure SMS_8
Respectively the power consumption and the heat production of the nth electric boiler in the time period t;
Figure SMS_9
Is the electric heat conversion efficiency of the nth electric boiler>
Figure SMS_10
Dividing into the minimum and maximum heating power of the nth electric boiler;
Figure SMS_11
And (4) starting and stopping states of the nth electric boiler in a time period t (1 represents starting, and 0 represents stopping). />
(2.3) gas turbine model
Figure SMS_12
Figure SMS_13
Figure SMS_14
Figure SMS_15
Figure SMS_16
Figure SMS_17
Figure SMS_18
In the formula: t is scheduling time; q is a gas turbine index;
Figure SMS_21
and &>
Figure SMS_22
Respectively representing the power generation power and the gas consumption power of the gas turbine; f (-) represents the gas turbine energy consumption curve;
Figure SMS_25
And &>
Figure SMS_20
Respectively representing the consumption of natural gas required by the startup and shutdown of the gas turbine; a is a q 、b q And c q Represents the gas coefficient of F (-); l is HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure SMS_24
And dividing the power into the minimum and maximum power generation power of the q gas turbines.
Figure SMS_27
Representing the power generation power of the qth gas turbine in the t period;
Figure SMS_28
The start-stop state (1 represents start-up and 0 represents stop) and the judgment is carried out on the qth gas turbine in the t period>
Figure SMS_19
Starting and stopping states of a qth gas turbine in a time period t;
Figure SMS_23
The upward climbing rate and the downward climbing rate of the qth gas turbine,
Figure SMS_26
for the continuous starting and stopping time of the qth gas turbine within the time period t-1, the pressure is adjusted>
Figure SMS_29
The minimum startup and shutdown time of the qth gas turbine in the time period t.
(2.4) Combined Heat and Power Unit model
Figure SMS_30
Figure SMS_31
Figure SMS_32
Figure SMS_33
Figure SMS_34
Figure SMS_35
Figure SMS_36
In the formula: t is scheduling time; p is an index of the cogeneration unit;
Figure SMS_39
and &>
Figure SMS_43
Respectively representing heat production power and gas consumption power of a combined heat and power generation unit (CHP);
Figure SMS_47
The generating power and the generating efficiency of the micro-combustion engine in the t period are based on the pressure value>
Figure SMS_38
For the heat dissipation loss rate, the device can be used for>
Figure SMS_44
And &>
Figure SMS_46
The heating coefficient and the flue gas recovery rate of the bromine refrigerator are respectively; l is HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure SMS_49
Dividing into the p (th) CHP unit minimum and maximum generating power;
Figure SMS_37
Representing the power generation power of the p-th cogeneration unit in the t period;
Figure SMS_42
The starting and stopping states of the pth CHP unit in the t period (1 represents starting, 0 represents stopping), and/or>
Figure SMS_45
The starting and stopping state of a P-th cogeneration unit in a time period t;
Figure SMS_48
The climbing rate and the descending rate of the pth CHP unit;
Figure SMS_40
For the continuous startup and shutdown time of the pth CHP unit in the t-1 time period,
Figure SMS_41
for the minimum startup of the p-th CHP unit in the time period tAnd downtime.
(2.5) energy storage device model
ESS(t)=(1-σ S )·ESS(t-1)+ES in (t)·η in -ES out (t)/η out
Figure SMS_50
Figure SMS_51
Figure SMS_52
Figure SMS_53
Figure SMS_54
Figure SMS_55
Figure SMS_56
Figure SMS_57
In the formula: t is scheduling time; ESS (t), ES in (t)、ES out (t) stored energy, stored energy power and released energy power of the energy storage device are respectively in a time period t; ESS (t-1) represents the stored energy of the energy storage device during the t-1 period; sigma S Is the self-consumption rate of the energy storage system; eta in 、η out Respectively storing energy and releasing energy efficiency for the energy storage equipment;
Figure SMS_58
gas storage, gas discharge power, GS, for a period of t in ,max 、GS out,max The maximum gas storage and gas release power of the gas storage device are respectively; c t GS The air storage quantity of the air storage device is greater or less for a time period t>
Figure SMS_59
The gas storage capacity of the gas storage equipment is t-1 time; c GS,min 、C GS,max Respectively the minimum and maximum gas storage capacity, eta, of the gas storage device CGS 、η GS,in 、η GS,out The self-consumption rate, the gas storage efficiency and the gas release efficiency of the gas storage equipment are improved.
Figure SMS_60
For heat-storage, heat-release power for period t, HS in,max 、HS out,max The maximum heat storage power and the maximum heat release power of the heat storage equipment are respectively;
Figure SMS_61
A reserve of heat storage apparatus for a time period t>
Figure SMS_62
The gas storage capacity of the heat storage equipment is t-1 time period; c HS,min 、C HS,max Minimum and maximum heat storage capacity, eta, of the heat storage apparatus CHS 、η HS,in 、η HS,out The self consumption rate, the heat storage efficiency and the heat release efficiency of the heat storage equipment are shown.
The demand response model in the step (3) is specifically as follows:
Figure SMS_63
Figure SMS_64
Figure SMS_65
Figure SMS_66
Figure SMS_67
Figure SMS_68
Figure SMS_69
in the formula: t is scheduling time;
Figure SMS_70
P t DR 、P t DR,inte 、P t DR,shif 、P t LD,max predicted value of electric load, demand side response electric load, demand side transfer electric load, maximum electric load allowed by system, P inte,max The maximum interruptible electrical load within the scheduled time.
Figure SMS_71
The maximum interruptible and extractable electrical load proportion for the t period.
Figure SMS_72
Positive means turning out of translatable load, and vice versa.
Figure SMS_73
Figure SMS_74
Figure SMS_75
In the formula: t is scheduling time;
Figure SMS_76
and &>
Figure SMS_77
Respectively obtaining a predicted value of the thermal load, a demand side response thermal load proportion, a maximum thermal load allowed by a system and a thermal load after considering demand response in a time period t; h DR,max The maximum interruptible thermal load within the scheduling time; n is a radical of hydrogen t The time period is scheduled for the whole time.
The two-stage adjustable robust model considering the combined heat and power demand response in the step (4) is specifically as follows:
(4.1) objective function
Figure SMS_78
Figure SMS_79
Figure SMS_80
Figure SMS_81
Figure SMS_82
Figure SMS_83
Figure SMS_84
Figure SMS_85
Figure SMS_86
Figure SMS_87
In the formula: t is scheduling time;
Figure SMS_90
respectively the electricity and gas purchase cost;
Figure SMS_93
The starting-up and shutdown costs of the pth CHP are respectively;
Figure SMS_96
Penalty fee for abandoning wind>
Figure SMS_89
Punishment of cost for light abandonment;
Figure SMS_94
To sell the electricity, C E,cc A compensation cost to interrupt the electrical load for demand side response;
Figure SMS_98
For purchasing electric power>
Figure SMS_100
In order to obtain the gas purchasing power,
Figure SMS_88
wind abandoning power and light abandoning power of the ith fan and the jth group of photovoltaic cells in the t period are respectively set;
Figure SMS_92
For selling electricity power, is selected>
Figure SMS_99
To needAnd the side is asked to interrupt the electric load. Δ t is the scheduling time interval, N t The number of time segments is scheduled.
Figure SMS_102
Respectively the unit prices of purchasing electricity, purchasing gas and selling electricity>
Figure SMS_91
Penalty prices in units of wind abandon and light abandon, respectively>
Figure SMS_95
Is a compensation price per unit demand side response to an interrupted electrical load.
Figure SMS_97
The start-stop state (1 represents the start-up and 0 represents the shutdown) of the pth CHP in the t period, and the judgment is performed>
Figure SMS_101
The cost of starting up and stopping the p-th CHP is respectively. N is a radical of hydrogen WT 、N PV 、N CHP The number of the fan, the photovoltaic cell, the electric boiler and the CHP are respectively.
(4.2) constraint Condition
Figure SMS_103
Figure SMS_104
Figure SMS_105
Figure SMS_106
Figure SMS_107
Figure SMS_108
Figure SMS_109
Figure SMS_110
Figure SMS_111
v 1t ,v 2t ,v 3t ,v 4t ≥0/>
Figure SMS_112
Figure SMS_113
Figure SMS_114
Figure SMS_115
Figure SMS_116
Figure SMS_117
Figure SMS_118
Figure SMS_119
Figure SMS_120
Figure SMS_121
Figure SMS_122
In the formula: t is scheduling time; (.) u The wind and light output is a corresponding adjusted variable under real-time change; v. of 1,t And v 2,t A slack variable is a power balance constraint; v. of 3,t And v 4,t Constraint relaxation variables for thermal balance; omega WT 、Ω PV Respectively representing uncertain collections of wind power output and photovoltaic output;
Figure SMS_123
deviation of wind power, photovoltaic output and predicted value respectively>
Figure SMS_124
The ratio of the wind power and photovoltaic output deviation value to the predicted value is obtained;
Figure SMS_125
Is a variable from 0 to 1 in the uncertain convergence; delta i 、Δ j Respectively obtaining wind power and photovoltaic output uncertainty precalculated values;
Figure SMS_126
Correcting the climbing and descending of the pth CHP unit;
Figure SMS_127
Correcting the climbing and descending of the qth gas turbine unit; and lambda (-) is a dual variable corresponding to the constraint condition.
The abstract expression of the two-stage adjustable robust optimization model of the park comprehensive energy system economic dispatch in the past is as follows:
Figure SMS_128
in the formula: x represents the starting and stopping states of units related to the CHP and the gas turbine, y and z represent the basic scene and the scheduling output of other units of the system adjusted according to the wind-light output transformation, u is an uncertain variable related to the uncertainty of the wind power output and the photovoltaic output, and c b 、c g A, B, B, F, h, C, D, E, F, G can be derived from the objective function and constraint conditions in 5.
The maximum and minimum subproblems identified in the worst scene in the step (6) are specifically as follows:
Figure SMS_129
Figure SMS_130
Figure SMS_131
the two-stage adjustable robust model set for solving and considering combined heat and power demand response by using the CCG method in the step (7) is as follows:
(7.1) Main problem
Figure SMS_132
Ax+By≤b
(7.2) sub-problem
Figure SMS_133
Figure SMS_134
(7.3) CCG solving step
Step 1: let iteration counter s =0 set the maximum value epsilon of security violation allowed by the system RO
Step 2: solving the main problem, if the main problem is solved, obtaining a system unit starting and stopping state x and a unit output arrangement y, and performing the step 3; otherwise, stopping iteration and outputting no solution.
And step 3: and (3) solving the maximum and minimum subproblems according to the x and y obtained by the solution in the step (2), and finding out the wind power and photovoltaic output under the worst scene which causes the maximum possibility of violating the safety specified value.
And 4, step 4: if the maximum possible violation safety-specifying value solved in step 3 is less than epsilon RO Then x and y are the final optimization solution and the iteration is stopped; otherwise, let s = s +1, according to the wind power and photovoltaic output value under the worst scene solved in step 3
Figure SMS_135
The CCG constraint shown below is added to the main question, returning to step 2.
f T v s ≤ε RO
Figure SMS_136
And (8) the data of the park integrated energy system also comprises the specific composition of the park integrated energy system, energy price, equipment parameters and values of each energy conversion equipment, demand response proportion, new energy output fluctuation conditions in a basic scene and a worst scene, maximum violation safety specified value and a predicted value of new energy output and load.
Compared with the prior art, the invention has the beneficial effects that:
(1) The optimization problem containing uncertainty is solved by adopting robust optimization, and a two-stage adjustable robust optimization model is established, so that a result close to reality can be obtained by only acquiring simple data in the scheduling process of the park comprehensive energy system, and the solving speed is high. The improved two-stage adjustable robust optimization model considers the wind and light uncertainty of the comprehensive energy system and simultaneously fully considers the influence of the energy coupling relation of various devices in the system and the response of the load side on the uncertainty of renewable energy. The day-ahead economic dispatching of the park comprehensive energy system aims at a basic scene with wind power and photovoltaic output values as predicted values, and when uncertainty occurs in operation, the park comprehensive energy system can adaptively and safely redistribute a generator set, a heat supply unit, P2G equipment, energy storage equipment and energy exchange between a park and a superior network. The two-stage robust optimization ensures that the system can meet safety constraints under any condition, minimizes the operation cost of a basic scene, and meets the requirements of system safety and economy.
(2) After the combined heat and power demand response is considered, the functions of gas energy storage, heat energy storage and electric heat load demand response are fully exerted, the coupling relation among various devices can be fully utilized by the combined heat and power demand response, the uncertainty capability of the system for coping with renewable energy sources is further improved, the wind and light curtailment is reduced, and the permeability of the renewable energy sources is increased. The comprehensive energy system model can adaptively adjust the output of the unit to adapt to the change of the generated energy of the renewable energy source through the energy conversion relation among the devices, thereby ensuring the safety of the operation of the system, promoting the consumption of the renewable energy source, improving the robustness of the system and improving the economic benefit of the park.
Drawings
FIG. 1 is a flow chart of the steps of the method of the present invention;
FIG. 2 is a diagram of the specific components of the park integrated energy system;
FIG. 3 is a CCG method solving flow chart;
FIG. 4 is a forecast of wind turbines, photovoltaic output, electrical load and thermal load for a typical winter day for a park multi-energy integrated system.
Detailed Description
In order to explain the technical solutions disclosed in the present invention in detail, the present invention is further explained with reference to the accompanying drawings and specific embodiments.
The invention discloses a park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty. The specific implementation step flow is shown in fig. 1, and the technical scheme of the invention comprises the following steps:
step 1: the specific composition of the park integrated energy system is determined, including the new energy form introduced and the specific equipment composition.
(1.1) the new energy form of the park integrated energy system is as follows: wind power and photovoltaic power generation;
(1.2) energy conversion equipment introduced into a park comprehensive energy system comprises: electricity changes gas equipment, electric boiler, gas turbine, combined heat and power generation unit, gas storage/heat-retaining equipment.
Step 2: and establishing various energy conversion equipment models in the park.
(2.1) model of electric gas-converting apparatus
The P2G technology can realize conversion from electric energy to natural gas, the natural gas is transmitted in a large-capacity, long-time and long-distance mode through a natural gas pipeline, powerful technical support is provided for renewable energy consumption, large-range and long-distance space-time transfer of wind power can be achieved, meanwhile, the P2G response is rapid, and the application prospect is strong. The relationship between the P2G equipment gas making power and the power consumption and the limitation of the gas making power are as follows:
Figure SMS_137
Figure SMS_138
in the formula: t is scheduling time; m is an index of the electric-to-gas equipment;
Figure SMS_139
respectively the gas making power, the electricity consumption power and the electricity conversion gas efficiency of the electricity conversion gas (P2G) equipment, L HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure SMS_140
The minimum and maximum pneumatic power of the mth station P2G.
(2.2) electric boiler model
The introduction of the electric boiler can break the electric-thermal coupling hard constraint of the CHP unit and change the traditional 'electricity by heat' scheduling mode. The electric boiler can coordinate the peak valley of the electric heating load, and the relation between the heating capacity and the power consumption and the heating capacity constraint are as follows:
Figure SMS_141
Figure SMS_142
in the formula: t is scheduling time; n is an electric boiler index;
Figure SMS_143
and &>
Figure SMS_144
Respectively the power consumption and the heat production of the nth electric boiler in the time period t;
Figure SMS_145
Is the electric heat conversion efficiency of the nth electric boiler>
Figure SMS_146
Dividing into the minimum and maximum heating power of the nth electric boiler;
Figure SMS_147
The start-stop state of the nth electric boiler in the time period t is shown (1 represents start-up, and 0 represents stop).
(2.3) gas turbine model
The gas turbine converts chemical energy in natural gas into electric energy, and the relation between the consumed natural gas power and the generated power, the generated power limit, the climbing constraint and the minimum on-off time constraint are as follows:
Figure SMS_148
Figure SMS_149
Figure SMS_150
Figure SMS_151
Figure SMS_152
Figure SMS_153
Figure SMS_154
in the formula: t is scheduling time; q is a gas turbine index;
Figure SMS_156
and &>
Figure SMS_158
Respectively representing the power generation power and the gas consumption power of the gas turbine; f (-) represents the gas turbine energy consumption curve;
Figure SMS_165
And &>
Figure SMS_155
Respectively representing the consumption of natural gas required by the startup and shutdown of the gas turbine; a is q 、b q And c q Represents the gas coefficient of F (-); l is HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure SMS_159
Dividing into the q gas turbine minimum and maximum power generation workAnd (4) rate.
Figure SMS_161
Representing the generated power of the qth gas turbine in the t period;
Figure SMS_163
The start-stop state (1 represents start-up and 0 represents stop) and the judgment is carried out on the qth gas turbine in the t period>
Figure SMS_157
Starting and stopping states of a qth gas turbine in a time period t;
Figure SMS_160
The upward climbing rate and the downward climbing rate of the qth gas turbine,
Figure SMS_162
for the continuous starting and stopping time of the qth gas turbine within the time period t-1, the pressure is adjusted>
Figure SMS_164
The minimum startup and shutdown time of the qth gas turbine in the time period t.
(2.4) Combined Heat and Power Unit model
Neglecting the influence of external environment change on the power generation and fuel combustion efficiency, the mathematical relationship of the thermoelectric relationship, the gas-electricity relationship, the relationship between the power generation power, the power generation power limit, the climbing constraint and the minimum on-off time constraint are as follows:
Figure SMS_166
Figure SMS_167
Figure SMS_168
Figure SMS_169
Figure SMS_170
Figure SMS_171
Figure SMS_172
in the formula: t is scheduling time; p is an index of the cogeneration unit;
Figure SMS_175
and &>
Figure SMS_179
Respectively representing heat production power and gas consumption power of a combined heat and power generation unit (CHP);
Figure SMS_181
The generating power and the generating efficiency of the micro-combustion engine in the t period are based on the pressure value>
Figure SMS_176
For heat dissipation loss rate>
Figure SMS_177
And &>
Figure SMS_183
The heating coefficient and the flue gas recovery rate of the bromine refrigerator are respectively; l is HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure SMS_185
Dividing into the p-th CHP unit minimum and maximum generating power;
Figure SMS_173
Representing the power generation power of the p-th cogeneration unit in the t period;
Figure SMS_180
Starting and stopping states (1 represents starting, 0 represents stopping) and/or on/off for the pth CHP unit in a t period>
Figure SMS_182
The starting and stopping state of a P-th cogeneration unit in a time period t;
Figure SMS_184
The climbing rate and the descending rate of the pth CHP unit;
Figure SMS_174
For the continuous startup and shutdown time of the pth CHP unit in the t-1 time period,
Figure SMS_178
the minimum startup and shutdown time of the p-th CHP unit in the time period t is obtained.
(2.5) energy storage device model
The relationship among the energy storage device model, the gas storage/release of the heat storage device, the constraint of the stored/released power and the gas storage at the time t, the constraint of the gas storage and the heat storage capacity of the heat storage device, the gas storage at the time t, the gas storage and the heat storage capacity of the heat storage device, the gas storage at the time t-1, the stored heat quantity, the stored/released gas at the time t and the stored/released power is as follows:
ESS(t)=(1-σ S )·ESS(t-1)+ES in (t)·η in -ES out (t)/η out
Figure SMS_186
Figure SMS_187
Figure SMS_188
Figure SMS_189
Figure SMS_190
Figure SMS_191
Figure SMS_192
Figure SMS_193
in the formula: t is scheduling time; ESS (t), ES in (t)、ES out (t) storing energy, stored energy power and released energy power of the energy storage device at a time period t, respectively; ESS (t-1) represents the stored energy of the energy storage device during the t-1 period; sigma S Is the self-consumption rate of the energy storage system; eta in 、η out Respectively storing energy and releasing energy efficiency for the energy storage equipment;
Figure SMS_194
gas storage, gas discharge power, GS, for a period of time t in ,max 、GS out,max The maximum gas storage and gas release power of the gas storage device are respectively;
Figure SMS_195
The air storage quantity of the air storage device is greater or less for a time period t>
Figure SMS_196
The gas storage capacity of the gas storage equipment is t-1 time; c GS,min 、C GS,max Respectively the minimum and maximum gas storage capacity, eta, of the gas storage device CGS 、η GS,in 、η GS,out The self-consumption rate, the gas storage efficiency and the gas release efficiency of the gas storage equipment are improved.
Figure SMS_197
For heat-storage, heat-release power for period t, HS in,max 、HS out,max The maximum heat storage power and the maximum heat release power of the heat storage equipment are respectively;
Figure SMS_198
A reserve of heat storage apparatus for a time period t>
Figure SMS_199
The gas storage capacity of the heat storage equipment is t-1 time period; c HS,min 、C HS,max Minimum and maximum heat storage capacity, eta, of the heat storage apparatus CHS 、η HS,in 、η HS,out The self consumption rate, the heat storage efficiency and the heat release efficiency of the heat storage equipment are obtained.
And step 3: and establishing a demand response model.
Dividing the demand side electric load capable of responding into a translatable electric load and an interruptible electric load, considering the interrupting cost of the interruptible load, wherein the sum of the demand side electric load capable of responding and the current time interval electric load needs to be less than the maximum electric load allowed in the current time interval, and the relation between the predicted value of the electric load and the electric load after the demand response is as follows:
Figure SMS_200
Figure SMS_201
Figure SMS_202
Figure SMS_203
Figure SMS_204
Figure SMS_205
Figure SMS_206
in the formula: t is scheduling time;
Figure SMS_207
P t DR 、P t DR,inte 、P t DR,shif 、P t LD,max predicted value of electric load, demand side response electric load, demand side transfer electric load, maximum electric load allowed by system, P inte,max The maximum interruptible electrical load during the scheduled time.
Figure SMS_208
The maximum interruptible and extractable electrical load proportion for the t period. P t DR,shif Positive means turning out of translatable load, and vice versa.
The demand side may respond to thermal load constraints as follows:
Figure SMS_209
Figure SMS_210
Figure SMS_211
in the formula: t is scheduling time;
Figure SMS_212
and &>
Figure SMS_213
The predicted value of the thermal load, the response thermal load of the demand side, the response thermal load proportion of the demand side and the system are respectivelyMaximum allowable thermal load and electrical load after considering demand response; h DR,max The maximum interruptible thermal load within the scheduling time; n is a radical of t The entire schedule period.
And 4, step 4: on the premise of meeting the system safety constraint, a two-stage adjustable robust model considering the combined heat and power demand response is established by taking the minimum running cost of a basic scene as an objective function.
(1) Objective function
The two-stage robust optimization model provided by the invention aims to minimize the running cost of a basic scene on the premise of meeting the system safety constraint, and because no load loss is allowed in the basic scene, an objective function does not contain a load loss related variable, and the objective function and a related equation are constrained as follows:
Figure SMS_214
Figure SMS_215
Figure SMS_216
Figure SMS_217
Figure SMS_218
Figure SMS_219
Figure SMS_220
Figure SMS_221
Figure SMS_222
in the formula: t is scheduling time;
Figure SMS_225
respectively the electricity and gas purchase cost;
Figure SMS_228
The starting-up and shutdown costs of the pth CHP are respectively;
Figure SMS_230
Penalty fee for wind abandon, based on the sum of the wind and the wind>
Figure SMS_226
Punishment of cost for light abandonment;
Figure SMS_229
To earn electricity sales, C E,cc A compensation cost to interrupt the electrical load for demand side response;
Figure SMS_232
For purchasing power, is turned on or off>
Figure SMS_234
In order to purchase the gas power, the gas-purchasing power,
Figure SMS_223
wind abandoning power and light abandoning power of the ith fan and the jth group of photovoltaic cells in a time period t are respectively set; p t out To sell electric power, P t DR,inte The electrical load is interrupted for the demand side. Δ t is the scheduling time interval, N t Is the number of scheduling periods.
Figure SMS_227
Respectively the unit purchase electricity, gas and electricity prices phi t WT 、φ t PV Punishment price of wind abandoning and light abandoning respectivelyShelf, or>
Figure SMS_231
Is a compensation price per unit demand side response to an interrupted electrical load.
Figure SMS_233
The start-stop state (1 represents the start-up and 0 represents the shutdown) of the pth CHP in the t period, and the judgment is performed>
Figure SMS_224
The cost of starting up and stopping the p-th CHP is respectively. N is a radical of WT 、N PV 、N CHP The number of the fan, the photovoltaic cell, the electric boiler and the CHP are respectively.
(2) Constraint conditions
In order to ensure the safe operation of the system, the section adopts a double-layer max-min model to identify the worst scene causing the least safe operation of the system, namely the maximum security violation specified value (security operation), wherein the maximum security violation specified value of the system in the worst scene is required to be smaller than a preset value epsilon RO ,ε RO The values are set in relation to predetermined system security levels to ensure safe operation of the campus energy complex. The method comprises the following steps of causing violation of the maximum value of system safety regulation, uncertainty set of wind power output and photovoltaic output, system energy balance constraint, constraint that a slack variable is constantly greater than zero, unit output constraint, unit error correction climbing constraint under uncertainty of renewable energy output, energy storage constraint, power exchange constraint with a superior network, demand side response constraint and wind and light abandoning constraint in the worst scene under an uncertain condition as follows:
Figure SMS_235
Figure SMS_236
Figure SMS_237
Figure SMS_238
Figure SMS_239
Figure SMS_240
Figure SMS_241
Figure SMS_242
v 1t ,v 2t ,v 3t ,v 4t ≥0
Figure SMS_243
Figure SMS_244
Figure SMS_245
Figure SMS_246
Figure SMS_247
Figure SMS_248
Figure SMS_249
Figure SMS_250
Figure SMS_251
Figure SMS_252
Figure SMS_253
Figure SMS_254
Figure SMS_255
Figure SMS_256
P t u,out ≤P out,min ,P t u,out ≤P out,max :(λ 9_3,t9_4,t )
Figure SMS_257
Figure SMS_258
Figure SMS_259
Figure SMS_260
in the formula: t is scheduling time;
Figure SMS_263
respectively comprise abandoned wind, abandoned light, electric load after demand response, P2G consumed electric power, electric boiler consumed electric power, park electricity selling power, CHP unit output and gas turbine output according to different wind power outputs>
Figure SMS_265
And photovoltaic output->
Figure SMS_267
An adjusted real-time value; (. Cndot.) u The wind and light output is a corresponding adjusted variable under real-time change; v. of 1,t And v 2,t A slack variable is a power balance constraint; v. of 3,t And v 4,t Constraint relaxation variables for thermal balance; omega WT 、Ω PV Respectively representing uncertain collections of wind power output and photovoltaic output;
Figure SMS_262
Deviation between wind power output and photovoltaic output and a predicted value, respectively>
Figure SMS_264
The ratio of the wind power and photovoltaic output deviation value to the predicted value is obtained;
Figure SMS_266
Is a variable from 0 to 1 in the uncertain convergence; delta of i 、Δ j Respectively calculating uncertainty precalculated values of wind power output and photovoltaic output;
Figure SMS_268
Correcting the climbing and descending of the pth CHP unit;
Figure SMS_261
Correcting the climbing and descending of the qth gas turbine unit; and lambda (-) is a dual variable corresponding to the constraint condition.
And 5, obtaining an abstract expression of a two-stage adjustable robust optimization model of the park comprehensive energy system day-ahead economic dispatching.
The proposed two-stage robust optimization model relates to three systems of electricity, gas and heat, involves more constraints, contains uncertain parameters, and is a nonlinear mixed integer programming problem. For ease of discussion, the two-stage robust optimization scheduling model proposed in this section may take the form of an abstract robust optimization model as follows:
Figure SMS_269
s.t.Ax+By≤b
Figure SMS_270
in the formula: x represents the starting and stopping states of units related to the CHP and the gas turbine, y and z represent the basic scene and the scheduling output of other units of the system adjusted according to the wind-light output transformation, u is an uncertain variable related to the uncertainty of the wind power output and the photovoltaic output, and c b 、c g A, B, B, C, D, E, F, G can be derived from the objective function and constraint conditions in 5.
And 6, establishing the maximum and minimum subproblems of the worst scene identification.
And 4, the sub-problem is the problem of identifying the worst scene, the scene causing the maximum violation of the safety specified value of the system is found through the double-layer max-min problem shown in the step 4, namely the specific value of the uncertain quantity in the worst scene is determined, and the double-layer max-min problem is converted into a single-layer bilinear maximum optimization sub-problem shown below through a dual theory.
Figure SMS_271
Figure SMS_272
Figure SMS_273
And 7, solving a two-stage adjustable robust model considering the joint thermoelectric demand response by utilizing a CCG (column and constraint generation) method.
(7.1) Main problem
Figure SMS_274
Ax+By≤b
(7.2) sub-problem
Figure SMS_275
Figure SMS_276
(7.3) CCG solving step
Step (7.3.1): let the iteration counter s =0 set the maximum value epsilon allowed by the system for violating the safety regulations RO
Step (7.3.2): solving the main problem, if the main problem is solved, obtaining a system unit starting and stopping state x and a unit output arrangement y, and performing the step (7.3.3); otherwise, stopping iteration and outputting no solution.
Step (7.3.3): and (4) solving the maximum and minimum subproblems according to the x and y obtained by the solution in the step (7.3.2), and finding out the wind power and photovoltaic output under the worst scene which causes the maximum possibility of violating the safety specified value.
Step (7.3.4): if the maximum possible violation safety provision resolved in step (7.3.3) is less than ε RO Then x and y are the final optimization solution and the iteration is stopped; conversely, let s = s +1, according to the worst scene solved in step (7.3.3)Lower wind and photovoltaic output value
Figure SMS_277
The CCG constraint shown below is added to the main question, returning to step (7.3.2).
f T v s ≤ε RO
Figure SMS_278
And 8, inputting the data of the park comprehensive energy system, and further comprising the specific composition, energy price, equipment parameters and values of each energy conversion equipment, demand response proportion, new energy output fluctuation conditions under a basic scene and a worst scene, maximum violation safety specified values, predicted values of new energy output and load and the like of the park comprehensive energy system, and solving the two-stage adjustable robust optimization operation model of the park comprehensive energy system by adopting a commercial solver Gurobi to obtain a robust optimization scheduling result.
The effects of the present invention will be described in detail below with reference to specific examples.
(1) Introduction to the examples.
Fig. 2 shows the specific components of the park energy system. And (4) selecting a multi-energy complementary park system model containing wind, light, gas, storage and consideration of electricity-to-gas and electricity-to-heat technologies in a simulation mode. The system comprises a gas turbine, a wind driven generator, an electric boiler, P2G equipment, heat storage equipment, gas storage equipment, two CHPs and a group of photovoltaic cells. Heating coefficient of bromine refrigerator
Figure SMS_279
And flue gas recovery ratio>
Figure SMS_280
Respectively 0.9 and 1.2, and the one-time startup and shutdown costs of the gas turbine, the CHP and the electric boiler are respectively as follows: 3.5, 1.94, 2.74 yuan. Assuming that the initial state of the CHP and the gas turbine is the off-stream state, the initial gas storage capacity of the gas storage equipment is 10m 3 The initial heat storage capacity of the heat storage equipment is 100 kW.h, and the self-consumption rate of the gas storage and heat storage equipment is 0.01. Multiple energy sources for parkThe predicted values of the fan, photovoltaic output, electrical load and thermal load of the integrated system on a typical winter day are shown in fig. 4.
(2) Description of embodiment scenarios.
In order to verify the effectiveness of the proposed two-stage robust adjustable optimization model considering the response of the combined heat and power demand, 7 scheduling operation modes shown in table 1 are set.
TABLE 1 7 scheduling operation modes
Figure SMS_281
Figure SMS_282
The wind power and photovoltaic predicted output is adjusted to be 2 times of the original output, and meanwhile uncertainty of wind and photovoltaic output is considered, and a new operation scheme 1-7 is obtained. The threshold value of the safety regulation violation is set to be 0, namely the system does not allow to lose load and overload under any scene, and the uncertain budget of the wind power output and the photovoltaic output is assumed to be 24.
And increasing the electrical load to 2.4 times of the original electrical load, and simultaneously considering the uncertainty of wind-solar output to obtain a new operation scheme 8-14, wherein the violation safety specified value is set to be 0.
(3) Examples analysis of results.
Table 2 gives the operating costs and violation of safety regulations for operating scenarios 1-7 at different uncertainty ratios, from which can be derived: along with the increase of the uncertain proportion of the wind power output and the photovoltaic output, the operation cost of the park is increased, and the maximum value of violating the safety regulations possibly appearing in the park is larger. However, with the continuous addition of energy conversion equipment such as heat storage equipment, P2G and gas storage equipment thereof, the park can continuously adjust the energy storage state in real time to cope with the real-time change of wind power photovoltaic output, the running risk of the system is reduced, the uncertainty of the system can be flexibly coped with by considering demand response, the combined thermoelectric demand response is considered, and the running cost under the basic scene of the park is further reduced.
TABLE 2 running costs and violations of safety regulations for running schemes 1-7 at different uncertainty ratios
Figure SMS_283
Figure SMS_284
Table 3 gives the number of iterations and violation of the safety specified values for the operating scenario 1 under different wind power and photovoltaic uncertainty budgets, from which can be derived: the larger the uncertainty budget value of the renewable energy source is, the fewer the number of iterations required for obtaining an optimization scheme meeting the system requirements is, and the larger the violation of the safety regulation value which may occur in the system is. The robustness of the system can be adjusted by adjusting the uncertainty budget value.
TABLE 3 iteration number and violation of safety regulation value of operation scheme 1 under different wind power and photovoltaic uncertainty budgets
Figure SMS_285
Figure SMS_286
Table 4 gives the operating costs and violation of safety regulations for operating scenarios 8-14 at different uncertainty ratios, from which one can derive: when the system load demand is obviously greater than the renewable energy output, the renewable energy output can be completely absorbed by the system, the gas-to-electricity conversion is superior to the electricity-to-gas conversion, the P2G equipment does not work, and the robustness of the system cannot be improved due to the introduction of the energy storage equipment. When the uncertain ratio of the output of the renewable energy is 0.2, compared with the schemes 8-10, the more comprehensive the equipment contained in the system is, the lower the running cost of the system in the basic scene is, and compared with the schemes 10-14, the introduction of the response of the demand side can further reduce the running cost of the system. Given the uncertain proportion of operating costs and violations of safety regulations for the different renewable energy sources of the comparison schemes 11-14, the combined heat and power demand response can reduce power load losses by reducing heat loads, greatly enhancing the robustness, flexibility and economy of the campus integrated energy system.
TABLE 4 running cost and violation of safety regulations for running schemes 8-14 at different uncertainty ratios
Figure SMS_287
Figure SMS_288
The above description is only an embodiment of the present invention, but not intended to limit the scope of the present invention, and all equivalent changes or substitutions made by using the contents of the present specification and the drawings, which are directly or indirectly applied to other related arts, should be included within the scope of the present invention.

Claims (1)

1. A park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty is characterized by comprising the following steps:
step 1: determining the specific composition of the comprehensive energy system of the park, including the introduced new energy form and the specific equipment composition;
step 2: establishing models of various energy conversion equipment in the park;
and 3, step 3: establishing a demand response model;
and 4, step 4: on the premise of meeting system safety constraints, establishing a two-stage adjustable robust model considering combined heat and power demand response by taking the minimum running cost of a basic scene as a target function;
step 5, obtaining an abstract expression of a two-stage adjustable robust optimization model of the park comprehensive energy system day-ahead economic dispatching;
step 6, establishing the maximum and minimum subproblems of the worst scene identification;
step 7, solving a two-stage adjustable robust model considering combined heat and power demand response by using a CCG method;
and step 8: inputting energy access, new energy output data, equipment parameters and operating parameters of the park integrated energy system, and solving a two-stage robust optimization model of the park integrated energy system day-ahead economic dispatching considering wind-solar uncertainty by adopting a commercial solver to obtain a dispatching strategy of the park integrated energy system day-ahead economic dispatching model;
the park comprehensive energy system in the step 1 specifically comprises the following components:
(1) The new energy form of the integrated energy system accessed to the park is as follows: wind power and photovoltaic power generation;
(2) The energy conversion equipment for introducing the park comprehensive energy system comprises: electricity-to-gas equipment, an electric boiler, a gas turbine, a cogeneration unit, and gas/heat storage equipment;
step 2, each energy conversion equipment model is as follows;
(1) Electric gas conversion equipment model
Figure QLYQS_1
Figure QLYQS_2
In the formula: t is scheduling time; m is an index of the electric-to-gas equipment;
Figure QLYQS_3
respectively gas making power, consumed power and electricity-to-gas efficiency, L, of an electricity-to-gas (P2G) facility HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure QLYQS_4
The minimum and maximum breathing power of the mth P2G;
(2) Electric boiler model
Figure QLYQS_5
Figure QLYQS_6
In the formula: t is scheduling time; n is an electric boiler index;
Figure QLYQS_7
and &>
Figure QLYQS_8
Respectively the power consumption and the heat production of the nth electric boiler in the t period;
Figure QLYQS_9
Is the electric heat conversion efficiency of the nth electric boiler>
Figure QLYQS_10
Dividing into the minimum and maximum heating power of the nth electric boiler;
Figure QLYQS_11
Starting and stopping the nth electric boiler at a time period t;
(3) Gas turbine model
Figure QLYQS_12
Figure QLYQS_13
Figure QLYQS_14
Figure QLYQS_15
Figure QLYQS_16
Figure QLYQS_17
Figure QLYQS_18
In the formula: t is scheduling time; q is a gas turbine index;
Figure QLYQS_20
and &>
Figure QLYQS_22
Respectively representing the power generation power and the gas consumption power of the gas turbine; f (-) represents the gas turbine energy consumption curve;
Figure QLYQS_25
And &>
Figure QLYQS_21
Respectively representing the consumption of natural gas required by the startup and shutdown of the gas turbine; a is q 、b q And c q A gas coefficient representing F ('); l is HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure QLYQS_24
Dividing the power into the minimum and maximum power generation powers of the q gas turbines;
Figure QLYQS_26
Representing the generated power of the qth gas turbine in the t-1 period;
Figure QLYQS_27
For the start-stop state of the qth gas turbine in the time period t, ->
Figure QLYQS_19
Starting and stopping states of a qth gas turbine in a t-1 time period;
Figure QLYQS_23
Is the climbing rate and the descending rate of the qth gas turbine>
Figure QLYQS_28
For the continuous start-up and shutdown time of the qth gas turbine in the t-1 time period, based on the preset time interval>
Figure QLYQS_29
The minimum startup and shutdown time of the qth gas turbine in the time period t is defined;
(4) Combined heat and power generation unit model
Figure QLYQS_30
Figure QLYQS_31
Figure QLYQS_32
Figure QLYQS_33
Figure QLYQS_34
Figure QLYQS_35
Figure QLYQS_36
In the formula: t is scheduling time; p is an index of the cogeneration unit;
Figure QLYQS_39
and &>
Figure QLYQS_45
Respectively representing heat production power and gas consumption power of a combined heat and power generation unit (CHP);
Figure QLYQS_48
The generating power and the generating efficiency of the micro-combustion engine in the t period are based on the pressure value>
Figure QLYQS_40
For the heat dissipation loss rate, the device can be used for>
Figure QLYQS_42
And &>
Figure QLYQS_46
The heating coefficient and the flue gas recovery rate of the bromine refrigerator are respectively; l is a radical of an alcohol HANG Taking 9.7kWh/m as low heating value of natural gas 3
Figure QLYQS_50
Dividing into the p (th) CHP unit minimum and maximum generating power;
Figure QLYQS_37
Representing the generated power of the p-th combined heat and power generation unit in a t-1 period;
Figure QLYQS_43
For the start-stop state of the pth CHP unit in the t-1 time period, the system is turned on or off>
Figure QLYQS_47
The starting and stopping state of the P-th cogeneration unit in the t-1 time period is shown;
Figure QLYQS_49
Figure QLYQS_38
The climbing rate and the descending rate of the pth CHP unit;
Figure QLYQS_41
for the continuous startup and shutdown time of the pth CHP unit in the t-1 period>
Figure QLYQS_44
The minimum startup and shutdown time of the pth CHP unit in the time period t is obtained;
(5) Energy storage equipment model
ESS(t)=(1-σ S )·ESS(t-1)+ES in (t)·η in -ES out (t)/η out
Figure QLYQS_51
Figure QLYQS_52
Figure QLYQS_53
Figure QLYQS_54
Figure QLYQS_55
Figure QLYQS_56
Figure QLYQS_57
Figure QLYQS_58
In the formula: t is scheduling time; ESS (t), ES in (t)、ES out (t) storing energy, stored energy power and released energy power of the energy storage device at a time period t, respectively; ESS (t-1) represents the stored energy of the energy storage device during the t-1 period; sigma S Is the self-consumption rate of the energy storage system; eta in 、η out Respectively storing energy and releasing energy efficiency for the energy storage equipment;
Figure QLYQS_59
gas storage, gas discharge power, GS, for a period of t in,max 、GS out,max The maximum gas storage and gas release power of the gas storage device are respectively;
Figure QLYQS_60
The air storage quantity of the air storage device is greater or less for a time period t>
Figure QLYQS_61
The gas storage capacity of the gas storage equipment is t-1 time; c GS,min 、C GS,max Respectively, the minimum and maximum gas storage capacities, eta, of the gas storage equipment CGS 、η GS,in 、η GS,out The self-consumption rate, the gas storage efficiency and the gas release efficiency of the gas storage equipment are improved;
Figure QLYQS_62
For the heat-storage and heat-release power of t time period, HS in ,max 、HS out,max The maximum heat storage power and the maximum heat release power of the heat storage equipment are respectively;
Figure QLYQS_63
For a gas reserve of the heat storage device for a period t>
Figure QLYQS_64
The gas storage capacity of the heat storage equipment is t-1 time period; c HS,min 、C HS,max Minimum and maximum heat storage capacity, eta, of the heat storage apparatus CHS 、η HS,in 、η HS,out The self consumption rate, the heat storage efficiency and the heat release efficiency of the heat storage equipment are obtained;
step 3, the demand response model is as follows:
Figure QLYQS_65
Figure QLYQS_66
Figure QLYQS_67
Figure QLYQS_68
Figure QLYQS_69
P t DR =P t DR,inte +P t DR,shif
Figure QLYQS_70
in the formula: t is scheduling time;
Figure QLYQS_71
P t DR 、P t DR,inte 、P t DR,shif 、P t LD,max predicted value of electric load, demand side response electric load, demand side interruptible electric load, and demand side transfer for t periodElectric load, maximum electric load allowed by the system, P inte,max The maximum interruptible electrical load within the scheduling time;
Figure QLYQS_72
A maximum interruptible and exportable electrical load ratio for the t period; p is t DR,shif Positive represents the load that can be shifted out, and vice versa represents the load that can be shifted in;
Figure QLYQS_73
Figure QLYQS_74
Figure QLYQS_75
in the formula: t is scheduling time;
Figure QLYQS_76
and &>
Figure QLYQS_77
Respectively a predicted value of the heat load, a response heat load of a demand side, a response heat load proportion of the demand side, the maximum heat load allowed by a system and the heat load after considering the demand response in a time period t; h DR,max The maximum interruptible thermal load within the scheduling time; n is a radical of t Scheduling the time period for the whole;
step 4 the two-stage tunable robust model considering the combined heat and power demand response is as follows:
(1) Objective function
Figure QLYQS_78
Figure QLYQS_79
Figure QLYQS_80
Figure QLYQS_81
Figure QLYQS_82
Figure QLYQS_83
Figure QLYQS_84
Figure QLYQS_85
Figure QLYQS_86
In the formula: t is scheduling time;
Figure QLYQS_88
the electricity and gas purchase costs are respectively;
Figure QLYQS_92
The starting-up and shutdown costs of the pth CHP are respectively;
Figure QLYQS_93
Penalty fee for wind abandon, based on the sum of the wind and the wind>
Figure QLYQS_89
Punishment of cost for light abandonment;
Figure QLYQS_91
To earn electricity sales, C E,cc A compensation cost to interrupt the electrical load for demand side response; p t in For purchasing power, is turned on or off>
Figure QLYQS_96
For purchasing air power, is selected>
Figure QLYQS_97
Wind abandoning power and light abandoning power of the ith fan and the jth group of photovoltaic cells in the t period are respectively set; p is t out To sell electric power, P t DR,inte For interrupting the electrical load on the demand side, Δ t is the scheduling time interval, N t For scheduling period number, ->
Figure QLYQS_87
The unit purchase price of electricity, gas and electricity, phi t WT 、φ t PV A penalty price for wind abandon and light abandon of the unit, respectively>
Figure QLYQS_90
For the compensation price of the unit demand side in response to an interrupted electrical load, based on the value of the compensation>
Figure QLYQS_94
The starting and stopping states of the pth CHP in the period t, 1 represents starting, 0 represents stopping,
Figure QLYQS_95
the cost of starting up and stopping the p-th CHP once respectively, N WT 、N PV 、N CHP The number of the fan, the photovoltaic cell, the electric boiler and the CHP are respectively;
(2) Constraint conditions
Figure QLYQS_98
Figure QLYQS_99
Figure QLYQS_100
Figure QLYQS_101
Figure QLYQS_102
Figure QLYQS_103
Figure QLYQS_104
Figure QLYQS_105
Figure QLYQS_106
Figure QLYQS_107
Figure QLYQS_108
Figure QLYQS_109
Figure QLYQS_110
P t u,DR =P t u,DR,inte +P t u,DR,shif :(λ 10_6,t )
Figure QLYQS_111
Figure QLYQS_112
Figure QLYQS_113
Figure QLYQS_114
In the formula: t is scheduling time; (. Cndot.) u Is a corresponding adjusted variable under the real-time change of wind and light output; v. of 1,t And v 2,t A power balance constraint relaxation variable; v. of 3,t And v 4,t Constraint relaxation variables for thermal balance; omega WT 、Ω PV Respectively representing uncertain collections of wind power output and photovoltaic output;
Figure QLYQS_115
deviation of wind power, photovoltaic output and predicted value respectively>
Figure QLYQS_122
The ratio of the wind power and photovoltaic output deviation value to the predicted value is obtained;
Figure QLYQS_117
Is a variable of 0 to 1 in the uncertain convergence set; delta i 、Δ j Respectively obtaining wind power and photovoltaic output uncertainty precalculated values;
Figure QLYQS_120
Correcting the climbing and descending of the pth CHP unit;
Figure QLYQS_121
Correcting the climbing and descending of the qth gas turbine unit; lambda (-) is a dual variable corresponding to a constraint condition>
Figure QLYQS_123
And &>
Figure QLYQS_116
Respectively taking the average values of the predicted values of the wind power output and the photovoltaic output;
Figure QLYQS_119
And &>
Figure QLYQS_124
Respectively the abandoned wind power and the abandoned light power after the uncertainty parameters are considered;
Figure QLYQS_125
Representing the wind abandoning proportion;
and 5, an abstract expression of the two-stage adjustable robust optimization model of the park comprehensive energy system economic dispatch in the day ahead is as follows:
Figure QLYQS_126
s.t.Ax+By≤b
Figure QLYQS_127
in the formula: x represents the starting and stopping states of units related to the CHP and the gas turbine, y and z represent the basic scene and the scheduling output of other units of the system adjusted according to the wind-light output transformation, u is an uncertain variable related to the uncertainty of the wind power output and the photovoltaic output, and c b 、c g A, B, B, F, h, C, D, E, F, G can be derived from the objective function and constraint conditions, C b And c g The method comprises the following steps that A, B, C, D, E, F and G are abstract coefficient vector matrixes of variables in inequality constraints respectively; b. h represents an abstract matrix of constants in the inequality constraint; f is an abstract coefficient vector corresponding to a max-min double-layer problem objective function;
step 6, the maximum and minimum subproblems of the worst scene recognition are as follows:
ΔD=ΔD 1 +ΔD 2
Figure QLYQS_128
Figure QLYQS_129
s.t.λ 9_1,t ≤0,λ 9_2,t ≤0,λ 9_3,t ≤0,λ 9_4,t ≤0,λ 9_5,t ≤0,λ 9_6,t ≤0,t∈N T
λ 6_2,q,t ≤0,λ 6_3,q,t ≤0,λ 6_4,q,t ≤0,λ 6_6,q,t ≤0,λ 6_7,q,t ≤0,q∈N GT ,t∈N T
λ 7_2,p,t ≤0,λ 7_3,p,t ≤0,λ 7_4,p,t ≤0,λ 7_6,p,t ≤0,λ 7_7,p,t ≤0,p∈N CHP ,t∈N T
λ 8_1,t ≤0,λ 8_2,t ≤0,λ 8_3,t ≤0,λ 8_4,t ≤0,t∈N T
λ 13_1,t ≤0,λ 13_2,t ≤0,λ 13_3,t ≤0,λ 13_4,t ≤0,t∈N T
λ 4_2,m,t ≤0,λ 4_3,m,t ≤0,t∈N T ,m∈N P2G ,t∈N T
λ 5_2,n,t ≤0,λ 5_3,n,t ≤0,t∈N T ,n∈N EB ,t∈N T
λ 10_1,t ≤0,λ 10_2,t ≤0,λ 10_3,t ≤0,λ 10_4,t ≤0,t∈N T
λ 11_3,t ≤0,λ 11_3,t ≤0,λ 12_1,t ≤0,λ 12_2,t ≤0,t∈N T
λ 10_7 ≤0,λ 10_8 ≤0,λ 11_4 ≤0,λ 11_5 ≤0,t∈N T
Figure QLYQS_130
Figure QLYQS_131
s.t.-λ 1,t6_1,q,t ·b q6_2,q,t6_3,q,t6_4,q,t+1 +
λ 6_4,q,t6_5,p,t+16_5,p,t6_6,q,t6_7,q,t ≤0
Figure QLYQS_132
Figure QLYQS_133
1,t9_1,t9_2,t ≤0,t∈N T
λ 1,t9_3,t9_4,t ≤0,t∈N T
λ 1,t10_5,t ≤0,t∈N T
10_1,t10_2,t10_3,t10_5,t10_710_8 =0,t∈N T
10_1,t10_4,t10_5,t10_5,t10_6 =0,t∈N T
λ 1,t12_1,i,t ≤0,t∈N T ,i∈N WT
λ 1,t12_2,j,t ≤0,t∈N T ,j∈N PV
2,t7_1,p,t ≤0,t∈N T ,p∈N CHP
2,t6_1,q,t ≤0,t∈N T ,q∈N GT
2,t4_1,m,t4_2,m,t4_3,m,t ≤0,t∈N T ,m∈N P2G
λ 2,t9_5,t9_6,t ≤0,t∈N T
λ 2,t8_2,t8_5,tGS,out ≤0,t∈N T
2,t8_1,t8_5,t ·η GS,in ≤0,t∈N T
8_3,t8_4,t8_5,t+1 ·η GS,in8_5,t ≤0,t∈N T3,t ·η h13_2,t13_5,tHS,out ≤0,t∈N T
λ 3,t ·η h13_1,t13_5,t ·η HS,in ≤0,t∈N T13_3,t13_4,t13_5,t+1 ·η HS,in13_5,t ≤0,t∈N T
3,t ·η h5_1,n,t5_2,n,t5_3,n,t ≤0,t∈N T ,n∈N EB3,t ·η h7_8,n,t ≤0,t∈N T ,n∈N EB
λ 3,t11_1,t ≤0,t∈N T λ 11_1,t11_2,t11_3,t11_411_5,t ≤0,t∈N T
-1≤λ 1,t ≤1,-1≤λ 3,t ≤1,t∈N T
and 7, solving a two-stage adjustable robust model considering combined heat and power demand response by using a CCG method as follows:
(1) Major problems
Figure QLYQS_134
Ax+By≤b
(2) Sub-problems
Figure QLYQS_135
Figure QLYQS_136
(3) CCG solving step
Step 1: let the iteration counter s =0 set the maximum value epsilon allowed by the system for violating the safety regulations RO
Step 2: solving the main problem, if the main problem is solved, obtaining a system unit starting and stopping state x and a unit output arrangement y, and performing the step 3; otherwise, stopping iteration and outputting no solution;
and 3, step 3: solving the maximum and minimum subproblems according to the x and y obtained by the solving in the step 2, and finding out the wind power and photovoltaic output under the worst scene which causes the maximum possibility of violating the safety specified value;
and 4, step 4: if solved in step 3Maximum possible violation of the safety prescribed value less than epsilon RO Then x and y are the final optimization solution and the iteration is stopped; otherwise, let s = s +1, according to the wind power and photovoltaic output value under the worst scene solved in step 3
Figure QLYQS_137
Adding CCG constraint shown below into the main problem, and returning to the step 2;
and 8, the park comprehensive energy system data further comprises the specific composition of the park comprehensive energy system, the energy price, the equipment parameters and values of each energy conversion equipment, the demand response proportion, the new energy output fluctuation situation under the basic scene and the worst scene, the maximum violation safety specified value and the predicted value of the new energy output and load.
CN202010773898.1A 2020-08-04 2020-08-04 Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty Active CN113256045B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010773898.1A CN113256045B (en) 2020-08-04 2020-08-04 Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010773898.1A CN113256045B (en) 2020-08-04 2020-08-04 Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty

Publications (2)

Publication Number Publication Date
CN113256045A CN113256045A (en) 2021-08-13
CN113256045B true CN113256045B (en) 2023-04-07

Family

ID=77220096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010773898.1A Active CN113256045B (en) 2020-08-04 2020-08-04 Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty

Country Status (1)

Country Link
CN (1) CN113256045B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113705991B (en) * 2021-08-17 2023-08-29 国网四川省电力公司技能培训中心 Establishment and low-carbon scheduling method for multi-energy park
CN114519449A (en) * 2021-12-01 2022-05-20 中国华能集团有限公司河北雄安分公司 Operation optimization method for park energy system
CN114784797B (en) * 2022-04-25 2024-01-19 东南大学溧阳研究院 Thermoelectric optimization day-ahead dispatching method for residence comprehensive energy system considering multiple uncertainties
CN114912714B (en) * 2022-06-17 2023-11-07 国网江苏省电力有限公司苏州供电分公司 Low-carbon economic dispatching method and system considering wind-light output uncertainty under lightning climate
CN115689242A (en) * 2022-11-15 2023-02-03 云南电网有限责任公司电力科学研究院 Scheduling method of multi-energy system, computer device and computer-readable storage medium
CN117039901B (en) * 2023-08-10 2024-07-05 国网山东省电力公司威海供电公司 Station area uncertainty on-site management scheduling method based on low-voltage flexible direct current interconnection

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108832665A (en) * 2018-07-04 2018-11-16 四川大学 A kind of probabilistic electric heating integrated system Robust distributed coordination optimization scheduling model of consideration wind-powered electricity generation
CN110009152A (en) * 2019-04-03 2019-07-12 东南大学 A kind of consideration electricity turns gas and probabilistic regional complex energy system operation robust Optimal methods
CN110555595A (en) * 2019-08-09 2019-12-10 四川大学 biogas-wind-light all-renewable energy system based on energy hub and method thereof
CN111355265A (en) * 2020-04-10 2020-06-30 中国人民解放军国防科技大学 Micro-grid energy two-stage robust optimization method and system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108832665A (en) * 2018-07-04 2018-11-16 四川大学 A kind of probabilistic electric heating integrated system Robust distributed coordination optimization scheduling model of consideration wind-powered electricity generation
CN110009152A (en) * 2019-04-03 2019-07-12 东南大学 A kind of consideration electricity turns gas and probabilistic regional complex energy system operation robust Optimal methods
CN110555595A (en) * 2019-08-09 2019-12-10 四川大学 biogas-wind-light all-renewable energy system based on energy hub and method thereof
CN111355265A (en) * 2020-04-10 2020-06-30 中国人民解放军国防科技大学 Micro-grid energy two-stage robust optimization method and system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Chuan He 等.Co-optimization scheduling of interdependent power and gas systems with electricity and gas uncertainties.Energy.2018,第159卷1003-1015. *
Jun Lu 等.Robust day-ahead coordinated scheduling of multi-energy systems with integrated heat-electricity demand response and high penetration of renewable energy.Renewable Energy.2021,第178卷466-482. *
单福州 等.基于改进两阶段鲁棒优化的区域综合能源系统经济调度.电测与仪表.2018,第55卷(第23期),103-108. *
税月 等.考虑风电不确定性的电气能源系统两阶段分布鲁棒协同调度.电力系统自动化.2018,第42卷(第13期),43-50,75. *

Also Published As

Publication number Publication date
CN113256045A (en) 2021-08-13

Similar Documents

Publication Publication Date Title
CN113256045B (en) Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty
Zhou et al. Optimized operation method of small and medium-sized integrated energy system for P2G equipment under strong uncertainty
CN109004686B (en) Cold, heat and power combined supply type micro-grid system considering ice storage air conditioner multi-mode
CN104734168B (en) Microgrid running optimization system and method based on power and heat combined dispatching
CN111950807B (en) Comprehensive energy system optimization operation method considering uncertainty and demand response
CN112701687B (en) Robust optimization operation method of gas-electricity distribution network system considering price type combined demand response
CN106786793A (en) A kind of supply of cooling, heating and electrical powers type microgrid operation method based on robust optimization
CN107358345B (en) Distributed combined cooling heating and power system optimization operation method considering demand side management
CN108090666B (en) AA-CAES-containing power grid electric energy and reserve capacity cooperative scheduling method
CN104537443A (en) Co-generation type micro-grid economy coordination and optimization dispatching method
CN113779783A (en) Multi-uncertainty-considered planning and operation joint optimization method for regional comprehensive energy system
CN112836882B (en) Regional comprehensive energy system operation optimization method considering equipment load rate change
CN110620403A (en) Day-ahead scheduling method and system for collaborative operation of energy system considering renewable energy
CN112446546B (en) Comprehensive energy system two-stage optimal configuration method considering energy reliability
CN109286208A (en) A kind of integrated energy system dispatching method and system
CN115099007B (en) Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve
CN113469430A (en) Multi-energy complementary capacity configuration method of integrated energy park
CN115860413A (en) Grid-connected combined cooling heating and power micro-grid system economic scheduling method based on load demand response and double-layer adjustable robust optimization
Jie et al. Study of energy scheduling and optimal cost management of a new structure CCHP system: A case study supplying energy for a chemical enterprise in Jiangsu Province
CN111030101A (en) Clean energy consumption linkage regulation and control method and system based on diversified big data
Gao et al. Fuzzy decision-based optimal energy dispatch for integrated energy systems with energy storage
CN113240205A (en) Regional energy utilization system substitution optimization method based on multi-energy comprehensive utilization
Cheng et al. A stochastic optimal model of micro energy internet contains rooftop PV and CCHP system
CN114091762A (en) Energy system double-layer operation optimization method and system based on stepped carbon trading
CN111769561A (en) RIES operation optimization method and device

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