CN112036613A - Park comprehensive energy optimization method and device based on ADMM alternating direction multiplier method - Google Patents

Park comprehensive energy optimization method and device based on ADMM alternating direction multiplier method Download PDF

Info

Publication number
CN112036613A
CN112036613A CN202010823959.0A CN202010823959A CN112036613A CN 112036613 A CN112036613 A CN 112036613A CN 202010823959 A CN202010823959 A CN 202010823959A CN 112036613 A CN112036613 A CN 112036613A
Authority
CN
China
Prior art keywords
optimization
natural gas
emergency
admm
output
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010823959.0A
Other languages
Chinese (zh)
Inventor
宋艺航
尚楠
黄国日
王刚
陈政
陈晖�
辜炜德
张翔
李沛
冷媛
杨再敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Energy Development Research Institute of China Southern Power Grid Co Ltd
Original Assignee
Energy Development Research Institute of China Southern Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Energy Development Research Institute of China Southern Power Grid Co Ltd filed Critical Energy Development Research Institute of China Southern Power Grid Co Ltd
Priority to CN202010823959.0A priority Critical patent/CN112036613A/en
Publication of CN112036613A publication Critical patent/CN112036613A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/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
    • 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
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P80/00—Climate change mitigation technologies for sector-wide applications
    • Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (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)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a park comprehensive energy optimization method and a device based on an ADMM alternating direction multiplier method, wherein the method comprises the following steps: acquiring declaration information of participants, and displaying real-time demand and output information on a public interface; according to the declaration information, performing pre-optimization processing by adopting an ADMM alternative direction multiplier method, decomposing two optimization sub-problems of an electric power system and a natural gas system by using the electric power and energy system, generating optimization parameters and optimization results by taking the minimum operation cost of the energy system of the whole park as a target, and sending the optimization parameters and the optimization results to enterprises in the park; performing safety check according to the optimized parameters; and if an emergency occurs, sending an emergency notice according to the emergency, and performing emergency optimization processing. The method can solve the problems of opaque information, complex model and the like in the centralized unified optimization method of the energy system, and realizes the cooperative optimization of the comprehensive energy system in the park through a small amount of information interaction between the power system and the natural gas system.

Description

Park comprehensive energy optimization method and device based on ADMM alternating direction multiplier method
Technical Field
The invention relates to the technical field of electric power systems, in particular to a park comprehensive energy optimization method and device based on an ADMM alternating direction multiplier method.
Background
With the improvement of environmental protection requirements and the development of clean energy technologies, research methods for integrating various energy sources such as natural gas and electric energy in areas and realizing coordinated planning and optimized operation among different energy sources are not developed in the background of comprehensive energy systems.
In the traditional optimized scheduling of the power-natural gas system, the scheduling of the power system and the scheduling of the natural gas system are not consistent, but have respective unique modes, and the interaction and coupling relationship between the power system and the natural gas system is often ignored, so that an over-ideal optimized result can be caused. In the current optimization research for the power-natural gas system, the power system and the natural gas system are generally considered to be operated uniformly by a joint dispatching center, which ignores the influence of information opaqueness which may be caused by different operators of the power system and the natural gas system. In case of emergency, the energy supplier and each participant lack a scheduling method capable of handling the emergency in time. In summary, for an integrated energy system with non-intercommunicating information, a distributed optimization method and a means for effectively handling emergency need to be provided.
Disclosure of Invention
The purpose of the invention is: the utility model provides a park energy integration optimization method and device based on ADMM alternative direction multiplier method, can solve the problem such as information opaqueness, model complicacy in the centralized unified optimization method of energy system, through a small amount of information interaction between electric power system and natural gas system, has realized the collaborative optimization of the energy integration system in the park.
In order to achieve the above object, the present invention provides a park integrated energy optimization method based on an ADMM alternative direction multiplier method, including:
acquiring declaration information of participants, and displaying real-time demand and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant;
according to the declaration information, performing pre-optimization processing by adopting an ADMM alternative direction multiplier method, decomposing two optimization sub-problems of an electric power system and a natural gas system by using the electric power and energy system, generating optimization parameters and optimization results by taking the minimum operation cost of the energy system of the whole park as a target, and sending the optimization parameters and the optimization results to enterprises in the park;
according to the optimized parameters, each enterprise in the park arranges the output and the demand of the next day and checks the safety; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check;
if an emergency occurs, sending an emergency notice according to the emergency, and performing emergency optimization processing; wherein the condition comprises: emergency forces, demand changes, and line faults.
Further, the power and energy system is decomposed into two optimization sub-problems of a power system and a natural gas system; the method specifically comprises the following steps: in the decomposition process, a variable which is common to an electric power system and a natural gas system and is the natural gas consumption of a gas turbine set is selected as a shared variable, and in the natural gas system, the mathematical expression of the consumption characteristic of the gas turbine set is as follows:
Figure BDA0002634148120000021
wherein a, b and c are consumption coefficients of the gas turbine unit; pGNActive power output by the gas turbine unit; gGNThe consumption of natural gas;
in the power system optimization sub-problem, inGNRepresenting the consumption of natural gas, the following mathematical expression is satisfied:
fGN＝gGN
according to the formula, the energy optimization of the electric power and natural gas system is carried out based on the ADMM alternative direction multiplier method.
Further, the mathematical expression of the power system optimization sub-problem objective function is as follows:
Figure BDA0002634148120000031
in the formula, WeIs the power system fuel price;
Figure BDA0002634148120000032
the consumption of the generator set at the minimum output is obtained; kGThe number of sections of the output of the generator set is represented; mjThe slope of the j-th piecewise linearization of the generator set and the gas generator set; delta PjActive output is the j section of the generator set; y is a Lagrange multiplier variable; rho is a constant penalty factor;
the constraint condition comprises system active power balance constraint; the unit is subjected to piecewise linearization output constraint; line transmission power constraints; and (4) carrying out piecewise linearization constraint on the consumption function of the gas turbine.
Further, the mathematical expression of the natural gas system optimization subproblem objective function is as follows:
Figure BDA0002634148120000033
in the formula, WgRepresenting the price of the gas source; gpIndicating the steam injection quantity of the gas source;
the constraint conditions comprise upper and lower limit constraints of gas sources and gas load steam injection quantity; the method comprises the steps of (1) performing piecewise linearization constraint on a natural gas pipeline airflow function; gas transmission pipeline node air pressure constraint; node air supply balance constraint; the air pressure ratio of the air compressor is restrained by compensating air pressure loss; and (5) natural gas network incidence matrix constraint.
Further, according to the declaration information, performing pre-optimization processing by using an ADMM alternating direction multiplier method, specifically:
setting an initial value comprising the original residual r(0)Dual residual s(0)Penalty factor rho, shared variable
Figure BDA0002634148120000034
Lagrange multiplier coefficient
Figure BDA0002634148120000035
Sequentially and respectively calculating an electric power energy flow optimization subproblem and a natural gas energy flow optimization subproblem, taking the kth iteration as an example: firstly, calculating a power energy flow optimization sub-problem, wherein a mathematical expression is as follows:
Figure BDA0002634148120000036
secondly, calculating a natural gas energy flow optimization sub-problem, wherein a mathematical expression is as follows:
Figure BDA0002634148120000041
and updating the Lagrange multiplier variable y after the kth iteration is completed, wherein the updating formula is as follows:
Figure BDA0002634148120000042
and judging convergence;
judging the convergence, and if the convergence is reached, outputting a pre-optimization result; if not, entering k +1 times of iteration, and judging a convergence formula as follows:
Figure BDA0002634148120000043
Figure BDA0002634148120000044
in the formula, r(k+1)And s(k+1)Respectively an original residual error and a dual residual error after the (k + 1) th iteration;prianddualrespectively the corresponding maximum allowed deviation;
and performing alternate iteration on the power energy flow optimization subproblem and the natural gas energy flow optimization subproblem until convergence, and outputting a pre-optimization result.
Further, if an emergency occurs, sending an emergency notice according to the emergency, and performing emergency optimization operation; wherein the condition comprises: emergency forces, demand changes, and line faults; the method specifically comprises the following steps: if the emergency output or the requirement changes, the enterprise directly issues emergency announcements, an operator and other enterprises receive the announcements at the same time, and the operator performs emergency optimization processing; if a line fault occurs, an operator directly issues an emergency notice and performs pre-optimization processing according to the ADMM alternative direction multiplier method again.
The embodiment of the invention also provides a park comprehensive energy optimization device based on the ADMM alternative direction multiplier method, which comprises the following steps: the system comprises an information acquisition unit, a pre-optimization processing unit, a safety check unit and an emergency processing unit, wherein the information acquisition unit is used for acquiring information;
the information acquisition unit is used for acquiring declaration information of participants and displaying real-time requirement and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant;
the pre-optimization processing unit is used for performing pre-optimization processing by adopting an ADMM alternative direction multiplier method according to the declaration information, decomposing the electric power and energy system into two optimization sub-problems of an electric power system and a natural gas system, generating optimization parameters and optimization results by taking the minimum operation cost of the energy system of the whole park as a target, and sending the optimization parameters and the optimization results to each enterprise of the park;
the safety checking unit is used for carrying out output and demand arrangement on each enterprise in the park on the next day according to the optimized parameters and carrying out safety checking; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check;
the emergency processing unit is used for sending an emergency notice according to the emergency and carrying out emergency optimization processing if the emergency occurs; wherein the condition comprises: emergency forces, demand changes, and line faults.
Further, the power and energy system is decomposed into two optimization sub-problems of a power system and a natural gas system; the method specifically comprises the following steps: in the decomposition process, a variable which is common to an electric power system and a natural gas system and is the natural gas consumption of a gas turbine set is selected as a shared variable, and in the natural gas system, the mathematical expression of the consumption characteristic of the gas turbine set is as follows:
Figure BDA0002634148120000051
wherein a, b and c are consumption coefficients of the gas turbine unit; pGNActive power output by the gas turbine unit; gGNIs the natural gas consumption.
In the power system optimization sub-problem, inGNRepresenting the consumption of natural gas, the following mathematical expression is satisfied:
fGN＝gGN
according to the formula, the energy optimization of the electric power and natural gas system is carried out based on the ADMM alternative direction multiplier method.
Further, the mathematical expression of the power system optimization sub-problem objective function is as follows:
Figure BDA0002634148120000061
in the formula, WeIs the power system fuel price;
Figure BDA0002634148120000062
the consumption of the generator set at the minimum output is obtained; kGThe number of sections of the output of the generator set is represented; mjThe slope of the j-th piecewise linearization of the generator set and the gas generator set; delta PjActive output is the j section of the generator set; y is a Lagrange multiplier variable; rho is a constant penalty factor;
the constraint condition comprises system active power balance constraint; the unit is subjected to piecewise linearization output constraint; line transmission power constraints; the gas turbine consumption function is subjected to piecewise linearization constraint;
the mathematical expression of the natural gas system optimization subproblem objective function is as follows:
Figure BDA0002634148120000063
in the formula, WgRepresenting the price of the gas source; gpIndicating the steam injection quantity of the gas source;
the constraint conditions comprise upper and lower limit constraints of gas sources and gas load steam injection quantity; the method comprises the steps of (1) performing piecewise linearization constraint on a natural gas pipeline airflow function; gas transmission pipeline node air pressure constraint; node air supply balance constraint; the air pressure ratio of the air compressor is restrained by compensating air pressure loss; and (5) natural gas network incidence matrix constraint.
An embodiment of the present invention further provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the method for optimizing the park integrated energy based on the ADMM alternative direction multiplier method according to any one of claims 1 to 6 is implemented.
Compared with the prior art, the method and the device for optimizing the comprehensive energy of the park based on the ADMM alternative direction multiplier method have the beneficial effects that:
the invention discloses a park comprehensive energy optimization method based on an ADMM alternating direction multiplier method, which comprises the following steps: acquiring declaration information of participants, and displaying real-time demand and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant; according to the declaration information, performing pre-optimization processing by adopting an ADMM alternative direction multiplier method, decomposing two optimization sub-problems of an electric power system and a natural gas system by using the electric power and energy system, generating optimization parameters and optimization results by taking the minimum operation cost of the energy system of the whole park as a target, and sending the optimization parameters and the optimization results to enterprises in the park; according to the optimized parameters, each enterprise in the park arranges the output and the demand of the next day and checks the safety; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check; if an emergency occurs, sending an emergency notice according to the emergency, and performing emergency optimization processing; wherein the condition comprises: emergency forces, demand changes, and line faults. The method can solve the problems of opaque information, complex model and the like in the centralized unified optimization method of the energy system, and realizes the cooperative optimization of the comprehensive energy system in the park through a small amount of information interaction between the power system and the natural gas system.
Drawings
Fig. 1 is a schematic flow chart of a campus comprehensive energy optimization method based on an ADMM alternative direction multiplier method according to a first embodiment of the present invention;
FIG. 2 is a schematic flow chart of the pre-optimization process in the optimization method of the park comprehensive energy based on the ADMM alternative direction multiplier method according to the present invention;
fig. 3 is a schematic structural diagram of a campus comprehensive energy optimization device based on the ADMM alternative direction multiplier method according to a second embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be understood that the step numbers used herein are for convenience of description only and are not intended as limitations on the order in which the steps are performed.
It is to be understood that the terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the specification of the present invention and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
The terms "comprises" and "comprising" indicate the presence of the described features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
The term "and/or" refers to and includes any and all possible combinations of one or more of the associated listed items.
The first embodiment of the present invention:
referring to fig. 1-2, a method for optimizing park comprehensive energy based on the ADMM alternative direction multiplier method according to an embodiment of the present invention at least includes the following steps:
s101, obtaining declaration information of a participant, and displaying real-time demand and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant;
in step S101, the declaration information includes: new participants participate in the application, fixed participants quit the application, the requirements of optimized participants, output prediction and the like. And the participants apply to the system operator through the background, are audited by the system operator, and display the real-time demand and output information on a public interface.
S102, according to the declaration information, performing pre-optimization processing by adopting an ADMM alternative direction multiplier method, decomposing an electric power and energy system into two optimization sub-problems of an electric power system and a natural gas system, generating optimization parameters and optimization results by taking the minimum operation cost of the whole park energy system as a target, and sending the optimization parameters and the optimization results to enterprises in the park;
it should be noted that, for step S102, the pre-optimization operation is performed based on the ADMM alternative direction multiplier method, the electric power-natural gas system is decomposed into two optimization sub-problems of the electric power system and the natural gas system, the objective is to minimize the operating cost of the energy system in the whole park, the distributed optimization is realized through a small amount of information transmission between the electric power operator and the natural gas operator, the optimization result preview is generated, the transmission constraint checking condition of the electric power line and the gas pipeline is displayed, the minimum output of each energy supply enterprise when the energy supply requirement of each enterprise at each moment is to be met, and finally the energy flow diagram at each moment is uploaded.
S103, according to the optimization parameters, each enterprise in the park arranges the output and the demand of the next day and checks the safety; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check;
it should be noted that, in step S103, each enterprise sets a notification according to the received optimization parameters, performs output and demand arrangement on the next day, and performs security check. And if the safety check fails, informing the enterprise to reduce the load or perform load peak shifting operation according to the enterprise credit, and performing optimization pre-operation based on the ADMM alternating direction multiplier method again until the safety check is met.
S104, if an emergency occurs, sending an emergency notice according to the emergency, and performing emergency optimization processing; wherein the condition comprises: emergency forces, demand changes, and line faults.
It should be noted that for step 104, if there is an emergency force and demand change, or a line fault. For the former, enterprises directly issue emergency announcements, operators and other enterprises receive the announcements at the same time, and the operators perform emergency optimization at the moment; for the latter, the operator will issue the emergency announcement directly, and the operator will again perform the pre-optimization operation.
Decomposing the power and energy system into two optimization sub-problems of a power system and a natural gas system in step S102; the method specifically comprises the following steps: in the decomposition process, a variable which is common to an electric power system and a natural gas system and is the natural gas consumption of a gas turbine set is selected as a shared variable, and in the natural gas system, the mathematical expression of the consumption characteristic of the gas turbine set is as follows:
Figure BDA0002634148120000091
wherein a, b and c are consumption coefficients of the gas turbine unit; pGNActive power output by the gas turbine unit; gGNThe consumption of natural gas;
in the power system optimization sub-problem, inGNRepresenting the consumption of natural gas, the following mathematical expression is satisfied:
fGN＝gGN
according to the formula, the energy optimization of the electric power and natural gas system is carried out based on the ADMM alternative direction multiplier method.
For step S102, the mathematical expression of the power system optimization sub-problem objective function is:
Figure BDA0002634148120000101
in the formula, WeIs the power system fuel price;
Figure BDA0002634148120000102
the consumption of the generator set at the minimum output is obtained; kGThe number of sections of the output of the generator set is represented; mjThe slope of the j-th piecewise linearization of the generator set and the gas generator set; delta PjActive output is the j section of the generator set; y is a Lagrange multiplier variable; rho is a constant penalty factor;
the constraint condition comprises system active power balance constraint; the unit is subjected to piecewise linearization output constraint; line transmission power constraints; and (4) carrying out piecewise linearization constraint on the consumption function of the gas turbine.
For step S102, the mathematical expression of the natural gas system optimization subproblem objective function is:
Figure BDA0002634148120000103
in the formula, WgRepresenting the price of the gas source; gpIndicating the steam injection quantity of the gas source;
the constraint conditions comprise upper and lower limit constraints of gas sources and gas load steam injection quantity; the method comprises the steps of (1) performing piecewise linearization constraint on a natural gas pipeline airflow function; gas transmission pipeline node air pressure constraint; node air supply balance constraint; the air pressure ratio of the air compressor is restrained by compensating air pressure loss; and (5) natural gas network incidence matrix constraint.
For step S102, performing pre-optimization processing by using an ADMM alternative direction multiplier method according to the declaration information, specifically:
s201, setting an initial value comprising an original residual error r(0)Dual residual s(0)Penalty factor rho, shared variable
Figure BDA0002634148120000111
Lagrange multiplier coefficient
Figure BDA0002634148120000112
S202, respectively calculating an electric power energy flow optimization sub-problem and a natural gas energy flow optimization sub-problem in sequence, taking the kth iteration as an example: firstly, calculating a power energy flow optimization sub-problem, wherein a mathematical expression is as follows:
Figure BDA0002634148120000113
secondly, calculating a natural gas energy flow optimization sub-problem, wherein a mathematical expression is as follows:
Figure BDA0002634148120000114
and updating the Lagrange multiplier variable y after the kth iteration is completed, wherein the updating formula is as follows:
Figure BDA0002634148120000115
and judging convergence;
s203, judging convergence, and if the convergence is judged, outputting a pre-optimization result; if not, entering k +1 times of iteration, and judging a convergence formula as follows:
Figure BDA0002634148120000116
Figure BDA0002634148120000117
in the formula, r(k+1)And s(k+1)Respectively an original residual error and a dual residual error after the (k + 1) th iteration;prianddualrespectively the corresponding maximum allowed deviation;
and S204, performing alternate iteration on the power energy flow optimization subproblem and the natural gas energy flow optimization subproblem until convergence, and outputting a pre-optimization result.
For step S104, if the emergency occurs, sending an emergency notice according to the emergency, and performing emergency optimization operation; wherein the condition comprises: emergency forces, demand changes, and line faults; the method specifically comprises the following steps: if the emergency output or the requirement changes, the enterprise directly issues emergency announcements, an operator and other enterprises receive the announcements at the same time, and the operator performs emergency optimization processing; if a line fault occurs, an operator directly issues an emergency notice and performs pre-optimization processing according to the ADMM alternative direction multiplier method again.
The embodiment of the invention provides a park comprehensive energy optimization method based on an ADMM alternative direction multiplier method, which comprises the following steps: acquiring declaration information of participants, and displaying real-time demand and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant; according to the declaration information, performing pre-optimization processing by adopting an ADMM alternative direction multiplier method, decomposing two optimization sub-problems of an electric power system and a natural gas system by using the electric power and energy system, generating optimization parameters and optimization results by taking the minimum operation cost of the energy system of the whole park as a target, and sending the optimization parameters and the optimization results to enterprises in the park; according to the optimized parameters, each enterprise in the park arranges the output and the demand of the next day and checks the safety; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check; if an emergency occurs, sending an emergency notice according to the emergency, and performing emergency optimization processing; wherein the condition comprises: emergency forces, demand changes, and line faults. The method can solve the problems of opaque information, complex model and the like in the centralized unified optimization method of the energy system, and realizes the cooperative optimization of the comprehensive energy system in the park through a small amount of information interaction between the power system and the natural gas system.
Second embodiment of the invention:
referring to fig. 3, an apparatus 300 for optimizing park integrated energy based on the ADMM alternative direction multiplier method according to an embodiment of the present invention includes: an information acquisition unit 301, a pre-optimization processing unit 302, a security check unit 303, and an emergency processing unit 304, wherein;
the information acquisition unit 301 is configured to acquire declaration information of a participant, and display real-time requirement and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant;
the pre-optimization processing unit 302 is configured to perform pre-optimization processing according to the declaration information by using an ADMM alternative direction multiplier method, decompose an electric power and energy system into two optimization sub-problems of an electric power system and a natural gas system, generate an optimization parameter and an optimization result with the objective that the operating cost of the energy system of the whole park is the minimum, and send the optimization parameter and the optimization result to each enterprise in the park;
the safety check unit 303 is configured to perform output and demand arrangement for each enterprise in the park on the next day according to the optimized parameters, and perform safety check; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check;
the emergency processing unit 304 is configured to send an emergency notice according to an emergency situation if the emergency situation occurs, and perform emergency optimization processing; wherein the condition comprises: emergency forces, demand changes, and line faults.
In one embodiment of the present invention, the power and energy system is divided into two optimization sub-problems of a power system and a natural gas system; the method specifically comprises the following steps: in the decomposition process, a variable which is common to an electric power system and a natural gas system and is the natural gas consumption of a gas turbine set is selected as a shared variable, and in the natural gas system, the mathematical expression of the consumption characteristic of the gas turbine set is as follows:
Figure BDA0002634148120000131
wherein a, b and c are consumption coefficients of the gas turbine unit; pGNActive power output by the gas turbine unit; gGNIs the natural gas consumption.
In the power system optimization sub-problem, inGNRepresenting the consumption of natural gas, the following mathematical expression is satisfied:
fGN＝gGN
according to the formula, the energy optimization of the electric power and natural gas system is carried out based on the ADMM alternative direction multiplier method.
In one embodiment of the present invention, the mathematical expression of the objective function of the power system optimization sub-problem is:
Figure BDA0002634148120000132
in the formula, WeIs the power system fuel price;
Figure BDA0002634148120000141
the consumption of the generator set at the minimum output is obtained; kGThe number of sections of the output of the generator set is represented; mjThe slope of the j-th piecewise linearization of the generator set and the gas generator set; delta PjActive output is the j section of the generator set; y is a Lagrange multiplier variable; rho is a constant penalty factor;
the constraint condition comprises system active power balance constraint; the unit is subjected to piecewise linearization output constraint; line transmission power constraints; the gas turbine consumption function is subjected to piecewise linearization constraint;
the mathematical expression of the natural gas system optimization subproblem objective function is as follows:
Figure BDA0002634148120000142
in the formula, WgRepresenting the price of the gas source; gpIndicating the steam injection quantity of the gas source;
the constraint conditions comprise upper and lower limit constraints of gas sources and gas load steam injection quantity; the method comprises the steps of (1) performing piecewise linearization constraint on a natural gas pipeline airflow function; gas transmission pipeline node air pressure constraint; node air supply balance constraint; the air pressure ratio of the air compressor is restrained by compensating air pressure loss; and (5) natural gas network incidence matrix constraint.
The embodiment of the invention provides a campus comprehensive energy optimization device 300 based on an ADMM alternative direction multiplier method, which comprises: an information acquisition unit 301, a pre-optimization processing unit 302, a security check unit 303, and an emergency processing unit 304, wherein; the information acquisition unit 301 is configured to acquire declaration information of a participant, and display real-time requirement and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant; the pre-optimization processing unit 302 is configured to perform pre-optimization processing according to the declaration information by using an ADMM alternative direction multiplier method, decompose an electric power and energy system into two optimization sub-problems of an electric power system and a natural gas system, generate an optimization parameter and an optimization result with the objective that the operating cost of the energy system of the whole park is the minimum, and send the optimization parameter and the optimization result to each enterprise in the park; the safety check unit 303 is configured to perform output and demand arrangement for each enterprise in the park on the next day according to the optimized parameters, and perform safety check; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check; the emergency processing unit 304 is configured to send an emergency notice according to an emergency situation if the emergency situation occurs, and perform emergency optimization processing; wherein the condition comprises: emergency forces, demand changes, and line faults. The device can solve the problems of opaque information, complex model and the like in the centralized unified optimization method of the energy system, and realizes the cooperative optimization of the comprehensive energy system in the park through a small amount of information interaction between the power system and the natural gas system.
Third embodiment of the invention:
an embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, wherein the computer program, when executed by a processor, implements the method for optimizing park integrated energy based on the ADMM alternative direction multiplier method according to any one of claims 1 to 6.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and substitutions can be made without departing from the technical principle of the present invention, and these modifications and substitutions should also be regarded as the protection scope of the present invention.

Claims (10)

1. A park comprehensive energy optimization method based on an ADMM alternating direction multiplier method is characterized by comprising the following steps:
acquiring declaration information of participants, and displaying real-time demand and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant;
according to the declaration information, performing pre-optimization processing by adopting an ADMM alternative direction multiplier method, decomposing two optimization sub-problems of an electric power system and a natural gas system by using the electric power and energy system, generating optimization parameters and optimization results by taking the minimum operation cost of the energy system of the whole park as a target, and sending the optimization parameters and the optimization results to enterprises in the park;
according to the optimized parameters, each enterprise in the park arranges the output and the demand of the next day and checks the safety; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check;
if an emergency occurs, sending an emergency notice according to the emergency, and performing emergency optimization processing; wherein the condition comprises: emergency forces, demand changes, and line faults.
2. The method of claim 1, wherein the decomposing of the power and energy system into two optimization sub-problems, power system and natural gas system; the method specifically comprises the following steps: in the decomposition process, a variable which is common to an electric power system and a natural gas system and is the natural gas consumption of a gas turbine set is selected as a shared variable, and in the natural gas system, the mathematical expression of the consumption characteristic of the gas turbine set is as follows:
Figure FDA0002634148110000011
wherein a, b and c are consumption coefficients of the gas turbine unit; pGNActive power output by the gas turbine unit; gGNThe consumption of natural gas;
in the power system optimization sub-problem, inGNRepresenting the consumption of natural gas, the following mathematical expression is satisfied:
fGN＝gGN
according to the formula, the energy optimization of the electric power and natural gas system is carried out based on the ADMM alternative direction multiplier method.
3. The ADMM alternating direction multiplier method-based campus comprehensive energy optimization method according to claim 2, wherein the mathematical expression of the power system optimization sub-problem objective function is:
Figure FDA0002634148110000021
in the formula, WeIs the power system fuel price;
Figure FDA0002634148110000022
the consumption of the generator set at the minimum output is obtained; kGThe number of sections of the output of the generator set is represented; mjThe slope of the j-th piecewise linearization of the generator set and the gas generator set; delta PjActive output is the j section of the generator set; y is a Lagrange multiplier variable; rho is a constant penalty factor;
the constraint condition comprises system active power balance constraint; the unit is subjected to piecewise linearization output constraint; line transmission power constraints; and (4) carrying out piecewise linearization constraint on the consumption function of the gas turbine.
4. The method of claim 2, wherein the mathematical expression of the natural gas system optimization subproblem objective function is:
Figure FDA0002634148110000031
in the formula, WgRepresenting the price of the gas source; gpIndicating the steam injection quantity of the gas source;
the constraint conditions comprise upper and lower limit constraints of gas sources and gas load steam injection quantity; the method comprises the steps of (1) performing piecewise linearization constraint on a natural gas pipeline airflow function; gas transmission pipeline node air pressure constraint; node air supply balance constraint; the air pressure ratio of the air compressor is restrained by compensating air pressure loss; and (5) natural gas network incidence matrix constraint.
5. The method for optimizing the comprehensive energy of the campus based on the ADMM alternative direction multiplier method as claimed in claim 1, wherein the pre-optimization processing is performed by the ADMM alternative direction multiplier method according to the declaration information, and specifically comprises:
setting an initial value comprising the original residual r(0)Dual residual s(0)Penalty factor rho, shared variable
Figure FDA0002634148110000032
Lagrange multiplier coefficient y(0)；
Sequentially and respectively calculating an electric power energy flow optimization subproblem and a natural gas energy flow optimization subproblem, taking the kth iteration as an example: firstly, calculating a power energy flow optimization sub-problem, wherein a mathematical expression is as follows:
Figure FDA0002634148110000033
secondly, calculating a natural gas energy flow optimization sub-problem, wherein a mathematical expression is as follows:
Figure FDA0002634148110000034
and updating the Lagrange multiplier variable y after the kth iteration is completed, wherein the updating formula is as follows:
Figure FDA0002634148110000035
and judging convergence;
judging the convergence, and if the convergence is reached, outputting a pre-optimization result; if not, entering k +1 times of iteration, and judging a convergence formula as follows:
Figure FDA0002634148110000036
Figure FDA0002634148110000041
in the formula, r(k+1)And s(k+1)Respectively an original residual error and a dual residual error after the (k + 1) th iteration;prianddualrespectively the corresponding maximum allowed deviation;
and performing alternate iteration on the power energy flow optimization subproblem and the natural gas energy flow optimization subproblem until convergence, and outputting a pre-optimization result.
6. The method of claim 1, wherein if an emergency occurs, an emergency notice is issued and an emergency optimization operation is performed according to the emergency; wherein the condition comprises: emergency forces, demand changes, and line faults; the method specifically comprises the following steps: if the emergency output or the requirement changes, the enterprise directly issues emergency announcements, an operator and other enterprises receive the announcements at the same time, and the operator performs emergency optimization processing; if a line fault occurs, an operator directly issues an emergency notice and performs pre-optimization processing according to the ADMM alternative direction multiplier method again.
7. A park comprehensive energy optimization device based on an ADMM alternative direction multiplier method is characterized by comprising the following steps: the system comprises an information acquisition unit, a pre-optimization processing unit, a safety check unit and an emergency processing unit, wherein the information acquisition unit is used for acquiring information;
the information acquisition unit is used for acquiring declaration information of participants and displaying real-time requirement and output information on a public interface; wherein the declaration information includes: the method comprises the following steps that a new participant participates in application, a fixed participant quits the application, and the demand and output prediction of each optimized participant;
the pre-optimization processing unit is used for performing pre-optimization processing by adopting an ADMM alternative direction multiplier method according to the declaration information, decomposing the electric power and energy system into two optimization sub-problems of an electric power system and a natural gas system, generating optimization parameters and optimization results by taking the minimum operation cost of the energy system of the whole park as a target, and sending the optimization parameters and the optimization results to each enterprise of the park;
the safety checking unit is used for carrying out output and demand arrangement on each enterprise in the park on the next day according to the optimized parameters and carrying out safety checking; if the safety check does not pass, informing the enterprise to reduce the load or carry out load peak shifting operation according to enterprise credit, and carrying out optimization processing according to the ADMM alternating direction multiplier method again until the output and demand arrangement of the enterprise meets the safety check;
the emergency processing unit is used for sending an emergency notice according to the emergency and carrying out emergency optimization processing if the emergency occurs; wherein the condition comprises: emergency forces, demand changes, and line faults.
8. The apparatus of claim 7, wherein the power and energy system is split into two optimization sub-problems, power system and natural gas system; the method specifically comprises the following steps: in the decomposition process, a variable which is common to an electric power system and a natural gas system and is the natural gas consumption of a gas turbine set is selected as a shared variable, and in the natural gas system, the mathematical expression of the consumption characteristic of the gas turbine set is as follows:
Figure FDA0002634148110000051
wherein a, b and c are consumption coefficients of the gas turbine unit; pGNActive power output by the gas turbine unit; gGNIs the natural gas consumption.
In the power system optimization sub-problem, inGNRepresenting the consumption of natural gas, the following mathematical expression is satisfied:
fGN＝gGN
according to the formula, the energy optimization of the electric power and natural gas system is carried out based on the ADMM alternative direction multiplier method.
9. The apparatus of claim 8, wherein the mathematical expression of the objective function of the optimization subproblem of the power system is:
Figure FDA0002634148110000061
in the formula, WeIs the power system fuel price;
Figure FDA0002634148110000062
the consumption of the generator set at the minimum output is obtained; kGThe number of sections of the output of the generator set is represented; mjThe slope of the j-th piecewise linearization of the generator set and the gas generator set; delta PjActive output is the j section of the generator set; y is a Lagrange multiplier variable; rho is a constant penalty factor;
the constraint condition comprises system active power balance constraint; the unit is subjected to piecewise linearization output constraint; line transmission power constraints; the gas turbine consumption function is subjected to piecewise linearization constraint;
the mathematical expression of the natural gas system optimization subproblem objective function is as follows:
Figure FDA0002634148110000063
in the formula, WgRepresenting the price of the gas source; gpIndicating the steam injection quantity of the gas source;
the constraint conditions comprise upper and lower limit constraints of gas sources and gas load steam injection quantity; the method comprises the steps of (1) performing piecewise linearization constraint on a natural gas pipeline airflow function; gas transmission pipeline node air pressure constraint; node air supply balance constraint; the air pressure ratio of the air compressor is restrained by compensating air pressure loss; and (5) natural gas network incidence matrix constraint.
10. A computer-readable storage medium storing a computer program, wherein the computer program when executed by a processor implements the method for energy optimization of a campus based on the ADMM alternative direction multiplier as recited in any one of claims 1 to 6.
CN202010823959.0A 2020-08-14 2020-08-14 Park comprehensive energy optimization method and device based on ADMM alternating direction multiplier method Pending CN112036613A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010823959.0A CN112036613A (en) 2020-08-14 2020-08-14 Park comprehensive energy optimization method and device based on ADMM alternating direction multiplier method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010823959.0A CN112036613A (en) 2020-08-14 2020-08-14 Park comprehensive energy optimization method and device based on ADMM alternating direction multiplier method

Publications (1)

Publication Number Publication Date
CN112036613A true CN112036613A (en) 2020-12-04

Family

ID=73577405

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010823959.0A Pending CN112036613A (en) 2020-08-14 2020-08-14 Park comprehensive energy optimization method and device based on ADMM alternating direction multiplier method

Country Status (1)

Country Link
CN (1) CN112036613A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112767148A (en) * 2021-01-08 2021-05-07 合肥工业大学 Electric power trading method and system based on block chain and distributed optimization
CN113507127A (en) * 2021-09-08 2021-10-15 国网江西省电力有限公司电力科学研究院 Power distribution network distributed voltage regulation method and device based on alternating direction multiplier method
CN113762632A (en) * 2021-09-10 2021-12-07 国网四川省电力公司经济技术研究院 Collaborative optimization operation method and system of electrical comprehensive energy system
CN116433225A (en) * 2023-06-12 2023-07-14 国网湖北省电力有限公司经济技术研究院 Multi-time scale fault recovery method, device and equipment for interconnected micro-grid

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150199606A1 (en) * 2014-01-14 2015-07-16 Mitsubishi Electric Research Laboratories, Inc. System and Method for Optimal Power Flow Analysis
CN107292456A (en) * 2017-08-01 2017-10-24 重庆大学 Electrical energy flow point cloth collaboration optimized calculation method based on alternating direction multiplier method
CN107508284A (en) * 2017-08-15 2017-12-22 华北电力大学 The micro-capacitance sensor distributed optimization dispatching method of meter and electrical interconnection
CN110866773A (en) * 2019-09-29 2020-03-06 西安交通大学 Electric-gas combined market distributed clearing method based on alternative direction multiplier method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150199606A1 (en) * 2014-01-14 2015-07-16 Mitsubishi Electric Research Laboratories, Inc. System and Method for Optimal Power Flow Analysis
CN107292456A (en) * 2017-08-01 2017-10-24 重庆大学 Electrical energy flow point cloth collaboration optimized calculation method based on alternating direction multiplier method
CN107508284A (en) * 2017-08-15 2017-12-22 华北电力大学 The micro-capacitance sensor distributed optimization dispatching method of meter and electrical interconnection
CN110866773A (en) * 2019-09-29 2020-03-06 西安交通大学 Electric-gas combined market distributed clearing method based on alternative direction multiplier method

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112767148A (en) * 2021-01-08 2021-05-07 合肥工业大学 Electric power trading method and system based on block chain and distributed optimization
CN112767148B (en) * 2021-01-08 2022-09-06 合肥工业大学 Electric power trading method and system based on block chain and distributed optimization
CN113507127A (en) * 2021-09-08 2021-10-15 国网江西省电力有限公司电力科学研究院 Power distribution network distributed voltage regulation method and device based on alternating direction multiplier method
CN113762632A (en) * 2021-09-10 2021-12-07 国网四川省电力公司经济技术研究院 Collaborative optimization operation method and system of electrical comprehensive energy system
CN116433225A (en) * 2023-06-12 2023-07-14 国网湖北省电力有限公司经济技术研究院 Multi-time scale fault recovery method, device and equipment for interconnected micro-grid
CN116433225B (en) * 2023-06-12 2023-08-29 国网湖北省电力有限公司经济技术研究院 Multi-time scale fault recovery method, device and equipment for interconnected micro-grid

Similar Documents

Publication Publication Date Title
CN112036613A (en) Park comprehensive energy optimization method and device based on ADMM alternating direction multiplier method
Mosa et al. Energy management system of low voltage dc microgrid using mixed-integer nonlinear programing and a global optimization technique
CN104102212B (en) Dispatching method, apparatus and system for gas and steam system in iron and steel enterprises
CN101179195B (en) Power distribution network planning scheme assistant decision system
Ameli et al. Coordinated operation strategies for natural gas and power systems in presence of gas‐related flexibilities
CN104851052A (en) Intelligent error preventing method and system used for outputting dispatching operation ticket
Xuan et al. Two-stage planning for electricity-gas coupled integrated energy system with carbon capture, utilization, and storage considering carbon tax and price uncertainties
CN111199354A (en) Hopfield neural network-based power communication network operation and maintenance work order scheduling method
Cunningham et al. Application of delay-time analysis via Monte Carlo simulation
Wu et al. Reserve constrained dynamic economic dispatch with valve-point effect: A two-stage mixed integer linear programming approach
Mi et al. Operation of natural gas and electricity networks with line pack
CN114169669A (en) Power generation industry carbon emission prediction method, platform, computing equipment and medium
CN115936722A (en) Data management method of intelligent gas call center, internet of things system and medium
Liang et al. A data-driven convex model for hybrid microgrid operation with bidirectional converters
Yang et al. A fast calculation method for long-term security-constrained unit commitment of large-scale power systems with renewable energy
CN116565859A (en) Power grid optimal scheduling system and method based on artificial intelligence
López-Lezama et al. A contingency-based security-constrained optimal power flow model for revealing the marginal cost of a blackout risk-equalizing policy in the colombian electricity market
Yuan et al. Fast solution method for TCUC with long time horizon based on horizon splitting
Baring-Gould et al. Hybrid2: the hybrid power system simulation model
Nepomuceno et al. Equivalent optimization model for loss minimization: a suitable analysis approach
Huang et al. A new method for dynamic economic emission dispatch problem
CN105488733A (en) Power business expansion engineering data processing apparatus and method
van Amelsvoort SG-rating–Putting values on smart grid architectures
Claeys Distribution network modeling: from simulation towards optimization
Blas et al. Improving the optimization of electric power systems through a discrete event based simulation model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20201204

RJ01 Rejection of invention patent application after publication