CN110361969B - Optimized operation method of cooling, heating and power comprehensive energy system - Google Patents

Abstract

The invention provides an optimized operation method of a cooling, heating and power comprehensive energy system, which comprises the following steps: 1) the method comprises the steps that the optimal economical efficiency of the overall operation of the cooling, heating and power comprehensive energy system is taken as a core, and a target function with the minimum total operation cost of the system is constructed by considering the multi-time scale characteristics of the cooling, heating and power comprehensive energy system; 2) considering the multi-time scale characteristics of the cooling, heating and power integrated energy system, establishing an equipment constraint model and a power balance constraint model as constraint conditions of a minimum objective function for the total running cost of the system; 3) and (3) solving the objective function with the minimum total running cost of the system by adopting a branch-and-bound method according to the constraint conditions in the step 2). The method provided by the invention aims at the complex structure and the operation mechanism of the cooling, heating and power comprehensive energy system, can improve the energy utilization efficiency, reduce the operation cost and realize the optimized operation of the cooling, heating and power comprehensive energy system.

Description

Optimized operation method of cooling, heating and power comprehensive energy system

Technical Field

The invention belongs to the field of comprehensive energy systems, and particularly relates to an optimized operation method of a cooling, heating and power comprehensive energy system.

Background

The combined cooling heating and power energy system is a combined production and supply system which is based on the concept of cascade utilization of energy and takes natural gas as primary energy to generate heat energy, electric energy and cold energy. The method takes natural gas as fuel, utilizes equipment such as a small gas turbine, a gas internal combustion engine, a micro-combustion engine and the like to combust the natural gas to obtain high-temperature flue gas which is firstly used for generating power and then utilizes waste heat to heat in winter; cooling in summer by driving the absorption refrigerator; meanwhile, domestic hot water can be provided, and exhaust heat is fully utilized. The utilization rate of primary energy can be improved to about 80 percent, and the primary energy is greatly saved.

The gas combined cooling heating and power system can be divided into a regional type and a building type according to the supply range. The regional system is mainly used for a cooling, heating and power energy supply center built in large regions such as various industrial, commercial or scientific parks. The equipment generally adopts a unit with larger capacity, an independent energy supply center is often required to be built, and external network equipment for supplying cold, heat and electricity is also required to be considered. The building type system is a cold and heat power supply system constructed for buildings with specific functions, such as office buildings, commercial buildings, hospitals and some comprehensive buildings, generally only needs a unit with smaller capacity, and machine rooms are usually arranged in the buildings without considering external network construction.

Compared with the traditional centralized power generation and remote power transmission modes, the gas combined cooling, heating and power supply can greatly improve the energy utilization efficiency: the generating efficiency of a large-scale power plant is generally 30-40%; the energy utilization efficiency of the cooling, heating and power comprehensive energy system is improved to 80-90%, and no power transmission loss exists.

The cooling, heating and power comprehensive energy system is essentially a cooling, heating and power multi-energy coupling system, is complex in structure and operation mechanism, has the characteristics of coexistence and interaction of various laws and variables, nonlinearity, uncertainty, multiple levels and the like, and is complex and diverse in system structure and working flow. At present, how to refer to the multi-time scale characteristics of a multi-energy coupling system of cold, heat and electricity, various energy sources are fully utilized in a gradient manner, efficient complementary supply of various energy sources such as cold, heat and electricity is realized, the energy utilization efficiency is improved, the operation cost is reduced, and the method is still a great problem in the operation process. Therefore, in order to solve the above problems, a specific solution needs to be provided to perfect the optimal operation of the cooling, heating and power integrated energy system.

Disclosure of Invention

Aiming at the problems, the invention provides a method for optimizing the operation of a cooling, heating and power comprehensive energy system, which comprises the following steps:

1) the optimization of the overall operation economy of the cooling, heating and power comprehensive energy system is taken as a core, and a system operation total cost minimum objective function is constructed on the basis of the multi-time scale characteristic of the cooling, heating and power comprehensive energy system;

2) establishing an equipment constraint model and a power balance constraint model based on the multi-time scale characteristics of the cooling, heating and power integrated energy system, wherein the equipment constraint model and the power balance constraint model are used as constraint conditions of a minimum objective function for the total running cost of the system;

3) and (3) solving the objective function with the minimum total running cost of the system by adopting a branch-and-bound method according to the constraint conditions in the step 2).

Wherein, the minimum objective function of the total running cost of the system in the step 1) is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; f_grid(t₁,t_2,i) Is at the t₁The system electricity purchasing cost in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; f_gas(t₁,t_2,i) Is at the t₁The cost for purchasing natural gas by the system in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; f_main(t₁,t_2,i) Is at the t₁Maintenance costs of system equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; f_poll(t₁,t_2,i) Is at the t₁The ith electric energy regulation period in the regulation period of the heat energy and the cold energyThe cost for discharging and treating the polluted gas in the device.

The electricity purchasing cost F of the system in the objective function with the minimum total running cost of the system_grid(t₁,t_2,i) Specifically, the following are shown:

F_grid(t₁,t_2,i)＝P_grid(t₁,t_2,i)gΔt₂gf_grid(t₁,t_2,i)

wherein, Δ t₂Time intervals for the power conditioning cycle; p_grid(t₁,t_2,i) Is at the t₁The electricity purchasing power of the system in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; f. of_grid(t₁,t_2,i) Is at the t₁Adjusting the real-time electricity price of the power grid in the ith electric energy adjusting period in the adjusting period of the heat energy and the cold energy;

the cost F for purchasing natural gas by the system in the objective function of minimum total operating cost of the system_gas(t₁,t_2,i) Specifically, the following are shown:

F_gas(t₁,t_2,i)＝V_gas(t₁,t_2,i)gΔt₂gf_gas(t₁,t_2,i)

wherein, V_gas(t₁,t_2,i) Is at the t₁The system of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy consumes the volume of the natural gas; Δ t₂Time intervals for the power conditioning cycle; f. of_gas(t₁,t_2,i) Is at the t₁The natural gas price of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy;

the system equipment maintenance cost F in the objective function of minimum total system operation cost_main(t₁,t_2,i) Specifically, the following are shown:

F_main(t₁,t_2,i)＝k_GE[P_GE(t₁,t_2,i)]gΔt₂gP_GE(t₁,t_2,i)+k_AP.cool[Q_AP.cool(t₁,t_2,i)]gΔt₂gQ_AP.cool(t₁,t_2,i)+k_AP.heat[Q_AP.heat(t₁,t_2,i)]gΔt₂gQ_AP.heat(t₁,t_2,i)+k_AC.heat[Q_AC.heat(t₁,t_2,i)]gΔt₂gQ_AC.heat(t₁,t_2,i)

wherein k is_GE[P_GE(t₁,t_2,i)]Is at the t₁Maintenance coefficients of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy under different output powers; p_GE(t₁,t_2,i) Is at the t₁The gas internal combustion engine outputs electric power in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; k is a radical of_AP.cool[Q_AP.cool(t₁,t_2,i)]Is at the t₁The cold power maintenance coefficient of the flue gas absorption heat pump equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.cool(t₁,t_2,i) Is at the t₁The flue gas of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the heat pump to output cold power; k is a radical of_AP.heat[Q_AP.heat(t₁,t_2,i)]Is at the t₁The thermal power maintenance coefficient of the flue gas absorption heat pump equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy;

is at the t₁The flue gas of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the heat power output by the heat pump; k is a radical of_AC.heat[Q_AC.heat(t₁,t_2,i)]Is at the t₁The maintenance coefficient of the absorption refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AC.heat(t₁,t_2,i) Is at the t₁The absorbed thermal power of the absorption refrigerator in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy;

said systemThe pollution gas emission abatement cost F in the objective function of minimum total system operating cost_poll(t₁,t_2,i) Specifically, the following are shown:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; λ is the number of pollutant emission types of the system, including: CO 2₂、SO₂、NO_x；_λTo comprise CO₂、SO₂、NO_xThe cost of abatement of various emissions therein; alpha is alpha_grid.λEmission coefficients for grid power versus different emissions; p_grid(t₁,t_2,i) Is at the t₁The electricity purchasing power of the system and the power grid in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy is obtained; alpha is alpha_GE.λEmission coefficients for different emissions for electric power of a gas internal combustion engine; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy.

The equipment constraint model in the step 2) comprises one or more models of a gas internal combustion engine constraint model, a cylinder sleeve water heat exchanger constraint model, an absorption refrigerator constraint model, an electric boiler constraint model, an electric refrigerator constraint model, a flue gas absorption heat pump equipment constraint model, an electric storage equipment constraint model, a heat storage equipment constraint model and a photovoltaic generator set constraint model.

The gas internal combustion engine constraint model comprises the following steps:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; generated power P of gas internal combustion engine_GE(t₁,t_2,i) At the t th₁The regulation period of 4 electric energy in the regulation period of the heat energy and the cold energy is constant; eta_GE.elec(t₁,t_2,i) Is at the t₁The power generation efficiency of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy is improved; p_maxThe rated power generation power of the gas internal combustion engine; q_GE.heat(t₁,t_2,i) Is at the t₁The thermal power output by the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy; thermal power Q output by gas internal combustion engine_GE.heat(t₁,t_2,i) At the t th₁The regulation period of 4 electric energy in the regulation period of the heat energy and the cold energy is constant; eta_LIs the inherent loss rate of the gas internal combustion engine; p_GE(t₁,t_2,i-1) the power generated by the gas combustion engine for the previous cycle of electrical energy regulation; p_GE.maxThe output gradient constraint of the gas internal combustion engine is carried out; LHV is the low calorific value of natural gas; eta_gasThe utilization rate of natural gas of the gas internal combustion engine is obtained; a is₃、a₂、a₁、a₀Respectively are fitting constants;

the constraint model of the flue gas absorption heat pump is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; t (T)₁,t_2,i) Is at the t₁Flue gas suction in the ith electric energy regulation period in the regulation period of heat energy and cold energyThe inlet temperature of the heat recovery pump; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_maxThe rated power generation power of the gas internal combustion engine; c_w(t₁,t_2,i) Is at the t₁The specific heat capacities of hot water with different temperatures in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy are respectively regulated; COP_AP(t₁,t_2,i) Is at the t₁The energy efficiency coefficient of the smoke absorption heat pump in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.heat(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the heating power of the heat pump;

is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the refrigeration power of the heat pump; q_AP.heat(t₁,t_2,i-1) adjusting the heating power of the cycle flue gas absorption heat pump for the last electrical energy; q_AP.cool(t₁,t_2,i-1) the refrigeration power of the flue gas absorption heat pump for the last electrical energy regulation cycle; heating power Q of flue gas absorption heat pump_AP.heat(t₁,t_2,i) And a refrigeration power Q_AP.cool(t₁,t_2,i) At the t th₁The regulation period of 4 electric energy in the regulation period of the heat energy and the cold energy is constant; lambda [ alpha ]_heat(t₁,t_2,i)、λ_cool(t₁,t_2,i) Are respectively the t-th₁The flue gas heating proportion and the refrigerating proportion of the flue gas absorption heat pump in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy are regulated; t is_heat、T_coolRespectively the hot water outlet temperature and the cold water outlet temperature; l is_heat(t₁,t_2,i)、L_cool(t₁,t_2,i) Are respectively the t-th₁The hot water and the cold water of the flue gas absorption heat pump in the ith electric energy regulation period in the regulation period of the heat energy and the cold energyWater flow rate; l is_heat.max、L_cool.maxMaximum heating and refrigerating flows are respectively; eta_AP.heat、η_AP.coolRespectively the heating efficiency and the refrigerating efficiency of the flue gas absorption heat pump; q_AP.heat.maxThe output gradient constraint of the heating power of the flue gas absorption heat pump is carried out; q_AP.cool.maxThe refrigeration power output gradient of the flue gas absorption heat pump is restrained; b₅、b₄、b₃、b₂、b₁、b₀Respectively are fitting constants;

the constraint model of the cylinder sleeve water heat exchanger is as follows:

Q_JW(t₁,t_2,i)＝η_JW gQ_GE.heat(t₁,t_2,i)

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_JW(t₁,t_2,i) Is at the t₁The cylinder sleeve water heat exchanger in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy outputs heat power; eta_JWThe heat exchange efficiency of the cylinder sleeve water heat exchanger is obtained; q_GE.heat(t₁,t_2,i) Is at the t₁The thermal power output by the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy;

the absorption refrigerator constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; q_ac.heat(t₁) Is at the t₁The absorption refrigerating machine absorbs heat power in a regulation period of heat energy and cold energy; q_ac.cool(t₁) Is at the t₁The cold power output by the absorption refrigerator with the regulation period of the heat energy and the cold energy; q_ac.cool(t₁-1) absorption chiller refrigeration power for the last heat and cold conditioning cycle; COP_acIs the energy efficiency coefficient of the absorption refrigerator; q_ac.heat.min、Q_ac.heat.maxRespectively the minimum and maximum heat power absorbed by the absorption refrigerator; q_ac.cool.maxIs the output gradient constraint of the absorption refrigerator;

the electric boiler constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EB(t₁,t_2,i) Is at the t₁Inputting electric power to the electric boiler in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_EB(t₁,t_2,i) Is at the t₁The electric boiler outputs heat power in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_EB(t₁,t_2,i-1) regulating the output thermal power of the electric boiler for the previous electric energy cycle; COP_EBThe energy production coefficient of the electric boiler; p_EB.min、P_EB.maxRespectively the minimum electric power and the maximum electric power of the electric boiler; q_EB.maxThe output gradient constraint of the electric boiler is carried out;

the electric refrigerator constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EC(t₁,t_2,i) Is at the t₁The input electric power of the electric refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_EC(t₁,t_2,i) Is at the t₁For the i-th electric regulation period within the regulation periods of heat and coldThe output cold power of the electric refrigerator; q_EC(t₁,t_2,i-1) adjusting the output cold power of the electric refrigerator for the last electric energy regulation cycle; COP_ECIs the energy efficiency coefficient of the electric refrigerator; p_EC.min、P_EC.maxRespectively the minimum electric power and the maximum electric power of the electric refrigerator; q_EC.maxIs the output gradient constraint of the electric refrigerator;

the photovoltaic generator set constraint model comprises the following steps:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_PV(t₁,t_2,i) Is at the t₁Adjusting the real-time power of the photovoltaic generator set in the ith electric energy adjusting period in the adjusting period of the heat energy and the cold energy; p_STCRated output of the photovoltaic generator set; g_ING(t₁,t_2,i) Is at the t₁Adjusting the real-time irradiation intensity of the ith electric energy adjusting period in the adjusting periods of the heat energy and the cold energy; g_STCThe rated irradiation intensity of the photovoltaic generator set; k is the power generation coefficient of the photovoltaic generator set; t is_out(t₁,t_2,i) Is at the t₁The external temperature of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; t is_sIs the reference temperature of the generator set;

the power storage equipment constraint model comprises the following steps:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; e_batt(t₁,t_2,i) Is at the t₁The real-time capacity of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; k is a radical of_LThe electric energy self-loss coefficient of the electricity storage equipment is obtained; eta_batt.chaThe charging efficiency of the electric storage device; eta_batt.disThe discharge efficiency of the electric storage device; p_batt.cha(t₁,t_2,i) Is at the t₁Charging power of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_batt.dis(t₁,t_2,i) Is at the t₁The discharge power of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_batt.dis.max、P_batt.dis.minThe maximum and small discharge power of the power storage equipment are respectively; p_batt.cha.max、P_batt.cha.minThe maximum charging power and the minimum charging power of the power storage equipment are respectively; e_batt.max、E_batt.minThe maximum and minimum electricity storage capacities of the electricity storage equipment are respectively set;

the heat storage equipment constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; b is_stor(t₁,t_2,i) Is at the t₁The real-time capacity of the heat storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy is regulated; k is a radical of_sThe heat energy self-loss coefficient of the heat storage equipment; eta_stor.chaThe heat absorption efficiency of the heat storage device; eta_stor.disThe heat release efficiency of the heat storage device; q_stor.cha(t₁,t_2,i) Is at the t₁The heat absorption power of the heat storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy is regulated; q_stor.dis(t₁,t_2,i) Is at the t₁Storage of the i-th electric energy regulation cycle within the regulation cycle of thermal and cold energyThe heat release power of the thermal device; q_stor.cha.max、Q_stor.cha.minThe maximum and minimum heat absorption power of the heat storage equipment are respectively; q_stor.dis.max、Q_stor.dis.minThe maximum and minimum heat release power of the heat storage equipment respectively; b is_stor.max、B_stor.minThe maximum and minimum heat storage capacities of the heat storage device are respectively.

The device power balance constraint model in the step 2) comprises an electric power balance constraint model, a thermal power balance constraint model and a cold power balance constraint model.

The electric power balance constraint model is as follows:

P_grid(t₁,t_2,i)+P_PV(t₁,t_2,i)+P_GE(t₁,t_2,i)+P_batt.dis(t₁,t_2,i)gD_batt.dis(t₁,t_2,i)＝P_batt.cha(t₁,t_2,i)gD_batt.cha(t₁,t_2,i)+P_ele(t₁,t_2,i)+P_EB(t₁,t_2,i)+P_EC(t₁,t_2,i)

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_grid(t₁,t_2,i) Is at the t₁The power of the power grid in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_PV(t₁,t_2,i) Is at the t₁The real-time power of the photovoltaic unit in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy is regulated; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_batt.dis(t₁,t_2,i)、P_batt.cha(t₁,t_2,i) Are respectively the t-th₁Discharge and charge power of the electricity storage device in the i-th electric energy regulation period within the regulation periods of thermal and cold energy, D_batt.dis(t₁,t_2,i)、D_batt.cha(t₁,t_2,i) Are respectively the t-th₁The discharge and charge variables of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy are changed; p_ele(t₁,t_2,i) Is at the t₁The electric load of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EB(t₁,t_2,i) Is at the t₁The electric power consumption of the electric boiler of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EC(t₁,t_2,i) Is at the t₁The electric refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy consumes electric power;

the thermal power balance constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_JW(t₁,t_2,i) Is at the t₁The output thermal power of the cylinder sleeve water heat exchanger in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy; q_AP.heat(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the output heat power of the heat pump; q_EB(t₁,t_2,i) Is at the t₁The output heat power of the electric boiler in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_stor.dis(t₁)、Q_stor.cha(t₁) Are respectively the t-th₁Heat-release and heat-absorption power of heat-storage devices for a controlled period of thermal and cold energy, D_stor.dis(t₁)、D_stor.cha(t₁) Are respectively the t-th₁The heat release and absorption variables of the heat storage equipment in the regulation period of the heat energy and the cold energy; q_heat(t₁) Is at the t₁Thermal load of the regulation cycle of the individual heat and cold energies; q_AC.heat(t₁) Is at the t₁The absorption refrigerating machine absorbs heat power in a regulation period of heat energy and cold energy;

the cold power balance constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AC.cool(t₁) Is at the t₁The cold power output by the absorption refrigerator with the regulation period of the heat energy and the cold energy; q_EC(t₁,t_2,i) Is at the t₁The cold power output by the electric refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.cool(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the cold power output by the heat pump; q_cool(t₁) Is at the t₁The system outputs cold power to the outside in a regulation period of heat energy and cold energy.

The solving in the step 3) specifically comprises the following steps:

inputting known parameters, relaxing constraint conditions, and decomposing the original problem into a plurality of sub-problems;

solving the subproblems, judging whether the solved subproblem solution is a feasible solution, and if so, ending the calculation process;

if not, setting the sub-problem solution as an original problem upper bound, setting the maximum feasible solution target as an original problem lower bound, and comparing the upper bound with the lower bound;

if the upper bound is larger than the lower bound, solving the subproblem again; if the upper bound is smaller than the lower bound, the original problem is solved, and the calculation process is finished.

The invention provides a set of specific solutions aiming at the structure and the operation mechanism of the cooling, heating and power comprehensive energy system. Based on the multi-time scale characteristic of the cooling, heating and power comprehensive energy system, various energy sources are fully and stepwisely utilized, efficient complementary supply of various energy sources such as cold, heat and power is achieved, the energy utilization efficiency is improved, the operation cost is reduced, and the optimized operation of the cooling, heating and power comprehensive energy system is achieved.

Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.

Fig. 1 shows a topological structure diagram of a cooling, heating and power integrated energy system according to an embodiment of the present invention;

fig. 2 shows a flowchart of an optimal operation constraint model solving method for a cooling, heating and power integrated energy system according to an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

The working principle of the cooling, heating and power comprehensive energy system is as follows:

fig. 1 is a topological structure diagram of a cooling, heating and power integrated energy system according to an embodiment of the present invention. As shown in fig. 1, the main devices of the cooling, heating and power integrated energy system include a gas internal combustion engine, a cylinder liner water heat exchanger, an absorption refrigerator, an electric boiler, an electric refrigerator, a flue gas absorption heat pump device, and two energy storage devices, namely an electricity storage device and a heat storage device, are provided, in order to improve the permeability of renewable energy, the system is further connected to a photovoltaic generator set and is connected to an electric network to ensure that sufficient electric energy is supplied to a power load.

The comprehensive energy system of cooling, heating and power is in micro-energy network level and takes a gas internal combustion engine as a core system. The gas internal combustion engine directly supplies electric power to a power load by consuming natural gas and generating electric energy. The hot steam generated by the gas internal combustion engine during working is converted into hot water through the cylinder liner water heat exchanger to supply a thermal load. Meanwhile, the flue gas generated during the combustion of the natural gas can be received by the flue gas absorption heat pump, and the flue gas absorption heat pump operates to convert the flue gas into heat energy and cold energy to be directly supplied to users.

When the heat energy supply is sufficient and the cold energy supply is insufficient, in order to compensate the cold energy supply insufficiency, partial heat energy absorbed by the absorption type refrigerating machine can be converted into cold energy to be supplied to a load for use. When the electric energy is sufficient to supply cold or insufficient heat, the electric boiler can absorb the electric energy to convert into heat energy or the electric refrigerator can absorb the electric energy to convert into cold energy for supplement. And electricity storage and heat storage equipment is added in the system to ensure that the system has enough power capacity margin so as to ensure the stability of the system. In addition, the active access of the photovoltaic generator set improves the permeability of new energy of the system and increases the environmental protection and economic benefits of the system. When the electric energy load demand is large, the system can interact with the power grid, but in order to reduce the construction cost and coordination cost of the system, the power grid information channel and the physical channel, the system in the embodiment adopts a principle of grid connection and no network access, and electric energy is purchased from the power grid to make up for the shortage of the electric energy of the system and ensure the stable operation of the system.

P in FIG. 1_eleRepresents the electric power output by the system; p_GERepresenting the power generation of the gas internal combustion engine; p_gridRepresenting an electric networkPower; p_PVRepresenting the generated power of the photovoltaic generator set; p_EBRepresenting the input electrical power of the electrical boiler; p_ECRepresenting the input electrical power of the electrical refrigerator; q_GEIndicating the output thermal power of the gas internal combustion engine; q_JWRepresenting the output thermal power of the cylinder liner water heat exchanger; q_EBRepresenting the output thermal power of the electric boiler; q_AP.heatThe output thermal power of the flue gas absorption heat pump is represented; q_heatRepresents the thermal power output by the system; q_AC.heatRepresenting the input thermal power of the absorption chiller; q_AC.coolRepresenting the output cold power of the absorption chiller; q_AP.coolThe output cold power of the flue gas absorption heat pump is represented; q_ECRepresenting the output cold power of the electric refrigerator; q_coolIndicating the cold power output by the system.

The optimization operation of the cooling, heating and power comprehensive energy system with multiple time scales is considered:

in order to improve the comprehensive energy efficiency of the cooling, heating and power comprehensive energy system, realize stable and balanced output of cooling, heating and power energy and reduce the operation cost, the invention provides an optimized operation method of the cooling, heating and power comprehensive energy system, which comprises the following contents:

because the time scales of cold, heat and electricity output by each device in the cooling, heating and power integrated energy system are different, the cooling, heating and power integrated energy system has the characteristic of multiple time scales. In the embodiment of the present invention, the adjustment period of the heat energy and the cold energy is set to 1 hour, and the adjustment period of the electric energy is set to 15 minutes, that is, in one adjustment period of the heat energy and the cold energy, 4 electric energy adjustment periods are included, the heat energy and the cold energy are respectively adjusted at the whole point of each hour, and the electric energy is adjusted at the 0 th minute, the 15 th minute, the 30 th minute and the 45 th minute of each hour.

The method provided by the invention combines the multi-time scale characteristic of the cooling, heating and power comprehensive energy system, and constructs the objective function with the minimum total running cost of the system by taking the optimal overall running economy of the cooling, heating and power comprehensive energy system as the core.

Meanwhile, the method combines the multi-time scale characteristic of the cooling, heating and power comprehensive energy system, and establishes an equipment constraint model for key equipment including a gas internal combustion engine, a cylinder sleeve water heat exchanger, an absorption refrigerator, an electric boiler, an electric refrigerator, a flue gas absorption heat pump device, an electricity storage device, a heat storage device and a photovoltaic generator set.

In addition, in order to meet the power balance of cold power, heat power and electric power in the system, the method establishes a power balance constraint model by combining the multi-time scale characteristic of the cooling, heating and power comprehensive energy system.

The method takes the equipment constraint model and the power balance constraint model as constraint conditions, adopts a branch-and-bound method to solve the minimum objective function of the total running cost of the system, the solved result is the minimum value of the total running cost of the system, and the equipment constraint conditions and the power balance constraint conditions which meet the minimum value of the total running cost of the system are the optimal running scheme of the system. This example will further illustrate the method in detail:

the system running total cost minimum objective function is as follows:

in the formula (1.1), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; f_grid(t₁,t_2,i) Is at the t₁The system electricity purchasing cost in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; f_gas(t₁,t_2,i) Is at the t₁The cost for purchasing natural gas by the system in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; f_main(t₁,t_2,i) Is at the t₁Maintenance costs of system equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; f_poll(t₁,t_2,i) Is at the t₁In the regulation cycle of heat energy and cold energyAnd (4) pollution gas emission treatment cost in i electric energy regulation periods.

In the formula (1.1), the system electricity purchasing cost is specifically expressed as follows:

F_grid(t₁,t_2,i)＝P_grid(t₁,t_2,i)gΔt₂gf_grid(t₁,t_2,i) (1.2)

in the formula (1.2), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; p_grid(t₁,t_2,i) Is at the t₁The electricity purchasing power of the system and the power grid in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy is regulated; f. of_grid(t₁,t_2,i) Is at the t₁And the real-time electricity price of the power grid in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy. The value at delta t can be obtained by solving the formula (1.2)₂The electricity purchase cost of the system in the meantime.

In the formula (1.1), the cost for purchasing natural gas by the system is specifically expressed as follows:

F_gas(t₁,t_2,i)＝V_gas(t₁,t_2,i)gΔt₂gf_gas(t₁,t_2,i) (1.3)

in the formula (1.3), V_gas(t₁,t_2,i) Is at the t₁The system consumption natural gas volume in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; f. of_gas(t₁,t_2,i) Is at the t₁The natural gas price in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy. The value at delta t can be obtained by solving the formula (1.3)₂The cost of purchasing natural gas from the system in the interim.

In the formula (1.1), the system equipment maintenance cost is specifically expressed as follows:

in the formula (1.4), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; k is a radical of_GE[P_GE(t₁,t_2,i)]Is at the t₁Maintenance coefficients of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy under different output powers; p_GE(t₁,t_2,i) Is at the t₁Outputting electric power by the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; k is a radical of_AP.cool[Q_AP.cool(t₁,t_2,i)]Is at the t₁The cold power maintenance coefficient of the flue gas absorption heat pump equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.cool(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the cold power output by the heat pump; k is a radical of_AP.heat[Q_AP.heat(t₁,t_2,i)]Is at the t₁The thermal power maintenance coefficient of the flue gas absorption heat pump equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.heat(t₁,t_2,i) Is at the t₁The flue gas of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the heat power output by the heat pump; k is a radical of_AC.heat[Q_AC.heat(t₁,t_2,i)]Is at the t₁The maintenance coefficient of the absorption refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AC.heat(t₁,t_2,i) Is at the t₁The electric energy in the ith regulation period of the heat energy and the cold energy regulates the heat power absorbed by the absorption refrigerator in the ith regulation period. The value at delta t can be obtained by solving the formula (1.4)₂The maintenance costs of the system equipment in the meantime.

In the formula (1.1), the pollution gas emission control cost is specifically expressed as follows:

in the formula (1.5), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; λ is the number of types of pollutant emissions of the system, said pollutant emissions comprising: CO 2₂、SO₂、NO_x；_λTo comprise CO₂、SO₂、NO_xThe cost of remediation of the pollutant emissions contained; alpha is alpha_grid.λEmission coefficients for grid power versus different emissions; p_grid(t₁,t_2,i) Is at the t₁The electricity purchasing power of the system and the power grid in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy is regulated; alpha is alpha_GE.λEmission coefficients for different emissions for electric power of a gas internal combustion engine; p_GE(t₁,t_2,i) Is at the t₁The power generated by the gas combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy. The value at delta t can be obtained by solving the formula (1.5)₂The pollution gas discharge treatment cost in the period.

The constraint condition of the objective function with minimum total running cost of the system mainly comprises equipment constraint model constraint and power balance constraint model constraint.

For the equipment constraint model, the method of the invention combines the multi-time scale characteristic of the cooling, heating and power comprehensive energy system, and establishes the equipment constraint model for key equipment including a gas internal combustion engine, a cylinder sleeve water heat exchanger, an absorption refrigerator, an electric boiler, an electric refrigerator, a flue gas absorption heat pump equipment, an electricity storage equipment, a heat storage equipment and a photovoltaic generator set in the system respectively. The constraint model of the device is specifically as follows:

the gas internal combustion engine constraint model:

in the constraint model formula (2.1) of the gas internal combustion engine, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; generated power P of gas internal combustion engine_GE(t₁,t_2,i) At the t th₁Constant for 4 electric energy regulation periods (namely one hour) in the regulation periods of the heat energy and the cold energy; eta_GE.elec(t₁,t_2,i) The power generation efficiency of the gas internal combustion engine; p_maxThe rated power generation power of the gas internal combustion engine; q_GE.heat(t₁,t_2,i) Is at the t₁The thermal power output by the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy; thermal power Q output by gas internal combustion engine_GE.heat(t₁,t_2,i) At the t th₁Constant for 4 electric energy regulation periods (namely one hour) in the regulation periods of the heat energy and the cold energy; eta_LIs the inherent loss rate of the gas internal combustion engine; p_GE(t₁,t_2,i-1) the power generated by the gas combustion engine for the previous cycle of electrical energy regulation; p_GE.maxThe output gradient constraint of the gas internal combustion engine is carried out; LHV is the low calorific value of natural gas; eta_gasThe utilization rate of natural gas of the gas internal combustion engine is obtained; a is₃、a₂、a₁、a₀Respectively, fitting constants.

The flue gas absorption heat pump constraint model is as follows:

in the constraint model formula (2.2) of the flue gas absorption heat pump, t is₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁In the regulation cycle of heat energy and cold energyi power regulation cycles; t (T)₁,t_2,i) Is at the t₁The inlet temperature of the flue gas absorption heat pump in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_maxThe rated power generation power of the gas internal combustion engine; c_w(t₁,t_2,i) Is at the t₁The specific heat capacities of hot water with different temperatures in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy are respectively regulated; COP_AP(t₁,t_2,i) Is at the t₁The energy efficiency coefficient of the smoke absorption heat pump in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.heat(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the heating power of the heat pump; q_AP.cool(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy absorbs the refrigeration power of the heat pump; q_AP.heat(t₁,t_2,i-1) adjusting the heating power of the cycle flue gas absorption heat pump for the last electrical energy; q_AP.cool(t₁,t_2,i-1) the refrigeration power of the flue gas absorption heat pump for the last electrical energy regulation cycle; heating power Q of flue gas absorption heat pump_AP.heat(t₁,t_2,i) And a refrigeration power Q_AP.cool(t₁,t_2,i) At the t th₁Constant for 4 electric energy regulation periods (namely one hour) in the regulation periods of the heat energy and the cold energy; lambda [ alpha ]_heat(t₁,t_2,i)、λ_cool(t₁,t_2,i) Are respectively the t-th₁The heating proportion and the refrigerating proportion of the flue gas absorption heat pump in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy are regulated; t is_heat、T_coolRespectively the hot water outlet temperature and the cold water outlet temperature; l is_heat(t₁,t_2,i)、L_cool(t₁,t_2,i) Are respectively the t-th₁The flue gas in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy absorbs the hot water and cold water flow of the heat pump; l is_heat.max、L_cool.maxMaximum heating and refrigerating flows are respectively; eta_AP.heat、η_AP.coolRespectively the heating efficiency and the refrigerating efficiency of the flue gas absorption heat pump; q_AP.heat.maxThe output gradient constraint of the heating power of the flue gas absorption heat pump is carried out; q_AP.cool.maxThe refrigeration power output gradient of the flue gas absorption heat pump is restrained; b₅、b₄、b₃、b₂、b₁、b₀Respectively, fitting constants.

Constraint model of cylinder liner water heat exchanger:

Q_JW(t₁,t_2,i)＝η_JWgQ_GE.heat(t₁,t_2,i) (2.3)

in the constraint model formula (2.3) of the cylinder liner water heat exchanger, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_JW(t₁,t_2,i) Is at the t₁Outputting thermal power by the cylinder sleeve water heat exchanger in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; eta_JWThe heat exchange efficiency of the cylinder sleeve water heat exchanger is obtained; q_GE.heat(t₁,t_2,i) Is at the t₁The electric energy in the ith regulation period of the heat energy and the cold energy regulates the thermal power output by the gas combustion engine in the ith regulation period.

Absorption chiller constraint model:

in the absorption chiller constraint model equation (2.4), t₁Represents the conditioning cycle of the heat and cold energy; q_ac.heat(t₁) Is at the t₁The absorption refrigerating machine absorbs heat power in a regulation period of heat energy and cold energy; q_ac.cool(t₁) Is at the t₁With a period of regulation of heat and coldThe cold power output by the absorption refrigerator; q_ac.cool(t₁-1) absorption chiller refrigeration power for the last heat and cold conditioning cycle; COP_acIs the energy efficiency coefficient of the absorption refrigerator; q_ac.heat.min、Q_ac.heat.maxRespectively the minimum and maximum heat power absorbed by the absorption refrigerator; q_ac.cool.maxIs the output gradient constraint of the absorption chiller.

Electric boiler constraint model:

in the constraint model formula (2.5) of the electric boiler, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EB(t₁,t_2,i) Is at the t₁Inputting electric power into an electric boiler in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_EB(t₁,t_2,i) Is at the t₁The electric boiler in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy outputs heat power; q_EB(t₁,t_2,i-1) regulating the output thermal power of the electric boiler for the previous electric energy cycle; COP_EBThe energy production coefficient of the electric boiler; p_EB.min、P_EB.maxRespectively the minimum electric power and the maximum electric power of the electric boiler; q_EB.maxIs the output gradient constraint of the electric boiler.

Electric refrigerator constraint model:

in the electric refrigerator constraint model formula (2.6), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EC(t₁,t_2,i) Is at the t₁Of heat and coldThe input electric power of the electric refrigerator in the ith electric energy regulation period in the regulation period; q_EC(t₁,t_2,i) Is at the t₁The output cold power of the electric refrigerator in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; q_EC(t₁,t_2,i-1) adjusting the output cold power of the electric refrigerator for the last electric energy regulation cycle; COP_ECIs the energy efficiency coefficient of the electric refrigerator; p_EC.min、P_EC.maxRespectively the minimum electric power and the maximum electric power of the electric refrigerator; q_EC.maxIs the output gradient constraint of the electric refrigerator.

Photovoltaic generator set restraint model:

in the constraint model formula (2.7) of the photovoltaic generator set, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_PV(t₁,t_2,i) Is at the t₁Adjusting the real-time power of the photovoltaic generator set in the ith electric energy adjusting period in the adjusting period of the heat energy and the cold energy; p_STCRated output of the photovoltaic generator set; g_ING(t₁,t_2,i) Is at the t₁Adjusting the real-time irradiation intensity in the ith electric energy adjusting period in the adjusting periods of the heat energy and the cold energy; g_STCThe rated irradiation intensity of the photovoltaic generator set; k is the power generation coefficient of the photovoltaic generator set; t is_out(t₁,t_2,i) Is at the t₁The external temperature in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; t is_sIs the reference temperature of the generator set.

The power storage equipment constraint model is as follows:

constraint model formula (2.8) of power storage equipment) In, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; e_batt(t₁,t_2,i) Is at the t₁The real-time capacity of the electricity storage equipment in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; k is a radical of_LThe electric energy self-loss coefficient of the electricity storage equipment is obtained; eta_batt.chaThe charging efficiency of the electric storage device; eta_batt.disThe discharge efficiency of the electric storage device; p_batt.cha(t₁,t_2,i) Is at the t₁Charging power of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_batt.dis(t₁,t_2,i) Is at the t₁The discharge power of the electricity storage equipment in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; p_batt.dis.max、P_batt.dis.minThe maximum and small discharge power of the power storage equipment are respectively; p_batt.cha.max、P_batt.cha.minThe maximum charging power and the minimum charging power of the power storage equipment are respectively; e_batt.max、E_batt.minThe maximum and minimum electric storage capacities of the electric storage device are respectively.

The heat storage equipment constraint model is as follows:

in the heat storage equipment constraint model (2.9), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; b is_stor(t₁,t_2,i) Is at the t₁The real-time capacity of the heat storage equipment in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy is regulated; k is a radical of_sThe heat energy self-loss coefficient of the heat storage equipment; eta_stor.chaThe heat absorption efficiency of the heat storage device; eta_stor.disThe heat release efficiency of the heat storage device; q_stor.cha(t₁,t_2,i) Is at the t₁The heat absorption power of the heat storage equipment in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy is regulated; q_stor.dis(t₁,t_2,i) Is at the t₁The heat release power of the heat storage equipment in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; q_stor.cha.max、Q_stor.cha.minThe maximum and minimum heat absorption power of the heat storage equipment are respectively; q_stor.dis.max、Q_stor.dis.minThe maximum and minimum heat release power of the heat storage equipment respectively; b is_stor.max、B_stor.minThe maximum and minimum heat storage capacities of the heat storage device are respectively.

The constraint condition of the objective function with the minimum total running cost of the system mainly comprises that besides the equipment constraint model, the method provided by the invention also establishes a power balance constraint model comprising an electric power balance constraint model, a thermal power balance constraint model and a cold power balance constraint model in order to satisfy the power balance of cold, heat and electricity. The power balance constraint model is specifically as follows:

electric power balance constraint model:

in the electric power balance constraint model formula (3.1), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_grid(t₁,t_2,i) Is at the t₁The power of the power grid in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; p_PV(t₁,t_2,i) Is at the t₁Adjusting the real-time power of the photovoltaic unit in the ith electric energy adjusting period in the adjusting period of the heat energy and the cold energy; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_batt.dis(t₁,t_2,i)、P_batt.cha(t₁,t_2,i) Are respectively the t-th₁Discharge and charge power of the electricity storage device in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy, D_batt.dis(t₁,t_2,i)、D_batt.cha(t₁,t_2,i) Are respectively the t-th₁The discharge and charge variables of the electricity storage equipment in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy are regulated; p_ele(t₁,t_2,i) Is a power load; p_EB(t₁,t_2,i) Is at the t₁The consumed electric power of the electric boiler in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EC(t₁,t_2,i) Is at the t₁The electric refrigerator in the ith electric power conditioning period in the conditioning periods of the heat energy and the cold energy consumes electric power.

Thermal power balance constraint model:

in the thermal power balance constraint model formula (3.2), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_JW(t₁,t_2,i) Is at the t₁The output thermal power of the cylinder sleeve water heat exchanger in the ith electric energy regulation period in the regulation period of the thermal energy and the cold energy; q_AP.heat(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy absorbs the output heat power of the heat pump; q_EB(t₁,t_2,i) Is at the t₁The output thermal power of the electric boiler in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_stor.dis(t₁)、Q_stor.cha(t₁) Are respectively the t-th₁Heat-release and heat-absorption power of the heat storage device during the individual cycle of regulation of thermal and cold energy, D_stor.dis(t₁)、D_stor.cha(t₁) Are respectively the t-th₁The heat release and absorption variables of the heat storage equipment in the regulation period of the heat energy and the cold energy; q_heat(t₁) Is at the t₁Thermal load during the conditioning cycle of the individual heat and cold energies; q_AC.heat(t₁) Is at the t₁The absorption refrigerating machine absorbs heat power in the regulation period of the heat energy and the cold energy.

Cold power balance constraint model:

in the cold power balance constraint model equation (3.3), t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AC.cool(t₁) Is at the t₁The cold power output by the absorption refrigerator in the regulation period of the heat energy and the cold energy; q_EC(t₁,t_2,i) Is at the t₁The cold power output by the electric refrigerator in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy;

is at the t₁The flue gas in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy absorbs the cold power output by the heat pump; q_cool(t₁) Is at the t₁The cold power output by the system in the regulation period of the heat energy and the cold energy is reduced.

The optimized operation constraint model of the cooling, heating and power comprehensive energy system provided by the invention has higher order and large dimension, and is difficult to calculate by a common solving method. Therefore, the optimization model is solved by adopting a branch-and-bound algorithm. Fig. 2 is a flowchart of a method for solving an optimized operation constraint model of a cooling, heating and power integrated energy system provided by the invention, and the specific calculation process is as follows:

firstly, inputting known parameters, relaxing constraint conditions and decomposing the original problem into a plurality of subproblems. Then, solving an optimal solution aiming at the subproblems, judging whether the solved subproblem optimal solution is a feasible solution or not, if so, taking the conclusion as the optimal solution, and ending the calculation process; if not, setting the optimal solution as an upper boundary of the original problem, setting the maximum target of the feasible solution as a lower boundary of the original problem, and comparing the upper boundary with the lower boundary. If the upper bound is larger than the lower bound, returning to the subproblem solving step, and solving the optimal solution of the subproblem again; if the upper bound is smaller than the lower bound, the original problem is solved, and the calculation process is finished.

The method of the invention provides a set of specific solutions aiming at the complex structure and the operation mechanism of the cooling, heating and power comprehensive energy system and based on the multi-time scale characteristic of the cooling, heating and power comprehensive energy system. The system has the advantages of fully and stepwisely utilizing various energy sources, realizing efficient complementary supply of various energy sources such as cold, heat, electricity and the like, improving the energy utilization efficiency, reducing the operation cost and realizing the optimized operation of the cooling, heating and power comprehensive energy system.

Although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (3)

1. An optimal operation method of a cooling, heating and power integrated energy system is characterized by comprising the following steps:

1) the optimization of the overall operation economy of the cooling, heating and power comprehensive energy system is taken as a core, and a system operation total cost minimum objective function is constructed on the basis of the multi-time scale characteristic of the cooling, heating and power comprehensive energy system;

2) establishing an equipment constraint model and a power balance constraint model based on the multi-time scale characteristics of the cooling, heating and power integrated energy system, wherein the equipment constraint model and the power balance constraint model are used as constraint conditions of a minimum objective function for the total running cost of the system;

the equipment constraint model in the step 2) comprises one or more models of a gas internal combustion engine constraint model, a cylinder sleeve water heat exchanger constraint model, an absorption refrigerator constraint model, an electric boiler constraint model, an electric refrigerator constraint model, a flue gas absorption heat pump equipment constraint model, an electric storage equipment constraint model, a heat storage equipment constraint model and a photovoltaic generator set constraint model;

the gas internal combustion engine constraint model comprises the following steps:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; generated power P of gas internal combustion engine_GE(t₁,t_2,i) At the t th₁The regulation period of 4 electric energy in the regulation period of the heat energy and the cold energy is constant; eta_GE.elec(t₁,t_2,i) Is at the t₁The power generation efficiency of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy is improved; p_maxThe rated power generation power of the gas internal combustion engine; q_GE.heat(t₁,t_2,i) Is at the t₁The thermal power output by the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy; thermal power Q output by gas internal combustion engine_GE.heat(t₁,t_2,i) At the t th₁The regulation period of 4 electric energy in the regulation period of the heat energy and the cold energy is constant; eta_LIs the inherent loss rate of the gas internal combustion engine; p_GE(t₁,t_2,i-1) the power generated by the gas combustion engine for the previous cycle of electrical energy regulation; p_GE.maxThe output gradient constraint of the gas internal combustion engine is carried out; LHV is the low calorific value of natural gas; eta_gasThe utilization rate of natural gas of the gas internal combustion engine is obtained; a is₃、a₂、a₁、a₀Respectively are fitting constants;

the constraint model of the flue gas absorption heat pump is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; t (T)₁,t_2,i) Is at the t₁The inlet temperature of the flue gas absorption heat pump in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_maxThe rated power generation power of the gas internal combustion engine; c_w(t₁,t_2,i) Is at the t₁The specific heat capacities of hot water with different temperatures in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy are respectively regulated; COP_AP(t₁,t_2,i) Is at the t₁The energy efficiency coefficient of the smoke absorption heat pump in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.heat(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the heating power of the heat pump; q_AP.cool(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the refrigeration power of the heat pump; q_AP.heat(t₁,t_2,i-1) adjusting the heating power of the cycle flue gas absorption heat pump for the last electrical energy; q_AP.cool(t₁,t_2,i-1) the refrigeration power of the flue gas absorption heat pump for the last electrical energy regulation cycle; heating power Q of flue gas absorption heat pump_AP.heat(t₁,t_2,i) And a refrigeration power Q_AP.cool(t₁,t_2,i) At the t th₁A heat energy andthe cold energy is constant in 4 electric energy regulation periods in the regulation period; lambda [ alpha ]_heat(t₁,t_2,i)、λ_cool(t₁,t_2,i) Are respectively the t-th₁The flue gas heating proportion and the refrigerating proportion of the flue gas absorption heat pump in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy are regulated; t is_heat、T_coolRespectively the hot water outlet temperature and the cold water outlet temperature; l is_heat(t₁,t_2,i)、L_cool(t₁,t_2,i) Are respectively the t-th₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the hot water and cold water flow of the heat pump; l is_heat.max、L_cool.maxMaximum heating and refrigerating flows are respectively; eta_AP.heat、η_AP.coolRespectively the heating efficiency and the refrigerating efficiency of the flue gas absorption heat pump; q_AP.heat.maxThe output gradient constraint of the heating power of the flue gas absorption heat pump is carried out; q_AP.cool.maxThe refrigeration power output gradient of the flue gas absorption heat pump is restrained; b₅、b₄、b₃、b₂、b₁、b₀Respectively are fitting constants;

the constraint model of the cylinder sleeve water heat exchanger is as follows:

Q_JW(t₁,t_2,i)＝η_JWgQ_GE.heat(t₁,t_2,i)

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_JW(t₁,t_2,i) Is at the t₁The cylinder sleeve water heat exchanger in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy outputs heat power; eta_JWThe heat exchange efficiency of the cylinder sleeve water heat exchanger is obtained; q_GE.heat(t₁,t_2,i) Is at the t₁The thermal power output by the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy;

the absorption refrigerator constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; q_ac.heat(t₁) Is at the t₁The absorption refrigerating machine absorbs heat power in a regulation period of heat energy and cold energy; q_ac.cool(t₁) Is at the t₁The cold power output by the absorption refrigerator with the regulation period of the heat energy and the cold energy; q_ac.cool(t₁-1) absorption chiller refrigeration power for the last heat and cold conditioning cycle; COP_acIs the energy efficiency coefficient of the absorption refrigerator; q_ac.heat.min、Q_ac.heat.maxRespectively the minimum and maximum heat power absorbed by the absorption refrigerator; q_ac.cool.maxIs the output gradient constraint of the absorption refrigerator;

the electric boiler constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EB(t₁,t_2,i) Is at the t₁Inputting electric power to the electric boiler in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_EB(t₁,t_2,i) Is at the t₁The electric boiler outputs heat power in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_EB(t₁,t_2,i-1) regulating the output thermal power of the electric boiler for the previous electric energy cycle; COP_EBThe energy production coefficient of the electric boiler; p_EB.min、P_EB.maxRespectively the minimum electric power and the maximum electric power of the electric boiler; q_EB.maxThe output gradient constraint of the electric boiler is carried out;

the electric refrigerator constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EC(t₁,t_2,i) Is at the t₁The input electric power of the electric refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_EC(t₁,t_2,i) Is at the t₁The output cold power of the electric refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_EC(t₁,t_2,i-1) adjusting the output cold power of the electric refrigerator for the last electric energy regulation cycle; COP_ECIs the energy efficiency coefficient of the electric refrigerator; p_EC.min、P_EC.maxRespectively the minimum electric power and the maximum electric power of the electric refrigerator; q_EC.maxIs the output gradient constraint of the electric refrigerator;

the photovoltaic generator set constraint model comprises the following steps:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_PV(t₁,t_2,i) Is at the t₁Adjusting the real-time power of the photovoltaic generator set in the ith electric energy adjusting period in the adjusting period of the heat energy and the cold energy; p_STCRated output of the photovoltaic generator set; g_ING(t₁,t_2,i) Is at the t₁Adjusting the real-time irradiation intensity of the ith electric energy adjusting period in the adjusting periods of the heat energy and the cold energy; g_STCThe rated irradiation intensity of the photovoltaic generator set; k is the power generation coefficient of the photovoltaic generator set; t is_out(t₁,t_2,i) Is at the t₁During a conditioning cycle of heat and coldThe outside temperature of the ith electric energy regulation period; t is_sIs the reference temperature of the generator set;

the power storage equipment constraint model comprises the following steps:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; e_batt(t₁,t_2,i) Is at the t₁The real-time capacity of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; k is a radical of_LThe electric energy self-loss coefficient of the electricity storage equipment is obtained; eta_batt.chaThe charging efficiency of the electric storage device; eta_batt.disThe discharge efficiency of the electric storage device; p_batt.cha(t₁,t_2,i) Is at the t₁Charging power of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_batt.dis(t₁,t_2,i) Is at the t₁The discharge power of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_batt.dis.max、P_batt.dis.minThe maximum and small discharge power of the power storage equipment are respectively; p_batt.cha.max、P_batt.cha.minThe maximum charging power and the minimum charging power of the power storage equipment are respectively; e_batt.max、E_batt.minThe maximum and minimum electricity storage capacities of the electricity storage equipment are respectively set;

the heat storage equipment constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation period of the heat energy and the cold energy；Δt₂Time intervals for the power conditioning cycle; b is_stor(t₁,t_2,i) Is at the t₁The real-time capacity of the heat storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy is regulated; k is a radical of_sThe heat energy self-loss coefficient of the heat storage equipment; eta_stor.chaThe heat absorption efficiency of the heat storage device; eta_stor.disThe heat release efficiency of the heat storage device; q_stor.cha(t₁,t_2,i) Is at the t₁The heat absorption power of the heat storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy is regulated; q_stor.dis(t₁,t_2,i) Is at the t₁The heat release power of the heat storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_stor.cha.max、Q_stor.cha.minThe maximum and minimum heat absorption power of the heat storage equipment are respectively; q_stor.dis.max、Q_stor.dis.minThe maximum and minimum heat release power of the heat storage equipment respectively; b is_stor.max、B_stor.minMaximum and minimum heat storage capacities of the heat storage device are respectively set;

the power balance constraint model in the step 2) comprises an electric power balance constraint model, a thermal power balance constraint model and a cold power balance constraint model;

the electric power balance constraint model is as follows:

P_grid(t₁,t_2,i)+P_PV(t₁,t_2,i)+P_GE(t₁,t_2,i)+P_batt.dis(t₁,t_2,i)gD_batt.dis(t₁,t_2,i)＝

P_batt.cha(t₁,t_2,i)gD_batt.cha(t₁,t_2,i)+P_ele(t₁,t_2,i)+P_EB(t₁,t_2,i)+P_EC(t₁,t_2,i)

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_grid(t₁,t_2,i) Is at the t₁The power of the power grid in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_PV(t₁,t_2,i) Is at the t₁The real-time power of the photovoltaic unit in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy is regulated; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_batt.dis(t₁,t_2,i)、P_batt.cha(t₁,t_2,i) Are respectively the t-th₁Discharge and charge power of the electricity storage device in the i-th electric energy regulation period within the regulation periods of thermal and cold energy, D_batt.dis(t₁,t_2,i)、D_batt.cha(t₁,t_2,i) Are respectively the t-th₁The discharge and charge variables of the electricity storage equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy are changed; p_ele(t₁,t_2,i) Is at the t₁The electric load of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EB(t₁,t_2,i) Is at the t₁The electric power consumption of the electric boiler of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; p_EC(t₁,t_2,i) Is at the t₁The electric refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy consumes electric power;

the thermal power balance constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_JW(t₁,t_2,i) Is at the t₁The output thermal power of the cylinder sleeve water heat exchanger in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy; q_AP.heat(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the output heat power of the heat pump; q_EB(t₁,t_2,i) Is at the t₁The output heat power of the electric boiler in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_stor.dis(t₁)、Q_stor.cha(t₁) Are respectively the t-th₁Heat-release and heat-absorption power of heat-storage devices for a controlled period of thermal and cold energy, D_stor.dis(t₁)、D_stor.cha(t₁) Are respectively the t-th₁The heat release and absorption variables of the heat storage equipment in the regulation period of the heat energy and the cold energy; q_heat(t₁) Is at the t₁Thermal load of the regulation cycle of the individual heat and cold energies; q_AC.heat(t₁) Is at the t₁The absorption refrigerating machine absorbs heat power in a regulation period of heat energy and cold energy;

the cold power balance constraint model is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AC.cool(t₁) Is at the t₁The cold power output by the absorption refrigerator with the regulation period of the heat energy and the cold energy; q_EC(t₁,t_2,i) Is at the t₁The cold power output by the electric refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.cool(t₁,t_2,i) Is at the t₁The flue gas in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the cold power output by the heat pump; q_cool(t₁) Is at the t₁The cold power output by the system is regulated according to the regulation period of the heat energy and the cold energy;

3) solving the objective function with the minimum total running cost of the system by adopting a branch-and-bound method according to the constraint conditions in the step 2); the solving specifically comprises the following steps:

inputting known parameters, relaxing constraint conditions, and decomposing the original problem into a plurality of sub-problems;

solving the subproblems, judging whether the solved subproblem solution is a feasible solution, and if so, ending the calculation process;

if not, setting the sub-problem solution as an original problem upper bound, setting the maximum feasible solution target as an original problem lower bound, and comparing the upper bound with the lower bound;

if the upper bound is larger than the lower bound, solving the subproblem again; if the upper bound is smaller than the lower bound, the original problem is solved, and the calculation process is finished.

2. The method for optimizing the operation of a system according to claim 1, wherein the objective function for minimizing the total cost of operating the system in step 1) is as follows:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; f_grid(t₁,t_2,i) Is at the t₁The system electricity purchasing cost in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; f_gas(t₁,t_2,i) Is at the t₁The cost for purchasing natural gas by the system in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; f_main(t₁,t_2,i) Is at the t₁Maintenance costs of system equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; f_poll(t₁,t_2,i) Is at the t₁Contaminated gas exhaust in the ith electric power conditioning cycle of the individual thermal and cold power conditioning cyclesAnd (5) controlling the treatment cost.

3. The method of claim 2, wherein the system electricity purchase cost F in the objective function of minimizing the total system operation cost is the system electricity purchase cost_grid(t₁,t_2,i) Specifically, the following are shown:

F_grid(t₁,t_2,i)＝P_grid(t₁,t_2,i)gΔt₂gf_grid(t₁,t_2,i)

wherein, Δ t₂Time intervals for the power conditioning cycle; p_grid(t₁,t_2,i) Is at the t₁The electricity purchasing power of the system in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy; f. of_grid(t₁,t_2,i) Is at the t₁Adjusting the real-time electricity price of the power grid in the ith electric energy adjusting period in the adjusting period of the heat energy and the cold energy;

the cost F for purchasing natural gas by the system in the objective function of minimum total operating cost of the system_gas(t₁,t_2,i) Specifically, the following are shown:

F_gas(t₁,t_2,i)＝V_gas(t₁,t_2,i)gΔt₂gf_gas(t₁,t_2,i)

wherein, V_gas(t₁,t_2,i) Is at the t₁The system of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy consumes the volume of the natural gas; Δ t₂Time intervals for the power conditioning cycle; f. of_gas(t₁,t_2,i) Is at the t₁The natural gas price of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy;

the system equipment maintenance cost F in the objective function of minimum total system operation cost_main(t₁,t_2,i) Specifically, the following are shown:

F_main(t₁,t_2,i)＝k_GE[P_GE(t₁,t_2,i)]gΔt₂gP_GE(t₁,t_2,i)+k_AP.cool[Q_AP.cool(t₁,t_2,i)]gΔt₂gQ_AP.cool(t₁,t_2,i)+k_AP.heat[Q_AP.heat(t₁,t_2,i)]gΔt₂gQ_AP.heat(t₁,t_2,i)+k_AC.heat[Q_AC.heat(t₁,t_2,i)]gΔt₂gQ_AC.heat(t₁,t_2,i)

wherein k is_GE[P_GE(t₁,t_2,i)]Is at the t₁Maintenance coefficients of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy under different output powers; p_GE(t₁,t_2,i) Is at the t₁The gas internal combustion engine outputs electric power in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; k is a radical of_AP.cool[Q_AP.cool(t₁,t_2,i)]Is at the t₁The cold power maintenance coefficient of the flue gas absorption heat pump equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.cool(t₁,t_2,i) Is at the t₁The flue gas of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the heat pump to output cold power; k is a radical of_AP.heat[Q_AP.heat(t₁,t_2,i)]Is at the t₁The thermal power maintenance coefficient of the flue gas absorption heat pump equipment in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AP.heat(t₁,t_2,i) Is at the t₁The flue gas of the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy absorbs the heat power output by the heat pump; k is a radical of_AC.heat[Q_AC.heat(t₁,t_2,i)]Is at the t₁The maintenance coefficient of the absorption refrigerator in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; q_AC.heat(t₁,t_2,i) Is at the t₁The absorbed thermal power of the absorption refrigerator in the ith electric energy regulation period in the regulation periods of the thermal energy and the cold energy;

the above-mentionedThe pollution gas emission abatement cost F in the objective function of minimum total system operation cost_poll(t₁,t_2,i) Specifically, the following are shown:

wherein, t₁Represents the conditioning cycle of the heat and cold energy; t is t_2,iRepresents the t th₁The ith electric energy regulation period in the regulation periods of the heat energy and the cold energy; Δ t₂Time intervals for the power conditioning cycle; λ is the number of pollutant emission types of the system, including: CO 2₂、SO₂、NO_x；_λTo comprise CO₂、SO₂、NO_xThe cost of abatement of various emissions therein; alpha is alpha_grid.λEmission coefficients for grid power versus different emissions; p_grid(t₁,t_2,i) Is at the t₁The electricity purchasing power of the system and the power grid in the ith electric energy regulation period in the regulation period of the heat energy and the cold energy is obtained; alpha is alpha_GE.λEmission coefficients for different emissions for electric power of a gas internal combustion engine; p_GE(t₁,t_2,i) Is at the t₁The power generation power of the gas internal combustion engine in the ith electric energy regulation period in the regulation periods of the heat energy and the cold energy.

Images