CN113065729A - Comprehensive energy system optimization scheduling method and system considering exergy loss - Google Patents

Comprehensive energy system optimization scheduling method and system considering exergy loss Download PDF

Info

Publication number
CN113065729A
CN113065729A CN202010003464.3A CN202010003464A CN113065729A CN 113065729 A CN113065729 A CN 113065729A CN 202010003464 A CN202010003464 A CN 202010003464A CN 113065729 A CN113065729 A CN 113065729A
Authority
CN
China
Prior art keywords
energy system
heat supply
supply network
node
integrated energy
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
CN202010003464.3A
Other languages
Chinese (zh)
Inventor
张新鹤
李克成
刘铠诚
何桂雄
钟鸣
闫华光
黄伟
张志刚
王旭东
霍现旭
李树鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Tianjin Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Tianjin Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Tianjin Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202010003464.3A priority Critical patent/CN113065729A/en
Publication of CN113065729A publication Critical patent/CN113065729A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • 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

Landscapes

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

Abstract

The invention relates to a meter
Figure DDA0002354310160000011
The optimal scheduling method and system for the damaged comprehensive energy system comprise the following steps: according to the integrated energy system
Figure DDA0002354310160000012
Determining the optimal voltage of each energy supply node of a power supply network in the comprehensive energy system by loss; determining the optimal transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system according to the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system; and adjusting the transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system to be the optimal transformer transformation ratio. The invention providesThe technical scheme is as follows
Figure DDA0002354310160000013
As a standard for measuring the energy supply capacity and quality of the system, a power network and a heat supply network of a calculation comprehensive energy system are provided
Figure DDA0002354310160000014
Method of energy loss and integrated energy system
Figure DDA0002354310160000015
The minimum loss is used as a target to optimally schedule the comprehensive energy system, so that the comprehensive energy system can work outwards to the maximum extent, energy quality matching of the energy utilization of the comprehensive energy system is realized, and the energy utilization efficiency of the comprehensive energy system is improved.

Description

Comprehensive energy system optimization scheduling method and system considering exergy loss
Technical Field
The invention relates to energySource and thermodynamics field, in particular to a meter
Figure BDA0002354310140000012
A comprehensive energy system optimization scheduling method and system are disclosed.
Background
Energy with different forms and properties is transmitted in the comprehensive energy system, and a dissipation phenomenon is necessarily existed in the conversion process, so that the comprehensive energy conversion efficiency of the system is reduced, and losses generated in the energy transmission process of users with different energy levels, such as users facing cities, parks, district ballast levels and the like, are different, so that the energy loss in the process must be analyzed, and meanwhile, a related network optimization scheme is formulated to reduce the losses.
On the one hand, when a system is spontaneous, the direction in which the reaction occurs must be changed along the direction in which the quality of the total energy of the system is reduced, i.e. the principle of degradation of the energy. The traditional optimization scheduling method is provided on the basis that the energy utilization efficiency and the energy loss meet the first law of thermodynamics (no matter what conversion occurs in the energy form, the energy entering and exiting the multi-energy flow network is surely conserved in quantity), but the optimization scheduling method can only determine the quantity change trend of the energy in the transmission process, but cannot determine the energy quality transformation trend in the transmission process (the gradual reduction of the system energy quality is equal to the gradual reduction of the external acting capacity of the system), and the accuracy of the optimization scheduling method needs to be improved;
on the other hand, the comprehensive energy system covers energy flow sub-networks with different energy flow properties, such as a gas network, a power network, a heat supply network and the like, however, evaluation systems for energy efficiency are different among different systems, and the optimization result of a single subsystem in the traditional optimization scheduling method is mostly completed by sacrificing the benefits of other subsystems, and the applicability of the optimization scheduling method needs to be improved.
Disclosure of Invention
In view of the deficiencies of the prior art, the object of the present invention is a meter
Figure BDA0002354310140000013
Optimal scheduling method of damaged integrated energy system and method thereof
Figure BDA0002354310140000014
As a standard for measuring the energy supply capacity and quality of the system, a power network and a heat supply network of a calculation comprehensive energy system are provided
Figure BDA0002354310140000015
Method of energy loss and integrated energy system
Figure BDA0002354310140000016
The minimum loss is used as a target to optimally schedule the comprehensive energy system, so that the comprehensive energy system can work outwards to the maximum extent, energy quality matching of the energy utilization of the comprehensive energy system is realized, and the energy utilization efficiency of the comprehensive energy system is improved.
The purpose of the invention is realized by adopting the following technical scheme:
the invention provides a meter
Figure BDA0002354310140000017
In a method for optimal scheduling of a damaged integrated energy system, the improvement comprising:
according to the integrated energy system
Figure BDA0002354310140000018
Determining the optimal voltage of each energy supply node of a power supply network in the comprehensive energy system by loss;
determining the optimal transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system according to the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system;
and adjusting the transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system to be the optimal transformer transformation ratio.
Preferably, said system is based on an integrated energy system
Figure BDA0002354310140000022
The loss amount is determined to be the optimal voltage of each energy supply node of the power supply network in the integrated energy system, and the method comprises the following steps:
by integrated energy systems
Figure BDA0002354310140000023
Establishing an objective function of an optimized dispatching model of the comprehensive energy system for the purpose of minimum loss;
and solving the objective function of the comprehensive energy system optimization scheduling model based on the constraint condition corresponding to the objective function of the comprehensive energy system optimization scheduling model, and obtaining the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system.
Further, an objective function of the comprehensive energy system optimization scheduling model is determined according to the following formula:
Figure BDA0002354310140000021
in the formula, ClossFor integrated energy systems
Figure BDA0002354310140000024
Loss, c1Is electricity
Figure BDA0002354310140000025
Is/are as follows
Figure BDA0002354310140000026
Mass coefficient, c2Is heat
Figure BDA0002354310140000027
Is/are as follows
Figure BDA0002354310140000028
Mass coefficient, w is pressure
Figure BDA0002354310140000029
And electricity
Figure BDA00023543101400000210
The conversion coefficient of (a) is,
Figure BDA00023543101400000211
for electricity of i-th transmission line of power supply network in integrated energy system
Figure BDA00023543101400000212
The loss amount of the waste water is reduced,
Figure BDA00023543101400000213
for pressing of kth heat supply pipeline of heat supply network in comprehensive energy system
Figure BDA00023543101400000214
The loss amount of the waste water is reduced,
Figure BDA00023543101400000215
for heating of kth heat supply pipeline of heat supply network in integrated energy system
Figure BDA00023543101400000216
Loss, k is equal to [ 1-S ∈ ]hl]，ShlFor the number of heat supply pipelines of a heat supply network in an integrated energy system, i belongs to [ 1-S ]el]，SelThe number of transmission lines of the power supply network in the comprehensive energy system.
Further, the constraint conditions of the objective function of the integrated energy system optimization scheduling model include: equality constraint conditions, inequality constraint conditions of running output of a generator set, a gas turbine set and a circulating water pump in a power supply network, tolerance constraint conditions, transmission capacity constraint conditions, transformer tap gear constraint conditions and pipeline temperature constraint conditions.
The invention provides a meter
Figure BDA00023543101400000217
In an integrated energy system optimization scheduling system, the improvement comprising:
a first determination module for determining the system based on the integrated energy system
Figure BDA00023543101400000218
Determining the optimal voltage of each energy supply node of a power supply network in the comprehensive energy system by loss;
the second determination module is used for determining the optimal transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system according to the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system;
and the adjusting module is used for adjusting the transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system to the optimal transformer transformation ratio.
Preferably, the first determining module includes:
building units for integrated energy systems
Figure BDA00023543101400000219
Establishing an objective function of an optimized dispatching model of the comprehensive energy system for the purpose of minimum loss;
and the obtaining unit is used for solving the objective function of the comprehensive energy system optimization scheduling model based on the constraint condition corresponding to the objective function of the comprehensive energy system optimization scheduling model, and obtaining the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system.
Further, an objective function of the comprehensive energy system optimization scheduling model is determined according to the following formula:
Figure BDA0002354310140000031
in the formula, ClossFor integrated energy systems
Figure BDA0002354310140000032
Loss, c1Is electricity
Figure BDA0002354310140000033
Is/are as follows
Figure BDA0002354310140000034
The mass coefficient of the light beam is measured,c2is heat
Figure BDA0002354310140000035
Is/are as follows
Figure BDA0002354310140000036
Mass coefficient, w is pressure
Figure BDA0002354310140000039
And electricity
Figure BDA0002354310140000038
The conversion coefficient of (a) is,
Figure BDA00023543101400000316
for electricity of i-th transmission line of power supply network in integrated energy system
Figure BDA0002354310140000037
The loss amount of the waste water is reduced,
Figure BDA00023543101400000317
for pressing of kth heat supply pipeline of heat supply network in comprehensive energy system
Figure BDA00023543101400000310
The loss amount of the waste water is reduced,
Figure BDA00023543101400000318
for heating of kth heat supply pipeline of heat supply network in integrated energy system
Figure BDA00023543101400000311
Loss, k is equal to [ 1-S ∈ ]hl]，ShlFor the number of heat supply pipelines of a heat supply network in an integrated energy system, i belongs to [ 1-S ]el]，SelThe number of transmission lines of the power supply network in the comprehensive energy system.
Further, the constraint conditions of the objective function of the integrated energy system optimization scheduling model include: equality constraint conditions, inequality constraint conditions of running output of a generator set, a gas turbine set and a circulating water pump in a power supply network, tolerance constraint conditions, transmission capacity constraint conditions, transformer tap gear constraint conditions and pipeline temperature constraint conditions.
Compared with the closest prior art, the invention has the following beneficial effects:
the technical scheme provided by the invention is based on an integrated energy system
Figure BDA00023543101400000312
Determining the optimal voltage of each energy supply node of a power supply network in the comprehensive energy system by loss; determining the optimal transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system according to the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system; adjusting the transformer transformation ratio of each energy supply node of a power supply network in the comprehensive energy system to be the optimal transformer transformation ratio; the scheme is to
Figure BDA00023543101400000314
As a standard for measuring the energy supply capacity and quality of the system to integrate the energy system
Figure BDA00023543101400000313
The minimum loss is used as a target to optimally schedule the comprehensive energy system, so that the comprehensive energy system can work outwards to the maximum extent, energy quality matching of the energy utilization of the comprehensive energy system is realized, and the energy utilization efficiency of the comprehensive energy system is improved.
Drawings
FIG. 1 is a meter and
Figure BDA00023543101400000315
a flow chart of a comprehensive energy system optimization scheduling method of the loss;
FIG. 2 is a schematic diagram of an island integrated energy system according to an embodiment of the present invention;
FIG. 3 is a graph of algorithm iterations in an embodiment of the present invention;
FIG. 4 is a graph showing voltage amplitude variations in cooling and heating modes of an island integrated energy system according to an embodiment of the present invention;
FIG. 5 is a graph showing the change in mass flow of the pipeline in cooling and heating modes for an island energy system according to an embodiment of the present invention;
FIG. 6(a) is a diagram showing the variation of the heating temperature of each node of an island integrated energy system before and after optimization according to an embodiment of the present invention;
fig. 6(b) is a diagram illustrating changes of cooling temperatures at nodes of an island integrated energy system before and after optimization according to an embodiment of the present invention;
FIG. 7 is a graph comparing the cooling/heating capacity of an island integrated energy system according to an embodiment of the present invention;
FIG. 8 shows an embodiment of the invention before and after optimization of an island integrated energy system
Figure BDA0002354310140000042
Loss contrast graph;
FIG. 9 is a meter and
Figure BDA0002354310140000043
the structure diagram of the comprehensive energy system optimization scheduling system is lost.
Detailed Description
The following describes embodiments of the present invention in further detail with reference to the accompanying drawings.
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention provides a meter
Figure BDA0002354310140000044
As shown in fig. 1, the method for optimal scheduling of a damaged integrated energy system includes:
101, according to the comprehensive energy system
Figure BDA0002354310140000045
Determining the optimal voltage of each energy supply node of a power supply network in the comprehensive energy system by loss;
102, determining the optimal transformer transformation ratio of each energy supply node of the power supply network in the integrated energy system according to the optimal voltage of each energy supply node of the power supply network in the integrated energy system;
and 103, regulating the transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system to be the optimal transformer transformation ratio.
Specifically, the step 101 includes:
step a: by integrated energy systems
Figure BDA0002354310140000046
Establishing an objective function of an optimized dispatching model of the comprehensive energy system for the purpose of minimum loss;
step b: and solving the objective function of the comprehensive energy system optimization scheduling model based on the constraint condition corresponding to the objective function of the comprehensive energy system optimization scheduling model, and obtaining the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system.
Further, an objective function of the comprehensive energy system optimization scheduling model is determined according to the following formula:
Figure BDA0002354310140000041
in the formula, ClossFor integrated energy systems
Figure BDA00023543101400000411
Loss, c1Is electricity
Figure BDA00023543101400000410
Is/are as follows
Figure BDA0002354310140000049
Mass coefficient, c2Is heat
Figure BDA0002354310140000048
Is/are as follows
Figure BDA0002354310140000047
Mass coefficient, w is pressure
Figure BDA0002354310140000054
And electricity
Figure BDA0002354310140000055
The conversion coefficient of (a) is,
Figure BDA00023543101400000513
for electricity of i-th transmission line of power supply network in integrated energy system
Figure BDA0002354310140000058
The loss amount of the waste water is reduced,
Figure BDA00023543101400000514
for pressing of kth heat supply pipeline of heat supply network in comprehensive energy system
Figure BDA0002354310140000057
The loss amount of the waste water is reduced,
Figure BDA00023543101400000515
for heating of kth heat supply pipeline of heat supply network in integrated energy system
Figure BDA0002354310140000056
Loss, k is equal to [ 1-S ∈ ]hl]，ShlFor the number of heat supply pipelines of a heat supply network in an integrated energy system, i belongs to [ 1-S ]el]，SelThe number of transmission lines of a power supply network in the comprehensive energy system;
wherein the electricity of the ith transmission line of the power supply network in the integrated energy system is determined according to the following formula
Figure BDA0002354310140000059
Loss of volume
Figure BDA00023543101400000516
Figure BDA0002354310140000051
In the formula (I), the compound is shown in the specification,
Figure BDA00023543101400000517
for the electricity flowing from the x-th transmission line connected with the head end node of the i-th transmission line of the power supply network in the integrated energy system into the i-th transmission line of the power supply network in the integrated energy system
Figure BDA00023543101400000510
Figure BDA00023543101400000518
For electricity flowing from the i-th transmission line of the power supply network in the integrated energy system to the c-th transmission line connected with the tail end node thereof
Figure BDA00023543101400000511
RiResistance, P, of the i-th transmission line of a power supply network in an integrated energy systemi,sFor the electric power flowing into the head-end node of the i-th transmission line of the supply network in the integrated energy system, Ui,sThe voltage of the head end node of the ith transmission line of a power supply network in the comprehensive energy system is x ∈ [ 1-S ]x]，SxC belongs to [ 1-S ] is the sum of the transmission lines connected with the head end node of the ith transmission line of the power supply network in the integrated energy systemc]，ScThe total sum of the transmission lines connected with the terminal node of the ith transmission line of the power supply network in the comprehensive energy system;
in the preferred embodiment of the present invention, the transfer pipe is made of a uniform material, and the first end and the last end of the transfer pipe are connected only
When the pipeline is transmitted, the voltage of the head end node and the tail end node of the transmission line of the power supply network in the comprehensive energy system meets the following formula:
Figure BDA0002354310140000052
in the formula of UBFor the voltage of the head-end node, U, of the transmission line of a power supply network in an integrated energy systemCIs the voltage of the end node of the transmission line of the supply network in the integrated energy system, R is the resistance of the transmission line of the supply network in the integrated energy system, PBInputting power for a head-end node of a transmission line of a power supply network in the integrated energy system;
therefore, the electric power of the head and tail end nodes of the transmission line of the power supply network in the integrated energy system satisfies the following formula:
Figure BDA0002354310140000053
in the formula, PCOutputting power for a terminal node of a transmission line of a power supply network in the integrated energy system;
therefore, the first and the last end nodes of the transmission line of the power supply network in the integrated energy system are electrified
Figure BDA00023543101400000512
Satisfies the following formula:
Figure BDA0002354310140000061
in the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000067
inputting electricity for a head-end node of a transmission line of a power supply network in an integrated energy system
Figure BDA0002354310140000066
Figure BDA0002354310140000068
For power supply networks in integrated energy systemsEnd node output of transmission line
Figure BDA0002354310140000069
U0The quiescent voltage of a node of a power supply network in the comprehensive energy system is taken as 0;
when the head end and the tail end of the transmission pipeline are popularized to be connected with a plurality of transmission pipelines, the head end and the tail end of the transmission pipeline of the power supply network in the comprehensive energy system are powered by the nodes
Figure BDA00023543101400000610
Satisfies the following formula:
Figure BDA0002354310140000062
in the formula, PBzInputting electricity for a head-end node of a transmission line of a power supply network in an integrated energy system
Figure BDA00023543101400000611
To sum up, i.e.
Figure BDA0002354310140000063
Figure BDA00023543101400000612
Electricity for transmission lines of power supply networks in integrated energy systems
Figure BDA00023543101400000613
The amount of change.
Determining the pressure of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formula
Figure BDA00023543101400000614
Loss of volume
Figure BDA00023543101400000615
Figure BDA0002354310140000064
In the formula, RkIs the fluid flow resistance, x, of the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe volume flow of the kth heat supply pipeline of a heat supply network in the comprehensive energy system;
in the preferred embodiment of the present invention, assuming that the transmission pipelines are uniform in texture, when the head end and the tail end of the transmission pipeline are both connected to only one transmission pipeline, the voltage of the head end node and the tail end node of the heat supply pipeline of the heat supply network in the integrated energy system satisfies the following formula:
pC＝pB-RyxQ
in the formula, pBPressure, p, for the head-end node of the heat supply pipe of the heat supply network in the integrated energy systemCFor the pressure intensity, R, of the tail end node of the heat supply pipeline in the comprehensive energy systemyIs the fluid flow resistance, x, of the heat supply pipeline in the integrated energy systemQThe volume flow of the heat supply pipeline in the comprehensive energy system;
therefore, the pressure power of the head end node and the tail end node of the heat supply pipeline of the heat supply network in the comprehensive energy system satisfies the following formula:
PCy＝PBy-RyxQ 2
in the formula, PByThe node voltage power of the head end of a heat supply pipeline of a heat supply network in the comprehensive energy system is obtained; pCyThe pressure power of the tail end node of the heat supply pipeline of the heat supply network in the comprehensive energy system is obtained;
therefore, the first end node and the tail end node of the heat supply pipeline of the heat supply network in the integrated energy system are pressed
Figure BDA00023543101400000616
Satisfies the following formula:
Figure BDA0002354310140000065
in the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000078
head end node pressure of heat supply pipeline for heat supply network in comprehensive energy system
Figure BDA0002354310140000075
PCyTerminal node pressure of heat supply pipeline for heat supply network in comprehensive energy system
Figure BDA0002354310140000076
When the head end and the tail end of the transmission pipeline are popularized to be connected with a plurality of transmission pipelines, the head end node pressure of the heat supply pipeline of the heat supply network in the comprehensive energy system is realized
Figure BDA0002354310140000077
Satisfies the following formula:
Figure BDA0002354310140000071
in the formula, xQzFor the total volume flow input of the heating pipes of the heating network in the integrated energy system, i.e.
Figure BDA0002354310140000079
Figure BDA00023543101400000710
For pressing of heat supply pipelines of heat supply network in integrated energy system
Figure BDA00023543101400000711
The amount of change.
Determining heat of kth heat supply pipeline of heat supply network in comprehensive energy system according to following formula
Figure BDA00023543101400000712
Loss of volume
Figure BDA00023543101400000713
Figure BDA0002354310140000072
In the formula (I), the compound is shown in the specification,
Figure BDA00023543101400000714
the heat of the kth heat supply pipeline flowing into the kth heat supply pipeline of the heat supply network in the integrated energy system from the kth heat supply pipeline connected with the head end node of the kth heat supply pipeline of the heat supply network in the integrated energy system
Figure BDA00023543101400000715
Figure BDA00023543101400000716
The heat flowing from the kth heat supply pipeline of the heat supply network in the integrated energy system to the kth heat supply pipeline connected with the tail end node of the heat supply network
Figure BDA00023543101400000717
Rho is the density of the heat supply pipeline transmission working medium, caThe specific heat capacity of the working medium is transmitted for the heat supply pipeline,
Figure BDA00023543101400000718
for the temperature, X, of the head end node of the kth heat supply pipeline of a heat supply network in an integrated energy systemT0For the value of the temperature silence, epsilon, of the heating network in the integrated energy systemkThe composite heat transfer coefficient of the kth heat supply pipeline of a heat supply network in an integrated energy system is shown, wherein z belongs to [ 1-S ]z]，SzThe number of heat supply pipelines connected with the head end node of the kth heat supply pipeline of the heat supply network in the integrated energy system is b belongs to [ 1-S ]b]，SbThe number of the heat supply pipelines connected with the terminal nodes of the kth heat supply pipeline of the heat supply network in the comprehensive energy system.
In the preferred embodiment of the present invention, assuming that the transmission pipelines are uniform in texture, when the head end and the tail end of the transmission pipeline are both connected to only one transmission pipeline, the temperature of the node at the head end and the tail end of the heat supply pipeline of the heat supply network in the integrated energy system satisfies the following formula: :
Figure BDA0002354310140000073
Figure BDA0002354310140000074
heat transfer coefficient, tau, for insulated pipes with n layers of insulationiI is the heat transfer coefficient of each layer of heat insulating material from inside to outside, diI belongs to n and is the radius of each layer of material; x is the number ofQIs the volume flow of the heat supply pipeline of the heat supply network in the comprehensive energy system; rho, c are the density and specific heat capacity of the hot/cold working medium in turn, XT0Is a measure of the intensity of the silence of the heat energy.
Therefore, the heat of the head and tail end nodes of the heat supply pipeline of the heat supply network in the integrated energy system
Figure BDA0002354310140000087
Satisfies the following formula:
Figure BDA0002354310140000081
in the formula (I), the compound is shown in the specification,
Figure BDA00023543101400000810
for head end node heat of heat supply pipeline of heat supply network in comprehensive energy system
Figure BDA0002354310140000088
Figure BDA00023543101400000811
For terminal node heat of heat supply pipeline of heat supply network in integrated energy system
Figure BDA0002354310140000089
When the head end and the tail end of the transmission pipeline are both only connected with a plurality of transmission pipelines, the heat supply network in the comprehensive energy system
Head end node of heat supply pipelinePress and press
Figure BDA00023543101400000812
Satisfies the following formula:
Figure BDA0002354310140000082
in the formula, xQzFor the total volume flow input of the heating pipes of the heating network in the integrated energy system, i.e.
Figure BDA0002354310140000083
Figure BDA00023543101400000813
For pressing of heat supply pipelines of heat supply network in integrated energy system
Figure BDA00023543101400000814
The amount of change.
Further, electricity is determined as follows
Figure BDA00023543101400000815
Is/are as follows
Figure BDA00023543101400000816
Mass coefficient c1：
Figure BDA0002354310140000084
In the formula, phi1As a factor for evaluating the energy level of electric energy, μ1Is an energy level factor of the electrical energy;
heat was determined as follows
Figure BDA00023543101400000817
Is/are as follows
Figure BDA00023543101400000818
Mass coefficient c2：
Figure BDA0002354310140000085
In the formula, phi2Is an energy level evaluation factor of thermal energy, mu2An energy level factor that is thermal energy;
determining the resistance R of the ith transmission line of a power supply network in an integrated energy system according to the following formulai：
Figure BDA0002354310140000086
In the formula, LiLength of i-th transmission line of power supply network in integrated energy system, AiThe transmission channel area K of the ith transmission line of a power supply network in an integrated energy systemiThe resistivity of a conductor material of the ith transmission line of a power supply network in the comprehensive energy system;
determining the fluid flow resistance R of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formulak：
Figure BDA0002354310140000091
In the formula, mukFluid viscosity, L, for the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe length r of the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe inner radius of the kth heat supply pipeline of a heat supply network in the comprehensive energy system;
determining the volume flow x of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formulak：
Figure BDA0002354310140000092
In the formula, pk,sPressure intensity p of head end node of kth heat supply pipeline of heat supply network in integrated energy systemk,mFor supplying heat in an integrated energy systemThe pressure of the tail end node of the kth heat supply pipeline of the network;
determining the composite heat transfer coefficient epsilon of the kth heat supply pipeline of the heat supply network in the comprehensive energy system according to the following formulak：
Figure BDA0002354310140000093
In the formula (d)k(σ+1)Radius of the material of the sigma +1 layer of the kth heat supply pipeline of the heat supply network in the integrated energy system, dkσRadius, gamma, of the material of the sigma-th layer of the kth heat supply pipeline of a heat supply network in an integrated energy systemkσThe heat transfer coefficient of the material of the sigma layer of the kth heat supply pipeline of the heat supply network in the comprehensive energy system is sigma ∈ (1-psi-1), and psi is the layer number of the heat transfer material of the kth heat supply pipeline of the heat supply network in the comprehensive energy system;
further, the constraint conditions of the objective function of the integrated energy system optimization scheduling model include: equality constraint conditions, inequality constraint conditions of running output of a generator set, a gas turbine set and a circulating water pump in a power supply network, tolerance constraint conditions, transmission capacity constraint conditions, transformer tap gear constraint conditions and pipeline temperature constraint conditions.
Further, the equality constraints are determined as follows:
Figure BDA0002354310140000101
in the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000103
for the input electrical power of the q-th supply node of the supply network in the integrated energy system,
Figure BDA0002354310140000104
for the electrical demand of the f-th load node of the power supply network in the integrated energy system,
Figure BDA0002354310140000106
for the current of the ith transmission line of the power supply network in the integrated energy system,
Figure BDA0002354310140000105
for the voltage variation of the ith transmission line of the power supply network in the integrated energy system,
Figure BDA0002354310140000107
for the input thermal power of the d-th heating node of the heating network in the integrated energy system,
Figure BDA0002354310140000108
for the thermal demand of the r-th load node of the heating network in the integrated energy system,
Figure BDA0002354310140000109
is the entropy of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure BDA00023543101400001010
the temperature variation of the kth heat supply pipeline of a heat supply network in an integrated energy system is shown, and q belongs to [ 1-S ]q]，SqIs the total number of power supply nodes of a power supply network in an integrated energy system, f belongs to [ 1-S ]f]，SfThe total number of load nodes of a power supply network in an integrated energy system is d ∈ 1-Sd]，SdThe total number of heat supply nodes of a heat supply network in the comprehensive energy system is r belongs to [ 1-S ]r]The total number of load nodes of a heat supply network in the comprehensive energy system;
determining the inequality constraint conditions of the running output of the generator set, the gas generator set and the circulating water pump in the power supply network according to the following formula:
Figure BDA0002354310140000102
in the formula, PFDFor the operating output, P, of a generator set in a power supply network in an integrated energy systemFD,maxFor comprehensive energy systemUpper limit value of operating output, P, of generator set in medium-voltage networkFD,minIs the lower limit value of the running output, P, of the generator set in the power supply network in the comprehensive energy systemRQFor the running output, P, of gas turbine units in the heat supply network in the integrated energy systemRQ,maxThe upper limit value P of the operating output of the gas turbine unit in the heat supply network in the comprehensive energy systemRQ,minThe lower limit value P of the running output of the gas unit in the heat supply network in the comprehensive energy systemXHThe running output, P, of a circulating water pump in a heat supply network in the comprehensive energy systemXH,maxThe upper limit value P of the running output of a circulating water pump in a heat supply network in the comprehensive energy systemXH,minThe lower limit value of the running output of a circulating water pump in a heat supply network in the comprehensive energy system;
the tolerance constraint is determined as follows:
Figure BDA0002354310140000111
in the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000113
for the voltage at the h-th node in the supply network in the integrated energy system,
Figure BDA0002354310140000114
the maximum tolerance value of the h node in the power supply network in the integrated energy system,
Figure BDA0002354310140000115
for the minimum tolerance value of the h node in the power supply network in the integrated energy system,
Figure BDA0002354310140000116
for the temperature of the jth node in the heating network in the integrated energy system,
Figure BDA0002354310140000119
the maximum tolerance value of the temperature of the J-th node in the heat supply network in the integrated energy system,
Figure BDA0002354310140000118
for the minimum tolerance value of the temperature of the J-th node in the heat supply network in the integrated energy system,
Figure BDA0002354310140000117
is the working medium pressure of the kth heat supply pipeline in a heat supply network in the comprehensive energy system,
Figure BDA00023543101400001110
the maximum tolerance value of the working medium pressure of the kth heat supply pipeline in the heat supply network in the comprehensive energy system,
Figure BDA00023543101400001111
the minimum tolerance value of the working medium pressure of the kth heat supply pipeline in the heat supply network in the comprehensive energy system,
Figure BDA00023543101400001112
is the air pressure of the tth gas turbine set in the heat supply network in the comprehensive energy system,
Figure BDA00023543101400001113
is the maximum tolerance value of the air pressure of the tth gas turbine set in the heat supply network in the comprehensive energy system,
Figure BDA00023543101400001114
the minimum tolerance value of the air pressure of the tth gas turbine set in a heat supply network in the comprehensive energy system is J epsilon (1-S)J)，SJH is the total number of nodes in the heat supply network in the integrated energy system and belongs to (1-S)h)，ShIs the total number of nodes in a power supply network in the integrated energy system, and belongs to (1-S)t)，StThe total number of the gas units in the heat supply network in the comprehensive energy system;
determining the transmission capacity constraint as follows:
Figure BDA0002354310140000112
in the formula (I), the compound is shown in the specification,
Figure BDA00023543101400001115
for the current transfer capacity of the i-th transmission line of the power supply network in the integrated energy system,
Figure BDA00023543101400001116
for the minimum current transfer capacity of the i-th transmission line of the power supply network in the integrated energy system,
Figure BDA00023543101400001117
for the maximum current transfer capacity of the i-th transmission line of the power supply network in the integrated energy system,
Figure BDA00023543101400001118
the volume flow transmission capacity of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure BDA00023543101400001119
the maximum volume flow transmission capacity of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure BDA00023543101400001120
the minimum volume flow transmission capacity of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure BDA0002354310140000124
the heat energy transmission capacity of the kth heat supply pipeline of the heat supply network in the comprehensive energy system,
Figure BDA0002354310140000125
the maximum heat energy transmission capacity of the kth heat supply pipeline of the heat supply network in the comprehensive energy system,
Figure BDA0002354310140000126
the minimum heat energy transmission capacity of the kth heat supply pipeline of a heat supply network in the comprehensive energy system is obtained;
determining the transformer tap gear constraints as follows:
Figure BDA0002354310140000121
in the formula, TapδFor the delta-th transformer tap position of the supply network in an integrated energy system,
Figure BDA0002354310140000127
the lower limit value can be adjusted for the delta-th transformer tap position of the power supply network in the integrated energy system,
Figure BDA0002354310140000128
the upper limit value of the delta-th transformer tap gear of a power supply network in an integrated energy system can be adjusted, and delta belongs to [ 1-S ]Tap]，STapThe number of transformer taps of a power supply network in the integrated energy system;
determining the pipeline temperature constraint as follows:
Figure BDA0002354310140000122
in the formula, CsA first coefficient matrix, C, for the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy systemrIs a first coefficient matrix, X, of the node temperature between the water return pipelines of the heat supply network in the integrated energy systemTsIs a temperature matrix of the load heat flow inlet, XTrIs a temperature matrix of the load heat flux outlet, XT0A matrix of silence values for temperature, bsA second coefficient matrix of the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy system, brAnd the second coefficient matrix is the temperature of the nodes between the water return pipelines of the heat supply network in the comprehensive energy system.
Further, a first coefficient matrix C of the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy system is determined according to the following formulas：
Figure BDA0002354310140000123
In the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000129
for a node upsilon in a heat supply network in an integrated energy systemaTemperature and node gammaaA first coefficient between temperatures, va，γa∈(γN)，γNThe method comprises the steps of collecting nodes among heat supply pipelines in a heat supply network in an integrated energy system;
second coefficient matrix b for determining temperature of nodes between heat supply pipelines of heat supply network in comprehensive energy system according to following formulas：
Figure BDA0002354310140000131
In the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000139
for a node upsilon in a heat supply network in an integrated energy systemaA second coefficient of temperature, T being the transposed sign;
wherein the hot working medium is supplied from the node gammaaThrough the pipe MaFlow direction node vaThe method comprises the following steps:
if node vaIs not only processed by node gammaaSupplied thermal medium and node gammaaIs a load node, then:
Figure BDA0002354310140000132
Figure BDA0002354310140000133
if node vaIs not only processed by node gammaaSupplied thermal medium and node gammaaIs notA load node, then:
Figure BDA0002354310140000134
Figure BDA0002354310140000135
Figure BDA0002354310140000136
if node vaIs only composed of node gammaaA hot working medium is supplied to
Figure BDA0002354310140000137
Figure BDA0002354310140000138
In the above formula, the first and second carbon atoms are,
Figure BDA00023543101400001310
heat supply pipeline M for heat supply network in comprehensive energy systemaThe volume flow rate of (a) is,
Figure BDA00023543101400001311
heat supply pipeline M for heat supply network in comprehensive energy systemaThe composite heat transfer coefficient of (a) is,
Figure BDA00023543101400001312
heat supply pipeline M for heat supply network in comprehensive energy systemaThe length of (a) of (b),
Figure BDA00023543101400001313
Figure BDA00023543101400001314
to synthesizeTo node upsilon in heat supply network in energy systemaTotal number of heat supply pipes for supplying hot working medium, XT0Is a value of the silence state for the temperature,
Figure BDA00023543101400001315
for node gamma in heat supply network in comprehensive energy systemaThe temperature of the mixture is controlled by the temperature,
Figure BDA00023543101400001316
for a node upsilon in a heat supply network in an integrated energy systemaA first coefficient between the temperature and the temperature of the node of the temperature;
determining a first coefficient matrix C of the node temperature between water return pipelines of a heat supply network in an integrated energy system according to the following formular：
Figure BDA0002354310140000141
In the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000148
for node eta in heat supply network in comprehensive energy systemaTemperature and node phiaFirst coefficient between temperatures, phia,ηa∈(1～χN)，χNThe total number of nodes between water return pipelines in a heat supply network in the comprehensive energy system;
determining a second coefficient matrix b of the node temperature between the water return pipelines of the heat supply network in the comprehensive energy system according to the following formular
Figure BDA0002354310140000142
In the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000149
for node eta in heat supply network in comprehensive energy systemaA second coefficient of temperature;
wherein the hot working medium is supplied from a node phiaThroughPipeline HaFlow direction node ηaThe method comprises the following steps:
if node etaaNot only by the node phiaThe supplied hot working medium is as follows:
Figure BDA0002354310140000143
Figure BDA0002354310140000144
Figure BDA0002354310140000145
if node etaaOnly by node phiaWhen supplying hot working medium, then:
Figure BDA0002354310140000146
Figure BDA0002354310140000147
in the above formula, the first and second carbon atoms are,
Figure BDA00023543101400001410
for node phi in heat supply network in comprehensive energy systemaThe outlet temperature of (a) is set,
Figure BDA00023543101400001411
for node phi in heat supply network in comprehensive energy systemaThe volume flow rate of (a) is,
Figure BDA00023543101400001412
heat supply pipeline H for heat supply network in comprehensive energy systemaThe composite heat transfer coefficient of (a) is,
Figure BDA00023543101400001413
heat supply pipeline H for heat supply network in comprehensive energy systemaThe length of (a) of (b),
Figure BDA00023543101400001414
for node eta in heat supply network in comprehensive energy systemaA first coefficient between the temperature and its own node temperature,
Figure BDA0002354310140000154
Figure BDA0002354310140000155
for a heat supply network in an integrated energy system to a node etaaTotal number of return pipes supplying hot working medium.
Specifically, the step 102 includes:
determining an optimal transformer transformation ratio for the h-th energy supply node of a power supply network in an integrated energy system according to the following formula
Figure BDA0002354310140000156
Figure BDA0002354310140000151
In the formula of Ueh,bThe voltage which is accessed by the bus for the h-th energy supply node of the power supply network in the integrated energy system,
Figure BDA0002354310140000157
the optimal voltage of the h energy supply node of the power supply network in the integrated energy system is h epsilon (1-S)h)，ShIs the total number of nodes in the power supply network in the integrated energy system.
In an embodiment of the present invention, an analysis is performed on an integrated energy system providing power and heat at an island, and fig. 2 is a topological structure diagram of the integrated energy system at the island, where the loads of power network nodes and the loads of heat supply network nodes are shown in table 1 and table 2, respectively. The power network and the heat supply network are interconnected through three CHP units, a circulating water pump and the like to realize energy coupling.
TABLE 1
Power network node numbering 1 2 3 4 5 6
load/(MW) 0.2 0 0.5 0.5 0.2 0.2
TABLE 2
Figure BDA0002354310140000152
In fig. 2, the units G1, G2, G3 are a gas turbine unit, an extraction steam turbine unit and a reciprocating internal combustion engine unit, respectively, the power network comprises 9 buses and 5 loads, and the heating network comprises 32 nodes and 32 pipelines. The state variable and the controlled variable of the examples are shown in tables 3 and 4, respectively.
TABLE 3
Figure BDA0002354310140000153
Figure BDA0002354310140000161
TABLE 4
Figure BDA0002354310140000162
In this example, the dimension of the particle swarm algorithm is 7, the learning factors c1 and c2 of the particle swarm algorithm are set to be 2, the inertia factor w is set to be 1, the iteration number of the algorithm is 60, and the particle number is 40.
Iteration of the algorithm: in the calculation example, the selection of fitness function is caused by system line transmission
Figure BDA0002354310140000163
The reciprocal of the loss amount can be used to obtain the convergence characteristic by describing the optimal fitness value of each generation in the population, and fig. 3 is an iteration curve of a particle swarm algorithm in two modes of cooling and heating for 60 times.
As can be seen from fig. 3, the particle swarm algorithm is used for calculation, the fitness value of each generation of particles is shown as small real points, and the small real points are connected by lines to form an algorithm convergence characteristic curve, which reflects the convergence trend of the algorithm; the fitness curve of the example basically reaches the optimum after 30 iterations, so that the method has strong convergence capability and high convergence speed.
And (4) comparing the results: the optimal value of the objective function and the optimal value of the control variable can be obtained through a particle swarm algorithm, under the optimal particle result, the voltage amplitude of the comprehensive energy system in the cooling and heating modes is shown in fig. 4, the pipeline mass flow of the comprehensive energy system in the cooling and heating modes is shown in fig. 5, the changes of the heating temperature of each node of the comprehensive energy system before and after optimization are respectively shown in fig. 6(a), and the changes of the cooling temperature of each node of the comprehensive energy system before and after optimization are respectively shown in fig. 6 (b);
as can be seen in FIG. 4, the voltage amplitude of each node in the optimized power system calculated by the particle swarm optimization is improved, and the voltage amplitude of each power node after optimization is obviously reduced compared with that before optimization, which is beneficial to reducing system power
Figure BDA0002354310140000164
Loss of (2). From data analysis, the optimization results of the regional comprehensive energy system on the power network in the heating mode and the cooling mode are very close.
As can be seen from fig. 5, in the heat supply mode, the mass flow in the heat supply system pipeline after being optimized by the particle swarm algorithm is reduced to a certain extent as compared with that before being optimized, but the reduction range is small, and the mass flow in the cooling system pipeline after being optimized by the particle swarm algorithm is obviously reduced as compared with that before being optimized. Thereby, the mass flow rate of the regional integrated energy system in the cooling/heating mode is reduced, thereby reducing the system pressure
Figure BDA0002354310140000165
Loss of loss, while also contributing to improved system safety.
As can be seen in fig. 6(a) and 6(b), the node temperature curves of the two modes before and after optimization are approximately consistent. However, since the mass flow rate of each node is reduced after optimization, the quality of energy supply of the system must be improved in order to ensure that the cold/heat load requirements of the system are met, so that the temperature of a plurality of thermodynamic nodes of the system in the heating mode after optimization is slightly increased compared with the temperature of the thermodynamic nodes before optimization, and the temperature of each thermodynamic node in the cooling mode is averagely reduced by 2 ℃ compared with the temperature of each thermodynamic node before optimization. Under different modes, the system is ensured by adjusting the temperature of the nodes of the thermodynamic network
Figure BDA0002354310140000166
The total amount of (A) is stable.
From the optimization results of the cooling mode and the heating mode of the regional integrated energy system, the cooling mode transmits the cold energy station under the same power conditionThe mass flow rate is much greater than the mass flow rate required for thermal energy transfer in the heating mode. This is because under certain ambient temperature conditions, the cooling temperature of the cold source node is typically a few degrees celsius to a dozen degrees celsius, while the heating temperature of the heat source system node is typically sixty to ninety degrees celsius, which results in cooling per mass flow rate
Figure BDA0002354310140000172
Generally less than heat
Figure BDA0002354310140000173
That is, the working capacity of cold water is generally less than that of hot water. As shown in fig. 7, when the temperature of the source node is constant, the mass flow rate increase caused by the increase of the cold demand in the cooling mode is much larger than the mass flow rate increase caused by the increase of the heat demand in the heating mode when the cold/heat demand is gradually increased. In this case, when the capacity of the pipe is constant, the cooling capacity of the cooling system is very limited, which limits the radiation range of the cooling system.
It can be seen that the adoption meter
Figure BDA0002354310140000174
After the comprehensive energy system operation optimization method is adopted, the mass flow rate of the system is reduced in different modes, so that the safety margin of the system is effectively increased, and the safety and reliability of the system are enhanced; under the condition that the system parameters are not changed much, as shown in fig. 8, under the heat supply mode of the original island comprehensive energy system,
Figure BDA0002354310140000177
the loss is 0.0306MW, after optimization
Figure BDA0002354310140000175
The loss is 0.0157MW, total system
Figure BDA0002354310140000176
The loss rate of the product is 48.91%; in the cold supply mode, the air conditioner is arranged,
Figure BDA0002354310140000178
the loss is 0.0534MW, after being optimized
Figure BDA00023543101400001710
The loss is 0.0349MW, total system
Figure BDA0002354310140000179
The loss reduction ratio reaches 34.59 percent, the energy utilization efficiency of the comprehensive energy system is effectively improved, the energy quality matching principle of the system is realized, and the effectiveness of the method is proved.
The invention provides a meter
Figure BDA00023543101400001711
As shown in fig. 9, the system for optimizing and scheduling integrated energy system includes:
a first determination module for determining the system based on the integrated energy system
Figure BDA00023543101400001712
Determining the optimal voltage of each energy supply node of a power supply network in the comprehensive energy system by loss;
the second determination module is used for determining the optimal transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system according to the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system;
and the adjusting module is used for adjusting the transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system to the optimal transformer transformation ratio.
Specifically, the first determining module includes:
building units for integrated energy systems
Figure BDA00023543101400001713
Establishing an objective function of an optimized dispatching model of the comprehensive energy system for the purpose of minimum loss;
and the solving unit is used for solving the objective function of the comprehensive energy system optimization scheduling model based on the constraint condition corresponding to the objective function of the comprehensive energy system optimization scheduling model, and acquiring the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system.
Further, an objective function of the comprehensive energy system optimization scheduling model is determined according to the following formula:
Figure BDA0002354310140000171
in the formula, ClossFor integrated energy systems
Figure BDA00023543101400001714
Loss, c1Is electricity
Figure BDA00023543101400001715
Is/are as follows
Figure BDA00023543101400001716
Mass coefficient, c2Is heat
Figure BDA00023543101400001717
Is/are as follows
Figure BDA00023543101400001718
Mass coefficient, w is pressure
Figure BDA00023543101400001720
And electricity
Figure BDA00023543101400001721
The conversion coefficient of (a) is,
Figure BDA00023543101400001722
for electricity of i-th transmission line of power supply network in integrated energy system
Figure BDA00023543101400001719
The loss amount of the waste water is reduced,
Figure BDA0002354310140000184
for pressing of kth heat supply pipeline of heat supply network in comprehensive energy system
Figure BDA0002354310140000185
The loss amount of the waste water is reduced,
Figure BDA0002354310140000186
for heating of kth heat supply pipeline of heat supply network in integrated energy system
Figure BDA0002354310140000187
Loss, k is equal to [ 1-S ∈ ]hl]，ShlFor the number of heat supply pipelines of a heat supply network in an integrated energy system, i belongs to [ 1-S ]el]，SelThe number of transmission lines of a power supply network in the comprehensive energy system;
wherein the electricity of the ith transmission line of the power supply network in the integrated energy system is determined according to the following formula
Figure BDA0002354310140000188
Loss of volume
Figure BDA0002354310140000189
Figure BDA0002354310140000181
In the formula (I), the compound is shown in the specification,
Figure BDA00023543101400001810
for the electricity flowing from the x-th transmission line connected with the head end node of the i-th transmission line of the power supply network in the integrated energy system into the i-th transmission line of the power supply network in the integrated energy system
Figure BDA00023543101400001811
Figure BDA00023543101400001812
For the current flowing from the ith transmission line of the power supply network in the integrated energy system to be connected with the tail end node of the ith transmission lineElectricity of the c-th transmission line
Figure BDA00023543101400001813
RiResistance, P, of the i-th transmission line of a power supply network in an integrated energy systemi,sFor the electric power flowing into the head-end node of the i-th transmission line of the supply network in the integrated energy system, Ui,sThe voltage of the head end node of the ith transmission line of a power supply network in the comprehensive energy system is x ∈ [ 1-S ]x]，SxC belongs to [ 1-S ] is the sum of the transmission lines connected with the head end node of the ith transmission line of the power supply network in the integrated energy systemc]，ScThe total sum of the transmission lines connected with the terminal node of the ith transmission line of the power supply network in the comprehensive energy system;
determining the pressure of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formula
Figure BDA00023543101400001815
Loss of volume
Figure BDA00023543101400001814
Figure BDA0002354310140000182
In the formula, RkIs the fluid flow resistance, x, of the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe volume flow of the kth heat supply pipeline of a heat supply network in the comprehensive energy system;
determining heat of kth heat supply pipeline of heat supply network in comprehensive energy system according to following formula
Figure BDA00023543101400001816
Loss of volume
Figure BDA00023543101400001817
Figure BDA0002354310140000183
In the formula (I), the compound is shown in the specification,
Figure BDA00023543101400001818
the heat of the kth heat supply pipeline flowing into the kth heat supply pipeline of the heat supply network in the integrated energy system from the kth heat supply pipeline connected with the head end node of the kth heat supply pipeline of the heat supply network in the integrated energy system
Figure BDA00023543101400001819
Figure BDA00023543101400001820
The heat flowing from the kth heat supply pipeline of the heat supply network in the integrated energy system to the kth heat supply pipeline connected with the tail end node of the heat supply network
Figure BDA00023543101400001821
Rho is the density of the heat supply pipeline transmission working medium, caThe specific heat capacity of the working medium is transmitted for the heat supply pipeline,
Figure BDA0002354310140000195
for the temperature, X, of the head end node of the kth heat supply pipeline of a heat supply network in an integrated energy systemT0For the value of the temperature silence, epsilon, of the heating network in the integrated energy systemkThe composite heat transfer coefficient of the kth heat supply pipeline of a heat supply network in an integrated energy system is shown, wherein z belongs to [ 1-S ]z]，SzThe number of heat supply pipelines connected with the head end node of the kth heat supply pipeline of the heat supply network in the integrated energy system is b belongs to [ 1-S ]b]，SbThe number of the heat supply pipelines connected with the terminal nodes of the kth heat supply pipeline of the heat supply network in the comprehensive energy system.
Further, electricity is determined as follows
Figure BDA0002354310140000196
Is/are as follows
Figure BDA0002354310140000197
Mass coefficient c1：
Figure BDA0002354310140000191
In the formula, phi1As a factor for evaluating the energy level of electric energy, μ1Is an energy level factor of the electrical energy;
heat was determined as follows
Figure BDA0002354310140000198
Is/are as follows
Figure BDA0002354310140000199
Mass coefficient c2：
Figure BDA0002354310140000192
In the formula, phi2Is an energy level evaluation factor of thermal energy, mu2An energy level factor that is thermal energy;
determining the resistance R of the ith transmission line of a power supply network in an integrated energy system according to the following formulai：
Figure BDA0002354310140000193
In the formula, LiLength of i-th transmission line of power supply network in integrated energy system, AiThe transmission channel area K of the ith transmission line of a power supply network in an integrated energy systemiThe resistivity of a conductor material of the ith transmission line of a power supply network in the comprehensive energy system;
determining the fluid flow resistance R of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formulak：
Figure BDA0002354310140000194
In the formula, mukIn an integrated energy systemFluid viscosity, L, of the kth heating pipeline of a heating networkkThe length r of the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe inner radius of the kth heat supply pipeline of a heat supply network in the comprehensive energy system;
determining the volume flow x of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formulak：
Figure BDA0002354310140000201
In the formula, pk,sPressure intensity p of head end node of kth heat supply pipeline of heat supply network in integrated energy systemk,mThe pressure of the tail end node of the kth heat supply pipeline of the heat supply network in the comprehensive energy system is obtained;
determining the composite heat transfer coefficient epsilon of the kth heat supply pipeline of the heat supply network in the comprehensive energy system according to the following formulak：
Figure BDA0002354310140000202
In the formula (d)k(σ+1)Radius of the material of the sigma +1 layer of the kth heat supply pipeline of the heat supply network in the integrated energy system, dkσRadius, gamma, of the material of the sigma-th layer of the kth heat supply pipeline of a heat supply network in an integrated energy systemkσThe heat transfer coefficient of the material of the sigma layer of the kth heat supply pipeline of the heat supply network in the comprehensive energy system is sigma e (1-psi-1), and psi is the layer number of the heat transfer material of the kth heat supply pipeline of the heat supply network in the comprehensive energy system.
Further, the constraint conditions of the objective function of the integrated energy system optimization scheduling model include: equality constraint conditions, inequality constraint conditions of running output of a generator set, a gas turbine set and a circulating water pump in a power supply network, tolerance constraint conditions, transmission capacity constraint conditions, transformer tap gear constraint conditions and pipeline temperature constraint conditions.
Further, the equality constraints are determined as follows:
Figure BDA0002354310140000203
in the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000204
for the input electrical power of the q-th supply node of the supply network in the integrated energy system,
Figure BDA0002354310140000205
for the electrical demand of the f-th load node of the power supply network in the integrated energy system,
Figure BDA0002354310140000206
for the current of the ith transmission line of the power supply network in the integrated energy system,
Figure BDA0002354310140000207
for the voltage variation of the ith transmission line of the power supply network in the integrated energy system,
Figure BDA0002354310140000208
for the input thermal power of the d-th heating node of the heating network in the integrated energy system,
Figure BDA0002354310140000209
for the thermal demand of the r-th load node of the heating network in the integrated energy system,
Figure BDA00023543101400002010
is the entropy of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure BDA00023543101400002011
the temperature variation of the kth heat supply pipeline of a heat supply network in an integrated energy system is shown, and q belongs to [ 1-S ]q]，SqIs the total number of power supply nodes of a power supply network in an integrated energy system, f belongs to [ 1-S ]f]，SfFor comprehensive energy systemThe total number of load nodes of the medium power supply network is d from [1 to S ∈d]，SdThe total number of heat supply nodes of a heat supply network in the comprehensive energy system is r belongs to [ 1-S ]r]The total number of load nodes of a heat supply network in the comprehensive energy system;
determining the inequality constraint conditions of the running output of the generator set, the gas generator set and the circulating water pump in the power supply network according to the following formula:
Figure BDA0002354310140000211
in the formula, PFDFor the operating output, P, of a generator set in a power supply network in an integrated energy systemFD,maxIs the upper limit value of the operating output, P, of the generator set in the power supply network in the comprehensive energy systemFD,minIs the lower limit value of the running output, P, of the generator set in the power supply network in the comprehensive energy systemRQFor the running output, P, of gas turbine units in the heat supply network in the integrated energy systemRQ,maxThe upper limit value P of the operating output of the gas turbine unit in the heat supply network in the comprehensive energy systemRQ,minThe lower limit value P of the running output of the gas unit in the heat supply network in the comprehensive energy systemXHThe running output, P, of a circulating water pump in a heat supply network in the comprehensive energy systemXH,maxThe upper limit value P of the running output of a circulating water pump in a heat supply network in the comprehensive energy systemXH,minThe lower limit value of the running output of a circulating water pump in a heat supply network in the comprehensive energy system;
the tolerance constraint is determined as follows:
Figure BDA0002354310140000212
in the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000213
for the voltage at the h-th node in the supply network in the integrated energy system,
Figure BDA0002354310140000214
the maximum tolerance value of the h node in the power supply network in the integrated energy system,
Figure BDA0002354310140000215
for the minimum tolerance value of the h node in the power supply network in the integrated energy system,
Figure BDA0002354310140000216
for the temperature of the jth node in the heating network in the integrated energy system,
Figure BDA0002354310140000217
the maximum tolerance value of the temperature of the J-th node in the heat supply network in the integrated energy system,
Figure BDA0002354310140000218
for the minimum tolerance value of the temperature of the J-th node in the heat supply network in the integrated energy system,
Figure BDA0002354310140000219
is the working medium pressure of the kth heat supply pipeline in a heat supply network in the comprehensive energy system,
Figure BDA00023543101400002110
the maximum tolerance value of the working medium pressure of the kth heat supply pipeline in the heat supply network in the comprehensive energy system,
Figure BDA0002354310140000223
the minimum tolerance value of the working medium pressure of the kth heat supply pipeline in the heat supply network in the comprehensive energy system,
Figure BDA0002354310140000224
is the air pressure of the tth gas turbine set in the heat supply network in the comprehensive energy system,
Figure BDA0002354310140000225
is the maximum tolerance value of the air pressure of the tth gas turbine set in the heat supply network in the comprehensive energy system,
Figure BDA0002354310140000226
the minimum tolerance value of the air pressure of the tth gas turbine set in a heat supply network in the comprehensive energy system is J epsilon (1-S)J)，SJH is the total number of nodes in the heat supply network in the integrated energy system and belongs to (1-S)h)，ShIs the total number of nodes in a power supply network in the integrated energy system, and belongs to (1-S)t)，StThe total number of the gas units in the heat supply network in the comprehensive energy system;
determining the transmission capacity constraint as follows:
Figure BDA0002354310140000221
in the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000227
for the current transfer capacity of the i-th transmission line of the power supply network in the integrated energy system,
Figure BDA0002354310140000228
for the minimum current transfer capacity of the i-th transmission line of the power supply network in the integrated energy system,
Figure BDA0002354310140000229
for the maximum current transfer capacity of the i-th transmission line of the power supply network in the integrated energy system,
Figure BDA00023543101400002210
the volume flow transmission capacity of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure BDA00023543101400002211
the maximum volume flow transmission capacity of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure BDA00023543101400002212
for supplying heat in an integrated energy systemThe minimum volume flow transmission capacity of the kth heat supply pipeline of the network,
Figure BDA00023543101400002213
the heat energy transmission capacity of the kth heat supply pipeline of the heat supply network in the comprehensive energy system,
Figure BDA00023543101400002214
the maximum heat energy transmission capacity of the kth heat supply pipeline of the heat supply network in the comprehensive energy system,
Figure BDA00023543101400002215
the minimum heat energy transmission capacity of the kth heat supply pipeline of a heat supply network in the comprehensive energy system is obtained;
determining the transformer tap gear constraints as follows:
Figure BDA0002354310140000222
in the formula, TapδFor the delta-th transformer tap position of the supply network in an integrated energy system,
Figure BDA00023543101400002216
the lower limit value can be adjusted for the delta-th transformer tap position of the power supply network in the integrated energy system,
Figure BDA00023543101400002217
the upper limit value of the delta-th transformer tap gear of a power supply network in an integrated energy system can be adjusted, and delta belongs to [ 1-S ]Tap]，STapThe number of transformer taps of a power supply network in the integrated energy system;
determining the pipeline temperature constraint as follows:
Figure BDA0002354310140000231
in the formula, CsAs a comprehensive energy systemFirst coefficient matrix, C, of the junction temperatures between the heat supply lines of the in-system heat supply networkrIs a first coefficient matrix, X, of the node temperature between the water return pipelines of the heat supply network in the integrated energy systemTsIs a temperature matrix of the load heat flow inlet, XTrIs a temperature matrix of the load heat flux outlet, XT0A matrix of silence values for temperature, bsA second coefficient matrix of the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy system, brAnd the second coefficient matrix is the temperature of the nodes between the water return pipelines of the heat supply network in the comprehensive energy system.
Further, a first coefficient matrix C of the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy system is determined according to the following formulas：
Figure BDA0002354310140000232
In the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000235
for a node upsilon in a heat supply network in an integrated energy systemaTemperature and node gammaaA first coefficient between temperatures, va，γa∈(γN)，γNThe method comprises the steps of collecting nodes among heat supply pipelines in a heat supply network in an integrated energy system;
second coefficient matrix b for determining temperature of nodes between heat supply pipelines of heat supply network in comprehensive energy system according to following formulas：
Figure BDA0002354310140000233
In the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000236
for a node upsilon in a heat supply network in an integrated energy systemaA second coefficient of temperature, T being the transposed sign;
wherein the hot working medium is supplied from the node gammaaThrough the tubeWay MaFlow direction node vaThe method comprises the following steps:
if node vaIs not only processed by node gammaaSupplied thermal medium and node gammaaIs a load node, then:
Figure BDA0002354310140000234
Figure BDA0002354310140000241
if node vaIs not only processed by node gammaaSupplied thermal medium and node gammaaIf the node is not a load node, then:
Figure BDA0002354310140000242
Figure BDA0002354310140000243
Figure BDA0002354310140000244
if node vaIs only composed of node gammaaA hot working medium is supplied to
Figure BDA0002354310140000245
Figure BDA0002354310140000246
In the above formula, the first and second carbon atoms are,
Figure BDA0002354310140000248
heat supply pipeline M for heat supply network in comprehensive energy systemaThe volume flow rate of (a) is,
Figure BDA0002354310140000249
heat supply pipeline M for heat supply network in comprehensive energy systemaThe composite heat transfer coefficient of (a) is,
Figure BDA00023543101400002410
heat supply pipeline M for heat supply network in comprehensive energy systemaThe length of (a) of (b),
Figure BDA00023543101400002411
Figure BDA00023543101400002412
for supplying heat to nodes upsilon in heat supply network in comprehensive energy systemaTotal number of heat supply pipes for supplying hot working medium, XT0Is a value of the silence state for the temperature,
Figure BDA00023543101400002413
for node gamma in heat supply network in comprehensive energy systemaThe temperature of the mixture is controlled by the temperature,
Figure BDA00023543101400002414
for a node upsilon in a heat supply network in an integrated energy systemaA first coefficient between the temperature and the temperature of the node of the temperature;
determining a first coefficient matrix C of the node temperature between water return pipelines of a heat supply network in an integrated energy system according to the following formular：
Figure BDA0002354310140000247
In the formula (I), the compound is shown in the specification,
Figure BDA00023543101400002415
for node eta in heat supply network in comprehensive energy systemaTemperature and node phiaFirst coefficient between temperatures, phia,ηa∈(1～χN)，χNFor heating network in comprehensive energy systemThe total number of nodes between the middle water return pipelines;
determining a second coefficient matrix b of the node temperature between the water return pipelines of the heat supply network in the comprehensive energy system according to the following formular
Figure BDA0002354310140000251
In the formula (I), the compound is shown in the specification,
Figure BDA0002354310140000258
for node eta in heat supply network in comprehensive energy systemaA second coefficient of temperature;
wherein the hot working medium is supplied from a node phiaThrough a pipeline HaFlow direction node ηaThe method comprises the following steps:
if node etaaNot only by the node phiaThe supplied hot working medium is as follows:
Figure BDA0002354310140000252
Figure BDA0002354310140000253
Figure BDA0002354310140000254
if node etaaOnly by node phiaWhen supplying hot working medium, then:
Figure BDA0002354310140000255
Figure BDA0002354310140000256
in the above formula, the first and second carbon atoms are,
Figure BDA0002354310140000259
for node phi in heat supply network in comprehensive energy systemaThe outlet temperature of (a) is set,
Figure BDA00023543101400002510
for node phi in heat supply network in comprehensive energy systemaThe volume flow rate of (a) is,
Figure BDA00023543101400002511
heat supply pipeline H for heat supply network in comprehensive energy systemaThe composite heat transfer coefficient of (a) is,
Figure BDA00023543101400002512
heat supply pipeline H for heat supply network in comprehensive energy systemaThe length of (a) of (b),
Figure BDA00023543101400002513
for node eta in heat supply network in comprehensive energy systemaA first coefficient between the temperature and its own node temperature,
Figure BDA00023543101400002514
Figure BDA00023543101400002515
for a heat supply network in an integrated energy system to a node etaaTotal number of return pipes supplying hot working medium.
Specifically, the second determining module is configured to:
determining an optimal transformer transformation ratio for the h-th energy supply node of a power supply network in an integrated energy system according to the following formula
Figure BDA00023543101400002516
Figure BDA0002354310140000257
In the formula of Ueh,bFor supplying in an integrated energy systemThe voltage that the h-th energy supply node of the electrical network is connected to by the busbar,
Figure BDA00023543101400002517
the optimal voltage of the h energy supply node of the power supply network in the integrated energy system is h epsilon (1-S)h)，ShIs the total number of nodes in the power supply network in the integrated energy system.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.

Claims (18)

1. Watch and
Figure FDA0002354310130000019
the optimal scheduling method of the damaged comprehensive energy system is characterized by comprising the following steps:
according to the integrated energy system
Figure FDA00023543101300000110
Determining the optimal voltage of each energy supply node of a power supply network in the comprehensive energy system by loss;
determining the optimal transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system according to the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system;
and adjusting the transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system to be the optimal transformer transformation ratio.
2. The method of claim 1, wherein the method is based on an integrated energy system
Figure FDA00023543101300000111
The loss amount is determined to be the optimal voltage of each energy supply node of the power supply network in the integrated energy system, and the method comprises the following steps:
by integrated energy systems
Figure FDA00023543101300000112
Establishing an objective function of an optimized dispatching model of the comprehensive energy system for the purpose of minimum loss;
and solving the objective function of the comprehensive energy system optimization scheduling model based on the constraint condition corresponding to the objective function of the comprehensive energy system optimization scheduling model, and obtaining the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system.
3. The method of claim 2, wherein the objective function of the integrated energy system optimization scheduling model is determined as follows:
Figure FDA0002354310130000011
in the formula, ClossFor integrated energy systems
Figure FDA00023543101300000113
Loss, c1Is electricity
Figure FDA00023543101300000114
Is/are as follows
Figure FDA00023543101300000115
Mass coefficient, c2Is heat
Figure FDA00023543101300000116
Is/are as follows
Figure FDA00023543101300000117
Mass coefficient, w is pressure
Figure FDA00023543101300000118
And electricity
Figure FDA00023543101300000119
The conversion coefficient of (a) is,
Figure FDA0002354310130000012
for electricity of i-th transmission line of power supply network in integrated energy system
Figure FDA00023543101300000120
The loss amount of the waste water is reduced,
Figure FDA0002354310130000013
for pressing of kth heat supply pipeline of heat supply network in comprehensive energy system
Figure FDA00023543101300000121
The loss amount of the waste water is reduced,
Figure FDA0002354310130000014
for heating of kth heat supply pipeline of heat supply network in integrated energy system
Figure FDA00023543101300000122
Loss, k is equal to [ 1-S ∈ ]hl]，ShlFor the number of heat supply pipelines of a heat supply network in an integrated energy system, i belongs to [ 1-S ]el]，SelThe number of transmission lines of the power supply network in the comprehensive energy system.
4. A method according to claim 3, characterized in that the electricity of the i-th transmission line of the supply network in the integrated energy system is determined according to the following equation
Figure FDA00023543101300000123
Loss of volume
Figure FDA0002354310130000015
Figure FDA0002354310130000016
In the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000017
for the electricity flowing from the x-th transmission line connected with the head end node of the i-th transmission line of the power supply network in the integrated energy system into the i-th transmission line of the power supply network in the integrated energy system
Figure FDA00023543101300000124
，
Figure FDA0002354310130000018
For electricity flowing from the i-th transmission line of the power supply network in the integrated energy system to the c-th transmission line connected with the tail end node thereof
Figure FDA00023543101300000125
，RiResistance, P, of the i-th transmission line of a power supply network in an integrated energy systemi,sFor the electric power flowing into the head-end node of the i-th transmission line of the supply network in the integrated energy system, Ui,sThe voltage of the head end node of the ith transmission line of a power supply network in the comprehensive energy system is x ∈ [ 1-S ]x]，SxC belongs to [ 1-S ] is the sum of the transmission lines connected with the head end node of the ith transmission line of the power supply network in the integrated energy systemc]，ScThe total sum of the transmission lines connected with the terminal node of the ith transmission line of the power supply network in the comprehensive energy system;
determining the pressure of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formula
Figure FDA0002354310130000029
Loss of volume
Figure FDA0002354310130000021
Figure FDA0002354310130000022
In the formula, RkIs the fluid flow resistance, x, of the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe volume flow of the kth heat supply pipeline of a heat supply network in the comprehensive energy system;
determining heat of kth heat supply pipeline of heat supply network in comprehensive energy system according to following formula
Figure FDA00023543101300000210
Loss of volume
Figure FDA0002354310130000023
Figure FDA0002354310130000024
In the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000025
the heat of the kth heat supply pipeline flowing into the kth heat supply pipeline of the heat supply network in the integrated energy system from the kth heat supply pipeline connected with the head end node of the kth heat supply pipeline of the heat supply network in the integrated energy system
Figure FDA00023543101300000211
，
Figure FDA0002354310130000026
The heat flowing from the kth heat supply pipeline of the heat supply network in the integrated energy system to the kth heat supply pipeline connected with the tail end node of the heat supply network
Figure FDA00023543101300000212
Rho is the density of the heat supply pipeline transmission working medium, caThe specific heat capacity of the working medium is transmitted for the heat supply pipeline,
Figure FDA0002354310130000027
for the temperature, X, of the head end node of the kth heat supply pipeline of a heat supply network in an integrated energy systemT0For the value of the temperature silence, epsilon, of the heating network in the integrated energy systemkThe composite heat transfer coefficient of the kth heat supply pipeline of a heat supply network in an integrated energy system is shown, wherein z belongs to [ 1-S ]z]，SzThe number of heat supply pipelines connected with the head end node of the kth heat supply pipeline of the heat supply network in the integrated energy system is b belongs to [ 1-S ]b]，SbThe number of the heat supply pipelines connected with the terminal nodes of the kth heat supply pipeline of the heat supply network in the comprehensive energy system.
5. The method of claim 4, wherein determining electricity is performed as follows
Figure FDA00023543101300000213
Is/are as follows
Figure FDA00023543101300000214
Mass coefficient c1：
Figure FDA0002354310130000028
In the formula, phi1As a factor for evaluating the energy level of electric energy, μ1Is an energy level factor of the electrical energy;
heat was determined as follows
Figure FDA00023543101300000215
Is/are as follows
Figure FDA00023543101300000216
Mass coefficient c2：
Figure FDA0002354310130000031
In the formula, phi2Is an energy level evaluation factor of thermal energy, mu2An energy level factor that is thermal energy;
determining the resistance R of the ith transmission line of a power supply network in an integrated energy system according to the following formulai：
Figure FDA0002354310130000032
In the formula, LiLength of i-th transmission line of power supply network in integrated energy system, AiThe transmission channel area K of the ith transmission line of a power supply network in an integrated energy systemiThe resistivity of a conductor material of the ith transmission line of a power supply network in the comprehensive energy system;
determining the fluid flow resistance R of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formulak：
Figure FDA0002354310130000033
In the formula, mukFluid viscosity, L, for the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe length r of the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe inner radius of the kth heat supply pipeline of a heat supply network in the comprehensive energy system;
determining the volume flow x of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formulak：
Figure FDA0002354310130000034
In the formula, pk,sThe pressure intensity of the head end node of the kth heat supply pipeline of the heat supply network in the comprehensive energy system,pk,mthe pressure of the tail end node of the kth heat supply pipeline of the heat supply network in the comprehensive energy system is obtained;
determining the composite heat transfer coefficient epsilon of the kth heat supply pipeline of the heat supply network in the comprehensive energy system according to the following formulak：
Figure FDA0002354310130000035
In the formula (d)k(σ+1)Radius of the material of the sigma +1 layer of the kth heat supply pipeline of the heat supply network in the integrated energy system, dkσRadius, gamma, of the material of the sigma-th layer of the kth heat supply pipeline of a heat supply network in an integrated energy systemkσThe heat transfer coefficient of the material of the sigma layer of the kth heat supply pipeline of the heat supply network in the comprehensive energy system is sigma e (1-psi-1), and psi is the layer number of the heat transfer material of the kth heat supply pipeline of the heat supply network in the comprehensive energy system.
6. The method of claim 3, wherein the constraints of the objective function of the integrated energy system optimization scheduling model comprise: equality constraint conditions, inequality constraint conditions of running output of a generator set, a gas turbine set and a circulating water pump in a power supply network, tolerance constraint conditions, transmission capacity constraint conditions, transformer tap gear constraint conditions and pipeline temperature constraint conditions.
7. The method of claim 6, wherein the equality constraint is determined as follows:
Figure FDA0002354310130000041
in the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000042
for the input electrical power of the q-th supply node of the supply network in the integrated energy system,
Figure FDA0002354310130000043
for the electrical demand of the f-th load node of the power supply network in the integrated energy system,
Figure FDA0002354310130000044
for the current of the ith transmission line of the power supply network in the integrated energy system,
Figure FDA0002354310130000045
for the voltage variation of the ith transmission line of the power supply network in the integrated energy system,
Figure FDA0002354310130000046
for the input thermal power of the d-th heating node of the heating network in the integrated energy system,
Figure FDA0002354310130000047
for the thermal demand of the r-th load node of the heating network in the integrated energy system,
Figure FDA0002354310130000048
is the entropy of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure FDA0002354310130000049
the temperature variation of the kth heat supply pipeline of a heat supply network in an integrated energy system is shown, and q belongs to [ 1-S ]q]，SqIs the total number of power supply nodes of a power supply network in an integrated energy system, f belongs to [ 1-S ]f]，SfThe total number of load nodes of a power supply network in an integrated energy system is d ∈ 1-Sd]，SdThe total number of heat supply nodes of a heat supply network in the comprehensive energy system is r belongs to [ 1-S ]r]The total number of load nodes of a heat supply network in the comprehensive energy system;
determining the pipeline temperature constraint as follows:
Figure FDA00023543101300000410
in the formula, CsA first coefficient matrix, C, for the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy systemrIs a first coefficient matrix, X, of the node temperature between the water return pipelines of the heat supply network in the integrated energy systemTsIs a temperature matrix of the load heat flow inlet, XTrIs a temperature matrix of the load heat flux outlet, XT0A matrix of silence values for temperature, bsA second coefficient matrix of the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy system, brAnd the second coefficient matrix is the temperature of the nodes between the water return pipelines of the heat supply network in the comprehensive energy system.
8. The method of claim 7, wherein the first coefficient matrix C for the node temperatures between heating pipes of the heating network in the integrated energy system is determined as followss：
Figure FDA0002354310130000051
In the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000052
for a node upsilon in a heat supply network in an integrated energy systemaTemperature and node gammaaA first coefficient between temperatures, va，γa∈(γN)，γNThe method comprises the steps of collecting nodes among heat supply pipelines in a heat supply network in an integrated energy system;
second coefficient matrix b for determining temperature of nodes between heat supply pipelines of heat supply network in comprehensive energy system according to following formulas：
Figure FDA0002354310130000053
In the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000054
for a node upsilon in a heat supply network in an integrated energy systemaA second coefficient of temperature, T being the transposed sign;
wherein the hot working medium is supplied from the node gammaaThrough the pipe MaFlow direction node vaThe method comprises the following steps:
if node vaIs not only processed by node gammaaSupplied thermal medium and node gammaaIs a load node, then:
Figure FDA0002354310130000055
Figure FDA0002354310130000056
if node vaIs not only processed by node gammaaSupplied thermal medium and node gammaaIf the node is not a load node, then:
Figure FDA0002354310130000057
Figure FDA0002354310130000058
Figure FDA0002354310130000059
if node vaIs only composed of node gammaaA hot working medium is supplied to
Figure FDA0002354310130000061
Figure FDA0002354310130000062
In the above formula, the first and second carbon atoms are,
Figure FDA0002354310130000063
heat supply pipeline M for heat supply network in comprehensive energy systemaThe volume flow rate of (a) is,
Figure FDA0002354310130000064
heat supply pipeline M for heat supply network in comprehensive energy systemaThe composite heat transfer coefficient of (a) is,
Figure FDA0002354310130000065
heat supply pipeline M for heat supply network in comprehensive energy systemaThe length of (a) of (b),
Figure FDA0002354310130000066
Figure FDA0002354310130000067
for supplying heat to nodes upsilon in heat supply network in comprehensive energy systemaTotal number of heat supply pipes for supplying hot working medium, XT0Is a value of the silence state for the temperature,
Figure FDA0002354310130000068
for node gamma in heat supply network in comprehensive energy systemaThe temperature of the mixture is controlled by the temperature,
Figure FDA0002354310130000069
for a node upsilon in a heat supply network in an integrated energy systemaA first coefficient between the temperature and the temperature of the node of the temperature;
determining a first coefficient matrix C of the node temperature between water return pipelines of a heat supply network in an integrated energy system according to the following formular：
Figure FDA00023543101300000610
In the formula (I), the compound is shown in the specification,
Figure FDA00023543101300000611
for node eta in heat supply network in comprehensive energy systemaTemperature and node phiaFirst coefficient between temperatures, phia,ηa∈(1～χN)，χNThe total number of nodes between water return pipelines in a heat supply network in the comprehensive energy system;
determining a second coefficient matrix b of the node temperature between the water return pipelines of the heat supply network in the comprehensive energy system according to the following formular
Figure FDA00023543101300000612
In the formula (I), the compound is shown in the specification,
Figure FDA00023543101300000613
for node eta in heat supply network in comprehensive energy systemaA second coefficient of temperature;
wherein the hot working medium is supplied from a node phiaThrough a pipeline HaFlow direction node ηaThe method comprises the following steps:
if node etaaNot only by the node phiaThe supplied hot working medium is as follows:
Figure FDA00023543101300000614
Figure FDA0002354310130000071
Figure FDA0002354310130000072
if node etaaOnly by node phiaWhen supplying hot working medium, then:
Figure FDA0002354310130000073
Figure FDA0002354310130000074
in the above formula, the first and second carbon atoms are,
Figure FDA0002354310130000075
for node phi in heat supply network in comprehensive energy systemaThe outlet temperature of (a) is set,
Figure FDA0002354310130000076
for node phi in heat supply network in comprehensive energy systemaThe volume flow rate of (a) is,
Figure FDA0002354310130000077
heat supply pipeline H for heat supply network in comprehensive energy systemaThe composite heat transfer coefficient of (a) is,
Figure FDA0002354310130000078
heat supply pipeline H for heat supply network in comprehensive energy systemaThe length of (a) of (b),
Figure FDA0002354310130000079
for node eta in heat supply network in comprehensive energy systemaA first coefficient between the temperature and its own node temperature,
Figure FDA00023543101300000710
Figure FDA00023543101300000711
for a heat supply network in an integrated energy system to a node etaaTotal number of return pipes supplying hot working medium.
9. The method of claim 1, wherein determining the optimal transformer transformation ratio for each energy supply node of the power supply network in the integrated energy system based on the optimal voltage for each energy supply node of the power supply network in the integrated energy system comprises:
determining an optimal transformer transformation ratio for the h-th energy supply node of a power supply network in an integrated energy system according to the following formula
Figure FDA00023543101300000712
Figure FDA00023543101300000713
In the formula of Ueh,bThe voltage which is accessed by the bus for the h-th energy supply node of the power supply network in the integrated energy system,
Figure FDA00023543101300000714
the optimal voltage of the h energy supply node of the power supply network in the integrated energy system is h epsilon (1-S)h)，ShIs the total number of nodes in the power supply network in the integrated energy system.
10. Watch and
Figure FDA00023543101300000715
comprehensive energy system optimization scheduling system that decreases, its characterized in that, the system includes:
a first determination module for determining the system based on the integrated energy system
Figure FDA00023543101300000716
Determining the maximum of the energy supply nodes of an energy supply network in an integrated energy systemOptimizing voltage;
the second determination module is used for determining the optimal transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system according to the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system;
and the adjusting module is used for adjusting the transformer transformation ratio of each energy supply node of the power supply network in the comprehensive energy system to the optimal transformer transformation ratio.
11. The system of claim 10, wherein the first determination module comprises:
building units for integrated energy systems
Figure FDA0002354310130000089
Establishing an objective function of an optimized dispatching model of the comprehensive energy system for the purpose of minimum loss;
and the obtaining unit is used for solving the objective function of the comprehensive energy system optimization scheduling model based on the constraint condition corresponding to the objective function of the comprehensive energy system optimization scheduling model, and obtaining the optimal voltage of each energy supply node of the power supply network in the comprehensive energy system.
12. The system of claim 11, wherein the objective function of the integrated energy system optimization scheduling model is determined as follows:
Figure FDA0002354310130000081
in the formula, ClossFor integrated energy systems
Figure FDA00023543101300000810
Loss, c1Is electricity
Figure FDA00023543101300000811
Is/are as follows
Figure FDA00023543101300000812
Mass coefficient, c2Is heat
Figure FDA00023543101300000813
Is/are as follows
Figure FDA00023543101300000814
Mass coefficient, w is pressure
Figure FDA00023543101300000815
And electricity
Figure FDA00023543101300000816
The conversion coefficient of (a) is,
Figure FDA0002354310130000082
for electricity of i-th transmission line of power supply network in integrated energy system
Figure FDA00023543101300000817
The loss amount of the waste water is reduced,
Figure FDA0002354310130000083
for pressing of kth heat supply pipeline of heat supply network in comprehensive energy system
Figure FDA00023543101300000818
The loss amount of the waste water is reduced,
Figure FDA0002354310130000084
for heating of kth heat supply pipeline of heat supply network in integrated energy system
Figure FDA00023543101300000819
Loss, k is equal to [ 1-S ∈ ]hl]，ShlFor the number of heat supply pipelines of a heat supply network in an integrated energy system, i belongs to [ 1-S ]el]，SelThe number of transmission lines of the power supply network in the comprehensive energy system.
13. The system of claim 12, wherein the electricity for the ith transmission line of the power supply network in the integrated energy system is determined according to the following equation
Figure FDA00023543101300000820
Loss of volume
Figure FDA0002354310130000085
Figure FDA0002354310130000086
In the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000087
for the electricity flowing from the x-th transmission line connected with the head end node of the i-th transmission line of the power supply network in the integrated energy system into the i-th transmission line of the power supply network in the integrated energy system
Figure FDA00023543101300000821
，
Figure FDA0002354310130000088
For electricity flowing from the i-th transmission line of the power supply network in the integrated energy system to the c-th transmission line connected with the tail end node thereof
Figure FDA00023543101300000822
，RiResistance, P, of the i-th transmission line of a power supply network in an integrated energy systemi,sFor the electric power flowing into the head-end node of the i-th transmission line of the supply network in the integrated energy system, Ui,sThe voltage of the head end node of the ith transmission line of a power supply network in the comprehensive energy system is x ∈ [ 1-S ]x]，SxFor transmission line connected to head end node of ith transmission line of power supply network in integrated energy systemAnd c is from [1 to S ∈c]，ScThe total sum of the transmission lines connected with the terminal node of the ith transmission line of the power supply network in the comprehensive energy system;
determining the pressure of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formula
Figure FDA00023543101300000910
Loss of volume
Figure FDA0002354310130000091
Figure FDA0002354310130000092
In the formula, RkIs the fluid flow resistance, x, of the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe volume flow of the kth heat supply pipeline of a heat supply network in the comprehensive energy system;
determining heat of kth heat supply pipeline of heat supply network in comprehensive energy system according to following formula
Figure FDA00023543101300000911
Loss of volume
Figure FDA0002354310130000093
Figure FDA0002354310130000094
In the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000095
the heat of the kth heat supply pipeline flowing into the kth heat supply pipeline of the heat supply network in the integrated energy system from the kth heat supply pipeline connected with the head end node of the kth heat supply pipeline of the heat supply network in the integrated energy system
Figure FDA00023543101300000912
，
Figure FDA0002354310130000096
The heat flowing from the kth heat supply pipeline of the heat supply network in the integrated energy system to the kth heat supply pipeline connected with the tail end node of the heat supply network
Figure FDA00023543101300000913
Rho is the density of the heat supply pipeline transmission working medium, caThe specific heat capacity of the working medium is transmitted for the heat supply pipeline,
Figure FDA0002354310130000097
for the temperature, X, of the head end node of the kth heat supply pipeline of a heat supply network in an integrated energy systemT0For the value of the temperature silence, epsilon, of the heating network in the integrated energy systemkThe composite heat transfer coefficient of the kth heat supply pipeline of a heat supply network in an integrated energy system is shown, wherein z belongs to [ 1-S ]z]，SzThe number of heat supply pipelines connected with the head end node of the kth heat supply pipeline of the heat supply network in the integrated energy system is b belongs to [ 1-S ]b]，SbThe number of the heat supply pipelines connected with the terminal nodes of the kth heat supply pipeline of the heat supply network in the comprehensive energy system.
14. The system of claim 13, wherein the electricity is determined as follows
Figure FDA00023543101300000914
Is/are as follows
Figure FDA00023543101300000915
Mass coefficient c1：
Figure FDA0002354310130000098
In the formula, phi1As a factor for evaluating the energy level of electric energy, μ1As energy of electric energyA level factor;
heat was determined as follows
Figure FDA00023543101300000916
Is/are as follows
Figure FDA00023543101300000917
Mass coefficient c2：
Figure FDA0002354310130000099
In the formula, phi2Is an energy level evaluation factor of thermal energy, mu2An energy level factor that is thermal energy;
determining the resistance R of the ith transmission line of a power supply network in an integrated energy system according to the following formulai：
Figure FDA0002354310130000101
In the formula, LiLength of i-th transmission line of power supply network in integrated energy system, AiThe transmission channel area K of the ith transmission line of a power supply network in an integrated energy systemiThe resistivity of a conductor material of the ith transmission line of a power supply network in the comprehensive energy system;
determining the fluid flow resistance R of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formulak：
Figure FDA0002354310130000102
In the formula, mukFluid viscosity, L, for the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe length r of the kth heat supply pipeline of a heat supply network in an integrated energy systemkThe inner radius of the kth heat supply pipeline of a heat supply network in the comprehensive energy system;
determining the volume flow x of the kth heat supply pipeline of a heat supply network in an integrated energy system according to the following formulak：
Figure FDA0002354310130000103
In the formula, pk,sPressure intensity p of head end node of kth heat supply pipeline of heat supply network in integrated energy systemk,mThe pressure of the tail end node of the kth heat supply pipeline of the heat supply network in the comprehensive energy system is obtained;
determining the composite heat transfer coefficient epsilon of the kth heat supply pipeline of the heat supply network in the comprehensive energy system according to the following formulak：
Figure FDA0002354310130000104
In the formula (d)k(σ+1)Radius of the material of the sigma +1 layer of the kth heat supply pipeline of the heat supply network in the integrated energy system, dkσRadius, gamma, of the material of the sigma-th layer of the kth heat supply pipeline of a heat supply network in an integrated energy systemkσThe heat transfer coefficient of the material of the sigma layer of the kth heat supply pipeline of the heat supply network in the comprehensive energy system is sigma e (1-psi-1), and psi is the layer number of the heat transfer material of the kth heat supply pipeline of the heat supply network in the comprehensive energy system.
15. The system of claim 12, wherein the constraints of the objective function of the integrated energy system optimization scheduling model include: equality constraint conditions, inequality constraint conditions of running output of a generator set, a gas turbine set and a circulating water pump in a power supply network, tolerance constraint conditions, transmission capacity constraint conditions, transformer tap gear constraint conditions and pipeline temperature constraint conditions.
16. The system of claim 15, wherein the equality constraint is determined as follows:
Figure FDA0002354310130000111
in the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000112
for the input electrical power of the q-th supply node of the supply network in the integrated energy system,
Figure FDA0002354310130000113
for the electrical demand of the f-th load node of the power supply network in the integrated energy system,
Figure FDA0002354310130000114
for the current of the ith transmission line of the power supply network in the integrated energy system,
Figure FDA0002354310130000115
for the voltage variation of the ith transmission line of the power supply network in the integrated energy system,
Figure FDA0002354310130000116
for the input thermal power of the d-th heating node of the heating network in the integrated energy system,
Figure FDA0002354310130000117
for the thermal demand of the r-th load node of the heating network in the integrated energy system,
Figure FDA0002354310130000118
is the entropy of the kth heat supply pipeline of a heat supply network in the comprehensive energy system,
Figure FDA0002354310130000119
the temperature variation of the kth heat supply pipeline of a heat supply network in an integrated energy system is shown, and q belongs to [ 1-S ]q]，SqFor supplying power in an integrated energy systemThe total number of power supply nodes of the network, f belongs to [ 1-S ]f]，SfThe total number of load nodes of a power supply network in an integrated energy system is d ∈ 1-Sd]，SdThe total number of heat supply nodes of a heat supply network in the comprehensive energy system is r belongs to [ 1-S ]r]The total number of load nodes of a heat supply network in the comprehensive energy system;
determining the pipeline temperature constraint as follows:
Figure FDA00023543101300001110
in the formula, CsA first coefficient matrix, C, for the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy systemrIs a first coefficient matrix, X, of the node temperature between the water return pipelines of the heat supply network in the integrated energy systemTsIs a temperature matrix of the load heat flow inlet, XTrIs a temperature matrix of the load heat flux outlet, XT0A matrix of silence values for temperature, bsA second coefficient matrix of the temperature of the nodes between the heat supply pipelines of the heat supply network in the integrated energy system, brAnd the second coefficient matrix is the temperature of the nodes between the water return pipelines of the heat supply network in the comprehensive energy system.
17. The system of claim 16, wherein the first coefficient matrix C for determining the temperature of the nodes between the heating pipes of the heating network in the integrated energy system is determined according tos：
Figure FDA0002354310130000121
In the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000122
for a node upsilon in a heat supply network in an integrated energy systemaTemperature and node gammaaA first coefficient between temperatures, va，γa∈(γN)，γNFor a collection of nodes between heat supply pipelines in a heat supply network in an integrated energy system, CsInitially, all elements in the matrix are 0;
second coefficient matrix b for determining temperature of nodes between heat supply pipelines of heat supply network in comprehensive energy system according to following formulas：
Figure FDA0002354310130000123
In the formula (I), the compound is shown in the specification,
Figure FDA0002354310130000124
for a node upsilon in a heat supply network in an integrated energy systemaSecond coefficient of temperature, T being transposed sign, bsInitially, all elements in the matrix are 0;
wherein the hot working medium is supplied from the node gammaaThrough the pipe MaFlow direction node vaThe method comprises the following steps:
if node vaIs not only processed by node gammaaSupplied thermal medium and node gammaaIs a load node, then:
Figure FDA0002354310130000125
Figure FDA0002354310130000126
if node vaIs not only processed by node gammaaSupplied thermal medium and node gammaaIf the node is not a load node, then:
Figure FDA0002354310130000127
Figure FDA0002354310130000128
Figure FDA0002354310130000129
if node vaIs only composed of node gammaaA hot working medium is supplied to
Figure FDA0002354310130000131
Figure FDA0002354310130000132
In the above formula, the first and second carbon atoms are,
Figure FDA0002354310130000133
heat supply pipeline M for heat supply network in comprehensive energy systemaThe volume flow rate of (a) is,
Figure FDA0002354310130000134
heat supply pipeline M for heat supply network in comprehensive energy systemaThe composite heat transfer coefficient of (a) is,
Figure FDA0002354310130000135
heat supply pipeline M for heat supply network in comprehensive energy systemaThe length of (a) of (b),
Figure FDA0002354310130000136
Figure FDA0002354310130000137
for supplying heat to nodes upsilon in heat supply network in comprehensive energy systemaTotal number of heat supply pipes for supplying hot working medium, XT0Is a value of the silence state for the temperature,
Figure FDA0002354310130000138
for node gamma in heat supply network in comprehensive energy systemaThe temperature of the mixture is controlled by the temperature,
Figure FDA0002354310130000139
for a node upsilon in a heat supply network in an integrated energy systemaA first coefficient between the temperature and the temperature of the node of the temperature;
determining a first coefficient matrix C of the node temperature between water return pipelines of a heat supply network in an integrated energy system according to the following formular：
Figure FDA00023543101300001310
In the formula (I), the compound is shown in the specification,
Figure FDA00023543101300001311
for node eta in heat supply network in comprehensive energy systemaTemperature and node phiaFirst coefficient between temperatures, phia,ηa∈(1～χN)，χNIs the total number of nodes C between water return pipelines in a heat supply network in an integrated energy systemrInitially, all elements in the matrix are 0;
determining a second coefficient matrix b of the node temperature between the water return pipelines of the heat supply network in the comprehensive energy system according to the following formular
Figure FDA00023543101300001312
In the formula (I), the compound is shown in the specification,
Figure FDA00023543101300001313
for node eta in heat supply network in comprehensive energy systemaSecond coefficient of temperature, brInitially, all elements in the matrix are 0;
wherein the hot working medium is supplied from a node phiaThrough a pipeline HaFlow direction node ηaThe method comprises the following steps:
if node etaaNot only by the node phiaThe supplied hot working medium is as follows:
Figure FDA00023543101300001314
Figure FDA0002354310130000141
Figure FDA0002354310130000142
if node etaaOnly by node phiaWhen supplying hot working medium, then:
Figure FDA0002354310130000143
Figure FDA0002354310130000144
in the above formula, the first and second carbon atoms are,
Figure FDA0002354310130000145
for node phi in heat supply network in comprehensive energy systemaThe outlet temperature of (a) is set,
Figure FDA0002354310130000146
for node phi in heat supply network in comprehensive energy systemaThe volume flow rate of (a) is,
Figure FDA0002354310130000147
heat supply pipeline H for heat supply network in comprehensive energy systemaThe composite heat transfer coefficient of (a) is,
Figure FDA0002354310130000148
heat supply pipeline H for heat supply network in comprehensive energy systemaThe length of (a) of (b),
Figure FDA0002354310130000149
for node eta in heat supply network in comprehensive energy systemaA first coefficient between the temperature and its own node temperature,
Figure FDA00023543101300001410
Figure FDA00023543101300001411
for a heat supply network in an integrated energy system to a node etaaTotal number of return pipes supplying hot working medium.
18. The system of claim 10, wherein the root second determination module is to:
determining an optimal transformer transformation ratio for the h-th energy supply node of a power supply network in an integrated energy system according to the following formula
Figure FDA00023543101300001412
Figure FDA00023543101300001413
In the formula of Ueh,bThe voltage which is accessed by the bus for the h-th energy supply node of the power supply network in the integrated energy system,
Figure FDA00023543101300001414
the optimal voltage of the h energy supply node of the power supply network in the integrated energy system is h epsilon (1-S)h)，ShIs the total number of nodes in the power supply network in the integrated energy system.
CN202010003464.3A 2020-01-02 2020-01-02 Comprehensive energy system optimization scheduling method and system considering exergy loss Pending CN113065729A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010003464.3A CN113065729A (en) 2020-01-02 2020-01-02 Comprehensive energy system optimization scheduling method and system considering exergy loss

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010003464.3A CN113065729A (en) 2020-01-02 2020-01-02 Comprehensive energy system optimization scheduling method and system considering exergy loss

Publications (1)

Publication Number Publication Date
CN113065729A true CN113065729A (en) 2021-07-02

Family

ID=76558346

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010003464.3A Pending CN113065729A (en) 2020-01-02 2020-01-02 Comprehensive energy system optimization scheduling method and system considering exergy loss

Country Status (1)

Country Link
CN (1) CN113065729A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113821004A (en) * 2021-08-23 2021-12-21 南方电网科学研究院有限责任公司 Optimization method, device and equipment for building energy management
CN114970964A (en) * 2022-04-24 2022-08-30 国网福建省电力有限公司 Electric heating and gas comprehensive energy system optimal scheduling method based on universal exergy loss model

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113821004A (en) * 2021-08-23 2021-12-21 南方电网科学研究院有限责任公司 Optimization method, device and equipment for building energy management
CN114970964A (en) * 2022-04-24 2022-08-30 国网福建省电力有限公司 Electric heating and gas comprehensive energy system optimal scheduling method based on universal exergy loss model

Similar Documents

Publication Publication Date Title
CN108921727B (en) Regional comprehensive energy system reliability assessment method considering thermal load dynamic characteristics
CN109063890B (en) Thermal load distribution method based on thermal power plant whole-plant peak shaving capacity maximization
CN111061986B (en) Thermoelectric comprehensive energy system tide calculation method with multiple operation modes
CN107730129A (en) Consider the electrical heat interacted system methods of risk assessment of photo-thermal cogeneration of heat and power and electric boiler
Teng et al. Optimal operation strategy for combined heat and power system based on solid electric thermal storage boiler and thermal inertia
CN113065729A (en) Comprehensive energy system optimization scheduling method and system considering exergy loss
CN111832138B (en) Regional pipe network topology optimization method
CN111555285B (en) Energy flow decoupling analysis and calculation method for distributed combined cooling heating and power comprehensive energy system
CN116341881B (en) Robust advanced scheduling method and system for electric-thermal system considering flexibility of heat supply network
CN112241923B (en) Distribution network power balance method based on comprehensive energy system source load equivalent external characteristics
CN113240204A (en) Energy station capacity optimal configuration method and system considering renewable energy consumption area
CN112150024B (en) Multi-scene energy efficiency evaluation method for comprehensive energy system
CN113036819A (en) Electric heating comprehensive energy system robust scheduling method considering source-load bilateral uncertainty
CN114077934A (en) Comprehensive energy microgrid interconnection system and scheduling method thereof
CN115099007A (en) Comprehensive energy system optimized operation method based on comprehensive cost-energy consumption curve
CN113469464B (en) Regional comprehensive energy distribution optimization method considering different heating modes of building
CN113139675B (en) Comprehensive electrical load interval prediction method for park comprehensive energy system
CN105373638A (en) Energy consumption evaluation method and electronic device
CN114662938A (en) Energy efficiency evaluation method and device of comprehensive energy system, terminal and storage medium
CN113761727A (en) Method for constructing optimal scheduling model of thermoelectric combined system with distributed electric heat pumps
CN113298164A (en) Building energy utilization optimization method based on photovoltaic output uncertainty
CN112070274A (en) Efficiency evaluation method and system of comprehensive energy system
CN113077125B (en) Energy efficiency-considered integrated energy system typical scene generation method
CN114675690B (en) Temperature regulation and control method and device for heat exchanger equipment
Rusovs et al. Hybrid heating system sizing and operation control with Neural Network based simulation

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