CN111768036B - Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid - Google Patents
Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid Download PDFInfo
- Publication number
- CN111768036B CN111768036B CN202010608816.8A CN202010608816A CN111768036B CN 111768036 B CN111768036 B CN 111768036B CN 202010608816 A CN202010608816 A CN 202010608816A CN 111768036 B CN111768036 B CN 111768036B
- Authority
- CN
- China
- Prior art keywords
- power
- distribution system
- node
- constraint
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 108
- 230000002452 interceptive effect Effects 0.000 title claims abstract description 48
- 238000005457 optimization Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000003993 interaction Effects 0.000 claims description 15
- 238000010521 absorption reaction Methods 0.000 claims description 10
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000020169 heat generation Effects 0.000 claims description 4
- 238000010248 power generation Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims 1
- 230000033228 biological regulation Effects 0.000 abstract description 8
- 239000007789 gas Substances 0.000 description 37
- 238000010586 diagram Methods 0.000 description 10
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 238000011160 research Methods 0.000 description 5
- 239000003345 natural gas Substances 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000004378 air conditioning Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a power optimization method for interactive operation of a comprehensive energy distribution system and an upper power grid, which comprises the steps of obtaining basic data, establishing a model, determining constraint conditions, determining an adjustment domain, solving the power optimization model for interactive operation of the comprehensive energy distribution system and the upper power grid by utilizing an MATLAB platform and a YALMIP tool box, outputting related results and the like. According to the power optimization method for the interactive operation of the comprehensive energy distribution system and the upper power grid, provided by the invention, on the premise of ensuring each load demand, an energy distribution scheme of an energy station for the interactive operation power optimization with the upper power grid can be found; from the aspect of meeting the regulation and control instructions of the upper power grid, the energy supply of the energy station is established by taking the deviation between the actual power and the expected power which are minimized to interact with the upper power grid as an objective function and taking the power distribution system trend constraint and the like as constraint conditions as a power optimization model for the interactive operation of the comprehensive energy distribution system and the upper power grid, and the energy supply of the energy station is determined by solving the model based on the MATLAB platform and the YALMIP tool box.
Description
Technical Field
The invention belongs to the technical field of comprehensive energy distribution systems, and particularly relates to a power optimization method for interactive operation of a comprehensive energy distribution system and an upper power grid.
Background
The society development has a strong dependence on energy, and exhaustion of fossil energy and pollution to the environment are disadvantageous factors such as the fact that people have to seek clean energy to replace traditional fossil energy. With the continuous maturing of researches on central air conditioning, gas turbines and the like, the coupling degree of each energy system is increased, the traditional energy systems are independently planned and operated, the coupling condition among the energy sources is not fully considered, and the comprehensive energy distribution system (integrated energy distribution system, IEDS) considers the coupling relation among the energy sources on the basis of the traditional energy systems, and relates to links of energy source production, transportation, distribution, conversion and the like, and the systems are uniformly planned and processed, so that the whole system is optimized, and a remarkable effect is produced on improving the utilization rate of the energy sources. The IEDS may also improve the reliability of the system, and when one energy system fails, it may be powered by another energy system, resulting in an overall improvement in reliability of the system. In addition, innovative technologies based on natural gas have also led to the conversion of energy structures supplied by electric energy as a main energy source in the past, and research on IEDS for energy distribution by using electric power and a natural gas network as a main energy source network has been paid great attention.
The energy systems are in different grades and have complex differences in space-time, and the characteristic brings challenges to the analysis and research of the multi-energy system in the aspects of energy conversion, storage, distribution and the like. The energy station is an energy input/output multiport model, can describe energy supply and energy demand with high abstraction, can adjust the energy supply of various devices by adjusting the interrelationship among different energies, further realizes multi-energy complementation, and has various advantages of economy, environmental protection, energy efficiency and the like. The energy station is able to integrate any number of energy carriers, thereby providing a high degree of flexibility in system modeling.
The upper grid supplies power to the IEDS, and in order to ensure safe and stable operation of the upper grid, the upper grid may wish to consume more or less power by the IEDS. The traditional power grid can only achieve the aim by increasing or decreasing the load, but the adjustment of the load involves a user side, and certain difficulty exists in achieving the aim, and the IEDS can achieve the adjustment of the interactive power of the comprehensive energy distribution system and the upper power grid by adjusting the distribution relation of each energy of the energy stations, so that the electric energy consumption requirement of the upper power grid on the IEDS is met. However, the related technology in this aspect has few researches, and the problems of safe operation of a natural gas network and feasible regulation of the interactive power of the comprehensive energy distribution system and the upper power grid are not considered in the related researches, so that the situation that the difference from the expected result is large easily occurs.
Disclosure of Invention
In order to realize the adjustment of the interactive power of the comprehensive energy distribution system and the upper power grid, the invention provides a power optimization method for the interactive operation of the comprehensive energy distribution system and the upper power grid.
The invention provides a power optimization method for the interactive operation of a comprehensive energy distribution system and an upper power grid, which comprises the following steps:
s1, collecting basic data of a comprehensive energy distribution system to be researched;
s2, establishing a power optimization model for the interactive operation of the comprehensive energy distribution system and the upper power grid through basic data, and taking the difference value between the actual value and the expected value of the interactive power of the minimum comprehensive energy distribution system and the upper power grid as an objective function;
s3, determining constraint conditions, wherein the set constraint conditions comprise: power flow constraint of a power distribution system, operation voltage level constraint, branch current limitation, pipeline flow constraint of a gas distribution system, gas pressure level constraint, pipeline flow constraint, energy station energy flow constraint and equipment output constraint;
s4, determining an adjustable range of the interaction power of the comprehensive energy distribution system and the upper power grid, and setting an expected value of the interaction power of the comprehensive energy distribution system and the upper power grid according to the adjustable range;
and S5, solving a power optimization model of the interactive operation of the comprehensive energy distribution system and the upper power grid to obtain the actual value of the interactive power of the comprehensive energy distribution system and the upper power grid and the corresponding value of the energy balance relation of each device of the energy station.
Preferably, the base data comprises: the grid structure, the electric load level and the electric parameters of the comprehensive energy distribution system, the network topology, the air load level, the pipeline parameters of the air distribution system, the multi-energy load of the energy station and the equipment efficiency.
Preferably, the expression of the objective function is:
Obj＝min|P actual -P expect | (1)
wherein P is actual To interact with the upper grid for actual power, P expect Power is desired for interaction with the upper grid.
Preferably, the expression of the power flow constraint of the power distribution system is:
wherein i epsilon u (j) is a set adjacent to node j and node i is a line head node; k epsilon v (j) is a set adjacent to node j and node k is a line end node; p (P) ij Active power flow at the head end of the line; p (P) jk Active power flow for the end of the line; p (P) j Net injection of active power for node j; q (Q) ij Reactive power flow is the head end of the line; q (Q) jk Reactive power flow at the end of the line; q (Q) j Net injection of reactive power for node j; v (V) i 、V j The voltage amplitudes of the head and the tail of the line are respectively; i ij For the magnitude of the current flowing through the line; r is (r) ij And x ij The resistance value and the reactance value of the circuit are respectively;
the expression of the operating voltage level constraint is:
V i，min ≤V i ≤V i，max (3)
wherein V is i For node i voltage, V i,max And V i,min Respectively the upper limit and the lower limit of the voltage amplitude of the node i;
the expression of the branch current limit is:
I ij ≤I ij,max (4)
wherein I is ij I for the amplitude of the current flowing through the line ij,max The maximum allowable value of the current flowing through the line;
the expression of the flow constraint of the gas distribution system is as follows:
wherein D is the diameter of the pipeline; l is the length of the pipeline; s is the gas density; f is the coefficient of friction; q mn Flow rate for the flow through the pipe; p is p m And p n The pressure values of the head and the tail nodes of the pipeline are respectively; matrix A is a branch node association matrix, a ij The value of (2) depends on the relation between the node i and the branch j, when the node i and the branch j are not associated, the value is 0, when the node i and the branch j are associated and the gas flow in the pipeline flows into the node, the value is 1, and when the node i and the gas flow in the pipeline flows out of the node, the value is-1; q (Q) branch A vector consisting of branch flows; q (Q) node Vector composed of the traffic load of each node;
the expression of the air pressure level constraint is as follows:
p m，min ≤p m ≤p m,max (6)
wherein p is m For the node pressure value, p m,max And p m,min Respectively the upper limit and the lower limit of the node pressure value;
the expression of the pipeline flow constraint is as follows:
q mn，min ≤q mn ≤q mn,max (7)
wherein q mn For the flow value of the pipeline, q mn,max And q mn,min The upper limit and the lower limit of the pipeline flow value are respectively;
the expression of the energy station energy flow constraint is as follows:
wherein L is el 、L hl 、L cl Respectively electric, thermal and cold loads; p (P) el Interacting power for the energy station and the power distribution system; η (eta) GT_el 、η GT_hl 、η GBl 、η ECl 、η ACl The power generation efficiency of the gas turbine, the heat generation efficiency of the gas turbine, the efficiency of the gas boiler, the efficiency of the electric refrigerator and the efficiency of the absorption refrigerator are respectively; p (P) GT_inl 、P GT_el 、P GT_hl The power of gas input, electric output and heat output of the gas turbine are respectively; p (P) GB_inl 、P GB_outl The power of gas input and heat output of the gas boiler are respectively; p (P) EC_inl 、P EC_outl The power of the electric input and the cold output of the electric refrigerator are respectively; p (P) AC_inl 、P AC_outl The power of the heat input and the cold output of the absorption refrigerator are respectively;
the expression of the equipment output constraint is as follows:
wherein P is GT_inl,min 、P GT_inl,max The upper and lower limits of the output of the gas turbine are respectively; p (P) GB_inl,min 、P GB_inl,max The upper and lower limits of the output of the gas boiler are respectively; p (P) EC_inl,min 、P EC_inl,max The upper and lower limits of the output of the electric refrigerator are respectively; p (P) AC_inl,min 、P AC_inl,max The upper and lower limits of the output of the absorption refrigerator are respectively set.
Preferably, the adjustable range of the interaction power with the upper power grid is the upper and lower limits of the interaction power with the upper power grid.
Preferably, the power optimization model of the interactive operation of the comprehensive energy distribution system and the upper power grid is solved through the MATLAB platform and the YALMIP toolbox.
The invention provides a power optimization method for interactive operation of a comprehensive energy distribution system and an upper power grid, which is used for modeling the power optimization problem of the interactive operation of the comprehensive energy distribution system and the upper power grid, taking the deviation between the actual power and the expected power of the interaction of the minimum comprehensive energy distribution system and the upper power grid as an objective function and taking the power flow constraint of the distribution system, the operation voltage level constraint, the branch current constraint, the power flow constraint of a distribution system, the air pressure level constraint, the pipeline flow constraint, the energy station energy flow constraint, the equipment output constraint and the like as constraint conditions. The method is characterized in that the method is realized through energy distribution of an energy station, a reference range of a regulation instruction is provided for an upper power grid according to a feasible regulation domain of the interaction power of the comprehensive energy distribution system and the upper power grid, and then the deviation between the actual power and the expected power of the interaction of the comprehensive energy distribution system and the upper power grid is minimized as an objective function, so that the power optimization of the interaction operation of the comprehensive energy distribution system and the upper power grid is realized. In the aspect of model establishment, the method comprehensively considers the models of the power distribution system and the air distribution system, and introduces the models of the energy stations to describe the energy supply and the energy demand.
Drawings
FIG. 1a is a schematic diagram of an IEEE33 node power distribution system in an IEDS embodiment;
FIG. 1b is a schematic diagram of an 11-node gas distribution system in an IEDS embodiment;
FIG. 1c is a schematic diagram of an energy station 1 in an IEDS embodiment;
FIG. 1d is a schematic diagram of an energy station 2 in an IEDS embodiment;
FIG. 1e is a schematic diagram of two energy stations in an IEDS embodiment;
FIG. 2 is a flow chart of a power optimization method for the interactive operation of the integrated energy distribution system and the upper power grid according to the invention;
FIG. 3 is a graph showing upper and lower limits of the integrated energy distribution system and the upper grid interactive power regulation;
FIG. 4 is a graph of optimizing the power balance of the output and the consumed electric power of each device of the obtained energy station 1;
FIG. 5 is a graph of optimizing the thermal power balance of the output and consumption of each device of the resulting energy station 1;
FIG. 6 is a graph of optimizing the resulting power plant 1 output versus the cold power consumed;
FIG. 7 is a graph of optimizing the power balance of the output and consumption of each device of the resulting energy station 2;
FIG. 8 is a graph of optimizing the thermal power balance of the output and consumption of each device of the resulting energy station 2;
FIG. 9 is a graph of optimizing the resulting power plant 2 output versus consumed cold power balance;
FIG. 10 is a comparison of actual values and expected values obtained from power optimization of the integrated energy distribution system and the superior grid interoperation.
Detailed Description
The invention provides a power optimization method for the interactive operation of a comprehensive energy distribution system and an upper power grid, which can find an energy distribution scheme of a power optimized energy station for the interactive operation of the comprehensive energy distribution system and the upper power grid on the premise of ensuring the requirements of various loads. The invention considers the interactive power optimization of the upper power grid and the comprehensive energy distribution system by the upper power grid regulation and control instruction, establishes a power optimization model for the interactive operation of the comprehensive energy distribution system and the upper power grid by taking the deviation between the actual power and the expected power of the minimized interaction of the comprehensive energy distribution system and the upper power grid as an objective function and taking the power flow constraint of the distribution system, the operation voltage level constraint, the branch current constraint, the power flow constraint of the distribution system, the air pressure level constraint, the pipeline flow constraint, the energy station energy flow constraint, the equipment output constraint and the like as constraint conditions, and solves the model based on a MATLAB platform and a YALMIP tool box to determine the energy station energy flow balance relation.
The power optimization method for the interactive operation of the comprehensive energy distribution system and the upper power grid is further described by the detailed description of the preferred embodiment with reference to the accompanying drawings.
In the embodiment of the invention, the IEDS example consists of an IEEE33 node power distribution system, an 11 node gas distribution system and two energy stations, fig. 1a is a schematic diagram of the IEEE33 node power distribution system in the IEDS example, fig. 1b is a schematic diagram of the 11 node gas distribution system in the IEDS example, fig. 1c is a schematic diagram of the energy station 1 in the IEDS example, fig. 1d is a schematic diagram of the energy station 2 in the IEDS example, and fig. 1e is a schematic diagram of the two energy stations in the IEDS example.
First, basic data of the IEDS calculation example shown in fig. 1 is obtained, including basic data of a grid structure, an electric load level, an electric parameter of the integrated energy power distribution system, and network topology, an air load level, pipeline parameters, a multi-energy load of an energy station, equipment efficiency and the like of the air distribution system. The access position of the energy station is set, in this embodiment, 6 nodes and 13 nodes of the power distribution system are used as the electric energy input positions of the energy station, and N4 nodes and N9 nodes of the power distribution system are used as the natural gas energy input positions of the energy station. The reference voltage of the system is set to be 12.66kV, and the air pressure value of the air source point N1 is 75mbar.
Secondly, establishing a power optimization model for the interactive operation of the comprehensive energy distribution system and the upper power grid by using the acquired basic data, setting a difference value between an actual power value and an expected value of the minimum interaction of the comprehensive energy distribution system and the upper power grid as an objective function, wherein the expression of the objective function is as follows:
Obj＝min|P actual -P expect | (1)
wherein P is actual For the actual power interacting with the upper grid, P expect Is the desired power to interact with the upper grid. In the embodiment of the invention, the power interacted between the comprehensive energy power distribution system and the upper power grid is the power value of the power distribution system node 1 connected with the upper power grid.
Secondly, determining constraint conditions, and setting constraint conditions including power distribution system power flow constraint, operation voltage level constraint, branch current constraint, gas distribution system power flow constraint, gas pressure level constraint, pipeline flow constraint, energy station energy flow constraint and equipment output constraint. The power distribution system tide constraint expression is as follows:
wherein i epsilon u (j) is a set adjacent to node j and node i is a line head node; k epsilon v (j) is a set adjacent to node j and node k is a line end node; p (P) ij Active power flow at the head end of the line; p (P) jk Active power flow for the end of the line; p (P) j Net injection of active power for node j; q (Q) ij Reactive power flow is the head end of the line; q (Q) jk Reactive power flow at the end of the line; q (Q) j Net injection of reactive power for node j; v (V) i 、V j The voltage amplitudes of the head and the tail of the line are respectively; i ij For the magnitude of the current flowing through the line; r is (r) ij And x ij The resistance value and the reactance value of the line, respectively.
The expression of the operating voltage level constraint is:
V i，min ≤V i ≤V i，max (3)
wherein V is i For node i voltage, V i,max And V i,min The upper and lower limits of the voltage amplitude of the node i are respectively set.
The expression of the branch current limit is:
I ij ≤I ij,max (4)
wherein I is ij I for the amplitude of the current flowing through the line ij,max Is the maximum allowable value of the current flowing through the line.
The expression of the flow constraint of the gas distribution system is as follows:
wherein D is the diameter of the pipeline; l is the length of the pipeline; s is the gas density; f is the coefficient of friction; qmn is the flow rate flowing through the pipe; p is p m And p n The pressure values of the head and the tail nodes of the pipeline are respectively; matrix A is a branch node association matrix, a ij The value of (1) depends on the relation between the node i and the branch j, when the node i and the branch j are not related, the value is 0, when the node i and the branch j are related and the gas flow in the pipeline flows into the node, the value is 1, and when the node i and the branch j are related and the gas flow in the pipeline flows outWhen the node is connected, the value is-1; q (Q) branch A vector consisting of branch flows; q (Q) node And a vector formed by the traffic loads of all the nodes.
The expression of the air pressure level constraint is as follows:
p m，min ≤p m ≤p m,max (6)
wherein p is m For the node pressure value, p m,max And p m,min Respectively the upper limit and the lower limit of the node pressure value.
The expression of the pipeline flow constraint is as follows:
q mn，min ≤q mn ≤q mn,max (7)
wherein q mn For the flow value of the pipeline, q mn,max And q mn,min The upper and lower limits of the pipeline flow value are respectively set.
The expression of the energy station energy flow constraint is as follows:
wherein L is el 、L hl 、L cl Respectively electric, thermal and cold loads; p (P) el Interacting power for the energy station and the power distribution system; η (eta) GT_el 、η GT_hl 、η GBl 、η ECl 、η ACl The power generation efficiency of the gas turbine, the heat generation efficiency of the gas turbine, the efficiency of the gas boiler, the efficiency of the electric refrigerator and the efficiency of the absorption refrigerator are respectively; p (P) GT_inl 、P GT_el 、P GT_hl The power of gas input, electric output and heat output of the gas turbine are respectively; p (P) GB_inl 、P GB_outl The power of gas input and heat output of the gas boiler are respectively; p (P) EC_inl 、P EC_outl The power of the electric input and the cold output of the electric refrigerator are respectively; p (P) AC_inl 、P AC_outl The power of the heat input and the cold output of the absorption refrigerator are respectively.
The expression of the equipment output constraint is as follows:
wherein P is GT_inl,min 、P GT_inl,max The upper and lower limits of the output of the gas turbine are respectively; p (P) GB_inl,min 、P GB_inl,max The upper and lower limits of the output of the gas boiler are respectively; p (P) EC_inl,min 、P EC_inl,max The upper and lower limits of the output of the electric refrigerator are respectively; p (P) AC_inl,min 、P AC_inl,max The upper and lower limits of the output of the absorption refrigerator are respectively set.
In the embodiment of the invention, the power generation efficiency of the gas turbine in the energy station is set to be 0.3, the heat generation efficiency of the gas turbine is set to be 0.4, the efficiency of the gas boiler is set to be 0.9, the efficiency of the electric refrigerator is set to be 4, and the efficiency of the absorption refrigerator is set to be 1.7.
And secondly, determining an adjustable range of the interactive power of the comprehensive energy distribution system and the upper power grid, namely, the upper limit and the lower limit of the interactive power of the comprehensive energy distribution system and the upper power grid, and providing a reference for setting an expected value of the interactive power of the comprehensive energy distribution system and the upper power grid. In the embodiment of the present invention, the upper and lower limit values of the power distribution system node 1 under the condition that all constraint conditions are satisfied are determined, and the result is shown in fig. 3.
And finally, solving a power optimization model for the interactive operation of the comprehensive energy distribution system and the upper power grid by utilizing the MATLAB platform and the YALMIP tool box, and outputting a related result, wherein the related result comprises an actual value of the interactive power of the comprehensive energy distribution system and the upper power grid, a corresponding value of the energy balance relation of each device of the energy station and the like.
The optimal configuration results of the embodiment of the present invention are shown in fig. 4 to 10, where fig. 4, fig. 5, and fig. 6 are respectively the power balance relationship between the output and the consumed power of each device of the energy station 1, fig. 7, fig. 8, and fig. 9 are respectively the power balance relationship between the output and the consumed power of each device of the energy station 2, and fig. 10 is a comparison between the actual power and the expected power of the interaction between the integrated energy distribution system and the upper grid.
The obtained result can be used for obtaining that the power optimization model for the interactive operation of the comprehensive energy distribution system and the upper power grid provided by the invention has reasonability, so that the characteristic of the flexibility of the energy station is reflected. As can be seen from fig. 10, the power optimization method for the interactive operation of the integrated energy distribution system and the upper power grid can make the actual power interacted with the upper power grid very similar to the expected power, thereby meeting the requirement of the regulation and control instruction of the upper power grid. Figures 4-9 can provide a power output optimization configuration scheme for each device of the energy station. The power optimization method for the interactive operation of the comprehensive energy distribution system and the upper power grid can provide reference for the upper power grid to give out instructions, so that the instructions are in a reasonable regulation and control range, and an effective scheduling scheme is provided for scheduling personnel.
While the present invention has been described in detail through the foregoing description of the preferred embodiment, it should be understood that the foregoing description is not to be considered as limiting the invention. Many modifications and substitutions of the present invention will become apparent to those of ordinary skill in the art upon reading the foregoing. Accordingly, the scope of the invention should be limited only by the attached claims.
Claims (5)
1. The power optimization method for the interactive operation of the comprehensive energy distribution system and the upper power grid is characterized by comprising the following steps of:
s1, collecting basic data of a comprehensive energy distribution system to be researched;
s2, establishing a power optimization model for the interactive operation of the comprehensive energy distribution system and the upper power grid through basic data, and taking the difference value between the actual value and the expected value of the interactive power of the minimum comprehensive energy distribution system and the upper power grid as an objective function;
s3, determining constraint conditions, wherein the set constraint conditions comprise: power flow constraint of a power distribution system, operation voltage level constraint, branch current limitation, pipeline flow constraint of a gas distribution system, gas pressure level constraint, pipeline flow constraint, energy station energy flow constraint and equipment output constraint;
s4, determining an adjustable range of the interaction power of the comprehensive energy distribution system and the upper power grid, and setting an expected value of the interaction power of the comprehensive energy distribution system and the upper power grid according to the adjustable range;
s5, solving a power optimization model of the interactive operation of the comprehensive energy distribution system and the upper power grid to obtain an actual value of the interactive power of the comprehensive energy distribution system and the upper power grid and a corresponding value of the energy balance relation of each device of the energy station;
the power distribution system tide constraint expression is as follows:
wherein i epsilon u (j) is a set adjacent to node j and node i is a line head node; k epsilon v (j) is a set adjacent to node j and node k is a line end node; p (P) ij Active power flow at the head end of the line; p (P) jk Active power flow for the end of the line; p (P) j Net injection of active power for node j; q (Q) ij Reactive power flow is the head end of the line; q (Q) jk Reactive power flow at the end of the line; q (Q) j Net injection of reactive power for node j; v (V) i 、V j The voltage amplitudes of the head and the tail of the line are respectively; i ij For the magnitude of the current flowing through the line; r is (r) ij And x ij The resistance value and the reactance value of the circuit are respectively;
the expression of the operating voltage level constraint is:
V x，min ≤V x ≤V x，max (3)
wherein V is x For node x voltage, V x,max And V x,min Respectively the upper limit and the lower limit of the voltage amplitude of the node x; the expression of the branch current limit is:
I ij ≤I ij,max (4)
wherein I is ij I for the amplitude of the current flowing through the line ij,max The maximum allowable value of the current flowing through the line; the expression of the flow constraint of the gas distribution system is as follows:
wherein D is the diameter of the pipeline; l is the length of the pipeline; s is the gas density; f is the coefficient of friction; q mn Is a pipeline flow value; p is p m And p n The pressure values of the head and the tail nodes of the pipeline are respectively; matrix A is a branch node association matrix, a ij The value of (2) depends on the relation between the node i and the branch j, when the node i and the branch j are not associated, the value is 0, when the node i and the branch j are associated and the gas flow in the pipeline flows into the node, the value is 1, and when the node i and the gas flow in the pipeline flows out of the node, the value is-1; q (Q) branch A vector consisting of branch flows; q (Q) node Vector composed of the traffic load of each node;
the expression of the air pressure level constraint is as follows:
p x，min ≤p x ≤p x,max (6)
wherein p is x For the node pressure value, p x,max And p x,min Respectively the upper limit and the lower limit of the node pressure value;
the expression of the pipeline flow constraint is as follows:
q mn，min ≤q mn ≤q mn,max (7)
wherein q mn For the flow value of the pipeline, q mn,max And q mn,min The upper limit and the lower limit of the pipeline flow value are respectively;
the expression of the energy station energy flow constraint is as follows:
wherein L is el 、L hl 、L cl Respectively electric, thermal and cold loads; p (P) el Interacting power for the energy station and the power distribution system; η (eta) GT_el 、η GT_hl 、η GBl 、η ECl 、η ACl The power generation efficiency of the gas turbine, the heat generation efficiency of the gas turbine, the efficiency of the gas boiler, the efficiency of the electric refrigerator and the efficiency of the absorption refrigerator are respectively; p (P) GT_inl 、P GT_el 、P GT_hl Gas input and electric power transmission of gas turbineOutput power of heat output; p (P) GB_inl 、P GB_outl The power of gas input and heat output of the gas boiler are respectively; p (P) EC_inl 、P EC_outl The power of the electric input and the cold output of the electric refrigerator are respectively; p (P) AC_inl 、P AC_outl The power of the heat input and the cold output of the absorption refrigerator are respectively;
the expression of the equipment output constraint is as follows:
wherein P is GT_inl,min 、P GT_inl,max The upper and lower limits of the output of the gas turbine are respectively; p (P) GB_inl,min 、P GB_inl,max The upper and lower limits of the output of the gas boiler are respectively; p (P) EC_inl,min 、P EC_inl,max The upper and lower limits of the output of the electric refrigerator are respectively; p (P) AC_inl,min 、P AC_inl,max The upper and lower limits of the output of the absorption refrigerator are respectively set.
2. The method for optimizing power for interactive operation of an integrated energy distribution system with an upper grid according to claim 1, wherein said base data comprises: the grid structure, the electric load level and the electric parameters of the comprehensive energy distribution system, the network topology, the air load level, the pipeline parameters of the air distribution system, the multi-energy load of the energy station and the equipment efficiency.
3. The power optimization method for interactive operation of the integrated energy distribution system and the upper power grid according to claim 1, wherein the expression of the objective function is:
Obj＝min|P actual -P expect | (1)
wherein P is actual To interact with the upper grid for actual power, P expect Power is desired for interaction with the upper grid.
4. The method for optimizing power of an integrated energy distribution system for interactive operation with a higher grid according to claim 1, wherein the adjustable range of the interactive power with the higher grid is the upper and lower limits of the interactive power with the higher grid.
5. The method for optimizing power for interactive operation of an integrated energy distribution system and an upper power grid according to claim 1, wherein the power optimization model for interactive operation of the integrated energy distribution system and the upper power grid is solved through an MATLAB platform and a YALMIP tool kit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010608816.8A CN111768036B (en) | 2020-06-29 | 2020-06-29 | Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010608816.8A CN111768036B (en) | 2020-06-29 | 2020-06-29 | Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111768036A CN111768036A (en) | 2020-10-13 |
| CN111768036B true CN111768036B (en) | 2023-11-03 |
Family
ID=72724107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010608816.8A Active CN111768036B (en) | 2020-06-29 | 2020-06-29 | Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111768036B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112488535A (en) * | 2020-12-01 | 2021-03-12 | 中国科学院地理科学与资源研究所 | Energy network security measure method, system and data platform |
| CN114548956B (en) * | 2022-04-26 | 2022-09-16 | 国网天津市电力公司电力科学研究院 | Integrated energy power distribution system based on honeycomb topology and operation method thereof |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105446163A (en) * | 2014-08-29 | 2016-03-30 | 国家电网公司 | Micro-grid versatile flow modeling approach based on different types of energy |
| CN106960282A (en) * | 2017-03-27 | 2017-07-18 | 武汉大学 | A kind of coordinated operation method of pneumoelectric comprehensive energy distribution network systems |
| CN107769215A (en) * | 2018-01-19 | 2018-03-06 | 国网天津市电力公司 | Garden energy mix system optimization dispatching method based on energy hub |
| CN107910881A (en) * | 2017-12-06 | 2018-04-13 | 重庆大学 | A kind of ADMM control methods based on network load contingency management |
| CN108898265A (en) * | 2018-05-02 | 2018-11-27 | 广州供电局有限公司 | A kind of integrated energy system integration planing method |
| CN108964041A (en) * | 2018-07-23 | 2018-12-07 | 燕山大学 | A kind of control method of electricity-interconnection integrated energy system peak load deliverability |
| CN110232217A (en) * | 2019-05-16 | 2019-09-13 | 天津大学 | A kind of comprehensive energy distribution system operation domain modeling method |
| CN110263387A (en) * | 2019-05-29 | 2019-09-20 | 国网福建省电力有限公司 | The integrated energy system optimal operation model and its linear processing methods that power distribution network, natural gas grid based on energy hub are established |
| CN110443398A (en) * | 2019-05-13 | 2019-11-12 | 华北电力大学（保定） | A kind of regional complex energy resource system optimizing operation method based on repeated game model |
| CN110458353A (en) * | 2019-08-08 | 2019-11-15 | 上海交通大学 | Electric heating combines microgrid energy step optimization method and system |
| CN110705779A (en) * | 2019-09-27 | 2020-01-17 | 河海大学 | Electric power-traffic network multi-period cooperative scheduling method considering traffic flow time domain coupling |
| WO2020093295A1 (en) * | 2018-11-06 | 2020-05-14 | 国电南瑞科技股份有限公司 | Multi-period power flow optimization method for electricity-gas interconnection integrated energy system |
-
2020
- 2020-06-29 CN CN202010608816.8A patent/CN111768036B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105446163A (en) * | 2014-08-29 | 2016-03-30 | 国家电网公司 | Micro-grid versatile flow modeling approach based on different types of energy |
| CN106960282A (en) * | 2017-03-27 | 2017-07-18 | 武汉大学 | A kind of coordinated operation method of pneumoelectric comprehensive energy distribution network systems |
| CN107910881A (en) * | 2017-12-06 | 2018-04-13 | 重庆大学 | A kind of ADMM control methods based on network load contingency management |
| CN107769215A (en) * | 2018-01-19 | 2018-03-06 | 国网天津市电力公司 | Garden energy mix system optimization dispatching method based on energy hub |
| CN108898265A (en) * | 2018-05-02 | 2018-11-27 | 广州供电局有限公司 | A kind of integrated energy system integration planing method |
| CN108964041A (en) * | 2018-07-23 | 2018-12-07 | 燕山大学 | A kind of control method of electricity-interconnection integrated energy system peak load deliverability |
| WO2020093295A1 (en) * | 2018-11-06 | 2020-05-14 | 国电南瑞科技股份有限公司 | Multi-period power flow optimization method for electricity-gas interconnection integrated energy system |
| CN110443398A (en) * | 2019-05-13 | 2019-11-12 | 华北电力大学（保定） | A kind of regional complex energy resource system optimizing operation method based on repeated game model |
| CN110232217A (en) * | 2019-05-16 | 2019-09-13 | 天津大学 | A kind of comprehensive energy distribution system operation domain modeling method |
| CN110263387A (en) * | 2019-05-29 | 2019-09-20 | 国网福建省电力有限公司 | The integrated energy system optimal operation model and its linear processing methods that power distribution network, natural gas grid based on energy hub are established |
| CN110458353A (en) * | 2019-08-08 | 2019-11-15 | 上海交通大学 | Electric heating combines microgrid energy step optimization method and system |
| CN110705779A (en) * | 2019-09-27 | 2020-01-17 | 河海大学 | Electric power-traffic network multi-period cooperative scheduling method considering traffic flow time domain coupling |
Non-Patent Citations (2)
| Title |
|---|
| 综合能源系统潮流及最优潮流计算模型与方法综述;陈思宇;柴庆宣;李延松;胡林献;;热力发电(第07期);全文 * |
| 面向能源互联网的多能流静态安全分析方法;潘昭光;孙宏斌;郭庆来;;电网技术(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111768036A (en) | 2020-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106786603B (en) | A kind of regional complex energy resource system multiobjective optimization mixed current algorithm | |
| CN108154309B (en) | Energy internet economic dispatching method considering multi-load dynamic response of cold, heat and electricity | |
| CN110110913B (en) | Optimal configuration method for energy station of large-scale park comprehensive energy system | |
| CN104751246B (en) | A kind of active power distribution network planing method based on Stochastic Chance-constrained | |
| CN104734155B (en) | A kind of method obtaining electrical interconnection energy resource system available transmission capacity | |
| CN109409705B (en) | Multi-objective optimization scheduling method for regional comprehensive energy system | |
| CN106786753B (en) | The system and its adjusting method of the Regional Energy internet of multi-user | |
| CN109063992A (en) | Consider the power distribution network Expansion Planning method of regional complex energy resource system optimization operation | |
| CN106096777B (en) | A kind of Optimization Scheduling of electric-gas coupling multipotency streaming system | |
| CN111768036B (en) | Power optimization method for interactive operation of comprehensive energy distribution system and superior power grid | |
| CN109217291A (en) | Consider the electrical interconnection system Multipurpose Optimal Method of peak load shifting model | |
| CN111952958A (en) | Power distribution network flexible soft switch optimal configuration method considering control mode conversion | |
| CN107196316A (en) | Multistage reactive voltage control method for coordinating in active distribution network | |
| CN106803130A (en) | Distributed power source accesses the planing method of power distribution network | |
| CN113570117A (en) | Electricity-gas comprehensive energy microgrid optimal scheduling method based on second-order cone relaxation conversion method | |
| CN106953361B (en) | Exchange micro-grid operational control method and device | |
| CN116402210A (en) | Multi-objective optimization method, system, equipment and medium for comprehensive energy system | |
| CN112836957B (en) | Regional comprehensive energy system planning method considering power supply reliability | |
| CN114077934A (en) | Comprehensive energy microgrid interconnection system and scheduling method thereof | |
| Mahdad et al. | Solving multi-objective optimal power flow problem considering wind-STATCOM using differential evolution | |
| CN109472716A (en) | A kind of production capacity end power distribution method multifactor based on the considerations of genetic algorithm | |
| CN111815111B (en) | Regional comprehensive energy expansion planning method considering pipeline risk level | |
| CN107317342A (en) | A kind of idle planning of distributing wind power plant and powerless control method | |
| Tai et al. | Design and optimisation of offshore grids in baltic sea for scenario year 2030 | |
| Balram et al. | Utilization of renewable energy sources in generation and distribution 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |