CN112734591A - Electric heating comprehensive coordination scheduling method and device, equipment and medium - Google Patents

Electric heating comprehensive coordination scheduling method and device, equipment and medium Download PDF

Info

Publication number
CN112734591A
CN112734591A CN202011357166.0A CN202011357166A CN112734591A CN 112734591 A CN112734591 A CN 112734591A CN 202011357166 A CN202011357166 A CN 202011357166A CN 112734591 A CN112734591 A CN 112734591A
Authority
CN
China
Prior art keywords
thermal
electric heating
unit
constructing
pressure cylinder
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
CN202011357166.0A
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.)
Inner Mongolia Fengtai Power Generation Co ltd
Tsinghua University
North United Power Co Ltd
Original Assignee
Inner Mongolia Fengtai Power Generation Co ltd
Tsinghua University
North United 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 Inner Mongolia Fengtai Power Generation Co ltd, Tsinghua University, North United Power Co Ltd filed Critical Inner Mongolia Fengtai Power Generation Co ltd
Priority to CN202011357166.0A priority Critical patent/CN112734591A/en
Publication of CN112734591A publication Critical patent/CN112734591A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/04Power grid distribution networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides an electric heating comprehensive coordination scheduling method, an electric heating comprehensive coordination scheduling device, electric heating comprehensive coordination scheduling equipment and an electric heating comprehensive coordination scheduling medium, wherein the method comprises the following steps: obtaining model parameters of an electric power system, model parameters of a thermal power plant and model parameters of a thermal system: constructing an electric heating comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters, wherein the electric heating comprehensive coordination scheduling model is a model of electric heating; and obtaining a solving result based on the electric heating comprehensive coordination scheduling model so as to carry out scheduling. By introducing a detailed thermal power plant model into the electric heating comprehensive coordination scheduling problem, the physical process in the thermal power plant can be more accurately described, the real physical law of electric heating system coupling is reflected, and scheduling is more reasonably performed.

Description

Electric heating comprehensive coordination scheduling method and device, equipment and medium
Technical Field
The invention relates to the technical field of energy utilization, in particular to an electric heating comprehensive coordination scheduling method, an electric heating comprehensive coordination scheduling device, electric heating comprehensive coordination scheduling equipment and an electric heating comprehensive coordination scheduling medium.
Background
At present, the problems of energy safety and environmental pollution are increasingly prominent, the use mode of the traditional fossil energy is changed, and the improvement of the energy utilization efficiency gradually becomes the urgent need of the human society. The cogeneration technology is a key technology for realizing the cascade utilization of energy, and the development of the cogeneration technology in China is very rapid in recent years. With the gradual development of thermal power plants, the coupling between the traditional power system and the thermodynamic system is increasingly tight, and an electric-heat comprehensive energy system is gradually developed. In order to fully exert the synergistic effect of the electric heating comprehensive energy system on realizing the optimal utilization of energy, the current electric heating comprehensive coordination scheduling problem becomes a hot spot problem which is widely concerned.
In the existing electric heating comprehensive coordination scheduling problem, the key electric heating coupling loop of the thermal power plant is usually simplified into a polygonal operation area of one or more thermoelectric power units, however, in fact, the thermal power plant usually has a more complex internal structure, and the adjustment of the loop involves various key parameters of a plurality of thermal components in the thermal power plant; on the other hand, the increasing abundance of the adjusting means of the thermal power plant causes the change of the unit operation interval, and the traditional modeling mode based on the operation interval cannot meet the analysis requirement. Then, it is necessary to model the thermal power plant in more detail in the electric heating comprehensive coordination scheduling problem, which is the basis for objectively reflecting the electric heating system coupling law.
Disclosure of Invention
The invention provides an electric heating comprehensive coordination scheduling method, an electric heating comprehensive coordination scheduling device and a medium.
In a first aspect, the present invention provides an electric heating comprehensive coordination scheduling method, including:
obtaining model parameters of an electric power system, model parameters of a thermal power plant and model parameters of a thermal system:
constructing an electric heating comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters, wherein the electric heating comprehensive coordination scheduling model is a model of electric heating;
and obtaining a solving result based on the electric heating comprehensive coordination scheduling model so as to carry out scheduling.
Further, the power system model parameters comprise a coal consumption characteristic function of a conventional thermal power generating unit, the maximum climbing speed of the conventional thermal power generating unit and the thermoelectric unit, and the maximum and minimum technical output of the conventional thermal power generating unit; the method comprises the following steps of (1) connecting a power bus and a power transmission line, active power transfer distribution factors from the power bus to the power transmission line, the maximum active transmission capacity of the power transmission line, an electric load prediction curve in a scheduling time interval and a wind power active output limit prediction curve;
the model parameters of the thermal power plant comprise topological connection relations among thermal parts, inlet and outlet temperatures, inlet and outlet pressures, steam inlet flow, ideal specific enthalpy drop, internal efficiency, maximum and minimum main steam flow of the thermal power plant, high-pressure cylinder bypass maximum flow, low-pressure cylinder minimum steam inlet flow, sliding pressure curve shape parameters, temperature-reduced water and boiler feed water specific enthalpy, generator efficiency, mechanical efficiency, boiler efficiency, physical property database, specific heat capacity of water, coal heat value, phase-change heat exchanger and performance coefficients of an absorption heat pump of the thermal power plant under rated working conditions;
the thermodynamic system model parameters comprise a node-pipeline connection relation of a primary thermodynamic pipeline, a temperature dynamic equation coefficient, transmission delay, mass flow, maximum and minimum node temperatures, soil temperature, minimum equivalent thermal resistance of a secondary thermodynamic pipeline network, maximum and minimum indoor temperatures, thermal capacity and thermal resistance of a building enclosure structure and an outdoor environment temperature prediction curve.
Further, the constructing of the electric-thermal comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters comprises:
constructing an objective function of an electric heating comprehensive coordination scheduling model, wherein the expression is as follows:
Figure RE-GDA0002966498820000031
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000032
the active power output of the conventional thermal power generating unit i in the time period t is shown;
Figure RE-GDA0002966498820000033
the coal consumption characteristic function of a conventional thermal power generating unit i is obtained;
Figure RE-GDA0002966498820000034
the coal consumption of the thermoelectric unit i in the time period t is shown;
constructing a power system constraint condition of an electric heating comprehensive coordination scheduling model;
constructing thermal power plant constraint conditions of an electric heating comprehensive coordination scheduling model; the thermal power plant comprises one or more of thermal parts or thermal links such as a thermal motor set, a phase-change heat exchanger, an absorption heat pump, an in-plant heat supply network and the like; the thermoelectric unit comprises one or more types of back pressure units or extraction condensing units;
constructing a thermodynamic system constraint condition of an electric heating comprehensive coordination scheduling model;
and constructing an abstract form of the electric heating comprehensive coordination scheduling model.
Further, the power system constraints include: the method comprises the following steps of (1) power system energy balance constraint, unit active output upper and lower limit constraint, unit climbing rate constraint and line active transmission capacity constraint;
the power system energy balance constraint is as follows:
Figure RE-GDA0002966498820000035
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000036
the active power of the thermoelectric unit i in the time period t;
Figure RE-GDA0002966498820000037
is the active power output of the wind farm w at time t;
Figure RE-GDA0002966498820000038
is a predicted value of the electrical load e at the time period t;
the active output upper and lower limits of the unit are constrained as follows:
Figure RE-GDA0002966498820000039
wherein the content of the first and second substances,
Figure RE-GDA00029664988200000310
and
Figure RE-GDA00029664988200000311
respectively the maximum and minimum technical output of a conventional thermal power generating unit i;
Figure RE-GDA00029664988200000312
the active output limit predicted value of the wind power plant w in the time period t is obtained;
the unit climbing rate constraint is as follows:
Figure RE-GDA00029664988200000313
wherein the content of the first and second substances,
Figure RE-GDA00029664988200000314
and
Figure RE-GDA00029664988200000315
the maximum climbing speed and the maximum climbing speed of the conventional thermal power generating unit i are respectively set;
Figure RE-GDA0002966498820000041
and
Figure RE-GDA0002966498820000042
the maximum climbing speed and the maximum climbing speed of the thermoelectric unit i are respectively set; Δ t is the length of the scheduling period;
the active transmission capacity constraint of the line is as follows:
Figure RE-GDA0002966498820000043
wherein, SFulIs the active power transfer distribution factor from the power bus u to the power transmission line l;
Figure RE-GDA0002966498820000044
is the maximum active transmission capacity of the power transmission line l.
Further, the thermal power plant constraint conditions for constructing the electric-thermal comprehensive coordination scheduling model comprise:
constructing operation constraints of high, medium and low pressure cylinders of the thermoelectric unit; for the high pressure cylinder, the middle pressure cylinder and the low pressure cylinder of the back pressure unit and the high pressure cylinder and the middle pressure cylinder of the extraction condensing unit, the following operation constraints are constructed:
Figure RE-GDA0002966498820000045
wherein, the superscript x can be HP, IP and LP, which respectively represent the relevant parameters of high, middle and low pressure cylinders; the superscript 0 represents the relevant parameter under the rated working condition;
Figure BDA0002802877070000045
the steam inlet flow of the high, medium and low pressure cylinders;
Figure BDA0002802877070000046
and
Figure BDA0002802877070000047
the inlet pressure and the outlet pressure of the high-pressure cylinder, the medium-pressure cylinder and the low-pressure cylinder are respectively set;
Figure BDA0002802877070000048
and
Figure BDA0002802877070000049
the inlet and outlet temperatures of the high, medium and low pressure cylinders are respectively;
Figure BDA00028028770700000410
the ideal temperature of the steam after isentropic adiabatic expansion in the high, medium and low pressure cylinders;
Figure BDA00028028770700000411
and
Figure BDA00028028770700000412
the actual specific enthalpy drop in the high-pressure cylinder, the intermediate-pressure cylinder and the low-pressure cylinder and the ideal specific enthalpy drop during isentropic adiabatic expansion are respectively;
Figure BDA00028028770700000413
the internal efficiency of the high, medium and low pressure cylinders under rated working conditions; h (-) and s (-) are a specific enthalpy function and a specific entropy function, respectively, and the function values can be queried through a physical property database; the physical property database is an existing database;
for the low pressure cylinder of the extraction and condensation unit, the following operational constraints are established:
Figure RE-GDA0002966498820000051
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000052
and
Figure RE-GDA0002966498820000053
the actual dryness of the last-stage steam of the low-pressure cylinder and the ideal dryness of the steam of the low-pressure cylinder after isentropic adiabatic expansion are respectively obtained;
constructing coupling constraints among high, medium and low pressure cylinders of the thermoelectric unit:
for the back pressure unit, the following steam flow balance equation is constructed:
Figure RE-GDA0002966498820000054
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000055
the serial number of the phase change heat exchanger connected with the thermoelectric unit i;
Figure RE-GDA0002966498820000056
and
Figure RE-GDA0002966498820000057
respectively the main steam flow and the regenerative steam flow of the thermoelectric unit i in a time period t;
Figure RE-GDA0002966498820000058
is the heating steam flow of the phase change heat exchanger j in the time period t; alpha is alphai、βiAnd gammaiThe final steam flow of the high, medium and low pressure cylinders of the thermoelectric unit i accounts for the proportion of the steam inlet flow;
Figure RE-GDA0002966498820000059
and
Figure RE-GDA00029664988200000510
the maximum main steam flow and the minimum main steam flow of the thermoelectric unit i are respectively;
for the extraction condensing unit, the following steam flow balance equation is constructed:
Figure RE-GDA00029664988200000511
Figure RE-GDA00029664988200000512
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000061
is the serial number of the absorption heat pump connected with the thermoelectric unit i;
Figure RE-GDA0002966498820000062
and
Figure RE-GDA0002966498820000063
the bypass flow of a high-pressure cylinder and the steam extraction flow of a medium-pressure cylinder of the thermoelectric unit i in a time period t are shown;
Figure RE-GDA0002966498820000064
and
Figure RE-GDA0002966498820000065
the maximum flow of a high-pressure cylinder bypass and the minimum steam inlet flow of a low-pressure cylinder of the thermoelectric unit i are respectively;
Figure RE-GDA0002966498820000066
and
Figure RE-GDA0002966498820000067
the proportions of the high-pressure cylinder bypass desuperheating water and the medium-pressure cylinder bypass desuperheating water of the thermoelectric generator set i at the time t are respectively;
Figure RE-GDA0002966498820000068
is the specific enthalpy of the reduced temperature water of the thermoelectric unit i;
Figure RE-GDA0002966498820000069
is the driving steam flow of the absorption heat pump j in the time period t;
and for the back pressure unit and the extraction condensing unit, constructing a coupling relation of inlet and outlet temperatures and pressures of high, medium and low pressure cylinders:
Figure RE-GDA00029664988200000610
wherein the content of the first and second substances,
Figure RE-GDA00029664988200000611
the steam inlet temperature of the high and medium pressure cylinders under rated working conditions;
Figure RE-GDA00029664988200000612
aiand biThe shape parameters of a sliding pressure curve of the thermoelectric unit i are all the shape parameters;
and constructing an external characteristic equation for the back pressure unit and the extraction condensing unit:
Figure RE-GDA00029664988200000613
wherein the content of the first and second substances,
Figure RE-GDA00029664988200000614
and
Figure RE-GDA00029664988200000615
the generator efficiency and the mechanical efficiency of the thermoelectric unit i are respectively;
constructing coal consumption characteristic constraints of the thermoelectric unit;
Figure RE-GDA00029664988200000616
wherein the content of the first and second substances,
Figure RE-GDA00029664988200000617
is the boiler feed water specific enthalpy of the thermoelectric unit i;
Figure RE-GDA00029664988200000618
is the boiler efficiency of the thermoelectric unit i; HV is the calorific value of the coal;
constructing operation constraint of the phase change heat exchanger;
Figure RE-GDA00029664988200000619
wherein, cwIs the specific heat capacity of water;
Figure RE-GDA00029664988200000620
and
Figure RE-GDA00029664988200000621
the heat exchange quantity and the pipe network water flow of the phase change heat exchanger j in the time period t are respectively;
Figure RE-GDA0002966498820000071
and
Figure RE-GDA0002966498820000072
the temperature of the inlet and the outlet of the pipe network water of the phase change heat exchanger j at the time t is respectively;
Figure RE-GDA0002966498820000073
is the coefficient of performance of the phase change heat exchanger j;
constructing operation constraint of the absorption heat pump;
Figure RE-GDA0002966498820000074
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000075
and
Figure RE-GDA0002966498820000076
the heat produced by the absorption heat pump j in the time period t and the water flow of a pipe network are respectively;
Figure RE-GDA0002966498820000077
and
Figure RE-GDA0002966498820000078
the temperature of the water inlet and outlet of the pipe network of the absorption heat pump j in the time period t respectively;
Figure RE-GDA0002966498820000079
is the coefficient of performance of the absorption heat pump j;
and (3) constructing the operation constraint of the heat supply network in the plant:
Figure RE-GDA00029664988200000710
wherein the content of the first and second substances,
Figure RE-GDA00029664988200000711
and
Figure RE-GDA00029664988200000712
respectively a pipeline number set of heat supply network nodes n in the inflow factory and the outflow factory;
Figure RE-GDA00029664988200000713
is the flow of the thermal pipe p in the plant during the time period t;
Figure RE-GDA00029664988200000714
and
Figure RE-GDA00029664988200000715
the inlet temperature and the outlet temperature of the thermal pipeline p in the plant at a time t are respectively set;
Figure RE-GDA00029664988200000716
is the mixed temperature of the heating network node n in the plant in the time period t.
Further, the thermodynamic system constraint conditions for constructing the electric-thermal comprehensive coordination scheduling model include:
constructing a first-stage heating power pipe network operation constraint;
Figure RE-GDA00029664988200000717
wherein the content of the first and second substances,
Figure RE-GDA00029664988200000718
and
Figure RE-GDA00029664988200000719
the inlet and outlet temperatures of the primary thermal pipeline p in the time period t are respectively; kappap,1、κp,2And kappap,3The coefficients of the temperature dynamic equation of the primary thermal pipeline p are respectively;
Figure RE-GDA00029664988200000720
and
Figure RE-GDA00029664988200000721
the transmission delay and the mass flow of the primary thermal pipeline p are respectively; t issoilIs the soil temperature;
Figure RE-GDA00029664988200000722
is the post-mixing temperature of node n at time period t;
Figure RE-GDA00029664988200000723
and
Figure RE-GDA00029664988200000724
maximum and minimum post-mix temperatures for node n, respectively;
constructing a secondary heating power pipe network operation constraint;
Figure RE-GDA00029664988200000725
wherein p ishIs the number of the first-stage heating power pipeline connected with the heat load h;
Figure RE-GDA0002966498820000081
is the thermal power of the thermal load h over a period t;
Figure RE-GDA0002966498820000082
is the indoor temperature of the thermal load h over a period t;
Figure RE-GDA0002966498820000083
the minimum equivalent thermal resistance of the secondary heating power pipe network h;
building envelope operation constraints;
Figure RE-GDA0002966498820000084
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000085
is the temperature of the middle node of the building envelope h in the time period t;
Figure RE-GDA0002966498820000086
and
Figure RE-GDA0002966498820000087
respectively the heat capacity and the heat resistance of the building envelope h; t isindoor,maxAnd Tindoor,minMaximum and minimum indoor temperatures, respectively;
Figure RE-GDA0002966498820000088
is the outdoor ambient temperature prediction value for time period t.
Further, the constructing of the abstract form of the electric heating comprehensive coordination scheduling model comprises:
Figure RE-GDA0002966498820000089
wherein the decision vector X1The included decision variables are shown as formula (21); the other decision variables are represented as decision vector X2(ii) a The objective function is abstractly represented as F (X)1,X2) (ii) a Constraint equations (6), (7), (10), and (11) are abstractly expressed as X1 ═ G1 (X2); other constraints are abstractly represented as G2(X1, X2) ≦ 0;
Figure RE-GDA00029664988200000810
further, the obtaining of a solution result for scheduling based on the electrothermal synthesis coordination scheduling model comprises:
s1: initializing decision vector X1The value of (A) and the result of (B) are recorded as
Figure RE-GDA00029664988200000811
Setting the convergence error ε>0. The maximum iteration number N and the current iteration number N are 0;
s2: x in abstract form of electric heat comprehensive coordination scheduling model1Is arranged as
Figure RE-GDA00029664988200000812
Solving the remaining decision vector X using a commercial solver2To obtain X2Is expressed as
Figure RE-GDA00029664988200000813
The solver adopts one or more of a CPLEX solver or a GUROBI solver;
s3: according to the above
Figure RE-GDA0002966498820000091
Updating decision vector X1Is expressed as
Figure RE-GDA0002966498820000092
S4: determining
Figure RE-GDA0002966498820000093
Or n>Obtaining a solution result under the condition that N is satisfied
Figure RE-GDA0002966498820000094
Otherwise, let n be incremented by 1 and return to step S2 until a qualified solution is obtained.
In a second aspect, the present invention provides an electric heating comprehensive coordination scheduling device, including:
the parameter acquisition module is used for acquiring model parameters of the power system, model parameters of the thermal power plant and model parameters of the thermodynamic system:
the model construction module is used for constructing an electric heating comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters;
and the solving and scheduling module is used for obtaining a solving result based on the electric heating comprehensive coordination scheduling model so as to perform scheduling.
In a third aspect, the present invention further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of any of the above-mentioned electrothermal comprehensive coordinated scheduling methods when executing the program.
In a fourth aspect, the present invention also provides a non-transitory computer readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the steps of the electrothermal integrated coordinated scheduling method according to any one of the above.
The invention provides an electric heating comprehensive coordination scheduling method, an electric heating comprehensive coordination scheduling device and a medium.
Drawings
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a flow chart of an integrated electric heating and cooling coordination scheduling method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the electric heating comprehensive coordination scheduling device according to the present invention;
fig. 3 is a schematic structural diagram of an electronic device provided in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
An electric heating comprehensive coordination scheduling method of the present invention is described below with reference to fig. 1.
Fig. 1 is a flowchart of an electric heating comprehensive coordination scheduling method according to an embodiment of the present invention.
The embodiment of the invention provides an electric heating comprehensive coordination scheduling method, which comprises the following steps:
step 110: obtaining model parameters of an electric power system, model parameters of a thermal power plant and model parameters of a thermal system:
specifically, electric power system model parameters, thermal power plant model parameters and thermodynamic system model parameters are respectively obtained from operation management departments of the electric power system, the thermal power plant and the thermodynamic system: the power system model parameters comprise a coal consumption characteristic function of a conventional thermal power generating unit, the maximum climbing speed of the conventional thermal power generating unit and the thermoelectric power generating unit, and the maximum and minimum technical output of the conventional thermal power generating unit; the method comprises the following steps of (1) connecting a power bus and a power transmission line, active power transfer distribution factors from the power bus to the power transmission line, the maximum active transmission capacity of the power transmission line, an electric load prediction curve in a scheduling time interval and a wind power active output limit prediction curve;
the model parameters of the thermal power plant comprise topological connection relations among thermal parts, inlet and outlet temperatures, inlet and outlet pressures, steam inlet flow, ideal specific enthalpy drop, internal efficiency, maximum and minimum main steam flow of the thermal power plant, high-pressure cylinder bypass maximum flow, low-pressure cylinder minimum steam inlet flow, sliding pressure curve shape parameters, temperature-reduced water and boiler feed water specific enthalpy, generator efficiency, mechanical efficiency, boiler efficiency, physical property database, specific heat capacity of water, coal heat value, phase-change heat exchanger and performance coefficients of an absorption heat pump of the thermal power plant under rated working conditions;
the thermodynamic system model parameters comprise a node-pipeline connection relation of a primary thermodynamic pipeline, a temperature dynamic equation coefficient, transmission delay, mass flow, maximum and minimum node temperatures, soil temperature, minimum equivalent thermal resistance of a secondary thermodynamic pipeline network, maximum and minimum indoor temperatures, thermal capacity and thermal resistance of a building enclosure structure and an outdoor environment temperature prediction curve.
Step 120: constructing an electric heating comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters, wherein the electric heating comprehensive coordination scheduling model is a model of electric heating;
specifically, a series of objective functions, power system constraints, thermal power plant constraints, thermodynamic system constraints, and the like may be constructed. In order to realize that an electric heating comprehensive coordination scheduling model is constructed based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters, an objective function of the electric heating comprehensive coordination scheduling model can be constructed, the expression is as follows:
Figure RE-GDA0002966498820000111
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000112
the active power output of the conventional thermal power generating unit i in the time period t is shown;
Figure RE-GDA0002966498820000113
the coal consumption characteristic function of a conventional thermal power generating unit i is obtained;
Figure RE-GDA0002966498820000114
the coal consumption of the thermoelectric unit i in the time period t is shown; the objective function shown in the formula (1) is to minimize the total operation coal consumption of the electric heating integrated energy system, including the operation coal consumption of a conventional thermal power generating unit and the operation coal consumption of the thermal power generating unit.
Electric power system constraint conditions of the electric heating comprehensive coordination scheduling model can be established; the thermal power plant constraint conditions of the electric heating comprehensive coordination scheduling model can be established; the thermal power plant comprises one or more of thermal parts or thermal links such as a thermal motor set, a phase-change heat exchanger, an absorption heat pump, an in-plant heat supply network and the like; the thermoelectric unit comprises one or more types of back pressure units or extraction condensing units; a thermodynamic system constraint condition of the electric heating comprehensive coordination scheduling model can be established; an abstract form of the electric heating comprehensive coordination scheduling model can be constructed.
Further, the power system constraints include: the method comprises the following steps of (1) power system energy balance constraint, unit active output upper and lower limit constraint, unit climbing rate constraint and line active transmission capacity constraint;
the power system energy balance constraint is as follows:
Figure RE-GDA0002966498820000121
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000122
the active power of the thermoelectric unit i in the time period t;
Figure RE-GDA0002966498820000123
is the active power output of the wind farm w at time t;
Figure RE-GDA0002966498820000124
is a predicted value of the electrical load e at the time period t; in any scheduling period, the power system needs to meet energy balance;
the active output upper and lower limits of the unit are constrained as follows:
Figure RE-GDA0002966498820000125
wherein, the active power of the conventional thermal power plant and the wind power plant needs to be in a certain range,
Figure RE-GDA0002966498820000126
and
Figure RE-GDA0002966498820000127
respectively a conventional thermal engineMaximum, minimum technical output for group i;
Figure RE-GDA0002966498820000128
the active output limit predicted value of the wind power plant w in the time period t is obtained;
the unit climbing rate constraint is as follows:
Figure RE-GDA0002966498820000129
wherein the content of the first and second substances,
Figure RE-GDA00029664988200001210
and
Figure RE-GDA00029664988200001211
the maximum climbing speed and the maximum climbing speed of the conventional thermal power generating unit i are respectively set;
Figure RE-GDA00029664988200001212
and
Figure RE-GDA00029664988200001213
the maximum climbing speed and the maximum climbing speed of the thermoelectric unit i are respectively set; Δ t is the length of the scheduling period; here, the active regulation rate of the conventional thermal power generating unit and the thermoelectric power generating unit should not exceed a certain threshold value.
The active transmission capacity constraint of the line is as follows:
Figure RE-GDA0002966498820000131
wherein, SFulIs the active power transfer distribution factor from the power bus u to the power transmission line l;
Figure RE-GDA0002966498820000132
is the maximum active transmission capacity of the power transmission line l.
Further, the thermal power plant constraint conditions for constructing the electric-thermal comprehensive coordination scheduling model comprise:
constructing operation constraints of high, medium and low pressure cylinders of the thermoelectric unit; for the high pressure cylinder, the middle pressure cylinder and the low pressure cylinder of the back pressure unit and the high pressure cylinder and the middle pressure cylinder of the extraction condensing unit, the following operation constraints are constructed:
Figure RE-GDA0002966498820000133
the formula (6) sequentially gives a through-flow equation and a thermal power conversion equation in a high-pressure cylinder, a middle-pressure cylinder and a low-pressure cylinder of the back pressure unit and a high-pressure cylinder and a middle-pressure cylinder of the pumping and condensing unit, and at the moment, the final-stage steam of each cylinder is usually in an overheat state; wherein, the superscript x can be HP, IP and LP, which respectively represent the relevant parameters of high, middle and low pressure cylinders; the superscript 0 represents the relevant parameter under the rated working condition;
Figure RE-GDA0002966498820000134
the steam inlet flow of the high, medium and low pressure cylinders;
Figure RE-GDA0002966498820000135
and
Figure RE-GDA0002966498820000136
the inlet pressure and the outlet pressure of the high-pressure cylinder, the medium-pressure cylinder and the low-pressure cylinder are respectively set;
Figure RE-GDA0002966498820000137
and
Figure RE-GDA0002966498820000138
the inlet and outlet temperatures of the high, medium and low pressure cylinders are respectively;
Figure RE-GDA0002966498820000139
the ideal temperature of the steam after isentropic adiabatic expansion in the high, medium and low pressure cylinders;
Figure RE-GDA00029664988200001310
and
Figure RE-GDA00029664988200001311
the actual specific enthalpy drop in the high-pressure cylinder, the intermediate-pressure cylinder and the low-pressure cylinder and the ideal specific enthalpy drop during isentropic adiabatic expansion are respectively;
Figure RE-GDA00029664988200001312
the internal efficiency of the high, medium and low pressure cylinders under rated working conditions; h (-) and s (-) are a specific enthalpy function and a specific entropy function, respectively, and the function values can be queried through a physical property database; the physical property database is an existing database;
for the low pressure cylinder of the extraction and condensation unit, the following operational constraints are established:
Figure RE-GDA0002966498820000141
because the steam of the last stage of the low pressure cylinder of the extraction and condensation type unit is usually in a wet saturation state, the through-flow equation and the heat-work conversion equation of the low pressure cylinder of the extraction and condensation type unit are independently given by the formula (7); wherein the content of the first and second substances,
Figure RE-GDA0002966498820000142
and
Figure RE-GDA0002966498820000143
the actual dryness of the last-stage steam of the low-pressure cylinder and the ideal dryness of the steam of the low-pressure cylinder after isentropic adiabatic expansion are respectively obtained;
constructing coupling constraints among high, medium and low pressure cylinders of the thermoelectric unit:
for the back pressure unit, the following steam flow balance equation is constructed:
Figure RE-GDA0002966498820000144
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000145
the serial number of the phase change heat exchanger connected with the thermoelectric unit i;
Figure RE-GDA0002966498820000146
and
Figure RE-GDA0002966498820000147
respectively the main steam flow and the regenerative steam flow of the thermoelectric unit i in a time period t;
Figure RE-GDA0002966498820000148
is the heating steam flow of the phase change heat exchanger j in the time period t; alpha is alphai、βiAnd gammaiThe final steam flow of the high, medium and low pressure cylinders of the thermoelectric unit i accounts for the proportion of the steam inlet flow;
Figure RE-GDA0002966498820000149
and
Figure RE-GDA00029664988200001410
the maximum main steam flow and the minimum main steam flow of the thermoelectric unit i are respectively;
for the extraction condensing unit, the following steam flow balance equation is constructed:
Figure RE-GDA00029664988200001411
Figure RE-GDA0002966498820000151
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000152
is the serial number of the absorption heat pump connected with the thermoelectric unit i;
Figure RE-GDA0002966498820000153
and
Figure RE-GDA0002966498820000154
is the bypass flow of the high-pressure cylinder and the steam extraction flow of the medium-pressure cylinder of the thermoelectric unit i in the time period tAn amount;
Figure RE-GDA0002966498820000155
and
Figure RE-GDA0002966498820000156
the maximum flow of a high-pressure cylinder bypass and the minimum steam inlet flow of a low-pressure cylinder of the thermoelectric unit i are respectively;
Figure RE-GDA0002966498820000157
and
Figure RE-GDA0002966498820000158
the proportions of the high-pressure cylinder bypass desuperheating water and the medium-pressure cylinder bypass desuperheating water of the thermoelectric generator set i at the time t are respectively;
Figure RE-GDA0002966498820000159
is the specific enthalpy of the reduced temperature water of the thermoelectric unit i;
Figure RE-GDA00029664988200001510
is the driving steam flow of the absorption heat pump j in the time period t;
and for the back pressure unit and the extraction condensing unit, constructing a coupling relation of inlet and outlet temperatures and pressures of high, medium and low pressure cylinders:
Figure RE-GDA00029664988200001511
wherein the content of the first and second substances,
Figure RE-GDA00029664988200001512
the steam inlet temperature of the high and medium pressure cylinders under rated working conditions;
Figure RE-GDA00029664988200001513
aiand biThe shape parameters of a sliding pressure curve of the thermoelectric unit i are all the shape parameters;
and constructing an external characteristic equation for the back pressure unit and the extraction condensing unit:
Figure RE-GDA00029664988200001514
wherein the content of the first and second substances,
Figure RE-GDA00029664988200001515
and
Figure RE-GDA00029664988200001516
the generator efficiency and the mechanical efficiency of the thermoelectric unit i are respectively; in the formula (12), the power output of the thermoelectric unit is in direct proportion to the sum of the specific enthalpy drops in the high, medium and low pressure cylinders;
constructing coal consumption characteristic constraints of the thermoelectric unit;
Figure RE-GDA00029664988200001517
wherein the content of the first and second substances,
Figure RE-GDA00029664988200001518
is the boiler feed water specific enthalpy of the thermoelectric unit i;
Figure RE-GDA00029664988200001519
is the boiler efficiency of the thermoelectric unit i; HV is the calorific value of the coal; here, the operation coal consumption of the thermoelectric power generation unit includes coal consumption in the main steam heating process and coal consumption in the reheating process.
Constructing operation constraint of the phase change heat exchanger;
Figure RE-GDA0002966498820000161
wherein, cwIs the specific heat capacity of water;
Figure RE-GDA0002966498820000162
and
Figure RE-GDA0002966498820000163
the heat exchange quantity and the pipe network of the phase change heat exchanger j in the time period t respectivelyWater flow rate;
Figure RE-GDA0002966498820000164
and
Figure RE-GDA0002966498820000165
the temperature of the inlet and the outlet of the pipe network water of the phase change heat exchanger j at the time t is respectively;
Figure RE-GDA0002966498820000166
is the coefficient of performance of the phase change heat exchanger j;
constructing operation constraint of the absorption heat pump;
Figure RE-GDA0002966498820000167
wherein the content of the first and second substances,
Figure RE-GDA0002966498820000168
and
Figure RE-GDA0002966498820000169
the heat produced by the absorption heat pump j in the time period t and the water flow of a pipe network are respectively;
Figure RE-GDA00029664988200001610
and
Figure RE-GDA00029664988200001611
the temperature of the water inlet and outlet of the pipe network of the absorption heat pump j in the time period t respectively;
Figure RE-GDA00029664988200001612
is the coefficient of performance of the absorption heat pump j;
and (3) constructing the operation constraint of the heat supply network in the plant:
Figure RE-GDA00029664988200001613
wherein the content of the first and second substances,
Figure RE-GDA00029664988200001614
and
Figure RE-GDA00029664988200001615
respectively a pipeline number set of heat supply network nodes n in the inflow factory and the outflow factory;
Figure RE-GDA00029664988200001616
is the flow of the thermal pipe p in the plant during the time period t;
Figure RE-GDA00029664988200001617
and
Figure RE-GDA00029664988200001618
the inlet temperature and the outlet temperature of the thermal pipeline p in the plant at a time t are respectively set;
Figure RE-GDA00029664988200001619
is the mixed temperature of the heating network node n in the plant in the time period t. Here, equation (16) gives the flow balance equation at the heating network node in the plant and the energy balance equation.
Further, the thermodynamic system constraint conditions for constructing the electric-thermal comprehensive coordination scheduling model comprise:
constructing a first-stage heating power pipe network operation constraint;
Figure RE-GDA00029664988200001620
wherein the content of the first and second substances,
Figure RE-GDA00029664988200001621
and
Figure RE-GDA00029664988200001622
the inlet and outlet temperatures of the primary thermal pipeline p in the time period t are respectively; kappap,1、κp,2And kappap,3The coefficients of the temperature dynamic equation of the primary thermal pipeline p are respectively;
Figure RE-GDA00029664988200001623
and
Figure RE-GDA0002966498820000171
the transmission delay and the mass flow of the primary thermal pipeline p are respectively; t issoilIs the soil temperature;
Figure RE-GDA0002966498820000172
is the post-mixing temperature of node n at time period t;
Figure RE-GDA0002966498820000173
and
Figure RE-GDA0002966498820000174
maximum and minimum post-mix temperatures for node n, respectively; the formula (17) sequentially provides a temperature dynamic equation of the primary heat distribution pipeline, an energy conservation equation at a node in the primary heat distribution pipeline network, and upper and lower temperature limit constraints at the node in the primary heat distribution pipeline network.
Constructing a secondary heating power pipe network operation constraint;
Figure RE-GDA0002966498820000175
wherein p ishIs the number of the first-stage heating power pipeline connected with the heat load h;
Figure RE-GDA0002966498820000176
is the thermal power of the thermal load h over a period t;
Figure RE-GDA0002966498820000177
is the indoor temperature of the thermal load h over a period t;
Figure RE-GDA0002966498820000178
the minimum equivalent thermal resistance of the secondary heating power pipe network h;
building envelope operation constraints;
Figure RE-GDA0002966498820000179
wherein equation (19) first gives the discretized building envelope heat transfer equation, where
Figure RE-GDA00029664988200001710
Is the temperature of the middle node of the building envelope h in the time period t;
Figure RE-GDA00029664988200001711
and
Figure RE-GDA00029664988200001712
respectively the heat capacity and the heat resistance of the building envelope h; t isindoor,maxAnd Tindoor,minMaximum and minimum indoor temperatures, respectively;
Figure RE-GDA00029664988200001713
is the predicted outdoor ambient temperature value for time period t, and meanwhile, the indoor temperature needs to be within a certain range in order to ensure the heat supply comfort of the heat user.
It is worth to be noted that the abstract form for constructing the electric heating comprehensive coordination scheduling model includes:
Figure RE-GDA00029664988200001714
wherein the decision vector X1The included decision variables are shown as formula (21); the other decision variables are represented as decision vector X2(ii) a The objective function is abstractly represented as F (X)1,X2) (ii) a Constraint equations (6), (7), (10), and (11) are abstractly expressed as X1 ═ G1 (X2); other constraints are abstractly represented as G2(X1, X2) ≦ 0;
Figure RE-GDA0002966498820000181
step 130: and obtaining a solving result based on the electric heating comprehensive coordination scheduling model so as to carry out scheduling.
After the electric heating comprehensive coordination scheduling model is obtained, solving the electric heating comprehensive coordination scheduling model considering the internal structure of the thermal power plant, wherein the thermal power plant constraint equation in the electric heating comprehensive coordination scheduling model considering the internal structure of the thermal power plant is nonlinear and cannot be directly solved by adopting a commercial solver, so that the invention starts from the special structure of the formula (20) and adopts an iteration method to solve the problem, and the method comprises the following steps:
s1: initializing decision vector X1The value of (A) and the result of (B) are recorded as
Figure RE-GDA0002966498820000182
Setting the convergence error ε>0. The maximum iteration number N and the current iteration number N are 0; in this embodiment, let ε be 10-3, N=20;
S2: x in abstract form of electric heat comprehensive coordination scheduling model1Is arranged as
Figure RE-GDA0002966498820000183
Solving the remaining decision vector X using a commercial solver2To obtain X2Is expressed as
Figure RE-GDA0002966498820000184
The solver adopts one or more of a CPLEX solver or a GUROBI solver; here, the CPLEX solver or the GUROBI solver can efficiently process various optimization types such as linear programming, quadratic programming, mixed integer linear programming, mixed integer quadratic programming, second order cone programming, and the like;
s3: according to the above
Figure RE-GDA0002966498820000185
Updating decision vector X1Is expressed as
Figure RE-GDA0002966498820000186
S4: determining
Figure RE-GDA0002966498820000187
Or n>Obtaining a solution result under the condition that N is satisfied
Figure RE-GDA0002966498820000188
Otherwise, n is incremented by 1(n ═ n +1), and the process returns to step S2 until a satisfactory solution result is obtained.
It can be seen that the invention converts the electric heating comprehensive coordination scheduling problem which cannot be solved directly by a commercial solver into a series of problems which can be processed by the commercial solver, thereby facilitating the use of the invention in practical engineering.
The electric heating comprehensive coordination scheduling device provided by the invention has the advantages that the detailed thermal power plant model is introduced into the electric heating comprehensive coordination scheduling problem, the physical process in the thermal power plant can be more accurately described, the real physical law of electric heating system coupling is reflected, the scheduling is more reasonably carried out, meanwhile, in order to process the problem that the model is difficult to solve caused by the introduction of the detailed thermal power plant model, the practical iterative solution method is provided, the electric heating comprehensive coordination scheduling model can be converted into a series of models which can be processed by a commercial solver, and the application in practical engineering is facilitated.
The electric heating comprehensive coordination scheduling device provided by the invention is described below, and the electric heating comprehensive coordination scheduling device described below and the electric heating comprehensive coordination scheduling method described above can be referred to correspondingly.
Referring to fig. 2, fig. 2 is a schematic structural diagram of an electric-heating comprehensive coordination scheduling device according to the present invention.
The invention provides an electric heating comprehensive coordination scheduling device 200, comprising:
a parameter obtaining module 210, configured to obtain power system model parameters, thermal power plant model parameters, and thermodynamic system model parameters:
a model construction module 220, configured to construct an electric-thermal comprehensive coordination scheduling model based on the power system model parameters, the thermal power plant model parameters, and the thermodynamic system model parameters;
and the solution scheduling module 230 is used for obtaining a solution result based on the electrothermal comprehensive coordination scheduling model so as to perform scheduling.
The parameter acquisition module is configured to: obtaining model parameters from operation management departments of an electric power system, a thermal power plant and a thermodynamic system and sending the model parameters to a model construction module; the model parameters comprise electric power system model parameters, thermal power plant model parameters and thermodynamic system model parameters;
the model building module is configured to: according to the model parameters, sequentially constructing an objective function, an electric power system constraint condition, a thermal power plant constraint condition and a thermodynamic system constraint condition of the electric heating comprehensive coordination scheduling model, constructing an abstract form of the electric heating comprehensive coordination scheduling model and sending the abstract form to a solving scheduling module;
the solution scheduling module is configured to: and solving the electric heating comprehensive coordination scheduling model by adopting the method in the embodiment of the method, and outputting a solving result.
As can be seen, the invention introduces a detailed thermal power plant model in the electric heating comprehensive coordination scheduling problem, can more accurately describe the physical process in the thermal power plant, and reflects the real physical law of electric heating system coupling. Meanwhile, in order to solve the problem that a model is difficult to solve due to the introduction of a detailed thermal power plant model, the invention provides a practical iterative solution method, and the method can convert an electric heating comprehensive coordination scheduling model into a series of models which can be processed by a commercial solver, so that the application in practical engineering is facilitated.
The electric heating comprehensive coordination scheduling device provided by the invention has the advantages that the detailed thermal power plant model is introduced into the electric heating comprehensive coordination scheduling problem, the physical process in the thermal power plant can be more accurately described, the real physical law of electric heating system coupling is reflected, the scheduling is more reasonably carried out, meanwhile, in order to process the problem that the model is difficult to solve caused by the introduction of the detailed thermal power plant model, the practical iterative solution method is provided, the electric heating comprehensive coordination scheduling model can be converted into a series of models which can be processed by a commercial solver, and the application in practical engineering is facilitated.
Fig. 3 illustrates a physical structure diagram of an electronic device, which may include, as shown in fig. 3: a processor (processor)310, a communication Interface (communication Interface)320, a memory (memory)330 and a communication bus 340, wherein the processor 310, the communication Interface 320 and the memory 330 communicate with each other via the communication bus 340. The processor 310 may invoke logic instructions in the memory 330 to perform an electro-thermal integrated coordinated scheduling method comprising: obtaining model parameters of an electric power system, model parameters of a thermal power plant and model parameters of a thermal system: constructing an electric heating comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters, wherein the electric heating comprehensive coordination scheduling model is a model of electric heating; and obtaining a solving result based on the electric heating comprehensive coordination scheduling model so as to carry out scheduling.
In addition, the logic instructions in the memory 330 may be implemented in the form of software functional units and stored in a computer readable storage medium when the software functional units are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, the present invention also provides a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions, which when executed by a computer, enable the computer to perform the electric-thermal integrated coordination scheduling method provided by the above methods, the method comprising: obtaining model parameters of an electric power system, model parameters of a thermal power plant and model parameters of a thermal system: constructing an electric heating comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters, wherein the electric heating comprehensive coordination scheduling model is a model of electric heating; and obtaining a solving result based on the electric heating comprehensive coordination scheduling model so as to carry out scheduling.
In yet another aspect, the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored, the computer program being implemented by a processor to perform the electrothermal integrated coordination scheduling methods provided above, the method comprising: obtaining model parameters of an electric power system, model parameters of a thermal power plant and model parameters of a thermal system: constructing an electric heating comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters, wherein the electric heating comprehensive coordination scheduling model is a model of electric heating; and obtaining a solving result based on the electric heating comprehensive coordination scheduling model so as to carry out scheduling.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (10)

1. An electric heating comprehensive coordination scheduling method is characterized by comprising the following steps:
obtaining model parameters of an electric power system, model parameters of a thermal power plant and model parameters of a thermal system:
constructing an electric heating comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters, wherein the electric heating comprehensive coordination scheduling model is a model of electric heating;
and obtaining a solving result based on the electric heating comprehensive coordination scheduling model so as to carry out scheduling.
2. The electric-thermal integrated coordinated scheduling method of claim 1,
the power system model parameters comprise a coal consumption characteristic function of a conventional thermal power generating unit, the maximum climbing speed of the conventional thermal power generating unit and the thermoelectric power generating unit, and the maximum and minimum technical output of the conventional thermal power generating unit; the method comprises the following steps of (1) connecting a power bus and a power transmission line, active power transfer distribution factors from the power bus to the power transmission line, the maximum active transmission capacity of the power transmission line, an electric load prediction curve in a scheduling time interval and a wind power active output limit prediction curve;
the model parameters of the thermal power plant comprise topological connection relations among thermal parts, inlet and outlet temperatures, inlet and outlet pressures, steam inlet flow, ideal specific enthalpy drop, internal efficiency, maximum and minimum main steam flow of the thermal power plant, high-pressure cylinder bypass maximum flow, low-pressure cylinder minimum steam inlet flow, sliding pressure curve shape parameters, temperature-reduced water and boiler feed water specific enthalpy, generator efficiency, mechanical efficiency, boiler efficiency, physical property database, specific heat capacity of water, coal heat value, phase-change heat exchanger and performance coefficients of an absorption heat pump of the thermal power plant under rated working conditions;
the thermodynamic system model parameters comprise a node-pipeline connection relation of a primary thermodynamic pipeline, a temperature dynamic equation coefficient, transmission delay, mass flow, maximum and minimum node temperatures, soil temperature, minimum equivalent thermal resistance of a secondary thermodynamic pipeline network, maximum and minimum indoor temperatures, thermal capacity and thermal resistance of a building enclosure structure and an outdoor environment temperature prediction curve.
3. The electric-thermal integrated coordinated scheduling method of claim 1,
the step of constructing an electric-heat comprehensive coordination scheduling model based on the electric power system model parameters, the thermal power plant model parameters and the thermodynamic system model parameters comprises the following steps:
constructing an objective function of an electric heating comprehensive coordination scheduling model, wherein the expression is as follows:
Figure FDA0002802877060000011
wherein the content of the first and second substances,
Figure FDA0002802877060000021
the active power output of the conventional thermal power generating unit i in the time period t is shown;
Figure FDA0002802877060000022
the coal consumption characteristic function of a conventional thermal power generating unit i is obtained;
Figure FDA0002802877060000023
the coal consumption of the thermoelectric unit i in the time period t is shown;
constructing a power system constraint condition of an electric heating comprehensive coordination scheduling model;
constructing thermal power plant constraint conditions of an electric heating comprehensive coordination scheduling model; the thermal power plant comprises one or more of thermal parts or thermal links such as a thermal motor set, a phase-change heat exchanger, an absorption heat pump, an in-plant heat supply network and the like; the thermoelectric unit comprises one or more types of back pressure units or extraction condensing units;
constructing a thermodynamic system constraint condition of an electric heating comprehensive coordination scheduling model;
and constructing an abstract form of the electric heating comprehensive coordination scheduling model.
4. The electric-thermal integrated coordinated scheduling method of claim 3,
the power system constraints include: the method comprises the following steps of (1) power system energy balance constraint, unit active output upper and lower limit constraint, unit climbing rate constraint and line active transmission capacity constraint;
the power system energy balance constraint is as follows:
Figure FDA0002802877060000024
wherein the content of the first and second substances,
Figure FDA0002802877060000025
the active power of the thermoelectric unit i in the time period t;
Figure FDA0002802877060000026
is the active power output of the wind farm w at time t;
Figure FDA0002802877060000027
is a predicted value of the electrical load e at the time period t;
the active output upper and lower limits of the unit are constrained as follows:
Figure FDA0002802877060000028
wherein, Pi TU,maxAnd Pi TU,minRespectively the maximum and minimum technical output of a conventional thermal power generating unit i;
Figure FDA0002802877060000029
the active output limit predicted value of the wind power plant w in the time period t is obtained;
the unit climbing rate constraint is as follows:
Figure FDA00028028770600000210
wherein the content of the first and second substances,
Figure FDA00028028770600000211
and
Figure FDA00028028770600000212
the maximum climbing speed and the maximum climbing speed of the conventional thermal power generating unit i are respectively set;
Figure FDA00028028770600000213
and
Figure FDA00028028770600000214
the maximum climbing speed and the maximum climbing speed of the thermoelectric unit i are respectively set; Δ t is the length of the scheduling period;
the active transmission capacity constraint of the line is as follows:
Figure FDA0002802877060000031
wherein, SFulIs the active power transfer distribution factor from the power bus u to the power transmission line l; pl L,maxIs electric powerThe maximum active transmission capacity of the transmission line l.
5. The electric-thermal integrated coordinated scheduling method of claim 3,
the thermal power plant constraint conditions for constructing the electric heating comprehensive coordination scheduling model comprise:
constructing operation constraints of high, medium and low pressure cylinders of the thermoelectric unit; for the high pressure cylinder, the middle pressure cylinder and the low pressure cylinder of the back pressure unit and the high pressure cylinder and the middle pressure cylinder of the extraction condensing unit, the following operation constraints are constructed:
Figure FDA0002802877060000032
wherein, the superscript x can be HP, IP and LP, which respectively represent the relevant parameters of high, middle and low pressure cylinders; the superscript 0 represents the relevant parameter under the rated working condition;
Figure FDA0002802877060000033
the steam inlet flow of the high, medium and low pressure cylinders;
Figure FDA0002802877060000034
and
Figure FDA0002802877060000035
the inlet pressure and the outlet pressure of the high-pressure cylinder, the medium-pressure cylinder and the low-pressure cylinder are respectively set;
Figure FDA0002802877060000036
and
Figure FDA0002802877060000037
the inlet and outlet temperatures of the high, medium and low pressure cylinders are respectively;
Figure FDA0002802877060000038
the ideal temperature of the steam after isentropic adiabatic expansion in the high, medium and low pressure cylinders;
Figure FDA0002802877060000039
and
Figure FDA00028028770600000310
the actual specific enthalpy drop in the high-pressure cylinder, the intermediate-pressure cylinder and the low-pressure cylinder and the ideal specific enthalpy drop during isentropic adiabatic expansion are respectively;
Figure FDA00028028770600000311
the internal efficiency of the high, medium and low pressure cylinders under rated working conditions; h (-) and s (-) are a specific enthalpy function and a specific entropy function, respectively, and the function values can be queried through a physical property database; the physical property database is an existing database;
for the low pressure cylinder of the extraction and condensation unit, the following operational constraints are established:
Figure FDA0002802877060000041
wherein the content of the first and second substances,
Figure FDA0002802877060000042
and
Figure FDA0002802877060000043
the actual dryness of the last-stage steam of the low-pressure cylinder and the ideal dryness of the steam of the low-pressure cylinder after isentropic adiabatic expansion are respectively obtained;
constructing coupling constraints among high, medium and low pressure cylinders of the thermoelectric unit:
for the back pressure unit, the following steam flow balance equation is constructed:
Figure FDA0002802877060000044
wherein the content of the first and second substances,
Figure FDA0002802877060000045
the serial number of the phase change heat exchanger connected with the thermoelectric unit i;
Figure FDA0002802877060000046
and
Figure FDA0002802877060000047
respectively the main steam flow and the regenerative steam flow of the thermoelectric unit i in a time period t;
Figure FDA0002802877060000048
is the heating steam flow of the phase change heat exchanger j in the time period t; alpha is alphai、βiAnd gammaiThe final steam flow of the high, medium and low pressure cylinders of the thermoelectric unit i accounts for the proportion of the steam inlet flow;
Figure FDA0002802877060000049
and
Figure FDA00028028770600000410
the maximum main steam flow and the minimum main steam flow of the thermoelectric unit i are respectively;
for the extraction condensing unit, the following steam flow balance equation is constructed:
Figure FDA00028028770600000411
Figure FDA00028028770600000412
wherein the content of the first and second substances,
Figure FDA00028028770600000413
is the serial number of the absorption heat pump connected with the thermoelectric unit i;
Figure FDA00028028770600000414
and
Figure FDA00028028770600000415
the bypass flow of a high-pressure cylinder and the steam extraction flow of a medium-pressure cylinder of the thermoelectric unit i in a time period t are shown;
Figure FDA00028028770600000416
and
Figure FDA0002802877060000051
the maximum flow of a high-pressure cylinder bypass and the minimum steam inlet flow of a low-pressure cylinder of the thermoelectric unit i are respectively;
Figure FDA0002802877060000052
and
Figure FDA0002802877060000053
the proportions of the high-pressure cylinder bypass desuperheating water and the medium-pressure cylinder bypass desuperheating water of the thermoelectric generator set i at the time t are respectively;
Figure FDA0002802877060000054
is the specific enthalpy of the reduced temperature water of the thermoelectric unit i;
Figure FDA0002802877060000055
is the driving steam flow of the absorption heat pump j in the time period t;
and for the back pressure unit and the extraction condensing unit, constructing a coupling relation of inlet and outlet temperatures and pressures of high, medium and low pressure cylinders:
Figure FDA0002802877060000056
wherein, Ti in,0The steam inlet temperature of the high and medium pressure cylinders under rated working conditions;
Figure FDA0002802877060000057
aiand biThe shape parameters of a sliding pressure curve of the thermoelectric unit i are all the shape parameters;
and constructing an external characteristic equation for the back pressure unit and the extraction condensing unit:
Figure FDA0002802877060000058
wherein the content of the first and second substances,
Figure FDA0002802877060000059
and
Figure FDA00028028770600000510
the generator efficiency and the mechanical efficiency of the thermoelectric unit i are respectively;
constructing coal consumption characteristic constraints of the thermoelectric unit;
Figure FDA00028028770600000511
wherein the content of the first and second substances,
Figure FDA00028028770600000512
is the boiler feed water specific enthalpy of the thermoelectric unit i;
Figure FDA00028028770600000513
is the boiler efficiency of the thermoelectric unit i; HV is the calorific value of the coal;
constructing operation constraint of the phase change heat exchanger;
Figure FDA00028028770600000514
wherein, cwIs the specific heat capacity of water;
Figure FDA00028028770600000515
and
Figure FDA00028028770600000516
the heat exchange quantity and the pipe network water flow of the phase change heat exchanger j in the time period t are respectively;
Figure FDA00028028770600000517
and
Figure FDA00028028770600000518
the temperature of the inlet and the outlet of the pipe network water of the phase change heat exchanger j at the time t is respectively;
Figure FDA00028028770600000519
is the coefficient of performance of the phase change heat exchanger j;
constructing operation constraint of the absorption heat pump;
Figure FDA0002802877060000061
wherein the content of the first and second substances,
Figure FDA0002802877060000062
and
Figure FDA0002802877060000063
the heat produced by the absorption heat pump j in the time period t and the water flow of a pipe network are respectively;
Figure FDA0002802877060000064
and
Figure FDA0002802877060000065
the temperature of the water inlet and outlet of the pipe network of the absorption heat pump j in the time period t respectively;
Figure FDA0002802877060000066
is the coefficient of performance of the absorption heat pump j;
and (3) constructing the operation constraint of the heat supply network in the plant:
Figure FDA0002802877060000067
wherein the content of the first and second substances,
Figure FDA0002802877060000068
and
Figure FDA0002802877060000069
respectively a pipeline number set of heat supply network nodes n in the inflow factory and the outflow factory;
Figure FDA00028028770600000610
is the flow of the thermal pipe p in the plant during the time period t;
Figure FDA00028028770600000611
and
Figure FDA00028028770600000612
the inlet temperature and the outlet temperature of the thermal pipeline p in the plant at a time t are respectively set;
Figure FDA00028028770600000613
is the mixed temperature of the heating network node n in the plant in the time period t.
6. The electric-thermal integrated coordinated scheduling method of claim 3,
the thermodynamic system constraint conditions for constructing the electric heating comprehensive coordination scheduling model comprise:
constructing a first-stage heating power pipe network operation constraint;
Figure FDA00028028770600000614
wherein the content of the first and second substances,
Figure FDA00028028770600000615
and
Figure FDA00028028770600000616
the inlet and outlet temperatures of the primary thermal pipeline p in the time period t are respectively; kappap,1、κp,2And kappap,3The coefficients of the temperature dynamic equation of the primary thermal pipeline p are respectively;
Figure FDA00028028770600000617
and
Figure FDA00028028770600000618
the transmission delay and the mass flow of the primary thermal pipeline p are respectively; t issoilIs the soil temperature;
Figure FDA00028028770600000619
is the post-mixing temperature of node n at time period t;
Figure FDA00028028770600000620
and
Figure FDA00028028770600000621
maximum and minimum post-mix temperatures for node n, respectively;
constructing a secondary heating power pipe network operation constraint;
Figure FDA00028028770600000622
wherein p ishIs the number of the first-stage heating power pipeline connected with the heat load h;
Figure FDA00028028770600000623
is the thermal power of the thermal load h over a period t;
Figure FDA00028028770600000624
is the indoor temperature of the thermal load h over a period t;
Figure FDA00028028770600000625
the minimum equivalent thermal resistance of the secondary heating power pipe network h;
building envelope operation constraints;
Figure FDA0002802877060000071
wherein the content of the first and second substances,
Figure FDA0002802877060000072
is the temperature of the middle node of the building envelope h in the time period t;
Figure FDA0002802877060000073
and
Figure FDA0002802877060000074
respectively the heat capacity and the heat resistance of the building envelope h; t isindoor,maxAnd Tindoor,minMaximum and minimum indoor temperatures, respectively; t ist envIs the outdoor ambient temperature prediction value for time period t.
7. The electric-thermal integrated coordinated scheduling method of claim 3,
the abstract form for constructing the electric heating comprehensive coordination scheduling model comprises the following steps:
Figure FDA0002802877060000075
wherein the decision vector X1The included decision variables are shown as formula (21); the other decision variables are represented as decision vector X2(ii) a The objective function is abstractly represented as F (X)1,X2) (ii) a Constraint equations (6), (7), (10), and (11) are abstractly expressed as X1 ═ G1 (X2); other constraints are abstractly represented as G2(X1, X2) ≦ 0;
Figure FDA0002802877060000076
8. the electric-thermal integrated coordinated scheduling method of claim 7,
the obtaining of a solution result for scheduling based on the electrothermal integrated coordination scheduling model comprises:
s1: initializing decision vector X1The value of (A) and the result of (B) are recorded as
Figure FDA0002802877060000077
Setting the convergence error ε>0. The maximum iteration number N and the current iteration number N are 0;
s2: x in abstract form of electric heat comprehensive coordination scheduling model1Is arranged as
Figure FDA0002802877060000078
Solving the remaining decision vector X using a commercial solver2To obtain X2Is expressed as
Figure FDA0002802877060000079
The solver adopts one or more of a CPLEX solver or a GUROBI solver;
s3: according to the above
Figure FDA0002802877060000081
Updating decision vector X1Is expressed as
Figure FDA0002802877060000082
S4: determining
Figure FDA0002802877060000083
Or n>Obtaining a solution result under the condition that N is satisfied
Figure FDA0002802877060000084
Otherwise, let n be incremented by 1 and return to step S2 until a qualified solution is obtained.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the electrothermal integrated coordinated scheduling method according to any one of claims 1 to 7 when executing the program.
10. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the electro-thermal integrated coordinated scheduling method according to any one of claims 1 to 7.
CN202011357166.0A 2020-11-26 2020-11-26 Electric heating comprehensive coordination scheduling method and device, equipment and medium Pending CN112734591A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011357166.0A CN112734591A (en) 2020-11-26 2020-11-26 Electric heating comprehensive coordination scheduling method and device, equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011357166.0A CN112734591A (en) 2020-11-26 2020-11-26 Electric heating comprehensive coordination scheduling method and device, equipment and medium

Publications (1)

Publication Number Publication Date
CN112734591A true CN112734591A (en) 2021-04-30

Family

ID=75597843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011357166.0A Pending CN112734591A (en) 2020-11-26 2020-11-26 Electric heating comprehensive coordination scheduling method and device, equipment and medium

Country Status (1)

Country Link
CN (1) CN112734591A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113902040A (en) * 2021-11-15 2022-01-07 中国电力科学研究院有限公司 Method, system, equipment and storage medium for coordinating and optimizing electricity-heat comprehensive energy system
CN114370896A (en) * 2021-12-29 2022-04-19 贵州电网有限责任公司 Method for monitoring heating power generation capacity of heat storage tank of expansion power generation system

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113902040A (en) * 2021-11-15 2022-01-07 中国电力科学研究院有限公司 Method, system, equipment and storage medium for coordinating and optimizing electricity-heat comprehensive energy system
CN114370896A (en) * 2021-12-29 2022-04-19 贵州电网有限责任公司 Method for monitoring heating power generation capacity of heat storage tank of expansion power generation system
CN114370896B (en) * 2021-12-29 2024-03-19 贵州电网有限责任公司 Method for monitoring heating power generation capacity of heat storage tank of expansion power generation system

Similar Documents

Publication Publication Date Title
Woudstra et al. Thermodynamic evaluation of combined cycle plants
Zhang et al. Compressed air energy storage system with variable configuration for accommodating large-amplitude wind power fluctuation
Yang et al. Theoretical evaluation on the impact of heat exchanger in Advanced Adiabatic Compressed Air Energy Storage system
Wang et al. Combined heat and power control considering thermal inertia of district heating network for flexible electric power regulation
Spelling et al. Thermoeconomic optimization of a combined-cycle solar tower power plant
US9003763B2 (en) Compressed air energy system integrated with gas turbine
Kaviri et al. Modeling and multi-objective exergy based optimization of a combined cycle power plant using a genetic algorithm
Ghazi et al. Modeling and thermo-economic optimization of heat recovery heat exchangers using a multimodal genetic algorithm
CN110991857B (en) Method for evaluating wind power consumption capability of electric heating integrated energy system
CN112734591A (en) Electric heating comprehensive coordination scheduling method and device, equipment and medium
CN112347607B (en) Thermoelectric combined dispatching method based on convex relaxation
CN111191353B (en) Coordinated operation method of electric-thermal comprehensive energy system considering thermal inertia
Wang et al. Peak regulation performance study of the gas turbine combined cycle based combined heating and power system with gas turbine interstage extraction gas method
Attonaty et al. Thermodynamic analysis of a 200 MWh electricity storage system based on high temperature thermal energy storage
Gou et al. Dynamic simulation of a gas turbine for heat recovery at varying load and environment conditions
CN114962222B (en) Advanced adiabatic compressed air energy storage energy hub equipment and modeling method thereof
Li et al. Thermal-economy optimization for single/dual/triple-pressure HRSG of gas-steam combined cycle by multi-objective genetic algorithm
Wang et al. Integration of compressed air energy storage into combined heat and power plants: A solution to flexibility and economy
CN113886761A (en) Energy efficiency analysis and evaluation method for comprehensive energy system
Guo et al. Off-design characteristics of a novel integrated system of coal-fired cogeneration unit and compressed air energy storage
Guédez et al. Techno-economic comparative analysis of innovative combined cycle power plant layouts integrated with heat pumps and thermal energy storage
Srinivas et al. Parametric simulation of combined cycle power pant: A case study
CN113610305B (en) Optimized scheduling method of comprehensive energy system
Lee et al. A study on the variation of the performance and the cost of power generation in a combined heat and power plant with the change of the user facility’s return temperature
CN109815530B (en) Scheduling method of cogeneration system considering heat exchange saturation of waste heat boiler

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