CN111563634A - Day-ahead optimization scheduling method for comprehensive energy system considering heat supply network reconstruction - Google Patents

Day-ahead optimization scheduling method for comprehensive energy system considering heat supply network reconstruction Download PDF

Info

Publication number
CN111563634A
CN111563634A CN202010447996.6A CN202010447996A CN111563634A CN 111563634 A CN111563634 A CN 111563634A CN 202010447996 A CN202010447996 A CN 202010447996A CN 111563634 A CN111563634 A CN 111563634A
Authority
CN
China
Prior art keywords
power
scheduling
thermal
scheduling time
energy system
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
CN202010447996.6A
Other languages
Chinese (zh)
Other versions
CN111563634B (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
Original Assignee
Tsinghua University
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 filed Critical Tsinghua University
Priority to CN202010447996.6A priority Critical patent/CN111563634B/en
Publication of CN111563634A publication Critical patent/CN111563634A/en
Application granted granted Critical
Publication of CN111563634B publication Critical patent/CN111563634B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas 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
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention relates to a day-ahead optimization scheduling method for a comprehensive energy system considering heat supply network reconstruction, and belongs to the technical field of operation control of the comprehensive energy system. The scheduling method of the invention establishes an objective function which minimizes the running cost; establishing constraint conditions including the contents of a multi-class equipment output calculation equation, range constraint, heat supply network constraint, power grid constraint, power balance and the like in consideration of reconstruction; by solving the optimized scheduling model, a comprehensive energy system day-ahead optimized scheduling scheme of statistics and heat supply network reconstruction can be obtained. According to the scheduling method, the heat supply network reconstruction is used as an adjusting means in day-ahead optimized scheduling, the flexibility of heat supply network adjustment can be fully exerted, the aims of more consumption of renewable energy sources, reduction of transmission blockage of a power grid and the like are achieved, and the efficient operation of a comprehensive energy system is promoted.

Description

Day-ahead optimization scheduling method for comprehensive energy system considering heat supply network reconstruction
Technical Field
The invention relates to a day-ahead optimization scheduling method for a comprehensive energy system considering heat supply network reconstruction, and belongs to the technical field of operation control of the comprehensive energy system.
Technical Field
The continuous development of social economy leads to the increasing of power load, and the problem of blocking in a power grid has attracted the attention of research and application level. Meanwhile, the installed scale of the renewable energy is continuously enlarged, but due to the influence of the time-space characteristics of the renewable energy, the power grid regulation capacity, the power grid transmission limitation and other factors, wind and light abandonment occurs, and the effective application of the renewable energy cannot be realized.
The cogeneration unit takes natural gas as input energy, simultaneously produces electric power and heat, can realize the cascade utilization of energy, and is a more efficient energy utilization and production mode. The combined heat and power generation unit is widely applied, so that the coupling of a power network, a heat supply/cold supply network and a gas network is realized, and a comprehensive energy system is formed. Energy flows in the comprehensive energy system are mutually interwoven and mutually influenced. By reasonably adjusting coupling equipment such as a cogeneration unit and the like, adjustable resources in the traditional power grid can be increased, and the problems of transmission blockage, wind abandonment, light abandonment and the like in the power grid are expected to be relieved or even eliminated. At present, an optimized scheduling method of a comprehensive energy system researched and applied is mainly used for reasonably arranging a start-stop plan and an output plan of a cogeneration unit and effectively utilizing thermal inertia of a heat supply network to relieve the problems in the operation of a power grid. All of the scheduling methods are based on a fixed heat supply network topological structure, and the application of the heat supply network regulation capacity is limited. And (3) considering the optimized scheduling of the heat supply network reconstruction, namely, taking the heat supply network topology as a regulating variable in the optimized scheduling, so that new adjustable capacity is provided for comprehensive energy, and the problem of the power grid in operation is further solved.
Disclosure of Invention
The invention aims to provide a day-ahead optimization scheduling method of a comprehensive energy system considering heat supply network reconstruction, aiming at increasing the regulation flexibility of the comprehensive energy system, increasing the consumption of renewable energy and relieving the transmission blockage of a power grid through the heat supply network reconstruction.
The invention provides a day-ahead optimization scheduling method for a comprehensive energy system considering heat supply network reconstruction, which comprises the following steps of:
(1) establishing an objective function of day-ahead optimized scheduling of the heat supply network reconstruction comprehensive energy system:
Figure BDA0002506639130000021
where T is the sequence number of the scheduling time, T is the set of scheduling times in one day, i1、i2、i3、i4Respectively numbering a cogeneration unit, a heat supply boiler, a traditional generator and a wind turbine generator in the comprehensive energy system, wherein n is the numbering of a power bus in the comprehensive energy system; kCHPAs a collection of cogeneration units, KHBFor integration of boilers, KTUAs a collection of conventional generators, KWDAs a collection of wind turbines, KbusIs a collection of power buses;
Figure BDA0002506639130000022
and
Figure BDA0002506639130000023
respectively a cogeneration unit i1Active power and thermal power generated at the scheduling time t;
Figure BDA0002506639130000024
to characterize cogeneration units i1Binary variable of running state at scheduling time t, cogeneration unit i1Is in the run state at the scheduled time t,
Figure BDA0002506639130000025
value of 1, cogeneration unit i1Is in a power-off state at a scheduled time t,
Figure BDA0002506639130000026
is 0;
Figure BDA0002506639130000027
to characterize cogeneration units i1Binary variable of starting-up command at scheduling time t, cogeneration unit i1The power-on command is executed at the scheduled time t,
Figure BDA0002506639130000028
value of 1, cogeneration unit i1The power-on command is not executed at the scheduling time t,
Figure BDA0002506639130000029
is 0;
Figure BDA00025066391300000210
to characterize cogeneration units i1Binary variable of shutdown instruction at scheduling time t, cogeneration unit i1The shutdown instruction is executed at the scheduled time t,
Figure BDA00025066391300000211
value of 1, cogeneration unit i1The shutdown instruction is not executed at the scheduled time t,
Figure BDA00025066391300000212
is 0;
Figure BDA00025066391300000213
for heating boilers i2Thermal power generated at the scheduling time t;
Figure BDA00025066391300000214
for a conventional generator i3Active power generated at a scheduling time t;
Figure BDA00025066391300000215
to characterize a conventional generator i3Binary variable of running state at scheduling time t, conventional generator i3Is in the run state at the scheduled time t,
Figure BDA00025066391300000216
value of 1, conventional generator i3Is in a power-off state at a scheduled time t,
Figure BDA00025066391300000217
is 0;
Figure BDA00025066391300000218
to characterize a conventional generator i3Binary variable of start-up command at scheduling time t, conventional generator i3The power-on command is executed at the scheduled time t,
Figure BDA00025066391300000219
value of 1, conventional generator i3The power-on command is not executed at the scheduling time t,
Figure BDA00025066391300000220
is 0;
Figure BDA00025066391300000221
to characterize a conventional generator i3Binary variable of shutdown command at scheduling time t, conventional generator i3The shutdown instruction is executed at the scheduled time t,
Figure BDA00025066391300000222
value of 1, conventional generator i3The shutdown instruction is not executed at the scheduled time t,
Figure BDA00025066391300000223
is 0;
Figure BDA00025066391300000224
for wind turbine generator i4The active power generated at the scheduling instant t,
Figure BDA00025066391300000225
active power for load shedding of the power bus n at the scheduling time t;
Figure BDA00025066391300000226
for cogeneration units i1The operating cost function of (a);
Figure BDA0002506639130000031
for heating boilers i2The operating cost function of (a);
Figure BDA0002506639130000032
for a conventional generator i3The operating cost function of (a);
Figure BDA0002506639130000033
for wind turbine generator i4A cost function of operation;
Figure BDA0002506639130000034
an operating cost function for shedding the power bus n;
(2) establishing a constraint condition of day-ahead optimized scheduling of the heat supply network reconstruction comprehensive energy system:
(2-1) establishing active power and thermal power constraints of a cogeneration unit in the comprehensive energy system:
Figure BDA0002506639130000035
Figure BDA0002506639130000036
Figure BDA0002506639130000037
wherein k is the serial number of the polygon end point of the feasible region of the cogeneration unit,
Figure BDA0002506639130000038
for cogeneration units i1The total number of end points of the feasible region polygon,
Figure BDA0002506639130000039
for cogeneration units i1The weight of the k-th endpoint in the feasible domain,
Figure BDA00025066391300000310
for cogeneration units i1The active power value of the k-th endpoint of the polygon in the feasible region,
Figure BDA00025066391300000311
for cogeneration units i1A heat value of a k-th endpoint of the feasible region polygon;
(2-2) establishing heat power constraint of a heat supply boiler in the comprehensive energy system:
Figure BDA00025066391300000312
Figure BDA00025066391300000313
in the formula (I), the compound is shown in the specification,
Figure BDA00025066391300000314
for heating boilers i2The thermal efficiency of the heat-generating element,
Figure BDA00025066391300000315
for entering the heating boiler i at the scheduled time t2The power of the natural gas (es) of (c),
Figure BDA00025066391300000316
for heating boilers i2Maximum value of thermal power;
(2-3) establishing a heat power balance equation and a range constraint in the comprehensive energy system:
Figure BDA00025066391300000317
Figure BDA00025066391300000318
Figure BDA00025066391300000319
Figure BDA00025066391300000320
Figure BDA00025066391300000321
Figure BDA0002506639130000041
wherein j is the number of the heating power station consisting of the cogeneration unit and the heating boiler in the step (1), and KHSAs a collection of thermal power stations, Kj CHPIs the set of all cogeneration units in the heating power station j, Kj HBBeing the collection of all the heating boilers in the heating station j,
Figure BDA0002506639130000042
the thermal power provided for the thermal station j at the scheduling time t, b is the number of the thermal pipeline, KPIs the collection of all the thermal power conduits,
Figure BDA0002506639130000043
is the thermal power at the outlet of the thermal pipeline b,
Figure BDA0002506639130000044
the thermal power of the inlet of the thermal pipeline b,
Figure BDA0002506639130000045
the thermal power lost in the thermal power pipeline b, h is the number of the heat supply network nodes, KNFor the set of all the heat network nodes,
Figure BDA0002506639130000046
for the set of thermal conduits with the heat network node h as the end point,
Figure BDA0002506639130000047
for a collection of thermal conduits with the heat network node h as the starting point,
Figure BDA0002506639130000048
for a collection of thermal stations connected to a heat supply network node h, i5The numbers of the heat exchange stations in the comprehensive energy system,
Figure BDA0002506639130000049
being a collection of heat exchange stations connected to the heat network node h,
Figure BDA00025066391300000410
to heat exchange station i5At the moment of scheduling the thermal power at the instant t,
Figure BDA00025066391300000411
for cogeneration units i1The lower limit of the thermal power is,
Figure BDA00025066391300000412
for cogeneration units i1Upper limit of thermal power, ub,tIs a binary variable for judging whether the thermal pipeline b is communicated at the scheduling time t, the thermal pipeline b is in a communicated state at the scheduling time t, and u isb,tIs 1, the thermal pipe b is in the closed state at the scheduling instant t, ub,tThe value of (a) is 0,
Figure BDA00025066391300000413
the upper limit of the thermal power is delivered to the thermal pipeline b;
(2-4) establishing range constraint of the heat supply network reconstruction state in the comprehensive energy system:
Figure BDA00025066391300000414
Figure BDA00025066391300000415
Figure BDA00025066391300000416
Figure BDA00025066391300000417
in the formula ub,t-1A binary variable u indicating whether the thermal pipeline b is communicated at the scheduling time t-1b,t-1Value-taking rules and step (2)U in (3)b,tThe value-taking rules are the same; x is the number ofb,tTo characterize the binary variable of the opening instruction of thermal pipe b at scheduling time t, thermal pipe b executes the opening instruction at scheduling time t, xb,tIs 1, the thermal pipe b does not execute an opening command at the scheduling time t, xb,tIs 0; y isb,tTo characterize the binary variable of the closing instruction of thermal pipe b at scheduling time t, thermal pipe b executes the closing instruction at scheduling time t, yb,tIs 1, the thermal pipe b does not execute a closing command at the scheduling instant t, yb,tIs 0; x is the number ofb,σBinary variable, x, characterizing the command for opening the thermal pipe b at the scheduling instant sigmab,σValue rule of (1) and xb,tThe value-taking rules are the same; y isb,σTo characterize the binary variable of the closing command of the thermodynamic pipe b at the scheduling instant sigma, yb,σValue rule of (a) and yb,tThe value-taking rules are the same; MRbFor a minimum operating time interval of the electric valve on the thermal pipe b,
Figure BDA0002506639130000051
the initial state of the thermal pipeline b at the scheduling time t;
(2-5) establishing an active power balance equation in the comprehensive energy system:
Figure BDA0002506639130000052
in the formula, Dn,tThe active power consumed by the electric load on the electric bus n at the scheduling time t;
(2-6) establishing range constraint of active power transmission in the comprehensive energy system:
Figure BDA0002506639130000053
wherein l is the number of the power line in the integrated energy system, KLFor the aggregate of all power lines, SFl,nFor power transfer distribution factor, SFl,nIs obtained from the comprehensive energy management system of the energy,
Figure BDA0002506639130000054
being a collection of conventional generators connected on the power bus n,
Figure BDA0002506639130000055
for the set of cogeneration units connected to the power bus n,
Figure BDA0002506639130000056
for a collection of wind turbines connected to a power bus n, FlAn upper limit of the active power transmitted for the power line l;
(2-7) establishing range constraint of active power in the comprehensive energy system:
Figure BDA0002506639130000057
Figure BDA0002506639130000058
Figure BDA0002506639130000059
in the formula (I), the compound is shown in the specification,
Figure BDA00025066391300000510
and
Figure BDA00025066391300000511
respectively a cogeneration unit i1The lower and upper limits of the active power generated,
Figure BDA00025066391300000512
and
Figure BDA00025066391300000513
respectively a conventional generator i3The lower and upper limits of the active power generated,
Figure BDA00025066391300000514
for wind turbine generator i4An upper limit of active power generated;
(2-8) establishing the range constraint of active power climbing in the comprehensive energy system:
Figure BDA00025066391300000515
Figure BDA0002506639130000061
Figure BDA0002506639130000062
Figure BDA0002506639130000063
in the formula (I), the compound is shown in the specification,
Figure BDA0002506639130000064
to characterize cogeneration units i1A binary variable of the running state at the scheduling time t +1,
Figure BDA0002506639130000065
value-taking rule of (1)
Figure BDA0002506639130000066
The value-taking rules are the same;
Figure BDA0002506639130000067
and
Figure BDA0002506639130000068
respectively a cogeneration unit i1The power-off rate and the power-on rate of,
Figure BDA0002506639130000069
and
Figure BDA00025066391300000610
respectively a cogeneration uniti1The downward climbing rate and the upward climbing rate;
Figure BDA00025066391300000611
for cogeneration units i1Active power generated at scheduling time t + 1;
Figure BDA00025066391300000612
to characterize a conventional generator i3A binary variable of the running state at the scheduling time t +1,
Figure BDA00025066391300000613
value-taking rule of (1)
Figure BDA00025066391300000614
The value-taking rules are the same;
Figure BDA00025066391300000615
and
Figure BDA00025066391300000616
respectively a conventional generator i3The power-off rate and the power-on rate of,
Figure BDA00025066391300000617
and
Figure BDA00025066391300000618
respectively a conventional generator i3Downward climbing rate and upward climbing rate of
Figure BDA00025066391300000619
For a conventional generator i3Active power generated at scheduling time t + 1;
(2-9) establishing range constraint of active power rotation standby in the comprehensive energy system:
Figure BDA00025066391300000620
Figure BDA00025066391300000621
in the formula (I), the compound is shown in the specification,
Figure BDA00025066391300000622
for each scheduling time
Figure BDA00025066391300000623
Maximum value of, SRUtReserve capacity, SRD, for an integrated energy system at scheduled time ttotaitThe total downward reserve capacity of the comprehensive energy system at the scheduling time t;
(2-10) establishing unit start-stop logic constraint in the comprehensive energy system:
Figure BDA00025066391300000624
Figure BDA00025066391300000625
in the formula (I), the compound is shown in the specification,to characterize cogeneration units i1A binary variable of the running state at the scheduling time t-1,
Figure BDA00025066391300000627
value-taking rule of (1)
Figure BDA00025066391300000628
The value-taking rules are the same;
Figure BDA00025066391300000629
to characterize cogeneration units i1The binary variable of the shutdown instruction at scheduling time t-1,
Figure BDA00025066391300000630
value-taking rule of (1)
Figure BDA00025066391300000631
The value-taking rules are the same;
Figure BDA00025066391300000632
to characterize a conventional generator i3A binary variable of the running state at the scheduling time t-1,
Figure BDA00025066391300000633
value-taking rule of (1)
Figure BDA00025066391300000634
The value-taking rules are the same;
Figure BDA0002506639130000071
to characterize a conventional generator i3The binary variable of the shutdown instruction at scheduling time t-1,
Figure BDA0002506639130000072
value-taking rule of (1)
Figure BDA0002506639130000073
The value-taking rules are the same;
(2-11) establishing a range constraint of the minimum startup time and the minimum shutdown time in the comprehensive energy system:
Figure BDA0002506639130000074
Figure BDA0002506639130000075
Figure BDA0002506639130000076
Figure BDA0002506639130000077
in the formula (I), the compound is shown in the specification,
Figure BDA0002506639130000078
and
Figure BDA0002506639130000079
respectively a cogeneration unit i1The minimum on-time and the minimum off-time,
Figure BDA00025066391300000710
and
Figure BDA00025066391300000711
respectively a conventional generator i3Minimum turn-on time and minimum turn-off time;
(3) the method comprises the steps of (1) forming a comprehensive energy system day-ahead optimization scheduling model considering heat supply network reconstruction by the objective function of the heat supply network reconstruction comprehensive energy system day-ahead optimization scheduling established in the step (1) and the constraint conditions of the heat supply network reconstruction comprehensive energy system day-ahead optimization scheduling established in the step (2), solving the optimization model by adopting a branch-and-bound method to obtain the comprehensive energy system day-ahead optimization scheduling parameters considering the heat supply network reconstruction, and including a cogeneration unit i in the step (1)1Active power generated at scheduling time t
Figure BDA00025066391300000712
And thermal power
Figure BDA00025066391300000713
Characterization of cogeneration units i1Binary variable of running state at scheduling time t
Figure BDA00025066391300000714
Characterization of cogeneration units i1Binary variable of starting instruction at scheduling time t
Figure BDA00025066391300000715
Characterization of cogeneration units i1Binary variable of shutdown instruction at scheduling time t
Figure BDA00025066391300000716
Heating boiler i2Thermal power generated at scheduled time t
Figure BDA00025066391300000717
Conventional generator i3Active power generated at scheduling time t
Figure BDA00025066391300000718
Characterizing a conventional Generator i3Binary variable of running state at scheduling time t
Figure BDA00025066391300000719
Characterizing a conventional Generator i3Binary variable of starting instruction at scheduling time t
Figure BDA00025066391300000720
Characterizing a conventional Generator i3Binary variable of shutdown instruction at scheduling time t
Figure BDA00025066391300000721
Binary variable u for judging whether thermal pipeline b is communicated at scheduling time t or notb,tBinary variable x for representing opening instruction of thermal pipeline b at scheduling time tb,tBinary variable y for representing closing instruction of thermal pipeline b at scheduling time tb,t。
The day-ahead optimization scheduling method of the comprehensive energy system considering heat supply network reconstruction, provided by the invention, has the advantages that:
the invention relates to a day-ahead optimization scheduling method of a comprehensive energy system considering heat supply network reconstruction, which establishes an objective function with minimized running cost; establishing constraint conditions including the contents of a multi-class equipment output calculation equation, range constraint, heat supply network constraint, power grid constraint, power balance and the like in consideration of reconstruction; by solving the optimized scheduling model, a comprehensive energy system day-ahead optimized scheduling scheme of statistics and heat supply network reconstruction can be obtained. The method increases the regulation flexibility of the comprehensive energy system through heat supply network reconstruction, increases the consumption of renewable energy sources and relieves the transmission blockage of a power grid. According to the scheduling method, the heat supply network reconstruction is used as an adjusting means in day-ahead optimized scheduling, the flexibility of heat supply network adjustment can be fully exerted, the aims of more consumption of renewable energy sources, reduction of transmission blockage of a power grid and the like are achieved, and the efficient operation of a comprehensive energy system is promoted.
Detailed Description
The invention provides a day-ahead optimization scheduling method for a comprehensive energy system considering heat supply network reconstruction, which comprises the following steps of:
(1) establishing an objective function of day-ahead optimized scheduling of the heat supply network reconstruction comprehensive energy system:
Figure BDA0002506639130000081
where T is the sequence number of the scheduling time, T is the set of scheduling times in one day, i1、i2、i3、i4Respectively numbering a cogeneration unit, a heat supply boiler, a traditional generator and a wind turbine generator in the comprehensive energy system, wherein n is the numbering of a power bus in the comprehensive energy system; kCHPAs a collection of cogeneration units, KHBFor integration of boilers, KTUAs a collection of conventional generators, KWDAs a collection of wind turbines, KbusIs a collection of power buses;
Figure BDA0002506639130000082
and
Figure BDA0002506639130000083
respectively a cogeneration unit i1Active power and thermal power generated at the scheduling time t;
Figure BDA0002506639130000084
to characterize cogeneration units i1Binary variable of running state at scheduling time t, cogeneration unit i1Is in the run state at the scheduled time t,
Figure BDA0002506639130000085
value of 1, cogeneration unit i1At a scheduling time tIn the power-off state, the power-on state,
Figure BDA0002506639130000086
is 0;
Figure BDA0002506639130000087
to characterize cogeneration units i1Binary variable of starting-up command at scheduling time t, cogeneration unit i1The power-on command is executed at the scheduled time t,
Figure BDA0002506639130000088
value of 1, cogeneration unit i1The power-on command is not executed at the scheduling time t,
Figure BDA0002506639130000089
is 0;
Figure BDA00025066391300000810
to characterize cogeneration units i1Binary variable of shutdown instruction at scheduling time t, cogeneration unit i1The shutdown instruction is executed at the scheduled time t,
Figure BDA00025066391300000811
value of 1, cogeneration unit i1The shutdown instruction is not executed at the scheduled time t,
Figure BDA0002506639130000091
is 0;
Figure BDA0002506639130000092
for heating boilers i2Thermal power generated at the scheduling time t;
Figure BDA0002506639130000093
for a conventional generator i3Active power generated at a scheduling time t;
Figure BDA0002506639130000094
to characterize a conventional generator i3Binary of running state at scheduling time tSystem variable, conventional generator i3Is in the run state at the scheduled time t,
Figure BDA0002506639130000095
value of 1, conventional generator i3Is in a power-off state at a scheduled time t,
Figure BDA0002506639130000096
is 0;
Figure BDA0002506639130000097
to characterize a conventional generator i3Binary variable of start-up command at scheduling time t, conventional generator i3The power-on command is executed at the scheduled time t,
Figure BDA0002506639130000098
value of 1, conventional generator i3The power-on command is not executed at the scheduling time t,
Figure BDA0002506639130000099
is 0;
Figure BDA00025066391300000910
to characterize a conventional generator i3Binary variable of shutdown command at scheduling time t, conventional generator i3The shutdown instruction is executed at the scheduled time t,
Figure BDA00025066391300000911
value of 1, conventional generator i3The shutdown instruction is not executed at the scheduled time t,
Figure BDA00025066391300000912
is 0;
Figure BDA00025066391300000913
for wind turbine generator i4The active power generated at the scheduling instant t,
Figure BDA00025066391300000914
for the power bus n isActive power of load shedding at a scheduling moment t;
Figure BDA00025066391300000915
for cogeneration units i1The operating cost function of (a);
Figure BDA00025066391300000916
for heating boilers i2The operating cost function of (a);
Figure BDA00025066391300000917
for a conventional generator i3The operating cost function of (a);
Figure BDA00025066391300000918
for wind turbine generator i4A cost function of operation;
Figure BDA00025066391300000919
an operating cost function for shedding the power bus n;
(2) establishing a constraint condition of day-ahead optimized scheduling of the heat supply network reconstruction comprehensive energy system:
(2-1) establishing active power and thermal power constraints of a cogeneration unit in the comprehensive energy system:
Figure BDA00025066391300000920
Figure BDA00025066391300000921
Figure BDA00025066391300000922
wherein k is the serial number of the polygon end point of the feasible region of the cogeneration unit,
Figure BDA00025066391300000923
for cogeneration units i1The total number of end points of the feasible region polygon,
Figure BDA00025066391300000924
for cogeneration units i1The weight of the k-th endpoint in the feasible domain,
Figure BDA00025066391300000925
for cogeneration units i1The active power value of the k-th endpoint of the polygon in the feasible region,
Figure BDA00025066391300000926
for cogeneration units i1A heat value of a k-th endpoint of the feasible region polygon;
(2-2) establishing heat power constraint of a heat supply boiler in the comprehensive energy system:
Figure BDA00025066391300000927
Figure BDA0002506639130000101
in the formula (I), the compound is shown in the specification,
Figure BDA0002506639130000102
for heating boilers i2The thermal efficiency of the heat-generating element,
Figure BDA0002506639130000103
for entering the heating boiler i at the scheduled time t2The power of the natural gas (es) of (c),
Figure BDA0002506639130000104
for heating boilers i2Maximum value of thermal power;
(2-3) establishing a heat power balance equation and a range constraint in the comprehensive energy system:
Figure BDA0002506639130000105
Figure BDA0002506639130000106
Figure BDA0002506639130000107
Figure BDA0002506639130000108
Figure BDA0002506639130000109
Figure BDA00025066391300001010
wherein j is the number of the heating power station consisting of the cogeneration unit and the heating boiler in the step (1), and KHSAs a collection of thermal power stations, Kj CHPIs the set of all cogeneration units in the heating power station j, Kj HBBeing the collection of all the heating boilers in the heating station j,
Figure BDA00025066391300001011
the thermal power provided for the thermal station j at the scheduling time t, b is the number of the thermal pipeline, KPIs the collection of all the thermal power conduits,
Figure BDA00025066391300001012
is the thermal power at the outlet of the thermal pipeline b,
Figure BDA00025066391300001013
the thermal power of the inlet of the thermal pipeline b,
Figure BDA00025066391300001014
the thermal power lost in the thermal power pipeline b, h is the number of the heat supply network nodes, KNFor the set of all the heat network nodes,
Figure BDA00025066391300001015
for the set of thermal conduits with the heat network node h as the end point,
Figure BDA00025066391300001016
for a collection of thermal conduits with the heat network node h as the starting point,
Figure BDA00025066391300001017
for a collection of thermal stations connected to a heat supply network node h, i5The numbers of the heat exchange stations in the comprehensive energy system,
Figure BDA00025066391300001018
being a collection of heat exchange stations connected to the heat network node h,
Figure BDA00025066391300001019
to heat exchange station i5At the moment of scheduling the thermal power at the instant t,
Figure BDA00025066391300001020
for cogeneration units i1The lower limit of the thermal power is,
Figure BDA00025066391300001021
for cogeneration units i1Upper limit of thermal power, ubtIs a binary variable for judging whether the thermal pipeline b is communicated at the scheduling time t, the thermal pipeline b is in a communicated state at the scheduling time t, and u isb,tIs 1, the thermal pipe b is in the closed state at the scheduling instant t, ub,tThe value of (a) is 0,
Figure BDA00025066391300001022
the upper limit of the thermal power is delivered to the thermal pipeline b;
(2-4) establishing range constraint of the heat supply network reconstruction state in the comprehensive energy system:
Figure BDA0002506639130000111
Figure BDA0002506639130000112
Figure BDA0002506639130000113
Figure BDA0002506639130000114
in the formula ub,t-1A binary variable u indicating whether the thermal pipeline b is communicated at the scheduling time t-1b,t-1Value-taking rule of (2) and step (3)b,tThe value-taking rules are the same; x is the number ofb,tTo characterize the binary variable of the opening instruction of thermal pipe b at scheduling time t, thermal pipe b executes the opening instruction at scheduling time t, xb,tIs 1, the thermal pipe b does not execute an opening command at the scheduling time t, xb,tIs 0; y isb,tTo characterize the binary variable of the closing instruction of thermal pipe b at scheduling time t, thermal pipe b executes the closing instruction at scheduling time t, yb,tIs 1, the thermal pipe b does not execute a closing command at the scheduling instant t, yb,tIs 0; x is the number ofb,σBinary variable, x, characterizing the command for opening the thermal pipe b at the scheduling instant sigmab,σValue rule of (1) and xb,tThe value-taking rules are the same; y isb,σTo characterize the binary variable of the closing command of the thermodynamic pipe b at the scheduling instant sigma, yb,σValue rule of (a) and yb,tThe value-taking rules are the same; MRbFor a minimum operating time interval of the electric valve on the thermal pipe b,
Figure BDA0002506639130000115
the initial state of the thermal pipeline b at the scheduling time t;
(2-5) establishing an active power balance equation in the comprehensive energy system:
Figure BDA0002506639130000116
in the formula, Dn,tThe active power consumed by the electric load on the electric bus n at the scheduling time t;
(2-6) establishing range constraint of active power transmission in the comprehensive energy system:
Figure BDA0002506639130000117
wherein l is the number of the power line in the integrated energy system, KLFor the aggregate of all power lines, SFl,nFor power transfer distribution factor, SFl,nIs obtained from the comprehensive energy management system of the energy,
Figure BDA0002506639130000118
being a collection of conventional generators connected on the power bus n,
Figure BDA0002506639130000121
for the set of cogeneration units connected to the power bus n,
Figure BDA0002506639130000122
for a collection of wind turbines connected to a power bus n, FlAn upper limit of the active power transmitted for the power line l;
(2-7) establishing range constraint of active power in the comprehensive energy system:
Figure BDA0002506639130000123
Figure BDA0002506639130000124
Figure BDA0002506639130000125
in the formula (I), the compound is shown in the specification,
Figure BDA0002506639130000126
and
Figure BDA0002506639130000127
respectively a cogeneration unit i1The lower and upper limits of the active power generated,
Figure BDA0002506639130000128
and
Figure BDA0002506639130000129
respectively a conventional generator i3The lower and upper limits of the active power generated,
Figure BDA00025066391300001210
for wind turbine generator i4An upper limit of active power generated;
(2-8) establishing the range constraint of active power climbing in the comprehensive energy system:
Figure BDA00025066391300001211
Figure BDA00025066391300001212
Figure BDA00025066391300001213
Figure BDA00025066391300001214
in the formula (I), the compound is shown in the specification,
Figure BDA00025066391300001215
to characterize cogeneration units i1A binary variable of the running state at the scheduling time t +1,
Figure BDA00025066391300001216
value-taking rule of (1)
Figure BDA00025066391300001217
The value-taking rules are the same;
Figure BDA00025066391300001218
and
Figure BDA00025066391300001219
respectively a cogeneration unit i1The power-off rate and the power-on rate of,
Figure BDA00025066391300001220
and
Figure BDA00025066391300001221
respectively a cogeneration unit i1The downward climbing rate and the upward climbing rate;
Figure BDA00025066391300001222
for cogeneration units i1Active power generated at scheduling time t + 1;
Figure BDA00025066391300001223
to characterize a conventional generator i3A binary variable of the running state at the scheduling time t +1,
Figure BDA00025066391300001224
value-taking rule of (1)
Figure BDA00025066391300001225
The value-taking rules are the same;
Figure BDA00025066391300001226
and
Figure BDA00025066391300001227
respectively a conventional generator i3The power-off rate and the power-on rate of,
Figure BDA00025066391300001228
and
Figure BDA00025066391300001229
respectively a conventional generator i3The downward climbing rate and the upward climbing rate;
Figure BDA00025066391300001230
for a conventional generator i3Active power generated at scheduling time t + 1;
(2-9) establishing range constraint of active power rotation standby in the comprehensive energy system:
Figure BDA00025066391300001231
Figure BDA0002506639130000131
in the formula (I), the compound is shown in the specification,
Figure BDA0002506639130000132
for each scheduling time
Figure BDA0002506639130000133
Maximum value of, SRUtReserve capacity, SRD, for an integrated energy system at scheduled time ttotaitThe total downward reserve capacity of the comprehensive energy system at the scheduling time t;
(2-10) establishing unit start-stop logic constraint in the comprehensive energy system:
Figure BDA0002506639130000134
Figure BDA0002506639130000135
in the formula (I), the compound is shown in the specification,
Figure BDA0002506639130000136
to characterize cogeneration units i1A binary variable of the running state at the scheduling time t-1,
Figure BDA0002506639130000137
value-taking rule of (1)
Figure BDA0002506639130000138
Value ofThe rules are the same;
Figure BDA0002506639130000139
to characterize cogeneration units i1The binary variable of the shutdown instruction at scheduling time t-1,
Figure BDA00025066391300001310
value-taking rule of (1)
Figure BDA00025066391300001311
The value-taking rules are the same;
Figure BDA00025066391300001312
to characterize a conventional generator i3A binary variable of the running state at the scheduling time t-1,
Figure BDA00025066391300001313
value-taking rule of (1)
Figure BDA00025066391300001314
The value-taking rules are the same;
Figure BDA00025066391300001315
to characterize a conventional generator i3The binary variable of the shutdown instruction at scheduling time t-1,
Figure BDA00025066391300001316
value-taking rule of (1)
Figure BDA00025066391300001317
The value-taking rules are the same;
(2-11) establishing a range constraint of the minimum startup time and the minimum shutdown time in the comprehensive energy system:
Figure BDA00025066391300001318
Figure BDA00025066391300001319
Figure BDA00025066391300001320
Figure BDA00025066391300001321
in the formula (I), the compound is shown in the specification,
Figure BDA00025066391300001322
and
Figure BDA00025066391300001323
respectively a cogeneration unit i1The minimum on-time and the minimum off-time,
Figure BDA00025066391300001324
and
Figure BDA00025066391300001325
respectively a conventional generator i3Minimum turn-on time and minimum turn-off time;
(3) the method comprises the steps of (1) forming a comprehensive energy system day-ahead optimization scheduling model considering heat supply network reconstruction by the objective function of the heat supply network reconstruction comprehensive energy system day-ahead optimization scheduling established in the step (1) and the constraint conditions of the heat supply network reconstruction comprehensive energy system day-ahead optimization scheduling established in the step (2), solving the optimization model by adopting a branch-and-bound method to obtain the comprehensive energy system day-ahead optimization scheduling parameters considering the heat supply network reconstruction, and including a cogeneration unit i in the step (1)1Active power generated at scheduling time t
Figure BDA0002506639130000141
And thermal power
Figure BDA0002506639130000142
Characterization of cogeneration units i1Binary variable of running state at scheduling time t
Figure BDA0002506639130000143
Characterization of cogeneration units i1Binary variable of starting instruction at scheduling time t
Figure BDA0002506639130000144
Characterization of cogeneration units i1Binary variable of shutdown instruction at scheduling time t
Figure BDA0002506639130000145
Heating boiler i2Thermal power generated at scheduled time t
Figure BDA0002506639130000146
Conventional generator i3Active power generated at scheduling time t
Figure BDA0002506639130000147
Characterizing a conventional Generator i3Binary variable of running state at scheduling time t
Figure BDA0002506639130000148
Characterizing a conventional Generator i3Binary variable of starting instruction at scheduling time t
Figure BDA0002506639130000149
Characterizing a conventional Generator i3Binary variable of shutdown instruction at scheduling time t
Figure BDA00025066391300001410
Binary variable u for judging whether thermal pipeline b is communicated at scheduling time t or notb,tBinary variable x for representing opening instruction of thermal pipeline b at scheduling time tb,tBinary variable y for representing closing instruction of thermal pipeline b at scheduling time tb,t。
This optimization model may be solved by a Gurobi commercial solver in one embodiment of the invention.

Claims (1)

1. A day-ahead optimization scheduling method of an integrated energy system considering heat supply network reconstruction is characterized by comprising the following steps:
(1) establishing an objective function of day-ahead optimized scheduling of the heat supply network reconstruction comprehensive energy system:
Figure FDA0002506639120000011
where T is the sequence number of the scheduling time, T is the set of scheduling times in one day, i1、i2、i3、i4Respectively numbering a cogeneration unit, a heat supply boiler, a traditional generator and a wind turbine generator in the comprehensive energy system, wherein n is the numbering of a power bus in the comprehensive energy system; kCHPAs a collection of cogeneration units, KHBFor integration of boilers, KTUAs a collection of conventional generators, KWDAs a collection of wind turbines, KbusIs a collection of power buses;
Figure FDA0002506639120000012
and
Figure FDA0002506639120000013
respectively a cogeneration unit i1Active power and thermal power generated at the scheduling time t;
Figure FDA0002506639120000014
to characterize cogeneration units i1Binary variable of running state at scheduling time t, cogeneration unit i1Is in the run state at the scheduled time t,
Figure FDA0002506639120000015
value of 1, cogeneration unit i1Is in a power-off state at a scheduled time t,
Figure FDA0002506639120000016
is 0;
Figure FDA0002506639120000017
to characterize cogeneration units i1Binary variable of starting-up command at scheduling time t, cogeneration unit i1The power-on command is executed at the scheduled time t,
Figure FDA0002506639120000018
value of 1, cogeneration unit i1The power-on command is not executed at the scheduling time t,
Figure FDA0002506639120000019
is 0;
Figure FDA00025066391200000110
to characterize cogeneration units i1Binary variable of shutdown instruction at scheduling time t, cogeneration unit i1The shutdown instruction is executed at the scheduled time t,
Figure FDA00025066391200000111
value of 1, cogeneration unit i1The shutdown instruction is not executed at the scheduled time t,
Figure FDA00025066391200000112
is 0;
Figure FDA00025066391200000113
for heating boilers i2Thermal power generated at the scheduling time t;
Figure FDA00025066391200000114
for a conventional generator i3Active power generated at a scheduling time t;
Figure FDA00025066391200000115
to characterize a conventional generator i3Binary variable of running state at scheduling time t, conventional generator i3Is in the run state at the scheduled time t,
Figure FDA00025066391200000116
value of 1, conventional generator i3Is in a power-off state at a scheduled time t,
Figure FDA00025066391200000117
is 0;
Figure FDA00025066391200000118
to characterize a conventional generator i3Binary variable of start-up command at scheduling time t, conventional generator i3The power-on command is executed at the scheduled time t,
Figure FDA00025066391200000119
value of 1, conventional generator i3The power-on command is not executed at the scheduling time t,
Figure FDA00025066391200000120
is 0;
Figure FDA00025066391200000121
to characterize a conventional generator i3Binary variable of shutdown command at scheduling time t, conventional generator i3The shutdown instruction is executed at the scheduled time t,
Figure FDA0002506639120000021
value of 1, conventional generator i3The shutdown instruction is not executed at the scheduled time t,
Figure FDA0002506639120000022
is 0;
Figure FDA0002506639120000023
for wind turbine generator i4The active power generated at the scheduling instant t,
Figure FDA0002506639120000024
as a power busn active power at the scheduling time t for load shedding;
Figure FDA0002506639120000025
for cogeneration units i1The operating cost function of (a);
Figure FDA0002506639120000026
for heating boilers i2The operating cost function of (a);
Figure FDA0002506639120000027
for a conventional generator i3The operating cost function of (a);
Figure FDA0002506639120000028
for wind turbine generator i4A cost function of operation;
Figure FDA0002506639120000029
an operating cost function for shedding the power bus n;
(2) establishing a constraint condition of day-ahead optimized scheduling of the heat supply network reconstruction comprehensive energy system:
(2-1) establishing active power and thermal power constraints of a cogeneration unit in the comprehensive energy system:
Figure FDA00025066391200000210
Figure FDA00025066391200000211
Figure FDA00025066391200000212
wherein k is the serial number of the polygon end point of the feasible region of the cogeneration unit,
Figure FDA00025066391200000213
as a heat and power unitProduct unit i1The total number of end points of the feasible region polygon,
Figure FDA00025066391200000214
for cogeneration units i1The weight of the k-th endpoint in the feasible domain,
Figure FDA00025066391200000215
for cogeneration units i1The active power value of the k-th endpoint of the polygon in the feasible region,
Figure FDA00025066391200000216
for cogeneration units i1A heat value of a k-th endpoint of the feasible region polygon;
(2-2) establishing heat power constraint of a heat supply boiler in the comprehensive energy system:
Figure FDA00025066391200000217
Figure FDA00025066391200000218
in the formula (I), the compound is shown in the specification,
Figure FDA00025066391200000219
for heating boilers i2The thermal efficiency of the heat-generating element,
Figure FDA00025066391200000220
for entering the heating boiler i at the scheduled time t2The power of the natural gas (es) of (c),
Figure FDA00025066391200000221
for heating boilers i2Maximum value of thermal power;
(2-3) establishing a heat power balance equation and a range constraint in the comprehensive energy system:
Figure FDA00025066391200000222
Figure FDA00025066391200000223
Figure FDA0002506639120000031
Figure FDA0002506639120000032
Figure FDA0002506639120000033
Figure FDA0002506639120000034
wherein j is the number of the heating power station consisting of the cogeneration unit and the heating boiler in the step (1), and KHSAs a collection of thermal power stations, Kj CHPIs the set of all cogeneration units in the heating power station j, Kj HBBeing the collection of all the heating boilers in the heating station j,
Figure FDA0002506639120000035
the thermal power provided for the thermal station j at the scheduling time t, b is the number of the thermal pipeline, KPIs the collection of all the thermal power conduits,
Figure FDA0002506639120000036
is the thermal power at the outlet of the thermal pipeline b,
Figure FDA0002506639120000037
the thermal power of the inlet of the thermal pipeline b,
Figure FDA0002506639120000038
is heat powerThe thermal power lost in the pipeline b, h is the number of the heat supply network nodes, KNFor the set of all the heat network nodes,
Figure FDA0002506639120000039
for the set of thermal conduits with the heat network node h as the end point,
Figure FDA00025066391200000310
for a collection of thermal conduits with the heat network node h as the starting point,
Figure FDA00025066391200000311
for a collection of thermal stations connected to a heat supply network node h, i5The numbers of the heat exchange stations in the comprehensive energy system,
Figure FDA00025066391200000312
being a collection of heat exchange stations connected to the heat network node h,
Figure FDA00025066391200000313
to heat exchange station i5At the moment of scheduling the thermal power at the instant t,
Figure FDA00025066391200000314
for cogeneration units i1The lower limit of the thermal power is,
Figure FDA00025066391200000315
for cogeneration units i1Upper limit of thermal power, ub,tIs a binary variable for judging whether the thermal pipeline b is communicated at the scheduling time t, the thermal pipeline b is in a communicated state at the scheduling time t, and u isb,tIs 1, the thermal pipe b is in the closed state at the scheduling instant t, ub,tThe value of (a) is 0,
Figure FDA00025066391200000316
the upper limit of the thermal power is delivered to the thermal pipeline b;
(2-4) establishing range constraint of the heat supply network reconstruction state in the comprehensive energy system:
Figure FDA00025066391200000317
Figure FDA00025066391200000318
Figure FDA00025066391200000319
Figure FDA00025066391200000320
in the formula ub,t-1A binary variable u indicating whether the thermal pipeline b is communicated at the scheduling time t-1b,t-1Value-taking rule of (2) and step (3)b,tThe value-taking rules are the same; x is the number ofb,tTo characterize the binary variable of the opening instruction of thermal pipe b at scheduling time t, thermal pipe b executes the opening instruction at scheduling time t, xb,tIs 1, the thermal pipe b does not execute an opening command at the scheduling time t, xb,tIs 0; y isb,tTo characterize the binary variable of the closing instruction of thermal pipe b at scheduling time t, thermal pipe b executes the closing instruction at scheduling time t, yb,tIs 1, the thermal pipe b does not execute a closing command at the scheduling instant t, yb,tIs 0; x is the number ofb,σBinary variable, x, characterizing the command for opening the thermal pipe b at the scheduling instant sigmab,σValue rule of (1) and xb,tThe value-taking rules are the same; y isb,σTo characterize the binary variable of the closing command of the thermodynamic pipe b at the scheduling instant sigma, yb,σValue rule of (a) and yb,tThe value-taking rules are the same; MRbFor a minimum operating time interval of the electric valve on the thermal pipe b,
Figure FDA0002506639120000041
the initial state of the thermal pipeline b at the scheduling time t;
(2-5) establishing an active power balance equation in the comprehensive energy system:
Figure FDA0002506639120000042
in the formula, Dn,tThe active power consumed by the electric load on the electric bus n at the scheduling time t;
(2-6) establishing range constraint of active power transmission in the comprehensive energy system:
Figure FDA0002506639120000043
wherein l is the number of the power line in the integrated energy system, KLFor the aggregate of all power lines, SFl,nFor power transfer distribution factor, SFl,nIs obtained from the comprehensive energy management system of the energy,
Figure FDA0002506639120000044
being a collection of conventional generators connected on the power bus n,
Figure FDA0002506639120000045
for the set of cogeneration units connected to the power bus n,
Figure FDA0002506639120000046
for a collection of wind turbines connected to a power bus n, FlAn upper limit of the active power transmitted for the power line l;
(2-7) establishing range constraint of active power in the comprehensive energy system:
Figure FDA0002506639120000047
Figure FDA0002506639120000048
Figure FDA0002506639120000049
in the formula (I), the compound is shown in the specification,
Figure FDA00025066391200000410
and
Figure FDA00025066391200000411
respectively a cogeneration unit i1The lower and upper limits of the active power generated,
Figure FDA00025066391200000412
and
Figure FDA00025066391200000413
respectively a conventional generator i3The lower and upper limits of the active power generated,
Figure FDA0002506639120000051
for wind turbine generator i4An upper limit of active power generated;
(2-8) establishing the range constraint of active power climbing in the comprehensive energy system:
Figure FDA0002506639120000052
Figure FDA0002506639120000053
Figure FDA0002506639120000054
Figure FDA0002506639120000055
in the formula (I), the compound is shown in the specification,
Figure FDA0002506639120000056
to characterize cogeneration units i1A binary variable of the running state at the scheduling time t +1,
Figure FDA0002506639120000057
value-taking rule of (1)
Figure FDA0002506639120000058
The value-taking rules are the same;
Figure FDA0002506639120000059
and
Figure FDA00025066391200000510
respectively a cogeneration unit i1The power-off rate and the power-on rate of,
Figure FDA00025066391200000511
and
Figure FDA00025066391200000512
respectively a cogeneration unit i1The downward climbing rate and the upward climbing rate;
Figure FDA00025066391200000513
for cogeneration units i1Active power generated at scheduling time t + 1;
Figure FDA00025066391200000514
to characterize a conventional generator i3A binary variable of the running state at the scheduling time t +1,
Figure FDA00025066391200000515
value-taking rule of (1)
Figure FDA00025066391200000516
The value-taking rules are the same;
Figure FDA00025066391200000517
and
Figure FDA00025066391200000518
respectively a conventional generator i3The power-off rate and the power-on rate of,
Figure FDA00025066391200000519
and
Figure FDA00025066391200000520
respectively a conventional generator i3The downward climbing rate and the upward climbing rate;
Figure FDA00025066391200000521
for a conventional generator i3Active power generated at scheduling time t + 1;
(2-9) establishing range constraint of active power rotation standby in the comprehensive energy system:
Figure FDA00025066391200000522
Figure FDA00025066391200000523
in the formula (I), the compound is shown in the specification,
Figure FDA00025066391200000524
for each scheduling time
Figure FDA00025066391200000525
Maximum value of, SRUtReserve capacity, SRD, for an integrated energy system at scheduled time ttotaitThe total downward reserve capacity of the comprehensive energy system at the scheduling time t;
(2-10) establishing unit start-stop logic constraint in the comprehensive energy system:
Figure FDA00025066391200000526
Figure FDA00025066391200000527
in the formula (I), the compound is shown in the specification,
Figure FDA00025066391200000528
to characterize cogeneration units i1A binary variable of the running state at the scheduling time t-1,
Figure FDA00025066391200000529
value-taking rule of (1)
Figure FDA0002506639120000061
The value-taking rules are the same;
Figure FDA0002506639120000062
to characterize cogeneration units i1The binary variable of the shutdown instruction at scheduling time t-1,
Figure FDA0002506639120000063
value-taking rule of (1)
Figure FDA0002506639120000064
The value-taking rules are the same;
Figure FDA0002506639120000065
to characterize a conventional generator i3A binary variable of the running state at the scheduling time t-1,
Figure FDA0002506639120000066
value-taking rule of (1)
Figure FDA0002506639120000067
The value-taking rules are the same;
Figure FDA0002506639120000068
to characterize a conventional generator i3The binary variable of the shutdown instruction at scheduling time t-1,
Figure FDA0002506639120000069
value-taking rule of (1)
Figure FDA00025066391200000610
The value-taking rules are the same;
(2-11) establishing a range constraint of the minimum startup time and the minimum shutdown time in the comprehensive energy system:
Figure FDA00025066391200000611
Figure FDA00025066391200000612
Figure FDA00025066391200000613
Figure FDA00025066391200000614
in the formula (I), the compound is shown in the specification,
Figure FDA00025066391200000615
and
Figure FDA00025066391200000616
respectively a cogeneration unit i1The minimum on-time and the minimum off-time,
Figure FDA00025066391200000617
and
Figure FDA00025066391200000618
respectively a conventional generator i3Minimum turn-on time and minimum turn-off time;
(3) the method comprises the steps of (1) forming a comprehensive energy system day-ahead optimization scheduling model considering heat supply network reconstruction by the objective function of the heat supply network reconstruction comprehensive energy system day-ahead optimization scheduling established in the step (1) and the constraint conditions of the heat supply network reconstruction comprehensive energy system day-ahead optimization scheduling established in the step (2), solving the optimization model by adopting a branch-and-bound method to obtain the comprehensive energy system day-ahead optimization scheduling parameters considering the heat supply network reconstruction, and including a cogeneration unit i in the step (1)1Active power generated at scheduling time t
Figure FDA00025066391200000619
And thermal power
Figure FDA00025066391200000620
Characterization of cogeneration units i1Binary variable of running state at scheduling time t
Figure FDA00025066391200000621
Characterization of cogeneration units i1Binary variable of starting instruction at scheduling time t
Figure FDA00025066391200000622
Characterization of cogeneration units i1Binary variable of shutdown instruction at scheduling time t
Figure FDA00025066391200000623
Heating boiler i2Thermal power generated at scheduled time t
Figure FDA00025066391200000624
Conventional generator i3Active power generated at scheduling time t
Figure FDA00025066391200000625
Characterizing a conventional Generator i3Binary variable of running state at scheduling time t
Figure FDA0002506639120000071
Characterizing a conventional Generator i3Binary variable of starting instruction at scheduling time t
Figure FDA0002506639120000072
Characterizing a conventional Generator i3Binary variable of shutdown instruction at scheduling time t
Figure FDA0002506639120000073
Binary variable u for judging whether thermal pipeline b is communicated at scheduling time t or notb,tBinary variable x for representing opening instruction of thermal pipeline b at scheduling time tb,tBinary variable y for representing closing instruction of thermal pipeline b at scheduling time tb,t。
CN202010447996.6A 2020-05-25 2020-05-25 Day-ahead optimization scheduling method for comprehensive energy system considering heat supply network reconstruction Active CN111563634B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010447996.6A CN111563634B (en) 2020-05-25 2020-05-25 Day-ahead optimization scheduling method for comprehensive energy system considering heat supply network reconstruction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010447996.6A CN111563634B (en) 2020-05-25 2020-05-25 Day-ahead optimization scheduling method for comprehensive energy system considering heat supply network reconstruction

Publications (2)

Publication Number Publication Date
CN111563634A true CN111563634A (en) 2020-08-21
CN111563634B CN111563634B (en) 2022-04-08

Family

ID=72069681

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010447996.6A Active CN111563634B (en) 2020-05-25 2020-05-25 Day-ahead optimization scheduling method for comprehensive energy system considering heat supply network reconstruction

Country Status (1)

Country Link
CN (1) CN111563634B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106849188A (en) * 2017-01-23 2017-06-13 中国电力科学研究院 The combined heat and power optimization method and system of a kind of promotion wind electricity digestion
CN107067116A (en) * 2017-04-26 2017-08-18 燕山大学 A kind of multizone electric heating integrated system economic environment combined dispatching method for solving
CN109840638A (en) * 2019-02-28 2019-06-04 华南理工大学 A kind of meter and the probabilistic combined heat and power of heat supply network robust economic load dispatching method a few days ago
CN109978625A (en) * 2019-03-28 2019-07-05 河海大学 It is a kind of meter and electric heating gas network integrated energy system multiple target running optimizatin method
US20190369581A1 (en) * 2017-01-11 2019-12-05 Southeast University Integrated energy system operational optimization method considering thermal inertia of district heating networks and buildings

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190369581A1 (en) * 2017-01-11 2019-12-05 Southeast University Integrated energy system operational optimization method considering thermal inertia of district heating networks and buildings
CN106849188A (en) * 2017-01-23 2017-06-13 中国电力科学研究院 The combined heat and power optimization method and system of a kind of promotion wind electricity digestion
CN107067116A (en) * 2017-04-26 2017-08-18 燕山大学 A kind of multizone electric heating integrated system economic environment combined dispatching method for solving
CN109840638A (en) * 2019-02-28 2019-06-04 华南理工大学 A kind of meter and the probabilistic combined heat and power of heat supply network robust economic load dispatching method a few days ago
CN109978625A (en) * 2019-03-28 2019-07-05 河海大学 It is a kind of meter and electric heating gas network integrated energy system multiple target running optimizatin method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HONGBIN SUN 等: "Preventive control for combined heat and power system: A day-ahead security constrained economic dispatch model", 《IEEE》 *
赵博: "Optimal Scheduling Method for Electrical-Thermal Integrated Energy System Considering Heat Storage Characteristics of Heating Network", 《WEB OF SCIENCE》 *

Also Published As

Publication number Publication date
CN111563634B (en) 2022-04-08

Similar Documents

Publication Publication Date Title
CN108717594B (en) Economic optimization scheduling method for combined cooling heating and power type multi-microgrid system
CN110245878B (en) Distributed comprehensive energy demand response collaborative optimization method for intelligent building group
CN106960282B (en) Coordinated operation method of gas-electricity comprehensive energy distribution network system
CN108625911B (en) Thermodynamic system for improving electric output adjusting capacity of heat supply unit
WO2018130231A1 (en) Integrated energy system optimization method based on heating grid and home thermal inertia
CN106950840B (en) Power grid peak clipping-oriented hierarchical distributed coordination control method for comprehensive energy system
CN110361969B (en) Optimized operation method of cooling, heating and power comprehensive energy system
CN110535128A (en) Based on the multizone integrated energy system coordinated dispatching method with energy comfort level
CN113725915A (en) Rural electric heating comprehensive energy system operation optimization method considering renewable energy uncertainty and thermal inertia
CN105186583B (en) Energy router and its energy dispatching method based on multi-agent modeling
CN114154328A (en) Flexibility-improved two-stage random optimization scheduling method for electric heating comprehensive energy system
CN110224443B (en) Decision-making method for combined start-stop control of generator set and heat exchange station
Zhao et al. Integrated unit commitment and economic dispatch of combined heat and power system considering heat-power decoupling retrofit of CHP unit
CN105470957B (en) Power grid load modeling method for production simulation
CN111563634B (en) Day-ahead optimization scheduling method for comprehensive energy system considering heat supply network reconstruction
CN110991845B (en) Distributed cooperative scheduling method for electric-thermal coupling system
CN107132806A (en) A kind of control method of cold, heat and power triple supply system
CN113610305B (en) Optimized scheduling method of comprehensive energy system
CN115860412A (en) Collaborative planning method and terminal of electric heating coupling system
CN115341964A (en) Thermal power generating unit heat storage peak regulation power generation system and method with desalting water tank
CN108155644A (en) Meter and the dynamic heating system of ancillary service and electric system coupled scheduler method
CN111313410B (en) Equipment regulation and control device of comprehensive energy system
CN113251459B (en) Method and system for deep peak shaving of thermoelectric unit
Xi et al. Robust optimization of integrated energy system considering multiple thermal inertia and wind power uncertainty
CN205141679U (en) Comprehensive energy power supply 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
GR01 Patent grant