CN113128070A - Optimal configuration method for comprehensive energy system of intermittent distributed power supply - Google Patents
Optimal configuration method for comprehensive energy system of intermittent distributed power supply Download PDFInfo
- Publication number
- CN113128070A CN113128070A CN202110498724.3A CN202110498724A CN113128070A CN 113128070 A CN113128070 A CN 113128070A CN 202110498724 A CN202110498724 A CN 202110498724A CN 113128070 A CN113128070 A CN 113128070A
- Authority
- CN
- China
- Prior art keywords
- power
- energy
- heat
- gas
- rated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 14
- 238000013178 mathematical model Methods 0.000 claims abstract description 11
- 239000007789 gas Substances 0.000 claims description 58
- 238000002485 combustion reaction Methods 0.000 claims description 44
- 239000002918 waste heat Substances 0.000 claims description 39
- 230000005611 electricity Effects 0.000 claims description 32
- 238000010438 heat treatment Methods 0.000 claims description 28
- 238000010248 power generation Methods 0.000 claims description 26
- 238000010521 absorption reaction Methods 0.000 claims description 24
- 238000001816 cooling Methods 0.000 claims description 24
- 238000005057 refrigeration Methods 0.000 claims description 17
- 238000012423 maintenance Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 10
- 230000009194 climbing Effects 0.000 claims description 8
- 238000007906 compression Methods 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 230000007774 longterm Effects 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 239000003546 flue gas Substances 0.000 claims description 3
- 210000001061 forehead Anatomy 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 239000003345 natural gas Substances 0.000 claims description 3
- 239000004065 semiconductor Substances 0.000 claims description 3
- 230000003993 interaction Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 6
- 230000006835 compression Effects 0.000 description 4
- 238000005286 illumination Methods 0.000 description 4
- 239000003507 refrigerant Substances 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 238000009834 vaporization Methods 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Abstract
The invention discloses an optimal configuration method of an intermittent distributed power supply comprehensive energy system, which comprises the following three steps: (1) analyzing the flow directions and interaction modes of cold energy flow, heat energy flow, electric energy flow and gas flow of the CCHP system, and providing an optimal configuration method of the comprehensive energy system; (2) intermittent distributed energy sources such as wind, light and the like are considered, a mathematical model with uncertainty and randomness is constructed, and the mathematical model and other energy sources are subjected to multi-energy complementation, so that cascade utilization is realized, and the energy utilization rate is improved; (3) and (3) taking the lowest total investment cost of the life whole period as an objective function, optimizing and solving to obtain the optimal capacity, and comparing the optimal capacity with the configuration of the traditional comprehensive energy system to verify the superiority of the method. The invention considers wind, light and other intermittent distributed energy sources, constructs a mathematical model with uncertainty and randomness, performs multi-energy complementation with other energy sources, realizes cascade utilization, improves the energy utilization rate, takes the lowest total investment cost of the life whole period as a target function, and optimizes the solution to obtain the optimal capacity.
Description
Technical Field
The invention relates to an optimal configuration method of an integrated energy system of an intermittent distributed power supply, and belongs to the technical field of energy.
Background
People's life and production are closely related to energy, and a comprehensive energy system which considers the energy coupling characteristics, can practically improve the energy utilization rate and effectively utilize various renewable resources gradually becomes the core for solving the energy problem. A combined cooling, heating and power (CCHP) system based on distributed energy can achieve diversified targets such as high energy utilization rate, low energy cost and environmental friendliness according to the cascade utilization of energy, and becomes an important direction and a main form of the development of a regional comprehensive energy system. In order to effectively improve the energy utilization rate and promote sustainable development, a combined cooling heating and power type comprehensive energy system considering intermittent distributed power sources can be constructed to perform coordinated optimization control so as to obtain better economic benefit.
Disclosure of Invention
The invention aims to provide an optimized configuration method of an integrated energy system of an intermittent distributed power supply, which is used for analyzing the flowing directions and interaction modes of cold energy flow, heat energy flow, electric energy flow and gas flow of a CCHP system and providing the optimized configuration method of the integrated energy system; intermittent distributed energy sources such as wind, light and the like are considered, a mathematical model with uncertainty and randomness is constructed, and the mathematical model and other energy sources are subjected to multi-energy complementation, so that cascade utilization is realized, and the energy utilization rate is improved; and (4) taking the lowest total investment cost of the life whole period as an objective function, and optimizing and solving to obtain the optimal capacity.
The purpose of the invention is realized by the following technical scheme:
the optimal configuration method of the comprehensive energy system of the intermittent distributed power supply comprises the following steps:
firstly, modeling each unit as follows:
(1) intermittent renewable energy source
1) Photovoltaic power generation
According to the photovoltaic effect of the semiconductor material, the photovoltaic cell can directly convert solar energy into electric energy, the temperature coefficient epsilon of the photovoltaic cell is 0.2, and a mathematical model is as follows:
f(t)=I+ε(T-Trated) (2)
in the formula, PPVRepresenting the photovoltaic output power, rratedIndicating the nominal light intensity, PratedRepresenting the photovoltaic power rating, TratedRepresents the nominal temperature, and T represents the photovoltaic cell temperature;
2) wind power generation
The long-term distribution of wind speed follows the Weibull distribution, which is commonly represented by a three-constant model:
wherein k is a Weibull shape parameter, c is a Weibull scale parameter, and v is a wind speed;
the mathematical model of wind power generation is as follows:
wherein v isci、vco、vratedRespectively representing cut-in wind speed, cut-out wind speed, rated wind speed, PWIndicating fan output power, PratedRepresents a rated power;
(2) combined supply equipment
1) Gas turbine
Variable working condition operation model of gas turbine:
η=aW3+bW2+cW (5-1)
in the formula, a is 1.071, b is-2.816, c is 2.745, eta and W are the specific values of the power generation efficiency and the output power of the gas turbine to the rated value, namely per unit; the formula (5) is a function of 3 times, and piecewise linearization is needed in CPLEX calculation; in addition, the climbing and descending constraints of the gas turbine are also considered, namely the absolute value of the difference value of the output of the gas turbine between two load values before and after the load change in actual operation is less than or equal to 0.5 times of rated power;
2) internal combustion engine
The variable working condition operation model of the internal combustion engine is as shown in the formula (5-2).
η=aW3+bW2+cW (5-2)
In the formula, a is 1.071, b is-2.816, c is 2.745, eta and W are the specific values of the power generation efficiency and the output power of the internal combustion engine to the rated value, namely the per unit value;
(3) refrigeration device
1) Double-effect absorption refrigerator
The refrigeration COP in summer is 1.2, namely 1.2 units of waste heat can be refrigerated every time the refrigeration COP absorbs 1 unit, the heating COP in winter is 0.9, namely 0.9 unit of waste heat can be heated every time the refrigeration COP absorbs 1 unit, the COP is the conversion rate between the absorbed energy and the released heat, which is called the heating energy efficiency ratio for short, and the variable working condition characteristic is not considered;
2) electric refrigerator
The refrigerating COP in summer is 4, namely 4 units of refrigerating can be performed by absorbing 1 unit of electricity, and the variable working condition characteristic is not considered;
(4) auxiliary equipment
1) Gas boiler
When the combined supply equipment cannot meet the heat demand of a user, the supplementary combustion of a gas-fired boiler is needed, the balance of supply and demand of heat energy is realized, and the variable working condition characteristic is not considered;
2) waste heat boiler
The equipment for changing water into steam by utilizing the heat of the flue gas discharged by the gas turbine does not consider the variable working condition characteristic;
secondly, setting the target function and the constraint condition as follows:
(1) objective function
After a photovoltaic power generation system and a wind power generation system are established, a randomness model is introduced into a comprehensive energy system model, an internal combustion engine is used as a prime mover, the total investment cost of the life whole period is used as an objective function, the objective function comprises the construction cost, the operation cost and the maintenance cost of each device, and the optimal capacity of the main devices of the system is optimized and configured; the objective function is as follows:
C=∑CT+∑Com+∑Cup (6)
in the formula, CTFor construction costs, ComFor operating costs, CupFor maintenance cost, the specific expressions of each part are as follows:
1)CTincluding all equipment converted to annual installed costs:
in the formula, IntRateFor interest rate, t is the current year, RconverseTo convert annual costs to a factor of the current year, RrecoverFor depreciation rate, MjFor construction unit price, RM,jIs the capacity of the device;
2)Comincluding the gas cost and the electricity purchasing cost when the system generated energy is insufficient:
in the formula, Cg,iFor gas costs of internal combustion engines and gas boilers, Ce,IFor the purchase of electricity for the voltage-compression refrigerator and for the upper-level power grid, the peak-valley rate system, P, is usedbuyIs the power purchase, PECRIs the electric power of the electric compressor, VICE、VGFBVolume of gas consumed by internal combustion engines and gas boilers, KgIs the price of natural gas, KeIs the price of electricity, NdayHours to use;
3)Cupincluding maintenance costs for all equipment;
(2) constraint conditions
The comprehensive energy system needs to meet 4 supply and demand balances, namely 4 equality constraint conditions, namely, electric balance, heat balance, cold balance and waste heat balance; the electric balance, the heat balance and the cold balance refer to electric quantity, heat quantity and cold quantity which are generated by equipment and can meet the electric demand, the heat demand and the cold demand of a system; the waste heat balance means that waste heat generated by a waste heat boiler and a gas boiler needs to meet the waste heat requirement of the double-effect absorption refrigerator; inequality constraints of equipment operation, namely equipment power limitation and climbing constraint of the internal combustion engine, need to be met;
1) power balance constraint
In the formula, PPV、PW、PICERespectively the generated power of the photovoltaic, the fan and the internal combustion engine, PLoad_power、PLoad_heating、PLoad_coolingRespectively electric, thermal and cold load, PbuyIs the power purchase, PECRIs the electrical power of the electric compressor, PDAR_heating、PWHB、PGFBRespectively the heating power, P, of the double-effect absorption refrigerator, the exhaust-heat boiler and the gas boilerDAR_cooling、PECR_coolingRefrigeration powers of the double-effect absorption refrigerator and the electric compressor, respectively, EICEFor the residual heat power of internal combustion engines, EDAR、EWHBThe heat absorption power of the double-effect absorption refrigerator and the waste heat boiler is respectively;
2) plant operating constraints
The light intensity, wind speed, temperature and the like are greatly influenced by natural environment and have strong randomness, if the total capacity of the distributed power supplies is overlarge, the power supply quality and stability of the system are influenced, so that the installation capacity of the distributed power supplies is limited, and the constraint condition is that
SPV、SWRespectively the apparent power of the photovoltaic and the fan,respectively representing the apparent power limit values of the photovoltaic and the draught fan;
PPV、PW、PICE、Pbuy、PECR、PWHB、PGFB、PDAR_heating、PDAR_coolingthe upper and lower limits of the power of the respective devices are satisfied:
in addition, the internal combustion engine also satisfies the hill climbing constraint:
|Pi,j+1 ICE-Pi,j ICE|≤0.5 PICE, forehead (12)。
Compared with the prior art, the invention has the beneficial effects that: an optimized configuration model of a comprehensive energy system which comprehensively covers a fan, a photovoltaic system, an internal combustion engine, a boiler, refrigeration equipment, an electric load, a heat load, a cold load and the like is established, wind-light-electricity multi-energy complementation and cooling-heat-electricity multi-energy flow coordinated transfer in a region are comprehensively considered, and the cascade utilization of energy is promoted. And the lowest total investment cost of the life whole period is taken as an objective function, and the optimal capacity of the system is obtained by optimal solution.
Drawings
FIG. 1 is a diagram of an integrated energy system architecture;
FIG. 2 is a 1 month typical day heat load balance diagram;
FIG. 3 is a 1 month typical day electricity load balancing diagram;
FIG. 4 is a 7 month typical day cooling load balancing diagram;
FIG. 5 is a 7 month typical day heat load balance diagram;
fig. 6 is a 7-month typical daily electricity load balance diagram.
Detailed Description
The invention is further described with reference to the following figures and specific examples.
Fig. 1 is a structural view of a combined cooling heating and power type integrated energy system in which intermittent distributed energy is considered. The gas turbine and the internal combustion engine of the comprehensive energy system are used as prime movers to generate electricity and generate waste heat, renewable energy is used for generating electricity, the double-effect absorption type refrigerating machine can absorb the waste heat of a boiler to refrigerate or heat, the electric refrigerating machine absorbs electric power to refrigerate, the waste heat boiler absorbs the waste heat generated by the internal combustion engine to heat, the gas boiler consumes gas to heat, and the gas boiler is purchased from a higher-level power grid when the electric quantity is insufficient.
The invention discloses an optimal configuration method of an integrated energy system of an intermittent distributed power supply, wherein the integrated energy system of the intermittent distributed power supply is divided into intermittent renewable energy, combined supply equipment, refrigerating equipment and auxiliary equipment, and the optimal configuration method of the integrated energy system of the intermittent distributed power supply comprises the following steps:
firstly, modeling each unit as follows:
(1) intermittent renewable energy source
1) Photovoltaic power generation
According to the photovoltaic effect of the semiconductor material, the photovoltaic cell can directly convert solar energy into electric energy, the temperature coefficient epsilon of the photovoltaic cell is 0.2, and a mathematical model is as follows:
f(t)=1+ε(T-Trated) (2)
in the formula, PPVRepresenting the photovoltaic output power, rratedRepresenting rated light intensity, taking 1000W/m2,PratedRepresenting the rated power of photovoltaic, taking 100kW, TratedRepresenting rated temperature, taking 0 ℃ and T representing the temperature of the photovoltaic cell;
2) wind power generation
Wind in a natural state has randomness and intermittence, the size and the direction of wind power change every moment, and accurate prediction cannot be carried out, so that the wind power often fluctuates frequently. When a model of the distributed wind power system is established, wind power characteristics of the distributed wind power system under different time scales need to be considered. It is generally accepted that the long-term distribution of wind speed follows the Weibull distribution, which is commonly represented by a three-constant model:
wherein k is a Weibull shape parameter, c is a Weibull scale parameter, and v is a wind speed;
the mathematical model of wind power generation is as follows:
wherein v isci、vco、vratedThe cut-in wind speed, the cut-out wind speed and the rated wind speed are respectively expressed, and 3m/s, 20m/s and 10m/s are respectively adopted. PWIndicating fan output power, PratedRepresenting the rated power, 100kW is taken.
(2) Combined supply equipment
1) Gas turbine
The internal combustion engine utilizes continuously flowing gas to push an impeller to rotate rapidly, converts gas energy into mechanical energy, and has the main advantages of small size, light weight and low efficiency. The rated power generation efficiency of the gas turbine is 0.3, the available waste heat coefficient is 0.5, and the life cycle is 20 years. The following formula is a variable working condition operation model of the gas turbine:
η=aW3+bW2+cW (5-1)
in the formula, a is 1.071, b is-2.816, c is 2.745, eta and W are the specific values of the power generation efficiency and the output power of the gas turbine to the rated value, namely per unit; equation (5) is a 3-degree function, and piecewise linearization is required in the CPLEX calculation. In addition, the ramp-up and ramp-down constraints of the gas turbine are taken into account, i.e. the absolute value of the difference in the output of the gas turbine between the two load values before and after the load change during actual operation is less than or equal to 0.5 nominal power.
2) Internal combustion engine
It is an engine which directly converts the heat released by the combustion of fuel in the interior of the machine into kinetic energy. The internal combustion engine has the advantages of high thermal efficiency, compact structure, strong maneuverability, simple and convenient operation and maintenance and the like, the rated power generation efficiency of the internal combustion engine is 0.42, the available waste heat coefficient is 0.4, and the life cycle is 20 years. The variable working condition operation model of the internal combustion engine is as shown in the formula (5-2).
η=aW3+bW2+cW (5-2)
In the formula, a is 1.071, b is-2.816, c is 2.745, eta and W are the specific values of the power generation efficiency and the output power of the internal combustion engine to the rated value, namely the per unit value;
(3) refrigeration device
1) Double-effect absorption refrigerator
The principle is that after the waste heat of a gas boiler and a gas turbine is absorbed, a high-temperature liquid refrigerant is throttled and depressurized, and absorbs a large amount of heat during vaporization, so that refrigeration is performed. Refrigerating COP in summer is 1.2 (1.2 units can be refrigerated by absorbing 1 unit of waste heat); the winter heating COP is 0.9 (0.9 unit can be heated every 1 unit of waste heat is absorbed). Cop (coefficient of performance), which is the rate of conversion between absorbed energy and released heat, is called simply the heating energy efficiency ratio. Regardless of the behavior of the variable mode.
2) Electric refrigerator
The principle of the method is that an electrically driven compressor is used for discharging a refrigerant in a low-temperature and low-pressure gas state into a machine body for compression, the refrigerant is introduced into a condenser under the condition that the temperature and the pressure meet the requirement, the high-temperature and high-pressure refrigerant releases heat, and condensed water absorbs heat, so that refrigeration is realized. The summer refrigeration COP is 4 (4 units can be refrigerated per 1 unit of electricity absorbed), regardless of the variable-operating characteristics.
(4) Auxiliary equipment
1) Gas boiler
When the combined supply equipment can not meet the heat demand of a user, the supplementary combustion of a gas-fired boiler is needed, and the balance of supply and demand of heat energy is realized. The rated efficiency is 0.85, and the variable working condition characteristic is not considered.
2) Waste heat boiler
The equipment for changing water into steam by utilizing the heat of the flue gas discharged by the gas turbine has the rated efficiency of 0.9, and does not consider the variable working condition characteristic.
Secondly, setting the target function and the constraint condition as follows:
(1) objective function
After a photovoltaic power generation system and a wind power generation system are established, a randomness model is introduced into a comprehensive energy system model, an internal combustion engine is used as a prime mover, the total investment cost of the life whole period is used as an objective function, the objective function comprises the construction cost, the operation cost and the maintenance cost of each device, and the optimal capacity of the main devices of the system is optimized and configured; the objective function is as follows:
C=∑CT+∑Com+∑Cup (6)
in the formula, CTFor construction costs, ComFor operating costs, CupFor maintenance cost, the specific expressions of each part are as follows:
1)CTincluding all equipment converted to annual installed costs:
in the formula, IntRateFor interest rate, take 0.03, t is the current year, RconverseTo convert annual costs to a factor of the current year, RrecoverFor depreciation rate, MjFor construction unit price, RM,jIs the capacity of the device;
2)Comincluding the gas cost and the electricity purchasing cost when the system generated energy is insufficient:
in the formula, Cg,iFor gas costs of internal combustion engines and gas boilers, Ce,iFor the purchase of electricity for the voltage-compression refrigerator and for the upper-level power grid, the peak-valley rate system, P, is usedbuyIs the power purchase, PECRIs the electric power of the electric compressor, VICE、VGFBVolume of gas consumed by internal combustion engines and gas boilers, KgIs natural gas price, and takes 2.28 yuan/m3,KeIs the price of electricity, NdayHours to use;
3)Cupincluding maintenance costs for all equipment;
(2) constraint conditions
The comprehensive energy system needs to meet 4 supply and demand balances, namely 4 equality constraint conditions, namely, electric balance, heat balance, cold balance and waste heat balance; the electric balance, the heat balance and the cold balance refer to electric quantity, heat quantity and cold quantity which are generated by equipment and can meet the electric demand, the heat demand and the cold demand of a system; the waste heat balance means that waste heat generated by a waste heat boiler and a gas boiler needs to meet the waste heat requirement of the double-effect absorption refrigerator; besides, inequality constraints of equipment operation, namely equipment power limitation and climbing constraint of the internal combustion engine, need to be met;
2) power balance constraint
In the formula, PPV、PW、PICERespectively the generated power of the photovoltaic, the fan and the internal combustion engine, PLoad_power、PLoad_heating、PLoad_coolingRespectively electric, thermal and cold load, PbuyIs the power purchase, PECRIs the electrical power of the electric compressor, PDAR_heating、PWHB、PGFBRespectively the heating power, P, of the double-effect absorption refrigerator, the exhaust-heat boiler and the gas boilerDAR_cooling、PECR_coolingRefrigeration powers of the double-effect absorption refrigerator and the electric compressor, respectively, EICEFor the residual heat power of internal combustion engines, EDAR、EWHBThe heat absorption power of the double-effect absorption refrigerator and the waste heat boiler is respectively;
2) plant operating constraints
The light intensity, wind speed, temperature and the like are greatly influenced by natural environment and have strong randomness, if the total capacity of the distributed power supplies is overlarge, the power supply quality and stability of the system are influenced, so that the installation capacity of the distributed power supplies is limited, and the constraint condition is that
SPV、SWRespectively the apparent power of the photovoltaic and the fan,respectively representing the apparent power limit values of the photovoltaic and the draught fan;
PPV、PW、PICE、Pbuy、PECR、PWHB、PGFB、PDAR_heating、PDAR_coolingthe upper and lower limits of the power of the respective devices are satisfied:
in addition, the internal combustion engine also satisfies the hill climbing constraint:
|Pi,j+1 ICE-Pi,j ICE|≤0.5 PICE, forehead (12)。
In the following, taking a combined cooling heating and power type comprehensive energy system in a certain place as an example, simulation analysis is performed on the method.
(1) Model parameters
The model takes 4-10 months as summer, the illumination is strong, the temperature is high, the wind speed is low, the rest is winter, and no cold load exists in winter. In summer, the refrigerating equipment is double-effect absorption refrigerator and electric refrigerator. The heating equipment is a waste heat boiler and a gas boiler. In winter, the heating equipment is a double-effect absorption refrigerator, a waste heat boiler and a gas boiler. The planning problem is mixed platform integer planning, and the installed capacity of each device needs to be optimized and configured on an MATLAB platform by applying a yalcip tool box and a Cplex algorithm, so that the life whole-cycle cost of the system is the lowest.
TABLE 1 cost of equipment
(2) Typical day of winter and summer results analysis
And (4) carrying out optimization solution according to the optimization model, wherein the lowest cost for obtaining the life whole cycle is 1.2550e +8 yuan. The capacity of the equipment in the integrated energy system is shown in table 2:
TABLE 2 CCHP plant capacity planning results
1) Typical day of 1 month results analysis
A typical day of 1 month is selected for winter result analysis, the 1 month is winter, the system assumes no cold load in winter when modeling, and therefore the electric balance and thermal balance diagram of the typical day of 1 month is shown in attached figures 2 and 3.
As shown in the attached figure 2, the double-effect absorption refrigerator, the gas boiler and the waste heat boiler are used for heating, and the generated heat is supplied to a heat load to meet the heat balance. At night, the electric load is less, the power generation capacity of the internal combustion engine is less, the waste heat available in the system is less, the gas-fired boiler becomes the main heating power, the supplied heat accounts for a large proportion, and the gas-fired boiler bears most of the heat load particularly in the early morning (23-3 o' clock).
Fig. 3 illustrates that the generated energy and the purchased electric quantity of the internal combustion engine, the fan and the photovoltaic in the system are equal to the power demand and the electric load of the voltage compression refrigerator, so that the electric balance is met. In addition, the system can purchase electricity from the superior power grid only in the morning, and most of the electricity load at night is borne by the electricity purchasing quantity. Because there is no illumination night, the photovoltaic cell can not generate electricity, the wind power is small, the electricity generation is limited, the user load is not large, the electricity price is low, and the electricity purchasing from the superior power grid can easily meet the electricity demand in the system. When the electricity load is large in the daytime, most of the electricity demand of the system is borne by the internal combustion engine, and the small part of the electricity demand is borne by the photovoltaic and the fan, so that electricity is not purchased from a superior power grid. Particularly, at noon, the illumination intensity is the largest, the wind power is also large, and the power generation capacity of the distributed power supply is large relative to other times.
2)7 month typical day results analysis
A typical day of 7 months was selected as the summer outcome analysis with cold, heat and electrical balance as shown in figures 4, 5 and 6.
As can be seen from the attached figure 4, the cooling load is concentrated at 6-19 points in summer, no cooling load exists in other times, the cooling demand in the system in the time period is borne by the double-effect absorption refrigerator and the voltage compression refrigerator, and the voltage compression refrigerator is the main force of cooling. In summer, the weather is hot, the cooling load is large, and an electric compressor (COP is 4) with strong refrigerating capacity is required to bear the main refrigerating power, and in contrast, the COP of the double-effect absorption refrigerator is only 1.2, so that the electric compressor is reasonable to become the main refrigerating power of the system.
Fig. 5 illustrates that the heat load in summer and the heat supply of the gas boiler in winter are different. The heat load is less in summer, and the heat requirement of the system is met mainly by the waste heat boiler.
Fig. 6 illustrates that in summer, when the temperature and the illumination are greatly increased, particularly in the midday, the photovoltaic power generation capacity is sharply increased. The power load in summer is very large, particularly the power load in the middle of the day reaches the peak, at the moment, the system needs to purchase power from a higher-level power grid to meet the power demand of the system, and the power load balance diagram is obviously different from the power load balance diagram in 1 month. The electricity load balance diagram at night is the same as that in winter, and the electricity demand in the system is all from the power purchase of the superior power grid for the same reason.
(3) Compared with the traditional comprehensive energy system
A conventional integrated energy system employed herein includes: internal combustion engines, boilers, refrigeration equipment, electrical loads, thermal loads, cold loads, and the like. The CPLEX is adopted for modeling and optimizing calculation, the lowest cost of the obtained life whole cycle is 1.2747e +08 yuan, and compared with a combined cooling heating and power system considering intermittent distributed power sources, the cost is increased more. The comprehensive energy system planning system equipment capacity of wind-light-electricity multi-energy complementation and cooling-heating-electricity multi-energy flow coordinated transfer supply in the region is considered, the energy utilization rate is effectively improved, good economic benefits are obtained, and the operation investment is reduced. Table 3, table 4 compares the projected capacity and cost composition of the two integrated energy systems, respectively:
TABLE 3 CCHP vs. projected Capacity for conventional Integrated energy System devices
TABLE 4 CCHP vs. conventional Integrated energy System cost
As can be seen from table 4, the conventional integrated energy system uses no fan or photovoltaic cell, so that the number of devices is relatively small, and the construction cost and the device maintenance cost are lower than those of the combined cooling heating and power system considering the intermittent distributed power source. However, the traditional comprehensive energy system has high operation cost, because the resource utilization efficiency is low, the variable working condition characteristic of wind and light distributed energy is not considered, the wind, light and electricity complementation is lacked, a large amount of gas cost and electricity purchasing cost are needed, the total cost is very high, and from the long-term operation result, the comprehensive energy system has the advantages of higher economic benefit, high resource utilization rate, use of more clean energy, reduction of environmental pollution and capability of realizing sustainable development of resources and cascade utilization of energy.
In addition to the above embodiments, the present invention may have other embodiments, and any technical solutions formed by equivalent substitutions or equivalent transformations fall within the scope of the claims of the present invention.
Claims (1)
1. The optimized configuration method of the comprehensive energy system of the intermittent distributed power supply comprises the following steps of:
firstly, modeling each unit as follows:
(1) intermittent renewable energy source
1) Photovoltaic power generation
According to the photovoltaic effect of the semiconductor material, the photovoltaic cell can directly convert solar energy into electric energy, the temperature coefficient epsilon of the photovoltaic cell is 0.2, and a mathematical model is as follows:
f(t)=I+ε(T-Trated) (2)
in the formula, PPVRepresenting the photovoltaic output power, rratedIndicating the nominal light intensity, PratedRepresenting the photovoltaic power rating, TratedRepresents the nominal temperature, and T represents the photovoltaic cell temperature;
2) wind power generation
The long-term distribution of wind speed follows the Weibull distribution, which is commonly represented by a three-constant model:
wherein k is a Weibull shape parameter, c is a Weibull scale parameter, and v is a wind speed;
the mathematical model of wind power generation is as follows:
wherein v isci、vco、vratedRespectively representing cut-in wind speed, cut-out wind speed, rated wind speed, PWIndicating fan output power, PratedRepresents a rated power;
(2) combined supply equipment
1) Gas turbine
Variable working condition operation model of gas turbine:
η=aW3+bW2+cW (5-1)
in the formula, a is 1.071, b is-2.816, c is 2.745, eta and W are the specific values of the power generation efficiency and the output power of the gas turbine to the rated value, namely per unit; the formula (5) is a function of 3 times, and piecewise linearization is needed in CPLEX calculation; in addition, the climbing and descending constraints of the gas turbine are also considered, namely the absolute value of the difference value of the output of the gas turbine between two load values before and after the load change in actual operation is less than or equal to 0.5 times of rated power;
2) internal combustion engine
The variable working condition operation model of the internal combustion engine is as shown in the formula (5-2).
η=aW3+bW2+cW (5-2)
In the formula, a is 1.071, b is-2.816, c is 2.745, eta and W are the specific values of the power generation efficiency and the output power of the internal combustion engine to the rated value, namely the per unit value;
(3) refrigeration device
1) Double-effect absorption refrigerator
The refrigeration COP in summer is 1.2, namely 1.2 units of waste heat can be refrigerated every time the refrigeration COP absorbs 1 unit, the heating COP in winter is 0.9, namely 0.9 unit of waste heat can be heated every time the refrigeration COP absorbs 1 unit, the COP is the conversion rate between the absorbed energy and the released heat, which is called the heating energy efficiency ratio for short, and the variable working condition characteristic is not considered;
2) electric refrigerator
The refrigerating COP in summer is 4, namely 4 units of refrigerating can be performed by absorbing 1 unit of electricity, and the variable working condition characteristic is not considered;
(4) auxiliary equipment
1) Gas boiler
When the combined supply equipment cannot meet the heat demand of a user, the supplementary combustion of a gas-fired boiler is needed, the balance of supply and demand of heat energy is realized, and the variable working condition characteristic is not considered;
2) waste heat boiler
The equipment for changing water into steam by utilizing the heat of the flue gas discharged by the gas turbine does not consider the variable working condition characteristic;
secondly, setting the target function and the constraint condition as follows:
(1) objective function
After a photovoltaic power generation system and a wind power generation system are established, a randomness model is introduced into a comprehensive energy system model, an internal combustion engine is used as a prime mover, the total investment cost of the life whole period is used as an objective function, the objective function comprises the construction cost, the operation cost and the maintenance cost of each device, and the optimal capacity of the main devices of the system is optimized and configured; the objective function is as follows:
C=∑CT+∑Com+∑Cup (6)
in the formula, CTFor construction costs, ComFor operating costs, CupFor maintenance cost, the specific expressions of each part are as follows:
1)CTincluding all equipment converted to annual installed costs:
in the formula, IntRateFor interest rate, t is the current year, RconverseTo convert annual costs to a factor of the current year, RrecoverFor depreciation rate, MjFor construction unit price, RM,jIs the capacity of the device;
2)Comincluding the gas cost and the electricity purchasing cost when the system generated energy is insufficient:
in the formula, Cg,iFor gas costs of internal combustion engines and gas boilers, Ce,IFor the purchase of electricity for the voltage-compression refrigerator and for the upper-level power grid, the peak-valley rate system, P, is usedbuyIs the power purchase, PECRIs the electric power of the electric compressor, VICE、VGFBVolume of gas consumed by internal combustion engines and gas boilers, KgIs the price of natural gas, KeIs the price of electricity, NdayHours to use;
3)Cupincluding maintenance costs for all equipment;
(2) constraint conditions
The comprehensive energy system needs to meet 4 supply and demand balances, namely 4 equality constraint conditions, namely, electric balance, heat balance, cold balance and waste heat balance; the electric balance, the heat balance and the cold balance refer to electric quantity, heat quantity and cold quantity which are generated by equipment and can meet the electric demand, the heat demand and the cold demand of a system; the waste heat balance means that waste heat generated by a waste heat boiler and a gas boiler needs to meet the waste heat requirement of the double-effect absorption refrigerator; inequality constraints of equipment operation, namely equipment power limitation and climbing constraint of the internal combustion engine, need to be met;
1) power balance constraint
In the formula, PPV、PW、PICERespectively the generated power of the photovoltaic, the fan and the internal combustion engine, PLoad_power、PLoad_heating、PLoad_coolingRespectively electric, thermal and cold load, PbuyIs the power purchase, PECRIs the electrical power of the electric compressor, PDAR_heating、PWHB、PGFBRespectively the heating power, P, of the double-effect absorption refrigerator, the exhaust-heat boiler and the gas boilerDAR_cooling、PECR_coolingRefrigeration powers of the double-effect absorption refrigerator and the electric compressor, respectively, EICEFor the residual heat power of internal combustion engines, EDAR、EWHBThe heat absorption power of the double-effect absorption refrigerator and the waste heat boiler is respectively;
2) plant operating constraints
The light intensity, wind speed, temperature and the like are greatly influenced by natural environment and have strong randomness, if the total capacity of the distributed power supplies is overlarge, the power supply quality and stability of the system are influenced, so that the installation capacity of the distributed power supplies is limited, and the constraint condition is that
SPV、SWRespectively the apparent power of the photovoltaic and the fan,respectively representing the apparent power limit values of the photovoltaic and the draught fan;
PPV、PW、PICE、Pbuy、PECR、PWHB、PGFB、PDAR_heating、PDAR_coolingthe upper and lower limits of the power of the respective devices are satisfied:
in addition, the internal combustion engine also satisfies the hill climbing constraint:
|Pi,j+1 ICE-Pi,j ICE|≤0.5PICE, forehead (12)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110498724.3A CN113128070A (en) | 2021-05-08 | 2021-05-08 | Optimal configuration method for comprehensive energy system of intermittent distributed power supply |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110498724.3A CN113128070A (en) | 2021-05-08 | 2021-05-08 | Optimal configuration method for comprehensive energy system of intermittent distributed power supply |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113128070A true CN113128070A (en) | 2021-07-16 |
Family
ID=76781092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110498724.3A Pending CN113128070A (en) | 2021-05-08 | 2021-05-08 | Optimal configuration method for comprehensive energy system of intermittent distributed power supply |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113128070A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113610380A (en) * | 2021-08-02 | 2021-11-05 | 上海电气集团股份有限公司 | Multi-energy complementary energy planning method |
CN115062869A (en) * | 2022-08-04 | 2022-09-16 | 国网山东省电力公司东营供电公司 | Comprehensive energy scheduling method and system considering carbon emission |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100319354A1 (en) * | 2008-02-04 | 2010-12-23 | Gianfranco Guidati | Low carbon emissions combined cycle power plant and process |
CN110598313A (en) * | 2019-09-10 | 2019-12-20 | 国网河北省电力有限公司 | Comprehensive energy system optimization configuration method considering energy storage full-life cycle operation and maintenance |
-
2021
- 2021-05-08 CN CN202110498724.3A patent/CN113128070A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100319354A1 (en) * | 2008-02-04 | 2010-12-23 | Gianfranco Guidati | Low carbon emissions combined cycle power plant and process |
CN110598313A (en) * | 2019-09-10 | 2019-12-20 | 国网河北省电力有限公司 | Comprehensive energy system optimization configuration method considering energy storage full-life cycle operation and maintenance |
Non-Patent Citations (4)
Title |
---|
任娜等: "多能流分布式综合能源系统容量匹配优化与调度研究", 《电网技术》 * |
吴仁光等: "面向综合能源配电网的储能系统优化配置方法", 《广东电力》 * |
李琳等: "考虑随机性的区域综合能源系统多目标日前优化调度方法", 《电力建设》 * |
汤继刚等: "考虑间歇性分布式电源的综合能源系统优化配置研究", 《电力系统》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113610380A (en) * | 2021-08-02 | 2021-11-05 | 上海电气集团股份有限公司 | Multi-energy complementary energy planning method |
CN115062869A (en) * | 2022-08-04 | 2022-09-16 | 国网山东省电力公司东营供电公司 | Comprehensive energy scheduling method and system considering carbon emission |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111445090B (en) | Double-layer planning method for off-grid type comprehensive energy system | |
Yang et al. | Design and simulation of gas turbine-based CCHP combined with solar and compressed air energy storage in a hotel building | |
Su et al. | A new hybrid photovoltaic/thermal and liquid desiccant system for trigeneration application | |
Yang et al. | Load matching and techno-economic analysis of CSP plant with S–CO2 Brayton cycle in CSP-PV-wind hybrid system | |
Li et al. | An improved operation strategy for CCHP system based on high-speed railways station case study | |
CN111860937B (en) | Combined cooling heating and power type microgrid optimization method based on improved multi-target Husky algorithm | |
Wang et al. | Study on solar-assisted cascade refrigeration system | |
CN109543889A (en) | A kind of regional complex energy resource system cooperates with optimizing operation method a few days ago | |
Chen et al. | Multi-objective optimization of a combined cooling, heating, and power system with subcooled compressed air energy storage considering off-design characteristics | |
CN113128070A (en) | Optimal configuration method for comprehensive energy system of intermittent distributed power supply | |
CN111859683A (en) | Optimal configuration method of park comprehensive energy system based on dynamic energy concentrator | |
CN112446546B (en) | Comprehensive energy system two-stage optimal configuration method considering energy reliability | |
CN115994615A (en) | Comprehensive energy operation control method and system based on multi-energy complementation | |
Yu et al. | Complementary configuration research of new combined cooling, heating, and power system driven by renewable energy under energy management modes | |
Zheng et al. | Discharging strategy of adiabatic compressed air energy storage system based on variable load and economic analysis | |
CN110119850B (en) | Heat storage amount two-stage optimization scheduling method based on photo-thermal power generation adjustment | |
CN110783917A (en) | Configuration method of multi-energy hub containing new energy consumption | |
CN113313329B (en) | Optimal scheduling method for power distribution network containing comprehensive energy system | |
CN114997460A (en) | Regional micro-energy network operation optimization method considering maximum consumption of renewable energy | |
CN114825459A (en) | Energy scheduling method for load-intensive intelligent park users to participate in carbon emission control | |
CN114421536A (en) | Multi-energy interactive regulation and control method based on energy storage | |
Hengrui et al. | An effective planning approach for integrated energy systems considering equipment operating characteristics | |
Guo et al. | Multi-objective optimization of multi-energy complementary system based on cascade utilization of heat storage | |
Jinglin et al. | Multi-objective Collaborative Planning Method for Micro-energy Systems Considering Thermoelectric Coupling Clusters | |
He et al. | An optimization method of multiple energy flows for CCHP based on fuzzy theory and PSO |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210716 |