CN115682090A - Optimization control method and device for central heating system - Google Patents

Optimization control method and device for central heating system Download PDF

Info

Publication number
CN115682090A
CN115682090A CN202211108359.1A CN202211108359A CN115682090A CN 115682090 A CN115682090 A CN 115682090A CN 202211108359 A CN202211108359 A CN 202211108359A CN 115682090 A CN115682090 A CN 115682090A
Authority
CN
China
Prior art keywords
heat
pump
follows
mathematical model
constraint
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211108359.1A
Other languages
Chinese (zh)
Inventor
蒋利民
张庆
李�昊
张思瑞
王占博
卜凡鹏
成岭
马美秀
覃剑
郭京超
郭炳庆
张静
林晶怡
李文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Smart Grid Research Institute Co ltd
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Beijing Electric Power Co Ltd
Original Assignee
State Grid Smart Grid Research Institute Co ltd
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Beijing Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Smart Grid Research Institute Co ltd, State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Beijing Electric Power Co Ltd filed Critical State Grid Smart Grid Research Institute Co ltd
Priority to CN202211108359.1A priority Critical patent/CN115682090A/en
Publication of CN115682090A publication Critical patent/CN115682090A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B30/00Energy efficient heating, ventilation or air conditioning [HVAC]
    • Y02B30/70Efficient control or regulation technologies, e.g. for control of refrigerant flow, motor or heating

Landscapes

  • Steam Or Hot-Water Central Heating Systems (AREA)

Abstract

The invention relates to the technical field of control of electric concurrent heating equipment of a heat exchange station, and particularly provides an optimal control method and device of a central heating system, wherein the optimal control method comprises the following steps: solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result; based on the optimization result, obtaining an optimization control scheme of the electric heat compensation equipment in the central heating system; wherein the electric supplementary heating device comprises at least one of: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: heat pump output, gas boiler output, opening of a primary regulating valve and rotating speed of a circulating pump. According to the technical scheme provided by the invention, on one hand, the heating quality and efficiency can be improved, on the other hand, the heat exchange station can be used as a load aggregator to participate in an auxiliary service market and respond to electricity price by controlling the output of the heat pump and adjusting the opening of the valve, so that the heating is ensured, and the economic benefit of the heat exchange station is improved.

Description

Optimization control method and device for central heating system
Technical Field
The invention relates to the technical field of control of electric heat supplementing equipment of a heat exchange station, in particular to an optimal control method and device of a central heating system.
Background
The electric heating is a device for converting clean electric energy into heat energy, the characteristics of energy saving and comfort are more and more accepted by heating users, wherein a heat pump consumes a small amount of electric energy while utilizing the heat of the surrounding environment through a reverse Carnot principle, so that higher temperature is achieved, the heating of the users is realized, the working principle is similar to that of a compressor of a refrigerator, and the heat energy of three same units can be generated when the electric energy of one unit is consumed under the common condition. High efficiency heat pumps are therefore the primary means of heating consumers.
In the urban heating field, municipal administration central heating is the important means of guaranteeing user's heat comfort degree to gas boiler is as the heat source, conveys the heat to each district heat exchange station with the form of heat medium through a pipe network, and the heat exchange station carries out the heat exchange through the heat exchanger, converts the required domestic water of user or heating water into, satisfies the user demand. However, municipal central heating development is relatively lagged behind at present, on one hand, the problems of heat source shortage, high energy consumption and high pollution exist, on the other hand, due to the complexity of a heating pipe network, heat loss is large in the heat transfer process, and in addition, when any module of a heating system breaks down, the heating effect is affected, and the heating quality and efficiency are reduced.
Disclosure of Invention
In order to overcome the defects, the invention provides an optimization control method and device for a central heating system.
In a first aspect, there is provided a method for optimally controlling a central heating system, the method comprising:
solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result;
obtaining an optimized control scheme of the electric heat supplementing equipment in the central heating system based on the optimized result;
wherein the electric supplementary heating device comprises at least one of: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: the output of the heat pump, the output of the gas boiler, the opening of the primary regulating valve and the rotating speed of the circulating pump.
Preferably, the pre-constructed comprehensive benefit optimization model comprises: an objective function that aims at economic cost minimization, and constraints configured for optimal control of the central heating system.
Further, the mathematical model of the objective function is as follows:
min F=C grid +C gas +C maintain +C car -E peak -E allowance
in the above formula, F is the comprehensive benefit of the heating system of a single heat exchange station containing heat pump heat compensation, C grid To purchase electricity cost, C gas For cost of gas, C maintain For operating maintenance costs, C car For carbon emission costs, E peak Earnings to participate in the peak shaving aid service market, E allowance It is a heating patch for clean energy.
Further, the mathematical model of the electricity purchase cost is as follows:
Figure BDA0003842327070000021
the mathematical model of the gas cost is as follows:
Figure BDA0003842327070000022
the mathematical model of the operation and maintenance cost is as follows:
Figure BDA0003842327070000023
the mathematical model of the carbon emission cost is as follows:
Figure BDA0003842327070000024
the mathematical model of the income gained by participating in the peak shaving auxiliary service market is as follows:
Figure BDA0003842327070000025
the mathematical model of the clean energy heating subsidy is as follows:
Figure BDA0003842327070000026
in the above formula, the first and second carbon atoms are,
Figure BDA0003842327070000027
is the electricity price in a unit time period, P wp (t) the power consumption of the circulation pump at time t,
Figure BDA0003842327070000028
the electric power consumed by the ith heat pump at the moment t, N is the number of heat pumps in the heat exchange station, G gb (t) is the gas consumption of the gas boiler at time t, eta gb In order to achieve the conversion efficiency of the gas boiler,
Figure BDA0003842327070000029
is the price of natural gas, k, per unit time period gb Is the operating and maintenance cost coefficient, k, of the gas boiler hp For the coefficient of operating maintenance of the heat pump, k wp For the operating and maintenance cost factor, Q, of the circulating pump gb (t) is the output of the gas boiler at time t, P wp (t) power consumption of the circulation pump at time t, p car Is the cost per carbon emission, delta ec 、δ gb Carbon emission coefficients for electricity purchase and heat supply of a gas boiler are respectively, k is an auxiliary peak regulation time interval specified by a power grid,
Figure BDA0003842327070000031
subsidizing the electricity price for assisting peak shaving in a unit time period, f hp As a function of the total electricity consumption and government subsidy policies of the heat pump system of the heat exchange station, c clean Is a subsidy price.
Further, the constraint condition includes at least one of: user heat load demand constraint, heat pump operation constraint, heat exchange station heat pump cluster total power consumption constraint, gas boiler output constraint, primary regulating valve relative opening constraint, heat medium flow constraint and circulating pump power constraint.
Further, the mathematical model of the user thermal load demand constraint is as follows:
Q hp (t)+Q hes (t)=Q user (t)
θ in (t+1)=θ in (t)+(Q user (t)-Q loss (t))△t/C M M k
Q loss (t)=ξ lossin (t)-θ out (t))
Figure BDA0003842327070000032
in (t+1)-θ in (t)|≤θ ch
Q hes (t)=k hes Q gb (t)
Figure BDA0003842327070000033
in the above formula, Q hp (t) the heat-supplementing power of the heat pump system at time t, Q hes (t) heat exchange station exchanges heat through heat exchanger at t moment, and provides heat power Q for user user (t) is the thermal load demand of the user at time t,. DELTA.t is the time interval,. Theta. in (t) indoor temperature of the building at time t, Q loss (t) Heat losses associated with the physical structure of the building and the outdoor environment at time t, C M Specific heat capacity of indoor air, M k Is the indoor air quality xi loss To and from a buildingCoefficient relating to physical structure and outdoor environment, theta out (t) is the outdoor temperature of the building at time t,
Figure BDA0003842327070000034
the lower limit of the room temperature is,
Figure BDA0003842327070000035
at the upper limit of room temperature, [ theta ] ch Is the maximum fluctuation range of the room temperature within a unit time period, k hes For the heat exchange efficiency of heat exchangers, Q gb (t) is the heating power of the gas boiler at time t,
Figure BDA0003842327070000036
and N is the number of heat pumps in the heat exchange station, and COP is the coefficient of performance of the heat pumps.
Further, the mathematical model of the heat pump operation constraint is as follows:
Figure BDA0003842327070000037
in the above formula, the first and second carbon atoms are,
Figure BDA0003842327070000038
is the rated power of the heat pump operation.
Further, the mathematical model of the heat exchange station heat pump cluster total consumed electric power constraint is as follows:
Figure BDA0003842327070000041
the mathematical model of the output restriction of the gas boiler is as follows:
Q gb.min ≤Q gb (t)≤Q gb.max
in the above formula, the first and second carbon atoms are,
Figure BDA0003842327070000042
and
Figure BDA0003842327070000043
minimum and maximum allowable power consumption power, Q, of a node j connected with a heat pump cluster of the heat exchange station at the moment t gb.min And Q gb.max The upper and lower output limits of the gas boiler.
Further, the mathematical model of the constraint of the relative opening of the primary regulating valve is as follows:
Figure BDA0003842327070000044
0.1≤K V ≤0.9
Figure BDA0003842327070000045
Figure BDA0003842327070000046
in the above formula, K V Relative opening of primary valve, /) V For the stroke of the valve core at a certain opening of the regulating valve,/ V.max The valve core stroke when the regulating valve is fully opened, V is the heat medium flow in the pipeline, R V For adjustable ratio of regulating valves, V max Maximum flow rate of heating medium, V, controlled by regulating valve min The minimum flow of the heating medium which can be controlled by the regulating valve.
Further, the mathematical model of the heat medium flow constraint is as follows:
Q gb (t)=ρ×c p ×V(t)×(T s -T h )
V min ≤V(t)≤V max
in the above formula, c p Is the specific heat capacity of the primary heat supply network heat medium, rho is the primary heat supply network heat medium density, V (T) is the heat medium flow in the pipeline at the time T, T s For temperature of the water supply, T h The temperature of the return water is shown.
Further, the mathematical model of the circulation pump power constraint is as follows:
Figure BDA0003842327070000047
0≤P wp ≤P wp.max
in the above formula, V e To design the operating flow, n e In order to design the running rotating speed of the water pump, n is the actual running rotating speed of the water pump,
Figure BDA0003842327070000051
for designing the power consumption of the water pump, P wp For actual operation of the pump consuming power, P wp.max The maximum power consumption of the circulating pump.
In a second aspect, there is provided an optimization control device for a central heating system, including:
the first analysis module is used for solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result;
the second analysis module is used for obtaining an optimized control scheme of the electric heat compensation equipment in the central heating system based on the optimization result;
wherein the electric supplementary heating device comprises at least one of the following: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: the output of the heat pump, the output of the gas boiler, the opening of the primary regulating valve and the rotating speed of the circulating pump.
Preferably, the pre-constructed comprehensive benefit optimization model comprises: an objective function that aims at economic cost minimization, and constraints configured for optimal control of the central heating system.
Further, the mathematical model of the objective function is as follows:
min F=C grid +C gas +C maintain +C car -E peak -E allowance
in the above formula, F is the comprehensive benefit of the heating system of a single heat exchange station containing heat pump heat compensation, C grid To purchase electricity cost, C gas As a gasCost, C maintain For operating maintenance costs, C car For carbon emission costs, E peak Gain to participate in the peak shaving aid service market, E allowance It is a supplementary heating paste for clean energy.
Further, the mathematical model of the electricity purchasing cost is as follows:
Figure BDA0003842327070000052
the mathematical model of the gas cost is as follows:
Figure BDA0003842327070000053
the mathematical model of the operation and maintenance cost is as follows:
Figure BDA0003842327070000054
the mathematical model of the carbon emission cost is as follows:
Figure BDA0003842327070000061
the mathematical model of the income obtained by participating in the peak shaving auxiliary service market is as follows:
Figure BDA0003842327070000062
the mathematical model of the clean energy heating subsidy is as follows:
Figure BDA0003842327070000063
in the above formula, the first and second carbon atoms are,
Figure BDA0003842327070000064
is the electricity price in a unit time period, P wp (t) is the power consumption of the circulation pump at time t,
Figure BDA0003842327070000065
the electric power consumed by the ith heat pump at the moment t, N is the number of heat pumps in the heat exchange station, G gb (t) is the gas consumption of the gas boiler at time t, eta gb In order to achieve the conversion efficiency of the gas boiler,
Figure BDA0003842327070000066
is the price of natural gas, k, per unit time period gb Is the operating and maintenance cost coefficient, k, of the gas boiler hp For the operating and maintenance cost coefficient, k, of the heat pump wp For operating and maintenance cost factor, Q, of the circulating pump gb (t) is the output of the gas boiler at time t, P wp (t) power consumption of the circulation pump at time t, p car Is the cost per carbon emission, delta ec 、δ gb Carbon emission coefficients for electricity purchase and heat supply of a gas boiler are respectively, k is an auxiliary peak regulation time interval specified by a power grid,
Figure BDA0003842327070000067
subsidizing the electricity price for assisting peak shaving in a unit time period, f hp Function of heat pump system of heat exchange station with respect to total electricity consumption and government subsidy policy, c clean Is a subsidy price.
Further, the constraint condition includes at least one of: the method comprises the following steps of user heat load demand constraint, heat pump operation constraint, heat exchange station heat pump cluster total power consumption constraint, gas boiler output constraint, primary regulating valve relative opening constraint, heat medium flow constraint and circulating pump power constraint.
Further, the mathematical model of the user thermal load demand constraint is as follows:
Q hp (t)+Q hes (t)=Q user (t)
θ in (t+1)=θ in (t)+(Q user (t)-Q loss (t))△t/C M M k
Q loss (t)=ξ lossin (t)-θ out (t))
Figure BDA0003842327070000068
in (t+1)-θ in (t)|≤θ ch
Q hes (t)=k hes Q gb (t)
Figure BDA0003842327070000071
in the above formula, Q hp (t) the heat-supplementing power of the heat pump system at time t, Q hes (t) heat exchange station exchanges heat through heat exchanger at t moment, and provides heat power Q for user user (t) is the thermal load demand of the user at time t,. DELTA.t is the time interval,. Theta. in (t) indoor temperature of the building at time t, Q loss (t) Heat losses associated with the physical structure of the building and the outdoor environment at time t, C M Specific heat capacity of indoor air, M k Is indoor air quality xi loss Is a coefficient relating to the physical structure of the building and the outdoor environment, theta out (t) is the outdoor temperature of the building at time t,
Figure BDA0003842327070000072
the lower limit of the room temperature is,
Figure BDA0003842327070000073
at the upper limit of room temperature, [ theta ] ch Is the maximum fluctuation range of the room temperature within a unit time period, k hes For the heat exchange efficiency of heat exchangers, Q gb (t) is the heating power of the gas boiler at time t,
Figure BDA0003842327070000074
and N is the number of heat pumps in the heat exchange station, and COP is the coefficient of performance of the heat pumps.
Further, the mathematical model of the heat pump operation constraint is as follows:
Figure BDA0003842327070000075
in the above-mentioned formula, the compound has the following structure,
Figure BDA0003842327070000076
is the rated power of the heat pump operation.
Further, the mathematical model of the heat exchange station heat pump cluster total consumed electric power constraint is as follows:
Figure BDA0003842327070000077
the mathematical model of the output restriction of the gas boiler is as follows:
Q gb.min ≤Q gb (t)≤Q gb.max
in the above-mentioned formula, the compound has the following structure,
Figure BDA0003842327070000078
and
Figure BDA0003842327070000079
minimum and maximum allowable power consumption power, Q, of a node j connected with a heat pump cluster of the heat exchange station at the moment t gb.min And Q gb.max The upper and lower output limits of the gas boiler.
Further, the mathematical model of the constraint of the relative opening of the primary regulating valve is as follows:
Figure BDA00038423270700000710
0.1≤K V ≤0.9
Figure BDA00038423270700000711
Figure BDA0003842327070000081
in the above formula, K V Is the relative opening of the primary valve, /) V Is the valve core stroke when the regulating valve has a certain opening degree V.max The valve core stroke when the regulating valve is fully opened, V is the heat medium flow in the pipeline, R V For regulating the adjustable ratio of valves, V max Maximum flow rate of heating medium, V, controlled by regulating valve min The minimum flow of the heating medium which can be controlled by the regulating valve.
Further, the mathematical model of the heat medium flow constraint is as follows:
Q gb (t)=ρ×c p ×V(t)×(T s -T h )
V min ≤V(t)≤V max
in the above formula, c p Is the specific heat capacity of the primary heat supply network heat medium, rho is the primary heat supply network heat medium density, V (T) is the heat medium flow in the pipeline at the time T, T s For supply water temperature, T h The temperature of the return water is shown.
Further, the mathematical model of the circulation pump power constraint is as follows:
Figure BDA0003842327070000082
0≤P wp ≤P wp.max
in the above formula, V e To design the operating flow, n e In order to design the running rotating speed of the water pump, n is the actual running rotating speed of the water pump,
Figure BDA0003842327070000083
for designing the power consumption of the water pump, P wp For actual operation of the pump consuming power, P wp.max The maximum power consumption of the circulating pump.
In a third aspect, a computer device is provided, comprising: one or more processors;
the processor to store one or more programs;
the one or more programs, when executed by the one or more processors, implement the method for optimized control of a central heating system.
In a fourth aspect, a computer-readable storage medium is provided, on which a computer program is stored, which, when executed, implements the method for optimizing a central heating system.
One or more technical schemes of the invention at least have one or more of the following beneficial effects:
the invention provides an optimization control method and device of a central heating system, which comprises the following steps: solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result; based on the optimization result, obtaining an optimization control scheme of the electric heat compensation equipment in the central heating system; wherein the electric supplementary heating device comprises at least one of: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: the output of the heat pump, the output of the gas boiler, the opening of the primary regulating valve and the rotating speed of the circulating pump. According to the technical scheme provided by the invention, the pre-constructed comprehensive benefit optimization model considers the self constraints of the heat pump and the gas boiler on one hand and considers the constraints of a primary regulating valve and a circulating pump on the other hand, so that the actual operation condition is met, the heating requirement is met, and the comprehensive benefit optimization model has certain significance for large-scale popularization;
furthermore, the technical scheme provided by the invention not only improves the economic operation of the electric heat supplementing equipment, but also responds to the peak-to-valley electrovalence signal of the power grid, and reduces the peak-to-valley difference of the power grid to a certain extent.
Drawings
FIG. 1 is a flow chart illustrating the main steps of an optimization control method for a central heating system according to an embodiment of the present invention;
fig. 2 is a diagram of a cooperative control system of a heat pump device and a heat supply pipe network of an urban heat exchange station according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a cooperative optimization control strategy of the heat pump equipment and the heat supply pipe network of the urban heat exchange station according to the embodiment of the invention;
fig. 4 is an overall architecture diagram of the cooperative optimization control device for heat pump equipment and a heat supply pipe network of the urban heat exchange station according to the embodiment of the invention;
fig. 5 is a block diagram of the main configuration of an optimization controller of a central heating system according to an embodiment of the present invention.
Detailed Description
The following provides a more detailed description of embodiments of the present invention, with reference to the accompanying drawings.
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.
As disclosed in the background art, the characteristics of energy saving and comfort of electric heating are more and more accepted by heating users as a device for converting clean electric energy into heat energy, wherein a heat pump uses ambient heat and consumes a small amount of electric energy through a reverse carnot principle, thereby achieving a higher temperature and realizing heating of users. High efficiency heat pumps are therefore the primary means of heating consumers.
In the urban heating field, municipal central heating is the important means of guaranteeing user's heat comfort degree to gas boiler is as the heat source, conveys the heat to each district heat exchange station with the form of heat medium through a pipe network, and the heat exchange station carries out heat exchange through the heat exchanger, converts the required domestic water of user or heating water into, satisfies the user demand. However, municipal central heating development is relatively lagged behind at present, on one hand, the problems of heat source shortage, high energy consumption and high pollution exist, on the other hand, due to the complexity of a heating pipe network, heat loss is large in the heat transfer process, and in addition, when any module of a heating system breaks down, the heating effect is affected, and the heating quality and efficiency are reduced.
In order to solve the above problems, the present invention provides an optimization control method and device for a central heating system, including: solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result; based on the optimization result, obtaining an optimization control scheme of the electric heat compensation equipment in the central heating system; wherein the electric supplementary heating device comprises at least one of the following: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: the output of the heat pump, the output of the gas boiler, the opening of the primary regulating valve and the rotating speed of the circulating pump. The pre-constructed comprehensive benefit optimization model in the technical scheme provided by the invention considers the self constraints of the heat pump and the gas boiler on one hand and the constraints of a primary regulating valve and a circulating pump on the other hand, so that the heat supply network operation is in line with the actual operation condition, the heating requirement is met, and the large-scale popularization of the heat supply network is of certain significance. The above scheme is explained in detail below.
Example 1
Referring to fig. 1, fig. 1 is a flow chart illustrating main steps of an optimization control method for a central heating system according to an embodiment of the present invention. As shown in fig. 1, the optimization control method for a central heating system in the embodiment of the present invention mainly includes the following steps:
step S101: solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result;
step S102: based on the optimization result, obtaining an optimization control scheme of the electric heat compensation equipment in the central heating system;
wherein the electric supplementary heating device comprises at least one of the following: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: the output of the heat pump, the output of the gas boiler, the opening of the primary regulating valve and the rotating speed of the circulating pump.
In this embodiment, the pre-constructed comprehensive benefit optimization model includes: an objective function that aims at economic cost minimization, and constraints configured for optimal control of the central heating system.
In one embodiment, the heat pump responds to the peak-valley electricity price of the power grid and gains obtained by participating in peak shaving of the auxiliary service market, on the other hand, the running cost of the heat pump, the gas consumption and running maintenance cost of the gas boiler, the carbon emission cost, the running cost of the heat grid and other costs are considered, the economic cost is the minimum as an objective function, the heat load balance, the heat pump power, the gas boiler output, the opening of the primary regulating valve and the circulating pump power are taken as constraints, and the heat pump, the gas boiler output, the opening of the primary regulating valve and the rotating speed of the circulating pump are regulated through an optimization algorithm, so that the optimal running state is achieved. Thus, the mathematical model of the objective function is as follows:
min F=C grid +C gas +C maintain +C car -E peak -E allowance
in the above formula, F is the comprehensive benefit of the heating system of a single heat exchange station containing heat pump heat compensation, C grid To purchase electricity cost, C gas For cost of gas, C maintain For operating maintenance costs, C car For carbon emission costs, E peak Earnings to participate in the peak shaving aid service market, E allowance It is a supplementary heating paste for clean energy.
Wherein, the mathematical model of the electricity purchasing cost is as follows:
Figure BDA0003842327070000111
the mathematical model of the gas cost is as follows:
Figure BDA0003842327070000112
the mathematical model of the operation and maintenance cost is as follows:
Figure BDA0003842327070000113
the mathematical model of the carbon emission cost is as follows:
Figure BDA0003842327070000114
the mathematical model of the income gained by participating in the peak shaving auxiliary service market is as follows:
Figure BDA0003842327070000115
the mathematical model of the clean energy heating subsidy is as follows:
Figure BDA0003842327070000116
in the above-mentioned formula, the compound has the following structure,
Figure BDA0003842327070000117
is the electricity price in a unit time period, P wp (t) the power consumption of the circulation pump at time t,
Figure BDA0003842327070000118
the electric power consumed by the ith heat pump at the moment t, N is the number of heat pumps in the heat exchange station, G gb (t) is the gas consumption of the gas boiler at time t, eta gb In order to achieve the conversion efficiency of the gas boiler,
Figure BDA0003842327070000119
is the natural gas price, k, over a unit period of time gb For the operating maintenance cost coefficient, k, of the gas boiler hp For the coefficient of operating maintenance of the heat pump, k wp For the operating and maintenance cost factor, Q, of the circulating pump gb (t) is the output of the gas boiler at time t, P wp (t) power consumption of the circulation pump at time t, p car Is the cost per carbon emission, delta ec 、δ gb Carbon emission coefficients for electricity purchase and heat supply of a gas boiler are respectively, k is an auxiliary peak regulation period specified by a power grid,
Figure BDA00038423270700001110
subsidizing the electricity price for assisting peak shaving in a unit time period, f hp Function of heat pump system of heat exchange station with respect to total electricity consumption and government subsidy policy, c clean Is a subsidy price.
Further, the constraint condition includes at least one of: the method comprises the following steps of user heat load demand constraint, heat pump operation constraint, heat exchange station heat pump cluster total power consumption constraint, gas boiler output constraint, primary regulating valve relative opening constraint, heat medium flow constraint and circulating pump power constraint.
In one embodiment, the user thermal load demand constraint is specifically:
the heat pump is as electric concurrent heating equipment, and the heat that its heat of release and heat supply network pass through the heat exchanger exchange jointly supplies heat to the user, satisfies user's heat load demand:
Q hp (t)+Q hes (t)=Q user (t)
for the user heat load demand, assuming that the user indoor temperature is uniformly distributed, and the coefficient related to the physical structure of the building is a fixed value, the dynamic change process of the user indoor temperature can be expressed as:
θ in (t+1)=θ in (t)+(Q user (t)-Q loss (t))△t/C M M k
Q loss (t)=ξ lossin (t)-θ out (t))
the thermal comfort of the user determines the range of room temperature fluctuations of the building allowed by the heating user. The room temperature is not lower than the lowest heating temperature in the national heating season, otherwise, the requirement of the physiological comfort degree of a heat user cannot be met; in addition to not being higher than the physiological capacity of the human body, the room temperature is also limited by the maximum adjustable level of the system. Within any scheduling period, the room temperature fluctuation can be expressed as:
Figure BDA0003842327070000121
in (t+1)-θ in (t)|≤θ ch
for a heat exchange station, the supply of thermal power can be expressed as:
Q hes (t)=k hes Q gb (t)
for a heat pump, the heating power can be expressed as:
Figure BDA0003842327070000122
in the above formula, Q hp (t) the heat-supplementing power of the heat pump system at time t, Q hes (t) heat exchange station exchanges heat through heat exchanger at t moment, and provides heat power Q for user user (t) is the thermal load demand of the user at time t,. DELTA.t is the time interval,. Theta. in (t) indoor temperature of the building at time t, Q loss (t) Heat losses associated with the physical structure of the building and the outdoor environment at time t, C M Specific heat capacity of indoor air, M k Is the indoor air quality xi loss Is a coefficient related to the physical structure of the building and the outdoor environment, θ out (t) is the outdoor temperature of the building at time t,
Figure BDA0003842327070000123
the lower limit of the room temperature is,
Figure BDA0003842327070000124
at the upper limit of room temperature, θ ch Is the maximum fluctuation range of the room temperature within a unit time period, k hes For the heat exchange efficiency of the heat exchanger, in actual operation, the gas-fired boiler is used as a total heat source, heat is transferred to each heat exchange station through a primary heat supply network, in order to simplify the model, the heat exchange station containing heat pump heat compensation is supposed to be specially supplied by the single gas-fired boiler, and Q is realized gb (t) is the heating power of the gas boiler at the time t,
Figure BDA0003842327070000131
and N is the number of heat pumps in the heat exchange station, and COP is the coefficient of performance of the heat pumps, and the COP is the same as the model number of each heat pump in the heat exchange station.
In one embodiment, the performance of the heat pump unit is closely related to the power of the compressor and the water pump during operation, if the operation power of the heat pump unit exceeds the rated power, the compressor and the water pump can run in an overload mode, and the compression pump mechanism and the driving motor can be damaged in serious conditions; when the operating power is too low, the service life of the heat pump units can be reduced and energy is wasted, so that in order to achieve a better operating effect, the performance of the heat pump units needs to be constrained, the actual operating power of each heat pump unit cannot be less than 0.25 time of the rated power, and the mathematical model of the heat pump operation constraint is as follows:
Figure BDA0003842327070000132
in the above formula, the first and second carbon atoms are,
Figure BDA0003842327070000133
is the rated power of the heat pump operation.
In one embodiment, the heat pump cluster is connected to a power distribution network, and in order to ensure the stability of a power system, the voltage of a node connected with the heat pump cluster and the power of a branch circuit where the node is located need to meet the requirement of the power grid, so that a mathematical model of the total electric power consumption constraint of the heat exchange station heat pump cluster is as follows:
Figure BDA0003842327070000134
the mathematical model of the output restriction of the gas boiler is as follows:
Q gb.min ≤Q gb (t)≤Q gb.max
in the above formula, the first and second carbon atoms are,
Figure BDA0003842327070000135
and
Figure BDA0003842327070000136
minimum and maximum allowable power consumption power, Q, of a node j connected with a heat pump cluster of the heat exchange station at the moment t gb.min And Q gb.max The upper and lower output limits of the gas boiler.
In one embodiment, the thermodynamic system has 2 primary modes of regulation to meet the customer heat load demand, namely a mass regulation to change the temperature of the heating medium in the system while maintaining the system flow constant and a mass regulation to change the system heating medium flow while maintaining the system temperature constant. Because the quantity regulation is adopted, the heating speed is very high, and the heat delay and the heat loss are not generated, so the quantity regulation is adopted. When the heating system is regulated, the primary regulating valve plays an important role, and the flow passing through the heating system is controlled by regulating the relative opening degree of the valve, so that the mathematical model of the constraint of the relative opening degree of the primary regulating valve is as follows:
Figure BDA0003842327070000141
when the primary regulating valve works, in order to prevent the influence on the regulating performance and the economical efficiency caused by the larger caliber selection of the valve, the maximum opening degree of the primary regulating valve is generally not more than 90 percent; in order to prevent the valve core from being damaged and the valve flow characteristic from being deteriorated and even malfunctioning due to the serious erosion of fluid to the valve core and the valve seat when the valve core is in a small opening degree, the minimum opening degree is generally not less than 10 percent, and the expression is as follows:
0.1≤K V ≤0.9
according to the difference of flow characteristics, the valves are divided into a linear flow characteristic, an equal percentage flow characteristic and a quick opening flow characteristic, and four parabolic characteristic regulating valves are adopted. The expression is as follows:
Figure BDA0003842327070000142
Figure BDA0003842327070000143
in the above formula, K V Relative opening of primary valve, /) V Is the valve core stroke when the regulating valve has a certain opening degree V.max Is the valve core stroke when the regulating valve is fully opened, V is the heat medium flow in the pipeline, R V For adjustable ratio of regulating valves, V max Maximum flow rate of heating medium, V, controlled by regulating valve min The minimum flow of the heating medium which can be controlled by the regulating valve.
In one embodiment, the heat generated at the heat source is transmitted to each heat exchange station through the heat medium, and in order to avoid imbalance caused by too small or too large flow of the heat medium in the pipeline, the flow of the heat medium needs to be controlled within a reasonable range, so that the mathematical model of the heat medium flow constraint is as follows:
Q gb (t)=ρ×c p ×V(t)×(T s -T h )
V min ≤V(t)≤V max
in the above formula, c p Is the specific heat capacity of the primary heat supply network heat medium, rho is the primary heat supply network heat medium density, V (T) is the heat medium flow in the pipeline at the time T, T s For temperature of the water supply, T h The temperature of the return water is shown.
In one embodiment, the circulating pump is used as an important device in the central heating pipe network, the rotating speed of the circulating pump is adjusted by consuming a part of electric energy, pressure difference is generated in the water supply pipe network and the water return pipe network respectively, and heat medium is pushed to carry out heat circulation in the heating system. In order to prevent the overflow problem of the water pump, the maximum output limit of the circulating pump is met during operation, and therefore, the power constraint mathematical model of the circulating pump is as follows:
Figure BDA0003842327070000151
0≤P wp ≤P wp.max
in the above formula, V e To design the operating flow, n e In order to design the running rotating speed of the water pump, n is the actual running rotating speed of the water pump,
Figure BDA0003842327070000152
for designing the power consumption of the water pump, P wp For actual operation of the pump consuming power, P wp.max The maximum power consumption of the circulating pump.
In a specific embodiment, the solution method is a simplex method in linear optimization:
considering that the optimization target and the constraint condition of the comprehensive benefit function model of the heating system of the single heat exchange station containing the heat pump for heat compensation are linear equations, the invention adopts a simplex method to solve, and converts the benefit objective function into a standard form:
Figure BDA0003842327070000153
in the formula: optimizing variable x j Is an n-dimensional variable comprising the output of a heat pump and a gas boiler; z is a model optimization objective function; a is a ij Coefficients of various constraint conditions satisfied by the heat pump and the gas boiler (both the constraint of the primary regulating valve and the constraint of the circulating pump can be converted into the output constraint of the gas boiler); max is a function of the maximum value.
The solving process mainly comprises 3 steps:
(1) Finding an initial basic feasible solution;
(2) And judging whether the solution is the optimal solution or not, if so, ending the solution. Otherwise, turning to (3);
(3) And (3) searching a new basic feasible solution, and repeating the step (2).
In one application scenario, as shown in fig. 2, the control system includes a heat pump system, a large power grid, a gas boiler, a primary regulating valve, a circulation pump, and a user heat load. The method comprises the steps of predicting the outdoor environment temperature to release the heat load demand of a user, taking one day as an operation working period, taking the optimal economy as an objective function, calculating and controlling the heat pump, the output of the gas boiler and the opening of a primary regulating valve under the constraint conditions of meeting the heat load demand of the user, the heat pump operation, the output of the gas boiler, the relative opening of the primary regulating valve, the heat medium flow, the power of a circulating pump and the like, and solving by adopting a simplex method in linear optimization, thereby realizing an optimal economic working mode.
Specifically, fig. 3 is a schematic diagram of a collaborative optimization control strategy of heat pump equipment and a heat supply pipe network of the urban heat exchange station according to the present invention, and as shown in fig. 3, the output Q of each heat pump is initialized by introducing electricity price, heat price, peak regulation period, equipment parameters, and the like into the program hp And the output Q of the gas boiler gb Finding a basic feasible solution according to the constraint condition; judging whether the solution is the optimal solution or not, if so, outputting the running state of each device; otherwise, the next basic feasible solution is searched by adjusting the output of the heat pump and the gas boiler, the opening of the primary valve and the rotating speed of the circulating pump until the optimal solution is found. The state of the optimal solution is the optimal economic operation mode.
Fig. 4 is an overall architecture diagram of the cooperative optimization control device for the heat pump equipment and the heat supply pipe network of the urban heat exchange station provided by the invention. As shown in fig. 4, a power grid signal, such as a node voltage and a branch power, is acquired through a sensor; heat supply network signals such as primary network heat medium flow, supply and return water temperature, primary regulating valve opening, circulating pump rotating speed and the like; heating plant signals, such as heating power of gas boilers and heat pumps; the thermal load requirements of the user, such as the ambient temperature and room temperature. Through measurement and calculation, the control strategy of fig. 2 is implemented, and the operation states of the heat pump and the heat supply network are controlled, namely the output of the heat pump and the gas boiler, the opening of the primary valve and the rotating speed of the circulating pump are adjusted, so that the optimal economic operation mode is achieved.
Example 2
Based on the same inventive concept, the present invention further provides an optimization control device for a central heating system, as shown in fig. 5, the optimization control device for a central heating system includes:
the first analysis module is used for solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result;
the second analysis module is used for obtaining an optimization control scheme of the electric heat compensation equipment in the central heating system based on the optimization result;
wherein the electric supplementary heating device comprises at least one of: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: heat pump output, gas boiler output, opening of a primary regulating valve and rotating speed of a circulating pump.
Preferably, the pre-constructed comprehensive benefit optimization model comprises: an objective function that aims at economic cost minimization, and constraints configured for optimal control of the central heating system.
Further, the mathematical model of the objective function is as follows:
min F=C grid +C gas +C maintain +C car -E peak -E allowance
in the above formula, F is the comprehensive benefit of the heating system of a single heat exchange station containing heat pump heat compensation, C grid To purchase electricity cost, C gas For cost of gas, C maintain For operating maintenance costs, C car For carbon emission costs, E peak Earnings to participate in the peak shaving aid service market, E allowance It is a supplementary heating paste for clean energy.
Further, the mathematical model of the electricity purchasing cost is as follows:
Figure BDA0003842327070000161
the mathematical model of the gas cost is as follows:
Figure BDA0003842327070000171
the mathematical model of the operation and maintenance cost is as follows:
Figure BDA0003842327070000172
the mathematical model of the carbon emission cost is as follows:
Figure BDA0003842327070000173
the mathematical model of the income gained by participating in the peak shaving auxiliary service market is as follows:
Figure BDA0003842327070000174
the mathematical model of the clean energy heating subsidy is as follows:
Figure BDA0003842327070000175
in the above formula, the first and second carbon atoms are,
Figure BDA0003842327070000176
is the electricity price per unit time period, P wp (t) is the power consumption of the circulation pump at time t,
Figure BDA0003842327070000177
the electric power consumed by the ith heat pump at the moment t, N is the number of heat pumps in the heat exchange station, G gb (t) is the gas consumption of the gas boiler at time t, eta gb In order to achieve the conversion efficiency of the gas boiler,
Figure BDA0003842327070000178
is the natural gas price, k, over a unit period of time gb For the operating maintenance cost coefficient, k, of the gas boiler hp For the operating and maintenance cost coefficient, k, of the heat pump wp For the operating and maintenance cost factor, Q, of the circulating pump gb (t) is the output of the gas boiler at time t, P wp (t) power consumption of the circulation pump at time t, p car Is the cost per carbon emission, delta ec 、δ gb Carbon emission coefficients for electricity purchase and heat supply of a gas boiler are respectively, k is an auxiliary peak regulation period specified by a power grid,
Figure BDA0003842327070000179
subsidizing the electricity price for assisting peak shaving in a unit time period, f hp As a function of the total electricity consumption and government subsidy policies of the heat pump system of the heat exchange station, c clean Is a subsidy price.
Further, the constraint condition includes at least one of: the method comprises the following steps of user heat load demand constraint, heat pump operation constraint, heat exchange station heat pump cluster total power consumption constraint, gas boiler output constraint, primary regulating valve relative opening constraint, heat medium flow constraint and circulating pump power constraint.
Further, the mathematical model of the user thermal load demand constraint is as follows:
Q hp (t)+Q hes (t)=Q user (t)
θ in (t+1)=θ in (t)+(Q user (t)-Q loss (t))△t/C M M k
Q loss (t)=ξ lossin (t)-θ out (t))
Figure BDA0003842327070000181
in (t+1)-θ in (t)|≤θ ch
Q hes (t)=k hes Q gb (t)
Figure BDA0003842327070000182
in the above formula, Q hp (t) the heat-supplementing power of the heat pump system at time t, Q hes (t) heat exchange station exchanges heat through heat exchanger at t moment, and provides heat power Q for user user (t) is the thermal load demand of the user at time t,. DELTA.t is the time interval,. Theta. in (t) indoor temperature of building at time t, Q loss (t) Heat losses associated with the physical structure of the building and the outdoor environment at time t, C M Specific heat capacity of indoor air, M k Is indoor air quality xi loss Is a coefficient related to the physical structure of the building and the outdoor environment, θ out (t) is the outdoor temperature of the building at time t,
Figure BDA0003842327070000183
the lower limit of the room temperature is,
Figure BDA0003842327070000184
at the upper limit of room temperature, θ ch Is the maximum fluctuation range of the room temperature within a unit time period, k hes For the heat exchange efficiency of heat exchangers, Q gb (t) is the heating power of the gas boiler at time t,
Figure BDA0003842327070000185
and N is the number of heat pumps in the heat exchange station, and COP is the coefficient of performance of the heat pumps.
Further, the mathematical model of the heat pump operation constraint is as follows:
Figure BDA0003842327070000186
in the above formula, the first and second carbon atoms are,
Figure BDA0003842327070000187
is the rated power for the heat pump to operate.
Further, the mathematical model of the heat exchange station heat pump cluster total consumed electric power constraint is as follows:
Figure BDA0003842327070000188
the mathematical model of the output restriction of the gas boiler is as follows:
Q gb.min ≤Q gb (t)≤Q gb.max
in the above-mentioned formula, the compound has the following structure,
Figure BDA0003842327070000189
and
Figure BDA00038423270700001810
minimum and maximum allowable power consumption power, Q, of a node j connected with a heat pump cluster of the heat exchange station at the moment t gb.min And Q gb.max The upper and lower output limits of the gas boiler.
Further, the mathematical model of the constraint of the relative opening of the primary regulating valve is as follows:
Figure BDA0003842327070000191
0.1≤K V ≤0.9
Figure BDA0003842327070000192
Figure BDA0003842327070000193
in the above formula, K V Is the relative opening of the primary valve, /) V For the stroke of the valve core at a certain opening of the regulating valve,/ V.max Is the valve core stroke when the regulating valve is fully opened, V is the heat medium flow in the pipeline, R V For adjustable ratio of regulating valves, V max Maximum flow rate of heating medium, V, controlled by regulating valve min The minimum flow of the heating medium which can be controlled by the regulating valve.
Further, the mathematical model of the heat medium flow constraint is as follows:
Q gb (t)=ρ×c p ×V(t)×(T s -T h )
V min ≤V(t)≤V max
in the above formula, c p Is the specific heat capacity of the primary heat supply network heat medium, rho is the primary heat supply network heat medium density, V (T) is the heat medium flow in the pipeline at the time T, T s For temperature of the water supply, T h The temperature of the return water is shown.
Further, the mathematical model of the circulation pump power constraint is as follows:
Figure BDA0003842327070000194
0≤P wp ≤P wp.max
in the above formula, V e To design the operating flow, n e In order to design the running rotating speed of the water pump, n is the actual running rotating speed of the water pump,
Figure BDA0003842327070000195
for designing the power consumption of the water pump, P wp For actual operation of the pump consuming power, P wp.max The maximum power consumption of the circulating pump.
Example 3
Based on the same inventive concept, the present invention also provides a computer apparatus comprising a processor and a memory, the memory being configured to store a computer program comprising program instructions, the processor being configured to execute the program instructions stored by the computer storage medium. The Processor may be a Central Processing Unit (CPU), and may also be other general purpose processors, digital Signal Processors (DSPs), application Specific Integrated Circuits (ASICs), field-Programmable gate arrays (FPGAs) or other Programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc., which are a computing core and a control core of the terminal, and are specifically adapted to load and execute one or more instructions in a computer storage medium so as to implement a corresponding method flow or a corresponding function, so as to implement the steps of the optimization control method of the Central heating system in the foregoing embodiments.
Example 4
Based on the same inventive concept, the present invention further provides a storage medium, in particular, a computer-readable storage medium (Memory), which is a Memory device in a computer device and is used for storing programs and data. It is understood that the computer readable storage medium herein can include both built-in storage media in the computer device and, of course, extended storage media supported by the computer device. The computer-readable storage medium provides a storage space storing an operating system of the terminal. Also, one or more instructions, which may be one or more computer programs (including program code), are stored in the memory space and are adapted to be loaded and executed by the processor. It should be noted that the computer readable storage medium may be a high-speed RAM memory, or a non-volatile memory (non-volatile memory), such as at least one disk memory. One or more instructions stored in the computer readable storage medium may be loaded and executed by a processor to implement the steps of the optimization control method for central heating system in the above embodiments.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: although the present invention has been described in detail with reference to the above embodiments, it should be understood by those skilled in the art that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.

Claims (24)

1. A method for optimizing control of a central heating system, the method comprising:
solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result;
obtaining an optimized control scheme of the electric heat supplementing equipment in the central heating system based on the optimized result;
wherein the electric supplementary heating device comprises at least one of: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: heat pump output, gas boiler output, opening of a primary regulating valve and rotating speed of a circulating pump.
2. The method of claim 1, wherein the pre-constructed synthetic benefit optimization model comprises: an objective function that aims at economic cost minimization, and constraints configured for optimal control of the central heating system.
3. The method of claim 2, wherein the mathematical model of the objective function is as follows:
minF=C grid +C gas +C maintain +C car -E peak -E allowance
in the above formula, F is the comprehensive benefit of the heating system of a single heat exchange station containing heat pump heat compensation, C grid To purchase electricity cost, C gas For cost of gas, C maintain For operating maintenance costs, C car For carbon emission costs, E peak Earnings to participate in the peak shaving aid service market, E allowance It is a supplementary heating paste for clean energy.
4. The method of claim 3, wherein the mathematical model of the electricity purchase cost is as follows:
Figure FDA0003842327060000011
the mathematical model of the gas cost is as follows:
Figure FDA0003842327060000012
the mathematical model of the operation and maintenance cost is as follows:
Figure FDA0003842327060000013
the mathematical model of the carbon emission cost is as follows:
Figure FDA0003842327060000014
the mathematical model of the income obtained by participating in the peak shaving auxiliary service market is as follows:
Figure FDA0003842327060000015
the mathematical model of the clean energy heating subsidy is as follows:
Figure FDA0003842327060000021
in the above formula, the first and second carbon atoms are,
Figure FDA0003842327060000022
is the electricity price in a unit time period, P wp (t) is the power consumption of the circulation pump at time t,
Figure FDA0003842327060000023
the electric power consumed by the ith heat pump at the moment t, N is the number of heat pumps in the heat exchange station, G gb (t) is the gas consumption of the gas boiler at time t, eta gb In order to achieve the conversion efficiency of the gas boiler,
Figure FDA0003842327060000024
is the price of natural gas, k, per unit time period gb For the operating maintenance cost coefficient, k, of the gas boiler hp For the coefficient of operating maintenance of the heat pump, k wp For operating and maintenance cost factor, Q, of the circulating pump gb (t) is the output of the gas boiler at time t, P wp (t) power consumption of the circulation pump at time t, p car Is the cost per carbon emission, delta ec 、δ gb Respectively as electricity purchasing and gas burning potThe carbon emission coefficient of furnace heat supply, k is the auxiliary peak regulation time period specified by the power grid,
Figure FDA0003842327060000025
subsidizing the electricity price for assisting peak shaving in a unit time period, f hp As a function of the total electricity consumption and government subsidy policies of the heat pump system of the heat exchange station, c clean Is a subsidy price.
5. The method of claim 2, wherein the constraints comprise at least one of: user heat load demand constraint, heat pump operation constraint, heat exchange station heat pump cluster total power consumption constraint, gas boiler output constraint, primary regulating valve relative opening constraint, heat medium flow constraint and circulating pump power constraint.
6. The method of claim 5, wherein the mathematical model of the user thermal load demand constraint is as follows:
Q hp (t)+Q hes (t)=Q user (t)
θ in (t+1)=θ in (t)+(Q user (t)-Q loss (t))△t/C M M k
Q loss (t)=ξ lossin (t)-θ out (t))
Figure FDA0003842327060000027
in (t+1)-θ in (t)|≤θ ch
Q hes (t)=k hes Q gb (t)
Figure FDA0003842327060000026
in the above formula, Q hp (t) the heat-supplementing power of the heat pump system at time t, Q hes When (t) is tThe heat exchange station exchanges heat through the heat exchanger and provides heat power Q for users user (t) is the thermal load demand of the user at time t,. DELTA.t is the time interval,. Theta. in (t) indoor temperature of building at time t, Q loss (t) Heat losses associated with the physical structure of the building and the outdoor environment at time t, C M Specific heat capacity of indoor air, M k Is indoor air quality xi loss Is a coefficient related to the physical structure of the building and the outdoor environment, θ out (t) is the outdoor temperature of the building at time t,
Figure FDA0003842327060000031
the lower limit of the room temperature is,
Figure FDA0003842327060000032
at the upper limit of room temperature, [ theta ] ch Is the maximum fluctuation range of the room temperature within a unit time period, k hes For the heat exchange efficiency of heat exchangers, Q gb (t) is the heating power of the gas boiler at the time t,
Figure FDA0003842327060000033
and N is the number of heat pumps in the heat exchange station, and COP is the coefficient of performance of the heat pumps.
7. The method of claim 6, wherein the mathematical model of the heat pump operating constraints is as follows:
Figure FDA0003842327060000034
in the above-mentioned formula, the compound has the following structure,
Figure FDA0003842327060000035
is the rated power of the heat pump operation.
8. The method of claim 7, wherein the mathematical model of the heat exchange station heat pump cluster total power consumption constraint is as follows:
Figure FDA0003842327060000036
the mathematical model of the output restriction of the gas boiler is as follows:
Q gb.min ≤Q gb (t)≤Q gb.max
in the above formula, the first and second carbon atoms are,
Figure FDA0003842327060000037
and
Figure FDA0003842327060000038
minimum and maximum allowable power consumption power, Q, of a node j connected with a heat pump cluster of the heat exchange station at the moment t gb.min And Q gb.max The upper and lower output limits of the gas boiler.
9. The method of claim 7, wherein the mathematical model of the primary governor valve relative opening constraint is as follows:
Figure FDA0003842327060000039
0.1≤K V ≤0.9
Figure FDA00038423270600000310
Figure FDA00038423270600000311
in the above formula, K V Relative opening of primary valve, /) V For the stroke of the valve core at a certain opening of the regulating valve,/ V.max Is the valve core stroke when the regulating valve is fully opened, V is the heat medium flow in the pipeline, R V For regulating the adjustable ratio of valves, V max Maximum flow rate of heating medium, V, controlled by regulating valve min The minimum flow of the heating medium which can be controlled by the regulating valve.
10. The method of claim 9, wherein the mathematical model of the heat medium flow constraint is as follows:
Q gb (t)=ρ×c p ×V(t)×(T s -T h )
V min ≤V(t)≤V max
in the above formula, c p Is the specific heat capacity of the primary heat supply network heat medium, rho is the primary heat supply network heat medium density, V (T) is the heat medium flow in the pipeline at the time T, T s For temperature of the water supply, T h The temperature of the return water is shown.
11. The method of claim 10, wherein the mathematical model of the circulation pump power constraint is as follows:
Figure FDA0003842327060000041
0≤P wp ≤P wp.max
in the above formula, V e To design the operating flow, n e In order to design the running rotating speed of the water pump, n is the actual running rotating speed of the water pump,
Figure FDA0003842327060000042
for designing the power consumption of the water pump, P wp For actual operation of the pump, power consumption, P wp.max The maximum power consumption of the circulating pump.
12. An optimal control device for a central heating system, the device comprising:
the first analysis module is used for solving a pre-constructed comprehensive benefit optimization model by adopting a simplex method to obtain an optimization result;
the second analysis module is used for obtaining an optimized control scheme of the electric heat compensation equipment in the central heating system based on the optimization result;
wherein the electric supplementary heating device comprises at least one of the following: the system comprises a heat pump, a gas boiler, a primary regulating valve and a circulating pump, wherein the optimization result comprises at least one of the following components: the output of the heat pump, the output of the gas boiler, the opening of the primary regulating valve and the rotating speed of the circulating pump.
13. The apparatus of claim 12, wherein the pre-constructed synthetic benefit optimization model comprises: an objective function that aims at economic cost minimization, and constraints configured for optimal control of the central heating system.
14. The apparatus of claim 13, wherein the mathematical model of the objective function is as follows:
minF=C grid +C gas +C maintain +C car -E peak -E allowance
in the above formula, F is the comprehensive benefit of the heating system of a single heat exchange station containing heat pump heat compensation, C grid To purchase electricity cost, C gas For cost of gas, C maintain For operating maintenance costs, C car For carbon emission costs, E peak Earnings to participate in the peak shaving aid service market, E allowance It is a supplementary heating paste for clean energy.
15. The apparatus of claim 14, wherein the mathematical model of the electricity purchase cost is as follows:
Figure FDA0003842327060000043
the mathematical model of the gas cost is as follows:
Figure FDA0003842327060000051
the mathematical model of the operation and maintenance cost is as follows:
Figure FDA0003842327060000052
the mathematical model of the carbon emission cost is as follows:
Figure FDA0003842327060000053
the mathematical model of the income gained by participating in the peak shaving auxiliary service market is as follows:
Figure FDA0003842327060000054
the mathematical model of the clean energy heating subsidy is as follows:
Figure FDA0003842327060000055
in the above-mentioned formula, the compound has the following structure,
Figure FDA0003842327060000056
is the electricity price in a unit time period, P wp (t) is the power consumption of the circulation pump at time t,
Figure FDA0003842327060000057
the electric power consumed by the ith heat pump at the moment t, N is the number of heat pumps in the heat exchange station, G gb (t) is the gas consumption of the gas boiler at time t, eta gb In order to achieve the conversion efficiency of the gas boiler,
Figure FDA0003842327060000058
is the price of natural gas, k, per unit time period gb Operating and maintenance cost system for gas boilerNumber, k hp For the operating and maintenance cost coefficient, k, of the heat pump wp For the operating and maintenance cost factor, Q, of the circulating pump gb (t) is the output of the gas boiler at time t, P wp (t) power consumption of the circulation pump at time t, p car Is the cost per carbon emission, delta ec 、δ gb Carbon emission coefficients for electricity purchase and heat supply of a gas boiler are respectively, k is an auxiliary peak regulation period specified by a power grid,
Figure FDA0003842327060000059
subsidizing the electricity price for assisting peak shaving in a unit time period, f hp Function of heat pump system of heat exchange station with respect to total electricity consumption and government subsidy policy, c clean The price is subsidized.
16. The apparatus of claim 13, wherein the constraint comprises at least one of: the method comprises the following steps of user heat load demand constraint, heat pump operation constraint, heat exchange station heat pump cluster total power consumption constraint, gas boiler output constraint, primary regulating valve relative opening constraint, heat medium flow constraint and circulating pump power constraint.
17. The apparatus of claim 16, wherein the mathematical model of the user thermal load demand constraint is as follows:
Q hp (t)+Q hes (t)=Q user (t)
θ in (t+1)=θ in (t)+(Q user (t)-Q loss (t))△t/C M M k
Q loss (t)=ξ lossin (t)-θ out (t))
Figure FDA00038423270600000610
in (t+1)-θ in (t)|≤θ ch
Q hes (t)=k hes Q gb (t)
Figure FDA0003842327060000061
in the above formula, Q hp (t) the heat-supplementing power of the heat pump system at time t, Q hes (t) heat exchange station exchanges heat through heat exchanger at t moment, and heat power Q provided for user user (t) is the thermal load demand of the user at time t,. DELTA.t is the time interval,. Theta. in (t) indoor temperature of the building at time t, Q loss (t) Heat losses associated with the physical structure of the building and the outdoor environment at time t, C M Specific heat capacity of indoor air, M k Is indoor air quality xi loss Is a coefficient relating to the physical structure of the building and the outdoor environment, theta out (t) is the outdoor temperature of the building at time t,
Figure FDA0003842327060000062
the lower limit of the room temperature is,
Figure FDA0003842327060000063
at the upper limit of room temperature, [ theta ] ch Is the maximum fluctuation range of the room temperature within a unit time period, k hes For the heat exchange efficiency of heat exchangers, Q gb (t) is the heating power of the gas boiler at time t,
Figure FDA0003842327060000064
and N is the number of heat pumps in the heat exchange station, and COP is the coefficient of performance of the heat pumps.
18. The apparatus of claim 17, wherein the mathematical model of the heat pump operating constraints is as follows:
Figure FDA0003842327060000065
in the above formula, the first and second carbon atoms are,
Figure FDA0003842327060000066
is the rated power for the heat pump to operate.
19. The apparatus of claim 18, wherein a mathematical model of the heat exchange station heat pump cluster total power consumption constraint is as follows:
Figure FDA0003842327060000067
the mathematical model of the output restriction of the gas boiler is as follows:
Q gb.min ≤Q gb (t)≤Q gb.max
in the above formula, the first and second carbon atoms are,
Figure FDA0003842327060000068
and
Figure FDA0003842327060000069
minimum and maximum allowable power consumption power, Q, of a node j connected with a heat pump cluster of the heat exchange station at the moment t gb.min And Q gb.max The upper and lower output limits of the gas boiler.
20. The apparatus of claim 18, wherein the mathematical model of the primary regulation valve relative opening constraint is as follows:
Figure FDA0003842327060000071
0.1≤K V ≤0.9
Figure FDA0003842327060000072
Figure FDA0003842327060000073
in the above formula, K V Relative opening of primary valve, /) V For the stroke of the valve core at a certain opening of the regulating valve,/ V.max The valve core stroke when the regulating valve is fully opened, V is the heat medium flow in the pipeline, R V For adjustable ratio of regulating valves, V max Maximum flow rate of heating medium, V, controlled by regulating valve min The minimum flow of the heating medium which can be controlled by the regulating valve.
21. The apparatus of claim 20, wherein the mathematical model of the heating medium flow restriction is as follows:
Q gb (t)=ρ×c p ×V(t)×(T s -T h )
V min ≤V(t)≤V max
in the above formula, c p Is the specific heat capacity of the primary heat supply network heat medium, rho is the primary heat supply network heat medium density, V (T) is the heat medium flow in the pipeline at the time T, T s For supply water temperature, T h The temperature of the return water.
22. The apparatus of claim 21, wherein the mathematical model of the circulation pump power constraint is as follows:
Figure FDA0003842327060000074
0≤P wp ≤P wp.max
in the above formula, V e To design the operating flow, n e In order to design the running rotating speed of the water pump, n is the actual running rotating speed of the water pump,
Figure FDA0003842327060000075
for designing the power consumption of the water pump, P wp For actual operation of the pump consuming power, P wp.max The maximum power consumption of the circulating pump.
23. A computer device, comprising: one or more processors;
the processor to store one or more programs;
the one or more programs, when executed by the one or more processors, implement a method of optimizing control of a central heating system according to any one of claims 1 to 11.
24. A computer-readable storage medium, characterized in that it has stored thereon a computer program which, when executed, carries out a method for the optimized control of a central heating system according to any one of claims 1 to 11.
CN202211108359.1A 2022-09-13 2022-09-13 Optimization control method and device for central heating system Pending CN115682090A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211108359.1A CN115682090A (en) 2022-09-13 2022-09-13 Optimization control method and device for central heating system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211108359.1A CN115682090A (en) 2022-09-13 2022-09-13 Optimization control method and device for central heating system

Publications (1)

Publication Number Publication Date
CN115682090A true CN115682090A (en) 2023-02-03

Family

ID=85063146

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211108359.1A Pending CN115682090A (en) 2022-09-13 2022-09-13 Optimization control method and device for central heating system

Country Status (1)

Country Link
CN (1) CN115682090A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116717839A (en) * 2023-08-10 2023-09-08 陕西拓普索尔电子科技有限责任公司 Heating control method, control device and heating system

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116717839A (en) * 2023-08-10 2023-09-08 陕西拓普索尔电子科技有限责任公司 Heating control method, control device and heating system

Similar Documents

Publication Publication Date Title
CN112583021B (en) Comprehensive energy system optimal scheduling method and device considering comprehensive demand response
CN105240897B (en) A kind of accumulation of heat peak adjusting device for heating system
CN109711080B (en) Multi-time scale optimization operation method for combined cooling heating and power system
CN107276096B (en) A kind of distribution network voltage control method of photovoltaic and air conditioner load coordination optimization
CN109269021B (en) Energy-saving operation optimization scheduling method for air conditioning system
CN106593771B (en) Electric heating decoupled system and its control method in a kind of co-generation unit
CN112398164A (en) Micro-energy-source network group optimization operation and cost distribution method containing shared energy storage system
EA025665B1 (en) Method of regulating a plant comprising cogenerating installations and thermodynamic systems intended for air conditioning and/or heating
CN112434861B (en) Electric heating coordination optimization scheduling method and device considering feedback adjustment of heating power pipe network
CN110829408A (en) Multi-domain scheduling method considering energy storage power system based on power generation cost constraint
CN107860057A (en) A kind of cogeneration of heat and power heating system thermic load economic optimization dispatching method
CN113446661A (en) Intelligent and efficient heat supply network operation adjusting method
CN115682090A (en) Optimization control method and device for central heating system
CN116344863B (en) Combined heat and power thermal management system of multi-fuel cell system and control method thereof
CN111625961A (en) Comprehensive energy system collaborative optimization operation regulation and control method
CN116451945A (en) Optimal scheduling method for virtual electric heating comprehensive energy power plant by considering pipeline hydraulic dynamic
CN113610305B (en) Optimized scheduling method of comprehensive energy system
CN113537618B (en) Comprehensive energy system optimization scheduling method considering resident user demand response
CN109405055A (en) A kind of heat supply and accumulation of heat decouple operating system with heat source
CN113065759A (en) Comprehensive energy system optimal scheduling method and system considering energy quality
Xia et al. Economic Optimal Scheduling of Integrated Energy System Considering Integrated Demand Response
CN115795842B (en) Calculation method and system for building cluster heat balance
CN118112938B (en) Short-period control strategy optimization method for combined type ground source heat pump system
US20240240801A1 (en) Control Units and Methods for Controlling Operation of a Heat Generation Plant of an Energy System
CN118134059B (en) Method, system, equipment and storage medium for collaborative optimization of electric heating network and building user

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication