CN111324850B - Considering heat supply network characteristics and heat load schedulable electric-heat combined scheduling method - Google Patents

Considering heat supply network characteristics and heat load schedulable electric-heat combined scheduling method Download PDF

Info

Publication number
CN111324850B
CN111324850B CN202010099345.2A CN202010099345A CN111324850B CN 111324850 B CN111324850 B CN 111324850B CN 202010099345 A CN202010099345 A CN 202010099345A CN 111324850 B CN111324850 B CN 111324850B
Authority
CN
China
Prior art keywords
heat
network
pipeline
time period
water
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.)
Active
Application number
CN202010099345.2A
Other languages
Chinese (zh)
Other versions
CN111324850A (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.)
Nantong University
Original Assignee
Nantong 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 Nantong University filed Critical Nantong University
Priority to CN202010099345.2A priority Critical patent/CN111324850B/en
Publication of CN111324850A publication Critical patent/CN111324850A/en
Application granted granted Critical
Publication of CN111324850B publication Critical patent/CN111324850B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • 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
    • 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 discloses an electric-heating combined dispatching method considering heat supply network characteristics and heat load schedulability, which is suitable for the field of optimized dispatching of electric-heating combined systems. On the basis of considering the thermal characteristics and the thermal load schedulability of the district heating network, the thermal characteristics of the district heating network are modeled, and the heat exchange station and the thermal users thereof are equivalent to the thermal load. An electric heating combined system optimization scheduling model considering the conventional constraint of a power system, the thermal characteristic constraint of a district heating network and the thermal load schedulability constraint is established by taking the minimum daily coal consumption of the system as an optimization target. The model is a mixed integer nonlinear programming model, and the output of the thermal power generating unit, the cogeneration unit and the wind generating set in each scheduling time period is obtained after the model is solved, so that a more practical and economic scheduling scheme is provided. The invention fully utilizes the heat storage capacity and the schedulability of the heat load of the regional heat supply network, enhances the peak regulation capacity of the thermoelectric generator set, effectively reduces the daily coal consumption of the system and promotes the wind power consumption.

Description

Electric-heat combined scheduling method considering heat supply network characteristics and heat load schedulability
Technical Field
The invention relates to the field of electric-heating combined system optimization scheduling, in particular to an electric-heating combined scheduling method considering heat supply network characteristics and heat load schedulability.
Background
In recent years, the development of wind power resources in China is continuously increased, and as the end of 2019, the installed capacity of wind power in China reaches 1.98 hundred million kilowatts, the total generated energy is 2914 hundred million kilowatt hours, and the installed capacity is increased by 8.9 percent on a same scale, so that the wind power is the 3 rd-rd main power source of a power system. Different from a conventional energy source unit, the wind power generation unit is restricted by natural characteristics, and the wind power has intermittent and random fluctuation characteristics, so that the operation difficulty of a power grid is increased and the operation efficiency is reduced when large-scale wind power grid connection is carried out, and even the wind is abandoned for limiting the power. In 2019, 1-9 months, the 'wind abandon' electricity quantity in China is up to 128 hundred million kilowatt hours, wherein the 'wind abandon' electricity quantity in three provinces of inner Mongolia, gansu and Xinjiang is more than 10 million kilowatt hours, so that the 'wind abandon' electricity limiting situation in China is severe.
At present, areas with severe electricity limiting situations of 'wind abandon' are mainly concentrated in northeast, northwest and north China areas in the winter heating period, and in the electric power system of the areas, the output of the cogeneration units is generally higher. By the heating period in winter, the cogeneration unit is limited by a mode of 'fixing power with heat', the peak regulation capability of the cogeneration unit is difficult to give full play, and the peak regulation method is an important reason for limiting the electricity by 'wind abandoning' on a large scale. Aiming at the problem, the invention provides a new technical scheme for promoting wind power consumption by fully utilizing the heat supply network characteristic and the heat load schedulability.
The physical characteristics of a thermodynamic system and an electric system are quite different, and the thermodynamic system and the electric system have thermal inertia and transmission delay characteristics. In a thermodynamic system, a heat source, a heat supply network, a heat exchange station and heat users thereof can store a large amount of heat energy. The heat supply network characteristics are mainly expressed in two aspects of hot water transmission delay and temperature loss, are closely related to the operation mode of the heat supply network, can promote thermoelectric decoupling, and further enhance the wind power receiving capacity of the power system. Meanwhile, the heat exchange station and the heat users thereof have thermal inertia, and the small-range fluctuation of the room temperature can not influence the thermal comfort of the users, so that the heat load is more schedulable than the electric load. In the optimal scheduling of the electric heating combined system, the wind power consumption can be promoted by fully utilizing the schedulability of the heat load. However, the existing electric-heat combined dispatching method is simple in modeling process of the thermodynamic system, only the steady-state operation process of the thermodynamic system is considered, the influence of the heat supply network characteristic and the heat load schedulability is not considered, the acceptance of the electric power system to wind power is restricted to a certain extent, and the engineering practical value and the economical efficiency of dispatching results are reduced.
Disclosure of Invention
The invention provides an electric-heat combined dispatching method considering heat supply network characteristics and heat load schedulability aiming at the defects in the prior art. On the basis of considering the thermal characteristics and the thermal load schedulability of the district heating network, the invention models the thermal characteristics of the district heating network and equates the heat exchange station and the heat user thereof as the thermal load. On the basis, an electric-heat combined system optimization scheduling model considering the conventional constraint of the power system, the thermal characteristic constraint of the district heating network and the thermal load schedulability constraint is established by taking the minimum daily coal consumption of the system as an optimization target, and a more economic and reasonable scheduling strategy is provided.
In order to achieve the purpose, the invention adopts the following technical scheme:
an electric-heat combined scheduling method considering heat supply network characteristics and heat load schedulability comprises the following steps:
step 1, considering the electric heating characteristics of a heat source, the hydraulic characteristics and the thermal characteristics of a regional heat supply network and the schedulability of a heat load, and establishing a regional heat supply system model;
step 2, establishing an optimal scheduling model of the electric heating combined system based on the regional heating system model, wherein the optimal scheduling model comprises the following steps: setting the minimum daily coal consumption of the electric heating combined system as an optimization target, and considering conventional constraints of a power system, thermal characteristic constraints of a regional heating network and schedulability constraints of heat load;
and 3, obtaining an electric-heat combined system optimized dispatching result considering the heat supply network characteristics and the heat load schedulability based on the regional heat supply system model and the electric-heat combined system optimized dispatching model.
As a preferred technical scheme of the invention: the electrothermal property of the heat source in step 1 is expressed by formula (1):
Figure BDA0002386412370000021
in the formula (1), k is an index of the cogeneration unit; t is a scheduling period index; c m,k 、C v,k The thermoelectric ratios of the cogeneration unit k under the working conditions of back pressure and air inlet are respectively; p CHP,k,t Electric power for the cogeneration unit k at time t; q CHP,k,t The thermal power of the cogeneration unit k in the time period t is obtained;
Figure BDA0002386412370000022
the upper limit of the electric power of the cogeneration unit k; omega CHP Is a combined heat and power generation unit set; Γ is a set of scheduling periods; e k Is a constant;
wherein the heat power Q of the cogeneration unit k in the time period t CHP,k,t Expressed by equation (2):
Figure BDA0002386412370000023
in the formula (2), C p Is the specific heat capacity of the hot water; m is f,t The hot water flow of the branch where the heat source f is located in the time period t; t is sg,f,t 、T sh,f,t Respectively providing hot water temperatures of the branch where the heat source f is located in the water supply network and the return water network at a time period t; f is a heat source index; omega HS Is a set of heat sources;
wherein, the hot water temperature T of the branch where the heat source f is located in the water supply and return network in the time period T sg,f,t 、T sh,f,t The interval values of (a) and (b) are expressed by the following formulas (3) and (4), respectively:
Figure BDA0002386412370000024
Figure BDA0002386412370000025
in the formulas (3) and (4),
Figure BDA0002386412370000026
the lower limit and the upper limit of the hot water temperature of a branch where a heat source is located in the water supply network are respectively set; />
Figure BDA0002386412370000031
Respectively is the lower limit and the upper limit of the hot water temperature of the branch where the heat source is located in the return water network.
As a preferred technical scheme of the invention: the hydraulic characteristics and the thermal characteristics of the regional heat supply network in the step 1 comprise pipeline water pressure loss constraint, node temperature mixing constraint, node flow continuity constraint, pipeline flow range constraint, pipeline transmission delay constraint and temperature loss constraint; the hydraulic loss constraint of the pipeline is expressed by formulas (5) and (6):
Figure BDA0002386412370000032
Figure BDA0002386412370000033
in the formulas (5) and (6),
Figure BDA0002386412370000034
respectively as node n in water supply and return network 1 Water pressure loss at time t;
Figure BDA0002386412370000035
Figure BDA0002386412370000036
respectively as node n in water supply and return network 2 Water pressure loss at time t; mu.s i The water pressure loss coefficient of the pipeline i; m is pg,i,t 、m ph,i,t Respectively the hot water flow of the pipeline i in the water supply network and the return water network in the time period t; i is a pipeline index; omega pipe Is a pipeline set;
wherein, the node n in the water supply and return network 1 、n 2 Expressed by equation (7):
Figure BDA0002386412370000037
in the formula (7), N start,i 、N end,i Respectively indexing a first node and a tail node of the pipeline i;
wherein, the node n in the water supply and return network 1 Water pressure loss at time t
Figure BDA0002386412370000038
The interval value of (2) is expressed by a formula (8), and a node n in a water supply and return network 2 Loss of water pressure in time period t>
Figure BDA0002386412370000039
The interval value of (a) is expressed by equation (9):
Figure BDA00023864123700000310
Figure BDA00023864123700000311
in the formulas (8) and (9),
Figure BDA00023864123700000312
respectively is the lower limit and the upper limit of the water pressure loss of the nodes in the water supply network;
Figure BDA00023864123700000313
respectively the lower limit and the upper limit of the water pressure loss of the nodes in the backwater network;
the node temperature hybrid constraint includes equations (10) - (13):
Figure BDA00023864123700000314
Figure BDA00023864123700000315
in the formulas (10) and (11), Ω start,j 、Ω end,j Respectively taking a heat supply network node j as a starting point and a terminal point;
Figure BDA00023864123700000316
Figure BDA00023864123700000317
respectively supplying water and returning water to the tail end temperature of the pipeline i in the network at the time t; t is ng,j,t 、T nh,j,t Respectively the hot water temperature of the node j in the water supply network and the water return network in the time period t; j is a heat supply network node index; omega node Is a heat supply network node set;
Figure BDA0002386412370000041
Figure BDA0002386412370000042
in the formulas (12) and (13),
Figure BDA0002386412370000043
respectively providing the head end temperature of the pipeline i in the water supply network and the head end temperature of the pipeline i in the water return network in a time period t;
the node traffic continuity constraints include equations (14), (15):
Figure BDA0002386412370000044
Figure BDA0002386412370000045
the pipeline flow range constraint includes equations (16), (17):
Figure BDA0002386412370000046
Figure BDA0002386412370000047
in the formulas (16), (17),
Figure BDA0002386412370000048
the upper limit of the hot water flow of the pipeline i in the water supply network and the water return network respectively;
the pipe transmission delay constraint comprises the following formulas (18) and (19):
Figure BDA0002386412370000049
Figure BDA00023864123700000410
in the formulae (18) and (19), τ pg,i,t 、τ ph,i,t Respectively the transmission delay time of a pipeline i in a water supply network and a water return network; v. of pg,i,t 、v ph,i,t Respectively the average flow velocity of hot water of a pipeline i in a water supply network and a water return network; l is i Is the length of conduit i; Δ t is the scheduling interval;
the pipeline temperature loss constraint includes equations (20) - (23):
Figure BDA00023864123700000411
Figure BDA00023864123700000412
in the formulas (20) and (21), Δ T pg,i,t 、ΔT ph,i,t Respectively the temperature drop of the pipeline i in the water supply network and the water return network in the time period t; λ is the heat transfer efficiency per unit length of the pipe; t is po,i,t The ambient temperature of the pipeline i at time t;
Figure BDA00023864123700000413
Figure BDA00023864123700000414
in the formulas (22), (23),
Figure BDA00023864123700000415
respectively considering the terminal temperature of a pipeline i in a water supply network and a return network in the time period t;
wherein, the head end temperature of the pipeline i in the water supply and return network in the time period t
Figure BDA0002386412370000051
Is expressed by formulas (24) and (25), and the temperature of the pipeline i in the water supply and return network at the tail end of the time period t is greater than or equal to the temperature of the water supply and return network at the tail end of the time period t>
Figure BDA0002386412370000052
The interval value of (2) is expressed by the following equations (26) and (27):
Figure BDA0002386412370000053
Figure BDA0002386412370000054
Figure BDA0002386412370000055
Figure BDA0002386412370000056
in the formulas (24), (25), (26) and (27),
Figure BDA0002386412370000057
the lower limit and the upper limit of the temperature of the pipeline in the water supply network are respectively set; />
Figure BDA0002386412370000058
The lower limit and the upper limit of the temperature of the pipeline in the water return network are respectively set.
As a preferred technical scheme of the invention: the schedulability of the thermal load in step 1 is represented by equation (28):
Figure BDA0002386412370000059
in the formula (28), Q HES,g,t The thermal power of the heat exchange station g in the time period t; u shape 1 、U 2 Respectively is a standard thermal load down-regulation coefficient and an upward floating coefficient; z is the time window coefficient of the equivalent thermal load; h j,t The standard heat load requirement of the heat supply network node j in the time period t; g is a heat exchange station index; omega HES Is a set of heat exchange stations;
wherein the heat exchange station g has a thermal power Q in a time period t HES,g,t Expressed by equation (29):
Figure BDA00023864123700000510
in the formula (29), m g,t The hot water flow of the branch where the heat exchange station g is located in the time period t; t is eg,g,t 、T eh,g,t Respectively providing hot water temperatures of the branch where the heat exchange station g is located in the water supply network and the return network at a time period t;
wherein, the hot water temperature T of the branch where the heat exchange station g is located in the water supply and return network in the time period T eg,g,t 、T eh,g,t The interval value of (c) is expressed by the following equations (30) and (31):
Figure BDA00023864123700000511
Figure BDA00023864123700000512
in the formulas (30) and (31),
Figure BDA00023864123700000513
respectively is the lower limit and the upper limit of the hot water temperature of a branch where a heat exchange station is located in a water supply network; />
Figure BDA0002386412370000061
The lower limit and the upper limit of the hot water temperature of the branch where the heat exchange station is located in the water return network are respectively set.
As a preferred technical scheme of the invention: the daily coal consumption of the electric heating combined system in the step 2 is expressed by a formula (32):
C＝min(C TU +C CHP ) (32)
in the formula (32), C is the daily coal consumption of the electric heating combined system; c TU The daily coal consumption of a thermoelectric generator set in the system; c CHP The daily coal consumption of the cogeneration unit in the system;
wherein, the daily coal consumption C of the medium-voltage generator set of the system TU Expressed by equation (33):
Figure BDA0002386412370000062
in the formula (33), b 0,l 、b 1,l 、b 2,l Respectively a quadratic term coefficient, a primary term coefficient and a constant term of coal consumption of the thermal power generating unit; p l,t Electric power of the thermal power generating unit l in a time period t; l is a thermal power generating unit indexing; omega TU The method comprises the steps of (1) collecting thermal power generating units;
wherein, the daily coal consumption of the cogeneration unit in the system is C CHP Expressed by equation (34):
Figure BDA0002386412370000063
in the formula (34), a 0,k 、a 1,k 、a 2,k 、a 3,k 、a 4,k 、a 5,k Is the k coal consumption coefficient of the cogeneration unit.
As a preferred technical scheme of the invention: the conventional constraints of the power system in the step 2 comprise power balance constraint, operation constraint, wind power constraint, climbing constraint, rotation standby constraint and power flow constraint; the power balance constraint is expressed using equation (35):
Figure BDA0002386412370000064
in the formula (35), P W,m,t The wind power of the wind turbine generator m in a time period t; d n,t The electrical load demand of the grid node n in the time period t; n is a power grid node index; omega WIND The method comprises the steps of (1) collecting a wind turbine generator set; omega BUS The method comprises the steps of (1) collecting power grid nodes;
the operating constraints include equations (36), (37), and (38):
Figure BDA0002386412370000065
Figure BDA0002386412370000066
Figure BDA0002386412370000067
in the formulae (36), (37) and (38), P l min 、P l max Respectively representing the lower limit and the upper limit of electric power of the thermal power generating unit l;
Figure BDA0002386412370000068
the lower limit of the electric power of the cogeneration unit k; />
Figure BDA0002386412370000069
Respectively is the lower limit and the upper limit of the thermal power of the cogeneration unit k;
the wind power constraint is expressed by a formula (39):
Figure BDA0002386412370000071
in the formula (39), the first and second groups,
Figure BDA0002386412370000072
the wind power upper limit of the wind turbine generator m in the scheduling time period t is set;
the hill climbing constraint includes equations (40), (41):
Figure BDA0002386412370000073
Figure BDA0002386412370000074
in the formulas (40) and (41), P l down 、P l up The lower climbing speed and the upper climbing speed of the electric power of the thermal power generating unit are respectively;
Figure BDA0002386412370000075
the lower climbing speed and the upper climbing speed of the electric power of the cogeneration unit are respectively;
the rotational standby constraint includes equations (42), (43), and (44):
Figure BDA0002386412370000076
Figure BDA0002386412370000077
Figure BDA0002386412370000078
in the equations (42), (43) and (44),
Figure BDA0002386412370000079
respectively providing upward and downward rotation reserve capacities for the thermal power generating unit l in the scheduling time t; RE up 、RE down Respectively providing standby requirements for upward and downward rotation of the system;
the power flow constraint is expressed by the formula (45):
Figure BDA00023864123700000710
in the formula (45), SF q,n Injecting a power offset coefficient for the transmission line q to the node n; f q The transmission capacity of the transmission line q; q is a transmission line index; omega LINE Is a power transmission line set.
Compared with the prior art, the electric-heat combined dispatching method considering the heat supply network characteristics and the heat load schedulability has the following technical effects by adopting the technical scheme:
(1) The mode of 'fixing the electricity with the heat' in the traditional electric heating combined system is broken through, and the coordinated operation of a thermodynamic system and an electric power system is promoted. (2) The heat storage capacity and the schedulability of heat load of a regional heat supply network are fully utilized, the peak regulation capacity of the thermoelectric unit is enhanced, the daily coal consumption of the system is effectively reduced, and wind power consumption is promoted.
Drawings
FIG. 1 is a schematic flow chart of a combined electric and heat dispatching method considering heat supply network characteristics and heat load schedulability;
FIG. 2 is a schematic diagram of an electric heating combined system;
FIG. 3 is a system daily load prediction curve;
FIG. 4 is a wind power forecast power curve;
FIG. 5 is a system daily net load curve;
FIG. 6 shows the thermal power optimization results of the system;
FIG. 7 is a wind power consumption comparison.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings. It is emphasized that the following description is merely exemplary in nature and is not intended to limit the scope of the invention or its application.
As shown in fig. 1, the present invention provides an electric-thermal combined scheduling method considering heat supply network characteristics and heat load schedulability, which comprises the following specific steps:
step 1, considering the electric heating characteristics of a heat source, the hydraulic characteristics and the thermal characteristics of a regional heat supply network and the schedulability of a heat load, and establishing a regional heat supply system model;
the modeling of the district heating system comprises the modeling of the electric heating characteristic of a heat source, the modeling of the hydraulic and thermal characteristics of a district heating network and the modeling of the schedulability of a heat load. Wherein the electrothermal property of the heat source is represented by the following formula:
Figure BDA0002386412370000081
in the formula, k is an index of a cogeneration unit; t is a scheduling period index; c m,k 、C v,k The thermoelectric ratios of the cogeneration unit k under the working conditions of back pressure and air inlet are respectively; p CHP,k,t Electric power of the cogeneration unit k in a time period t; q CHP,k,t The thermal power of the cogeneration unit k in the time period t is obtained;
Figure BDA0002386412370000082
the upper limit of the electric power of the cogeneration unit k; omega CHP Is a combined heat and power generation unit set; f is a scheduling time interval set; e k Is a constant;
wherein the heat power Q of the cogeneration unit k in the time period t CHP,k,t Represented by the formula:
Figure BDA0002386412370000083
in the formula, C p Is the specific heat capacity of the hot water; m is f,t The hot water flow of the branch where the heat source f is located in the time period t; t is sg,f,t 、T sh,f,t Respectively providing hot water temperatures of the branch where the heat source f is located in the water supply network and the return water network at a time period t; f is a heat source index; omega HS Is a set of heat sources;
wherein, the hot water temperature T of the branch where the heat source f is located in the water supply and return network in the time period T sg,f,t 、T sh,f,t The interval values of (a) and (b) are expressed by the following formulas (3) and (4), respectively:
Figure BDA0002386412370000084
Figure BDA0002386412370000085
in the formula (I), the compound is shown in the specification,
Figure BDA0002386412370000091
the lower limit and the upper limit of the hot water temperature of a branch where a heat source is located in the water supply network are respectively set;
Figure BDA0002386412370000092
the lower limit and the upper limit of the hot water temperature of a branch where a heat source is located in a return water network are respectively set;
the hydraulic characteristics and the thermal characteristics of the district heating network comprise pipeline water pressure loss constraint, node temperature mixing constraint, node flow continuity constraint, pipeline flow range constraint, pipeline transmission delay constraint and temperature loss constraint. The hydraulic loss constraint of the pipeline is expressed by formulas (5) and (6):
Figure BDA0002386412370000093
Figure BDA0002386412370000094
in the formula (I), the compound is shown in the specification,
Figure BDA0002386412370000095
respectively as node n in water supply and return network 1 Water pressure loss at time t;
Figure BDA0002386412370000096
respectively as node n in water supply and return network 2 Water pressure loss at time t; mu.s i The water pressure loss coefficient of the pipeline i; m is pg,i,t 、m ph,i,t Respectively the hot water flow of the pipeline i in the water supply network and the return water network in the time period t; i is a pipeline index; omega pipe Is a pipeline set;
wherein, the node n in the water supply and return network 1 、n 2 Represented by the formula:
Figure BDA0002386412370000097
in the formula, N start,i 、N end,i Respectively indexing a first node and a tail node of the pipeline i;
wherein, the node n in the water supply and return network 1 Loss of water pressure at time t
Figure BDA0002386412370000098
The interval value of (2) is expressed by a formula (8), and a node n in a water supply and return network 2 Loss of water pressure in time period t>
Figure BDA0002386412370000099
The interval value of (a) is expressed by equation (9): />
Figure BDA00023864123700000910
Figure BDA00023864123700000911
In the formula (I), the compound is shown in the specification,
Figure BDA00023864123700000912
respectively is the lower limit and the upper limit of the water pressure loss of the nodes in the water supply network; />
Figure BDA00023864123700000916
Respectively is the lower limit and the upper limit of the water pressure loss of the nodes in the backwater network;
the node temperature mixing constraint is expressed by the following formula (10) - (13):
Figure BDA00023864123700000913
Figure BDA00023864123700000914
in the formula, omega start,j 、Ω end,j Respectively taking a heat supply network node j as a starting point and a terminal point;
Figure BDA00023864123700000915
respectively supplying water and returning water to the tail end temperature of the pipeline i in the network at the time t; t is ng,j,t 、T nh,j,t Respectively the hot water temperature of the node j in the water supply network and the water return network in the time period t; j is a heat supply network node index; omega node Is a heat supply network node set;
Figure BDA0002386412370000101
Figure BDA0002386412370000102
in the formula (I), the compound is shown in the specification,
Figure BDA0002386412370000103
respectively providing the head end temperature of the pipeline i in the water supply network and the head end temperature of the pipeline i in the water return network in a time period t;
the node flow continuity constraint is expressed by the formulas (14) and (15):
Figure BDA0002386412370000104
Figure BDA0002386412370000105
the pipeline flow range constraint is expressed by the following formulas (16) and (17):
Figure BDA0002386412370000106
Figure BDA0002386412370000107
in the formula (I), the compound is shown in the specification,
Figure BDA0002386412370000108
the upper limit of the hot water flow of the pipeline i in the water supply network and the water return network respectively;
the pipeline transmission delay constraint is expressed by the formulas (18) and (19):
Figure BDA0002386412370000109
Figure BDA00023864123700001010
in the formula, τ pg,i,t 、τ ph,i,t Respectively the transmission delay time of a pipeline i in a water supply network and a water return network; v. of pg,i,t 、v ph,i,t Respectively the average flow velocity of hot water of a pipeline i in a water supply network and a water return network; l is i Is the length of conduit i; Δ t is the scheduling interval;
the pipeline temperature loss constraint includes equations (20) - (23):
Figure BDA00023864123700001011
/>
Figure BDA00023864123700001012
in the formula,. DELTA.T pg,i,t 、ΔT ph,i,t Respectively the temperature drop of the pipeline i in the water supply and return network in the time period t(ii) a λ is the heat transfer efficiency per unit length of the pipe; t is po,i,t The ambient temperature of the pipeline i at time t;
Figure BDA00023864123700001013
Figure BDA00023864123700001014
in the formula (I), the compound is shown in the specification,
Figure BDA00023864123700001015
the tail end temperatures of the pipeline i in the time period t in the water supply network and the water return network considering time delay are respectively considered;
wherein, the head end temperature of the pipeline i in the water supply and return network in the time period t
Figure BDA0002386412370000111
Is expressed by formulas (24) and (25), and the temperature of the pipeline i in the water supply and return network at the tail end of the time period t is greater than or equal to the temperature of the water supply and return network at the tail end of the time period t>
Figure BDA0002386412370000112
The interval value of (2) is expressed by the following equations (26) and (27):
Figure BDA0002386412370000113
Figure BDA0002386412370000114
Figure BDA0002386412370000115
Figure BDA0002386412370000116
in the formula (I), the compound is shown in the specification,
Figure BDA0002386412370000117
the lower limit and the upper limit of the temperature of the pipeline in the water supply network are respectively set; />
Figure BDA0002386412370000118
Respectively is the lower limit and the upper limit of the temperature of the pipeline in the water return network;
the schedulability of thermal load is represented by the following equation:
Figure BDA0002386412370000119
in the formula, Q HES,g,t The thermal power of the heat exchange station g in the time period t; u shape 1 、U 2 Respectively is a standard thermal load down-regulation coefficient and an upward floating coefficient; z is the time window coefficient of the equivalent thermal load; h j,t The standard heat load requirement of the heat supply network node j in the time period t; g is a heat exchange station index; omega HES Is a set of heat exchange stations;
wherein the heat exchange station g has a thermal power Q in a time period t HES,g,t Represented by the formula:
Figure BDA00023864123700001110
in the formula, m g,t The hot water flow of the branch where the heat exchange station g is located in the time period t; t is eg,g,t 、T eh,g,t Respectively providing hot water temperatures of the branch where the heat exchange station g is located in the water supply network and the return network at a time period t;
wherein, the hot water temperature T of the branch where the heat exchange station g is located in the water supply and return network in the time period T eg,g,t 、T eh,g,t The interval value of (c) is expressed by the following equations (30) and (31):
Figure BDA00023864123700001111
Figure BDA00023864123700001112
in the formula (I), the compound is shown in the specification,
Figure BDA00023864123700001113
respectively is the lower limit and the upper limit of the hot water temperature of a branch where a heat exchange station is located in a water supply network;
Figure BDA00023864123700001114
the lower limit and the upper limit of the hot water temperature of the branch where the heat exchange station is located in the water return network are respectively set. />
Step 2, establishing an optimal scheduling model of the electric heating combined system based on the regional heating system model, wherein the optimal scheduling model comprises the following steps: setting the minimum daily coal consumption of the electric heating combined system as an optimization target, and considering conventional constraints of the power system, thermal characteristic constraints of a regional heating network and schedulability constraints of heat load; wherein the daily coal consumption of the electric heating combined system is represented by the following formula:
C＝min(C TU +C CHP ) (32)
in the formula, C is the daily coal consumption of the electric heating combined system; c TU The daily coal consumption of a thermoelectric generator set in the system; c CHP The daily coal consumption of the cogeneration unit in the system;
wherein, the daily coal consumption C of the medium-voltage generator set of the system TU Represented by the formula:
Figure BDA0002386412370000121
in the formula, b 0,l 、b 1,l 、b 2,l Respectively a quadratic term coefficient, a primary term coefficient and a constant term of coal consumption of the thermal power generating unit; p l,t Electric power of the thermal power generating unit l in a time period t; l is a thermal power generating unit index; omega TU The method comprises the steps of (1) collecting thermal power generating units;
wherein, the daily coal consumption of the cogeneration unit in the system is C CHP Represented by the formula:
Figure BDA0002386412370000122
in the formula, a 0,k 、a 1,k 、a 2,k 、a 3,k 、a 4,k 、a 5,k The k coal consumption coefficient of the cogeneration unit;
the conventional constraints of the power system comprise power balance constraint, operation constraint, wind power constraint, climbing constraint, rotation standby constraint and power flow constraint. Wherein the power balance constraint is represented by:
Figure BDA0002386412370000123
in the formula, P W,m,t The wind power of the wind turbine generator m in a time period t; d n,t The electrical load demand of the grid node n in the time period t; n is a power grid node index; omega WIND The method comprises the steps of (1) collecting a wind turbine generator set; omega BUS The method comprises the steps of (1) collecting power grid nodes;
the operating constraints are expressed by equations (36) - (38):
Figure BDA0002386412370000124
Figure BDA0002386412370000125
Figure BDA0002386412370000126
in the formula, P l min 、P l max Respectively representing the lower limit and the upper limit of electric power of the thermal power generating unit l;
Figure BDA0002386412370000127
the lower limit of the electric power of the cogeneration unit k; />
Figure BDA0002386412370000128
Respectively is the lower limit and the upper limit of the thermal power of the cogeneration unit k;
the wind power constraint is represented by:
Figure BDA0002386412370000131
in the formula (I), the compound is shown in the specification,
Figure BDA0002386412370000132
the wind power upper limit of the wind turbine generator m in the scheduling time period t is set;
the climbing constraint is expressed by the formulas (40) and (41):
Figure BDA0002386412370000133
Figure BDA0002386412370000134
in the formula, P l down 、P l up The lower climbing speed and the upper climbing speed of the electric power of the thermal power generating unit are respectively;
Figure BDA0002386412370000135
the lower climbing speed and the upper climbing speed of the electric power of the cogeneration unit are respectively;
the rotational standby constraint is expressed by equations (42) - (44):
Figure BDA0002386412370000136
Figure BDA0002386412370000137
Figure BDA0002386412370000138
in the formula (I), the compound is shown in the specification,
Figure BDA0002386412370000139
respectively providing upward and downward rotation reserve capacities for the thermal power generating unit l in the scheduling time t; RE up 、RE down Respectively providing standby requirements for upward and downward rotation of the system;
the power flow constraint is represented by:
Figure BDA00023864123700001310
in the formula, SF q,n Injecting a power offset coefficient for the transmission line q to the node n; f q The transmission capacity of the transmission line q; q is a transmission line index; omega LINE Is a power transmission line set.
And 3, obtaining an electric-heat combined system optimized dispatching result considering the heat supply network characteristics and the heat load schedulability based on the regional heat supply system model and the electric-heat combined system optimized dispatching model.
The combined electric and heat dispatching method considering the heat supply network characteristics and the heat load schedulability proposed by the invention is described by a preferred example.
As shown in FIG. 2, the combined heat and power system consists of a 6-node power grid and a 6-node heat supply network. Comprises the following steps: the two conventional thermal power generating units are respectively a G1 positioned at a power grid node 1 and a G2 positioned at a power grid node 2; the two cogeneration units are respectively CHP 1 and CHP2 positioned at a power grid node 6; the wind turbine WF is also located at the grid node 6. The specific parameters of each unit of the system are shown in tables 1-3, and the parameters of the district heating system are shown in table 4.
TABLE 1 System output parameters of each unit
Figure BDA0002386412370000141
TABLE 2 Cogeneration Unit parameters
Figure BDA0002386412370000142
TABLE 3 coal consumption parameters of thermal power generating unit
Figure BDA0002386412370000143
TABLE 4 district heating System parameters
Figure BDA0002386412370000144
TABLE 5 pipeline parameters
Figure BDA0002386412370000145
As shown in fig. 2, the heat supply network is composed of 2 heat sources, 6 nodes and 10 heat supply pipelines, wherein the pipeline numbers 1, 2, 3, 4 and 5 are water supply pipelines, the pipeline numbers 6, 7, 8, 9 and 10 are water return pipelines, and the specific pipeline parameters are shown in table 5.
Note that the scheduling period interval in this embodiment is 15 minutes. Fig. 3 shows a system daily load prediction curve, and a wind power prediction power result is shown in fig. 4. In addition, a system daily net load curve (i.e., a value of the electrical load power minus the wind power predicted power) is shown in fig. 5.
To prove the rationality and effectiveness of the proposed model, the present embodiment sets two scenarios for comparison:
scene 1: and determining the thermal power of the cogeneration unit by system thermal load balance constraint without considering the heat supply network characteristic and the thermal load schedulability, wherein the specific expression is shown as the following formula:
Figure BDA0002386412370000151
scene 2: and (4) obtaining an optimal scheduling strategy of the electric-heating combined system by considering the heat supply network characteristics and the heat load schedulability.
The above data are adopted for simulation, and the daily scheduling results of the electric power and the thermal power in the system are respectively shown in a table 6 and a figure 6. As can be seen from table 6: the system can completely meet the requirements of electricity and heat load in operation, and the cogeneration unit needs to supply power and heat simultaneously, and the output electric power of the cogeneration unit is generally greater than that of the thermal power unit.
TABLE 6 System electric Power optimization results
Figure BDA0002386412370000152
Figure BDA0002386412370000161
As shown in fig. 6, the real-time thermal power of the cogeneration unit is not strictly equal to the standard thermal load. Because the heat supply network characteristic and the heat load schedulability are considered, the district heating system has certain heat storage capacity, which is specifically represented as: (1) When the thermal power of the cogeneration unit is greater than the standard thermal load, the district heating network stores a certain amount of thermal energy; (2) When the thermal power of the cogeneration unit is less than the standard thermal load, the thermal energy originally stored in the district heating network is released to satisfy the heating demand of the heat user and maintain the comfort of the user. As can be seen from fig. 7, from 8. During the two periods 0 to 8 and 17 to 0. In addition, due to the temperature loss characteristic of the district heating network, the total output thermal power of the cogeneration unit is slightly higher than the equivalent thermal load.
Table 7 shows the comparison of the operating economy of the combined heat and power system in different scenarios. As can be seen from the table, the daily coal consumption considering the heat supply network characteristics and the heat load schedulability is 7357.114 tons, the daily coal consumption not considering the heat supply network characteristics and the heat load schedulability is 7484.331 tons, the daily coal consumption of the system is reduced by 127.217 tons, and the economic benefit is remarkably improved.
TABLE 7 comparison of operating economics
Figure BDA0002386412370000171
As shown in fig. 7 and table 7, the "wind abandon" phenomenon is significantly improved in consideration of the heat supply network characteristics and the heat load schedulability, the total amount of the "wind abandon" in consideration of the heat supply network characteristics and the heat load schedulability is 259.581MWh, the total amount of the "wind abandon" in consideration of the heat supply network characteristics and the heat load schedulability is 387.803MWh, the total amount of the "wind abandon" is reduced by 128.222MWh, the wind abandon rate is reduced by 6.125%, and the wind power absorption is effectively promoted.
The invention provides an electric-heat combined dispatching method considering heat supply network characteristics and heat load schedulability, and compared with the prior art, the technical scheme has the following technical effects: (1) The mode of 'fixing the electricity with the heat' in the traditional electric heating combined system is broken through, and the coordinated operation of a thermodynamic system and an electric power system is promoted. (2) The heat storage capacity and the schedulability of heat load of a regional heat supply network are fully utilized, the peak regulation capacity of the thermoelectric unit is enhanced, the daily coal consumption of the system is effectively reduced, and wind power consumption is promoted.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only illustrative of the present invention, and are not intended to limit the scope of the present invention, and any person skilled in the art should understand that equivalent changes and modifications made without departing from the concept and principle of the present invention should fall within the protection scope of the present invention.

Claims (1)

1. An electric-heat combined scheduling method considering heat supply network characteristics and heat load schedulability is characterized by comprising the following steps of:
step 1, considering the electric heating characteristics of a heat source, the hydraulic characteristics and the thermal characteristics of a regional heat supply network and the schedulability of a heat load, and establishing a regional heat supply system model;
step 2, establishing an optimal scheduling model of the electric heating combined system based on the regional heating system model, wherein the optimal scheduling model comprises the following steps: setting the minimum daily coal consumption of the electric heating combined system as an optimization target, and considering conventional constraints of a power system, thermal characteristic constraints of a regional heating network and schedulability constraints of heat load;
step 3, obtaining an optimized dispatching result of the electric-heat combined system considering the heat supply network characteristics and the heat load schedulability based on the regional heat supply system model and the optimized dispatching model of the electric-heat combined system;
the electrothermal property of the heat source in step 1 is expressed by formula (1):
Figure FDA0004071403970000011
in the formula (1), k is an index of the cogeneration unit; t is a scheduling period index; c m,k 、C v,k The thermoelectric ratios of the cogeneration unit k under the working conditions of back pressure and air inlet are respectively set; p CHP,k,t Electric power of the cogeneration unit k in a time period t; q CHP,k,t The thermal power of the cogeneration unit k in the time period t is obtained;
Figure FDA0004071403970000012
the upper limit of the electric power of the cogeneration unit k; omega CHP Is a combined heat and power generation unit set; Γ is a set of scheduling periods; e k Is a constant;
wherein the heat power Q of the cogeneration unit k in the time period t CHP,k,t Expressed by equation (2):
Figure FDA0004071403970000013
in the formula (2), C p Is the specific heat capacity of the hot water; m is f,t Is a branch of the heat source fThe hot water flow in the time period t; t is sg,f,t 、T sh,f,t Respectively providing hot water temperatures of the branch where the heat source f is located in the water supply network and the return water network at a time period t; f is a heat source index; omega HS Is a set of heat sources;
wherein, the hot water temperature T of the branch where the heat source f is located in the water supply and return network in the time period T sg,f,t 、T sh,f,t The interval values of (a) and (b) are expressed by the following formulas (3) and (4), respectively:
Figure FDA0004071403970000014
Figure FDA0004071403970000015
in the formulas (3) and (4),
Figure FDA0004071403970000016
the lower limit and the upper limit of the hot water temperature of a branch where a heat source is located in the water supply network are respectively set;
Figure FDA0004071403970000017
the lower limit and the upper limit of the hot water temperature of a branch where a heat source is located in a return water network are respectively set;
the hydraulic characteristics and the thermal characteristics of the regional heat supply network in the step 1 comprise pipeline water pressure loss constraint, node temperature mixing constraint, node flow continuity constraint, pipeline flow range constraint, pipeline transmission delay constraint and temperature loss constraint; the hydraulic loss constraint of the pipeline is expressed by formulas (5) and (6):
Figure FDA0004071403970000021
Figure FDA0004071403970000022
in the formulas (5) and (6),
Figure FDA0004071403970000023
respectively as node n in water supply and return network 1 Water pressure loss at time t;
Figure FDA0004071403970000024
Figure FDA0004071403970000025
respectively as node n in water supply and return network 2 Water pressure loss at time t; mu.s i The water pressure loss coefficient of the pipeline i; m is pg,i,t 、m ph,i,t Respectively the hot water flow of the pipeline i in the water supply network and the return water network in the time period t; i is a pipeline index; omega pipe Is a pipeline set; />
Wherein, the node n in the water supply and return network 1 、n 2 Expressed by equation (7):
Figure FDA0004071403970000026
in the formula (7), N start,i 、N end,i Respectively indexing a first node and a tail node of the pipeline i;
wherein, the node n in the water supply and return network 1 Water pressure loss at time t
Figure FDA0004071403970000027
The interval value of (2) is expressed by a formula (8), and a node n in a water supply and return network 2 Loss of water pressure in time period t>
Figure FDA0004071403970000028
The interval value of (a) is expressed by equation (9):
Figure FDA0004071403970000029
Figure FDA00040714039700000210
in the formulae (8) and (9), Y g min 、Y g max Respectively is the lower limit and the upper limit of the water pressure loss of the nodes in the water supply network; y is h min 、Y h max Respectively is the lower limit and the upper limit of the water pressure loss of the nodes in the backwater network;
the nodal temperature mixing constraint includes equations (10) - (13):
Figure FDA00040714039700000211
Figure FDA00040714039700000212
in the formulas (10) and (11), Ω start,j 、Ω end,j Respectively taking a heat supply network node j as a starting point and a terminal point;
Figure FDA00040714039700000213
Figure FDA00040714039700000214
respectively supplying water and returning water to the tail end temperature of the pipeline i in the network at the time t; t is ng,j,t 、T nh,j,t Respectively the hot water temperature of the node j in the water supply network and the water return network in the time period t; j is a heat supply network node index; omega node Is a heat supply network node set;
Figure FDA00040714039700000215
Figure FDA00040714039700000216
in the formulas (12) and (13),
Figure FDA0004071403970000031
respectively providing the head end temperature of the pipeline i in the water supply network and the head end temperature of the pipeline i in the water return network in a time period t;
the node traffic continuity constraints include equations (14), (15):
Figure FDA0004071403970000032
Figure FDA0004071403970000033
the pipeline flow range constraint includes equations (16), (17):
Figure FDA0004071403970000034
Figure FDA0004071403970000035
in the formulas (16), (17),
Figure FDA0004071403970000036
the upper limit of the hot water flow of the pipeline i in the water supply network and the water return network respectively;
the pipe transmission delay constraint comprises the following formulas (18) and (19):
Figure FDA0004071403970000037
Figure FDA0004071403970000038
in the formulae (18) and (19), τ pg,i,t 、τ ph,i,t Respectively the transmission delay time of a pipeline i in a water supply network and a water return network; v. of pg,i,t 、v ph,i,t Respectively the average flow velocity of hot water of a pipeline i in a water supply network and a water return network; l is i Is the length of conduit i; Δ t is the scheduling interval;
the pipeline temperature loss constraint includes equations (20) - (23):
Figure FDA0004071403970000039
Figure FDA00040714039700000310
in the formulas (20) and (21), Δ T pg,i,t 、ΔT ph,i,t Respectively the temperature drop of the pipeline i in the water supply network and the water return network in the time period t; λ is the heat transfer efficiency per unit length of the pipe; t is po,i,t The ambient temperature of the pipeline i at time t;
Figure FDA00040714039700000311
Figure FDA00040714039700000312
in the formulas (22), (23),
Figure FDA00040714039700000313
respectively considering the terminal temperature of a pipeline i in a water supply network and a return network in the time period t;
wherein, the head end temperature of the pipeline i in the water supply and return network in the time period t
Figure FDA00040714039700000314
Is expressed by formulas (24) and (25), and the temperature of the pipeline i in the water supply and return network at the tail end of the time period t is greater than or equal to the temperature of the water supply and return network at the tail end of the time period t>
Figure FDA0004071403970000041
The interval value of (2) is expressed by the following equations (26) and (27):
Figure FDA0004071403970000042
Figure FDA0004071403970000043
Figure FDA0004071403970000044
Figure FDA0004071403970000045
in the formulas (24), (25), (26) and (27),
Figure FDA0004071403970000046
the lower limit and the upper limit of the temperature of the pipeline in the water supply network are respectively set;
Figure FDA0004071403970000047
respectively is the lower limit and the upper limit of the temperature of the pipeline in the water return network;
the schedulability of the thermal load in step 1 is represented by equation (28):
Figure FDA0004071403970000048
in the formula (28), Q HES,g,t The thermal power of the heat exchange station g in the time period t; u shape 1 、U 2 Respectively is a standard thermal load down-regulation coefficient and an upward floating coefficient; z is the time window coefficient of the equivalent thermal load; h j,t A standard thermal load demand for heat supply network node j at time period t; g is a heat exchange station index; omega HES Is a set of heat exchange stations;
wherein the heat exchange station g has a thermal power Q in a time period t HES,g,t Expressed by equation (29):
Figure FDA0004071403970000049
in the formula (29), m g,t The hot water flow of the branch where the heat exchange station g is located in the time period t; t is eg,g,t 、T eh,g,t Respectively providing hot water temperatures of the branch where the heat exchange station g is located in the water supply network and the return network at a time period t;
wherein, the hot water temperature T of the branch where the heat exchange station g is located in the water supply and return network in the time period T eg,g,t 、T eh,g,t The interval value of (c) is expressed by the following equations (30) and (31):
Figure FDA00040714039700000410
/>
Figure FDA00040714039700000411
in the formulas (30) and (31),
Figure FDA00040714039700000412
the lower limit and the upper limit of the hot water temperature of a branch where a heat exchange station is located in a water supply network are respectively set; />
Figure FDA00040714039700000413
The lower limit and the upper limit of the hot water temperature of a branch where a heat exchange station in a water return network is located are respectively set;
the daily coal consumption of the electric heating combined system in the step 2 is expressed by a formula (32):
C＝min(C TU +C CHP ) (32)
in the formula (32), C is the daily coal consumption of the electric heating combined system; c TU The daily coal consumption of a thermoelectric generator set in the system; c CHP The daily coal consumption of the cogeneration unit in the system;
wherein, the daily coal consumption C of the medium-voltage generator set of the system TU Expressed by equation (33):
Figure FDA0004071403970000051
in the formula (33), b 0,l 、b 1,l 、b 2,l Respectively a quadratic term coefficient, a primary term coefficient and a constant term of coal consumption of the thermal power generating unit; p l,t Electric power of the thermal power generating unit l in a time period t; l is a thermal power generating unit index; omega TU The method comprises the steps of (1) collecting thermal power generating units;
wherein, the daily coal consumption of the cogeneration unit in the system is C CHP Expressed by equation (34):
Figure FDA0004071403970000052
in the formula (34), a 0,k 、a 1,k 、a 2,k 、a 3,k 、a 4,k 、a 5,k The k coal consumption coefficient of the cogeneration unit;
the conventional constraints of the power system in the step 2 comprise power balance constraint, operation constraint, wind power constraint, climbing constraint, rotating standby constraint and power flow constraint; the power balance constraint is expressed using equation (35):
Figure FDA0004071403970000053
in the formula (35), P W,m,t For wind powerThe wind power of the unit m in the time period t; d n,t The electrical load demand of the grid node n in the time period t; n is a power grid node index; omega WIND The method comprises the steps of (1) collecting a wind turbine generator set; omega BUS The method comprises the steps of (1) collecting power grid nodes;
the operating constraints include equations (36), (37), and (38):
Figure FDA0004071403970000054
Figure FDA0004071403970000055
Figure FDA0004071403970000056
in the formulae (36), (37) and (38), P l min 、P l max Respectively representing the lower limit and the upper limit of electric power of the thermal power generating unit l;
Figure FDA0004071403970000057
the lower limit of the electric power of the cogeneration unit k; />
Figure FDA0004071403970000058
Respectively is the lower limit and the upper limit of the thermal power of the cogeneration unit k;
the wind power constraint is expressed by a formula (39):
Figure FDA0004071403970000059
in the formula (39), the first and second groups,
Figure FDA00040714039700000510
the wind power upper limit of the wind turbine generator m in the scheduling time period t is set;
the hill climbing constraint includes equations (40), (41):
Figure FDA0004071403970000061
Figure FDA0004071403970000062
in the formulas (40) and (41), P l down 、P l up The lower climbing speed and the upper climbing speed of the electric power of the thermal power generating unit are respectively;
Figure FDA0004071403970000063
the lower climbing speed and the upper climbing speed of the electric power of the cogeneration unit are respectively; />
The rotational standby constraint includes equations (42), (43), and (44):
Figure FDA0004071403970000064
Figure FDA0004071403970000065
Figure FDA0004071403970000066
in the equations (42), (43) and (44),
Figure FDA0004071403970000067
respectively providing upward and downward rotation reserve capacities for the thermal power generating unit l in the scheduling time t; RE up 、RE down Respectively providing standby requirements for upward and downward rotation of the system;
the power flow constraint is expressed by the formula (45):
Figure FDA0004071403970000068
in the formula (45), SF q,n Injecting a power offset coefficient for the transmission line q to the node n; f q The transmission capacity of the transmission line q; q is a transmission line index; omega LINE Is a power transmission line set.
CN202010099345.2A 2020-02-18 2020-02-18 Considering heat supply network characteristics and heat load schedulable electric-heat combined scheduling method Active CN111324850B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010099345.2A CN111324850B (en) 2020-02-18 2020-02-18 Considering heat supply network characteristics and heat load schedulable electric-heat combined scheduling method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010099345.2A CN111324850B (en) 2020-02-18 2020-02-18 Considering heat supply network characteristics and heat load schedulable electric-heat combined scheduling method

Publications (2)

Publication Number Publication Date
CN111324850A CN111324850A (en) 2020-06-23
CN111324850B true CN111324850B (en) 2023-04-07

Family

ID=71167087

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010099345.2A Active CN111324850B (en) 2020-02-18 2020-02-18 Considering heat supply network characteristics and heat load schedulable electric-heat combined scheduling method

Country Status (1)

Country Link
CN (1) CN111324850B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111682537B (en) * 2020-06-30 2022-06-10 广西大学 Electric heating combined system considering heat exchange delay characteristic and scheduling method thereof
CN111985165B (en) * 2020-07-22 2021-07-20 河海大学 Electric heating comprehensive energy system scheduling method considering heat storage characteristics of intelligent building
CN113190941A (en) * 2021-03-26 2021-07-30 国网辽宁省电力有限公司电力科学研究院 Centralized heating pipe network modeling method suitable for electric-heat joint scheduling
CN113761727B (en) * 2021-08-23 2024-02-23 中国人民解放军海军工程大学 Method for constructing optimal scheduling model of combined heat and power system containing distributed electric heat pump
CN113689043B (en) * 2021-08-25 2024-03-08 国网黑龙江省电力有限公司电力科学研究院 Electrothermal joint scheduling method considering start and stop of unit

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106849188B (en) * 2017-01-23 2020-03-06 中国电力科学研究院 Combined heat and power optimization method and system for promoting wind power consumption
CN110232640A (en) * 2019-01-15 2019-09-13 华北电力大学 It is a kind of towards wind electricity digestion the considerations of thermic load elasticity and heat supply network characteristic electric heating integrated distribution model

Also Published As

Publication number Publication date
CN111324850A (en) 2020-06-23

Similar Documents

Publication Publication Date Title
CN111324850B (en) Considering heat supply network characteristics and heat load schedulable electric-heat combined scheduling method
CN108258679B (en) Electric-thermal comprehensive energy system optimization scheduling method considering heat storage characteristics of heat supply network
CN106998079B (en) Modeling method of combined heat and power optimization scheduling model
CN108154309B (en) Energy internet economic dispatching method considering multi-load dynamic response of cold, heat and electricity
CN109190785B (en) Operation optimization method for electric-thermal coupling comprehensive energy system
CN109659927B (en) Energy storage capacity configuration method of comprehensive energy microgrid considering energy storage participation degree
CN111324849B (en) Electric heating combined system optimal scheduling method considering heat supply network characteristics
CN112736939B (en) Optimized capacity configuration method for hydrogen production and storage device of hydrogen-doped natural gas comprehensive energy system
CN110188492B (en) Combined cooling heating and power micro-grid optimized scheduling method considering heat supply network characteristics
CN110930073B (en) Day-ahead scheduling method for wind-light-photo-thermal combined power generation system considering price type demand response
CN110535128A (en) Based on the multizone integrated energy system coordinated dispatching method with energy comfort level
CN110797917A (en) Scheduling model of electric heating combined system
CN110620403A (en) Day-ahead scheduling method and system for collaborative operation of energy system considering renewable energy
CN113190941A (en) Centralized heating pipe network modeling method suitable for electric-heat joint scheduling
CN113725915A (en) Rural electric heating comprehensive energy system operation optimization method considering renewable energy uncertainty and thermal inertia
CN110991845B (en) Distributed cooperative scheduling method for electric-thermal coupling system
CN111476394B (en) Robust operation optimization method suitable for multi-energy systems such as electric heating gas system
CN114139379A (en) Simulation method of electric power-thermal steam coupling energy system
CN113313305A (en) Non-cooperative game-based comprehensive energy system optimization scheduling method
Zhao et al. Optimal scheduling method for electrical-thermal integrated energy system considering heat storage characteristics of heating network
CN113852139A (en) Wind-storage combined heating system optimal scheduling strategy considering electric heating demand response
CN112508730A (en) Wind power consumption strategy of comprehensive energy virtual power plant
Xi et al. Robust optimization of integrated energy system considering multiple thermal inertia and wind power uncertainty
Wang et al. Research on the electricity-gas coupling system with P2G to absorb surplus hydropower
CN112100828B (en) Electric heating system control method considering load quasi-dynamic characteristic of heating power network

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