CN114997457A - Flexibility optimization scheduling method for cogeneration unit based on heat storage characteristic of centralized heating system - Google Patents

Flexibility optimization scheduling method for cogeneration unit based on heat storage characteristic of centralized heating system Download PDF

Info

Publication number
CN114997457A
CN114997457A CN202210286544.3A CN202210286544A CN114997457A CN 114997457 A CN114997457 A CN 114997457A CN 202210286544 A CN202210286544 A CN 202210286544A CN 114997457 A CN114997457 A CN 114997457A
Authority
CN
China
Prior art keywords
heat
temperature
chp
unit
output
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210286544.3A
Other languages
Chinese (zh)
Inventor
戴梦迪
柯德平
郑景文
宋琳
徐箭
郭雨
吴笑民
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd
Original Assignee
Wuhan University WHU
Electric Power Research Institute of State Grid Hubei 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 Wuhan University WHU, Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd filed Critical Wuhan University WHU
Priority to CN202210286544.3A priority Critical patent/CN114997457A/en
Publication of CN114997457A publication Critical patent/CN114997457A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313—Resource planning in a project environment
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E20/00—Combustion technologies with mitigation potential
    • Y02E20/14—Combined heat and power generation [CHP]

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Power Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Heat-Pump Type And Storage Water Heaters (AREA)

Abstract

The invention relates to a flexibility optimal scheduling method of a cogeneration unit based on heat storage characteristics of a central heating system. Step 1, collecting the temperature of each node of a pipeline of a heat supply network, the water supply temperature and the water return temperature of the pipeline, collecting the water flow of the pipeline, and obtaining the upper and lower output limits of a CHP unit based on a thermal system scheduling objective function and constraint conditions of the CHP unit; and 2, taking the minimum wind curtailment amount of the wind power as an optimization target, comprehensively considering the power balance, the power system constraint and the CHP unit multi-time scale flexible output constraint, optimizing a target function based on the total cost, and solving by using a commercial solver to obtain the specific output value of the CHP unit. Therefore, the invention can fully excavate the adjustment flexibility of the CHP unit, expand the wind power online space, increase the utilization rate of renewable energy, realize the complementary mutual benefit of electric heating comprehensive energy and realize the optimal economy.

Description

Flexibility optimization scheduling method for cogeneration unit based on heat storage characteristics of centralized heating system
Technical Field
The invention relates to a flexibility optimization scheduling method for a cogeneration unit, in particular to a flexibility optimization scheduling method for a cogeneration unit based on the heat storage characteristic of a centralized heating system.
Background
The principle of Combined Heat and Power (CHP) is that it can generate electric energy and heat energy at the same time, and it is a reliable and effective way to generate electric energy and heat energy by a single fuel source, and can raise the operation efficiency of the machine set and reduce the cost.
The cogeneration units are mainly divided into an extraction type and a back pressure type, and most of the thermoelectric units applied in China are large extraction type units nowadays.
The back pressure type thermoelectric unit supplies heat through waste heat discharged by back pressure, the thermoelectric comprehensive efficiency is high and can reach 80% -85%, but the thermoelectric ratio is fixed, the electric output is completely determined by the thermal output, the regulation capability is not available, and the installed capacity ratio in a power system is small.
Compared with a backpressure unit, the steam extraction unit has stronger flexibility and regulating capacity, can regulate the electric output and the thermal output in a larger range, but has lower thermoelectric comprehensive efficiency, and is generally suitable for thermal power plants with more frequent load change and larger change amplitude.
However, the larger the extraction steam quantity is, the smaller the proportion of steam for the condensed steam power generation is, so that the regulating range of the unit is smaller. A large amount of heat loads need to be met in winter, so that the condition necessarily leads to the situation that the minimum output of the unit is increased because the unit is constrained to be fixed by heat during the off-peak period of the electric load at night. Therefore, the peak regulation capability of the system is poor, a large amount of wind power cannot be connected to the grid, and the phenomenon of serious wind abandon is caused.
Disclosure of Invention
The technical problem of the invention is mainly solved by the following technical scheme:
a flexibility optimal scheduling method for a cogeneration unit based on the heat storage characteristic of a centralized heating system is characterized by comprising the following steps:
step 1, acquiring the temperature of each node of a pipeline of a heat supply network, the water supply temperature and the water return temperature of the pipeline, acquiring the water flow of the pipeline, taking the output of a conventional unit and a CHP unit and the standby of each unit as decision variables by a scheduling model on the basis of considering the multi-time scale flexibility constraint of the CHP unit, and obtaining the output upper and lower limits of the CHP unit based on a thermodynamic system scheduling objective function and the constraint conditions of the CHP unit;
step 2, taking the minimum wind curtailment amount of wind power as an optimization target in terms of total operation and standby cost, comprehensively considering power balance and power system constraints and CHP unit multi-time scale flexible output constraints, and based on a total cost optimization target function, solving the problem by using a commercial solver; finally obtaining the flexible scheduling of the CHP unit; namely, the specific output value of the CHP unit is obtained.
In the foregoing method for optimizing flexibility of a cogeneration unit based on heat storage characteristics of a central heating system, in step 1, a scheduling objective function of a thermodynamic system of a CHP unit is based on the following formula:
Figure BDA0003558476470000021
wherein T is the total time segment number,
Figure BDA0003558476470000022
are respectively provided withThe cost of the thermal output fuel of the cogeneration unit and other heat generating devices,
Figure BDA0003558476470000023
in order to store heat and lose cost,
Figure BDA0003558476470000024
and
Figure BDA0003558476470000025
and the penalty cost is respectively the upper bound deviation and the lower bound deviation of the CHP unit thermal output.
In the above method for optimizing and scheduling flexibility of a cogeneration unit based on the heat storage characteristics of a central heating system, in step 1,
Figure BDA0003558476470000031
Figure BDA0003558476470000032
Figure BDA0003558476470000033
μ b and mu c The fuel costs of the heating plant and the CHP unit respectively,
Figure BDA0003558476470000034
and
Figure BDA0003558476470000035
fuel consumption rate, p, of heating equipment and CHP units, respectively s The heat loss coefficient of the heat storage device in storage;
Figure BDA0003558476470000036
Figure BDA0003558476470000037
wherein, mu + And mu - As a deviation penalty function for the CHP unit,
Figure BDA0003558476470000038
is an intermediate variable of the construct.
In the foregoing method for optimizing flexibility of a cogeneration unit based on a heat storage characteristic of a central heating system, in step 1, the constraint condition includes:
temperature restraint:
Figure BDA0003558476470000039
Figure BDA00035584764700000310
Figure BDA00035584764700000311
Figure BDA00035584764700000312
Figure BDA00035584764700000313
Figure BDA00035584764700000314
Figure BDA00035584764700000315
Figure BDA00035584764700000316
force restraint:
Figure BDA00035584764700000317
Figure BDA00035584764700000318
Figure BDA00035584764700000319
Figure BDA00035584764700000320
Figure BDA0003558476470000041
Figure BDA0003558476470000042
Figure BDA0003558476470000043
S s,T ＝S s,0
wherein S is s,t Is the amount of heat stored at time t, S s,0 The heat storage amount at time 0;
Figure BDA0003558476470000044
and
Figure BDA0003558476470000045
the decision variables composed of the upper bound and the lower bound of the operation interval of the heat output of the CHP unit,
Figure BDA0003558476470000046
indicating the temperature at time t of the water supply node,
Figure BDA0003558476470000047
the temperature at time t of the water return node is shown,
Figure BDA0003558476470000048
which represents the inlet temperature of the pipe,
Figure BDA0003558476470000049
which represents the outlet temperature of the pipe or pipes,
Figure BDA00035584764700000410
representing the actual temperature of the water supply node at actual operation,
Figure BDA00035584764700000411
represents the inlet temperature of the pipe at actual operation,
Figure BDA00035584764700000412
represents the outlet temperature of the pipe at actual operation,
Figure BDA00035584764700000413
the temperature change value for measuring the heat storage quantity of the heat supply network pipeline is represented; all the subscripts min and max on the belt represent the minimum value and the maximum value of the corresponding parameters, and can be set manually;
Figure BDA00035584764700000414
the thermal output of the CHP plant is indicated,
Figure BDA00035584764700000415
and
Figure BDA00035584764700000416
respectively representing the upper and lower bounds of the preceding paragraph;
Figure BDA00035584764700000417
is an intermediate variable;
Figure BDA00035584764700000418
indicating the heat supply of the radiator, S s,t The heat storage amount at the moment t;
Figure BDA00035584764700000419
the heat output of other heating equipment set for the system.
In the foregoing method for optimizing and scheduling flexibility of a cogeneration unit based on the heat storage characteristic of a central heating system, in step 2, the total cost optimization objective function is based on the following formula:
Figure BDA00035584764700000420
wherein C is the total running cost of system coordination, C CON (t) and C CHP (t) operating costs of the thermal power generating unit and the thermoelectric power generating unit at the time period t, respectively, C CUT (t) wind curtailment penalty cost for t period, C s And (t) the heat dissipation cost of the radiator in the period of t.
In the flexibility optimization scheduling method of the cogeneration unit based on the heat storage characteristic of the central heating system,
the operating cost of the thermal power generating unit is as follows:
C CON (t)＝μ coal (aP con (t) 2 +bP con (t)+c)
wherein, mu coal Is coal value, P con (t) the power output of the thermal power generating unit in a period of t, and a, b and c are fitting coefficients of the operating cost of the thermal power generating unit;
C CHP (t)＝μ coal (f 0 D2+f 1 P CHP D+f 2 P C 2 HP +f 3 D+f 4 P CHP +f 5 )
f 0 to f 5 Is the fitting coefficient
P CHP Is the CHP unit electrical power; dThe air extraction amount is;
C CUT (t)＝[P WT,y(t) -P WT,s(t) ]·ε
wherein, P WT,y(t) And P WT,s(t) Respectively a predicted value and an actual output power of the wind power, wherein epsilon is a wind abandon punishment coefficient;
C s (t)＝c·Q s,t
wherein c is a heat dissipation cost coefficient; q s,t Heat released for the heat sink.
In the above method for optimizing and scheduling flexibility of a cogeneration unit based on the heat storage characteristics of a central heating system, the constraint conditions include:
and (3) power constraint:
P CHP (t)+P CON (t)+P WT,S (t)＝P load (t)
heat restraint:
Q CHP (t)+Q HS (t)＝Q load (t)
heat dissipation and building temperature constraints:
Figure BDA0003558476470000051
at is the scheduling time interval and,
Figure BDA0003558476470000061
and
Figure BDA0003558476470000062
indoor and outdoor temperatures at time t, k 1 、k 2 、k 3 Are the corresponding coefficients; c is the specific heat capacity; s is the heat supply area, is the indoor heat loss coefficient, Q s Indicating the heat supply amount of the radiator;
pipeline thermal delay and temperature drop constraints:
Figure BDA0003558476470000063
Figure BDA0003558476470000064
wherein, Δ T t For the temperature loss at the time t,
Figure BDA0003558476470000065
the temperature of the water supplied to the head end of the pipeline at the time t,
Figure BDA0003558476470000066
is the external temperature of the pipeline at time t, k loss Lambda is the heat transfer efficiency per unit length of the pipe;
the temperature at the end of the pipe at time t can be found to be:
Figure BDA0003558476470000067
unit output restraint:
Figure BDA0003558476470000068
temperature restraint:
Figure BDA0003558476470000069
Figure BDA00035584764700000610
Figure BDA00035584764700000611
Figure BDA00035584764700000612
wherein, P CHP Is represented by CElectric power of HP unit, P CON Representing the electrical output, P, of the thermal power generating unit WT Representing the electric output of the wind turbine; other bands with min sum and max are the upper and lower limits of corresponding parameters, and the sizes can be set manually by an actual model;
Figure BDA0003558476470000071
indicating the temperature at time t of the water supply node,
Figure BDA0003558476470000072
the temperature at time t of the water return node is shown,
Figure BDA0003558476470000073
which represents the inlet temperature of the pipe,
Figure BDA0003558476470000074
indicating the outlet temperature of the pipe.
Therefore, the invention has the following advantages: 1. the adjustment flexibility of the CHP unit can be fully excavated, the wind power internet space is enlarged, and the utilization rate of renewable energy sources is increased. 2. The complementary mutual benefit of the electric heating comprehensive energy can be realized, and the optimal economical efficiency is realized.
Drawings
FIG. 1 shows the working interval of the CHP unit of the present invention.
FIG. 2 is a flow chart of a method of the present invention.
Detailed Description
The technical scheme of the invention is further specifically described by the following embodiments and the accompanying drawings.
Example (b):
in the embodiment of the present invention, a set of inequalities may be first expressed for the CHP operation:
P gt ＞rQ gt
Figure BDA0003558476470000075
0≤Q gt ≤Q gmax
where g denotes a CHP unit, t denotes a scheduling time, P gt For the electrical output of the unit, Q gt The thermal output is r is the elastic coefficient of the electrical output and the thermal output under the working condition of back pressure, Q gmax Is the maximum value of the thermal output, beta e And beta h The specific fuel consumption of the electrical output and the thermal output are respectively expressed,Fand
Figure BDA0003558476470000076
representing the maximum and minimum fuel consumption respectively. Namely: r is the slope of the line segment BC, and the slopes of the line segments AB and CD are beta h /β e 。
In the conventional mode of using heat to fix power, when CHP is operated at Q 0 While the output power range of the electric power is floated at E 1 E 2 . At the moment, the electric output of the unit is P 0 Greater than the minimum electrical output P of the unit min If the proportion of the CHP unit is large, the capacity of the power system for wind power consumption is seriously influenced, and the flexibility of system operation is also influenced.
For the area with a large wind power grid-connected proportion, the traditional cogeneration unit is limited by the operation constraint of 'fixing the power with the heat' and has poor regulating capability, so that the wind power is insufficient in the internet space, and a large amount of abandoned wind is generated. In order to increase the adjusting capacity of the cogeneration unit and increase the wind power online space, the constraint of 'fixing the power with heat' of the cogeneration unit needs to be decoupled, and the configuration of heat storage for the cogeneration unit becomes an effective measure for solving the problem.
For example, after the heat storage and heat accumulation device is added, heat supply of the CHP can be reduced, and heat is supplemented by the heat accumulation device, so that the electric output range of the CHP can be enlarged, the minimum electric output is reduced, and the wind power consumption capability can be improved. The CHP is operated at Q 1 Electric power output range is E 3 E 6 Compared with the traditional mode, the unit obtains a part of space E for dispatching electric power 5 E 6 And the method can be used for wind power consumption, and the part can be called as scheduling flexibility of CHP. When the unit actually operates, the regulating output space E in a spare form can be obtained by coordinating with the heat supply network system 9 E 10 At Q of 2 。
In the dispatching stage, the dispatching heat balance of the CHP unit at the time t comprises a heat source exchange node and a heat load exchange node:
Figure BDA0003558476470000081
Figure BDA0003558476470000082
Figure BDA0003558476470000083
Figure BDA0003558476470000084
Figure BDA0003558476470000085
wherein the content of the first and second substances,
Figure BDA0003558476470000086
is the thermal output of the CHP unit,
Figure BDA0003558476470000087
the heat output of other heat supply equipment is provided,
Figure BDA0003558476470000088
is the thermal load of the i-node,
Figure BDA0003558476470000091
and
Figure BDA0003558476470000092
respectively the inlet temperature and the outlet temperature of the heat load exchange node,
Figure BDA0003558476470000093
and
Figure BDA0003558476470000094
respectively the water supply temperature and the water return temperature of the heat source exchange node. Lambda [ alpha ] ij Temperature transfer factor, η, for the water supply node j to the heat load inflow node i ij And the transfer factor of the heat exchange node i to the return water temperature j of the heat source is obtained.
Figure BDA0003558476470000095
Indicating the number of heat source exchange stations connected to other thermal power equipment,
Figure BDA0003558476470000096
the number of heat source exchange stations connected to the CHP plant is shown. Omega hes Representing a set of heat load nodes.
The heat supply network system has certain heat storage characteristics. The heat storage characteristic of the heat supply network system means that the heat source output fluctuation can be compensated by adjusting the overall temperature change of the heat supply network system. Under the quality adjusting mode, the pipeline node temperature of the heat supply network system can be divided into two types, one type is the node temperature of a water supply pipeline connected with a heat source node, the type of temperature is controllable under the influence of a heat source, the other type is the temperature of other water supply and return water nodes, and the type of temperature can be changed in a certain temperature interval. The matching between the heat supply network storage characteristics and the CHP units relates to the relationship between the heat supply network temperatures, and has the following three characteristics:
(1) the heat output is in direct proportion to the difference between the water delivery temperature and the return water temperature;
the change of the thermal output only affects the water delivery temperature;
(2) the water supply temperature is also influenced by the water return temperature.
Therefore, when the CHP unit participates in real-time regulation and control, the actual heat output is assumed to be
Figure BDA0003558476470000097
Regardless of the effect of the return water temperature, the heat balance at time t can be expressed as:
Figure BDA0003558476470000098
wherein the content of the first and second substances,
Figure BDA0003558476470000099
the supply water temperature is influenced by the output of the CHP unit. The heat imbalance of the thermodynamic system is caused by the change of the thermal output, and the unbalanced heat is assumed by the temperature change of the whole heat supply network, and the temperature change reflected to each heat source node can be represented as:
Figure BDA00035584764700000910
Figure BDA00035584764700000911
wherein the content of the first and second substances,
Figure BDA00035584764700000912
is the overall temperature change of the heat source exchange node,
Figure BDA00035584764700000913
and
Figure BDA00035584764700000914
the actual water supply temperature and the actual water return temperature are respectively, so the world heat balance equation at the time t becomes:
Figure BDA0003558476470000101
thereby can obtain
Figure BDA0003558476470000102
Figure BDA0003558476470000103
Figure BDA0003558476470000104
Figure BDA0003558476470000105
Wherein the content of the first and second substances,
Figure BDA0003558476470000106
and
Figure BDA0003558476470000107
the actual inlet and outlet temperatures of the heat load exchange node, respectively.
The above formula shows the influence of the heat output variation of the CHP unit on the node temperature of the heat supply network.
Figure BDA0003558476470000108
The heat storage characteristics of the heat supply network represent the coordination of the random change of the CHP unit output, and the heat storage level of the heat supply network for flexibly releasing the CHP unit can be quantitatively supported. At the same time, by controlling
Figure BDA0003558476470000109
The flexibility of the CHP unit and the heat storage and release capacity of the heat supply network can be ensured not to be excessively utilized. It can be seen that the temperature of the heat source node of other non-CHP equipment is also influenced by the output change of the CHP unit, so that the temperature of each node is controlled
Figure BDA00035584764700001010
And the influence of the CHP unit participating in power system scheduling and control on other heat source equipment can be reduced.
Because the allowable heat output interval of the CHP unit is determined by the heat storage characteristic, the scheduling flexibility and the adjustment flexibility release of the CHP unit both depend on the cooperation with the heat storage characteristic of a heat supply network, but the heat storage characteristic of the heat supply network cannot be quantized, the maximum allowable heat output operation interval of the CHP unit is difficult to determine, and the scheduling and the allocation of the multi-time scale flexibility of the CHP unit are directly influenced.
The scheduling targets of the thermodynamic system are:
Figure BDA00035584764700001011
wherein T is the total time segment number,
Figure BDA00035584764700001012
respectively the heat output fuel cost of the cogeneration unit and other heat generating devices,
Figure BDA0003558476470000111
in order to store heat and lose cost,
Figure BDA0003558476470000112
and
Figure BDA0003558476470000113
and the penalty cost is respectively the upper bound deviation and the lower bound deviation of the CHP unit thermal output.
Wherein the content of the first and second substances,
Figure BDA0003558476470000114
Figure BDA0003558476470000115
Figure BDA0003558476470000116
μ b and mu c The fuel costs of the heating plant and the CHP plant respectively,
Figure BDA0003558476470000117
and
Figure BDA0003558476470000118
fuel consumption rate, p, of heating equipment and CHP units, respectively s The heat loss coefficient of the heat storage device in storage.
Figure BDA0003558476470000119
Figure BDA00035584764700001110
Wherein, mu + And mu - As a bias penalty function for the CHP unit,
Figure BDA00035584764700001111
is an intermediate variable of the construct.
Temperature restraint:
Figure BDA00035584764700001112
Figure BDA00035584764700001113
Figure BDA00035584764700001114
Figure BDA00035584764700001115
Figure BDA00035584764700001116
Figure BDA00035584764700001117
Figure BDA00035584764700001118
Figure BDA00035584764700001119
force restraint:
Figure BDA0003558476470000121
Figure BDA0003558476470000122
Figure BDA0003558476470000123
Figure BDA0003558476470000124
Figure BDA0003558476470000125
Figure BDA0003558476470000126
Figure BDA0003558476470000127
S s,T ＝S s,0
wherein S is s,t Is the amount of heat stored at time t, S s,0 The amount of heat stored at time 0.
Figure BDA0003558476470000128
And
Figure BDA0003558476470000129
and the decision variables are formed by the upper bound and the lower bound of the operation interval of the heat output of the CHP unit.
By solving the above problem, the upper and lower output limits of the CHP unit can be obtained.
On the basis of considering the CHP unit multi-time scale flexibility constraint, the dispatching model takes the output of a conventional unit and the CHP unit and the standby of each unit as decision variables, takes the total operation and standby cost and the minimum wind curtailment amount of wind power as optimization targets, and comprehensively considers the power balance, the power system constraint and the CHP unit multi-time scale flexibility output constraint.
The electric heat coordination optimization under multiple thermal inertia is to fully utilize the characteristics of multiple inertias of the heat supply system, reasonably adjust the output of each unit, realize the economic operation of the combined heat and power system, and the optimization target can be expressed as:
Figure BDA00035584764700001210
wherein C is the total running cost of system coordination, C CON (t) and C CHP (t) operating costs of the thermal power generating unit and the thermoelectric power generating unit at a time t, respectively, C CUT (t) wind curtailment penalty cost for t period, C s And (t) is the heat dissipation cost of the radiator in the period of t.
The operating cost of the thermal power generating unit is as follows:
C CON (t)＝μ coal (aP con (t) 2 +bP con (t)+c)
wherein, mu coal Is the price of coal, P con And (t) is the electric output of the thermal power generating unit in the time period t, and a, b and c are the operation cost fitting coefficients of the thermal power generating unit.
Figure BDA0003558476470000131
f 0 To f 5 Is the fitting coefficient
P CHP Is the CHP plant electrical power. D is the air extraction amount.
C CUT (t)＝[P WT,y(t) -P WT,s(t) ]·ε
Wherein, P WT,y(t) And P WT,s(t) The predicted value and the actual output power of the wind power are respectively, and epsilon is a wind abandon punishment coefficient.
C s (t)＝c·Q s,t
Wherein c is a heat dissipation cost coefficient. Q s,t Heat released for the heat sink.
And (3) power constraint:
P CHP (t)+P CON (t)+P WT,S (t)＝P load (t)
heat restraint:
Q CHP (t)+Q HS (t)＝Q load (t)
heat dissipation and building temperature constraints:
Figure BDA0003558476470000132
at is the scheduling time interval and,
Figure BDA0003558476470000133
and
Figure BDA0003558476470000134
is the indoor and outdoor temperature at time t, k 1 、k 2 、k 3 Are the corresponding coefficients. C is specific heat capacity. And s is the heat supply area and the indoor heat loss coefficient.
Pipeline thermal delay and temperature drop constraints:
Figure BDA0003558476470000141
Figure BDA0003558476470000142
wherein, Delta T t Is the temperature loss at the time of the t,
Figure BDA0003558476470000143
the temperature of the water supplied to the head end of the pipeline at the time t,
Figure BDA0003558476470000144
is the external temperature of the pipeline at time t, k loss λ is the heat transfer efficiency per unit length of the pipe, which is the temperature loss coefficient.
The temperature at the end of the pipe at time t can be found to be:
Figure BDA0003558476470000145
unit output restraint:
Figure BDA0003558476470000146
temperature restraint:
Figure BDA0003558476470000147
Figure BDA0003558476470000148
Figure BDA0003558476470000149
Figure BDA00035584764700001410
and finally, converting the problem into a quadratic programming problem, and solving the quadratic programming problem by utilizing a commercial solver. And finally, obtaining the flexible scheduling of the CHP unit. Namely, the specific output value of the CHP unit is obtained.
The specific embodiments described herein are merely illustrative of the spirit of the invention. Various modifications or additions may be made to the described embodiments, or alternatives may be employed, by those skilled in the art, without departing from the spirit or ambit of the invention as defined in the appended claims.

Claims (7)

1. A flexibility optimization scheduling method for a cogeneration unit based on heat storage characteristics of a central heating system is characterized by comprising the following steps:
step 1, acquiring the temperature of each node of a pipeline of a heat supply network, the water supply temperature and the water return temperature of the pipeline, acquiring the water flow of the pipeline, taking the output of a conventional unit and a CHP unit and the standby of each unit as decision variables by a scheduling model on the basis of considering the multi-time scale flexibility constraint of the CHP unit, and obtaining the output upper and lower limits of the CHP unit based on a thermodynamic system scheduling objective function and the constraint conditions of the CHP unit;
step 2, taking the minimum wind curtailment amount of wind power as an optimization target in terms of total operation and standby cost, comprehensively considering power balance and power system constraints and CHP unit multi-time scale flexible output constraints, and based on a total cost optimization target function, solving the problem by using a commercial solver; finally obtaining the flexible scheduling of the CHP unit; namely, the specific output value of the CHP unit is obtained.
2. The flexibility optimization scheduling method of the cogeneration unit based on the heat storage characteristic of the central heating system according to claim 1, wherein in the step 1, the scheduling objective function of the thermodynamic system of the CHP unit is based on the following formula:
Figure FDA0003558476460000011
wherein T is the total time segment number,
Figure FDA0003558476460000012
respectively the heat output fuel cost of the cogeneration unit and other heat generating devices,
Figure FDA0003558476460000013
in order to store the heat and lose the cost,
Figure FDA0003558476460000014
and
Figure FDA0003558476460000015
and the penalty cost is respectively the upper bound deviation and the lower bound deviation of the CHP unit thermal output.
3. The flexibility optimal scheduling method for the cogeneration unit based on the heat storage characteristic of the central heating system according to claim 2, wherein in the step 1,
Figure FDA0003558476460000021
Figure FDA0003558476460000022
Figure FDA0003558476460000023
μ b and mu c The fuel costs of the heating plant and the CHP unit respectively,
Figure FDA0003558476460000024
and
Figure FDA0003558476460000025
fuel consumption rate, p, of heating equipment and CHP units, respectively s For heat storage apparatusHeat loss coefficients in storage;
Figure FDA0003558476460000026
Figure FDA0003558476460000027
wherein, mu + And mu - As a bias penalty function for the CHP unit,
Figure FDA0003558476460000028
is an intermediate variable of the construct.
4. The flexibility optimization scheduling method for the cogeneration unit based on the heat storage characteristic of the central heating system according to claim 3, wherein in the step 1, the constraint condition comprises:
temperature restraint:
Figure FDA0003558476460000029
Figure FDA00035584764600000210
Figure FDA00035584764600000211
Figure FDA00035584764600000212
Figure FDA00035584764600000213
Figure FDA00035584764600000214
Figure FDA00035584764600000215
Figure FDA00035584764600000216
force restraint:
Figure FDA00035584764600000217
Figure FDA00035584764600000218
Figure FDA00035584764600000219
Figure FDA00035584764600000220
Figure FDA0003558476460000031
Figure FDA0003558476460000032
Figure FDA0003558476460000033
S s,T ＝S s,0
wherein S is s,t Is the amount of heat stored at time t, S s,0 The heat storage amount at time 0;
Figure FDA0003558476460000034
and
Figure FDA0003558476460000035
the decision variables composed of the upper bound and the lower bound of the operation interval of the heat output of the CHP unit,
Figure FDA0003558476460000036
indicating the temperature at time t of the water supply node,
Figure FDA0003558476460000037
the temperature at time t of the water return node is shown,
Figure FDA0003558476460000038
which represents the inlet temperature of the pipe,
Figure FDA0003558476460000039
which represents the outlet temperature of the pipe or pipes,
Figure FDA00035584764600000310
representing the actual temperature of the water supply node at actual operation,
Figure FDA00035584764600000311
represents the inlet temperature of the pipe at actual operation,
Figure FDA00035584764600000312
represents the outlet temperature of the pipe at actual operation,
Figure FDA00035584764600000313
the temperature change value for measuring the heat storage quantity of the heat supply network pipeline is represented; all the upper subscripts min and max represent the minimum and maximum values of the corresponding parameters, and can be set manually;
Figure FDA00035584764600000314
the thermal output of the CHP plant is indicated,
Figure FDA00035584764600000315
and
Figure FDA00035584764600000316
respectively representing the upper and lower bounds of the preceding paragraph;
Figure FDA00035584764600000317
is an intermediate variable;
Figure FDA00035584764600000318
indicating the heat supply of the radiator, S s,t The heat storage amount at the moment t;
Figure FDA00035584764600000319
the heat output of other heating equipment set for the system.
5. The flexibility optimization scheduling method for the cogeneration unit based on the heat storage characteristic of the central heating system according to claim 1, wherein in the step 2, the total cost optimization objective function is based on the following formula:
Figure FDA00035584764600000320
wherein C is the total running cost of system coordination, C CON (t) and C CHP (t) operating costs of the thermal power generating unit and the thermoelectric power generating unit at a time t, respectively, C CUT (t) isPenalty cost of wind curtailment at t time period, C s And (t) the heat dissipation cost of the radiator in the period of t.
6. The flexibility optimization scheduling method of cogeneration units based on heat storage characteristics of central heating system as claimed in claim 5,
the operating cost of the thermal power generating unit is as follows:
C CON (t)＝μ coal (aP con (t) 2 +bP con (t)+c)
wherein, mu coal Is the price of coal, P con (t) the power output of the thermal power generating unit in a period of t, and a, b and c are the operation cost fitting coefficients of the thermal power generating unit;
Figure FDA0003558476460000041
f 0 to f 5 Is the fitting coefficient
P CHP Is the CHP unit electrical power; d is the air extraction amount;
C CUT (t)＝[P WT,y(t) -P WT,s(t) ]·ε
wherein, P WT,y(t) And P WT,s(t) Respectively a predicted value and an actual output power of the wind power, wherein epsilon is a wind abandon punishment coefficient;
C s (t)＝c·Q s,t
wherein c is a heat dissipation cost coefficient; q s,t Heat released for the heat sink.
7. The flexibility optimal scheduling method for the cogeneration unit based on the heat storage characteristic of the central heating system according to claim 5, wherein the constraint condition comprises:
and (3) power constraint:
P CHP (t)+P CON (t)+P WT,S (t)＝P load (t)
heat restraint:
Q CHP (t)+Q HS (t)＝Q load (t)
heat dissipation and building temperature constraints:
Figure FDA0003558476460000042
at is the scheduling time interval and,
Figure FDA0003558476460000051
and
Figure FDA0003558476460000052
is the indoor and outdoor temperature at time t, k 1 、k 2 、k 3 Are the corresponding coefficients; c is the specific heat capacity; s is the heat supply area, is the indoor heat loss coefficient, Q s Indicating the heat supply amount of the radiator;
pipeline thermal delay and temperature drop constraints:
Figure FDA0003558476460000053
Figure FDA0003558476460000054
wherein, Delta T t For the temperature loss at the time t,
Figure FDA0003558476460000055
the temperature of the water supplied to the head end of the pipeline at the time t,
Figure FDA0003558476460000056
is the external temperature of the pipeline at time t, k loss Lambda is the heat transfer efficiency per unit length of the pipe;
the temperature at the end of the pipe at time t can be found to be:
Figure FDA0003558476460000057
unit output constraint:
Figure FDA0003558476460000058
temperature restraint:
Figure FDA0003558476460000059
Figure FDA00035584764600000510
Figure FDA00035584764600000511
Figure FDA00035584764600000512
wherein, P CHP Indicating the electrical output of the CHP unit, P CON Representing the electrical output, P, of the thermal power generating unit WT Representing the electric output of the wind turbine; other bands with min sum and max are the upper and lower limits of corresponding parameters, and the sizes can be set manually by an actual model;
Figure FDA0003558476460000061
indicating the temperature at time t of the water supply node,
Figure FDA0003558476460000062
the temperature at time t of the water return node is shown,
Figure FDA0003558476460000063
which represents the inlet temperature of the pipe,
Figure FDA0003558476460000064
indicating the outlet temperature of the pipe.
CN202210286544.3A 2022-03-22 2022-03-22 Flexibility optimization scheduling method for cogeneration unit based on heat storage characteristic of centralized heating system Pending CN114997457A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210286544.3A CN114997457A (en) 2022-03-22 2022-03-22 Flexibility optimization scheduling method for cogeneration unit based on heat storage characteristic of centralized heating system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210286544.3A CN114997457A (en) 2022-03-22 2022-03-22 Flexibility optimization scheduling method for cogeneration unit based on heat storage characteristic of centralized heating system

Publications (1)

Publication Number Publication Date
CN114997457A true CN114997457A (en) 2022-09-02

Family

ID=83024245

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210286544.3A Pending CN114997457A (en) 2022-03-22 2022-03-22 Flexibility optimization scheduling method for cogeneration unit based on heat storage characteristic of centralized heating system

Country Status (1)

Country Link
CN (1) CN114997457A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116680935A (en) * 2023-07-31 2023-09-01 天津六百光年智能科技有限公司 Two-network water supply prediction model construction method based on piecewise linear function

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116680935A (en) * 2023-07-31 2023-09-01 天津六百光年智能科技有限公司 Two-network water supply prediction model construction method based on piecewise linear function
CN116680935B (en) * 2023-07-31 2023-10-13 天津六百光年智能科技有限公司 Two-network water supply prediction model construction method based on piecewise linear function

Similar Documents

Publication Publication Date Title
CN109767105B (en) Wind, light, water, fire and storage combined system-based multi-energy complementary coordinated power generation scheduling method
CN107039990B (en) A kind of multi-heat source coordinated scheduling policy development method improving wind electricity digestion capability
CN109494784A (en) A kind of large capacity heat reservoir optimal control method improving wind electricity digestion
CN107332286B (en) A kind of cogeneration of heat and power containing heat accumulation and wind-powered electricity generation coordinated scheduling method
CN106593771B (en) Electric heating decoupled system and its control method in a kind of co-generation unit
CN110676847A (en) Optimal scheduling method considering wind power-heat storage unit-electric boiler combined operation
CN108009683A (en) A kind of high-capacity direct control load improves the dispatching method that wind-powered electricity generation receives ability
CN110458353A (en) Electric heating combines microgrid energy step optimization method and system
CN104179575B (en) Combined cycle gas-steam turbine unit is Optimal Scheduling and method a few days ago
CN102563959A (en) Integrated energy matching system and control method thereof
CN108628176A (en) It is a kind of meter and heat supply network combined heat and power Optimization Scheduling flexible
CN110165665A (en) A kind of source-lotus-storage dispatching method based on improvement multi-objective particle swarm algorithm
CN114154328A (en) Flexibility-improved two-stage random optimization scheduling method for electric heating comprehensive energy system
CN113190941A (en) Centralized heating pipe network modeling method suitable for electric-heat joint scheduling
CN112434915A (en) Regional comprehensive energy system flexibility optimization configuration method for curtailment consumption
CN111429020A (en) Multi-time scale economic scheduling method of electric heating system considering heat storage characteristics of heat supply network
CN114997457A (en) Flexibility optimization scheduling method for cogeneration unit based on heat storage characteristic of centralized heating system
CN112671040A (en) Day-ahead optimal scheduling method of multi-energy complementary system considering maximum new energy consumption
CN111428962A (en) Electric heating cooperative scheduling method based on heat supply network dynamic characteristic port equivalent model
CN111146819B (en) Electric heating combined system auxiliary service scheduling method considering heat supply network characteristics
Wu et al. Day-ahead optimal dispatch with CHP and wind turbines based on room temperature control
CN108960564A (en) Supply of cooling, heating and electrical powers type microgrid distribution system dispatching method based on centralization modeling
CN108494014A (en) A kind of energy mix cogeneration of heat and power economy optimum management method
CN117035202A (en) Double-layer collaborative expansion planning method for electric heating comprehensive energy system considering demand response
CN113762643B (en) Energy storage capacity optimization configuration method of regional comprehensive energy 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