CN104616208A - Model predication control based cooling heating and power generation type micro-grid operation method - Google Patents
Model predication control based cooling heating and power generation type micro-grid operation method Download PDFInfo
- Publication number
- CN104616208A CN104616208A CN201510059153.8A CN201510059153A CN104616208A CN 104616208 A CN104616208 A CN 104616208A CN 201510059153 A CN201510059153 A CN 201510059153A CN 104616208 A CN104616208 A CN 104616208A
- Authority
- CN
- China
- Prior art keywords
- represent
- unit
- power
- period
- tst
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001816 cooling Methods 0.000 title claims abstract description 36
- 238000010438 heat treatment Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000010248 power generation Methods 0.000 title abstract 5
- 238000005096 rolling process Methods 0.000 claims abstract description 25
- 238000012544 monitoring process Methods 0.000 claims abstract description 14
- 239000007789 gas Substances 0.000 claims description 224
- 238000005338 heat storage Methods 0.000 claims description 174
- 239000003463 adsorbent Substances 0.000 claims description 66
- 238000009825 accumulation Methods 0.000 claims description 48
- 238000005057 refrigeration Methods 0.000 claims description 38
- 238000007599 discharging Methods 0.000 claims description 37
- 230000005611 electricity Effects 0.000 claims description 37
- 238000005457 optimization Methods 0.000 claims description 34
- 238000012423 maintenance Methods 0.000 claims description 28
- 238000003860 storage Methods 0.000 claims description 24
- 238000004146 energy storage Methods 0.000 claims description 19
- 238000012937 correction Methods 0.000 claims description 10
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 10
- 239000002918 waste heat Substances 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 238000004378 air conditioning Methods 0.000 claims description 6
- 239000003345 natural gas Substances 0.000 claims description 5
- 239000000446 fuel Substances 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011022 operating instruction Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- -1 sulphuric dioxide Substances 0.000 description 1
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
- 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems 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)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a model predication control based cooling heating and power generation type micro-grid operation method. The method comprises the steps of building a predication model to predicate the wind power, photovoltaic power and cooling heating power generation load power data within a control time domain in the future; solving a rolling optimizing model at every moment according to the latest predication result and the real-time operation state of each device; calculating the force output of each device at the following time quantum; re-calculating the rolling optimizing model at the next time according to the latest predication result, monitoring the actual wind power value, photovoltaic value, cooling load value, heating load value and power generation load value on real time; updating the historic data; solving and feeding back to correct the model at every 5 minutes to obtain the adjustment of each device; distributing to each device for adjusting until reaching the next rolling optimizing time. With the adoption of the method, the influence of inaccurate predication on the optimal operation of the system can be removed well, thus the system operation risk is reduced, and the system operation stability is improved; meanwhile, the economy of the operation of a cooling heating and power generation type micro-grid can be increased.
Description
Technical field
The invention belongs to cold, heat and power triple supply system field, specifically, relate to a kind of supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control.
Background technology
At present, take natural gas as the clean energy resource of representative, use with regenerative resource and be combined with each other, and utilize generator waste heat overbottom pressure, together constitute supply of cooling, heating and electrical powers type micro-capacitance sensor, become the outer fast-developing cleaning new energy industry of Now Domestic.Supply of cooling, heating and electrical powers type micro-grid system is arranged on user side nearby with various ways such as small-scale, low capacity, modularization and distributings, the energy supply system of user to cold energy, heat energy, electrical energy demands can be met simultaneously, save tradition and concentrate the initial outlay cost of energy induction system and reduce energy loss in transmitting procedure.Supply of cooling, heating and electrical powers type micro-grid system take clean energy resource as main fuel, can reduce the pollutant emissions such as oxides of nitrogen, sulphuric dioxide, solid particle, greatly reduces the pressure that environment worsens.
Supply of cooling, heating and electrical powers type micro-capacitance sensor optimizing operation method is basis operation plan instruction operation a few days ago mostly, can realize the economical operation of system to a certain extent.But, exert oneself due to renewable power supply and there is the feature such as undulatory property, intermittence, add that the electricity consumption behavior of user has more randomness, often there is deviation in the predicted value of load and renewable power supply power and actual value, cause actual conditions often to depart from operation plan a few days ago, have a strong impact on supply of cooling, heating and electrical powers type micro-capacitance sensor optimizing operation.
Summary of the invention
Technical matters: technical matters to be solved by this invention is: a kind of supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control is provided, the impact because forecasting inaccuracy causes system cloud gray model can be eliminated well, ensure what system operating instructions always calculated in the optimization of up-to-date information, reduce the risk of system cloud gray model, improve the stability of system cloud gray model, improve the economy that supply of cooling, heating and electrical powers type micro-capacitance sensor runs simultaneously.
Technical scheme: for solving the problems of the technologies described above, a kind of supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control that the present invention proposes, comprises the following steps:
Step 10) set up parametric prediction model, comprise wind-powered electricity generation forecast model, photovoltaic forecast model, cooling load prediction model, heat load prediction model and electric load forecast model; According to each parametric prediction model, in the k moment, according to n the historical data gathered before the k moment, utilize each parametric prediction model, following wind power, photovoltaic power, refrigeration duty power, thermal load power and the electric load power controlled in time domain M of prediction; Historical data comprises wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power;
Step 20) according to step 10) control in time domain M wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power in future of obtaining, set up such as formula the rolling optimization model shown in (1) to formula (12):
First determine with micro-grid system operating cost minimum for objective function, shown in (1):
In formula, C represents system operation cost; K represents current time; M represents control time domain;
the fuel cost of expression system t period;
the operation expense of expression system t period;
represent t period and the mutual cost of electrical network;
Then determine constraint condition, comprise cold energy balance, thermal energy balance, electric energy balance and equipment and run constraint condition:
Cold energy equilibrium constraint is determined according to formula (5):
In formula,
represent the power input of adsorbent refrigerator t period, unit: kW; COP
acrepresent the coefficient of refrigerating performance of adsorbent refrigerator;
represent the power input of electric refrigerating machine t period, unit: kW; COP
ecrepresent the coefficient of refrigerating performance of electric refrigerating machine;
represent t cooling load of the air-conditioning system power, unit: kW;
Thermal energy equilibrium constraint is determined according to formula (6):
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; η
mtrepresent the efficiency of miniature gas turbine; η
lossrepresent the gas turbine energy proportion of goods damageds; η
hrrepresent waste-heat recoverer efficiency;
represent the power of gas fired-boiler t period, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat release power of heat storage tank t period, unit: kW;
represent the power input of adsorbent refrigerator t period, unit: kW;
represent the thermal load power of t period system, unit: kW; η
herepresent effectiveness of heat exchanger;
Electric energy balance constraint condition is determined according to formula (7):
In formula,
represent the predicted value of photovoltaic t period;
represent the predicted value of wind-powered electricity generation t period;
represent the electric power of miniature gas turbine t period, unit: kW;
expression system t period and the mutual power of major network, unit: kW;
represent the charge power of accumulator t period, unit: kW;
represent accumulator t period discharge power, unit: kW;
represent the power input of electric refrigerating machine t period, unit: kW;
represent t period system electric load power, unit: kW;
Determine that miniature gas turbine runs constraint condition according to formula (801) and formula (802), its Chinese style (801) is miniature gas turbine running status constraint condition; Formula (802) represents gas turbine unit Climing constant condition, comprises Unit Commitment Climing constant and runs Climing constant continuously:
In formula (801),
represent miniature gas turbine t period running status variable,
represent that miniature gas turbine runs,
represent that miniature gas turbine is shut down;
represent the lower limit that miniature gas turbine is exerted oneself, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW;
represent the upper limit that miniature gas turbine is exerted oneself, unit: kW;
represent that unit falls power, in maximum when continuous running status of micro-gas-turbine unit: kW;
represent miniature gas turbine t-1 period running status variable;
to represent when micro-gas-turbine unit is shut down maximum falls power, unit: kW;
represent the electric power of miniature gas turbine t-1 period, unit: kW;
represent most the increase power of micro-gas-turbine unit when continuous running status, unit: kW;
represent the power that increases most when micro-gas-turbine unit starts, unit: kW;
The constraint condition that gas fired-boiler runs is determined according to formula (9):
In formula,
represent the lower limit that gas fired-boiler is exerted oneself, unit: kW;
represent the power of gas fired-boiler t period, unit: kW;
represent the upper limit that gas fired-boiler is exerted oneself, unit: kW;
The constraint condition of the mutual power of electrical network is determined according to formula (10):
In formula,
the lower limit of expression system and the mutual power of major network, unit kW;
the upper limit of expression system and the mutual power of major network, unit: kW;
expression system t period and the mutual power of major network, unit: kW;
The constraint condition that accumulator runs is determined according to formula (11):
In formula,
represent the discharge condition of accumulator t period,
represent battery discharging;
represent that accumulator does not charge also not discharge;
represent charge in batteries power upper limit, unit: kW;
represent accumulator t period discharge power, unit: kW;
represent the charge power of accumulator t period, unit: kW;
represent the charged state of accumulator t period,
represent charge in batteries;
represent that accumulator does not charge also not discharge;
represent battery discharging power upper limit, unit: kW;
represent the energy of t+1 period in accumulator, unit: kWh;
represent the energy of t period in accumulator, unit: kWh; σ
btrepresent the self-energy proportion of goods damageds of accumulator; η
c, btrepresent the charge efficiency of accumulator; η
disc, btrepresent battery discharging efficiency;
represent the lower limit of accumulator storage power, unit: kWh;
represent the upper limit of accumulator storage power, unit: kWh; Δ t represents the time interval;
The constraint condition that heat storage tank runs is determined according to formula (12):
In formula,
represent the heat release state of heat storage tank t period,
represent heat storage tank heat release,
represent heat storage tank not heat release also not accumulation of heat;
represent the accumulation of heat power upper limit of heat storage tank, unit kW;
represent the heat release power of heat storage tank t period, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat storage state of heat storage tank t period,
represent heat storage tank accumulation of heat,
represent heat storage tank not heat release also not accumulation of heat;
represent the heat release power upper limit of heat storage tank, unit: kW;
represent the energy of t+1 period in heat storage tank, unit: kWh;
represent the energy of t period in heat storage tank, unit: kWh; σ
tstrepresent the self-energy proportion of goods damageds of heat storage tank; η
c, tstrepresent the heat storage efficiency of heat storage tank; η
disc, tstrepresent the efficiency of heat storage tank release heat;
represent the lower limit of heat storage tank storage power, unit: kWh;
represent the upper limit of heat storage tank storage power, unit: kWh;
In each moment, according to step 10) the up-to-date refrigeration duty that obtains, thermal load and electric load predict the outcome and miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, the running status that energy storage device is real-time, Yalmip optimization tool is adopted to solve rolling optimization model, measuring and calculating miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, energy storage device will control exerting oneself of time domain M period in future, by the miniature gas turbine of first period in the following M period, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, each equipment is delivered in the instruction of exerting oneself of energy storage device,
Step 30) monitoring real time data, and upgrade historical data: monitoring is obtained the wind power in each moment, photovoltaic power, refrigeration duty power, thermal load power and electric load power actual value, replaced the wind power in a upper moment, photovoltaic power, refrigeration duty power, thermal load power and electric load power actual value;
Step 40) set up feedback compensation model, the actual of equipment such as real-time correction miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, energy storage device are exerted oneself: the actual value of Real-Time Monitoring wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power, and the historical data upgrading wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load and electric load; Feedback compensation model is such as formula shown in (13) to formula (21):
Set up such as formula shown in (13) with the minimum correction function for objective function of relative adjustment amount:
In formula, AD represents overall relative adjustment amount; △ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
represent the rated power that miniature gas turbine exports, unit: kW; △ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW;
represent the rated power of electric refrigerating machine machine input, unit: kW; △ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent the rated power of charge in batteries, unit: kW; △ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the rated power of battery discharging, unit: kW; △ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
expression system and the mutual rated power of major network, unit: kW; △ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent that gas fired-boiler exports rated power, unit: kW; △ P
acrepresent the adjustment amount of adsorbent refrigerator power input, unit: kW;
represent the rated power of adsorbent refrigerator input, unit: kW; △ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the rated power of heat storage tank accumulation of heat, unit: kW; △ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the rated power of heat storage tank heat release, unit: kW; w
1represent the weight coefficient relevant to electric energy, w
2represent the weight coefficient relevant to cold and hot energy;
Set up such as formula the cold energy equilibrium constraint shown in (14):
In formula,
represent the power input of adsorbent refrigerator t period, unit: kW; △ P
acrepresent the adjustment amount of adsorbent refrigerator power input; COP
acrepresent the coefficient of refrigerating performance of adsorbent refrigerator;
represent the power input of electric refrigerating machine t period, unit: kW; △ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW; COP
ecrepresent the coefficient of refrigerating performance of electric refrigerating machine; Q
cfor real-time cooling load of the air-conditioning system power, unit: kW;
Set up such as formula the thermal energy equilibrium constraint shown in (15):
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; △ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW; η
mtrepresent the efficiency of miniature gas turbine; η
lossrepresent the gas turbine energy proportion of goods damageds; η
hrrepresent waste-heat recoverer efficiency;
represent the power of gas fired-boiler t period, unit: kW; △ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kW; △ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the heat release power of heat storage tank t period, unit: kW; △ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the power input of adsorbent refrigerator t period, unit: kW; △ P
acrepresent the adjustment amount of adsorbent refrigerator power input; η
herepresent effectiveness of heat exchanger; Q
hfor real-time system heat load power, unit: kW;
Set up such as formula the electric energy balance constraint condition shown in (16):
In formula, P
pvrepresent real-time photovoltaic power, unit: kW; P
wtrepresent real-time wind power, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW; △ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
expression system t period and the mutual power of major network, unit: kW; △ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
represent the charge power of accumulator t period, unit: kW; △ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent accumulator t period discharge power, unit: kW; △ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the power input of electric refrigerating machine t period, unit: kW; △ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW; P
elrepresent real-time system electric load power, unit: kW;
Set up the constraint condition run such as formula the miniature gas turbine shown in (17):
In formula,
represent miniature gas turbine t period running status variable,
represent that miniature gas turbine runs,
represent that miniature gas turbine does not run;
represent the lower limit that miniature gas turbine is exerted oneself, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW; △ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
represent the upper limit that miniature gas turbine is exerted oneself, unit: kW;
represent that unit falls power, in maximum in short-term when continuous running status of micro-gas-turbine unit: kW;
represent in short-term most the increase power of micro-gas-turbine unit when continuous running status, unit: kW;
Set up the constraint condition run such as formula the gas fired-boiler shown in (18):
In formula,
represent the lower limit that gas fired-boiler is exerted oneself, unit: kW;
represent the power of gas fired-boiler t period, unit: kW; △ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent the upper limit that gas fired-boiler is exerted oneself, unit: kW;
Set up the constraint condition such as formula the mutual power of electrical network shown in (19):
formula (19)
In formula,
the lower limit of expression system and the mutual power of major network, unit: kW;
expression system t period and the mutual power of major network, unit: kW; △ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
the upper limit of expression system and the mutual power of major network, unit: kW;
Set up the constraint condition run such as formula the accumulator shown in (20):
In formula,
represent the discharge condition of accumulator t period,
represent battery discharging,
represent that accumulator does not charge also not discharge;
represent charge in batteries power upper limit, unit: kW;
represent accumulator t period discharge power, unit: kW; △ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the charge power of accumulator t period, unit: kW; △ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent the charged state of accumulator t period,
represent charge in batteries,
represent that accumulator does not charge also not discharge;
represent battery discharging power upper limit, unit: kW;
energy before expression feedback compensation in accumulator, unit: kWh;
represent the energy in the accumulator after feedback compensation, unit: kWh; σ
btrepresent the self-energy proportion of goods damageds of accumulator; η
c, btrepresent the charge efficiency of accumulator; η
disc, btrepresent battery discharging efficiency;
represent the lower limit of accumulator storage power, unit: kWh;
represent the upper limit of accumulator storage power, unit: kWh; Δ t '=1/12h;
Set up the constraint condition run such as formula the heat storage tank shown in (21):
In formula,
represent the heat release state of heat storage tank t period,
represent heat storage tank heat release,
represent heat storage tank not heat release also not accumulation of heat;
represent charge in batteries power upper limit, unit: kW;
represent the heat release power of heat storage tank t period, unit: kW; △ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh; △ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the heat storage state of heat storage tank t period,
represent heat storage tank accumulation of heat,
represent heat storage tank not heat release also not accumulation of heat;
represent the heat release power upper limit of heat storage tank, unit: kW;
energy before expression feedback compensation in heat storage tank, unit: kWh;
represent the energy in the heat storage tank after feedback compensation, unit: kWh; σ
tstrepresent the self-energy proportion of goods damageds of heat storage tank; η
c, tstrepresent the heat storage efficiency of heat storage tank; η
disc, tstrepresent the efficiency of heat storage tank release heat;
represent the lower limit of heat storage tank storage power, unit: kWh;
represent the upper limit of heat storage tank storage power, unit: kWh;
Finally, Yalmip optimization tool is adopted to solve feedback compensation model, obtain the adjustment amount that miniature gas turbine is exerted oneself, gas fired-boiler is exerted oneself, adsorbent refrigerator power input, electric refrigerating machine power input, accumulator cell charging and discharging power, heat storage tank store the mutual power of heat release power, system and electrical network, these adjustment amounts are issued to respectively miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, accumulator, heat storage tank equipment adjusts; Within every 5 minutes, perform a step 30) and step 40), until be finished in control cycle Δ t.
Step 50) enter subsequent time, return step 10), until micro-capacitance sensor out of service.
Further, described step 20) in,
In formula, △ t represents the time interval; R
ngrepresent Gas Prices, unit: $/m
3;
represent the electric power of miniature gas turbine t period, unit: kW; η
mtrepresent the efficiency of miniature gas turbine;
represent the power of gas fired-boiler t period, unit: kW; η
brepresent the efficiency of gas fired-boiler; H
ngrepresent heating value of natural gas;
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; K
om, mtrepresent miniature gas turbine operation and maintenance cost, unit: $/kWh;
represent the power of gas fired-boiler t period, unit: kW; K
om, brepresent gas fired-boiler operation and maintenance cost, unit: $/kWh;
represent the thermal load power of t period system, unit: kW; η
herepresent effectiveness of heat exchanger; K
om, herepresent heat exchanger operation and maintenance cost, unit: $/kWh;
represent the power input of adsorbent refrigerator t period, unit: kW; K
om, acrepresent adsorbent refrigerator operation and maintenance cost, unit: $/kWh;
represent the power input of electric refrigerating machine t period, unit: kW; K
om, ecrepresent electric refrigerating machine operation and maintenance cost, unit: $/kWh;
represent the predicted value of photovoltaic t period; K
om, pvrepresent photovoltaic cell maintenance cost unit: $/kWh;
represent the predicted value of wind-powered electricity generation t period; K
om, wtrepresent blower fan maintenance cost unit: $/kWh;
represent the charge power of accumulator t period, unit: kW;
represent accumulator t period discharge power, unit: kW; K
om, btrepresent accumulator operation and maintenance cost, unit: $/kWh;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat release power of heat storage tank t period, unit: kW; K
om, tstrepresent heat storage tank operation and maintenance cost, unit: $/kWh;
In formula,
the expression system t period from the price of major network power purchase, unit: $/kWh;
the expression system t period is to the price of major network sale of electricity; α be 0 or 1, α=1 represent that micro-capacitance sensor can to major network sale of electricity, α=0 represents and does not allow micro-capacitance sensor to major network sale of electricity;
expression system t period and the mutual power of major network, unit: kW,
represent from major network power purchase,
represent to major network sale of electricity.
Beneficial effect: compared with prior art, the present invention has the following advantages:
The supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control that the present invention proposes, first the corresponding forecast model of wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load and electric load is set up, according to up-to-date wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load, electric load historical data, prediction will control wind power in time domain M, photovoltaic power, cool and thermal power load power data future, each moment, according to the up-to-date refrigeration duty obtained, thermal load and electric load predict the outcome and miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, the running status that energy storage device etc. are real-time, Yalmip optimization tool is adopted to solve this rolling optimization model, calculate miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, the equipment such as energy storage device are exerted oneself the following M period, but miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, the equipment such as energy storage device only perform the instruction of first period, predict the outcome according to up-to-date again to subsequent time, recalculate rolling optimization model, circulation like this is carried out, this ensure that system call instruction is all calculating on the basis of up-to-date information, thus make each equipment run the actual optimum of maintenance, Real-Time Monitoring wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load, the actual value of electric load, and upgrade wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load, the historical data of electric load, Yalmip optimization tool within every 5 minutes, is adopted to solve feedback compensation model, obtain miniature gas turbine to exert oneself, gas fired-boiler is exerted oneself, adsorbent refrigerator power input, electric refrigerating machine power input, accumulator cell charging and discharging power, heat storage tank stores heat release power, the adjustment amount of system and the mutual power of electrical network etc., these adjustment amounts are issued to miniature gas turbine respectively, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, accumulator, the equipment such as heat storage tank adjust, until the next rolling optimization moment.The inventive method can eliminate the impact because forecasting inaccuracy causes system optimized operation well, ensure what system operating instructions always calculated in the optimization of up-to-date information, reduce the risk of system cloud gray model, improve the stability of system cloud gray model, improve the economy that supply of cooling, heating and electrical powers type micro-capacitance sensor runs simultaneously.
Accompanying drawing explanation
Fig. 1 is supply of cooling, heating and electrical powers type micro-capacitance sensor structural representation in the present invention.
Fig. 2 is Systematical control block diagram of the present invention.
Fig. 3 is process flow diagram of the present invention.
Embodiment
In order to make object of the present invention, technical scheme and advantage clearly understand, below in conjunction with accompanying drawing and case study on implementation, the present invention is in depth described in detail.Should be appreciated that concrete case study on implementation described herein is only in order to explain the present invention, and be not used in restriction invention.
The present invention relates to a kind of supply of cooling, heating and electrical powers type micro-capacitance sensor optimizing operation method based on Model Predictive Control, belong to cold, heat and power triple supply system field, is supply of cooling, heating and electrical powers type micro-capacitance sensor structural representation as shown in Figure 1.Supply of cooling, heating and electrical powers type micro-capacitance sensor comprises photovoltaic cell, blower fan, accumulator, miniature gas turbine, gas fired-boiler, electric refrigerating machine, Absorption Refrigerator, heat interaction device and heat storage tank, there is cold energy stream in depositing in system, thermal energy stream, electric flux stream and natural gas flow.The inventive method mainly comprises three parts: prediction module, rolling optimization module and feedback compensation module.
As shown in Figure 3, the supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control of the present invention, comprises the following steps:
Step 10) set up parametric prediction model, comprise wind-powered electricity generation forecast model, photovoltaic forecast model, cooling load prediction model, heat load prediction model and electric load forecast model; According to each parametric prediction model, in the k moment, according to n the historical data gathered before the k moment, utilize each parametric prediction model, following wind power, photovoltaic power, refrigeration duty power, thermal load power and the electric load power controlled in time domain M of prediction; Historical data comprises wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power.
Step 20) according to step 10) control in time domain M wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power in future of obtaining, set up such as formula the rolling optimization model shown in (1) to formula (12).
First determine with micro-grid system operating cost minimum for objective function, shown in (1):
In formula, C represents system operation cost; K represents current time; M represents control time domain;
the fuel cost of expression system t period;
the operation expense of expression system t period;
represent t period and the mutual cost of electrical network.
In formula, △ t represents the time interval; R
ngrepresent Gas Prices, unit: $/m
3;
represent the electric power of miniature gas turbine t period, unit: kW; η
mtrepresent the efficiency of miniature gas turbine;
represent the power of gas fired-boiler t period, unit: kW; η
brepresent the efficiency of gas fired-boiler; H
ngrepresent heating value of natural gas, H
ng=9.78kWh/m
3.As preferably, △ t=0.25h.
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; K
om, mtrepresent miniature gas turbine operation and maintenance cost, unit: $/kWh;
represent the power of gas fired-boiler t period, unit: kW; K
om, brepresent gas fired-boiler operation and maintenance cost, unit: $/kWh;
represent the thermal load power of t period system, unit: kW; η
herepresent effectiveness of heat exchanger; K
om, herepresent heat exchanger operation and maintenance cost, unit: $/kWh;
represent the power input of adsorbent refrigerator t period, unit: kW; K
om, acrepresent adsorbent refrigerator operation and maintenance cost, unit: $/kWh;
represent the power input of electric refrigerating machine t period, unit: kW; K
om, ecrepresent electric refrigerating machine operation and maintenance cost, unit: $/kWh;
represent the predicted value of photovoltaic t period; K
om, pvrepresent photovoltaic cell maintenance cost unit: $/kWh;
represent the predicted value of wind-powered electricity generation t period; K
om, wtrepresent blower fan maintenance cost unit: $/kWh;
represent the charge power of accumulator t period, unit: kW;
represent accumulator t period discharge power, unit: kW; K
om, btrepresent accumulator operation and maintenance cost, unit: $/kWh;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat release power of heat storage tank t period, unit: kW; K
om, tstrepresent heat storage tank operation and maintenance cost, unit: $/kWh.
In formula,
the expression system t period from the price of major network power purchase, unit: $/kWh;
the expression system t period is to the price of major network sale of electricity; α be 0 or 1, α=1 represent that micro-capacitance sensor can to major network sale of electricity, α=0 represents and does not allow micro-capacitance sensor to major network sale of electricity;
expression system t period and the mutual power of major network, unit: kW,
represent from major network power purchase,
represent to major network sale of electricity.
Then determine constraint condition, comprise cold energy balance, thermal energy balance, electric energy balance and equipment and run constraint condition.
Cold energy equilibrium constraint is determined according to formula (5):
In formula,
represent the power input of adsorbent refrigerator t period, unit: kW; COP
acrepresent the coefficient of refrigerating performance of adsorbent refrigerator;
represent the power input of electric refrigerating machine t period, unit: kW; COP
ecrepresent the coefficient of refrigerating performance of electric refrigerating machine;
represent t cooling load of the air-conditioning system power, unit: kW.
Thermal energy equilibrium constraint is determined according to formula (6):
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; η
mtrepresent the efficiency of miniature gas turbine; η
lossrepresent the gas turbine energy proportion of goods damageds; η
hrrepresent waste-heat recoverer efficiency;
represent the power of gas fired-boiler t period, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat release power of heat storage tank t period, unit: kW;
represent the power input of adsorbent refrigerator t period, unit: kW;
represent the thermal load power of t period system, unit: kW; η
herepresent effectiveness of heat exchanger.
Electric energy balance constraint condition is determined according to formula (7):
In formula,
represent the predicted value of photovoltaic t period;
represent the predicted value of wind-powered electricity generation t period;
represent the electric power of miniature gas turbine t period, unit: kW;
expression system t period and the mutual power of major network, unit: kW;
represent the charge power of accumulator t period, unit: kW;
represent accumulator t period discharge power, unit: kW;
represent the power input of electric refrigerating machine t period, unit: kW;
represent t period system electric load power, unit: kW.
Determine that miniature gas turbine runs constraint condition according to formula (801) and formula (802), its Chinese style (801) is miniature gas turbine running status constraint condition; Formula (802) represents gas turbine unit Climing constant condition, comprises Unit Commitment Climing constant and runs Climing constant continuously:
In formula (801),
represent miniature gas turbine t period running status variable,
represent that miniature gas turbine runs,
represent that miniature gas turbine is shut down;
represent the lower limit that miniature gas turbine is exerted oneself, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW;
represent the upper limit that miniature gas turbine is exerted oneself, unit: kW;
represent that unit falls power, in maximum when continuous running status of micro-gas-turbine unit: kW;
represent miniature gas turbine t-1 period running status variable;
to represent when micro-gas-turbine unit is shut down maximum falls power, unit: kW;
represent the electric power of miniature gas turbine t-1 period, unit: kW;
represent most the increase power of micro-gas-turbine unit when continuous running status, unit: kW;
represent the power that increases most when micro-gas-turbine unit starts, unit: kW.
The constraint condition that gas fired-boiler runs is determined according to formula (9):
In formula,
represent the lower limit that gas fired-boiler is exerted oneself, unit: kW;
represent the power of gas fired-boiler t period, unit: kW;
represent the upper limit that gas fired-boiler is exerted oneself, unit: kW.
The constraint condition of the mutual power of electrical network is determined according to formula (10):
In formula,
the lower limit of expression system and the mutual power of major network, unit kW;
the upper limit of expression system and the mutual power of major network, unit: kW;
expression system t period and the mutual power of major network, unit: kW.
The constraint condition that accumulator runs is determined according to formula (11):
In formula,
represent the discharge condition of accumulator t period,
represent battery discharging;
represent that accumulator does not charge also not discharge;
represent charge in batteries power upper limit, unit: kW;
represent accumulator t period discharge power, unit: kW;
represent the charge power of accumulator t period, unit: kW;
represent the charged state of accumulator t period,
represent charge in batteries;
represent that accumulator does not charge also not discharge;
represent battery discharging power upper limit, unit: kW;
represent the energy of t+1 period in accumulator, unit: kWh;
represent the energy of t period in accumulator, unit: kWh; σ
btrepresent the self-energy proportion of goods damageds of accumulator; η
c, btrepresent the charge efficiency of accumulator; η
disc, btrepresent battery discharging efficiency;
represent the lower limit of accumulator storage power, unit: kWh;
represent the upper limit of accumulator storage power, unit: kWh; Δ t represents the time interval.
The constraint condition that heat storage tank runs is determined according to formula (12):
In formula,
represent the heat release state of heat storage tank t period,
represent heat storage tank heat release,
represent heat storage tank not heat release also not accumulation of heat;
represent the accumulation of heat power upper limit of heat storage tank, unit kW;
represent the heat release power of heat storage tank t period, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat storage state of heat storage tank t period,
represent heat storage tank accumulation of heat,
represent heat storage tank not heat release also not accumulation of heat;
represent the heat release power upper limit of heat storage tank, unit: kW;
represent the energy of t+1 period in heat storage tank, unit: kWh;
represent the energy of t period in heat storage tank, unit: kWh; σ
tstrepresent the self-energy proportion of goods damageds of heat storage tank; η
c, tstrepresent the heat storage efficiency of heat storage tank; η
disc, tstrepresent the efficiency of heat storage tank release heat;
represent the lower limit of heat storage tank storage power, unit: kWh;
represent the upper limit of heat storage tank storage power, unit: kWh.
In each moment, according to step 10) the up-to-date refrigeration duty that obtains, thermal load and electric load predict the outcome and miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, the running status that energy storage device is real-time, Yalmip optimization tool is adopted to solve rolling optimization model, measuring and calculating miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, energy storage device will control exerting oneself of time domain M period in future, by the miniature gas turbine of the first period in the following M period, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, each equipment is delivered in the instruction of exerting oneself of energy storage device.
Step 30) monitoring real time data, and upgrade historical data: monitoring is obtained the wind power in each moment, photovoltaic power, refrigeration duty power, thermal load power and electric load power actual value, replaced the wind power in a upper moment, photovoltaic power, refrigeration duty power, thermal load power and electric load power actual value.
Step 40) set up feedback compensation model, the actual of equipment such as real-time correction miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, energy storage device are exerted oneself: the actual value of Real-Time Monitoring wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power, and the historical data upgrading wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load and electric load; Feedback compensation model is such as formula shown in (13) to formula (21).
Set up such as formula shown in (13) with the minimum correction function for objective function of relative adjustment amount:
In formula, AD represents overall relative adjustment amount; △ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
represent the rated power that miniature gas turbine exports, unit: kW; △ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW;
represent the rated power of electric refrigerating machine machine input, unit: kW; △ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent the rated power of charge in batteries, unit: kW; △ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the rated power of battery discharging, unit: kW; △ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
expression system and the mutual rated power of major network, unit: kW; △ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent that gas fired-boiler exports rated power, unit: kW; △ P
acrepresent the adjustment amount of adsorbent refrigerator power input, unit: kW;
represent the rated power of adsorbent refrigerator input, unit: kW; △ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the rated power of heat storage tank accumulation of heat, unit: kW; △ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the rated power of heat storage tank heat release, unit: kW; w
1represent the weight coefficient relevant to electric energy, w
2represent the weight coefficient relevant to cold and hot energy.
Set up such as formula the cold energy equilibrium constraint shown in (14):
In formula,
represent the power input of adsorbent refrigerator t period, unit: kW; △ P
acrepresent the adjustment amount of adsorbent refrigerator power input; COP
acrepresent the coefficient of refrigerating performance of adsorbent refrigerator;
represent the power input of electric refrigerating machine t period, unit: kW; △ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW; COP
ecrepresent the coefficient of refrigerating performance of electric refrigerating machine; Q
cfor real-time cooling load of the air-conditioning system power, unit: kW.
Set up such as formula the thermal energy equilibrium constraint shown in (15):
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; △ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW; η
mtrepresent the efficiency of miniature gas turbine; η
lossrepresent the gas turbine energy proportion of goods damageds; η
hrrepresent waste-heat recoverer efficiency;
represent the power of gas fired-boiler t period, unit: kW; △ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kW; △ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the heat release power of heat storage tank t period, unit: kW; △ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the power input of adsorbent refrigerator t period, unit: kW; △ P
acrepresent the adjustment amount of adsorbent refrigerator power input; η
herepresent effectiveness of heat exchanger; Q
hfor real-time system heat load power, unit: kW.
Set up such as formula the electric energy balance constraint condition shown in (16):
In formula, P
pvrepresent real-time photovoltaic power, unit: kW; P
wtrepresent real-time wind power, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW; △ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
expression system t period and the mutual power of major network, unit: kW; △ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
represent the charge power of accumulator t period, unit: kW; △ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent accumulator t period discharge power, unit: kW; △ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the power input of electric refrigerating machine t period, unit: kW; △ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW; P
elrepresent real-time system electric load power, unit: kW.
Set up the constraint condition run such as formula the miniature gas turbine shown in (17):
In formula,
represent miniature gas turbine t period running status variable,
represent that miniature gas turbine runs,
represent that miniature gas turbine does not run;
represent the lower limit that miniature gas turbine is exerted oneself, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW; △ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
represent the upper limit that miniature gas turbine is exerted oneself, unit: kW;
represent that unit falls power, in maximum in short-term when continuous running status of micro-gas-turbine unit: kW;
represent in short-term most the increase power of micro-gas-turbine unit when continuous running status, unit: kW.
Set up the constraint condition run such as formula the gas fired-boiler shown in (18):
In formula,
represent the lower limit that gas fired-boiler is exerted oneself, unit: kW;
represent the power of gas fired-boiler t period, unit: kW; △ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent the upper limit that gas fired-boiler is exerted oneself, unit: kW.
Set up the constraint condition such as formula the mutual power of electrical network shown in (19):
formula (19)
In formula,
the lower limit of expression system and the mutual power of major network, unit: kW;
expression system t period and the mutual power of major network, unit: kW; △ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
the upper limit of expression system and the mutual power of major network, unit: kW.
Set up the constraint condition run such as formula the accumulator shown in (20):
In formula,
represent the discharge condition of accumulator t period,
represent battery discharging,
represent that accumulator does not charge also not discharge;
represent charge in batteries power upper limit, unit: kW;
represent accumulator t period discharge power, unit: kW; △ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the charge power of accumulator t period, unit: kW; △ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent the charged state of accumulator t period,
represent charge in batteries,
represent that accumulator does not charge also not discharge;
represent battery discharging power upper limit, unit: kW;
energy before expression feedback compensation in accumulator, unit: kWh;
represent the energy in the accumulator after feedback compensation, unit: kWh; σ
btrepresent the self-energy proportion of goods damageds of accumulator; η
c, btrepresent the charge efficiency of accumulator; η
disc, btrepresent battery discharging efficiency;
represent the lower limit of accumulator storage power, unit: kWh;
represent the upper limit of accumulator storage power, unit: kWh; Δ t '=1/12h.
Set up the constraint condition run such as formula the heat storage tank shown in (21):
In formula,
represent the heat release state of heat storage tank t period,
represent heat storage tank heat release,
represent heat storage tank not heat release also not accumulation of heat;
represent charge in batteries power upper limit, unit: kW;
represent the heat release power of heat storage tank t period, unit: kW; △ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh; △ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the heat storage state of heat storage tank t period,
represent heat storage tank accumulation of heat,
represent heat storage tank not heat release also not accumulation of heat;
represent the heat release power upper limit of heat storage tank, unit: kW;
energy before expression feedback compensation in heat storage tank, unit: kWh;
represent the energy in the heat storage tank after feedback compensation, unit: kWh; σ
tstrepresent the self-energy proportion of goods damageds of heat storage tank; η
c, tstrepresent the heat storage efficiency of heat storage tank; η
disc, tstrepresent the efficiency of heat storage tank release heat;
represent the lower limit of heat storage tank storage power, unit: kWh;
represent the upper limit of heat storage tank storage power, unit: kWh.
Finally, Yalmip optimization tool is adopted to solve feedback compensation model, obtain the adjustment amount that miniature gas turbine is exerted oneself, gas fired-boiler is exerted oneself, adsorbent refrigerator power input, electric refrigerating machine power input, accumulator cell charging and discharging power, heat storage tank store the mutual power of heat release power, system and electrical network, these adjustment amounts are issued to respectively miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, accumulator, heat storage tank equipment adjusts; Within every 5 minutes, perform a step 30) and step 40), until be finished in control cycle Δ t.The time interval of rolling optimization is Δ t=0.25h, and the time interval of feedback compensation is Δ t '=1/12h, and feedback compensation is the correction to rolling optimization result, will do 3 feedback compensations within the time interval of a rolling optimization.
Step 50) enter subsequent time, return step 10), until micro-capacitance sensor out of service.
Step 10 of the present invention) in, wind-powered electricity generation forecast model and photovoltaic forecast model adopt document " the wind energy turbine set short term power forecast model of Kalman filtering correction " (Zhao Pan, Dai Yiping, Xia Junrong etc. the wind energy turbine set short term power forecast model [J] of Kalman filtering correction. XI AN JIAOTONG UNIVERSITY Subject Index, 2011, (5): 47-51.) in forecast model.Cooling load prediction model, heat load prediction model and electric load forecast model all adopt document " application of Nonparametric Autoregressive method in short-term electric load prediction " (Zhao Yuan, Zhang Xiafei, Xie Kaigui. the application of Nonparametric Autoregressive method in short-term electric load prediction [J]. High-Voltage Technology, 2011, (2): 429-435.) in forecast model.
In step 20) in, at each moment k, according to step 10) the up-to-date refrigeration duty that obtains, thermal load and electric load predict the outcome and miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, the running status that energy storage device is real-time, Yalmip optimization tool is adopted to solve rolling optimization model, measuring and calculating miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, energy storage device will control exerting oneself of time domain M period in future, by the miniature gas turbine of first period in the following M period, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, each equipment is delivered in the instruction of exerting oneself of energy storage device.Predict the outcome according to up-to-date to subsequent time again, recalculate rolling optimization model, so circulation is carried out.This ensure that system call instruction is all calculating on the basis of up-to-date information.
Step 30) in, the actual value of Real-Time Monitoring wind power, photovoltaic power, refrigeration duty, thermal load and electric load, and the historical data upgrading wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load and electric load.The predicted value of wind power, photovoltaic power, refrigeration duty, thermal load and electric load and actual value are often different, need to pass through feedback compensation, the actual of equipment such as real-time correction miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, energy storage device are exerted oneself, feedback compensation is the small adjustment to rolling optimal dispatching instruction, in order to revise the error brought because ex ante forecasting is inaccurate, feedback compensation model is the Optimized model of a quiet hour section.Exerting oneself on predicted data basis of each equipment that rolling optimization obtains.When actual motion, predicted value is different with actual value.If each equipment runs according to the result of rolling optimization, just system cloud gray model cannot be met, so will revise exerting oneself of each equipment according to actual value, namely feedback compensation.
Fig. 2 is the Systematical control block diagram of the inventive method, X in figure
kthe predicted value of kth moment to wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load, the following M time domain of electric load; X is wind-powered electricity generation, photovoltaic, refrigeration duty, heat are carved, the Real-Time Monitoring value of electricity; U
k+1be that each facilities plan of first period of kth moment rolling optimization result is exerted oneself instruction, comprise Gas Turbine Output, gas fired-boiler is exerted oneself, adsorbent refrigerator, electric refrigerating machine power input, accumulator cell charging and discharging power, heat storage tank stores heat release power, system and the mutual power of electrical network etc.; △ u is the Gas Turbine Output after Real-time Feedback corrects, and gas fired-boiler is exerted oneself, adsorbent refrigerator, electric refrigerating machine power input, accumulator cell charging and discharging power, and heat storage tank stores heat release power, the adjustment amount of system and the mutual power of electrical network;
for each Gas Turbine Output, gas fired-boiler is exerted oneself, adsorbent refrigerator, electric refrigerating machine power input, accumulator cell charging and discharging power, and heat storage tank stores heat release power, system and real-time the exerting oneself of the mutual power of electrical network.
First method of the present invention sets up the corresponding forecast model of wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load and electric load, and according to up-to-date historical data, prediction is following controls wind power, photovoltaic power, cool and thermal power load power data in time domain M; Each moment, according to up-to-dately predicting the outcome and each equipment real-time running state of obtaining, Yalmip optimization tool is adopted to solve this rolling optimization model, calculate each equipment exerting oneself in the following M period, but only perform the instruction of first period, predict the outcome according to up-to-date again to subsequent time, recalculate rolling optimization model, so circulate; The actual value of Real-Time Monitoring wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load, electric load, and upgrade the historical data of wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load, electric load, within every 5 minutes, solve feedback compensation model, obtain the adjustment amount of each equipment, and be issued to each equipment respectively and adjust, until the next rolling optimization moment.The inventive method can eliminate the impact because forecasting inaccuracy causes system optimized operation well, ensure what system operating instructions always calculated in the optimization of up-to-date information, reduce the risk of system cloud gray model, improve the stability of system cloud gray model, improve the economy that supply of cooling, heating and electrical powers type micro-capacitance sensor runs simultaneously.
The above is only the preferred embodiment of the present invention; it is noted that for those skilled in the art, under the premise without departing from the principles of the invention; can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.
Claims (3)
1., based on a supply of cooling, heating and electrical powers type micro-capacitance sensor operation method for Model Predictive Control, it is characterized in that, this operation method comprises the following steps:
Step 10) set up parametric prediction model, comprise wind-powered electricity generation forecast model, photovoltaic forecast model, cooling load prediction model, heat load prediction model and electric load forecast model; According to each parametric prediction model, in the k moment, according to n the historical data gathered before the k moment, utilize each parametric prediction model, following wind power, photovoltaic power, refrigeration duty power, thermal load power and the electric load power controlled in time domain M of prediction; Historical data comprises wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power;
Step 20) according to step 10) control in time domain M wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power in future of obtaining, set up such as formula the rolling optimization model shown in (1) to formula (12):
First determine with micro-grid system operating cost minimum for objective function, shown in (1):
In formula, C represents system operation cost; K represents current time; M represents control time domain;
the fuel cost of expression system t period;
the operation expense of expression system t period;
represent t period and the mutual cost of electrical network;
Then determine constraint condition, comprise cold energy balance, thermal energy balance, electric energy balance and equipment and run constraint condition:
Cold energy equilibrium constraint is determined according to formula (5):
In formula,
represent the power input of adsorbent refrigerator t period, unit: kW; COP
acrepresent the coefficient of refrigerating performance of adsorbent refrigerator;
represent the power input of electric refrigerating machine t period, unit: kW; COP
ecrepresent the coefficient of refrigerating performance of electric refrigerating machine;
represent t cooling load of the air-conditioning system power, unit: kW;
Thermal energy equilibrium constraint is determined according to formula (6):
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; η
mtrepresent the efficiency of miniature gas turbine; η
lossrepresent the gas turbine energy proportion of goods damageds; η
hrrepresent waste-heat recoverer efficiency;
represent the power of gas fired-boiler t period, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat release power of heat storage tank t period, unit: kW;
represent the power input of adsorbent refrigerator t period, unit: kW;
represent the thermal load power of t period system, unit: kW; η
herepresent effectiveness of heat exchanger;
Electric energy balance constraint condition is determined according to formula (7):
In formula,
represent the predicted value of photovoltaic t period;
represent the predicted value of wind-powered electricity generation t period;
represent the electric power of miniature gas turbine t period, unit: kW;
expression system t period and the mutual power of major network, unit: kW;
represent the charge power of accumulator t period, unit: kW;
represent accumulator t period discharge power, unit: kW;
represent the power input of electric refrigerating machine t period, unit: kW;
represent t period system electric load power, unit: kW;
Determine that miniature gas turbine runs constraint condition according to formula (801) and formula (802), its Chinese style (801) is miniature gas turbine running status constraint condition; Formula (802) represents gas turbine unit Climing constant condition, comprises Unit Commitment Climing constant and runs Climing constant continuously:
In formula (801),
represent miniature gas turbine t period running status variable,
represent that miniature gas turbine runs,
represent that miniature gas turbine is shut down;
represent the lower limit that miniature gas turbine is exerted oneself, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW;
represent the upper limit that miniature gas turbine is exerted oneself, unit: kW;
represent that unit falls power, in maximum when continuous running status of micro-gas-turbine unit: kW;
represent miniature gas turbine t-1 period running status variable;
to represent when micro-gas-turbine unit is shut down maximum falls power, unit: kW;
represent the electric power of miniature gas turbine t-1 period, unit: kW;
represent most the increase power of micro-gas-turbine unit when continuous running status, unit: kW;
represent the power that increases most when micro-gas-turbine unit starts, unit: kW;
The constraint condition that gas fired-boiler runs is determined according to formula (9):
In formula,
represent the lower limit that gas fired-boiler is exerted oneself, unit: kW;
represent the power of gas fired-boiler t period, unit: kW;
represent the upper limit that gas fired-boiler is exerted oneself, unit: kW;
The constraint condition of the mutual power of electrical network is determined according to formula (10):
In formula,
the lower limit of expression system and the mutual power of major network, unit kW;
the upper limit of expression system and the mutual power of major network, unit: kW;
expression system t period and the mutual power of major network, unit: kW;
The constraint condition that accumulator runs is determined according to formula (11):
In formula,
represent the discharge condition of accumulator t period,
represent battery discharging;
represent that accumulator does not charge also not discharge;
represent charge in batteries power upper limit, unit: kW;
represent accumulator t period discharge power, unit: kW;
represent the charge power of accumulator t period, unit: kW;
represent the charged state of accumulator t period,
represent charge in batteries;
represent that accumulator does not charge also not discharge;
represent battery discharging power upper limit, unit: kW;
represent the energy of t+1 period in accumulator, unit: kWh;
represent the energy of t period in accumulator, unit: kWh; σ
btrepresent the self-energy proportion of goods damageds of accumulator; η
c, btrepresent the charge efficiency of accumulator; η
disc, btrepresent battery discharging efficiency;
represent the lower limit of accumulator storage power, unit: kWh;
represent the upper limit of accumulator storage power, unit: kWh; Δ t represents the time interval;
The constraint condition that heat storage tank runs is determined according to formula (12):
In formula,
represent the heat release state of heat storage tank t period,
represent heat storage tank heat release,
represent heat storage tank not heat release also not accumulation of heat;
represent the accumulation of heat power upper limit of heat storage tank, unit kW;
represent the heat release power of heat storage tank t period, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat storage state of heat storage tank t period,
represent heat storage tank accumulation of heat,
represent heat storage tank not heat release also not accumulation of heat;
represent the heat release power upper limit of heat storage tank, unit: kW;
represent the energy of t+1 period in heat storage tank, unit: kWh;
represent the energy of t period in heat storage tank, unit: kWh; σ
tstrepresent the self-energy proportion of goods damageds of heat storage tank; η
c, tstrepresent the heat storage efficiency of heat storage tank; η
disc, tstrepresent the efficiency of heat storage tank release heat;
represent the lower limit of heat storage tank storage power, unit: kWh;
represent the upper limit of heat storage tank storage power, unit: kWh;
In each moment, according to step 10) the up-to-date refrigeration duty that obtains, thermal load and electric load predict the outcome and miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, the running status that energy storage device is real-time, Yalmip optimization tool is adopted to solve rolling optimization model, measuring and calculating miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, energy storage device will control exerting oneself of time domain M period in future, by the miniature gas turbine of first period in the following M period, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, each equipment is delivered in the instruction of exerting oneself of energy storage device,
Step 30) monitoring real time data, and upgrade historical data: monitoring is obtained the wind power in each moment, photovoltaic power, refrigeration duty power, thermal load power and electric load power actual value, replaced the wind power in a upper moment, photovoltaic power, refrigeration duty power, thermal load power and electric load power actual value;
Step 40) set up feedback compensation model, the actual of equipment such as real-time correction miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, energy storage device are exerted oneself: the actual value of Real-Time Monitoring wind power, photovoltaic power, refrigeration duty power, thermal load power and electric load power, and the historical data upgrading wind-powered electricity generation, photovoltaic, refrigeration duty, thermal load and electric load; Feedback compensation model is such as formula shown in (13) to formula (21):
Set up such as formula shown in (13) with the minimum correction function for objective function of relative adjustment amount:
In formula, AD represents overall relative adjustment amount; Δ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
represent the rated power that miniature gas turbine exports, unit: kW; Δ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW;
represent the rated power of electric refrigerating machine machine input, unit: kW; Δ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent the rated power of charge in batteries, unit: kW; Δ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the rated power of battery discharging, unit: kW; Δ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
expression system and the mutual rated power of major network, unit: kW; Δ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent that gas fired-boiler exports rated power, unit: kW; Δ P
acrepresent the adjustment amount of adsorbent refrigerator power input, unit: kW;
represent the rated power of adsorbent refrigerator input, unit: kW; Δ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the rated power of heat storage tank accumulation of heat, unit: kW; Δ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the rated power of heat storage tank heat release, unit: kW; w
1represent the weight coefficient relevant to electric energy, w
2represent the weight coefficient relevant to cold and hot energy;
Set up such as formula the cold energy equilibrium constraint shown in (14):
In formula,
represent the power input of adsorbent refrigerator t period, unit: kW; Δ P
acrepresent the adjustment amount of adsorbent refrigerator power input; COP
acrepresent the coefficient of refrigerating performance of adsorbent refrigerator;
represent the power input of electric refrigerating machine t period, unit: kW; Δ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW; COP
ecrepresent the coefficient of refrigerating performance of electric refrigerating machine; Q
cfor real-time cooling load of the air-conditioning system power, unit: kW;
Set up such as formula the thermal energy equilibrium constraint shown in (15):
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; Δ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW; η
mtrepresent the efficiency of miniature gas turbine; η
lossrepresent the gas turbine energy proportion of goods damageds; η
hrrepresent waste-heat recoverer efficiency;
represent the power of gas fired-boiler t period, unit: kW; Δ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kW; Δ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the heat release power of heat storage tank t period, unit: kW; Δ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the power input of adsorbent refrigerator t period, unit: kW; Δ P
acrepresent the adjustment amount of adsorbent refrigerator power input; η
herepresent effectiveness of heat exchanger; Q
hfor real-time system heat load power, unit: kW;
Set up such as formula the electric energy balance constraint condition shown in (16):
In formula, P
pvrepresent real-time photovoltaic power, unit: kW; P
wtrepresent real-time wind power, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW; Δ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
expression system t period and the mutual power of major network, unit: kW; Δ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
represent the charge power of accumulator t period, unit: kW; Δ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent accumulator t period discharge power, unit: kW; Δ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the power input of electric refrigerating machine t period, unit: kW; Δ P
ecrepresent the adjustment amount of electric refrigerating machine machine power input, unit: kW; P
elrepresent real-time system electric load power, unit: kW;
Set up the constraint condition run such as formula the miniature gas turbine shown in (17):
In formula,
represent miniature gas turbine t period running status variable,
represent that miniature gas turbine runs,
represent that miniature gas turbine does not run;
represent the lower limit that miniature gas turbine is exerted oneself, unit: kW;
represent the electric power of miniature gas turbine t period, unit: kW; Δ P
mtrepresent the adjustment amount that miniature gas turbine is exerted oneself, unit: kW;
represent the upper limit that miniature gas turbine is exerted oneself, unit: kW;
represent that unit falls power, in maximum in short-term when continuous running status of micro-gas-turbine unit: kW;
represent in short-term most the increase power of micro-gas-turbine unit when continuous running status, unit: kW;
Set up the constraint condition run such as formula the gas fired-boiler shown in (18):
In formula,
represent the lower limit that gas fired-boiler is exerted oneself, unit: kW;
represent the power of gas fired-boiler t period, unit: kW; Δ P
brepresent the adjustment amount that gas fired-boiler is exerted oneself, unit: kW;
represent the upper limit that gas fired-boiler is exerted oneself, unit: kW;
Set up the constraint condition such as formula the mutual power of electrical network shown in (19):
formula (19)
In formula,
the lower limit of expression system and the mutual power of major network, unit: kW;
expression system t period and the mutual power of major network, unit: kW; Δ P
gthe adjustment amount of expression system and the mutual power of major network, unit: kW;
the upper limit of expression system and the mutual power of major network, unit: kW;
Set up the constraint condition run such as formula the accumulator shown in (20):
In formula,
represent the discharge condition of accumulator t period,
represent battery discharging,
represent that accumulator does not charge also not discharge;
represent charge in batteries power upper limit, unit: kW;
represent accumulator t period discharge power, unit: kW; Δ P
disc, btrepresent the adjustment amount of battery discharging power, unit: kW;
represent the charge power of accumulator t period, unit: kW; Δ P
c, btrepresent the adjustment amount of charge in batteries power, unit: kW;
represent the charged state of accumulator t period,
represent charge in batteries,
represent that accumulator does not charge also not discharge;
represent battery discharging power upper limit, unit: kW;
energy before expression feedback compensation in accumulator, unit: kWh;
represent the energy in the accumulator after feedback compensation, unit: kWh; σ
btrepresent the self-energy proportion of goods damageds of accumulator; η
c, btrepresent the charge efficiency of accumulator; η
disc, btrepresent battery discharging efficiency;
represent the lower limit of accumulator storage power, unit: kWh;
represent the upper limit of accumulator storage power, unit: kWh; Δ t '=1/12h;
Set up the constraint condition run such as formula the heat storage tank shown in (21):
In formula,
represent the heat release state of heat storage tank t period,
represent heat storage tank heat release,
represent heat storage tank not heat release also not accumulation of heat;
represent charge in batteries power upper limit, unit: kW;
represent the heat release power of heat storage tank t period, unit: kW; Δ P
disc, tstrepresent the adjustment amount of heat storage tank heat release power, unit: kW;
represent the accumulation of heat power of heat storage tank t period, unit: kWh; Δ P
c, tstrepresent the adjustment amount of heat storage tank accumulation of heat power, unit: kW;
represent the heat storage state of heat storage tank t period,
represent heat storage tank accumulation of heat,
represent heat storage tank not heat release also not accumulation of heat;
represent the heat release power upper limit of heat storage tank, unit: kW;
energy before expression feedback compensation in heat storage tank, unit: kWh;
represent the energy in the heat storage tank after feedback compensation, unit: kWh; σ
tstrepresent the self-energy proportion of goods damageds of heat storage tank; η
c, tstrepresent the heat storage efficiency of heat storage tank; η
disc, tstrepresent the efficiency of heat storage tank release heat;
represent the lower limit of heat storage tank storage power, unit: kWh;
represent the upper limit of heat storage tank storage power, unit: kWh;
Finally, Yalmip optimization tool is adopted to solve feedback compensation model, obtain the adjustment amount that miniature gas turbine is exerted oneself, gas fired-boiler is exerted oneself, adsorbent refrigerator power input, electric refrigerating machine power input, accumulator cell charging and discharging power, heat storage tank store the mutual power of heat release power, system and electrical network, these adjustment amounts are issued to respectively miniature gas turbine, gas fired-boiler, adsorbent refrigerator, electric refrigerating machine, accumulator, heat storage tank equipment adjusts; Within every 5 minutes, perform a step 30) and step 40), until be finished in control cycle Δ t;
Step 50) enter subsequent time, return step 10), until micro-capacitance sensor out of service.
2., according to the supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control according to claim 1, it is characterized in that, described step 20) in,
In formula, Δ t represents the time interval; R
ngrepresent Gas Prices, unit: $/m
3;
represent the electric power of miniature gas turbine t period, unit: kW; η
mtrepresent the efficiency of miniature gas turbine; P
b trepresent the power of gas fired-boiler t period, unit: kW; η
brepresent the efficiency of gas fired-boiler; H
ngrepresent heating value of natural gas;
In formula,
represent the electric power of miniature gas turbine t period, unit: kW; K
om, mtrepresent miniature gas turbine operation and maintenance cost, unit: $/kWh;
represent the power of gas fired-boiler t period, unit: kW; K
om, brepresent gas fired-boiler operation and maintenance cost, unit: $/kWh;
represent the thermal load power of t period system, unit: kW; η
herepresent effectiveness of heat exchanger; K
om, herepresent heat exchanger operation and maintenance cost, unit: $/kWh;
represent the power input of adsorbent refrigerator t period, unit: kW; K
om, acrepresent adsorbent refrigerator operation and maintenance cost, unit: $/kWh;
represent the power input of electric refrigerating machine t period, unit: kW; K
om, ecrepresent electric refrigerating machine operation and maintenance cost, unit: $/kWh;
represent the predicted value of photovoltaic t period; K
om, pvrepresent photovoltaic cell maintenance cost unit: $/kWh;
represent the predicted value of wind-powered electricity generation t period; K
om, wtrepresent blower fan maintenance cost unit: $/kWh;
represent the charge power of accumulator t period, unit: kW;
represent accumulator t period discharge power, unit: kW; K
om, btrepresent accumulator operation and maintenance cost, unit: $/kWh;
represent the accumulation of heat power of heat storage tank t period, unit: kWh;
represent the heat release power of heat storage tank t period, unit: kW; K
om, tstrepresent heat storage tank operation and maintenance cost, unit: $/kWh;
In formula,
the expression system t period from the price of major network power purchase, unit: $/kWh;
the expression system t period is to the price of major network sale of electricity; α be 0 or 1, α=1 represent that micro-capacitance sensor can to major network sale of electricity, α=0 represents and does not allow micro-capacitance sensor to major network sale of electricity;
expression system t period and the mutual power of major network, unit: kW,
represent from major network power purchase,
represent to major network sale of electricity.
3., according to the supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control according to claim 2, it is characterized in that, described step 20) in, Δ t=0.25h, H
ng=9.78kWh/m
3.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510059153.8A CN104616208B (en) | 2015-02-04 | 2015-02-04 | A kind of supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510059153.8A CN104616208B (en) | 2015-02-04 | 2015-02-04 | A kind of supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104616208A true CN104616208A (en) | 2015-05-13 |
CN104616208B CN104616208B (en) | 2017-10-13 |
Family
ID=53150641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510059153.8A Expired - Fee Related CN104616208B (en) | 2015-02-04 | 2015-02-04 | A kind of supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104616208B (en) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105305472A (en) * | 2015-09-25 | 2016-02-03 | 清华大学 | Transformer substation capacity optimization method based on multi-energy collaborative energy supply |
CN105356521A (en) * | 2015-12-14 | 2016-02-24 | 东南大学 | AC and Dc mixed micro-grid operation optimization method based on time-domain rolling control |
CN105375479A (en) * | 2015-12-14 | 2016-03-02 | 东南大学 | Model predicative control based energy management method of distributed energy resource system |
CN105374126A (en) * | 2015-11-09 | 2016-03-02 | 珠海格力电器股份有限公司 | Distribution cabinet and control method and device for same |
CN105391061A (en) * | 2015-12-01 | 2016-03-09 | 贵州电网有限责任公司 | Active power distribution network optimization system and optimization method of cooling, heating and power wind-photovoltaic storage joint optimization |
CN105676824A (en) * | 2016-03-02 | 2016-06-15 | 山东大学 | Optimized energy dispatching system and method for renewable-energy-source-based combined supply of cooling, heating and power |
CN105869075A (en) * | 2016-04-19 | 2016-08-17 | 东南大学 | Economic optimization scheduling method for cold, heat and electricity combined supply type miniature energy grid |
CN105930986A (en) * | 2016-06-05 | 2016-09-07 | 朱明增 | High-voltage power grid overload accident prevention diagnosis processing system and method |
CN106410855A (en) * | 2016-11-29 | 2017-02-15 | 国网安徽省电力公司芜湖供电公司 | Independent micro energy network and optimal configuration method thereof |
CN106505634A (en) * | 2016-12-14 | 2017-03-15 | 东南大学 | Based on two benches coordination optimization and the supply of cooling, heating and electrical powers type microgrid operation method for controlling |
CN106647268A (en) * | 2016-12-21 | 2017-05-10 | 东南大学 | MGT-CCHP layered optimal control system based on model prediction control |
CN106786793A (en) * | 2016-12-14 | 2017-05-31 | 东南大学 | A kind of supply of cooling, heating and electrical powers type microgrid operation method based on robust optimization |
CN107703745A (en) * | 2017-09-21 | 2018-02-16 | 东南大学 | MGT CCHP control systems based on economic forecasting control |
CN108008629A (en) * | 2016-10-29 | 2018-05-08 | 南京理工大学 | The complementary optimizing operation method for utilizing system of one kind of multiple energy |
CN108596442A (en) * | 2018-03-30 | 2018-09-28 | 天津大学 | The integrated energy system economic load dispatching method of meter and Conditional Lyapunov ExponentP |
CN109286208A (en) * | 2018-10-09 | 2019-01-29 | 中国电力科学研究院有限公司 | A kind of integrated energy system dispatching method and system |
CN109726906A (en) * | 2018-12-21 | 2019-05-07 | 国网浙江省电力有限公司电力科学研究院 | Co-generation unit dispatching method a few days ago based on the constraint of heat supply network partial differential equation |
CN109768567A (en) * | 2018-12-20 | 2019-05-17 | 清华大学 | A kind of Optimization Scheduling coupling multi-energy complementation system |
CN110837916A (en) * | 2019-09-30 | 2020-02-25 | 国创新能源汽车能源与信息创新中心(江苏)有限公司 | Energy rolling optimization management method applied to home villa scene |
CN111930055A (en) * | 2020-09-29 | 2020-11-13 | 国网(天津)综合能源服务有限公司 | Comprehensive energy sensing device with optimized control |
CN111953018A (en) * | 2020-08-06 | 2020-11-17 | 许继集团有限公司 | Distributed multi-energy complementary energy supply system and power distribution method |
CN113356952A (en) * | 2021-06-25 | 2021-09-07 | 西安交通大学 | Combined cooling and power system capable of pre-cooling air at inlet of gas turbine and operation method thereof |
CN113554211A (en) * | 2021-05-27 | 2021-10-26 | 浙江工业大学 | ALIF-GRU-MTL source-load joint prediction-based comprehensive energy system optimization scheduling method |
CN113757924A (en) * | 2021-09-09 | 2021-12-07 | 广州珠江新城能源有限公司 | Refrigerating capacity prediction method and control system for centralized cooling system |
RU2783369C2 (en) * | 2021-02-26 | 2022-11-11 | Общество С Ограниченной Ответственностью «Русэнергопроект» | Method for optimization of energy consumption and greenhouse gas emission for technological object and system for its implementation |
CN115693793A (en) * | 2022-10-11 | 2023-02-03 | 国网浙江省电力有限公司 | Energy optimization control method for regional micro-grid |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102494430A (en) * | 2011-10-23 | 2012-06-13 | 西安交通大学 | Cold-electricity cogeneration system comprising wind power and gas combined cycle unit and method for scheduling cold-electricity cogeneration system |
CN103065197A (en) * | 2012-12-12 | 2013-04-24 | 中国能源建设集团广东省电力设计研究院 | Optimal configuration method of distributed combined cooling heating and power system |
US20130166043A1 (en) * | 2011-07-20 | 2013-06-27 | Nec Laboratories America, Inc. | Optimal energy management of a microgrid system using multi-objective optimization |
CN103455850A (en) * | 2013-08-07 | 2013-12-18 | 东南大学 | Online optimization method of grid-connected operation of distributed cool-heat-electricity cogeneration system |
CN103617460A (en) * | 2013-12-06 | 2014-03-05 | 天津大学 | Double-layer optimization planning and designing method for combined cooling, heating and power micro-grid system |
-
2015
- 2015-02-04 CN CN201510059153.8A patent/CN104616208B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130166043A1 (en) * | 2011-07-20 | 2013-06-27 | Nec Laboratories America, Inc. | Optimal energy management of a microgrid system using multi-objective optimization |
CN102494430A (en) * | 2011-10-23 | 2012-06-13 | 西安交通大学 | Cold-electricity cogeneration system comprising wind power and gas combined cycle unit and method for scheduling cold-electricity cogeneration system |
CN103065197A (en) * | 2012-12-12 | 2013-04-24 | 中国能源建设集团广东省电力设计研究院 | Optimal configuration method of distributed combined cooling heating and power system |
CN103455850A (en) * | 2013-08-07 | 2013-12-18 | 东南大学 | Online optimization method of grid-connected operation of distributed cool-heat-electricity cogeneration system |
CN103617460A (en) * | 2013-12-06 | 2014-03-05 | 天津大学 | Double-layer optimization planning and designing method for combined cooling, heating and power micro-grid system |
Non-Patent Citations (2)
Title |
---|
吴雄 等: "含热电联供系统的微网经济运行", 《电力自动化设备》 * |
王锐 等: "含可再生能源的热电联供型微网经济运行优化", 《电力系统自动化》 * |
Cited By (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105305472B (en) * | 2015-09-25 | 2018-10-23 | 清华大学 | A kind of substation capacity optimization method based on multiple-energy-source collaboration energy supply |
CN105305472A (en) * | 2015-09-25 | 2016-02-03 | 清华大学 | Transformer substation capacity optimization method based on multi-energy collaborative energy supply |
CN105374126A (en) * | 2015-11-09 | 2016-03-02 | 珠海格力电器股份有限公司 | Distribution cabinet and control method and device for same |
CN105374126B (en) * | 2015-11-09 | 2017-12-19 | 珠海格力电器股份有限公司 | Dispense cabinet and its control method and device |
CN105391061A (en) * | 2015-12-01 | 2016-03-09 | 贵州电网有限责任公司 | Active power distribution network optimization system and optimization method of cooling, heating and power wind-photovoltaic storage joint optimization |
CN105375479B (en) * | 2015-12-14 | 2017-11-17 | 东南大学 | A kind of distributed energy energy management method based on Model Predictive Control |
CN105356521A (en) * | 2015-12-14 | 2016-02-24 | 东南大学 | AC and Dc mixed micro-grid operation optimization method based on time-domain rolling control |
CN105375479A (en) * | 2015-12-14 | 2016-03-02 | 东南大学 | Model predicative control based energy management method of distributed energy resource system |
CN105356521B (en) * | 2015-12-14 | 2018-02-23 | 东南大学 | A kind of alternating current-direct current mixing micro-capacitance sensor running optimizatin method that control is rolled based on time domain |
CN105676824A (en) * | 2016-03-02 | 2016-06-15 | 山东大学 | Optimized energy dispatching system and method for renewable-energy-source-based combined supply of cooling, heating and power |
CN105676824B (en) * | 2016-03-02 | 2018-05-01 | 山东大学 | A kind of energy-optimised scheduling System and method for of regenerative resource supply of cooling, heating and electrical powers |
CN105869075A (en) * | 2016-04-19 | 2016-08-17 | 东南大学 | Economic optimization scheduling method for cold, heat and electricity combined supply type miniature energy grid |
CN105930986A (en) * | 2016-06-05 | 2016-09-07 | 朱明增 | High-voltage power grid overload accident prevention diagnosis processing system and method |
CN108008629A (en) * | 2016-10-29 | 2018-05-08 | 南京理工大学 | The complementary optimizing operation method for utilizing system of one kind of multiple energy |
CN106410855A (en) * | 2016-11-29 | 2017-02-15 | 国网安徽省电力公司芜湖供电公司 | Independent micro energy network and optimal configuration method thereof |
CN106786793A (en) * | 2016-12-14 | 2017-05-31 | 东南大学 | A kind of supply of cooling, heating and electrical powers type microgrid operation method based on robust optimization |
CN106505634A (en) * | 2016-12-14 | 2017-03-15 | 东南大学 | Based on two benches coordination optimization and the supply of cooling, heating and electrical powers type microgrid operation method for controlling |
CN106786793B (en) * | 2016-12-14 | 2019-04-09 | 东南大学 | A kind of supply of cooling, heating and electrical powers type microgrid operation method based on robust optimization |
CN106505634B (en) * | 2016-12-14 | 2019-03-19 | 东南大学 | The supply of cooling, heating and electrical powers type microgrid operation method with control is coordinated and optimized based on two stages |
CN106647268A (en) * | 2016-12-21 | 2017-05-10 | 东南大学 | MGT-CCHP layered optimal control system based on model prediction control |
CN107703745A (en) * | 2017-09-21 | 2018-02-16 | 东南大学 | MGT CCHP control systems based on economic forecasting control |
CN107703745B (en) * | 2017-09-21 | 2020-11-03 | 东南大学 | MGT-CCHP control system based on economic predictive control |
CN108596442A (en) * | 2018-03-30 | 2018-09-28 | 天津大学 | The integrated energy system economic load dispatching method of meter and Conditional Lyapunov ExponentP |
CN109286208A (en) * | 2018-10-09 | 2019-01-29 | 中国电力科学研究院有限公司 | A kind of integrated energy system dispatching method and system |
CN109768567A (en) * | 2018-12-20 | 2019-05-17 | 清华大学 | A kind of Optimization Scheduling coupling multi-energy complementation system |
CN109726906A (en) * | 2018-12-21 | 2019-05-07 | 国网浙江省电力有限公司电力科学研究院 | Co-generation unit dispatching method a few days ago based on the constraint of heat supply network partial differential equation |
CN110837916A (en) * | 2019-09-30 | 2020-02-25 | 国创新能源汽车能源与信息创新中心(江苏)有限公司 | Energy rolling optimization management method applied to home villa scene |
CN111953018A (en) * | 2020-08-06 | 2020-11-17 | 许继集团有限公司 | Distributed multi-energy complementary energy supply system and power distribution method |
CN111930055B (en) * | 2020-09-29 | 2021-01-15 | 国网(天津)综合能源服务有限公司 | Comprehensive energy sensing device with optimized control |
CN111930055A (en) * | 2020-09-29 | 2020-11-13 | 国网(天津)综合能源服务有限公司 | Comprehensive energy sensing device with optimized control |
RU2783369C2 (en) * | 2021-02-26 | 2022-11-11 | Общество С Ограниченной Ответственностью «Русэнергопроект» | Method for optimization of energy consumption and greenhouse gas emission for technological object and system for its implementation |
CN113554211A (en) * | 2021-05-27 | 2021-10-26 | 浙江工业大学 | ALIF-GRU-MTL source-load joint prediction-based comprehensive energy system optimization scheduling method |
CN113356952A (en) * | 2021-06-25 | 2021-09-07 | 西安交通大学 | Combined cooling and power system capable of pre-cooling air at inlet of gas turbine and operation method thereof |
CN113356952B (en) * | 2021-06-25 | 2022-12-09 | 西安交通大学 | Combined cooling and power system capable of pre-cooling air at inlet of gas turbine and operation method thereof |
CN113757924A (en) * | 2021-09-09 | 2021-12-07 | 广州珠江新城能源有限公司 | Refrigerating capacity prediction method and control system for centralized cooling system |
CN115693793A (en) * | 2022-10-11 | 2023-02-03 | 国网浙江省电力有限公司 | Energy optimization control method for regional micro-grid |
CN115693793B (en) * | 2022-10-11 | 2024-05-17 | 国网浙江省电力有限公司 | Regional micro-grid energy optimization control method |
Also Published As
Publication number | Publication date |
---|---|
CN104616208B (en) | 2017-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104616208A (en) | Model predication control based cooling heating and power generation type micro-grid operation method | |
CN106505634B (en) | The supply of cooling, heating and electrical powers type microgrid operation method with control is coordinated and optimized based on two stages | |
CN108154309B (en) | Energy internet economic dispatching method considering multi-load dynamic response of cold, heat and electricity | |
CN110417006A (en) | Consider the integrated energy system Multiple Time Scales energy dispatching method of multipotency collaboration optimization | |
CN111445090A (en) | Double-layer planning method for off-grid type comprehensive energy system | |
CN107807523A (en) | Consider the Regional Energy internet multi-source coordination optimization operation reserve of tou power price | |
CN109858759A (en) | A kind of industrial park comprehensive energy balance dispatching method | |
CN111027846B (en) | Electricity-hydrogen comprehensive energy system considering heat-hydrogen co-production and capacity configuration method thereof | |
CN106786793A (en) | A kind of supply of cooling, heating and electrical powers type microgrid operation method based on robust optimization | |
CN106950840A (en) | Towards the integrated energy system layered distribution type control method for coordinating of power network peak clipping | |
CN109409595B (en) | Garden multi-energy complementary system day-ahead scheduling method | |
CN111144620A (en) | Electricity-hydrogen comprehensive energy system considering seasonal hydrogen storage and robust planning method thereof | |
CN104951899A (en) | Multi-time-scale optimal scheduling method for power distribution company containing large-scale renewable energy sources | |
CN105225022A (en) | A kind of economy optimizing operation method of cogeneration of heat and power type micro-capacitance sensor | |
CN111737884B (en) | Multi-target random planning method for micro-energy network containing multiple clean energy sources | |
CN112600253B (en) | Park comprehensive energy collaborative optimization method and equipment based on optimal energy utilization efficiency | |
CN111160636B (en) | CCHP type micro-grid scheduling optimization method | |
CN112836882B (en) | Regional comprehensive energy system operation optimization method considering equipment load rate change | |
CN108596442A (en) | The integrated energy system economic load dispatching method of meter and Conditional Lyapunov ExponentP | |
CN112446546B (en) | Comprehensive energy system two-stage optimal configuration method considering energy reliability | |
CN110311371A (en) | A kind of photovoltaic refrigeration storage system and its load active control method based on virtual energy storage | |
CN108008629A (en) | The complementary optimizing operation method for utilizing system of one kind of multiple energy | |
Zhao et al. | Optimization and analysis of an integrated energy system based on wind power utilization and on-site hydrogen refueling station | |
Luo et al. | A two-stage energy management strategy for CCHP microgrid considering house characteristics | |
Xu et al. | Optimal intraday rolling operation strategy of integrated energy system with multi-storage |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171013 |
|
CF01 | Termination of patent right due to non-payment of annual fee |