CN111768036A - Power optimization method for interactive operation of comprehensive energy power distribution system and upper-level power grid - Google Patents
Power optimization method for interactive operation of comprehensive energy power distribution system and upper-level power grid Download PDFInfo
- Publication number
- CN111768036A CN111768036A CN202010608816.8A CN202010608816A CN111768036A CN 111768036 A CN111768036 A CN 111768036A CN 202010608816 A CN202010608816 A CN 202010608816A CN 111768036 A CN111768036 A CN 111768036A
- Authority
- CN
- China
- Prior art keywords
- power
- energy
- distribution system
- node
- 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.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 104
- 230000002452 interceptive effect Effects 0.000 title claims abstract description 40
- 238000005457 optimization Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000003993 interaction Effects 0.000 claims description 23
- 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
- 230000008569 process Effects 0.000 claims description 2
- 230000033228 biological regulation Effects 0.000 abstract description 8
- 239000007789 gas Substances 0.000 description 39
- 238000010586 diagram Methods 0.000 description 15
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 4
- 239000003345 natural gas Substances 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 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
- 238000005265 energy consumption Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000020169 heat generation 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
- 238000010248 power generation Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/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 an integrated energy power distribution system and a superior 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 integrated energy power distribution system and the superior power grid by using an MATLAB platform and a YALMIP tool box, outputting related results and the like. The power optimization method for interactive operation of the comprehensive energy power distribution system and the superior power grid can find an energy distribution scheme of the energy station for power optimization of interactive operation with the superior power grid on the premise of ensuring each load requirement; considering the energy supply of the energy station from the perspective of meeting the regulation and control instruction of the upper-level power grid, establishing a power optimization model of interactive operation of the comprehensive energy power distribution system and the upper-level power grid, which takes the deviation between the actual power and the expected power interacted with the upper-level power grid as an objective function and takes power flow constraint of the power distribution system and the like as constraint conditions, and solving the model based on an MATLAB platform and a YALMIP tool box to determine the energy supply of the energy station.
Description
Technical Field
The invention belongs to the technical field of comprehensive energy power distribution systems, and particularly relates to a power optimization method for interactive operation of a comprehensive energy power distribution system and a superior power grid.
Background
Social development has a strong dependence on energy, and due to unfavorable factors such as exhaustion of fossil energy and pollution to the environment, people have to seek clean energy to replace the traditional fossil energy. With the continuous maturity of research in the aspects of central air conditioners, gas turbines and the like, the coupling degree of each energy system is also enhanced day by day, the traditional energy systems are independently planned and respectively operated, the coupling condition among the energy sources is not fully considered, an Integrated Energy Distribution System (IEDS) considers the coupling relationship among the energy sources on the basis of the traditional energy systems, links such as production, transmission, distribution and conversion of the energy sources are involved, the system is uniformly planned and processed, the whole system is optimized, and the utilization rate of the energy sources is improved. 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 improved reliability of the system. In addition, innovative technologies developed based on natural gas also shift the energy structure of the conventional power supply using electric energy as the main energy source, and research on IEDS for distributing energy using power and natural gas networks as the main energy source networks has attracted attention.
The energy systems are in different levels and have complex differences in space-time, and the characteristics bring challenges to the analysis and research of the multi-energy systems in the aspects of energy conversion, storage, distribution and the like. The energy station is an energy input and output multi-port model, can describe energy supply and energy demand in a highly abstract way, can adjust the energy supply of various devices by adjusting the interrelation 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 capable of integrating any number of energy carriers, thereby providing a high degree of flexibility in system modeling.
The upper-level power grid supplies electric energy to the IEDS, and in order to guarantee safe and stable operation of the upper-level power grid, the upper-level power grid can expect that the IEDS can consume more or less electric energy. The traditional power grid can achieve the purpose only by increasing or decreasing loads, but the adjustment of the loads relates to a user side, so that certain difficulty exists in achieving the purpose, and the IEDS can achieve the adjustment of the interaction power of the comprehensive energy power distribution system and the superior power grid by adjusting the distribution relation of the energy station, so that the requirement of the superior power grid on the energy consumption of the IEDS is met. However, related technologies in this aspect are less researched, and in the related research, the problems of safe operation of the natural gas network and feasible power regulation domain interaction between the integrated energy power distribution system and the upper-level power grid are less considered, and a situation that the difference from an expected result is large easily occurs.
Disclosure of Invention
In order to realize the adjustment of the interactive power of the comprehensive energy power distribution system and the superior power grid, the invention provides a power optimization method for interactive operation of the comprehensive energy power distribution system and the superior power grid.
The invention provides a power optimization method for interactive operation of an integrated energy power distribution system and a superior power grid, which comprises the following processes:
s1, collecting basic data of the comprehensive energy distribution system to be researched;
s2, establishing a power optimization model of interactive operation of the comprehensive energy power distribution system and a superior power grid through basic data, and taking the difference between the actual value and the expected value of the interactive power of the comprehensive energy power distribution system and the superior power grid as a target function;
and S3, determining constraint conditions, wherein the set constraint conditions comprise: the method comprises the following steps of power distribution system flow constraint, operation voltage level constraint, branch current constraint, gas distribution system pipeline flow constraint, air 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 power distribution system and a superior power grid, and setting an expected value of the interaction power of the comprehensive energy power distribution system and the superior power grid according to the adjustable range;
and S5, solving the power optimization model of the interactive operation of the comprehensive energy power distribution system and the superior power grid to obtain the actual interactive power value of the comprehensive energy power distribution system and the superior power grid and the corresponding value of the energy flow balance relation of each device of the energy source station.
Preferably, the base data comprises: the comprehensive energy distribution system comprises a grid structure, an electric load level and electric parameters of the comprehensive energy distribution system, and a network topology, an air load level, pipeline parameters, a multi-energy load of an energy station and equipment efficiency of the gas distribution system.
Preferably, the expression of the objective function is:
Obj＝min|Pactual-Pexpect| (1)
wherein, PactualFor interaction of actual power with the superordinate grid, PexpectThe desired power is exchanged with the superordinate network.
Preferably, the expression of the power flow constraint of the power distribution system is as follows:
wherein i ∈ u (j) is a set adjacent to node j and node i is a first node of the line, k ∈ v (j) is a set adjacent to node j and node k is a last node of the line, PijActive power flow for the head end of the line; pjkActive power flow for the end of the line; pjNet injection of active power for node j; qijIs the reactive power flow of the head end of the line; qjkIs the reactive power flow at the end of the line; qjNet injection of reactive power for node j; vi、VjRespectively the voltage amplitudes of the head end and the tail end of the circuit; i isijIs the magnitude of the current flowing through the line; r isijAnd xijRespectively the resistance value and the reactance value of the circuit;
the operating voltage level constraint is expressed as:
Vi，min≤Vi≤Vi，max(3)
wherein, ViIs the voltage of node i, Vi,maxAnd Vi,minRespectively the upper limit and the lower limit of the voltage amplitude of the node i;
the expression of the branch current limit is:
Iij≤Iij,max(4)
wherein, IijFor the amplitude of the current flowing through the line, Iij,maxThe maximum allowable value of the current flowing through the line;
the flow constraint expression 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 friction coefficient; q. q.smnThe flow rate flowing through the pipeline; p is a radical ofmAnd pnRespectively taking the pressure values of the head node and the tail node of the pipeline; the matrix A is a branch node incidence matrix, aijThe value of (a) depends on the relationship 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 branch j are associated and the gas flow in the pipeline flows out of the node, the value is-1; qbranchVector composed of branch flow; qnodeVector composed of flow load of each node;
the expression for the barometric pressure level constraint is:
pm，min≤pm≤pm,max(6)
wherein p ismIs a nodal pressure value, pm,maxAnd pm,minRespectively representing the upper limit and the lower limit of the node pressure value;
the expression of the pipeline flow constraint is as follows:
qmn，min≤qmn≤qmn,max(7)
wherein q ismnIs the value of the pipe flow, qmn,maxAnd q ismn,minRespectively the upper limit and the lower limit of the pipeline flow value;
the energy station energy flow constraint expression is as follows:
wherein L isel、Lhl、LclRespectively electric, hot and cold loads; pelPower interaction between energy station and distribution system ηGT_el、ηGT_hl、ηGBl、ηECl、ηAClRespectively the efficiency of the gas turbine for generating electric energy, the efficiency of the gas turbine for generating heat, the efficiency of a gas boiler, the efficiency of an electric refrigerator and the efficiency of an absorption refrigerator; pGT_inl、PGT_el、PGT_hlPower for gas input, electrical output, thermal output, respectively, of the gas turbine; pGB_inl、PGB_outlThe power of gas input and heat output of the gas boiler respectively; pEC_inl、PEC_outlPower of electric input and cold output of the electric refrigerator respectively; pAC_inl、PAC_outlThe power of the heat input and the cold output of the absorption refrigerator respectively;
the expression of the equipment output constraint is as follows:
wherein, PGT_inl,min、PGT_inl,maxRespectively the upper and lower output limits of the gas turbine; pGB_inl,min、PGB_inl,maxRespectively the upper and lower output limits of the gas boiler; pEC_inl,min、PEC_inl,maxRespectively the upper and lower output limits of the electric refrigerator; pAC_inl,min、PAC_inl,maxRespectively the upper and lower limits of the output of the absorption refrigerator.
Preferably, the adjustable range of the power interacting with the upper-level power grid is the upper limit and the lower limit of the power interacting with the upper-level power grid.
Preferably, the power optimization model of the interactive operation of the comprehensive energy power distribution system and the superior power grid is solved through an MATLAB platform and a YALMIP toolbox.
The invention provides a power optimization method for interactive operation of an integrated energy power distribution system and a superior power grid, which models the power optimization problem of interactive operation of the integrated energy power distribution system and the superior power grid, takes the deviation of the actual power and the expected power of the interaction of the integrated energy power distribution system and the superior power grid as a target function, and takes power flow constraint, operation voltage level constraint, branch current constraint, gas distribution system power flow constraint, air pressure level constraint, pipeline flow constraint, energy station energy flow constraint, equipment output constraint and the like of the power distribution system as constraint conditions. The method is different from the previous method for optimizing the interaction power with the superior power grid in that the method provided by the invention is realized through energy distribution of the energy station, firstly, a reference range of a regulation instruction is provided for the superior power grid according to a feasible regulation domain of the interaction power of the comprehensive energy power distribution system and the superior power grid, and then the power optimization of the interaction operation of the comprehensive energy power distribution system and the superior power grid is realized by taking the minimum deviation between the actual power and the expected power of the interaction of the comprehensive energy power distribution system and the superior power grid as a target function. In the aspect of model establishment, the method comprehensively considers the models of the power distribution system and the gas distribution system, and introduces the model of the energy station to describe energy supply and energy demand.
Drawings
FIG. 1a is a schematic diagram of an IEEE33 node power distribution system in an example of the IEDS algorithm;
FIG. 1b is a schematic diagram of an 11-node gas distribution system in an IEDS calculation example;
FIG. 1c is a schematic diagram of an energy station 1 in an IEDS algorithm;
FIG. 1d is a schematic diagram of an energy station 2 in an IEDS algorithm;
FIG. 1e is a schematic diagram of two energy stations in an example of IEDS calculation;
fig. 2 is a flow chart of a power optimization method for interactive operation of the integrated energy power distribution system and a superior power grid according to the present invention;
FIG. 3 shows the upper and lower limits of the interactive power regulation of the integrated energy power distribution system and the upper power grid;
FIG. 4 is a diagram illustrating the balance of output and consumed electric power of each device of the optimized energy station 1;
FIG. 5 is a diagram illustrating the balance of the output power and the consumed thermal power of each device of the optimized energy station 1;
FIG. 6 is a diagram illustrating the balance of the output and the consumed cold power of each device of the optimized energy station 1;
FIG. 7 is a diagram illustrating the balance of output and consumed electric power for optimizing the devices of the resulting energy station 2;
FIG. 8 is a graph of the thermal power balance of the output and consumption of each device of the optimized resulting energy station 2;
FIG. 9 is a diagram illustrating the balance of output and consumed cold power for each device of the resulting energy station 2;
fig. 10 is a comparison of actual values and expected values obtained from power optimization of the integrated energy distribution system operating interactively with the upper grid.
Detailed Description
The invention provides a power optimization method for interactive operation of an integrated energy power distribution system and a superior power grid, which can find an energy station energy distribution scheme for power optimization of interactive operation of the integrated energy power distribution system and the superior power grid on the premise of ensuring each load requirement. The invention considers the upper-level power grid regulation and control instruction to optimize the interactive power of the comprehensive energy power distribution system and the upper-level power grid, establishes a power optimization model of interactive operation of the comprehensive energy power distribution system and the upper-level power grid, which takes the deviation between the actual power and the expected power of the interaction of the minimized comprehensive energy power distribution system and the upper-level power grid as an objective function and takes power flow constraint, operation voltage level constraint, branch current constraint, distribution system power flow constraint, air pressure level constraint, pipeline flow constraint, energy station energy flow constraint, equipment output constraint and the like as constraint conditions, and solves the model based on an MATLAB platform and a YALMIP tool box to determine the energy flow balance relation of the energy station.
The method for optimizing power of interactive operation between the integrated energy power distribution system and the upper-level power grid according to the present invention is further described in detail by referring to the preferred embodiments in conjunction with the accompanying drawings.
An IEDS example in the embodiment of the present invention includes an IEEE33 node power distribution system, an 11 node air distribution system, and two energy stations, where fig. 1a is a schematic diagram of an IEEE33 node power distribution system in the IEDS example, fig. 1b is a schematic diagram of an 11 node air distribution system in the IEDS example, fig. 1c is a schematic diagram of an energy station 1 in the IEDS example, fig. 1d is a schematic diagram of an energy station 2 in the IEDS example, and fig. 1e is a schematic diagram of structures of two energy stations in the IEDS example.
First, basic data of the IEDS algorithm shown in fig. 1, including a grid structure, an electrical load level, and electrical parameters of the integrated energy distribution system, and basic data of a network topology, an air load level, pipeline parameters, a multi-energy load of an energy station, and equipment efficiency of the gas distribution system, are obtained. The access position of the energy station is set, in the embodiment, the 6 nodes and 13 nodes of the power distribution system are used as the electric energy input positions of the energy station, and the N4 nodes and N9 nodes of the gas distribution system are used as the natural gas energy input positions of the energy station. The reference voltage of the system was set at 12.66kV and the gas pressure at source point N1 was set at 75 mbar.
Secondly, establishing a power optimization model of interactive operation of the comprehensive energy power distribution system and a superior power grid by using the acquired basic data, and setting a difference value between an actual power value and an expected power value of the interaction of the comprehensive energy power distribution system and the superior power grid as an objective function, wherein the expression of the objective function is as follows:
Obj＝min|Pactual-Pexpect| (1)
wherein, PactualFor actual power interaction with the superordinate grid, PexpectIs the desired power to interact with the upper grid. In the embodiment of the invention, the power of the interaction between the comprehensive energy power distribution system and the upper-level power grid is the power value of the power distribution system node 1 connected with the upper-level 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, air pressure level constraint, pipeline flow constraint, energy station energy flow constraint and equipment output constraint. The power distribution system power flow constraint expression is as follows:
wherein i ∈ u (j) is a set adjacent to node j and node i is a first node of the line, k ∈ v (j) is a set adjacent to node j and node k is a last node of the line, PijActive power flow for the head end of the line; pjkActive power flow for the end of the line; pjNet injection of active power for node j; qijIs the reactive power flow of the head end of the line; qjkIs the reactive power flow at the end of the line; qjNet injection of reactive power for node j; vi、VjRespectively the voltage amplitudes of the head end and the tail end of the circuit; i isijIs the magnitude of the current flowing through the line; r isijAnd xijRespectively the resistance value and the reactance value of the line.
The operating voltage level constraint is expressed as:
Vi，min≤Vi≤Vi，max(3)
wherein, ViIs the voltage of node i, Vi,maxAnd Vi,minRespectively the upper and lower limits of the voltage amplitude of the node i.
The expression of the branch current limit is:
Iij≤Iij,max(4)
wherein, IijFor the amplitude of the current flowing through the line, Iij,maxThe maximum allowable value of the current flowing through the line.
The flow constraint expression 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 friction coefficient; qmn is the flow rate in the pipeline; p is a radical ofmAnd pnRespectively taking the pressure values of the head node and the tail node of the pipeline; the matrix A is a branch node incidence matrix, aijThe value of (a) depends on the relationship 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 branch j are associated and the gas flow in the pipeline flows out of the node, the value is-1; qbranchVector composed of branch flow; qnodeAnd forming a vector for the traffic load of each node.
The expression for the barometric pressure level constraint is:
pm，min≤pm≤pm,max(6)
wherein p ismIs a nodal pressure value, pm,maxAnd pm,minRespectively an upper limit and a lower limit of the node pressure value.
The expression of the pipeline flow constraint is as follows:
qmn，min≤qmn≤qmn,max(7)
wherein q ismnIs the value of the pipe flow, qmn,maxAnd q ismn,minRespectively the upper and lower limits of the pipeline flow value.
The energy station energy flow constraint expression is as follows:
wherein L isel、Lhl、LclRespectively electric, hot and cold loads; pelPower interaction between energy station and distribution system ηGT_el、ηGT_hl、ηGBl、ηECl、ηAClRespectively the efficiency of the gas turbine for generating electric energy, the efficiency of the gas turbine for generating heat, the efficiency of a gas boiler, the efficiency of an electric refrigerator and the efficiency of an absorption refrigerator; pGT_inl、PGT_el、PGT_hlPower for gas input, electrical output, thermal output, respectively, of the gas turbine; pGB_inl、PGB_outlThe power of gas input and heat output of the gas boiler respectively; pEC_inl、PEC_outlPower of electric input and cold output of the electric refrigerator respectively; pAC_inl、PAC_outlRespectively of absorption refrigeratorsHeat input, cold output power.
The expression of the equipment output constraint is as follows:
wherein, PGT_inl,min、PGT_inl,maxRespectively the upper and lower output limits of the gas turbine; pGB_inl,min、PGB_inl,maxRespectively the upper and lower output limits of the gas boiler; pEC_inl,min、PEC_inl,maxRespectively the upper and lower output limits of the electric refrigerator; pAC_inl,min、PAC_inl,maxRespectively the upper and lower limits of the output of the absorption refrigerator.
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 interaction power between the comprehensive energy power distribution system and the upper-level power grid, namely upper and lower limits of the interaction power between the comprehensive energy power distribution system and the upper-level power grid, and providing reference for setting an expected value of the interaction power between the comprehensive energy power distribution system and the upper-level power grid. In the embodiment of the present invention, the result is shown in fig. 3, namely, the result is to determine the upper and lower limit values of the power value of the power distribution system node 1 under the condition that each constraint condition is satisfied.
And finally, solving the power optimization model of the interactive operation of the comprehensive energy power distribution system and the superior power grid by using an MATLAB platform and a YALMIP tool box, and outputting related results including actual values of the interactive power of the comprehensive energy power distribution system and the superior power grid, corresponding values of energy flow balance relations of all equipment of the energy station and the like.
Fig. 4 to 10 show the optimized configuration results of the embodiment of the present invention, where fig. 4, fig. 5, and fig. 6 respectively show the power output and consumed power, heat, and cold power balance relationship of each device of the energy station 1, fig. 7, fig. 8, and fig. 9 respectively show the power output and consumed power, heat, and cold power balance relationship of each device of the energy station 2, and fig. 10 shows the actual power and expected power of the integrated energy power distribution system interacting with the upper-level power grid.
According to the obtained result, the power optimization model for interactive operation of the comprehensive energy power distribution system and the superior power grid has rationality, so that the characteristic of 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 power distribution system and the upper-level power grid can make the actual power interacted with the upper-level power grid very close to the expected power, so as to meet the requirement of the regulation and control instruction of the upper-level power grid. Fig. 4-9 may provide a configuration scheme for optimizing the output of each device of the energy source station. The power optimization method for the interactive operation of the comprehensive energy power distribution system and the superior power grid can provide reference for the superior power grid to issue 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 with reference to the preferred embodiments, it should be understood that the above description should not be taken as limiting the invention. Various modifications and alterations to this invention will become apparent to those skilled in the art upon reading the foregoing description. Accordingly, the scope of the invention should be determined from the following claims.
Claims (6)
1. A power optimization method for interactive operation of an integrated energy power distribution system and a superior power grid is characterized by comprising the following processes:
s1, collecting basic data of the comprehensive energy distribution system to be researched;
s2, establishing a power optimization model of interactive operation of the comprehensive energy power distribution system and a superior power grid through basic data, and taking the difference between the actual value and the expected value of the interactive power of the comprehensive energy power distribution system and the superior power grid as a target function;
and S3, determining constraint conditions, wherein the set constraint conditions comprise: the method comprises the following steps of power distribution system flow constraint, operation voltage level constraint, branch current constraint, gas distribution system pipeline flow constraint, air 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 power distribution system and a superior power grid, and setting an expected value of the interaction power of the comprehensive energy power distribution system and the superior power grid according to the adjustable range;
and S5, solving the power optimization model of the interactive operation of the comprehensive energy power distribution system and the superior power grid to obtain the actual interactive power value of the comprehensive energy power distribution system and the superior power grid and the corresponding value of the energy flow balance relation of each device of the energy source station.
2. The method of claim 1, wherein the basic data comprises: the comprehensive energy distribution system comprises a grid structure, an electric load level and electric parameters of the comprehensive energy distribution system, and a network topology, an air load level, pipeline parameters, a multi-energy load of an energy station and equipment efficiency of the gas distribution system.
3. The method of claim 1, wherein the objective function is expressed as:
Obj＝min|Pactual-Pexpect| (1)
wherein, PactualFor interaction of actual power with the superordinate grid, PexpectThe desired power is exchanged with the superordinate network.
4. The method for optimizing power of interactive operation of an integrated energy power distribution system and a superior power grid according to claim 1, wherein the power flow constraint of the power distribution system is expressed by:
wherein i ∈ u (j) is a set adjacent to node j and node i is a first node of the line, k ∈ v (j) is a set adjacent to node j and node k is a last node of the line, PijActive power flow for the head end of the line; pjkActive power flow for the end of the line; pjNet injection of active power for node j; qijReactive power flow for line head end；QjkIs the reactive power flow at the end of the line; qjNet injection of reactive power for node j; vi、VjRespectively the voltage amplitudes of the head end and the tail end of the circuit; i isijIs the magnitude of the current flowing through the line; r isijAnd xijRespectively the resistance value and the reactance value of the circuit;
the operating voltage level constraint is expressed as:
Vi，min≤Vi≤Vi，max(3)
wherein, ViIs the voltage of node i, Vi,maxAnd Vi,minRespectively the upper limit and the lower limit of the voltage amplitude of the node i; the expression of the branch current limit is:
Iij≤Iij,max(4)
wherein, IijFor the amplitude of the current flowing through the line, Iij,maxThe maximum allowable value of the current flowing through the line; the flow constraint expression 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 friction coefficient; q. q.smnThe flow rate flowing through the pipeline; p is a radical ofmAnd pnRespectively taking the pressure values of the head node and the tail node of the pipeline; the matrix A is a branch node incidence matrix, aijThe value of (a) depends on the relationship 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 branch j are associated and the gas flow in the pipeline flows out of the node, the value is-1; qbranchVector composed of branch flow; qnodeVector composed of flow load of each node;
the expression for the barometric pressure level constraint is:
pm，min≤pm≤pm,max(6)
wherein p ismIs a nodal pressure value, pm,maxAnd pm,minRespectively representing the upper limit and the lower limit of the node pressure value;
the expression of the pipeline flow constraint is as follows:
qmn，min≤qmn≤qmn,max(7)
wherein q ismnIs the value of the pipe flow, qmn,maxAnd q ismn,minRespectively the upper limit and the lower limit of the pipeline flow value;
the energy station energy flow constraint expression is as follows:
wherein L isel、Lhl、LclRespectively electric, hot and cold loads; pelPower interaction between energy station and distribution system ηGT_el、ηGT_hl、ηGBl、ηECl、ηAClRespectively the efficiency of the gas turbine for generating electric energy, the efficiency of the gas turbine for generating heat, the efficiency of a gas boiler, the efficiency of an electric refrigerator and the efficiency of an absorption refrigerator; pGT_inl、PGT_el、PGT_hlPower for gas input, electrical output, thermal output, respectively, of the gas turbine; pGB_inl、PGB_outlThe power of gas input and heat output of the gas boiler respectively; pEC_inl、PEC_outlPower of electric input and cold output of the electric refrigerator respectively; pAC_inl、PAC_outlAre respectively of absorption type
Power of heat input and cold output of the refrigerator;
the expression of the equipment output constraint is as follows:
wherein, PGT_inl,min、PGT_inl,maxRespectively the upper and lower output limits of the gas turbine; pGB_inl,min、PGB_inl,maxRespectively the upper and lower output limits of the gas boiler; pEC_inl,min、PEC_inl,maxRespectively the upper and lower output limits of the electric refrigerator; pAC_inl,min、PAC_inl,maxRespectively the upper and lower limits of the output of the absorption refrigerator.
5. The method as claimed in claim 1, wherein the adjustable range of the interactive power with the upper grid is upper and lower limits of the interactive power with the upper grid.
6. The method according to claim 1, wherein the power optimization model for the interactive operation of the integrated energy power distribution system and the upper power grid is solved through a MATLAB platform and a yalmap toolkit.
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 true CN111768036A (en) | 2020-10-13 |
| CN111768036B 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) |
Cited By (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 |
| CN114548956A (en) * | 2022-04-26 | 2022-05-27 | 国网天津市电力公司电力科学研究院 | 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 |
|---|
| 潘昭光;孙宏斌;郭庆来;: "面向能源互联网的多能流静态安全分析方法", 电网技术, no. 06 * |
| 陈思宇;柴庆宣;李延松;胡林献;: "综合能源***潮流及最优潮流计算模型与方法综述", 热力发电, no. 07 * |
Cited By (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 |
| CN114548956A (en) * | 2022-04-26 | 2022-05-27 | 国网天津市电力公司电力科学研究院 | Integrated energy power distribution system based on honeycomb topology and operation method thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111768036B (en) | 2023-11-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106786603B (en) | A kind of regional complex energy resource system multiobjective optimization mixed current algorithm | |
| CN105576657B (en) | Random fuzzy continuous power flow modeling method for active power distribution network and power transmission network integrated system | |
| JP2022550225A (en) | Electric heat pump - regulation method and system for optimizing cogeneration systems | |
| CN109063992A (en) | Consider the power distribution network Expansion Planning method of regional complex energy resource system optimization operation | |
| CN108565863A (en) | A kind of regional complex energy resource system multiple target tide optimization method considering randomness | |
| CN109102104B (en) | Multi-energy system collaborative low-carbon planning method based on carbon emission flow | |
| CN109409705B (en) | Multi-objective optimization scheduling method for regional comprehensive energy system | |
| CN108304972B (en) | Active power distribution network frame planning method based on supply and demand interaction and DG (distributed generation) operation characteristics | |
| CN108023364A (en) | Power distribution network distributed generation resource maximum access capability computational methods based on convex difference planning | |
| CN112035984B (en) | Collaborative planning method for comprehensive energy system of electricity-gas-storage area | |
| CN109217291A (en) | Consider the electrical interconnection system Multipurpose Optimal Method of peak load shifting model | |
| CN111768036A (en) | Power optimization method for interactive operation of comprehensive energy power distribution system and upper-level power grid | |
| CN114077934A (en) | Comprehensive energy microgrid interconnection system and scheduling method thereof | |
| CN113240204A (en) | Energy station capacity optimal configuration method and system considering renewable energy consumption area | |
| CN109301882A (en) | A kind of electric system distributing economic load dispatching method based on ADMM | |
| CN108594658A (en) | A kind of electric-gas coupled system maximum probability load margin Model for Multi-Objective Optimization and its method for solving | |
| CN110535120A (en) | Consider the interconnection system distributed Optimization Scheduling of electric-gas of air pollution diffusion | |
| CN116707023A (en) | Active power distribution network layering and partitioning comprehensive optimization method based on source-load correlation clustering | |
| Abdelaziz et al. | Energy hub optimization using modified firefly algorithm | |
| CN111815111B (en) | Regional comprehensive energy expansion planning method considering pipeline risk level | |
| CN109390970B (en) | Island microgrid distributed control method and system based on multi-Agent communication network | |
| CN107069703A (en) | A kind of alternating current-direct current distribution planing method of meter and new energy access | |
| CN112288216A (en) | Cooperative game-based capacity planning method and system for electric gas conversion device | |
| Kou | Particle Swarm Optimization Based Reactive Power Dispatch for Power Networks with Distributed Generation | |
| CN105262107A (en) | Wind power grid connected reactive optimization method |
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 |