CN113191638A - Electric heating coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network - Google Patents

Electric heating coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network Download PDF

Info

Publication number
CN113191638A
CN113191638A CN202110481506.9A CN202110481506A CN113191638A CN 113191638 A CN113191638 A CN 113191638A CN 202110481506 A CN202110481506 A CN 202110481506A CN 113191638 A CN113191638 A CN 113191638A
Authority
CN
China
Prior art keywords
heat
network
scheduling
temperature
user
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110481506.9A
Other languages
Chinese (zh)
Other versions
CN113191638B (en
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.)
Tsinghua University
Metering Center of State Grid Jibei Electric Power Co Ltd
Original Assignee
Tsinghua University
Metering Center of State Grid Jibei 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 Tsinghua University, Metering Center of State Grid Jibei Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN202110481506.9A priority Critical patent/CN113191638B/en
Priority claimed from CN202110481506.9A external-priority patent/CN113191638B/en
Publication of CN113191638A publication Critical patent/CN113191638A/en
Application granted granted Critical
Publication of CN113191638B publication Critical patent/CN113191638B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Educational Administration (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides an electric heating coordination scheduling method and device considering heat transmission dynamic characteristics of a heat supply network, wherein the method comprises the following steps: establishing an electric heating coordination scheduling model; establishing constraint conditions based on the heat transmission dynamic characteristics of a heat supply network and a multi-stage heat exchange process; obtaining the current value of a decision variable, and calculating through an electric heating coordinated scheduling model according to the current value of the decision variable and a constraint condition to obtain the current value of a target variable; determining a target interval of a target variable based on the set scheduling condition; feeding back and adjusting decision variables according to the relative relation between the current value of the target variable and the target interval, thereby realizing electric heating coordination scheduling; the influence of the transmission characteristics of the heat supply network on the target variable is definitely considered; according to the invention, heat transmission constraint is established based on a multistage heat exchange process, consideration factors of an electric heat coordination scheduling model are refined, and the model optimization result is more reliable and the model compatibility is better.

Description

Electric heating coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network
Technical Field
The invention relates to the technical field of electric heating coordination, in particular to an electric heating coordination scheduling method and device considering heat transmission dynamic characteristics of a heat supply network.
Background
The centralized heating network in the electric heating comprehensive energy system has passive heat storage capacity due to thermal inertia and can provide regulation capacity for the power system. In order to quantitatively evaluate the passive heat storage performance of the heat supply network and improve the consumption effect of renewable energy, an electric heating coordination optimization scheduling model containing the heat transmission dynamic characteristic of the heat supply network needs to be established. Taking the wind power consumption problem in the northeast China as an example, a transmission dynamic characteristic model of a pipe network is brought into an electric heating coordination optimization model, so that a relatively accurate heat supply network transmission delay rule can be obtained, a heat source is quantitatively and optimally controlled, the heat load of a user at the moment is provided in advance before delay, and the heat load of the user at the moment with heavy wind power generation and high heat load at night is also provided for a period of time in advance. Compared with the situation without considering the transmission delay, the method has the advantages that the heat supply amount and the output of the thermoelectric generator set are lower in the night time with large wind power generation and high heat load, and larger space is provided for wind power grid-connected consumption. For example, the heat load is highest and the wind power output is maximum at 3 am, if the transmission delay of a pipe network is not considered, the heat supply amount (namely, the heat output of the cogeneration unit) is highest at 3 am, the power output of the thermoelectric unit is maximum, and the wind power grid-connected space is possibly insufficient, so that the maximum wind abandon is achieved; if the transmission delay is considered, the maximum heat supply amount needs to be arranged in advance (for example, 2 am) by the heat source, so that the peak time of 3 am with the maximum wind power can be staggered, and the wind abandon amount of 3 am is reduced.
According to the method, the temperature dynamic distribution rule of each node of the pipe network can be obtained only by solving the heat transmission dynamic characteristic of the heat supply network, so that the real-time quantitative relation between the heat supply temperature of the heat source and the hot water temperature reaching a user is determined, the time of the heat supply network supplying heat in advance is determined, and accurate constraint is provided for electric heat coordination optimization.
Three problems exist in the existing electric heating coordination optimization considering the passive heat storage of a pipe network. First, a dynamic model of a heat supply network with both accuracy and compatibility with the power system is lacking. The traditional thermodynamic model cannot be directly substituted into electric heating coordination optimization, and the compatibility is poor. Second, many documents do not adequately consider the operational characteristics of components that are closely related to the operational characteristics of the heat supply network, based on the inaccuracy of the heat supply network model. The Jiangxianbao of Zhejiang university considers the time delay characteristic of the heat supply network and the variable flow variable temperature regulation into electric heat coordination optimization, but the adopted heat supply network model simplifies the temperature into the average heat supply temperature at the inlet and the outlet, does not consider the logarithmic average temperature difference in the heat exchange constraint, and cannot reflect the heat transmission and conversion process of each node in the thermodynamic system. Thirdly, heat exchange constraints among the primary network, the secondary network and the heat consumers are not considered, real-time change of indoor temperature is not obtained, or only heat exchange balance is considered and a logarithmic mean temperature difference equation is ignored, so that the obtained theoretical heat supply amount is inconsistent with the reality, and accurate constraints can not be provided for making heat supply and power supply plans.
In the above-mentioned problem of the electric heat coordination optimization model considering the dynamic characteristics of the heat supply network, the heat supply network model determines the real-time mapping relationship between the heat supply temperature of the heat source and the temperature of the heat medium reaching the user. The change rule of the primary network and the secondary network in the thermodynamic system, the change rule of the secondary network and the indoor temperature of the user also greatly depends on the precision of the heat supply network dynamic model. Therefore, analysis of the heat transmission dynamic characteristics of the pipe network is the basis for realizing dynamic adjustment and electric heat coordination optimization of the heating system. Only a heat transmission dynamic model which is good in precision and compatible with an electric power system algorithm is established, a pipe network thermal dynamic topological model can be established, and the heat exchanger and the operation characteristics of heat users are further combined, so that the dynamic response relation between a heat source and each user in a complex pipe network is accurately mastered, the purpose of quantitatively adjusting heat supply and generator set operation parameters according to a renewable energy consumption target is realized, and a theoretical basis is provided for exploring and utilizing the natural adjustment potential of the pipe network in an electric heating comprehensive energy system.
In summary, the method and the device for electric heat coordination scheduling considering the heat transmission dynamic characteristics of the heat supply network are provided, based on the calculation of the multi-stage heat exchange process, the problems that in the prior art, the compatibility/precision is insufficient, the heat transmission and conversion process of each node in the thermodynamic system cannot be reflected, and the heat exchange constraint between the primary network/the secondary network/the heat users is not considered can be solved, and the method and the device have high practical value and significance.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides an electric heating coordination scheduling method and device considering the heat transmission dynamic characteristics of a heat supply network.
The invention provides an electric heating coordination scheduling method considering heat transmission dynamic characteristics of a heat supply network, which comprises the following steps:
establishing an electric heating coordination scheduling model;
establishing heat supply network constraint based on the heat supply network heat transmission dynamic characteristic and the multistage heat exchange process as constraint conditions;
obtaining the current value of a decision variable, and calculating through an electric heating coordinated scheduling model according to the current value of the decision variable and a constraint condition to obtain the current value of a target variable;
determining a target interval of a target variable or a target value of the target variable based on the set scheduling condition;
feeding back and adjusting the decision variable according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, thereby realizing electric heating coordinated scheduling;
the electric heating coordination scheduling model is a model for adjusting a target variable through a decision variable under a set constraint condition;
the decision variables comprise any one or any combination of more of hourly heat source heat supply temperature, heat supply pipe network initial temperature, heat supply quantity of the thermoelectric unit, electric output, thermal power unit output and hourly wind power output;
the target variable comprises any one or any combination of more of air abandonment amount, heat deviation amount, heat supply amount, heat loss, user heat gain amount, user indoor temperature and wind power consumption amount;
the multistage heat exchange process comprises a heat transmission process and/or a heat exchange process of a heat supply pipe network; the heat transmission process of the heat supply pipe network comprises a heat transmission process of a heat medium supply pipeline and/or a heat transmission process of a heat medium recovery pipeline; the heat exchange process comprises any one or any combination of a heat exchange process between a primary network and a secondary network, a heat exchange process between the secondary network and a user, and a heat exchange process between the user and an outdoor environment.
According to the electric heat coordination scheduling method considering the heat transmission dynamic characteristics of the heat supply network, the heat supply network constraint comprises a first function corresponding to a heat exchange process:
Figure BDA0003049429870000041
Figure BDA0003049429870000042
Figure BDA0003049429870000043
in the formula:
Tc,in,tthe inlet temperature of the heat medium of the secondary network at the time t; t isoutdoor,tIs the outdoor atmospheric temperature at time t; t ish,in,tThe inlet temperature of the heat medium of the primary network at the time t;
Th,out,tthe outlet temperature of the primary network heat medium at the moment t; gcMass flow of the secondary network heat medium; ghMass flow rate of the primary net heat medium;
Tuser,tthe indoor temperature of the user at the moment t;
sh,c、sevl、sc,user、s12hare respectively:
Figure BDA0003049429870000051
Figure BDA0003049429870000052
Figure BDA0003049429870000053
wherein e is the base of the natural logarithm; k is a radical ofh,cThe heat exchange coefficient of the primary network heat exchanger and the secondary network heat exchanger is obtained; a. theh,cThe heat exchange area of the primary net and the secondary net heat exchanger is shown; c is the specific heat capacity of the heat medium; k is a radical ofevlAverage heat exchange coefficient of the building envelope for the user; a. theevlAverage heat exchange area of the building envelope for the user; k is a radical ofc,userThe heat exchange coefficient of the radiator between the secondary network and the user is shown; a. thec,userThe heat exchange area of the radiator between the secondary network and the user.
According to the electric heating coordinated scheduling method considering the heat transmission dynamic characteristics of the heat supply network, the heat supply network constraint comprises a second function corresponding to the heat transmission process of the heat supply network:
Figure BDA0003049429870000054
Figure BDA0003049429870000055
in the formula, phim,nThe outlet excess temperature is the outlet excess temperature with a space node of m and a time node of n; i. j is a space node serial number and a time node serial number respectively; the surplus temperature refers to the difference between the temperature of the heat medium and the temperature of the soil;
parameter An、BnAnd A'nSatisfies the following conditions:
Figure BDA0003049429870000061
Figure BDA0003049429870000062
Figure BDA0003049429870000063
wherein, C is a permutation and combination operator;
the parameters a and r satisfy:
Figure BDA0003049429870000064
Figure BDA0003049429870000065
in the formula, rho and c are respectively the density of the heat medium and the specific heat capacity of the heat medium; A. d is the cross-sectional area and the diameter of the pipeline respectively; v is the fluid flow rate; k is the heat leakage coefficient between the fluid in the pipeline and the soil; x and tau are respectively a space variable and a time variable.
According to the electric heat coordination scheduling method considering the heat transmission dynamic characteristic of the heat supply network, the constraint condition further comprises the constraint of a cogeneration unit;
the constraint of the cogeneration unit comprises a thermoelectric ratio constraint, a climbing constraint, a unit limit output constraint, a heat source water supply temperature constraint and a scheduling period constraint;
the thermoelectric ratio constraints are:
Figure BDA0003049429870000071
in the formula, epsilon is a set thermoelectric ratio of the cogeneration unit; qCHP,tThe thermal output of the cogeneration unit at the moment t; pCHP,tThe electric output of the cogeneration unit at the moment t;
the climbing restriction is as follows:
Figure BDA0003049429870000072
in the formula, PCHP,up、PCHP,downRespectively restricting the upward and downward climbing of the set cogeneration unit; pCHP,t-1The power output of the cogeneration unit at the time of t-1; pthp,up、 Pthp,downRespectively performing upward and downward climbing constraints on the set thermal power cogeneration unit; pthp,t、 Pthp,t-1The thermal power cogeneration unit power output at the time t and the time t-1 respectively; wherein the t-1 moment is the last moment of the t moment;
the unit ultimate output restriction package is as follows:
Figure BDA0003049429870000073
in the formula, PCHP,max、PCHP,minThe maximum power output and the minimum power output of the cogeneration unit are respectively; pthp,max、Pthp,minRespectively the maximum and minimum power output of the thermal power cogeneration unit;
the heat source water supply temperature constraint is as follows:
40℃≤Tsource,t≤90℃
in the formula, Tsource,tThe temperature of the heat source water supply at the time t;
the scheduling period constraint means that the thermal scheduling period is adjusted based on a set power scheduling period value, so that the electric heating coordination scheduling period, the thermal scheduling period and the power scheduling period are consistent and are set power scheduling period values.
According to the electric heating coordination scheduling method considering the heat transmission dynamic characteristic of the heat supply network, the constraint condition further comprises the constraint of the power grid and the wind generating set;
the power grid and wind power generator set constraints comprise power grid power balance constraints, wind power limit output constraints and line capacity constraints;
the power balance constraint of the power grid is as follows:
PCHP,t+Pwind,t=Puser,t+Ppump,t
in the formula, PCHP,tThe electric output of the cogeneration unit at the moment t; pwind,tThe actual wind power output at the moment t is obtained; puser,tThe load is the residential electricity load at the moment t; ppump,tThe power consumption of the circulating pump at the time t;
the wind power limit output constraint is as follows:
Pwind,t≤Pwind,max,t
in the formula, Pwind,max,tThe wind power limit output at the time t is obtained;
the line capacity constraints are:
Ptr,min≤Ptr,t≤Ptr,max
in the formula, Ptr,tThe transmission capacity of the power grid at the moment t; ptr,min、Ptr,maxThe minimum and maximum transmission capacities of the power grid at the moment t are respectively.
According to the electric heating coordination scheduling method considering the heat transmission dynamic characteristic of the heat supply network, the scheduling condition is that the air abandoning amount at all times in the whole day is minimum; the target variable corresponding to the scheduling condition with the minimum air abandon amount at all times in the whole day is the sum of the air abandon amounts, and the corresponding target value is 0.
According to the electric heat coordination scheduling method considering the heat transmission dynamic characteristics of the heat supply network, the scheduling condition is that the absolute value of the total heat deviation of the whole day is minimum, the target variable corresponding to the scheduling condition with the minimum absolute value of the total heat deviation of the whole day is the sum of the absolute values of the heat deviations, and the corresponding target value is 0.
The invention also provides an electric heating coordination scheduling device considering the heat transmission dynamic characteristics of the heat supply network, which comprises a model establishing module, a constraint condition setting module, a calculating module, a scheduling module and a feedback adjusting module;
the model establishing module can establish an electric heating coordination scheduling model;
the constraint condition setting module can establish heat transmission constraint based on a multi-stage heat exchange process to serve as a constraint condition;
the calculation module can obtain the current value of the decision variable, and according to the current value of the decision variable and the constraint condition, the current value of the target variable is calculated and obtained through an electric heating coordination scheduling model;
the scheduling module can determine a target interval of a target variable or a target value of the target variable based on a set scheduling condition;
the feedback regulation can feed back and regulate decision variables according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, so that the electric heating coordinated dispatching is realized;
the electric heating coordination scheduling model is a model for adjusting a target variable through a decision variable under a set constraint condition;
the decision variables comprise any one or any combination of more of hourly heat source heat supply temperature, heat supply pipe network initial temperature, heat supply quantity of the thermoelectric unit, electric output, thermal power unit output and hourly wind power output;
the target variable comprises any one or any combination of more of air abandonment amount, heat deviation amount, heat supply amount, heat loss, user heat gain amount, user indoor temperature and wind power consumption amount;
the multistage heat exchange process comprises a heat transmission process and/or a heat exchange process of a heat supply pipe network; the heat transmission process of the heat supply pipe network comprises a heat transmission process of a heat medium supply pipeline and/or a heat transmission process of a heat medium recovery pipeline; the heat exchange process comprises any one or any combination of a heat exchange process between a primary network and a secondary network, a heat exchange process between the secondary network and a user, and a heat exchange process between the user and an outdoor environment.
The invention further provides an electronic device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein when the processor executes the program, the steps of the electric heating coordination scheduling method considering the dynamic characteristics of heat transmission of the heat supply network are realized.
The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when being executed by a processor, implements the steps of the method for coordinated electric and heat scheduling taking into account the thermal transfer dynamics of a heat network as described in any of the above.
According to the electric heating coordination scheduling method and device considering the heat transmission dynamic characteristics of the heat supply network, the heat transmission constraint is established based on the multi-stage heat exchange process, the heat transmission and conversion process of each node in a thermodynamic system is reflected, and the heat exchange constraint between a primary network/a secondary network/a heat user is established, so that the consideration factors of an electric heating coordination scheduling model are refined, the model optimization result is more reliable, and the model compatibility is better.
Drawings
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a schematic flow chart of an electric heat coordination scheduling method taking account of heat transmission dynamic characteristics of a heat supply network provided by the present invention;
FIG. 2 is a schematic diagram of a heat exchange process of multiple loops of a heat supply pipe network according to an embodiment of the present invention;
FIG. 3 is a schematic view of a buried pipe provided by an embodiment of the present invention;
FIG. 4 is a schematic diagram of a preliminary thermoelectric comparison of a heat supply network according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of an equivalent thermal circuit of a heat supply network provided by an embodiment of the present invention;
FIG. 6 is a schematic diagram of an electric heating integrated energy system for a certain area according to an embodiment of the present invention;
fig. 7 is a schematic diagram of a time-by-time electrical load and wind power limit output curve of a single user node case provided in the embodiment of the present invention;
FIG. 8 is a schematic diagram of a time-by-time outdoor temperature curve for a single user node case according to an embodiment of the present invention
FIG. 9 is a schematic view of a time-lapse heating temperature curve of various comparative examples provided in the examples of the present invention;
FIG. 10 is a schematic graph of a time-lapse indoor temperature profile for various comparative examples provided by examples of the present invention;
FIG. 11 is a schematic view of a time-wise wind abandon curve of different comparative examples provided by the embodiment of the present invention;
fig. 12 is a schematic structural diagram of an electronic device provided in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, 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 method and apparatus for coordinating and scheduling electric heating taking into account the thermal transmission dynamic characteristics of the heat supply network according to the present invention will be described with reference to fig. 1-11.
The definition of 'taking account of the multi-stage heat exchange process' in the subject matters of the invention and the embodiments refers to that in the modeling process of the electric heating coordination scheduling model and the electric heating coordination scheduling process, the multi-stage heat exchange process of the heat supply pipe network part is considered, namely an equation is established on the basis of the multi-stage heat exchange process to restrict the electric heating coordination scheduling.
As shown in fig. 1, an embodiment of the present invention provides an electric heating coordination scheduling method considering heat transfer dynamic characteristics of a heat supply network, including:
step 1, establishing an electric heating coordination scheduling model;
step 2, establishing heat supply network constraint based on the heat supply network heat transmission dynamic characteristic and the multi-stage heat exchange process as constraint conditions;
step 3, obtaining the current value of the decision variable, and calculating to obtain the current value of the target variable through an electric heating coordination scheduling model according to the current value of the decision variable and the constraint condition;
step 4, determining a target interval of a target variable or a target value of the target variable based on the set scheduling condition;
step 5, feeding back and adjusting decision variables according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, thereby realizing electric heating coordination scheduling;
the electric heating coordination scheduling model is a model for adjusting a target variable through a decision variable under a set constraint condition;
the decision variables comprise any one or any combination of more of hourly heat source heat supply temperature, heat supply pipe network initial temperature, heat supply quantity of the thermoelectric unit, electric output, thermal power unit output and hourly wind power output;
the target variable comprises any one or any combination of more of air abandonment amount, heat deviation amount, heat supply amount, heat loss, user heat gain amount, user indoor temperature and wind power consumption amount;
the multistage heat exchange process comprises a heat transmission process and/or a heat exchange process of a heat supply pipe network; the heat transmission process of the heat supply pipe network comprises a heat transmission process of a heat medium supply pipeline and/or a heat transmission process of a heat medium recovery pipeline; the heat exchange process comprises any one or any combination of a heat exchange process between a primary network and a secondary network, a heat exchange process between the secondary network and a user, and a heat exchange process between the user and an outdoor environment.
The beneficial effect of this embodiment lies in:
the heat transmission constraint is established based on the multi-stage heat exchange process, the heat transmission and conversion process of each node in the thermodynamic system is reflected, and the heat exchange constraint between the primary network/the secondary network/the heat users is established, so that the consideration factors of the electric heat coordination scheduling model are refined, the model optimization result is more reliable, and the model compatibility is better.
According to the above embodiment, in the present embodiment:
the heat grid constraints include a first function corresponding to a heat exchange process:
Figure BDA0003049429870000121
Figure BDA0003049429870000122
Figure BDA0003049429870000131
in the formula:
Tc,in,tthe inlet temperature of the heat medium of the secondary network at the time t; t isoutdoor,tIs the outdoor atmospheric temperature at time t; t ish,in,tThe inlet temperature of the heat medium of the primary network at the time t;
Th,out,tthe outlet temperature of the primary network heat medium at the moment t; gcMass flow of the secondary network heat medium; ghMass flow rate of the primary net heat medium;
Tuser,tthe indoor temperature of the user at the moment t;
sh,c、sevl、sc,user、s12hare respectively:
Figure BDA0003049429870000132
Figure BDA0003049429870000133
Figure BDA0003049429870000134
wherein e is the base of the natural logarithm; k is a radical ofh,cThe heat exchange coefficient of the primary network heat exchanger and the secondary network heat exchanger is obtained; a. theh,cThe heat exchange area of the primary net and the secondary net heat exchanger is shown; c is the specific heat capacity of the heat medium; k is a radical ofevlAverage heat exchange coefficient of the building envelope for the user; a. theevlAverage heat exchange area of the building envelope for the user; k is a radical ofc,userThe heat exchange coefficient of the radiator between the secondary network and the user is shown; a. thec,userThe heat exchange area of the radiator between the secondary network and the user.
The derivation of the first function is as follows.
As shown in fig. 2, the heat transfer from the primary network heat medium to the user involves three layers of heat exchange processes, namely, a primary network and a secondary network, a secondary network heat medium and a heat user, and a heat user and the external environment; the heat exchange equations of the three heat exchange loops are listed in sequence according to the schematic diagram of the heat exchange process of the multiple loops of the heat supply network shown in fig. 2.
(1) Heat exchange process between primary net and secondary net
The formula (1a), the formula (1b), the formula (1c) and the formula (1d) are respectively a heat exchange quantity equation, a logarithmic mean temperature difference equation, a heat balance equation of a primary network heat medium and a heat balance equation of a secondary network heat medium between the primary network and the secondary network. The heat exchanger is arranged in a countercurrent mode, the four modes are combined, and the expression of the outlet temperature of the primary network and the outlet temperature of the secondary network with respect to the inlet temperature of the primary network and the secondary network can be solved.
Qh,c,t=kh,cAh,cΔTm,h,c,t (1a)
Figure BDA0003049429870000141
Qh,c,t=cGh(Th,in,t-Th,out,t) (1c)
Qh,c,t=cGc(Tc,out,t-Tc,in,t) (1d)
In the formula, Qh,c,tIs the heat exchange amount between the primary and secondary nets, kh,cIs the heat exchange coefficient of the heat exchanger between the primary and secondary networks, Ah,cIs the heat exchange area, delta T, of the primary and secondary internetwork heat exchangersm,h,c,tIs the logarithmic mean temperature difference, G, of the heat exchanger between the primary and secondary networksh、GcRespectively represents the mass flow rate, T, of the primary network heat medium water and the secondary network heat medium waterh,in,t、Th,out,tInlet and outlet temperatures, T, of primary network heating medium water respectivelyc,in,t、Tc,out,tThe inlet temperature and the outlet temperature of the secondary network heating medium water are respectively.
(2) Heat exchange process between secondary network and heat consumer
And (2) analyzing the heat exchange process of the heat user as a research object, namely the heat release amount of the secondary network to the user is equal to the heat release amount of the user and the outdoor, wherein the formulas (1e), (1f), (1g) and (1h) are a heat exchange amount equation, a logarithmic mean temperature difference equation, a secondary network heat medium water-heat balance equation and a user and external heat exchange balance equation between the secondary network and the user respectively.
Qc,user,t=kc,userAc,userΔTm,c,user,t (1e)
Figure BDA0003049429870000151
Qc,user,t=cGc(Tcs,t-Tcr,t) (1g)
Qc,user,t=kevlAevl(Tuser,t-Toutdoor,t) (1h)
Wherein Q isc,user,tIndicates the heat exchange amount between the secondary network heat medium water and the user, kc,user、Ac,userRespectively represents the heat exchange coefficient and the heat exchange area, delta T, of the radiator between the heat medium water of the secondary network and the userm,c,user,tIs the logarithmic mean temperature difference, T, of the heat exchanger between the primary and secondary networkscs,t、Tcr,tThe temperature of the secondary network heat medium water supply and the return water at the time t are respectively; k is a radical ofevlRepresents the average heat exchange coefficient of the building envelope, AevlRepresenting the average heat exchange area of the building envelope; t isuser,tFor the user's indoor temperature, Toutdoor,tIs the outdoor atmospheric temperature.
(3) Primary network, secondary network and user heat exchange process simultaneous solving method
For the multi-loop heat exchange process, many intermediate parameters exist in calculation, such as the temperature of a heat medium water inlet and outlet of an intermediate loop. Taking the heat exchange process shown in fig. 2 as an example, the supply and return water temperature of the secondary network is determined by both the primary network heat exchange equation and the user heat exchange equation. Therefore, the relation between the primary network water supply temperature and the indoor temperature and other parameters is solved through the joint formula (1a) to the formula (1 h).
Supply temperature T of hot media water at user sidecs,tNamely the outlet temperature T of the secondary network heating medium water of the heating plantc,out,tReturn temperature T of heat medium water at user sidecr,tNamely the inlet temperature T of the secondary network heat medium water of the heat supply stationc,in,tAs shown in formula (1 i).
Figure BDA0003049429870000161
Therefore, the inlet and outlet temperature difference (T) of the secondary network heat medium water at the heating station sidec,out,t-Tc,in,t) Temperature difference (T) of water supply and return of secondary network at user sidecs,t-Tcr,t) Expressed as shown in formula (1 j):
Tc,out,t-Tc,in,t=Tcs,t-Tcr,t (1j)
will Tc,in,tSubstituting to obtain the indoor temperature Tuser,tAbout the temperature T of a heat medium water inlet of a primary network of a heat exchange stationh,in,tOutdoor atmospheric temperature Toutdoor,tAnd the relation with the heat exchange parameters of the heat exchanger of the heat exchange station, the heat exchange parameters of the radiator and the heat exchange parameters of the user enclosure structure is shown as a formula (1 k).
Figure BDA0003049429870000171
The dynamically changed parameter in the formula (1k) is only the temperature T of the heat exchange station primary network heat medium water inleth,in,tAnd outdoor atmospheric temperature Toutdoor,tAnd the variable intermediate variable related to the secondary network is not contained, so that the variable intermediate variable can be directly used in an electric heating coordination optimization model.
Further simplifying the equation (1k) and establishing the equation (1L).
Figure BDA0003049429870000181
The formula (1k) can be rewritten as the formula (1.1), the formula (1.2) and the formula (1.3)
Figure BDA0003049429870000182
Figure BDA0003049429870000183
Figure BDA0003049429870000191
The beneficial effect of this embodiment lies in:
the heat transmission and heat exchange process in the pipe network comprises five processes: the heat transmission process comprises heat transmission of a water supply pipeline, heat exchange of a primary network and a secondary network, heat exchange of the secondary network and a user, heat exchange of the user and an outdoor environment and heat transmission of a water return pipeline. In the embodiment, a step-by-step heat exchange equation is established for each node, and a pipe network thermodynamic dynamic model considering pipe network heat transmission is established, so that in the prior art, even if heat exchange constraint is considered in electric heat coordination optimization considering a pipe network, a secondary network is not considered generally, but only a simple heat exchange process between a primary network and the secondary network is considered, for example, heat exchange between heat medium water of the secondary network and users is ignored, and the like. In the embodiment, a three-layer heat exchange process is considered comprehensively, a heat balance and heat exchange model containing three-layer heat exchange is established aiming at the problem of multi-loop heat exchange, and an expression of indoor temperature about primary network heat supply temperature and outdoor temperature is directly obtained.
According to any of the embodiments described above, in this embodiment:
the heat supply network constraints include a second function corresponding to a heat transfer process of the heat supply network:
Figure BDA0003049429870000192
Figure BDA0003049429870000193
in the formula, phim,nThe outlet excess temperature is the outlet excess temperature with a space node of m and a time node of n; i. j is a space node serial number and a time node serial number respectively; the surplus temperature refers to the difference between the temperature of the heat medium and the temperature of the soil;
parameter An、BnAnd A'nSatisfies the following conditions:
Figure BDA0003049429870000201
Figure BDA0003049429870000202
Figure BDA0003049429870000203
wherein, C is a permutation and combination operator;
the parameters a and r satisfy:
Figure BDA0003049429870000204
Figure BDA0003049429870000205
in the formula, rho and c are respectively the density of the heat medium and the specific heat capacity of the heat medium; A. d is the cross-sectional area and the diameter of the pipeline respectively; v is the fluid flow rate; k is the heat leakage coefficient between the fluid in the pipeline and the soil; x and tau are respectively a space variable and a time variable.
Specifically, the derivation process of the second function is as follows.
Referring to the directly-buried pipe shown in fig. 3, under a set equivalent thermal circuit condition, for a fluid infinitesimal with a length dx in a heat supply pipe network pipeline in a soil environment, considering the change of water temperature along with the extension direction x and time τ of the pipeline, a one-dimensional energy conservation equation is established and is expressed as formula (2 a):
Figure BDA0003049429870000211
in the formula (2a), rho and c are respectively the density of the heat medium and the specific heat capacity of the heat medium; a. D is the cross-sectional area and the diameter of the pipeline respectively; t is the temperature of the fluid infinitesimal; v is the fluid flow rate; k is the heat leakage coefficient between the fluid in the pipeline and the soil; t issoilIs the soil temperature; x and tau are respectively a space variable and a time variable;
the equivalent thermal circuit conditions comprise constant soil temperature, same fluid temperature on the same cross section and positive direction taking the extending direction of the heat supply pipe network pipeline as an x axis;
for each term on the left side of the medium number in the formula (2a), integrating in the pipeline extending direction x and time τ, and applying a step difference format, a first-order windward difference format and an implicit difference format to the first term, the second term and the third term on the left side of the medium number in the formula (2a) respectively to obtain a formula (2 b):
Figure BDA0003049429870000212
in the formula (2b), Ti,jRepresenting the temperature of the ith pipeline space node in the numerical calculation at the jth moment;
referring to the preliminary thermoelectric analogy diagram of the heat supply network shown in fig. 4, based on equation (2b), a linear equivalent thermal circuit is established:
regarding each fluid micro element as an equivalent thermoelectric node in an equivalent thermal circuit, and marking the equivalent thermoelectric nodes which are not end points and are positioned at the ith space node and the jth time node as T'i,jThen T'i,jThe equivalent thermal circuit is connected with other parts of the equivalent thermal circuit in at least four directions, and the directions are respectively marked as a first direction, a second direction, a third direction and a fourth direction;
let equivalent thermoelectric node T'i,jTemperature of itself is Ti,j,T′i,jThe self internal energy of the corresponding fluid infinitesimal element can be equivalent to one end and T'i,jCapacitor C connected and grounded at one endi,jThen the first direction is from T'i,jWarp Ci,jThe direction of grounding;
T′i,jthe heat convection between the corresponding fluid micro-element and the immediately adjacent fluid micro-element can be equivalent to one end and T'i,jConnected, one end of which is connected with T'i-1,jConnected power supply Ei,jThen the second direction is from T'i,jWarp Ei,jGet T'i-1,jThe direction of (a);
temperature T of soilsoilCan be equivalent to a constant potential point Ts,T′i,jThe thermal resistance of the heat exchange between the corresponding fluid micro element and the soil can be equivalent to one end and T'i,jConnection, one end and constant potential point TsConnected resistor Ri,jThen the third direction is from T'i,jWarp Ri,jConnect TSThe direction of (a);
T′i,jthe thermal convection between the corresponding fluid infinitesimal and the next adjacent fluid infinitesimal can be equivalent to one end and T'i,jConnected, one end of which is connected with T'i+1jConnected power supply Ei+1,jThen the fourth direction is from T'i,jWarp Ei+1,jGet T'i+1,jThe direction of (a);
referring to fig. 5, a linear equivalent thermal circuit is simplified:
prepared from T'i,jThe self internal energy of the corresponding fluid element is recorded as a vector q in the first direction1(ii) a Prepared from T'i,jThe thermal convection between the corresponding fluid infinitesimal and the immediately preceding adjacent fluid infinitesimal is denoted as a vector q in the second direction2(ii) a Prepared from T'i,j+1 Heat loss between the corresponding fluid infinitesimal element and the soil is recorded as the third partyUpward vector q3
Q is to be2Decomposed into vectors q2aAnd vector q2bAnd let q be2bAnd q is1Equal in size and opposite in direction, the simplified Ti,j+1Is Ti-1,jThe function of (a) is recorded as the recursion relation of the equivalent thermal circuit;
eliminating the temperature measurement of the intermediate node through the recursive relationship of the equivalent thermal circuit, establishing an equivalent thermal circuit equation and recording as a second function;
the second function includes equations (2.1) and (2.2):
Figure BDA0003049429870000221
Figure BDA0003049429870000231
the beneficial effect of this embodiment lies in:
a piecewise function of the outlet excess temperature (excess temperature, i.e., the difference between the fluid temperature and the soil temperature) is established with time, inlet temperature, and initial temperature as arguments. Due to An、BnThe model is a linear equation and can be directly used for electric heating coordination optimization. The problem of pipeline export temperature equation and electric heat coordination system incompatible in the prior art is solved.
According to any of the embodiments described above, in this embodiment:
the constraint condition further comprises a cogeneration unit constraint;
the constraint of the cogeneration unit comprises a thermoelectric ratio constraint, a climbing constraint, a unit limit output constraint, a heat source water supply temperature constraint and a scheduling period constraint;
the thermoelectric ratio constraints are:
Figure BDA0003049429870000232
in the formula, epsilon is a set thermoelectric ratio of the cogeneration unit; qCHP,tThe thermal output of the cogeneration unit at the moment t; pCHP,tThe electric output of the cogeneration unit at the moment t;
the climbing restriction is as follows:
Figure BDA0003049429870000233
in the formula, PCHP,up、PCHP,downRespectively restricting the upward and downward climbing of the set cogeneration unit; pCHP,t-1The power output of the cogeneration unit at the time of t-1; pthp,up、 Pthp,downRespectively performing upward and downward climbing constraints on the set thermal power cogeneration unit; pthp,t、 Pthp,t-1The thermal power cogeneration unit power output at the time t and the time t-1 respectively; wherein the t-1 moment is the last moment of the t moment;
the unit ultimate output restriction package is as follows:
Figure BDA0003049429870000241
in the formula, PCHP,max、PCHP,minThe maximum power output and the minimum power output of the cogeneration unit are respectively; pthp,max、Pthp,minRespectively the maximum and minimum power output of the thermal power cogeneration unit;
the heat source water supply temperature constraint is as follows:
40℃≤tsource,t≤90℃
in the formula, Tsource,tThe temperature of the heat source water supply at the time t;
the scheduling period constraint means that the thermal scheduling period is adjusted based on a set power scheduling period value, so that the electric heating coordination scheduling period, the thermal scheduling period and the power scheduling period are consistent and are set power scheduling period values.
The beneficial effect of this embodiment lies in:
the embodiment is different from the prior art, the optimization of the power system is not taken as a main body, the operation characteristics and the evaluation indexes of the thermodynamic system are comprehensively considered, and the compatibility of the electric heating coordination scheduling method considering the heat transmission dynamic characteristics of the heat supply network can be further improved.
According to any of the embodiments described above, in this embodiment:
the constraint conditions also comprise the constraints of a power grid and a wind generating set;
the power grid and wind power generator set constraints comprise power grid power balance constraints, wind power limit output constraints and line capacity constraints;
the power balance constraint of the power grid is as follows:
PCHP,t+Pwind,t=Puser,t+Ppump,t
in the formula, PCHP,tThe electric output of the cogeneration unit at the moment t; pwind,tThe actual wind power output at the moment t is obtained; puser,tThe load is the residential electricity load at the moment t; ppump,tThe power consumption of the circulating pump at the time t;
the wind power limit output constraint is as follows:
Pwind,t≤Pwind,max,t
in the formula, Pwind,max,tThe wind power limit output at the time t is obtained;
the line capacity constraints are:
Ptr,min≤Ptr,t≤Ptr,max
in the formula, Ptr,tThe transmission capacity of the power grid at the moment t; ptr,min、Ptr,maxThe minimum and maximum transmission capacities of the power grid at the moment t are respectively.
The beneficial effect of this embodiment lies in:
the embodiment is different from the prior art, the optimization of the power system is not taken as a main body, the operation characteristics and the evaluation indexes of the thermodynamic system are comprehensively considered, and the compatibility of the electric heating coordination scheduling method considering the heat transmission dynamic characteristics of the heat supply network can be further improved.
According to any of the embodiments described above, in this embodiment:
the scheduling condition is that the air abandoning amount at all the time in the whole day is minimum; the target variable corresponding to the scheduling condition with the minimum air abandon amount at all times in the whole day is the sum of the air abandon amounts, and the corresponding target value is 0.
The beneficial effect of this embodiment lies in:
electric heating coordination scheduling is carried out aiming at the abandoned wind quantity, and wind power consumption is maximized on the aspect of an electric power system.
According to any of the embodiments described above, in this embodiment:
the scheduling condition is that the absolute value of the total thermal deviation of the whole day is minimum, the target variable corresponding to the scheduling condition with the minimum absolute value of the total thermal deviation of the whole day is the sum of the absolute values of the thermal deviations, and the corresponding target value is 0.
The beneficial effect of this embodiment lies in:
and electric heat coordination scheduling is carried out according to the heat deviation, so that the overall thermal comfort of a user is promoted on the aspect of a thermodynamic system.
According to any of the embodiments described above, in this embodiment:
in order to quantitatively evaluate the heat storage performance of the heat supply network in the electric heating integrated energy system and improve the effect of renewable energy consumption, an electric heating coordination optimization scheduling model considering the heat transmission dynamic characteristics of the heat supply network needs to be established.
At present, a common semi-theoretical semi-empirical model of a heat supply network in an electric heating coordination optimization model cannot accurately calculate the dynamic temperature change rule of the heat supply network and cannot obtain the accurate mapping relation between the heat supply temperature of a heat source and the indoor temperature of a heat user; it is difficult to provide accurate constraints for the heating plan and the power supply plan, and large deviation is brought to quantitative evaluation of the effect of improving renewable energy sources of the heat supply network.
In view of the above problems, the present embodiment provides an electric heating coordination optimization scheduling model considering dynamic characteristics of a heat supply network based on an existing heat supply network linearization matrix model. The optimization model comprehensively considers the operation characteristics of a heat supply network, a heat exchanger closely related to the operation of the heat supply network, a heat user and the like, comprehensively considers the operation evaluation indexes of an electric power system and a thermodynamic system by taking the maximum optimal thermal comfort and wind power consumption as a two-stage optimization target, and provides a basis for analyzing and evaluating the influence of the thermal inertia of the heat supply network on the electric heat coordination optimization operation effect.
Finally, a typical change example and a comparative example are selected, a heat supply plan and a power generation plan of the electric heating comprehensive energy system are optimized in real time and quantitatively according to heat transmission characteristics of a heat supply network, and a conclusion that influences of dynamic characteristics of heat transmission of the heat supply network on wind power consumption and thermal comfort are taken into consideration is obtained.
Specifically, the embodiment provides a linearization model considering the heat transfer dynamic characteristic of the heat supply network on the basis of the existing heat supply network linearization model, and has the requirements of precision and compatibility; and physical characteristic models of thermodynamic system components closely related to the operation of a water supply network, such as a primary network and a secondary network, a secondary network and a heat consumer, a heat consumer and a three-level heat exchange model of an outdoor environment and the like, are established. On the basis, an electric heating coordination optimization scheduling model considering a heat transmission dynamic model of a heat supply network and a heat exchange process is constructed, the hourly heat supply temperature of a heat source is taken as a decision variable, and the maximum optimal thermal comfort and wind power consumption is an optimal target. And finally, selecting a typical variation example and a comparative example, and quantitatively researching the influence of the operation characteristics of the heat supply network on the electric heating coordination optimization operation effect.
The heat transmission and heat exchange process in the pipe network comprises five processes: the heat transmission process comprises heat transmission of a water supply pipeline, heat exchange of a primary network and a secondary network, heat exchange of the secondary network and a user, heat exchange of the user and an outdoor environment and heat transmission of a water return pipeline. In the embodiment, a step-by-step heat exchange equation is mainly established for each node, and a pipeline heat transmission-based pipeline network thermodynamic dynamic model is established.
The calculation process of the thermodynamic dynamic model of the whole pipe network comprises the following steps:
writing all pipe sections as an expression of the outlet temperature with respect to the inlet temperature and the heat transfer contribution matrix; then listing a heat balance equation of the pipe sections connected end to end at the intersection node and a mixed temperature equation of the mixed node column; finally, a heat transmission equation of each pipe section is combined with a heat balance and mixing temperature equation, and the inlet temperature (namely the heat supply temperature of a heat source) of the pipe network and the flow rate of the heat medium water are given; the thermodynamic characteristics such as the temperature dynamic change rule of each user node, the pipeline section along the way and the like can be obtained, and the expression of the characteristics is shown in an expression (2.1) and an expression (2.2).
The heat is transferred from the primary network heat medium water to the users, which relates to three heat exchange processes, namely, a primary network and a secondary network, a secondary network heat medium water and heat users, and a heat user and an external environment (neglecting the heat exchange process of the secondary network and the environment). The heat balance and heat exchange model comprising three layers of heat exchange is established aiming at the problem of multi-loop heat exchange by comprehensively considering the three-layer heat exchange process in the chapter, and the expression of the indoor temperature about the primary network heat supply temperature and the outdoor temperature is directly obtained.
According to the schematic diagram of the multi-loop heat exchange process shown in FIG. 2, the heat exchange equations of three heat exchange loops, i.e., the equations (1a) to (1h), are listed in sequence.
After the solutions of the joint vertical type (1a) to the formula (1h) are solved and simplified, the indoor temperature T can be obtaineduser,tAbout the temperature T of a heat medium water inlet of a primary network of a heat exchange stationh,in,tOutdoor atmospheric temperature Toutdoor,tThe relationship with the heat exchange parameters of the heat exchanger of the heat exchange station, the heat exchange parameters of the radiator and the heat exchange parameters of the user enclosure structure is shown in the formula (1.1), the formula (1.2) and the formula (1.3)
Electric heating coordination optimization modeling considering the overall operation of the pipe network and the system is performed below.
By taking the problem of wind abandonment in wind power consumption in winter in northeast China as a background, a regional electric heating comprehensive energy system comprising a centralized heating system, a wind power plant, a thermal power generating unit and a cogeneration unit is taken as an example, and an electric heating coordination optimization model considering the overall operation of a pipe network and a system is established. The optimization scheduling model does not only take power system optimization as a main body, but comprehensively considers the operating characteristics and evaluation indexes of the thermodynamic system.
Firstly, the regional electric heating integrated energy system and the electric heating coordination optimization model are summarized as follows.
Fig. 6 is a structure of a regional electric heating integrated energy system case in a certain area, which mainly comprises a generator set, a power grid, a heat supply pipe network, a heat exchanger, an electric user, a heat user and the like, wherein the generator set comprises a gas turbine cogeneration unit, a thermal power unit and a wind power unit, and the sum of the electric power output is equal to the electric load of the user at any moment so as to ensure the power balance of the power system. The thermal load of the residents is provided by the cogeneration unit.
The electric heating coordination optimization model is mainly composed of three elements of operation constraint, optimization target and decision variable. Because the power system and the heating system are coupled together through the cogeneration unit, the operation constraint and decision variable elements are divided into three types, namely a power system, a heating system and an electric-heating coupling. The operation constraint is composed of operation constraints of a power grid, a thermal power generating unit, a wind power generating unit and the like in the power system, pipe network heat transmission constraints, heat exchange constraints and cogeneration unit constraints in the heat supply system, and the decision variables comprise hourly heat source heat supply temperature, thermoelectric unit heat output, thermoelectric unit output, thermal power generating unit output and wind power output. The constraint and decision variables of the cogeneration unit belong to both an electric power system and a heat supply system, and are an electric-thermal coupling part. The optimization target is a two-stage optimization target considering the minimum wind power curtailment amount all day long and the best user thermal comfort.
The system operation constraints of the present embodiment will be further explained below.
The electric heat comprehensive energy system comprises a gas turbine cogeneration unit, a wind power plant, a centralized heating system, a power grid and the like, and the operation constraint of each part is specifically described as follows.
The constraint conditions comprise the constraint of a cogeneration unit, the constraint of a power grid and a wind generating unit besides the heat transmission constraint, namely a first function and a second function;
in the embodiment, the power dispatching system regulates the output of each power plant unit every 15 minutes, while the thermal power pipe network regulation period in the prior art is longer than the power dispatching period, and usually the temperature is regulated for tens of minutes or even hours. The key point of this embodiment is to incorporate regulation of the heat pipe network into the power scheduling framework, so in order to unify the time scale of operation regulation of each part of the integrated electric heating energy system, this embodiment establishes a heat transmission equation of a 15-minute regulation period for the heat pipe network.
In the embodiment, the transmission capacity of the related lines is assumed to be within the limit value, and the influence of the transmission capacity of the power grid on the regulation of the thermodynamic system and the coordinated operation of the electric heating comprehensive energy system is not considered.
The optimization objectives and decision variables of the present embodiment will be explained below.
In order to realize the optimal operation of the heating power system and the electric power system, the system operation indexes of the heating power system and the electric power system need to be considered. And selecting the total thermal comfort degree of the user at each moment as an evaluation index from the aspect of a thermodynamic system, and representing the degree of the thermal comfort degree by adopting the thermal deviation. And selecting the maximum wind power total consumption at each moment as an optimization target from the aspect of the power system, namely selecting the minimum wind power total wind abandon amount as the optimization target. The minimum air loss amount at all times of the whole day is taken as a primary optimization target, and the minimum absolute value of the total heat deviation amount of the whole day is taken as a secondary optimization target, which is shown in the following formula.
Figure BDA0003049429870000291
Wherein i represents the node number, Pi curw,tShowing the air curtailment quantity of the ith node at the time T, Ti dev,tThe thermal deviation of the ith node at the moment t is represented, alpha is a weight coefficient, and the value is a minimum value, such as 0.001, and the aim is to convert into an optimal optimization problem only considering thermal comfort when the air abandoning amount is 0.
The decision variables are the hourly heat source heat supply temperature, the heat supply amount and the electric output of the thermoelectric unit, the output of the thermal power unit and the hourly wind power output.
The method is an electric heating coordination optimization model considering the overall operation of the pipe network and the system, and the heat transmission dynamic characteristic of the pipe network, the temperature feedback regulation mechanism of the pipe network, the operation characteristic of the heat exchanger and the operation constraint of the power system are comprehensively considered. The model can provide an accurate operation boundary for an electric heating comprehensive energy system, and realizes system coordination optimization operation considering performance indexes of an electric power system and a thermodynamic system.
The following describes the operation effect of the coordinated electric heating optimization considering the dynamic characteristics of the heat supply network, and the beneficial effect compared with the prior art, with reference to the variation of the present embodiment and the comparative example of the prior art.
The method comprises the steps that a user node which is 1.26km away from a heat source is selected as a research object, a power generation system is composed of 2 gas turbine cogeneration units (one for one), and a wind power plant, a time-by-time output curve of the wind power unit adopts data in documents, the maximum power generation power and the minimum power generation power of the cogeneration units are 23MW and 8MW respectively, and the climbing power is 10 MW/h.
The heat exchange area and coefficient of the building envelope structure of the user node and the area and coefficient of the heat exchanger are shown in tables 1 and 2.
Table 1: single user node user side heat exchange parameter table
Figure BDA0003049429870000301
Table 2: single user node user side heat exchange parameter table
Figure BDA0003049429870000302
And taking historical electric load data of a certain area obtained by investigation as original data. The time-by-time electrical load curve and the wind power limit output curve of the user are shown in figure 7.
The outdoor temperature is taken as an example of data of 2016, 11, 16, 24 hours all day in Beijing, and the outdoor temperature curve of every hour all day is shown in figure 8.
The length of a pipeline from a heat source to a single-user node is 1.26km, the pipe diameter is 0.5m, the flow speed is 0.68m/s, the grid ratio is 0.25, the space node step length is 24m, and the time step length is 6 s. According to the 3.3-section heat transfer formula, a composite heat transfer contribution matrix of the user node can be obtained, wherein B represents an influence coefficient matrix of the inlet temperature on the outlet temperature, and a represents an influence coefficient matrix of the initial temperature on the outlet temperature, as shown in the following formula.
Figure BDA0003049429870000311
On this basis, this example provides 3 comparative examples for analysis:
comparative example 1, namely scheme 1, carries out electric heating coordinated scheduling based on the embodiment, and considers a pipe network temperature feedback regulation mechanism and a pipe network heat transmission accurate model, namely comparative example 1 establishes heat transmission constraint based on a multistage heat exchange process, and simultaneously considers the constraint of a first function and a second function on electric heating coordinated scheduling. The electric heating coordinated scheduling of the comparative example 1 takes the maximum wind power total consumption as a primary optimization target, takes the best thermal comfort as a secondary optimization target, and takes the lowest indoor temperature of more than or equal to 20 ℃ at each moment as an operation constraint.
Comparative example 2, namely scheme 2, is based on the prior art (Liu x., Jenkins n., Wu J., et al. Combined analysis of electric and heat networks [ J ]. Energy procedure, 2014, 61: 155-.
Comparative example 3, namely scheme 3, does not consider the flexible fluctuation of the transmission dynamic characteristic of the heat medium water in the pipeline and the thermal comfort of users, but considers the feedback adjustment of the temperature, and the adjustment standard is that the indoor temperature reaches the optimal value of 20 ℃.
The total air loss in all days of each scheme obtained by the optimization calculation is shown in table 3, and it can be seen that the total air loss is as follows from large to small: comparative example 3, comparative example 2, comparative example 1. The data in the table can provide reference for subsequent analysis.
Table 3: all-day total air abandon rate of each proportion
Figure BDA0003049429870000321
Fig. 9 and 10 are graphs of the heating temperature of the heat source and the indoor temperature of comparative example 1, comparative example 2, and comparative example 3. The three schemes all consider a pipe network temperature feedback regulation mechanism. And the influence of the heat transmission characteristic of the pipe network and the accurate modeling on the system operation is considered in the comparative analysis. As can be seen from the figure, the variation trends of the heating temperatures of all the schemes are basically consistent, namely the heating temperature at night is higher than that at daytime, because the temperature at night is low, in order to ensure the heat load of the user, the heating temperature at night is increased. Comparing the heat supply temperatures under different schemes, wherein the heat supply temperature of the comparative example 1 is obviously lower than that of the comparative example 2 and that of the comparative example 3 from 1 morning to 8 morning and from 14 afternoon to 16 afternoon, because the wind power limit output is larger in the above time period, the wind power is even larger than the total power load at night, the maximum wind power consumption of the comparative example 1 is a first-level optimization target, and in order to maximally consume the wind power, the heat supply temperature is correspondingly adjusted down in the time period, and the thermoelectric output is reduced; in a period of time between 1 am and 13 pm, the heating temperature of the comparative example 1 is significantly higher than that of the other two schemes, and the purpose is to improve the heating temperature in advance by using the transmission delay characteristic of the heat supply network so as to meet the heat load in the subsequent period of time (avoid insufficient heat load caused by reducing the heating load in the subsequent period of time); in other periods, the heating temperature of comparative example 1 was slightly lower than that of the other solutions, but the difference was small and was substantially equal. The results are in accordance with the scheduling principles of comparative example 1. Comparing the heating temperatures of comparative example 2 and comparative example 3, it can be seen that the heating temperature sequence of comparative example 2 is out of phase with comparative example 3, i.e., the heating temperature at the previous moment of comparative example 2 is approximately equal to the temperature at the moment of comparative example 3, because comparative example 2 considers the heat transfer delay and the heating amount is provided before a transfer delay time.
As can be seen from the indoor temperature graph 10, the room temperature in the initial period from the early 7 th to the early 9 th of the scheme 1 is low, because the heat source temperature is not transmitted to the heat user, and the initial temperature of the pipe network is not enough to make the room temperature reach the expected target; the indoor temperature is higher than 20 c from about 11 to 15 c and from night 22 to early morning 1c, because the wind power is about to increase in the period after these times, so increasing the heat supply amount in advance can suitably reduce the thermal output in the period of wind power generation, so the heat supply temperature at these times is higher, resulting in the room temperature exceeding the expected 20 c, and the conclusion is in agreement with fig. 9. At other times, the room temperature of comparative example 1 was stabilized at 20 ℃. The room temperature of comparative examples 2 and 3 was lower than 20 ℃ at a higher time than that of comparative example 1. For comparison 2, the reason why the temperature is frequently lower than 20 ℃ is that the pipe network model is not accurate enough, and the pipe network model is difficult to be directly incorporated into the electric heating coordination optimization for calculation, and an accurate operation boundary cannot be provided. The reason why the room temperature is frequently lower than 20 c for comparative example 3 is that the feedback regulation mechanism has a certain hysteresis, and a relatively stable room temperature cannot be secured by quantitatively optimizing the time-by-time heating temperature in advance as in comparative example 1. The total sum of absolute values of all-day thermal deviation of the three schemes of the comparative example 1, the comparative example 2 and the comparative example 3 is 31.91, 30.46 and 37.39 respectively, and the lower the total sum is, the better the thermal comfort is, so the thermal comfort of the comparative example 3 is the worst, and the comparative example 1 and the comparative example 2 are basically consistent, which shows that the consideration of the transmission characteristic of the heat network is beneficial to improving the thermal comfort of users to a certain extent.
Fig. 11 reflects the nighttime wind abandoning curves of comparative example 1, comparative example 2 and comparative example 3, only the wind abandoning data from 1 to 5 in the morning are selected, and the wind abandoning effects of different schemes at the nighttime wind power outbreak period are mainly analyzed because the differences of the wind abandoning amounts of the three schemes in other periods are not large. It can be seen that the hourly air rejection rate of comparative example 1 is much lower than that of comparative examples 2 and 3 at most times, and the air rejection rates of comparative examples 2 and 3 are alternately increased at each period. According to the table 3, the total daily air loss of the comparative example 1 is 167.57MW & h, which is reduced by 5.9% compared with the comparative example 3(178.10 MW & h), and the effect of promoting wind power consumption by considering the transmission characteristic of a pipe network in electric heat coordination optimization is shown; the total air curtailment rate is reduced by 4.9% compared with the ratio 2(176.28MW & h), which indicates that the accurate pipe network heat transmission dynamic model adopted in the comparative example 1 can promote wind power consumption better than the accurate model adopted in the ratio 2, and the accurate pipe network model can provide an accurate heat supply operation boundary for electric heat coordination optimization, so that the quantitative and coordination optimization of system operation is realized, and the aim of reducing the air curtailment rate is fulfilled.
From the above comparative example analysis, the following conclusions can be drawn.
The influence of each element such as pipe network characteristics on the electric heating coordination optimization operation effect is quantitatively evaluated. Taking a single hot user as an example, the calculation result shows that: the air volume of the system is reduced by 5.9% by considering the heat transmission characteristic of the pipe network.
The total air curtailment of the comparative example 1 in the whole day is 167.57MW & h, which is reduced by 5.9% compared with that of the comparative example 3(178.10 MW & h), and the effect of promoting wind power consumption can be achieved by considering the transmission characteristic of a pipe network in the electric heating coordination optimization; the total air curtailment rate is reduced by 4.9% compared with the ratio 2(176.28MW & h), which indicates that the accurate pipe network heat transmission dynamic model adopted in the comparative example 1 can promote wind power consumption better than the accurate model adopted in the ratio 2, and the accurate pipe network model can provide an accurate heat supply operation boundary for electric heat coordination optimization, so that the quantitative and coordination optimization of system operation is realized, and the aim of reducing the air curtailment rate is fulfilled.
The beneficial effect of this embodiment lies in:
an electric heating coordination optimization model considering the overall operation of the pipe network and the system is provided, and the heat transmission dynamic characteristic of the pipe network, the temperature feedback regulation mechanism of the pipe network, the operation characteristic of a heat exchanger and the operation constraint of the power system are comprehensively considered. The model can provide an accurate operation boundary for the electric heating comprehensive energy system so as to realize the coordination optimization calculation simulation considering the performance indexes of the electric power and the thermodynamic system.
The electric heating coordination scheduling device considering the heat transmission dynamic characteristic of the heat supply network provided by the invention is described below, and the electric heating coordination scheduling device considering the heat transmission dynamic characteristic of the heat supply network described below and the electric heating coordination scheduling method considering the heat transmission dynamic characteristic of the heat supply network described above can be referred to correspondingly.
The embodiment of the invention also provides an electric heating coordination scheduling device considering the heat transmission dynamic characteristics of the heat supply network, which comprises a model establishing module, a constraint condition setting module, a calculating module, a scheduling module and a feedback adjusting module;
the model establishing module can establish an electric heating coordination scheduling model;
the constraint condition setting module can establish heat transmission constraint based on a multi-stage heat exchange process to serve as a constraint condition;
the calculation module can obtain the current value of the decision variable, and according to the current value of the decision variable and the constraint condition, the current value of the target variable is calculated and obtained through an electric heating coordination scheduling model;
the scheduling module can determine a target interval of a target variable or a target value of the target variable based on a set scheduling condition;
the feedback regulation can feed back and regulate decision variables according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, so that the electric heating coordinated dispatching is realized;
the electric heating coordination scheduling model is a model for adjusting a target variable through a decision variable under a set constraint condition;
the decision variables comprise any one or any combination of more of hourly heat source heat supply temperature, heat supply pipe network initial temperature, heat supply quantity of the thermoelectric unit, electric output, thermal power unit output and hourly wind power output;
the target variable comprises any one or any combination of more of air abandonment amount, heat deviation amount, heat supply amount, heat loss, user heat gain amount, user indoor temperature and wind power consumption amount;
the multistage heat exchange process comprises a heat transmission process and/or a heat exchange process of a heat supply pipe network; the heat transmission process of the heat supply pipe network comprises a heat transmission process of a heat medium supply pipeline and/or a heat transmission process of a heat medium recovery pipeline; the heat exchange process comprises any one or any combination of a heat exchange process between a primary network and a secondary network, a heat exchange process between the secondary network and a user, and a heat exchange process between the user and an outdoor environment.
The beneficial effect of the embodiment is that
The heat transmission constraint is established based on the multi-stage heat exchange process, the heat transmission and conversion process of each node in the thermodynamic system is reflected, and the heat exchange constraint between the primary network/the secondary network/the heat users is established, so that the consideration factors of the electric heat coordination scheduling model are refined, the model optimization result is more reliable, and the model compatibility is better.
Fig. 12 illustrates a physical structure diagram of an electronic device, which may include, as shown in fig. 12: a processor (processor)1210, a communication Interface (Communications Interface)1220, a memory (memory)1230, and a communication bus 1240, wherein the processor 1210, the communication Interface 1220, and the memory 1230 communicate with each other via the communication bus 1240. Processor 1210 may invoke logic instructions in memory 1230 to perform a method of thermoelectric coordinated scheduling that accounts for thermal network thermal transfer dynamics, the method comprising: establishing an electric heating coordination scheduling model; establishing heat transmission constraint based on a multi-stage heat exchange process as a constraint condition; obtaining the current value of a decision variable, and calculating through an electric heating coordinated scheduling model according to the current value of the decision variable and a constraint condition to obtain the current value of a target variable; determining a target interval of a target variable or a target value of the target variable based on the set scheduling condition; and feeding back and adjusting the decision variable according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, thereby realizing the electric heating coordinated scheduling.
In addition, the logic instructions in the memory 1230 may be implemented in software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as a stand-alone product. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, the present invention also provides a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions, which when executed by a computer, enable the computer to perform a method of electric heating coordinated scheduling taking into account thermal transport dynamics of a heat grid, the method comprising: establishing an electric heating coordination scheduling model; establishing heat transmission constraint based on a multi-stage heat exchange process as a constraint condition; obtaining the current value of a decision variable, and calculating through an electric heating coordinated scheduling model according to the current value of the decision variable and a constraint condition to obtain the current value of a target variable; determining a target interval of a target variable or a target value of the target variable based on the set scheduling condition; and feeding back and adjusting the decision variable according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, thereby realizing the electric heating coordinated scheduling.
In yet another aspect, the present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program, which when executed by a processor, implements the method for coordinating and scheduling electric heat and heat in consideration of heat transfer dynamics of a heat supply network, the method comprising: establishing an electric heating coordination scheduling model; establishing heat supply network dynamic constraints based on a heat supply network linearization heat transfer model and a multi-stage heat exchange process, wherein the heat supply network dynamic constraints are used as constraint conditions; obtaining the current value of a decision variable, and calculating through an electric heating coordinated scheduling model according to the current value of the decision variable and a constraint condition to obtain the current value of a target variable; determining a target interval of a target variable or a target value of the target variable based on the set scheduling condition; and feeding back and adjusting the decision variable according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, thereby realizing the electric heating coordinated scheduling.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; 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 (10)

1. An electric heating coordination scheduling method considering heat transmission dynamic characteristics of a heat supply network is characterized by comprising the following steps:
establishing an electric heating coordination scheduling model;
establishing heat supply network constraint based on the heat supply network heat transmission dynamic characteristic and the multistage heat exchange process as constraint conditions;
obtaining the current value of a decision variable, and calculating through an electric heating coordinated scheduling model according to the current value of the decision variable and a constraint condition to obtain the current value of a target variable;
determining a target interval of a target variable or a target value of the target variable based on the set scheduling condition;
feeding back and adjusting the decision variable according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, thereby realizing electric heating coordinated scheduling;
the electric heating coordination scheduling model is a model for adjusting a target variable through a decision variable under a set constraint condition;
the decision variables comprise any one or any combination of more of hourly heat source heat supply temperature, heat supply pipe network initial temperature, heat supply quantity of the thermoelectric unit, electric output, thermal power unit output and hourly wind power output;
the target variable comprises any one or any combination of more of air abandonment amount, heat deviation amount, heat supply amount, heat loss, user heat gain amount, user indoor temperature and wind power consumption amount;
the multistage heat exchange process comprises a heat transmission process and/or a heat exchange process of a heat supply pipe network; the heat transmission process of the heat supply pipe network comprises a heat transmission process of a heat medium supply pipeline and/or a heat transmission process of a heat medium recovery pipeline; the heat exchange process comprises any one or any combination of a heat exchange process between a primary network and a secondary network, a heat exchange process between the secondary network and a user, and a heat exchange process between the user and an outdoor environment.
2. The method of claim 1, wherein the heat supply network constraints comprise a first function corresponding to a heat exchange process:
Figure FDA0003049429860000021
Figure FDA0003049429860000022
Figure FDA0003049429860000023
in the formula:
Tc,in,tthe inlet temperature of the heat medium of the secondary network at the time t; t isoutdoor,tIs the outdoor atmospheric temperature at time t; t ish,in,tThe inlet temperature of the heat medium of the primary network at the time t;
Th,out,tthe outlet temperature of the primary network heat medium at the moment t; gcMass flow of the secondary network heat medium; ghMass flow rate of the primary net heat medium;
Tuser,tthe indoor temperature of the user at the moment t;
sh,c、sevl、sc,user、s12hare respectively:
Figure FDA0003049429860000031
Figure FDA0003049429860000032
Figure FDA0003049429860000033
wherein e is the base of the natural logarithm; k is a radical ofh,cThe heat exchange coefficient of the primary network heat exchanger and the secondary network heat exchanger is obtained; a. theh,cThe heat exchange area of the primary net and the secondary net heat exchanger is shown; c is the specific heat capacity of the heat medium; k is a radical ofevlAverage heat exchange coefficient of the building envelope for the user; a. theevlAverage heat exchange area of the building envelope for the user; k is a radical ofc,userThe heat exchange coefficient of the radiator between the secondary network and the user is shown; a. thec,userThe heat exchange area of the radiator between the secondary network and the user.
3. The method of claim 1 or 2, wherein the heat supply network constraints comprise a second function corresponding to a heat supply network heat transfer process:
Figure FDA0003049429860000034
Figure FDA0003049429860000035
in the formula, phim,nThe outlet excess temperature is the outlet excess temperature with a space node of m and a time node of n; i. j is a space node serial number and a time node serial number respectively; the surplus temperature refers to the difference between the temperature of the heat medium and the temperature of the soil;
parameter An、BnAnd A'nSatisfies the following conditions:
Figure FDA0003049429860000041
Figure FDA0003049429860000042
Figure FDA0003049429860000043
wherein, C is a permutation and combination operator;
the parameters a and r satisfy:
Figure FDA0003049429860000044
Figure FDA0003049429860000045
in the formula, rho and c are respectively the density of the heat medium and the specific heat capacity of the heat medium; A. d is the cross-sectional area and the diameter of the pipeline respectively; v is the fluid flow rate; k is the heat leakage coefficient between the fluid in the pipeline and the soil; x and tau are respectively a space variable and a time variable.
4. The method of claim 1, wherein the constraints further include cogeneration unit constraints;
the constraint of the cogeneration unit comprises a thermoelectric ratio constraint, a climbing constraint, a unit limit output constraint, a heat source water supply temperature constraint and a scheduling period constraint;
the thermoelectric ratio constraints are:
Figure FDA0003049429860000051
in the formula, epsilon is a set thermoelectric ratio of the cogeneration unit; qCHP,tThe thermal output of the cogeneration unit at the moment t; pCHP,tThe electric output of the cogeneration unit at the moment t;
the climbing restriction is as follows:
Figure FDA0003049429860000052
in the formula, PCHP,up、PCHP,downRespectively restricting the upward and downward climbing of the set cogeneration unit; pCHP,t-1The power output of the cogeneration unit at the time of t-1; pthp,up、Pthp,downRespectively performing upward and downward climbing constraints on the set thermal power cogeneration unit; pthp,t、Pthp,t-1The thermal power cogeneration unit power output at the time t and the time t-1 respectively; wherein the t-1 moment is the last moment of the t moment;
the unit ultimate output restriction package is as follows:
Figure FDA0003049429860000053
in the formula, PCHP,max、PCHP,minThe maximum power output and the minimum power output of the cogeneration unit are respectively; pthp,max、Pthp,minRespectively the maximum and minimum power output of the thermal power cogeneration unit;
the heat source water supply temperature constraint is as follows:
40℃≤Tsource,t≤90℃
in the formula, Tsource,tThe temperature of the heat source water supply at the time t;
the scheduling period constraint means that the thermal scheduling period is adjusted based on a set power scheduling period value, so that the electric heating coordination scheduling period, the thermal scheduling period and the power scheduling period are consistent and are set power scheduling period values.
5. The method of claim 1, wherein the constraints further include grid and wind farm constraints;
the power grid and wind power generator set constraints comprise power grid power balance constraints, wind power limit output constraints and line capacity constraints;
the power balance constraint of the power grid is as follows:
PCHP,t+Owind,t=Ouser,t+Opump,t
in the formula, OCHP,tThe electric output of the cogeneration unit at the moment t; pwind,tThe actual wind power output at the moment t is obtained; puser,tThe load is the residential electricity load at the moment t; ppump,tThe power consumption of the circulating pump at the time t;
the wind power limit output constraint is as follows:
Pwind,t≤Pwind,max,t
in the formula, Pwind,max,tWind power limit at time tForce is exerted;
the line capacity constraints are:
Ptr,min≤Ptr,t≤Ptr,max
in the formula, Ptr,tThe transmission capacity of the power grid at the moment t; ptr,min、Ptr,maxThe minimum and maximum transmission capacities of the power grid at the moment t are respectively.
6. The electric heat coordination scheduling method considering heat transmission dynamic characteristics of the heat supply network according to claim 1, wherein the scheduling condition is that air abandonment amount at all time of the whole day is minimum; the target variable corresponding to the scheduling condition with the minimum air abandon amount at all times in the whole day is the sum of the air abandon amounts, and the corresponding target value is 0.
7. The electric heating coordinated scheduling method considering heat transmission dynamic characteristics of the heat supply network according to claim 1, wherein the scheduling condition is that the absolute value of total thermal deviation of the whole day is minimum, the target variable corresponding to the scheduling condition with the minimum absolute value of total thermal deviation of the whole day is the sum of the absolute values of the thermal deviation, and the corresponding target value is 0.
8. An electric heating coordination scheduling device considering heat transmission dynamic characteristics of a heat supply network is characterized by comprising a model establishing module, a constraint condition setting module, a calculating module, a scheduling module and a feedback adjusting module;
the model establishing module can establish an electric heating coordination scheduling model;
the constraint condition setting module can establish heat transmission constraint based on a multi-stage heat exchange process to serve as a constraint condition;
the calculation module can obtain the current value of the decision variable, and according to the current value of the decision variable and the constraint condition, the current value of the target variable is calculated and obtained through an electric heating coordination scheduling model;
the scheduling module can determine a target interval of a target variable or a target value of the target variable based on a set scheduling condition;
the feedback regulation can feed back and regulate decision variables according to the relative relation between the current value of the target variable and the target interval or according to the difference value between the current value of the target variable and the target value, so that the electric heating coordinated dispatching is realized;
the electric heating coordination scheduling model is a model for adjusting a target variable through a decision variable under a set constraint condition;
the decision variables comprise any one or any combination of more of hourly heat source heat supply temperature, heat supply pipe network initial temperature, heat supply quantity of the thermoelectric unit, electric output, thermal power unit output and hourly wind power output;
the target variable comprises any one or any combination of more of air abandonment amount, heat deviation amount, heat supply amount, heat loss, user heat gain amount, user indoor temperature and wind power consumption amount;
the multistage heat exchange process comprises a heat transmission process and/or a heat exchange process of a heat supply pipe network; the heat transmission process of the heat supply pipe network comprises a heat transmission process of a heat medium supply pipeline and/or a heat transmission process of a heat medium recovery pipeline; the heat exchange process comprises any one or any combination of a heat exchange process between a primary network and a secondary network, a heat exchange process between the secondary network and a user, and a heat exchange process between the user and an outdoor environment.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the method of any of claims 1 to 7 for electro-thermal coordinated scheduling taking into account thermal transport dynamics of a heat network.
10. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the method of any of claims 1 to 7 for electric heat coordinated scheduling taking into account thermal transport dynamics of a heat network.
CN202110481506.9A 2021-04-30 Electrothermal coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network Active CN113191638B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110481506.9A CN113191638B (en) 2021-04-30 Electrothermal coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110481506.9A CN113191638B (en) 2021-04-30 Electrothermal coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network

Publications (2)

Publication Number Publication Date
CN113191638A true CN113191638A (en) 2021-07-30
CN113191638B CN113191638B (en) 2024-10-29

Family

ID=

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114970987A (en) * 2022-05-11 2022-08-30 昆明理工大学 Electric heating comprehensive energy system optimization scheduling method and system based on thermoelectric simulation model
CN116207796A (en) * 2023-03-27 2023-06-02 山东大学 Electric-thermal comprehensive energy system scheduling method and system considering variable flow regulation

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180356774A1 (en) * 2017-02-22 2018-12-13 Tsinghua University Dispatch method and apparatus for combined heat and power system
CN109727158A (en) * 2019-01-25 2019-05-07 燕山大学 A kind of electric heating integrated energy system dispatching method based on the weak robust optimization of improvement
CN110007600A (en) * 2019-04-24 2019-07-12 四川大学 A kind of constrained multipotency stream coordinated scheduling aid decision-making system of tool
CN112434861A (en) * 2020-11-26 2021-03-02 清华大学 Electric heating coordination optimization scheduling method and device considering feedback adjustment of heating power pipe network

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180356774A1 (en) * 2017-02-22 2018-12-13 Tsinghua University Dispatch method and apparatus for combined heat and power system
CN109727158A (en) * 2019-01-25 2019-05-07 燕山大学 A kind of electric heating integrated energy system dispatching method based on the weak robust optimization of improvement
CN110007600A (en) * 2019-04-24 2019-07-12 四川大学 A kind of constrained multipotency stream coordinated scheduling aid decision-making system of tool
CN112434861A (en) * 2020-11-26 2021-03-02 清华大学 Electric heating coordination optimization scheduling method and device considering feedback adjustment of heating power pipe network

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114970987A (en) * 2022-05-11 2022-08-30 昆明理工大学 Electric heating comprehensive energy system optimization scheduling method and system based on thermoelectric simulation model
CN116207796A (en) * 2023-03-27 2023-06-02 山东大学 Electric-thermal comprehensive energy system scheduling method and system considering variable flow regulation
CN116207796B (en) * 2023-03-27 2023-10-27 山东大学 Electric-thermal comprehensive energy system scheduling method and system considering variable flow regulation

Similar Documents

Publication Publication Date Title
JP7261507B2 (en) Electric heat pump - regulation method and system for optimizing cogeneration systems
CN107301472B (en) Distributed photovoltaic planning method based on scene analysis method and voltage regulation strategy
CN105375479B (en) A kind of distributed energy energy management method based on Model Predictive Control
CN113190999B (en) Electric heating coordination method and device for improving wind power consumption by adjusting flow of heat supply pipe network
CN108321793B (en) Active power distribution network modeling and optimal scheduling method integrating flexible loads of intelligent building
KR102194002B1 (en) energy operation management system by using optimized physical learning model and machine learning methods
CN113190942B (en) Method and device for calculating virtual energy storage capacity of heat supply/cold system and electronic equipment
CN110707755B (en) Comprehensive energy system ultra-short-term scheduling method based on energy hub under consideration of non-ideal communication condition
Zhang et al. Combined heat and power scheduling: Utilizing building‐level thermal inertia for short‐term thermal energy storage in district heat system
CN110489915B (en) Electric-heat combined scheduling method and system considering comprehensive demand response
CN113159983B (en) Coordinated control method for off-grid phosgene combined heat and power generation comprehensive energy system
CN118199060B (en) Low-voltage flexible interconnection load balancing regulation and control system for distribution transformer
CN111898856B (en) Analysis method of physical-data fusion building based on extreme learning machine
CN111400937A (en) Load flow calculation method of comprehensive energy system
CN114048908A (en) Multi-time scale strengthening interval optimization method for unified power grid-distributed heat supply network system
CN111523697A (en) Comprehensive energy service cost allocation and pricing calculation method
CN114936762A (en) Comprehensive energy system expansion planning method considering flexible electric load
CN113869742A (en) Power dispatching system of comprehensive supply and demand side based on mobile home and critic network
CN113191638A (en) Electric heating coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network
Patteeuw et al. Demand response for residential heat pumps in interaction with the electricity generation system
CN114661011B (en) Building comprehensive energy system optimizing operation method, device and storage medium
CN113191638B (en) Electrothermal coordination scheduling method and device considering heat transmission dynamic characteristics of heat supply network
CN115758763A (en) Multi-energy flow system optimal configuration method and system considering source load uncertainty
CN114169620B (en) Consistency optimization method of electric-gas storage-gas coupling system considering line loss
Wu et al. Neural Network Based Fea sible Region Approximation Model for Optimal Operation of Integrated Electricity and Heating System

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant