CN113328443A - Energy consumption optimization control method for side of transformer area - Google Patents

Energy consumption optimization control method for side of transformer area Download PDF

Info

Publication number
CN113328443A
CN113328443A CN202110672510.3A CN202110672510A CN113328443A CN 113328443 A CN113328443 A CN 113328443A CN 202110672510 A CN202110672510 A CN 202110672510A CN 113328443 A CN113328443 A CN 113328443A
Authority
CN
China
Prior art keywords
energy
power
load
heat
constraint
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
CN202110672510.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 Jiangsu Electric Power Co ltd Marketing Service Center
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Jiangsu Electric Power Co ltd Marketing Service Center
State Grid Corp of China SGCC
State Grid Jiangsu 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 Jiangsu Electric Power Co ltd Marketing Service Center, State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Jiangsu Electric Power Co ltd Marketing Service Center
Priority to CN202110672510.3A priority Critical patent/CN113328443A/en
Publication of CN113328443A publication Critical patent/CN113328443A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/26Arrangements for eliminating or reducing asymmetry in polyphase networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/50Arrangements for eliminating or reducing asymmetry in polyphase networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E70/00Other energy conversion or management systems reducing GHG emissions
    • Y02E70/30Systems combining energy storage with energy generation of non-fossil origin
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Human Resources & Organizations (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Water Supply & Treatment (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Feedback Control In General (AREA)

Abstract

A platform side energy consumption optimization control method is characterized in that a multi-energy system model is established based on an energy hub; the method comprises the steps that the optimal 24-hour running cost of the multi-energy system is taken as a target function, and an energy utilization optimization model of the multi-energy system is constructed by combining the constraint of the running limit of the multi-energy system and the constraint of a running mechanism; solving an energy utilization optimization model of the multi-energy system based on a particle swarm algorithm to obtain output optimization data of each energy subsystem in the multi-energy system; the comprehensive optimization control method for the energy system with centralized and distributed cooperation and fusion of renewable energy and fossil energy is realized, the comprehensive utilization efficiency of multiple energy sources is improved, and the aims of energy conservation, environmental protection, economic operation and the like are fulfilled.

Description

Energy consumption optimization control method for side of transformer area
Technical Field
The invention relates to the technical field of energy coordination control, in particular to a station area side-oriented energy consumption optimization control method.
Background
Along with the development of economic society, the intimacy of energy and human life is continuously increased, and the scientificity and rationality of energy production, distribution, transportation and use not only concern the energy industry itself, but also have important influence on various aspects such as economic transformation development, environmental protection, social safety and the like. Although the construction of urban energy systems in China currently obtains great achievements, the problems of difficult renewable energy consumption, serious environmental pollution, low overall energy efficiency and the like are still serious. The traditional urban energy supply in China is mainly divided into several items such as electric power, gas and heat supply, and the items are respectively completed by corresponding energy departments. The energy supply mode which is divided into different classes and lacks of overall planning not only causes the waste of resources and facilities but also cannot exert the complementary advantages of the properties of various energy sources because the interconversion property among different energy sources is not considered, thereby influencing the improvement of the utilization efficiency of the whole energy sources.
In the prior art, energy-saving transformation is usually carried out on a single energy-using object, and although a certain energy-saving benefit is obtained, the energy-saving potential generated by the overall optimization of the multi-energy complementary energy network system is far greater than the energy-saving effect of the energy-saving transformation on the single energy-using object. Moreover, in the existing energy-saving control method of the multi-energy complementary energy network, the comprehensive consideration of energy utilization of a plurality of energy utilization and capacity objects is lacked, so that the system cannot be comprehensively optimized; in addition, the complex energy system theory and optimization method system for regional energy planning is weak, so that the corresponding multi-energy combined supply and coordinated interconnection optimization technologies are relatively few, and particularly, the system-oriented test and simulation capability is still to be improved.
Therefore, it is necessary to construct a comprehensive optimization control method for an energy system with centralized and distributed cooperation and fusion of renewable energy and fossil energy.
Disclosure of Invention
In order to solve the defects in the prior art, the invention aims to provide a district-oriented energy consumption optimization control method, realize a comprehensive optimization control method of an energy system integrating centralized and distributed cooperation and renewable energy and fossil energy, improve the comprehensive utilization efficiency of multiple energy sources, and realize the aims of energy conservation, environmental protection, economic operation and the like.
The invention adopts the following technical scheme.
The method for optimizing and controlling energy consumption of the platform area side comprises the following steps:
step 1, based on an energy hub, establishing a multi-energy system model by taking energy input data of each energy subsystem as input data and output data of each energy subsystem as output data;
step 2, based on a multi-energy system model, constructing an energy utilization optimization model of the multi-energy system by taking the optimal 24-hour running cost of the multi-energy system as a target function and combining the constraint of the running limit of the multi-energy system and the constraint of a running mechanism;
wherein the objective function satisfies the following relation:
Figure BDA0003119236720000021
wherein F is the 24-hour running cost of the multi-energy system, QG(t) the natural gas quantity, Q, purchased by the multi-energy system in the tthD(t) the amount of diesel oil purchased by the multi-energy system in the tth hour, CG(t) the purchase price of natural gas at the t hour, CD(t) the purchase price of the diesel oil in the tth hour;
the constraints of the operation limit of the multi-energy system comprise energy supply equipment power constraints and energy storage equipment capacity constraints; constraints on the operating mechanism of the multi-energy system include: electric power balance constraint, thermal power balance constraint, cold power balance constraint and power network power flow constraint;
and 3, solving the energy utilization optimization model of the multi-energy system by taking the energy input data and the operation constraint of various energy sources as input data based on the particle swarm optimization to obtain output optimization data of each energy subsystem in the multi-energy system.
Preferably, the first and second electrodes are formed of a metal,
in the step 1, each energy subsystem comprises a cogeneration unit, a diesel generator, a gas boiler unit, an ice storage system, a hot water storage tank system and a V2G battery replacement station.
In step 1, establishing a multi-energy system model comprises:
step 1.1, building an energy network, wherein the energy network comprises a power network for transmitting electric energy and a fluid network for transmitting heat energy; the power network and the fluid network are coupled by virtue of the pump and the total heat load balance, wherein the total heat load balance meets the following relational expression:
PT=PTh+PTe
in the formula, PThHeat, P, supplied to the user for the fluid networkTeHeat provided to the user for the power network;
step 1.2, based on a network theory, selecting a branch where a pump is located as a tree branch and the rest branches as connecting branches in a power network, and selecting a branch where the pump is located as a tree branch and the rest branches as connecting branches in a fluid network; numbering each branch, each node and each basic loop, and selecting the direction of each branch and each basic loop; each branch is a section of line without branches in the network, each node is an end point of each branch, and each basic loop is a closed loop which comprises and only comprises one continuous branch; the direction of each branch is selected to be the same as the current, and the direction of each basic loop is the clockwise direction;
step 1.3, selecting common nodes of all loops as reference nodes, and establishing a multi-energy system model according to the generalized kirchhoff law by using the following relational expression:
Figure BDA0003119236720000031
in the formula (I), the compound is shown in the specification,
a is an energy network incidence matrix,
h is an extensive flow matrix, the extensive flow matrix is a charge flow matrix in the process of transmitting electric energy by the power network, the extensive flow matrix is an entropy flow matrix in the process of transmitting heat energy by the fluid network,
Bfin the form of a matrix of elementary loops,
and delta X is a matrix of the difference value of the intensities of the two ends of the line, the matrix of the difference value of the intensities of the two ends of the line is a matrix of the potential difference value in the process of transferring electric energy by the power network, and the matrix of the difference value of the intensities of the two ends of the line is a matrix of the temperature difference value of the two ends of the line in the process of transferring heat energy by the fluid network.
Preferably, the first and second electrodes are formed of a metal,
in step 2, the power constraint of the energy supply equipment comprises the following steps: the method comprises the following steps of (1) gas internal combustion engine power constraint, air conditioner power constraint, flue gas absorption type refrigerating unit power constraint, compression type electric refrigerating unit power constraint, gas boiler power constraint and waste heat boiler power constraint;
energy storage device capacity constraints include: the capacity of the storage battery is restricted, the capacity of the heat storage tank is restricted, the capacity of the ice storage is restricted, the capacity of the heat pump is restricted, and the capacity of the water chilling unit is restricted.
In step 2, the electric power balance constraint satisfies the following relational expression:
Figure BDA0003119236720000032
in the formula (I), the compound is shown in the specification,
Pbuyin order to purchase the electric quantity,
Psellin order to sell the electricity quantity,
PPV,kis the power generation amount of the kth photovoltaic cell, wherein NPVRepresents the total number of photovoltaic cells,
PGE,iis the power generation amount of the ith gas internal combustion engine, wherein NGEIndicates the total number of gas-fired internal combustion engines,
PBT,outis the output power of the storage battery,
PBT,inis the input power of the storage battery,
Pacfor the power consumption of the air-conditioning system at the time t, PecFor the power consumption of the compression-type electric refrigerator at the time t, Php_hIs the power consumption, P, of the heat pump for heating at time thp_cIs the power consumption, P, of the heat pump for cooling at the time tloadThe electric load of the system at the t moment;
the thermal power balance constraint comprises a space thermal load thermal power balance constraint and a hot water load thermal power balance constraint, wherein the space thermal load thermal power balance constraint satisfies the following relational expression:
Figure BDA0003119236720000041
the hot water load thermal power balance constraint satisfies the following relational expression:
Figure BDA0003119236720000042
in the formula (I), the compound is shown in the specification,
Hgb,j1for the space heat load generated by the jth gas boiler,
Hgb,j2for the hot water load produced by the jth gas boiler,
Hh,out1is the space heat load generated by the heat storage tank,
Hh,out2is the load of hot water generated by the heat storage tank,
Hhp1for the space heat load generated by the heat pump,
Hhp2for the hot water load generated by the heat pump,
Hb,out1for the space heat load generated by the waste heat boiler,
Hb,out2is the hot water load generated by the waste heat boiler,
Hload1for the space heat load of the multi-energy system at the t moment,
Hload2for the hot water load of the multi-energy system at time t,
Hh,inis the load of hot water consumed by the heat storage tank,
Ngbthe total number of the gas-fired boilers;
the cold power balance constraint comprises a space cold load cold power balance constraint and a cold water load cold power balance constraint, wherein the space cold load cold power balance constraint satisfies the following relational expression:
Qac+Qsar1+Qgar1+Qec1+Qc1+Qhp+Qic,out1=Qload1
the cold water load cold power balance constraint satisfies the following relational expression:
Qsar2+Qgar2+Qec2+Qic,out2+Qc2=Qic,in+Qload2
in the formula (I), the compound is shown in the specification,
Qacfor the space cooling load generated by the air conditioning system,
Qsar1for space cooling loads generated by vapor absorption chiller units,
Qsar2for the cold water load generated by a vapor absorption chiller,
Qgar1is the space cold load generated by the smoke absorption refrigerating unit,
Qgar2is the cold water load generated by the smoke absorption refrigerating unit,
Qec1the space cold load generated by the compression type electric refrigerator,
Qec2is the cold water load produced by the compression type electric refrigerator,
Qc1is the space cold load generated by the water chilling unit,
Qc2is the cold water load generated by the water chilling unit,
Qhpis the space cold load generated by the heat pump,
Qic,out1the space cold load generated for ice cold accumulation,
Qic,out2the cold water load generated for ice storage,
Qic,inthe cold water load consumed for ice storage,
Qload1for the space cooling load of the multi-energy system at the time t,
Qload2the cold water load of the multi-energy system at the t moment;
the power network power flow constraint comprises non-out-of-limit constraint on the interactive power flow of a user and a distribution network and upper and lower limit constraint on the node voltage of the whole network, wherein the non-out-of-limit constraint satisfies the following relational expression:
Figure BDA0003119236720000051
Figure BDA0003119236720000061
in the formula (I), the compound is shown in the specification,
Figure BDA0003119236720000062
for the line load active power of network nodes l to n at time t,
Figure BDA0003119236720000063
for the line-load active power minimum limit from network node l to node n,
Figure BDA0003119236720000064
for the maximum limit of the line load active power of the network node l to the node n,
Figure BDA0003119236720000065
for line-flow reactive power of network nodes l to n at time t,
Figure BDA0003119236720000066
for the line-flow reactive power minimum limit from network node l to node n,
Figure BDA0003119236720000067
for the line flow reactive power maximum limit from network node l to node n,
l,n∈nnode,nnodethe total number of all nodes in the network;
the upper and lower limit constraints satisfy the following relational expression:
Figure BDA0003119236720000068
in the formula (I), the compound is shown in the specification,
Figure BDA0003119236720000069
is the voltage at the node l and,
Figure BDA00031192367200000610
the upper and lower limits of the voltage at node l.
Preferably, the first and second electrodes are formed of a metal,
the step 3 comprises the following steps:
step 3.1, taking energy input data and operation constraints of various energy sources as input data, and taking hourly output data of each energy subsystem as an optimization variable;
step 3.2, defining each optimized variable as a random particle swarm, and initializing;
3.3, taking the random particle swarm as input data of iterative computation, taking a target function as a fitness function, and iteratively updating an individual extreme value and a global extreme value of the random particle swarm according to preset particle swarm parameters to determine the position and the speed of the particles; obtaining a particle swarm optimization solution after multiple iterations; the preset particle swarm parameters comprise the size of the particle swarm and the iteration times;
step 3.4, verifying whether the comprehensive utilization efficiency of the cogeneration unit reaches the lowest or not under the obtained particle swarm optimization solution: if the comprehensive utilization efficiency of the cogeneration unit reaches the minimum standard, the obtained particle swarm optimization solution is output optimization data of each energy subsystem in the multi-energy system every hour; and if the comprehensive utilization efficiency of the cogeneration unit does not reach the minimum standard, returning to the step 3.1, and adding a minimum heat-power ratio constraint condition to adjust the heat-power ratio of the cogeneration unit.
In step 3.4, the comprehensive utilization efficiency of the cogeneration unit meets the following relational expression:
Figure BDA0003119236720000071
in the formula, QeFor electric energy output by cogeneration units, QhHeat energy, Q, output by cogeneration unitsinThe sum of the heat values of diesel oil and natural gas is consumed for the cogeneration unit.
In step 3.4, the thermoelectric ratio is adjusted by the constraint condition of the minimum thermoelectric ratio represented by the following relational expression:
R′CHPmin=RCHPmin+ΔRCHP
in the formula (II), R'CHPminFor adjusted minimum thermoelectric ratio, RCHPminFor minimum thermoelectric ratio before adjustment, Δ RCHPIs the minimum thermoelectric ratio constraint value.
Compared with the prior art, the method has the advantages that the method realizes multi-energy combined supply and coordinated interconnection optimization technology by constructing a complex energy system theory and an optimization method system facing regional energy planning, plays a positive role in innovating an energy service mode, guaranteeing the balance of supply and demand of a power grid, improving the economic operation level of the power grid, improving the overall efficiency of a regional energy system, reducing environmental pollution, developing an intelligent power grid, reducing the total investment and operation management cost of a social energy system and optimizing an energy structure, and provides an effective technical means for greatly reducing the energy consumption in a short period.
Drawings
FIG. 1 is a flow chart of the energy consumption optimization control method for the side of the platform of the invention;
FIG. 2 is a schematic diagram of a simple energy network in the optimization control method for energy consumption at the platform side according to the present invention;
FIG. 3 is a schematic diagram illustrating numbering of branches, nodes and basic loops of an energy network in the method for controlling energy consumption optimization of a platform side;
FIG. 4 is a flow chart of a PSO optimization algorithm in the optimization control method for energy consumption of the platform side;
FIG. 5 is a graph of spring cold, hot, and electrical load curves for a campus oriented to one embodiment of the energy consumption optimization control method for the side of the campus;
FIG. 6 is a graph of the cool, hot and electrical load in summer of a campus of one embodiment of the optimization control method for energy consumption of the power-consumption management system of the invention;
FIG. 7 is a diagram of the cold, hot and electric load curves in autumn of a park according to one embodiment of the energy consumption optimization control method for the side of the park;
FIG. 8 is a graph of winter cold, hot, and electrical loads for a campus according to one embodiment of the optimal control method for energy consumption on the power-consumption side of the campus of the present invention;
FIG. 9 illustrates the natural gas prices for seasonal parks in accordance with one embodiment of the invention directed to a power utility optimization control method;
FIG. 10 is a graph of thermoelectric ratio versus energy conversion efficiency for one embodiment of the energy use optimization control method for the side of the platform of the present invention;
FIG. 11 is a schematic diagram of the output of the ice storage air conditioner in spring and the equivalent electrical load before and after the individual consideration of the cold load according to an embodiment of the optimal control method for energy consumption of the platform area side;
FIG. 12 is a schematic diagram of equivalent electrical loads before and after output of a summer ice storage air conditioner and independent consideration of cold load in an embodiment of an energy consumption optimization control method for a side of a platform area according to the present invention;
FIG. 13 is a schematic diagram of equivalent electrical loads before and after output of an autumn ice storage air conditioner and independent consideration of cold load, according to an embodiment of an energy consumption optimization control method for a platform area side;
FIG. 14 is a schematic diagram of equivalent electrical loads before and after the output of the ice storage air conditioner in winter and the independent consideration of the cold load according to an embodiment of the optimal control method for energy consumption of the side of the platform;
FIG. 15 is a graph of summer diesel engine output before and after consideration of cold load alone for an embodiment of the power consumption optimization control method for the side of the platform of the present invention;
FIG. 16 is a graph comparing the output of the hot water storage tank in spring and the gas price in the embodiment of the optimal control method for energy consumption facing the platform area side;
FIG. 17 is a diagram of the comparison between the output of the heat storage water tank in summer and the gas price in the embodiment of the optimal control method for energy consumption facing the platform area side;
FIG. 18 is a comparison graph of the output of the autumn heat storage water tank and the gas price facing to the platform area side energy consumption optimization control method;
FIG. 19 is a graph comparing the output of the hot water storage tank in winter and the gas price in the embodiment of the optimal control method for the energy consumption of the side of the platform;
FIG. 20 is a graph comparing the spring output and the gas price for an embodiment of the energy consumption optimization control method facing the platform area side of the present invention;
FIG. 21 is a diagram of the comparison of summer output and gas price for an embodiment of the optimal control method for energy consumption at the side of a platform;
FIG. 22 is a diagram of comparing autumn output and gas prices facing to an embodiment of the method for controlling energy consumption optimization of the platform area side;
fig. 23 is a graph comparing the output in winter and the gas price of an embodiment of the optimal energy consumption control method facing the platform area side of the invention.
Detailed Description
The present application is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present application is not limited thereby.
As shown in fig. 1, a method for controlling energy consumption optimization facing a platform area side includes:
step 1, based on an energy hub, taking energy input data of each energy subsystem as input data, and taking output data of each energy subsystem as output data, and establishing a multi-energy system model.
In particular, the amount of the solvent to be used,
in step 1, the energy subsystems include but are not limited to: the system comprises a heat and power cogeneration unit (CHP), a diesel generator, a gas boiler unit, an ice storage system, a hot water storage tank system and a V2G battery replacement station.
In the preferred embodiment, a CHP with a afterburning waste heat boiler is arranged in the multi-energy system facing the platform area side to change the thermoelectric ratio of the energy supply side, and the CHP can collect heat energy generated by the micro-combustion engine in the power generation process by using the waste heat boiler so as to improve the comprehensive energy efficiency; the load loss rate of the system can be obviously reduced by installing the diesel engine power generation system so as to improve the reliability of the system; the gas boiler can be used as auxiliary heat production equipment to work in cooperation with CHP to supply heat load together; the ice storage system can be used to provide cooling load for the park; the high electric energy storage installation cost is considered, so that the capacity configuration of a battery energy storage system is reduced, the excess capacity of a large-scale storage system is realized by taking heat storage as a main energy storage mode and through electric-thermal coupling, the load curve is finally stabilized, and the reliability of the system is improved; the V2G trades power station can utilize electric automobile's energy storage ability and gardens comprehensive energy system to interact to play the effect of buffer system load, with the stability of improvement system.
It is noted that those skilled in the art can select different types of energy subsystems according to engineering practice, and the preferred embodiment of the present invention is a non-limiting preferred choice.
Furthermore, in step 1, the output data of the CHP is the thermoelectric ratio, i.e. the ratio of the thermal power actually output by the CHP to the actual output electric power, and satisfies the following relation:
Figure BDA0003119236720000091
in the formula (I), the compound is shown in the specification,
RCHPis the heat-electricity ratio of the cogeneration unit,
PCHPhis the output thermal power of the cogeneration unit,
PCHPefor the output electric power of the cogeneration unit,
PCHPinis the input power of the cogeneration unit,
ηCHPhfor the heat conversion efficiency of the cogeneration unit,
ηCHPethe electricity conversion efficiency of the cogeneration unit;
in the preferred embodiment, in order to change the output thermoelectric ratio of the CHP, the diesel afterburning amount needs to be increased to improve the heat conversion efficiency of the CHP, which may increase the operation cost of the system, so the operating constraint condition of the cogeneration unit is defined, and the thermoelectric ratio is adjusted by the fuel afterburning amount, and the following relation is satisfied:
RCHP(t)=Ka1(Kb1P0(t)+Kb2)+Ka2
in the formula (I), the compound is shown in the specification,
P0(t) is the amount of diesel afterburning,
Ka1for the first enhancement factor of the thermoelectric ratio,
Ka2for the first enhancement factor of the thermoelectric ratio,
Kb1the first influencing factor is a supplementary firing of the fuel,
Kb2a second influencing factor for fuel oil afterburning;
in the preferred embodiment, Ka1、Ka2、Kb1、Kb2Set at 4.4483, 1.2857, 0.0051209, 0.00038, respectively. Therefore, the mathematical relationship between the post-combustion fuel amount and the heat conversion efficiency can be expressed as:
ηCHPh=ηCHPeRCHP=ηCHPe[Ka1(Kb1P0(t)+Kb2)+Ka2]
from the above formula, the electrical conversion efficiency of the CHP is almost constant, and the heat conversion efficiency of the CHP can be adjusted only by adding the post-combustion fuel oil, so as to adjust the heat-electricity ratio of the CHP.
Further, in step 1, the output data of the diesel generator is output power, and the following relational expression is satisfied:
PDE=ηDEPDEmin
in the formula (I), the compound is shown in the specification,
PDEis the output power of the diesel generator,
ηDEin order to achieve the energy conversion efficiency of the diesel generator,
PDEminthe input power of the diesel generator;
the operation constraint condition of the diesel generator meets the following relational expression:
PDEmin<<PDE<<PDEmax
in the formula (I), the compound is shown in the specification,
PDEminfor starting diesel generatorsThe minimum input power is the minimum input power,
PDEmaxto maintain maximum input power for diesel generator operation.
Further, in step 1, the output data of the gas-fired boiler unit is output power, and the following relational expression is satisfied:
PGB=ηGBPGBmin
in the formula (I), the compound is shown in the specification,
PGBis the output power of the gas boiler unit,
ηGBin order to improve the energy conversion efficiency of the gas boiler unit,
ηGB=(100-a1-a2-a3) %, wherein a1For heat loss of exhaust gas, a2Loss due to incomplete combustion of chemical species, a3Heat dissipation loss; in the preferred embodiment, a1、a2、a3Set to 9.87, 0.95 and 2.18, respectively.
The operation constraint conditions of the gas boiler unit comprise operation capacity constraint and operation time constraint, wherein the operation capacity constraint satisfies the following relational expression:
QGBHmin≤QGBH(t)≤QGBHmax
in the formula (I), the compound is shown in the specification,
QGBHminis the minimum operation capacity of the gas boiler unit,
QGBHmaxthe maximum operation capacity of the gas boiler unit;
the runtime constraint satisfies the following relationship:
TGBON(t)≥TGBONmin
TGBOFF(t)≥TGBOFFmin
in the formula (I), the compound is shown in the specification,
TGBON(t) is the accumulated opening time of the gas boiler unit at the t-th moment,
TGBOFF(t) is the accumulated closing time of the gas boiler unit at the t-th moment,
TGBONminfor the minimum operation time of the gas boiler plant,
TGBOFFminis the minimum off-time of the gas boiler plant.
Further, in step 1, the output data of the ice storage system is refrigeration power, and the following relational expression is satisfied:
PC=PCinEER
in the formula (I), the compound is shown in the specification,
PCis the refrigerating power of the ice cold storage system,
PCinis the input power of the ice storage system,
the EER is the energy efficiency ratio of the ice storage system, and in the preferred embodiment, the energy efficiency ratio is taken to be 3.1.
The ice storage and melting work of the ice storage system is mainly completed in the ice storage tank. The operation constraint condition of the ice cold storage system meets the following relational expression:
Figure BDA0003119236720000121
in the formula (I), the compound is shown in the specification,
QCice(t) is the ice storage capacity of the ice storage tank in the ice storage system,
QCiceminis the minimum ice storage capacity of an ice storage tank in the ice cold storage system,
QCicemaxis the maximum ice storage amount of an ice storage tank in the ice storage system,
Qw(t) is the ice melting amount of an ice storage tank in the ice storage system,
Qwminis the minimum ice melting amount of an ice storage tank in the ice cold storage system,
Qwmaxis the maximum ice melting amount of an ice storage tank in the ice storage system,
MCICE(t) setting an ice making operation mode of the dual-working-condition main machine,
MCwand (t) setting the ice melting operation mode of the dual-working-condition host machine.
Further, in step 1, the output data of the hot water storage tank system is heat energy, and the following relational expression is satisfied:
QWT(t)=(1-εWT)QWT(t-1)+QWTh(t)-QWTr(t)
in the formula (I), the compound is shown in the specification,
QWT(t) is the heat stored by the hot water storage tank system at the time of the t,
εWTfor the self-heat release rate of the hot water storage tank system,
QWTh(t) is the operating capacity of the hot water storage tank system to store heat,
QWTr(t) the heat release of the hot water storage tank system;
the operation constraint condition of the hot water storage tank system meets the following relational expression:
Figure BDA0003119236720000122
in the formula (I), the compound is shown in the specification,
QWThminfor the minimum operating capacity limit of the hot water storage tank system when storing heat,
QWThmaxis the maximum operating capacity limit of the hot water storage tank system when storing heat,
QWTrminis the minimum operation capacity limit value of the hot water storage tank system when releasing heat,
QWTrmaxis the maximum operation capacity limit value of the hot water storage tank system when releasing heat,
QWTminis the minimum limit value of the heat storage quantity of the hot water storage tank system,
QWTmaxthe maximum limit value of the heat storage quantity of the hot water storage tank system.
Further, in step 1, the output data of the V2G power conversion station is the charging and discharging time, and respectively satisfy the following relations:
Figure BDA0003119236720000131
in the formula (I), the compound is shown in the specification,
tscfor the charging start time of the V2G power conversion station,
tc1system negative for V2G power conversion stationAt the time of the valley period of the lotus seeds,
tc2at the end of the system load valley period of the V2G power change station,
tchgfor the continuous charging time of the V2G power conversion station,
k*is a random number between 0 and 1, is a parameter for stabilizing the system load fluctuation to the maximum extent and reducing the peak-valley difference of a power grid of a charging automobile which embodies the internal parameters of the power station and the demand response,
Tgthe valley period duration of the V2G power change station;
Figure BDA0003119236720000132
in the formula (I), the compound is shown in the specification,
tSDfor the starting discharge moment of the V2G power conversion station,
td1at the time of the system load peak section of the V2G power change station,
td2for the end time of the system load peak section of the V2G power change station,
tchdfor the continuous discharge time of the V2G power conversion station, TfThe peak-segment time length of the V2G power change station is obtained;
in addition, considering the load characteristic and the energy storage characteristic of the V2G power conversion station, the operation of the power conversion station is mainly influenced by the power conversion amount LEVAnd the limitation of the state of charge, therefore, the operation constraint condition of the V2G power conversion station meets the following relational expression:
Figure BDA0003119236720000133
in the formula (I), the compound is shown in the specification,
QEV(t) is the energy storage capacity of the V2G power conversion station at the t-th moment,
QEVfor the rated energy storage capacity of the V2G power conversion station,
SOCminin order to be at a minimum state of charge,
SOCmaxat maximum state of charge;
in the present preferred embodiment of the present invention,in order to enable the V2G power station to participate in system peak regulation to the maximum extent, the SOC is setminIs 10%, SOCmaxIs 100%.
The energy hub is an abstraction of a multi-energy system, the multi-energy system is considered to be a dual-port network with multi-form energy input and multi-energy output, and an intermediate energy conversion link is represented by a coupling matrix; the energy hub is utilized to model the multi-energy system, the collaborative optimization of energy sources in various forms can be realized, the established model is connected with different energy carriers to associate different inputs and outputs of the system, so that the supply path of the energy sources is increased, different kinds of energy sources can be coupled, and finally the economy, flexibility, safety and energy efficiency of the energy system are improved.
In particular, the amount of the solvent to be used,
in step 1, establishing a multi-energy system model comprises:
step 1.1, building an energy network, wherein the energy network comprises a power network for transmitting electric energy and a fluid network for transmitting heat energy;
the energy network comprises a plurality of forms of energy transfer, in the preferred embodiment, in a simple energy network involving only electrical and thermal energy, as shown in FIG. 2, the solid lines represent the power network; the dotted lines represent the fluid network.
The power and fluid networks rely on pumps TM and total heat load PTBalanced to achieve coupling, in which the total heat load PTThe balance satisfies the following relation:
PT=PTh+PTe
in the formula, PThHeat, P, supplied to the user for the fluid networkTeHeat provided to the user by the power network.
Step 1.2, based on a network theory, selecting a branch where a pump is located as a tree branch and the rest branches as connecting branches in a power network, and selecting a branch where the pump is located as a tree branch and the rest branches as connecting branches in a fluid network; numbering each branch, each node and each basic loop, and selecting the direction of each branch and each basic loop; each branch is a section of line without branches in the network, each node is an end point of each branch, and each basic loop is a closed loop which comprises and only comprises one continuous branch; the direction of each branch is selected to be the same as the current, and the direction of each basic loop is the clockwise direction;
in the preferred embodiment, the simple energy network shown in fig. 2 is represented by fig. 3, based on network theory.
Step 1.3, selecting common nodes of all loops as reference nodes, and establishing a multi-energy system model according to the generalized kirchhoff law by using the following relational expression:
Figure BDA0003119236720000141
in the formula (I), the compound is shown in the specification,
a is an energy network incidence matrix,
h is an extensive flow matrix, the extensive flow matrix is a charge flow matrix in the process of transmitting electric energy by the power network, the extensive flow matrix is an entropy flow matrix in the process of transmitting heat energy by the fluid network,
Bfin the form of a matrix of elementary loops,
and delta X is a matrix of the difference value of the intensities of the two ends of the line, the matrix of the difference value of the intensities of the two ends of the line is a matrix of the potential difference value in the process of transferring electric energy by the power network, and the matrix of the difference value of the intensities of the two ends of the line is a matrix of the temperature difference value of the two ends of the line in the process of transferring heat energy by the fluid network.
In the preferred embodiment, as shown in FIG. 3, n is selected1、n2For the reference node, the energy network association matrix is as follows:
Figure BDA0003119236720000151
the extensive volume flow matrix is as follows:
H=[H1 H2 H3 H4 H5 H6 H7 H8]T
the basic loop matrix is:
Figure BDA0003119236720000152
the matrix of the intensity difference values at the two ends of the line is as follows:
ΔX=[ΔX1 ΔX2 ΔX3 ΔX4 ΔX5 ΔX6 ΔX7 ΔX8]T
in the electric energy transfer process, the intensity quantity X is electric potential, the intensity quantity difference value is electric potential difference value delta X, and the extensive quantity H is electric charge. In the process of transferring heat energy, the strength X is temperature, the strength difference value is temperature difference value delta X, and the extension H is entropy.
For the model of the multi-energy system in step 1.3, it can be seen from fig. 3 that the power branch b is removed3Branch b in the fluid network where the pump is located5Branch b of the electric network in which the pump is located1And besides the electric heater branch, other branches can write a branch characteristic equation:
Figure BDA0003119236720000153
from fig. 2, TR is the heat exchanger at the heat source to transfer the heat energy to the hot water network, therefore, in the above formula, L is other electrical loads in the power network, R is the flow resistance of the hot water network, and SeFor power supply, each parameter subscript line indicates the number of the branch in which it is located.
Branch b of the fluid network in which the pump is located5And branch b in which the pump is located in the electrical network1The energy conservation relationship therebetween, so the following relationship can be listed:
ΔX5H5+ΔX1H1=0
based on the thermal load characteristics of the electric heater, the following relationships are listed for the heater branches:
P2e=ΔX4H4
in the formula, P2eIs the heat load power of the electric heater.
The preferred embodimentIs the operating state at which the calculation meets the user load demand, and thus the load is known, i.e. PThOr PTAre known. When P is presentThWhen known, the heat exchange satisfies the following relational expression:
PTh=ρcΔTH7
in the formula, rho is the density of the fluid, c is the specific heat capacity of the fluid, and delta T is the temperature difference between the inlet and the outlet of the heat exchanger.
And 2, based on the multi-energy system model, considering the purchase cost of natural gas and diesel oil, taking the optimal 24-hour running cost of the multi-energy system as a target function, and combining the constraint of the running limit of the multi-energy system and the constraint of the running mechanism to construct an energy utilization optimization model of the multi-energy system.
Wherein the objective function satisfies the following relation:
Figure BDA0003119236720000161
in the formula (I), the compound is shown in the specification,
f is the 24-hour running cost of the multi-energy system,
QG(t) is the natural gas quantity purchased by the multi-energy system in the tth hour,
QD(t) is the amount of diesel oil purchased by the multi-energy system in the tth hour,
CG(t) is the purchase price of natural gas at the tth hour,
CD(t) is the purchase price of the diesel oil at the t hour.
Constraints on the operating limits of the multi-energy system include energy supply device power constraints and energy storage device capacity constraints.
Further, the powered device power constraints include:
(1) the gas internal combustion engine power constraint is as follows:
PGEmin≤PGE≤PGEmax
in the formula, PGEFor power of gas internal combustion engines, PGEminIs the minimum limit of the power of the gas internal combustion engine, PGEmaxIs the maximum limit of the power of the gas internal combustion engine, and
Figure BDA0003119236720000171
(2) air conditioning power constraints are as follows:
Qacmin≤Qac≤Qacmax
in the formula, QacFor air conditioning system power, Qacmin、QacmaxThe minimum value and the maximum value of the power of the air conditioning system. (3) The power constraint of the smoke absorption refrigerating unit is as follows:
Qgarmin≤Qgar≤Qgarmax
in the formula, QgarFor power, Q, of a flue gas absorption refrigerating unitgarmin、QgarmaxThe power is the minimum value and the maximum value of the smoke absorption refrigerating unit.
(4) The power constraint of the compression type electric refrigerator is as follows:
Qecmin≤Qec≤Qecmax
in the formula, QecFor power, Q, of compression-type electric refrigeratorsecmax、QecminThe power of the compression type electric refrigerator is up and down limited.
(5) Gas boiler power constraints are as follows:
Hgbmin≤Hgb≤Hgbmax
in the formula, HgbIs the power of the gas boiler, Hgbmax、HgbminThe upper and lower power limits of the gas boiler.
(6) The power constraint of the waste heat boiler is as follows:
0≤Hb,out≤Hb,outmax
in the formula, Hb,outIs the maximum value of the power of the waste heat boiler, Hb,outmaxAnd the maximum value of the power of the waste heat boiler.
Further, the energy storage device capacity constraints include:
(1) battery capacity constraints are as follows:
EBT,min≤EBT≤EBT,max
the above relationship is expressed using state of charge, as follows:
30%≤SOC≤70%
the remaining capacity of the battery must be satisfied at all times:
Figure BDA0003119236720000172
and making an assumption that the initial value and the final value of each period of the storage battery are equal, namely in one period, the storage battery is not discharged or charged, and the following relational expression is satisfied:
Figure BDA0003119236720000181
(2) the heat storage tank capacity constraint is as follows:
Hh,inmin≤Hh,in≤Hh,inmax
Hh,outmin≤Hh,out≤Hh,outmax
in the formula, Hh,inAs input heat value of the heat storage tank, Hh,inmin、Hh,inmaxRespectively represent the upper and lower limits of the input heat value of the heat storage tank, Hh,outIs the output heat value of the heat storage tank, Hh,outmin、Hh,outmaxRespectively represent the upper and lower limits of the output heat value of the heat storage tank.
And an assumption is made that the initial value and the final value of each cycle of the heat storage tank are equal, i.e. in one cycle, the heat storage tank does not exert any force nor store heat.
(3) Ice storage capacity constraints are as follows:
Qic,inmin≤Qic,in≤Qic,inmax
Qic,outmin≤Qic,out≤Qic,outmax
Wic,min≤Wic≤Wic,max
in the formula, Qic,inFor storing iceLoad of cold water consumed, Qic,inmax、Qic,inminUpper and lower cold water load limits, Q, consumed for ice storageic,outCold water load, Q, generated for ice storageic,outmax、Qic,outminThe upper and lower limits of the cold water load generated by ice cold accumulation.
And an assumption is made that the initial value and the final value of each period of ice cold accumulation are equal, namely in one period, the ice cold accumulation has no output and no cold accumulation.
(4) Heat pump capacity constraints, as follows:
0≤Hhp≤Hhpmax
0≤Qhp≤Qhpmax
in the formula, HhpHeat power generated for heat pumps, HhpmaxMaximum value of thermal power, Q, generated by the heat pumphpSpace cooling load, Q, generated for heat pumpshpmaxThe maximum space cooling load generated by the heat pump.
(5) Chiller capacity constraints are as follows:
0≤Hc≤Hcmax
0≤Qc≤Qcmax
in the formula, HcThermal power, H, generated for water chilling unitscmaxMaximum value of thermal power, Q, generated by water chilling unitcCold power, Q, produced for water chiller unitscmaxThe maximum value of the cold power generated by the water chilling unit.
Constraints on the operating mechanism of the multi-energy system include: an electrical power balance constraint, a thermal power balance constraint, a cold power balance constraint, and a power network flow constraint.
Specifically, the electric power balance constraint satisfies the following relation:
Figure BDA0003119236720000191
in the formula (I), the compound is shown in the specification,
Pbuyin order to purchase the electric quantity,
Psellin order to sell the electricity quantity,
PPV,kis the power generation amount of the kth photovoltaic cell, wherein NPVRepresents the total number of photovoltaic cells,
PGE,iis the power generation amount of the ith gas internal combustion engine, wherein NGEIndicates the total number of gas-fired internal combustion engines,
PBT,outis the output power of the storage battery,
PBT,inis the input power of a battery, PacFor the power consumption of the air conditioning system at the time t,
Pecfor the power consumption of the compression-type electric refrigerator at the t-th time,
Php_his the power consumption of the heat pump for heating at time t,
Php_cis the power consumption of the heat pump for cooling at time t,
Ploadthe electric load of the system at the t moment;
specifically, the thermal power balance constraint comprises a space thermal load thermal power balance constraint and a hot water load thermal power balance constraint, wherein the space thermal load thermal power balance constraint satisfies the following relational expression:
Figure BDA0003119236720000192
the hot water load thermal power balance constraint satisfies the following relational expression:
Figure BDA0003119236720000193
in the formula (I), the compound is shown in the specification,
Hgb,j1for the space heat load generated by the jth gas boiler,
Hgb,j2for the hot water load produced by the jth gas boiler,
Hh,out1is the space heat load generated by the heat storage tank,
Hh,out2hot water produced for heat storage tankThe load is applied to the workpiece to be processed,
Hhp1for the space heat load generated by the heat pump,
Hhp2for the hot water load generated by the heat pump,
Hb,out1for the space heat load generated by the waste heat boiler,
Hb,out2is the hot water load generated by the waste heat boiler,
Hload1for the space heat load of the multi-energy system at the t moment,
Hload2for the hot water load of the multi-energy system at time t,
Hh,inis the load of hot water consumed by the heat storage tank,
Ngbthe total number of the gas-fired boilers;
the heat energy consumed by the steam absorption refrigerator is smaller than the heat energy of the smoke generated by the internal combustion engine, the heat energy of the smoke generated by the internal combustion engine and the cylinder sleeve water is larger than the heat energy consumed by the steam absorption refrigerator, the smoke absorption refrigerator, the waste heat boiler and the water chilling unit, and the heat energy is specifically shown as the following two formulas:
Hgas≥Hgar_i
Hgas+Hwater≥Hgar_i+Hsar_i+Hb,in+Hc
in the formula (I), the compound is shown in the specification,
Hgasthe available calorific value of the flue gas discharged by the internal combustion engine at the t moment,
Hwateris the available calorific value of jacket cooling water of the internal combustion engine at the time t,
Hsar_ithe heat consumed by the vapor absorption chiller unit at time t,
Hgar_ithe heat consumed by the smoke absorption refrigerating unit at the t moment,
Hb,inthe heat consumed by the waste heat boiler at the time tt is shown.
Specifically, the cold power balance constraint includes a space cold load cold power balance constraint and a cold water load cold power balance constraint, wherein the space cold load cold power balance constraint satisfies the following relational expression:
Qac+Qsar1+Qgar1+Qec1+Qc1+Qhp+Qic,out1=Qload1
the cold water load cold power balance constraint satisfies the following relational expression:
Qsar2+Qgar2+Qec2+Qic,out2+Qc2=Qic,in+Qload2
in the formula (I), the compound is shown in the specification,
Qacfor the space cooling load generated by the air conditioning system,
Qsar1for space cooling loads generated by vapor absorption chiller units,
Qsar2for the cold water load generated by a vapor absorption chiller,
Qgar1is the space cold load generated by the smoke absorption refrigerating unit,
Qgar2is the cold water load generated by the smoke absorption refrigerating unit,
Qec1the space cold load generated by the compression type electric refrigerator,
Qec2is the cold water load produced by the compression type electric refrigerator,
Qc1is the space cold load generated by the water chilling unit,
Qc2is the cold water load generated by the water chilling unit,
Qhpis the space cold load generated by the heat pump,
Qic,out1the space cold load generated for ice cold accumulation,
Qic,out2the cold water load generated for ice storage,
Qic,inthe cold water load consumed for ice storage,
Qload1for the space cooling load of the multi-energy system at the time t,
Qload2the cold water load of the multi-energy system at the t moment;
specifically, the power network power flow constraints comprise non-out-of-limit constraints on the interactive power flow of the user and the distribution network and upper and lower limit constraints on the node voltage of the whole network, wherein the non-out-of-limit constraints satisfy the following relational expression:
Figure BDA0003119236720000211
Figure BDA0003119236720000212
in the formula (I), the compound is shown in the specification,
Figure BDA0003119236720000213
drawing the line tidal current active power from the network node l to the node n at the t-th moment,
Figure BDA0003119236720000214
the minimum limit value of the line tide active power between the nodes l and n of the hook network,
Figure BDA0003119236720000215
for the maximum limit of the line load active power of the network node l to the node n,
Figure BDA0003119236720000221
for line-flow reactive power of network nodes l to n at time t,
Figure BDA0003119236720000222
for the line-flow reactive power minimum limit from network node l to node n,
Figure BDA0003119236720000223
for the line tide reactive power from a network node l to a node nThe value of the large limit value is,
l,n∈nnode,nnodethe total number of all nodes in the network;
the upper and lower limit constraints satisfy the following relational expression:
Figure BDA0003119236720000224
in the formula (I), the compound is shown in the specification,
Figure BDA0003119236720000225
is the voltage at the node l and,
Figure BDA0003119236720000226
the upper and lower limits of the voltage at node l.
And 3, solving the energy utilization optimization model of the multi-energy system by taking the energy input data and the operation constraint of various energy sources as input data based on the particle swarm optimization to obtain output optimization data of each energy subsystem in the multi-energy system.
In particular, the amount of the solvent to be used,
as shown in fig. 4, step 3 includes:
step 3.1, taking energy input data and operation constraints of various energy sources as input data, and taking hourly output data of each energy subsystem as an optimization variable;
step 3.2, defining each optimized variable as a random particle swarm, and initializing;
3.3, taking the random particle swarm as input data of iterative computation, taking a target function as a fitness function, and iteratively updating an individual extreme value and a global extreme value of the random particle swarm according to preset particle swarm parameters to determine the position and the speed of the particles; obtaining a particle swarm optimization solution after multiple iterations; the preset particle swarm parameters comprise the scale of the particle swarm and the iteration times;
step 3.4, verifying whether the comprehensive utilization efficiency of the cogeneration unit reaches the lowest or not under the obtained particle swarm optimization solution: if the comprehensive utilization efficiency of the cogeneration unit reaches the minimum standard, the obtained particle swarm optimization solution is output optimization data of each energy subsystem in the multi-energy system every hour; and if the comprehensive utilization efficiency of the cogeneration unit does not reach the minimum standard, returning to the step 3.1, and adding a minimum heat-power ratio constraint condition to adjust the heat-power ratio of the cogeneration unit.
Further, in step 3.4, the comprehensive utilization efficiency of the cogeneration unit satisfies the following relational expression:
Figure BDA0003119236720000227
in the formula, QeFor electric energy output by cogeneration units, QhHeat energy, Q, output by cogeneration unitsinThe sum of the heat values of diesel oil and natural gas is consumed for the cogeneration unit.
Further, in step 3.4, the adjustment of the thermoelectric ratio is realized by the constraint condition of the minimum thermoelectric ratio expressed by the following relation:
R′CHPmin=RCHPmin+ΔRCHP
in the formula (II), R'CHPminFor adjusted minimum thermoelectric ratio, RCHPminFor minimum thermoelectric ratio before adjustment, Δ RCHPIs the minimum thermoelectric ratio constraint value.
Example 1.
In the preferred embodiment, taking an off-grid operation park as an example, the energy utilization optimization strategy of the multi-source system facing the platform area side is subjected to simulation verification.
The combined heat and power generation device arranged in the garden contains a afterburning waste heat boiler and has the characteristic of adjustable heat-power ratio. The influence of seasonal factors on the energy demand of the comprehensive energy system is analyzed by introducing a CHP containing an after-burning boiler and adjusting the heat-electricity ratio of CHPs in different seasons, the influence of cold load on the operation of the system is separately considered, and in addition, the optimization strategy also considers the optimization degree of real-time electricity price on the operation cost of the system.
In a park, when the rated power of a cogeneration unit is 3MW, the rated electric efficiency is 33%, and the heat efficiency is limited by the afterburning amount. The rated power of the diesel engine set is 10MW, and the efficiency is 40%. In order to fully utilize the capacity of the gas boiler and the ice storage unit, a 5MW gas boiler and two sets of ice storage air conditioning systems with rated power of 5MW and 4MW are configured as cold and hot energy storage equipment. In addition, the invention assumes that the capacity of the heat storage water tank configured in the park can completely collect the waste heat, considering that the installation cost of the heat storage water tank is low and the waste heat recovery is beneficial to improving the comprehensive energy efficiency. However, due to the high cost of the batteries and the limited number of the charging cars, the maximum storage capacity of the V2G power station is 900kWh, and 500kWh power is required to be supplied every day for charging cars in the park. The cold, hot and electric loads of the park throughout the year are shown in fig. 5 to 8.
Considering that the price of diesel oil is less influenced by load and is relatively low in price. Thus, the price of diesel in the preferred embodiment is only considered to be seasonally affected, and not to be time-varying. In four seasons of spring, summer, autumn and winter, the price of the diesel oil is 1.034, 0.983, 1.025 and 1.051$/L respectively. In addition, with small-scale distributed cogeneration and increase of installed capacity of a gas boiler, the consumption of natural gas is increased day by day, and finally, the phenomenon of insufficient supply of natural gas at the peak time of gas utilization occurs. For this reason, the preferred embodiment simulates time-of-use electricity prices, and optimizes the system energy cost by adopting the time-of-use natural gas prices according to different cold, heat and electricity demands every day. The natural gas prices in different seasons are shown in fig. 9.
(1) CHP thermoelectric ratio optimization
As can be seen from FIG. 10, the thermoelectric ratio R is variedCHPIs increased, eta0Initially, there is an increasing trend, with the energy conversion efficiency reaching a maximum of 78.2% when the thermoelectric ratio is 337%. After that, as the thermoelectric ratio continues to increase, the energy conversion efficiency gradually decreases. This is because the exhaust gas temperature of the boiler is higher and higher with the increase of the post combustion amount, so that the effect of the post combustion on improving the boiler efficiency is smaller and smaller, and further, the energy conversion efficiency of the CHP is gradually reduced.
Fig. 10 also shows the four season CHP heat-power ratio optimization results. In fig. 10, points a, b, c, and d represent the optimal thermoelectric ratios in spring, summer, fall, and winter, respectively, the four season CHP optimal thermal power ratios are 267%, 251%, 276%, and 304%, respectively, and the optimal fuel costs are 11275.9 $ 12847.6 $ 12226.9 $ 13701.2, respectively. Thus, if seasonal considerations are taken into account, the average daily energy cost for a year is $ 12512.9. In addition, the M point represents the annual energy cost optimum thermoelectric ratio and its corresponding overall CHP efficiency, regardless of seasonal factors. Studies have shown that when the thermoelectric ratio is 271%, the average daily energy cost is $ 13159.3.
From the above data it can be seen that:
first, the optimal heat-to-power ratio of the CHP can differ by up to 20% depending on the season, which indicates that seasonal factors have a greater impact on the form of energy demand;
secondly, by optimizing the heat-electricity ratio of the CHP by seasons, the daily average energy cost of the system can be remarkably reduced by 5.2% compared with that before optimization.
In conclusion, the seasonal factors have a large influence on the thermoelectric ratio of the CHP, and the thermoelectric ratio of the CHP is optimized by seasons, so that the thermoelectric load ratio of the CHP can be matched with that of the garden more, and the energy cost of the system is effectively reduced.
(2) Independent consideration of the influence of cold load on system operation
The optimization results in fig. 11 to 14 show that the individual analysis of the system cooling load demand helps to realize peak clipping and valley filling for the electrical load, and improve the system operation stability, and fig. 11, 12, 13 and 14 represent spring, summer, autumn and winter, respectively. P in FIGS. 11 to 14C1,PC2The running states of the two sets of ice storage air conditioning equipment are shown, and when the cold load is greater than zero, the ice storage air conditioning system is in a refrigeration working state; on the contrary, the ice storage air conditioning system is in an ice storage state. In addition, PT1The average daily equivalent electrical load curve (the cold load is converted into the electrical load) before the cold load is considered independently, PT2The electric load and other electric loads required by the ice storage air conditioning system after the cold load is considered independentlyAnd (4) summing.
The data shows that by analyzing the system cold load alone, the equivalent maximum electrical load in summer is reduced from 4.23MW to 3.25MW, which is reduced by 23.2%; the lowest electric load of the system is increased from 0.89MW to 1.27MW, which is increased by 29.9%; in addition, the peak-to-valley difference of the power system of the park is reduced from 3.34MW to 1.98MW, which is reduced by 40.7%. Although the effect of the independent consideration of the cold load on the peak-valley difference of the power-saving system is more obvious in summer, the independent consideration of the cold load of the system can also reduce the peak-valley difference of the power system in other seasons, so that the stability of the system is improved.
By introducing the ice storage device and separately considering the system cooling load, the installed capacity of the diesel generator can be reduced in addition to the stabilization of the electric load fluctuation. In the summer with the maximum annual electric load, P in FIG. 15DE1And PDE2The average daily curves of the output force of the diesel generator before and after the system cold load is considered independently.
As can be seen from fig. 15, the maximum output power of the diesel generator before and after considering the cooling load alone is 3.14MW and 1.98MW, respectively, which indicates that the capacity of the diesel generator can be theoretically reduced by 1.16MW by considering the cooling load alone. Therefore, by considering the system cooling load alone, the configuration capacity of the diesel generator can be greatly reduced. Meanwhile, the data also show that the ratio of the average output to the maximum output of the diesel generator in summer after optimization reaches 78%, which is far higher than 49% before optimization, so that the power generation capacity of the diesel generator is more fully utilized.
(3) Energy cost optimization taking into account real-time gas prices
In FIGS. 16 to 19, PWTWhen the temperature is more than zero, the heat storage water tank is in a heat release state; pWTWhen the temperature is less than zero, the heat storage water tank is in a heat storage state. In one day, the heat energy stored in the heat storage water tank is larger than the heat energy emitted by the heat storage water tank, because the heat loss of the heat storage water tank is 8-11% (different heat losses in different seasons). Fig. 16, 17, 18 and 19 represent spring, summer, fall and winter, respectively.
The optimization results in fig. 16 to 19 show that the energy cost of the system can be reduced by adjusting the heat charging and discharging strategy of the hot water storage tank based on the real-time gas price. Through the reasonable strategy of adjusting CHP, heat production of a gas boiler and heat energy charging and discharging of a hot water storage tank, the heat energy of 7.4MWh of natural gas in a high-price period and a low-price period can be transferred to the low-price period every year on average, so that the annual daily average energy cost is reduced by 7%.
Besides the hot water storage tank, the V2G also plays a role in responding to real-time gas price by transferring electric load, and energy cost is reduced. In FIGS. 20 to 23, when P isV2GWhen the voltage is more than zero, the V2G power changing station supplies electric load to the park; when P is presentV2GWhen the voltage is less than zero, the park power supply system charges the V2G power conversion station. Fig. 20, 21, 22 and 23 represent spring, summer, fall and winter, respectively.
As can be seen from fig. 20 to 23, during the period of low gas price, the park power supply system charges the V2G charging station; during peak gas prices, the V2G conversion station discharges the stored electric energy, and can transfer 1.3MW of electric load on average each day. In addition, the electric quantity discharged for the second time of the V2G power conversion station is smaller than that discharged for the first time in one day, because more electric vehicles need to be charged in a 17: 00-19: 00 time period in consideration of the running time of the electric vehicles. The capacity of the conventional V2G battery replacement station is small, the response real-time gas price of the battery replacement station is influenced, the energy cost of a system is reduced, and the V2G battery replacement station can play a greater role along with the rapid popularization of electric automobiles.
The research shows that: by optimizing the CHP thermoelectric ratio by seasons, the CHP thermoelectric ratio can be better matched with the thermoelectric ratio of the comprehensive energy system, so that the energy cost of the system is reduced by 5.2%. Meanwhile, for the comprehensive energy system running off the grid, the stability of the system running can be seriously influenced by a large peak-to-valley difference of the power load. And by considering the system cold load independently and using the ice storage system to participate in the power load demand response, the peak clipping and valley filling of the power system are indirectly realized. On the premise of meeting the system cold load, the peak-to-valley difference of the summer power load is reduced by 40.7%. In addition, the optimized operation strategy provided by the invention also fully utilizes the heat load transfer capacity of the heat storage water tank, improves the system stability, effectively responds to the real-time gas price and reduces the energy cost. Compared with the prior art, the method has the advantages that the method realizes multi-energy combined supply and coordinated interconnection optimization technology by constructing a complex energy system theory and an optimization method system facing regional energy planning, plays a positive role in innovating an energy service mode, guaranteeing the balance of supply and demand of a power grid, improving the economic operation level of the power grid, improving the overall efficiency of a regional energy system, reducing environmental pollution, developing an intelligent power grid, reducing the total investment and operation management cost of a social energy system and optimizing an energy structure, and provides an effective technical means for greatly reducing the energy consumption in a short period.
The present applicant has described and illustrated embodiments of the present invention in detail with reference to the accompanying drawings, but it should be understood by those skilled in the art that the above embodiments are merely preferred embodiments of the present invention, and the detailed description is only for the purpose of helping the reader to better understand the spirit of the present invention, and not for limiting the scope of the present invention, and on the contrary, any improvement or modification made based on the spirit of the present invention should fall within the scope of the present invention.

Claims (8)

1. An energy consumption optimization control method facing to a platform area side is characterized in that,
the method comprises the following steps:
step 1, based on an energy hub, establishing a multi-energy system model by taking energy input data of each energy subsystem as input data and output data of each energy subsystem as output data;
step 2, based on a multi-energy system model, constructing an energy utilization optimization model of the multi-energy system by taking the optimal 24-hour running cost of the multi-energy system as a target function and combining the constraint of the running limit of the multi-energy system and the constraint of a running mechanism;
wherein the objective function satisfies the following relation:
Figure FDA0003119236710000011
wherein F is the 24-hour running cost of the multi-energy system, QG(t) the natural gas quantity, Q, purchased by the multi-energy system in the tthD(t) the amount of diesel oil purchased by the multi-energy system in the tth hour, CG(t) the purchase price of natural gas at the t hour, CD(t) the purchase price of the diesel oil in the tth hour;
the constraints of the operation limit of the multi-energy system comprise energy supply equipment power constraints and energy storage equipment capacity constraints; the constraints of the operation mechanism of the multi-energy system comprise: electric power balance constraint, thermal power balance constraint, cold power balance constraint and power network power flow constraint;
and 3, solving the energy utilization optimization model of the multi-energy system by taking the energy input data and the operation constraint of various energy sources as input data based on the particle swarm optimization to obtain output optimization data of each energy subsystem in the multi-energy system.
2. The optimizing control method for energy consumption of platform side according to claim 1,
in the step 1, each energy subsystem comprises a cogeneration unit, a diesel generator, a gas boiler unit, an ice storage system, a hot water storage tank system and a V2G battery replacement station.
3. The optimizing control method for energy consumption of platform side according to claim 1,
in step 1, the establishing of the multi-energy system model includes:
step 1.1, building an energy network, wherein the energy network comprises a power network for transmitting electric energy and a fluid network for transmitting heat energy; the power network and the fluid network are coupled by virtue of the pump and the total heat load balance, wherein the total heat load balance meets the following relational expression:
PT=PTh+PTe
in the formula, PThHeat, P, supplied to the user for the fluid networkTeHeat provided to the user for the power network;
step 1.2, based on a network theory, selecting a branch where a pump is located as a tree branch and the rest branches as connecting branches in a power network, and selecting a branch where the pump is located as a tree branch and the rest branches as connecting branches in a fluid network; numbering each branch, each node and each basic loop, and selecting the direction of each branch and each basic loop; each branch is a section of line without branches in the network, each node is an end point of each branch, and each basic loop is a closed loop which comprises and only comprises one continuous branch; the direction of each branch is selected to be the same as the current, and the direction of each basic loop is the clockwise direction;
step 1.3, selecting common nodes of all loops as reference nodes, and establishing a multi-energy system model according to the generalized kirchhoff law by using the following relational expression:
Figure FDA0003119236710000021
in the formula (I), the compound is shown in the specification,
a is an energy network incidence matrix,
h is an extensive flow matrix, the extensive flow matrix is a charge flow matrix in the process of transmitting electric energy by the power network, the extensive flow matrix is an entropy flow matrix in the process of transmitting heat energy by the fluid network,
Bfin the form of a matrix of elementary loops,
and delta X is a matrix of the difference value of the intensities of the two ends of the line, the matrix of the difference value of the intensities of the two ends of the line is a matrix of the potential difference value in the process of transferring electric energy by the power network, and the matrix of the difference value of the intensities of the two ends of the line is a matrix of the temperature difference value of the two ends of the line in the process of transferring heat energy by the fluid network.
4. The optimizing control method for energy consumption of platform side according to claim 1,
in step 2, the power constraint of the energy supply equipment comprises: the method comprises the following steps of (1) gas internal combustion engine power constraint, air conditioner power constraint, flue gas absorption type refrigerating unit power constraint, compression type electric refrigerating unit power constraint, gas boiler power constraint and waste heat boiler power constraint;
the energy storage device capacity constraints include: the capacity of the storage battery is restricted, the capacity of the heat storage tank is restricted, the capacity of the ice storage is restricted, the capacity of the heat pump is restricted, and the capacity of the water chilling unit is restricted.
5. The optimizing control method for energy consumption of platform side according to claim 1,
in step 2, the electric power balance constraint satisfies the following relational expression:
Figure FDA0003119236710000022
in the formula (I), the compound is shown in the specification,
Pbuyin order to purchase the electric quantity,
Psellin order to sell the electricity quantity,
PPV,kis the power generation amount of the kth photovoltaic cell, wherein NPVRepresents the total number of photovoltaic cells,
PGE,iis the power generation amount of the ith gas internal combustion engine, wherein NGEIndicates the total number of gas-fired internal combustion engines,
PBT,outis the output power of the storage battery,
PBT,inis the input power of the storage battery,
Pacfor the power consumption of the air-conditioning system at the time t, PecFor the power consumption of the compression-type electric refrigerator at the time t, Php_hIs the power consumption, P, of the heat pump for heating at time thp_cIs the power consumption, P, of the heat pump for cooling at the time tloadThe electric load of the system at the t moment;
the thermal power balance constraint comprises a space thermal load thermal power balance constraint and a hot water load thermal power balance constraint, wherein the space thermal load thermal power balance constraint satisfies the following relational expression:
Figure FDA0003119236710000031
the hot water load thermal power balance constraint satisfies the following relational expression:
Figure FDA0003119236710000032
in the formula (I), the compound is shown in the specification,
Hgb,j1for the space heat load generated by the jth gas boiler,
Hgb,j2for the hot water load produced by the jth gas boiler,
Hh,out1is the space heat load generated by the heat storage tank,
Hh,out2is the load of hot water generated by the heat storage tank,
Hhp1for the space heat load generated by the heat pump,
Hhp2for the hot water load generated by the heat pump,
Hb,out1for the space heat load generated by the waste heat boiler,
Hb,out2is the hot water load generated by the waste heat boiler,
Hload1for the space heat load of the multi-energy system at the t moment,
Hload2for the hot water load of the multi-energy system at time t,
Hh,inis the load of hot water consumed by the heat storage tank,
Ngbthe total number of the gas-fired boilers;
the cold power balance constraint comprises a space cold load cold power balance constraint and a cold water load cold power balance constraint, wherein the space cold load cold power balance constraint satisfies the following relational expression:
Qac+Qsar1+Qgar1+Qec1+Qc1+Qhp+Qic,out1=Qload1
the cold water load cold power balance constraint satisfies the following relational expression:
Qsar2+Qgar2+Qec2+Qic,out2+Qc2=Qic,in+Qload2
in the formula (I), the compound is shown in the specification,
Qacfor the space cooling load generated by the air conditioning system,
Qsar1for space cooling loads generated by vapor absorption chiller units,
Qsar2for the cold water load generated by a vapor absorption chiller,
Qgar1is the space cold load generated by the smoke absorption refrigerating unit,
Qgar2is the cold water load generated by the smoke absorption refrigerating unit,
Qec1the space cold load generated by the compression type electric refrigerator,
Qec2is the cold water load produced by the compression type electric refrigerator,
Qc1is the space cold load generated by the water chilling unit,
Qc2is the cold water load generated by the water chilling unit,
Qhpis the space cold load generated by the heat pump,
Qic,out1the space cold load generated for ice cold accumulation,
Qic,out2the cold water load generated for ice storage,
Qic,inthe cold water load consumed for ice storage,
Qload1for the space cooling load of the multi-energy system at the time t,
Qload2the cold water load of the multi-energy system at the t moment;
the power network power flow constraint comprises non-out-of-limit constraint on user and distribution network interactive power flow and upper and lower limit constraint on the whole network node voltage, wherein the non-out-of-limit constraint satisfies the following relational expression:
Figure FDA0003119236710000051
Figure FDA0003119236710000052
in the formula (I), the compound is shown in the specification,
Figure FDA0003119236710000053
for the line load active power of network nodes l to n at time t,
Figure FDA0003119236710000054
for the line-load active power minimum limit from network node l to node n,
Figure FDA0003119236710000055
for the maximum limit of the line load active power of the network node l to the node n,
Figure FDA0003119236710000056
for line-flow reactive power of network nodes l to n at time t,
Figure FDA0003119236710000057
for the line-flow reactive power minimum limit from network node l to node n,
Figure FDA0003119236710000058
for the line flow reactive power maximum limit from network node l to node n,
l,n∈nnode,nnodethe total number of all nodes in the network;
the upper and lower limit constraints satisfy the following relational expression:
Figure FDA0003119236710000059
in the formula (I), the compound is shown in the specification,
Figure FDA00031192367100000510
is the voltage at the node l and,
Figure FDA00031192367100000511
the upper and lower limits of the voltage at the hook node l.
6. The optimizing control method for energy consumption of platform side according to claim 1,
the step 3 comprises the following steps:
step 3.1, taking energy input data and operation constraints of various energy sources as input data, and taking hourly output data of each energy subsystem as an optimization variable;
step 3.2, defining each optimized variable as a random particle swarm, and initializing;
3.3, taking the random particle swarm as input data of iterative computation, taking a target function as a fitness function, and iteratively updating an individual extreme value and a global extreme value of the random particle swarm according to preset particle swarm parameters to determine the position and the speed of the particles; obtaining a particle swarm optimization solution after multiple iterations; the preset particle swarm parameters comprise the size of the particle swarm and the iteration times;
step 3.4, verifying whether the comprehensive utilization efficiency of the cogeneration unit reaches the lowest or not under the obtained particle swarm optimization solution: if the comprehensive utilization efficiency of the cogeneration unit reaches the minimum standard, the obtained particle swarm optimization solution is output optimization data of each energy subsystem in the multi-energy system every hour; and if the comprehensive utilization efficiency of the cogeneration unit does not reach the minimum standard, returning to the step 3.1, and adding a minimum heat-power ratio constraint condition to adjust the heat-power ratio of the cogeneration unit.
7. The optimizing control method for energy consumption of platform side according to claim 6,
in step 3.4, the comprehensive utilization efficiency of the cogeneration unit meets the following relational expression:
Figure FDA0003119236710000061
in the formula, QeFor electric energy output by cogeneration units, QhHeat energy, Q, output by cogeneration unitsinThe sum of the heat values of diesel oil and natural gas is consumed for the cogeneration unit.
8. The optimizing control method for energy consumption of platform side according to claim 6,
in step 3.4, the thermoelectric ratio is adjusted by the constraint condition of the minimum thermoelectric ratio represented by the following relational expression:
R′CHPmin=RCHPmin+ΔRCHP
in the formula (II), R'CHPminFor adjusted minimum thermoelectric ratio, RCHPminFor minimum thermoelectric ratio before adjustment, Δ RCHPIs the minimum thermoelectric ratio constraint value.
CN202110672510.3A 2021-06-17 2021-06-17 Energy consumption optimization control method for side of transformer area Pending CN113328443A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110672510.3A CN113328443A (en) 2021-06-17 2021-06-17 Energy consumption optimization control method for side of transformer area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110672510.3A CN113328443A (en) 2021-06-17 2021-06-17 Energy consumption optimization control method for side of transformer area

Publications (1)

Publication Number Publication Date
CN113328443A true CN113328443A (en) 2021-08-31

Family

ID=77423771

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110672510.3A Pending CN113328443A (en) 2021-06-17 2021-06-17 Energy consumption optimization control method for side of transformer area

Country Status (1)

Country Link
CN (1) CN113328443A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115579923A (en) * 2022-10-18 2023-01-06 国网上海能源互联网研究院有限公司 Energy storage optimization operation method and device for multiple scenes and multiple targets

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105676646A (en) * 2016-03-11 2016-06-15 国网天津市电力公司 Linearization method for optimized operation of combined cooling heating and power supply system
CN105869075A (en) * 2016-04-19 2016-08-17 东南大学 Economic optimization scheduling method for cold, heat and electricity combined supply type miniature energy grid
CN110363363A (en) * 2019-07-31 2019-10-22 广东电网有限责任公司 A kind of dual-stage scheduling model optimization method for integrated energy system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105676646A (en) * 2016-03-11 2016-06-15 国网天津市电力公司 Linearization method for optimized operation of combined cooling heating and power supply system
CN105869075A (en) * 2016-04-19 2016-08-17 东南大学 Economic optimization scheduling method for cold, heat and electricity combined supply type miniature energy grid
CN110363363A (en) * 2019-07-31 2019-10-22 广东电网有限责任公司 A kind of dual-stage scheduling model optimization method for integrated energy system

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
崔秀清: "热电联产在混合储能系统中的优化配置", 《科技视野》 *
文俊中: "基于能量网络理论的分布式供能系统分析", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *
施泉生: "风电与电动汽车协同并网调度环境及经济模型", 《上海电力学院学报》 *
施锦月: "基于热电比可调模式的区域综合能源系统双层优化运行", 《电网技术》 *
罗蒙: "《绿色医院节能实用手册》", 31 July 2018, 上海交通大学出版社 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115579923A (en) * 2022-10-18 2023-01-06 国网上海能源互联网研究院有限公司 Energy storage optimization operation method and device for multiple scenes and multiple targets
CN115579923B (en) * 2022-10-18 2024-02-23 国网上海能源互联网研究院有限公司 Energy storage optimization operation method and device for multiple scenes and multiple targets

Similar Documents

Publication Publication Date Title
Yan et al. An integrated design for hybrid combined cooling, heating and power system with compressed air energy storage
Liu et al. Two-phase collaborative optimization and operation strategy for a new distributed energy system that combines multi-energy storage for a nearly zero energy community
CN104716644B (en) Renewable energy source cooling, heating and power microgrid system and control method
CN108631343A (en) One kind is provided multiple forms of energy to complement each other energy internet Optimization Scheduling
CN111860937B (en) Combined cooling heating and power type microgrid optimization method based on improved multi-target Husky algorithm
CN105869075A (en) Economic optimization scheduling method for cold, heat and electricity combined supply type miniature energy grid
Liu et al. Co-optimization of a novel distributed energy system integrated with hybrid energy storage in different nearly zero energy community scenarios
WO2019205561A1 (en) Cchp micro-grid structure including compressed air energy storage and operation method therefor
CN111737884B (en) Multi-target random planning method for micro-energy network containing multiple clean energy sources
CN109696891A (en) Micro- energy net system and its progress control method comprising air source heat pump and energy storage
CN111144707A (en) Multi-energy system collaborative planning modeling method based on energy hub
CN113762708A (en) Park level comprehensive energy system planning method considering multi-target cooperation
CN109184916A (en) The method of comprehensive energy router device and energy conversion
CN106779471A (en) A kind of multipotency interconnects alternating current-direct current mixing micro-capacitance sensor system and Optimal Configuration Method
CN112600253B (en) Park comprehensive energy collaborative optimization method and equipment based on optimal energy utilization efficiency
CN113033900A (en) Park level comprehensive energy system capacity optimal configuration method and system
Guo et al. Thermodynamic performance analyses and optimization design method of a novel distributed energy system coupled with hybrid-energy storage
CN115173470A (en) Comprehensive energy system scheduling method and system based on power grid peak shaving
CN204407890U (en) A kind of regenerative resource cool and thermal power micro-grid system
CN108197412B (en) Multi-energy coupling energy management system and optimization method
CN113328443A (en) Energy consumption optimization control method for side of transformer area
CN213783243U (en) Comprehensive energy system operation optimizing device for industrial park
Yang et al. Economic dispatching of CCHP microgrid considering grid-connected revenue
Fan et al. Research on collaborative optimization model of park-level integrated energy system participating in power peak shaving
Hui Urban multi-regional energy system planning considering integral energy operating optimization

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210831

RJ01 Rejection of invention patent application after publication