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

Abstract

The present invention provides a method for optimizing operation of a cooling, heating and power integrated energy system. The method includes the following steps: 1) Taking the overall economical optimization of the cooling, heating and power integration energy system as the core, considering the cooling, heating and power integration The multi-time-scale characteristics of the energy system construct the minimum objective function of the total operating cost of the system; 2) Considering the multi-time-scale characteristics of the integrated cooling, heating and power energy system, an equipment constraint model and a power balance constraint model are established, which are used as requirements for the operation of the system. Constraints of the total cost minimum objective function; 3) Using the branch and bound method, according to the constraints in step 2), the system operation total cost minimum objective function is solved. Aiming at the complex structure and operation mechanism of the cold-thermal-electricity integrated energy system, the method of the invention can improve the energy utilization efficiency, reduce the operation cost, and realize the optimal operation of the cold-thermal-electricity integrated 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. A method for optimizing operation of a cooling, heating and power integrated energy system, characterized in that the method comprises the steps: 1) Taking the optimization of the overall operation economy of the integrated cooling, heating and power system as the core, and constructing the minimum objective function of the total operating cost of the system based on the multi-time scale characteristics of the integrated cooling and heating energy system; 2) Based on the multi-time scale characteristics of the integrated cooling, heating and power system, an equipment constraint model and a power balance constraint model are established as constraints on the minimum objective function of the total operating cost of the system; The equipment constraint model in the step 2) includes a gas internal combustion engine constraint model, a cylinder jacket water heat exchanger constraint model, an absorption refrigerator constraint model, an electric boiler constraint model, an electric refrigerator constraint model, and a flue gas absorption heat pump device constraint model. one or more of a model, an electrical storage device constraint model, a thermal storage device constraint model, and a photovoltaic generator set constraint model; The gas engine constraint model is as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat energy and cold energy adjustment period; Δt ₂ is the time interval of the electric energy adjustment period; P _GE (t ₁ , _{t 2 , i} ) is the power generation of the gas-fired internal combustion engine in the i-th electric energy regulation cycle in the t ₁ -th thermal energy and cold energy regulation cycle _{; i} ) Constant in 4 electric energy regulation cycles in the t1th heat and cold energy regulation cycle; η _GE.elec (t ₁ ,t _2,i ) is in the _{t 1th} heat and cold energy regulation cycle The power generation efficiency of the gas-fired internal combustion engine in the i-th electric energy regulation cycle; P _max is the rated power generation of the gas-fired internal combustion engine; Q _GE.heat (t ₁ , t _{2 , i} ) is the t ₁ -th thermal and cold energy regulation cycle The thermal power output by the gas internal combustion engine in the i-th electric energy regulation cycle; the thermal power Q _GE.heat (t ₁ ,t _2,i ) output by the gas internal combustion engine in the t _1th thermal energy and cold energy regulation cycle 4 Constant in the electric energy adjustment cycle; η _L is the inherent loss rate of the gas engine; P _GE (t ₁ , t _{2, i} -1) is the power generated by the gas engine in the previous electric energy adjustment cycle; P _GE.max is the output of the gas engine slope constraint; LHV is the low calorific value of natural gas; η _gas is the natural gas utilization rate of the gas engine; a ₃ , a ₂ , a ₁ , and a ₀ are fitting constants respectively; The flue gas absorption heat pump constraint model is as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat energy and cold energy adjustment period; T(t ₁ ,t _2,i ) is the ith electric energy adjustment period The inlet temperature of the flue gas absorption heat pump in the ith electric energy regulation cycle in the _t1 heat and cold energy regulation cycle; P _GE (t ₁ ,t _2,i ) is the t _1th heat energy and cold energy regulation The power generation of the gas-fired internal combustion engine in the ith electric energy regulation cycle in the cycle; _Pmax is the rated power of the gas-fired internal combustion engine; _Cw (t ₁ , t _{2 , i} ) is the t ₁ heat and cold energy in the regulation cycle _The specific heat capacity of hot water at different temperatures in _{the ith} electric energy regulation cycle _of The energy efficiency coefficient of the absorption heat pump; Q _AP.heat (t ₁ ,t _2,i ) is the heating power of the flue gas absorption heat pump in the i-th electric energy regulation cycle in the t ₁ -th heat and cold energy regulation cycle ; Q _AP.cool (t ₁ , t _{2, i} ) is the cooling power of the flue gas absorption heat pump in the i-th electric energy regulation cycle in the t ₁ -th heat and cold energy regulation cycle; Q _AP.heat (t ₁ , t _{2, i} -1) is the heating power of the flue gas absorption heat pump in the previous electric energy adjustment cycle; Q _AP.cool (t ₁ , t _{2, i} -1) is the flue gas absorption heat pump in the previous electric energy adjustment cycle Cooling power; heating power Q _AP.heat (t ₁ ,t _2,i ) and cooling power Q _AP.cool (t ₁ , _{t 2,i} ) of the flue gas absorption heat pump It is constant in the 4 electric energy regulation cycles in the regulation cycle; λ _heat (t ₁ ,t _2,i ) and λ _cool (t ₁ ,t _2,i ) are respectively in the regulation cycle of the t _1th heat energy and cold energy The heating ratio and cooling ratio of the flue gas absorption heat pump in the ith electric energy regulation cycle; T _heat and T _cool are the hot water outlet temperature and the cold water outlet temperature respectively; L _heat (t ₁ ,t _2,i ), L _cool (t ₁ , t _{2 , i} ) are the hot and cold water flows of the flue gas absorption heat pump in the i-th electric energy regulation cycle in the t ₁ -th heat energy and cold energy regulation cycles, respectively; L _heat.max , L _{cool .max} is the maximum heating and cooling flow respectively; η _AP.heat and η _AP.cool are the heating and cooling efficiencies of the flue gas absorption heat pump respectively; Q _AP.heat.max is the heating power output gradient of the flue gas absorption heat pump constraints; Q _AP.cool.max is the cooling power output gradient constraint of the flue gas absorption heat pump; b ₅ , b ₄ , b ₃ , b ₂ , b ₁ , and b ₀ are fitting constants respectively; The constraint model of the cylinder jacket water heat exchanger is as follows: Q _JW (t ₁ ,t _2,i )=η _JW gQ _GE.heat (t ₁ ,t _2,i ) Among them, t ₁ represents the adjustment period of thermal energy and cold energy; t _2,i represents the ith electric energy adjustment period in the adjustment period of t ₁ th thermal energy and cold energy; Q _JW (t ₁ ,t _2,i ) is The output heat power of the liner water heat exchanger in the ith electric energy regulation cycle in the t _1th heat and cold energy regulation cycle; η _JW is the heat exchange efficiency of the liner water heat exchanger; Q _GE.heat (t ₁ , t _{2, i} ) is the thermal power output by the gas-fired internal combustion engine in the ith electric energy regulation cycle in the t ₁ th thermal energy and cold energy regulation cycle; The absorption chiller constraint model is as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; Q _ac.heat (t ₁ ) is the thermal power absorbed by the absorption chiller in the t _1th adjustment period of heat energy and cold energy; Q _ac.cool (t ₁ ) ) is the cooling power output by the absorption chiller in the t1th heat and cold energy regulation cycle; Q _ac.cool (t _{1 -1} ) is the absorption chiller cooling power in the last heat and cold energy regulation cycle; COP _ac is the energy efficiency coefficient of the absorption chiller; Q _ac.heat.min and Q _ac.heat.max are the minimum and maximum thermal power absorbed by the absorption chiller respectively; Q _ac.cool.max is the absorption chiller Output slope constraint; The electric boiler constraint model is as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat energy and cold energy adjustment period; P _EB (t ₁ ,t _2,i ) is The input electric power _of the electric boiler in the _ith electric energy adjustment cycle in the _{t1th heat} energy and cooling energy adjustment cycle _; The output thermal power of the electric boiler in the i-th electric energy regulation cycle; Q _EB (t ₁ , t _{2, i} -1) is the output thermal power of the electric boiler in the previous electric energy regulation cycle; COP _EB is the energy production coefficient of the electric boiler; P _EB.min and P _EB.max are the minimum and maximum electric power of the electric boiler respectively; Q _EB.max is the output gradient constraint of the electric boiler; The electric refrigerator constraint model is as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat energy and cold energy adjustment period; P _EC (t ₁ ,t _2,i ) is The input electric power of the electric refrigerator in the ith electric energy regulation cycle in the t1th heat and cold energy regulation cycle; Q _EC (t ₁ ,t _2,i ) is the _{t 1th} heat and cold energy regulation cycle The output cooling power _of the electric refrigerator in _{the ith electric} energy regulation cycle in energy efficiency coefficient; P _EC.min and P _EC.max are the minimum and maximum electric power of the electric refrigerator respectively; Q _EC.max is the output gradient constraint of the electric refrigerator; The photovoltaic generator set constraint model is as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat energy and cold energy adjustment period; P _PV (t ₁ ,t _2,i ) is The real-time power of the photovoltaic generator set in the ith electric energy adjustment cycle in the t1th thermal energy and cold energy adjustment cycle; P _STC is the rated output of the photovoltaic generator set; G _ING (t ₁ , t _{2 , i} ) is the ith t is the real-time irradiance intensity of the i _- th electric energy regulation cycle within the regulation cycle of thermal energy and cold energy; G _STC is the rated irradiance intensity of the photovoltaic generator set; k is the power generation coefficient of the photovoltaic generator set; T _out (t ₁ , t _{2, i} ) is the outside temperature of the ith electric energy regulation cycle in the t _1th thermal energy and cold energy regulation cycle; T _s is the reference temperature of the generator set; The energy storage device constraint model is as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat energy and cold energy adjustment period; Δt ₂ is the time interval of the electric energy adjustment period; E _batt (t ₁ , t _{2 , i} ) is the real-time capacity of the power storage device in the i-th electric energy regulation cycle in the t ₁ -th thermal energy and cold energy regulation cycle; k _L is the self-loss coefficient of the power storage device; η _{batt .cha} is the charging efficiency of the electrical storage device; η _batt.dis is the discharging efficiency of the electrical storage device; P _batt.cha (t ₁ , t _{2 , i} ) is the t ₁ th thermal and cold energy regulation cycle The charging power of the power storage device in the i electric energy adjustment cycle; P _batt.dis (t ₁ , t _{2, i} ) is the power storage device in the i th electric energy adjustment cycle in the t ₁ th heat energy and cold energy adjustment cycle P _{batt.dis.max and P batt.dis.min are the maximum and minimum discharge power of the power storage device respectively; P batt.cha.max and P batt.cha.min are the} maximum and minimum charging power of the power storage device respectively Power; E _batt.max and E _batt.min are the maximum and minimum storage capacity of the power storage device respectively; The thermal storage device constraint model is as follows:

Among them, t ₁ represents the adjustment period of thermal energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t ₁ -th thermal energy and cold energy adjustment period; Δt ₂ is the time interval of the electric energy adjustment period; B _stor (t ₁ , t _{2, i} ) is the real-time capacity of the heat storage device in the i-th electric energy regulation cycle in the _{t 1} -th heat and cold energy regulation cycle; ks is the thermal energy self-loss coefficient of the heat storage device; η _{stor .cha} is the heat absorption efficiency of the heat storage device; η _stor.dis is the heat release efficiency of the heat storage device; Q _stor.cha (t ₁ , t _{2 , i} ) is the t ₁ heat energy and cold energy regulation cycle in the The endothermic power of the heat storage device in the ith electric energy regulation cycle; Q _stor.dis (t ₁ ,t _2,i ) is the storage power of the i th electric energy regulation cycle in the t ₁ th thermal energy and cold energy regulation cycle The heat release power of the heat device; Q _stor.cha.max and Q _stor.cha.min are the maximum and minimum heat absorption power of the heat storage device respectively; Q _stor.dis.max and Q _stor.dis.min are the heat storage device The maximum and minimum heat release power of the equipment; B _stor.max and B _stor.min are the maximum and minimum heat storage capacity of the heat storage equipment respectively; The power balance constraint model in the step 2) includes 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 )+ _PE B (t ₁ ,t _2,i ) +P _EC (t ₁ ,t _2,i ) Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat and cold energy adjustment period; P _grid (t ₁ ,t _2,i ) is The grid power of the ith electric energy regulation cycle in the t1th heat and cold energy regulation cycle; P _PV (t ₁ ,t _2,i ) is the i th electricity in the _{t 1th} heat and cold energy regulation cycle The real-time power of photovoltaic units in each electric energy regulation cycle; P _GE (t ₁ , t _{2 , i} ) is the power generation power of the gas-fired internal combustion engine in the ith electric energy regulation cycle in the t ₁ th thermal energy and cold energy regulation cycle; P _{batt .dis} (t ₁ ,t _2,i ), P _batt.cha (t ₁ ,t _2,i ) are the power storage devices of the i-th electric energy regulation cycle in the t ₁ -th thermal energy and cold energy regulation cycles, respectively The discharge and charging power of , D _batt.dis (t ₁ , t _{2, i} ), D _batt.cha (t ₁ , t _{2, i} ) are the i-th regulation cycle of the t ₁ -th thermal energy and cold energy, respectively The discharge and charging variables of the power storage equipment in the electric energy regulation cycle; P _ele (t ₁ , t _{2 , i} ) is the power load of the i-th electric energy regulation cycle in the t ₁ -th thermal energy and cold energy regulation cycle; P _EB (t ₁ , t _{2, i} ) is the electric power consumption of the electric boiler in the ith electric energy regulation cycle in the t ₁ th thermal energy and cooling energy regulation cycle; P _EC (t ₁ , t _{2, i} ) is the ith electric power consumption t The electric power consumed by the electric refrigerator in the i _- th electric energy adjustment cycle within the adjustment cycle of heat energy and cold energy; The thermal power balance constraint model is as follows:

Among them, t ₁ represents the adjustment period of thermal energy and cold energy; t _2,i represents the ith electric energy adjustment period in the adjustment period of t ₁ th thermal energy and cold energy; Q _JW (t ₁ ,t _2,i ) is The output thermal power of the liner water heat exchanger in the ith electric energy regulation cycle in the t1th heat and cold energy regulation cycle; Q _AP.heat (t ₁ ,t _{2 ,i} ) is the t _1th thermal energy The output heat power of the flue gas absorption heat pump in the i-th electric energy regulation cycle in the regulation cycle of cold energy and cold energy; Q _EB (t ₁ , t _{2 , i} ) is the t The output thermal power of the electric boiler in i electric energy regulation cycles; Q _stor.dis (t ₁ ) and Q _stor.cha (t ₁ ) are the heat release of the heat storage device in the t ₁ -th regulation cycle of thermal energy and cooling energy, respectively , endothermic power, D _stor.dis (t ₁ ), D _stor.cha ( _{t 1} ) are the exothermic and endothermic variables of the heat storage device in the t1th heat energy and cold energy regulation cycle, respectively; Q _heat ( t ₁ ) is the thermal load of the t _1th regulation cycle of heat energy and cooling energy; Q _AC.heat (t ₁ ) is the thermal power absorbed by the absorption chiller in the t _1th regulation cycle of heat energy and cooling energy; The cold power balance constraint model is as follows:

Among them, t ₁ represents the regulation cycle of heat energy and cold energy; t _2,i represents the ith electric energy regulation cycle in the t ₁ th thermal energy and cold energy regulation cycle; Q _AC.cool (t ₁ ) is the t ₁ th electrical energy regulation cycle The cooling power output by the absorption chiller in the regulation cycle of heat and cold energy; Q _EC (t ₁ ,t _2,i ) is the ith electric energy regulation cycle in the regulation cycle of t ₁ heat and cold energy. Cooling power output by the electric refrigerator; Q _AP.cool (t ₁ , t _{2, i} ) is the cooling power output by the flue gas absorption heat pump in the i-th electric energy regulation cycle in the t ₁ -th heat and cold energy regulation cycle ; Q _cool (t ₁ ) is the cooling power output by the system in the t _1th regulation cycle of thermal energy and cold energy; 3) adopt the branch and bound method to solve the minimum objective function of the total operating cost of the system according to the constraints in the step 2); the solution specifically includes the following steps: Input known parameters, relax constraints, and decompose the original problem into many sub-problems; For the sub-problem solution, determine whether the obtained sub-problem solution is a feasible solution, if the determination result is yes, the calculation process ends; If the judgment is no, set the solution of the sub-problem as the upper bound of the original problem, set the maximum objective of the feasible solution as the lower bound of the original problem, and compare the upper bound with the lower bound; If the upper bound is greater than the lower bound, the sub-problem is re-solved; if the upper bound is less than the lower bound, the original problem has no solution, and the calculation process ends. 2. system optimization operation method according to claim 1, is characterized in that, in described step 1), system operation total cost minimum objective function is as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat energy and cold energy adjustment period; F _grid (t ₁ ,t _2,i ) is The electricity purchase cost of the system in the ith electric energy regulation cycle in the t1th heat and cold energy regulation cycle; F _gas (t ₁ ,t _2,i ) is the _{t 1th} heat and cold energy regulation cycle in the The cost of purchasing natural gas for the system in the i-th electric energy regulation cycle; F _main (t ₁ ,t _2,i ) is the system equipment maintenance in the i-th electric energy regulation cycle in the t ₁ -th heat energy and cooling energy regulation cycle Cost; F _poll (t ₁ , t _{2 , i} ) is the pollution gas emission control cost in the i-th electric energy regulation cycle in the t ₁ -th heat energy and cold energy regulation cycle. 3. The system optimization operation method according to claim 2, wherein the system electricity purchase cost F _grid (t ₁ , t _{2 , i} ) in the minimum objective function of the total system operation 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 ) Among them, Δt ₂ is the time interval of the electric energy regulation cycle; P _grid (t ₁ , t _{2 , i} ) is the power purchase of the system in the ith electric energy regulation cycle in the t ₁ th thermal energy and cold energy regulation cycle; f _grid (t ₁ , t _{2 , i} ) is the real-time electricity price of the power grid in the ith electric energy regulation period in the t ₁ th thermal energy and cold energy regulation period; In the minimum objective function of the total operating cost of the system, the system purchase natural gas cost F _gas (t ₁ , t _{2 , i} ) 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 ) Among them, V _gas (t ₁ , t _{2 , i} ) is the natural gas volume consumed by the system in the i-th electric energy regulation cycle in the t ₁ -th heat and cold energy regulation cycle; Δt ₂ is the time interval of the electric energy regulation cycle; f _gas (t ₁ , t _{2, i} ) is the price of natural gas in the ith electric energy regulation period within the t ₁ th thermal energy and cold energy regulation period; The system equipment maintenance cost F _main (t ₁ ,t _2,i ) in the minimum objective function of the total operating cost of the system is specifically expressed as follows: 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 ) Among them, k _GE [P _GE (t ₁ , t _{2 , i} )] is the maintenance coefficient of the gas-fired internal combustion engine under different output powers in the i-th electric energy regulation cycle in the t ₁ -th thermal energy and cold energy regulation cycle; P _GE (t ₁ ,t _2,i ) is the output electric power of the gas-fired internal combustion engine in the i-th electric energy regulation cycle in the t ₁ -th thermal energy and cold energy regulation cycle; k _AP.cool [Q _AP.cool (t ₁ ,t _2,i )] is the cooling power maintenance coefficient of the flue gas absorption heat pump equipment in the ith electric energy regulation cycle in the t1th heat and cold energy regulation cycle; Q _AP.cool (t ₁ ,t _{2 ,i} ) is the output cooling power of the flue gas absorption heat pump in the ith electric energy regulation cycle in the t1th heat and cold energy regulation cycle; k _AP.heat [Q _AP.heat (t ₁ ,t _{2 ,i} )] is the ith The thermal power maintenance coefficient of the flue gas absorption heat pump equipment in the ith electric energy adjustment cycle in the _t1 heat and cold energy adjustment cycle; Q _AP.heat (t ₁ ,t _2,i ) is the t _1th heat energy and The output heat power of the flue gas absorption heat pump in the i-th electric energy regulation cycle in the cooling energy regulation cycle; k _AC.heat [Q _AC.heat (t ₁ ,t _2,i )] is the t _1th heat energy and cold energy The maintenance coefficient _of the absorption _chiller in the _ith electric energy adjustment cycle in the adjustment cycle _of The thermal power absorbed by the absorption chiller of the electric energy regulation cycle; The pollution gas emission control cost F _poll (t ₁ ,t _2,i ) in the minimum objective function of the total operating cost of the system is specifically expressed as follows:

Among them, t ₁ represents the adjustment period of heat energy and cold energy; t _2,i represents the ith electric energy adjustment period in the t _1th heat energy and cold energy adjustment period; Δt ₂ is the time interval of the electric energy adjustment period; λ is the The number of types of pollutant emissions in the system, including: CO ₂ , SO ₂ , NO _x ; δ _λ is the treatment cost of different emissions including CO ₂ , SO ₂ , NO _x ; α _grid.λ is the difference in power grid power Emission coefficient of emissions; P _grid (t ₁ , t _{2, i} ) is the power purchased by the system and the grid in the i-th electric energy regulation cycle in the t ₁ -th thermal energy and cooling energy regulation cycle; α _GE.λ is Emission coefficient of the electric power of the gas internal combustion engine to different emissions; P _GE (t ₁ ,t _2,i ) is the power generation power of the gas internal combustion engine in the ith electric energy regulation cycle in the t _1th thermal energy and cold energy regulation cycle.

Images